AliPhysics  9b6b435 (9b6b435)
AliAnaParticleJetLeadingConeCorrelation.h
Go to the documentation of this file.
1 #ifndef ALIANAPARTICLEJETLEADINGCONECORRELATION_H
2 #define ALIANAPARTICLEJETLEADINGCONECORRELATION_H
3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4  * See cxx source for full Copyright notice */
5 
6 //_________________________________________________________________________
25 
26 
27 // --- ROOT system ---
28 class TH2F;
29 
30 //---- Analysis system ----
31 class AliAODTrack;
32 class AliVCluster;
33 class AliCaloTrackReader;
35 
37 
39 
40 public:
41 
43 
46 
48 
49  void InitParameters();
50 
51  void Print(const Option_t * opt) const;
52 
53  Bool_t AreJetsRecalculated() const { return fReMakeJet ; }
55  void SwitchOffJetsRecalculation() { fReMakeJet = kFALSE ; }
56 
57  Bool_t AreJetsOnlyInCTS() const { return fJetsOnlyInCTS ; }
59  void SwitchOffJetsOnlyInCTS() { fJetsOnlyInCTS = kFALSE ; }
60 
64 
65  Bool_t IsPbPb() const { return fPbPb ; }
66  void SetppCollisions() { fPbPb = kFALSE ; }
67  void SetPbPbCollisions() { fPbPb = kTRUE ; }
68 
73 
77 
79  { fPtTriggerSelectionCut = cut ; }
80  void SetJetSelectionMode(UInt_t sel) { fSelect= sel ; }
81 
82  Int_t GetJetNCones() const { return fJetNCone ; }
83  Int_t GetJetNPtThres() const { return fJetNPt ; }
84  Float_t GetJetCone() const { return fJetCone ; }
87  Float_t GetJetCones(Int_t i) const { return fJetCones[i] ; }
88  Float_t GetJetPtThreshold(Int_t i) const { return fJetPtThres[i] ; }
89  TString GetJetConeName(Int_t i) const { return fJetNameCones[i] ; }
91 
92 
95  void SetLeadingRatioCutRange(Double_t ratiomin, Double_t ratiomax)
96  { fLeadingRatioMaxCut =ratiomax; fLeadingRatioMinCut = ratiomin ; }
97 
98  void SetJetNCones(Int_t n) { fJetNCone = n ; }
99  void SetJetNPtThresholds(Int_t n) { fJetNPt = n ; }
100  void SetJetCones(Int_t i, Float_t cone, TString sc)
101  { fJetCones[i] = cone ; fJetNameCones[i] = sc ; }
102  void SetCone(Float_t cone) { fJetCone = cone ; }
106  { fJetPtThres[i] = pt ; fJetNamePtThres[i] = spt ; }
107 
108  void SetJetRatioCutRange(Double_t ratiomin, Double_t ratiomax)
109  { fJetRatioMaxCut = ratiomax ; fJetRatioMinCut = ratiomin ; }
110  void SetJetCTSRatioCutRange(Double_t ratiomin, Double_t ratiomax)
111  { fJetCTSRatioMaxCut = ratiomax ; fJetCTSRatioMinCut = ratiomin ; }
112 
113  Bool_t OnlyIsolated() const { return fSelectIsolated ; }
114  void SelectIsolated(Bool_t select) { fSelectIsolated = select ; }
115 
116  private:
117 
118  Double_t CalculateJetRatioLimit(Double_t ptTrig, const Double_t *param, const Double_t *x) const ;
119 
120  void FillJetHistos(AliCaloTrackParticleCorrelation * particle, const TLorentzVector jet,
121  const TString & type, const TString & lastname);
122 
123  TList * GetOutputContainer() const { return fOutCont ; }
124 
125  Bool_t IsJetSelected(Double_t ptTrig, Double_t ptjet) const ;
126  Bool_t IsParticleInJetCone(Double_t eta, Double_t phi, Double_t etal, Double_t phil) const ;
127 
131 
132  void MakeAnalysisFillAOD();
136 
137  Bool_t SelectCluster(AliVCluster * calo, Double_t *vertex, TLorentzVector & mom, Int_t & pdg) ;
138 
139  private:
140 
145 
146  // Leading particle selection parameters
151 
152  // Jet selection parameters
153  // Fixed cuts (old)
158 
159  // Cuts depending on jet pt
170 
183 
184  // Temporary kinematic containers
185  TVector3 fTrackVector;
186  TLorentzVector fBkgMom;
187  TLorentzVector fJetMom;
188  TLorentzVector fJetConstMom;
189  TLorentzVector fLeadingMom;
190  TLorentzVector fLeadingPi0Mom;
191  TLorentzVector fLeadingPhoMom1;
192  TLorentzVector fLeadingPhoMom2;
193  TLorentzVector fLeadingChargeMom;
194 
195  // Histograms
196  // Leading particle distributions
198 
206 
216 
221 
222  // Jet distributions
223  // Fixed cone and pt threshold
235 
247 
248  // Variable cone and pt threshold
249 
250  TH2F * fhJetPts [5][5];
251  TH2F * fhJetRatioPts [5][5];
257  TH2F * fhJetFFzs [5][5];
258  TH2F * fhJetFFxis[5][5];
259  TH2F * fhJetFFpts[5][5];
261 
262  TH2F * fhBkgPts [5][5];
263  TH2F * fhBkgRatioPts [5][5];
269  TH2F * fhBkgFFzs [5][5];
270  TH2F * fhBkgFFxis[5][5];
271  TH2F * fhBkgFFpts[5][5];
273 
276 
279 
283 
284 } ;
285 
286 #endif //ALIANAPARTICLEJETLEADINGCONECORRELATION_H
287 
288 
289 
TH2F * fhBkgFFz
! Accepted reconstructed background fragmentation function, z=ptjet/pttrig
Double_t fJetXMax2[6]
X Factor to set jet max limit for PbPb.
Int_t pdg
Double_t fJetPtThres[5]
Jet pT threshold under study(fSeveralConeAndPtCuts)
Double_t fJetSigma2[2]
Rec. sigma of jet energy parameters.
Int_t fJetNCone
Number of jet cones sizes, maximum 5.
Double_t fJetPtThreshold
Jet pT threshold under study(!fSeveralConeAndPtCuts)
void MakeAODJet(AliCaloTrackParticleCorrelation *particle)
double Double_t
Definition: External.C:58
TH2F * fhChargedLeadingDeltaPhiRatioPt30
! Difference of charged hadron and trigger phi as function of pT leading / trigger pT...
Definition: External.C:236
Double_t fJetXMin2[6]
X Factor to set jet min limit for PbPb.
Bool_t GetLeadingParticle(AliCaloTrackParticleCorrelation *particle)
Search Charged or Neutral leading particle, select the highest one and fill AOD.
TH2F * fhNeutralLeadingRatioPt
! Ratio of Pt leading neutral and trigger
Double_t fPtTriggerSelectionCut
Jet pt to change to low pt jets analysis.
Double_t fLeadingRatioMinCut
Leading/gamma Ratio cut minimum.
TH2F * fhJetLeadingRatioPts[5][5]
! Ratio of pt leading and pt jet
TH2F * fhBkgFFzs[5][5]
! Accepted reconstructed background fragmentation function, z=ptjet/pttrig
Double_t fJetXMin1[6]
X Factor to set jet min limit for pp.
TH2F * fhBkgFFxis[5][5]
! Accepted reconstructed background fragmentation function, xsi = ln(pttrig/ptjet) ...
Double_t fLeadingRatioMaxCut
Leading /gamma Ratio cut maximum.
TH2F * fhChargedLeadingPt
! Pt(Pt trigger) distribution of charged hadrons
TLorentzVector fLeadingPhoMom2
! leading particle momentum
TH2F * fhNeutralLeadingPhi
! Phi(Pt trigger) distribution of neutral hadrons
void MakeAnalysisFillHistograms()
Particle-Jet/Hadron Correlation Analysis, fill histograms.
TH2F * fhBkgNTracksInCones[5][5]
! Background multiplicity in cone
Double_t fJetSigma1[2]
Rec. sigma of jet energy parameters.
Bool_t fReMakeJet
Re-make the jet reconstruction from AODParticleCorrelation input.
TH2F * fhBkgPt
! Leading pt bakground vs pt trigger
Double_t fDeltaPhiMinCut
Maximum Delta Phi Gamma-Leading.
TH2F * fhJetFFpts[5][5]
! Jet particle pt distribution in cone
Double_t fJetCTSRatioMaxCut
Jet(CTS) /gamma Ratio cut maximum.
TH2F * fhBkgLeadingRatioPts[5][5]
! Ratio of pt leading and pt background
TH2F * fhNeutralLeadingDeltaPhiRatioPt30
! Difference of neutral hadron and trigger phi as function of pT leading / trigger pT...
TH2F * fhNeutralLeadingDeltaPt
! Difference of neutral hadron and trigger pT as function of trigger pT
TLorentzVector fLeadingChargeMom
! leading particle momentum
TH2F * fhNeutralLeadingXi
! Ln (pt leading neutral / pt trigger)
TH2F * fhBkgNTracksInCone
! Background multiplicity in cone
TH2F * fhChargedLeadingEta
! Eta(Pt trigger) distribution of charged hadrons
void GetLeadingCharge(AliCaloTrackParticleCorrelation *particle)
TH2F * fhNeutralLeadingDeltaEta
! Difference of charged particle and trigger eta as function of trigger pT
void InitParameters()
Initialize the parameters of the analysis.
TH2F * fhBkgLeadingRatioPt
! Ratio of pt leading and pt background
Daughter of AliCaloTrackParticle that includes correlation part.
TH2F * fhNeutralLeadingDeltaPhi
! Difference of neutral hadron and trigger phi as function of trigger pT
TH2F * fhBkgPts[5][5]
! Leading pt bakground vs pt trigger
Base class for CaloTrackCorr analysis algorithms.
TH2F * fhJetFFpt
! Jet particle pt distribution in cone
Bool_t IsJetSelected(Double_t ptTrig, Double_t ptjet) const
TH2F * fhJetFFz
! Accepted reconstructed jet fragmentation function, z=ptjet/pttrig
Int_t fJetNPt
Number of jet particle pT threshold, maximum 5.
TH2F * fhJetLeadingRatioPt
! Ratio of pt leading and pt jet
AliAnaParticleJetLeadingConeCorrelation()
Default constructor. Initialize parameters.
Double_t fJetCone
Jet cone sizes under study (!fSeveralConeAndPtCuts)
void MakeJetFromAOD(AliCaloTrackParticleCorrelation *particle)
int Int_t
Definition: External.C:63
TH2F * fhChargedLeadingPhi
! Phi(Pt trigger) distribution of charged hadrons
Bool_t SelectCluster(AliVCluster *calo, Double_t *vertex, TLorentzVector &mom, Int_t &pdg)
unsigned int UInt_t
Definition: External.C:33
TH2F * fhNeutralLeadingPt
! Pt(Pt trigger) distribution of neutral hadrons
float Float_t
Definition: External.C:68
Bool_t fSelectIsolated
Select only trigger particles isolated.
void FillJetHistos(AliCaloTrackParticleCorrelation *particle, const TLorentzVector jet, const TString &type, const TString &lastname)
Fill jet and background histograms.
Double_t fJetRatioMaxCut
Jet(EMCAL+CTS)/gamma Ratio cut maximum.
TH2F * fhChargedLeadingDeltaPt
! Difference of charged hadron and trigger pT as function of trigger p
TH2F * fhJetFFzs[5][5]
! Accepted reconstructed jet fragmentation function, z=ptjet/pttrig
Correlate trigger particle and reconstructed jet.
TH2F * fhChargedLeadingXi
! Ln (pt leading charge / pt trigger)
Base class for event, clusters and tracks filtering and preparation for the analysis.
AliAnaParticleJetLeadingConeCorrelation & operator=(const AliAnaParticleJetLeadingConeCorrelation &g)
Assignment operator not implemented.
TH2F * fhBkgRatioPt
! Ratio of pt background and pt trigger
TLorentzVector fJetConstMom
! jet constituent momentum
TH2F * fhChargedLeadingRatioPt
! Ratio of Pt leading charge and trigger
TH2F * fhChargedLeadingDeltaPhi
! Difference of charged hadron and trigger phi as function of trigger pT
TLorentzVector fLeadingPi0Mom
! leading particle momentum
TH2F * fhBkgFFxi
! Accepted reconstructed background fragmentation function, xsi = ln(pttrig/ptjet) ...
Double_t fDeltaPhiMaxCut
Minimum Delta Phi Gamma-Leading.
void SetDeltaPhiCutRange(Double_t phimin, Double_t phimax)
Double_t fJetPtThresPbPb
Jet pT threshold under study(!fSeveralConeAndPtCuts)
TH2F * fhJetFFxi
! Accepted reconstructed jet fragmentation function, xsi = ln(pttrig/ptjet)
TLorentzVector fLeadingPhoMom1
! leading particle momentum
TH2F * fhJetFFxis[5][5]
! Accepted reconstructed jet fragmentation function, xsi = ln(pttrig/ptjet)
void Print(const Option_t *opt) const
Print some relevant parameters set for the analysis.
Bool_t IsParticleInJetCone(Double_t eta, Double_t phi, Double_t etal, Double_t phil) const
void GetLeadingPi0(AliCaloTrackParticleCorrelation *particle)
TH2F * fhBkgLeadingDeltaPhis[5][5]
! Delta phi background-leading
TH2F * fhJetRatioPts[5][5]
! Ratio of pt jet and pt trigger
TH2F * fhJetNTracksInCones[5][5]
! jet multiplicity in cone
Double_t fJetCones[5]
Jet cone sizes under study (fSeveralConeAndPtCuts)
TH2F * fhChargedLeadingDeltaPhiRatioPt50
! Difference of charged hadron and trigger phi as function of pT leading / trigger pT...
UInt_t fSelect
kTRUE: Selects all jets, no limits.
Double_t fJetXMax1[6]
X Factor to set jet max limit for pp.
TString fJetNameCones[5]
String name of cone to append to histos.
Double_t CalculateJetRatioLimit(Double_t ptTrig, const Double_t *param, const Double_t *x) const
Bool_t fSeveralConeAndPtCuts
To play with the jet cone size and pt th.
Double_t fJetCTSRatioMinCut
Jet(CTS) /gamma Ratio cut maximum.
Class that contains methods to select candidate cluster pairs to neutral meson.
TH2F * fhNeutralLeadingDeltaPhiRatioPt50
! Difference of neutral hadron and trigger phi as function of pT leading / trigger pT...
const char Option_t
Definition: External.C:48
void SetLeadingRatioCutRange(Double_t ratiomin, Double_t ratiomax)
TH2F * fhBkgRatioPts[5][5]
! Ratio of pt background and pt trigger
TH2F * fhBkgLeadingDeltaEtas[5][5]
! Delta eta background-leading
TH2F * fhNeutralLeadingEta
! Eta(Pt trigger) distribution of neutral hadrons
TH2F * fhJetPts[5][5]
! Leading pt jet vs pt trigger
TString fJetNamePtThres[5]
String name of pt th to append to histos.
void SetJetCTSRatioCutRange(Double_t ratiomin, Double_t ratiomax)
bool Bool_t
Definition: External.C:53
void SetJetRatioCutRange(Double_t ratiomin, Double_t ratiomax)
Double_t fJetRatioMinCut
Jet(EMCAL+CTS)/gamma Ratio cut minimum.
TH2F * fhBkgDeltaPhis[5][5]
! Delta phi background-trigger
TLorentzVector fLeadingMom
! leading particle momentum
TH2F * fhBkgDeltaEtas[5][5]
! Delta eta background-trigger
TH2F * fhBkgFFpt
! Background particle pt distribution in cone
TH2F * fhBkgFFpts[5][5]
! Background particle pt distribution in cone
const Double_t phimin
TH2F * fhChargedLeadingDeltaEta
! Difference of charged particle and trigger eta as function of trigger pT