AliPhysics  421aab4 (421aab4)
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
AliJetResponseMaker.h
Go to the documentation of this file.
1 #ifndef ALIJETRESPONSEMAKER_H
2 #define ALIJETRESPONSEMAKER_H
3 /**************************************************************************
4 * Copyright(c) 1998-2009, ALICE Experiment at CERN, All rights reserved. *
5 * *
6 * Author: The ALICE Off-line Project. *
7 * Contributors are mentioned in the code where appropriate. *
8 * *
9 * Permission to use, copy, modify and distribute this software and its *
10 * documentation strictly for non-commercial purposes is hereby granted *
11 * without fee, provided that the above copyright notice appears in all *
12 * copies and that both the copyright notice and this permission notice *
13 * appear in the supporting documentation. The authors make no claims *
14 * about the suitability of this software for any purpose. It is *
15 * provided "as is" without express or implied warranty. *
16 **************************************************************************/
17 
18 //-----------------------------------------------------------------------
19 // Author : Salvatore Aiola, Yale University, salvatore.aiola@cern.ch
20 //-----------------------------------------------------------------------
21 
22 class TClonesArray;
23 class TH2;
24 class THnSparse;
25 class AliNamedArrayI;
26 
27 #include "AliEmcalJet.h"
29 #include "AliEmcalEmbeddingQA.h"
30 
32  public:
34  AliJetResponseMaker(const char *name);
35  virtual ~AliJetResponseMaker();
36 
40  kMCLabel = 2,
42  };
43 
45 
49  void SetMinJetMCPt(Float_t pt) { fMinJetMCPt = pt ; }
50  void SetHistoType(Int_t b) { fHistoType = b ; }
51  void SetDeltaPtAxis(Int_t b) { fDeltaPtAxis = b ; }
53  void SetNEFAxis(Int_t b) { fNEFAxis = b ; }
54  void SetZAxis(Int_t b) { fZAxis = b ; }
57  void SetZgAxis(Int_t b) { fZgAxis = b ; }
58  void SetdRAxis(Int_t b) { fdRAxis = b ; }
59  void SetPtgAxis(Int_t b) { fPtgAxis = b ; }
60  void SetDBCAxis(Int_t b) { fDBCAxis = b ; }
62 
63  protected:
64  void ExecOnce();
65  void DoJetLoop();
67  Bool_t Run();
69  void SetMatchingLevel(AliEmcalJet *jet1, AliEmcalJet *jet2, MatchingType matching);
70  void GetGeometricalMatchingLevel(AliEmcalJet *jet1, AliEmcalJet *jet2, Double_t &d) const;
71  void GetMCLabelMatchingLevel(AliEmcalJet *jet1, AliEmcalJet *jet2, Double_t &d1, Double_t &d2) const;
72  void GetSameCollectionsMatchingLevel(AliEmcalJet *jet1, AliEmcalJet *jet2, Double_t &d1, Double_t &d2) const;
73  void FillMatchingHistos(AliEmcalJet* jet1, AliEmcalJet* jet2, Double_t d, Double_t CE1, Double_t CE2);
74  void FillJetHisto(AliEmcalJet* jet, Int_t Set);
75  void AllocateTH2();
76  void AllocateTHnSparse();
77  Double_t GetRelativeEPAngle(Double_t jetAngle, Double_t epAngle) const;
78 
79  MatchingType fMatching; // matching type
80  Double_t fMatchingPar1; // matching parameter for jet1-jet2 matching
81  Double_t fMatchingPar2; // matching parameter for jet2-jet1 matching
82  Bool_t fUseCellsToMatch; // use cells instead of clusters to match jets (slower but sometimes needed)
83  Double_t fMinJetMCPt; // minimum jet MC pt
85  Int_t fHistoType; // histogram type (0=TH2, 1=THnSparse)
86  Int_t fDeltaPtAxis; // add delta pt axis in THnSparse (default=0)
87  Int_t fDeltaEtaDeltaPhiAxis; // add delta eta and delta phi axes in THnSparse (default=0)
88  Int_t fNEFAxis; // add NEF axis in matching THnSparse (default=0)
89  Int_t fZAxis; // add Z axis in matching THnSparse (default=0)
90  Int_t fFlavourZAxis; // add flavour Z axis in matching THnSparse (default=0)
91  Int_t fFlavourPtAxis; // add flavour pt axis in matching THnSparse (default=0)
92  Int_t fZgAxis; // add Zg axis in matching THnSparse (default=0)
93  Int_t fdRAxis; // add dR axis in matching THnSparse (default=0)
94  Int_t fPtgAxis; // add Ptg axis in matching THnSparse (default=0)
95  Int_t fDBCAxis; // add DBC (number of soft dropped branches) axis in matching THnSparse (default=0)
97 
100 
103 
104  // THnSparse
105  THnSparse *fHistJets1;
106  THnSparse *fHistJets2;
107  THnSparse *fHistMatching;
108 
109  // Jets 1
116 
117  // Jets 2
124 
125  // Jet1-Jet2 matching
134 
148 
160 
172 
173  private:
174  AliJetResponseMaker(const AliJetResponseMaker&); // not implemented
175  AliJetResponseMaker &operator=(const AliJetResponseMaker&); // not implemented
176 
177  ClassDef(AliJetResponseMaker, 29) // Jet response matrix producing task
178 };
179 #endif
TH2 * fHistJets2PtArea
phi-eta distribution of jets 2
TH2 * fHistRejectionReason1
whether the jet2 collection has to be average subtracted
TH2 * fHistJets2CorrPtArea
inclusive jet pt vs. area histogram 2
double Double_t
Definition: External.C:58
TH2 * fHistDeltaPtvsCommonEnergy1
delta pt between matched jets vs distance
TH2 * fHistJet2PtOverJet1PtvsJet2Pt
delta eta vs delta phi of matched jets
Bool_t FillHistograms()
Function filling histograms.
TH2 * fHistDeltaPtvsDeltaArea
delta pt between matched jets vs jet 2 area
void SetFlavourPtAxis(Int_t b)
TH2 * fHistJets2PhiEta
Constituent Pt over Jet Pt ratio vs. jet pt 1.
TH2 * fHistDeltaCorrPtvsArea2
delta pt corr between matched jets vs jet 1 area
void FillJetHisto(AliEmcalJet *jet, Int_t Set)
void SetDeltaPtAxis(Int_t b)
Int_t fJetRelativeEPAngle
add jet angle relative to the EP in matching THnSparse (default=0)
AliEmcalEmbeddingQA fEmbeddingQA
! Embedding QA hists (will only be added if embedding)
void GetSameCollectionsMatchingLevel(AliEmcalJet *jet1, AliEmcalJet *jet2, Double_t &d1, Double_t &d2) const
TH2 * fHistDeltaMCPtvsCommonEnergy2
jet 1 MC pt - jet2 pt vs common energy 1 (%)
void SetMinJetMCPt(Float_t pt)
TH2 * fHistJets1ZvsPt
Jet charged energy fraction vs. charged jet pt 1.
Double_t GetRelativeEPAngle(Double_t jetAngle, Double_t epAngle) const
THnSparse * fHistMatching
jet2 THnSparse
TH2 * fHistDeltaEtaDeltaPhi
common energy 1 (%) vs common energy 2 (%)
TH2 * fHistDeltaMCPtvsCommonEnergy1
jet 1 MC pt - jet2 pt vs distance
TH2 * fHistDeltaCorrPtOverJet2CorrPtvsJet2CorrPt
delta pt corr / jet 1 corr pt between matched jets vs jet 1 corr pt
TH2 * fHistCommonEnergy1vsCommonEnergy2
distance vs common energy 2 (%)
TH2 * fHistDeltaPtvsCommonEnergy2
delta pt between matched jets vs common energy 1 (%)
TH2 * fHistDeltaMCPtvsDistance
jet 1 MC pt - jet2 pt vs jet 2 pt
void SetJetRelativeEPAngleAxis(Int_t b)
QA Class for EMCal Embedding Framework.
TH2 * fHistJets1NEFvsPt
inclusive jet pt vs. area histogram 1
TH2 * fHistJets2ZvsPt
Jet charged energy fraction vs. charged jet pt 2.
Bool_t fIsJet2Rho
whether the jet1 collection has to be average subtracted
TH2 * fHistDeltaPtvsJet1Pt
jet 1 pt over jet 2 pt vs jet 1 pt
TH2 * fHistDeltaCorrPtvsCommonEnergy2
delta pt corr between matched jets vs common energy 1 (%)
TH2 * fHistDeltaPtvsJet2Pt
delta pt between matched jets vs jet 1 pt
TH2 * fHistJets1PhiEta
matching THnSparse
int Int_t
Definition: External.C:63
THnSparse * fHistJets1
Rejection reason vs. jet pt.
float Float_t
Definition: External.C:68
TH2 * fHistJet1PtOverJet2PtvsJet1Pt
jet 2 pt over jet 1 pt vs jet 2 pt
void SetPtHardBin(Int_t b)
TH2 * fHistDeltaCorrPtvsDistance
delta pt corr between matched jets vs jet 2 corr pt
void FillMatchingHistos(AliEmcalJet *jet1, AliEmcalJet *jet2, Double_t d, Double_t CE1, Double_t CE2)
void ExecOnce()
Perform steps needed to initialize the analysis.
TH2 * fHistJet1PtvsJet2Pt
delta pt between matched jets vs delta area
TH2 * fHistDeltaCorrPtvsJet1CorrPt
delta pt corr / jet 2 corr pt between matched jets vs jet 2 corr pt
TH2 * fHistDeltaMCPtOverJet2PtvsJet2Pt
jet 1 MC pt - jet2 pt / jet 1 MC pt vs jet 1 pt
TH2 * fHistDistancevsCommonEnergy2
distance vs common energy 1 (%)
void SetMatching(MatchingType t, Double_t p1=1, Double_t p2=1)
TH2 * fHistRejectionReason2
Rejection reason vs. jet pt.
TH2 * fHistDistancevsJet1Pt
common energy 2 (%) vs jet 2 pt
TH2 * fHistDeltaMCPtOverJet1MCPtvsJet1MCPt
correlation jet 1 corr pt vs jet 2 corr pt
Int_t fSelectPtHardBin
select one pt hard bin for analysis
TH2 * fHistDeltaPtOverJet1PtvsJet1Pt
delta pt between matched jets vs jet 2 pt
TH2 * fHistDeltaCorrPtvsArea1
delta pt corr between matched jets vs common energy 2 (%)
Bool_t Run()
Run function. This is the core function of the analysis and contains the user code. Therefore users have to implement this function.
TH2 * fHistCommonEnergy1vsJet1Pt
Constituent Pt over Jet Pt ratio vs. jet pt 2.
void SetMatchingLevel(AliEmcalJet *jet1, AliEmcalJet *jet2, MatchingType matching)
TH2 * fHistJets1CorrPtArea
inclusive jet pt vs. area histogram 1
TH2 * fHistDeltaMCPtvsJet1MCPt
jet 1 MC pt - jet2 pt / jet 2 pt vs jet 2 pt
TH2 * fHistDeltaMCPtvsArea2
jet 1 MC pt - jet2 pt vs jet 1 area
void SetHistoType(Int_t b)
TH2 * fHistDeltaMCPtvsJet2Pt
jet 1 MC pt - jet2 pt vs jet 1 MC pt
Definition: External.C:220
TH2 * fHistDistancevsJet2Pt
distance vs jet 1 pt
TH2 * fHistDeltaPtvsArea1
delta pt between matched jets vs common energy 2 (%)
TH2 * fHistDeltaCorrPtvsJet2CorrPt
delta pt corr between matched jets vs jet 1 corr pt
void GetGeometricalMatchingLevel(AliEmcalJet *jet1, AliEmcalJet *jet2, Double_t &d) const
THnSparse * fHistJets2
jet1 THnSparse
TH2 * fHistDeltaPtOverJet2PtvsJet2Pt
delta pt / jet 1 pt between matched jets vs jet 1 pt
TH2 * fHistDeltaPtvsArea2
delta pt between matched jets vs jet 1 area
TH2 * fHistJets2NEFvsPt
inclusive jet pt vs. area histogram 2
void SetDeltaEtaDeltaPhiAxis(Int_t b)
Base task in the EMCAL jet framework.
Represent a jet reconstructed using the EMCal jet framework.
Definition: AliEmcalJet.h:44
TH2 * fHistCommonEnergy2vsJet2Pt
common energy 1 (%) vs jet 1 pt
TH2 * fHistDeltaPtvsDistance
delta pt / jet 2 pt between matched jets vs jet 2 pt
TH2 * fHistJet1MCPtvsJet2Pt
jet 1 MC pt - jet2 pt vs delta area
TH2 * fHistJets2CEFvsCEFPt
Jet neutral energy fraction vs. jet pt 2.
TH2 * fHistDistancevsCommonEnergy1
distance vs jet 2 pt
TH2 * fHistJets1CEFvsCEFPt
Jet neutral energy fraction vs. jet pt 1.
bool Bool_t
Definition: External.C:53
void SetUseCellsToMatch(Bool_t i)
TH2 * fHistDeltaMCPtvsArea1
jet 1 MC pt - jet2 pt vs common energy 2 (%)
TH2 * fHistDeltaCorrPtOverJet1CorrPtvsJet1CorrPt
correlation jet 1 pt vs jet 2 pt
TH2 * fHistJets1PtArea
phi-eta distribution of jets 1
TH2 * fHistDeltaMCPtvsDeltaArea
jet 1 MC pt - jet2 pt vs jet 2 area
TH2 * fHistDeltaCorrPtvsCommonEnergy1
delta pt corr between matched jets vs distance
void GetMCLabelMatchingLevel(AliEmcalJet *jet1, AliEmcalJet *jet2, Double_t &d1, Double_t &d2) const
AliJetResponseMaker & operator=(const AliJetResponseMaker &)
void SetFlavourZAxis(Int_t b)
TH2 * fHistJet1CorrPtvsJet2CorrPt
delta pt corr between matched jets vs delta area
TH2 * fHistDeltaCorrPtvsDeltaArea
delta pt corr between matched jets vs jet 2 area