]>
Commit | Line | Data |
---|---|---|
96c271c1 | 1 | /************************************************************************* |
2 | * * | |
bce887f5 | 3 | * Task for Jet Chemistry Analysis in PWG-JE Jet Task Force Train * |
96c271c1 | 4 | * * |
5 | * * | |
bce887f5 | 6 | * contact: * |
7 | * Alice Zimmermann * | |
8 | * zimmermann@physi.uni-heidelberg.de * | |
9 | * * | |
10 | * * | |
11 | * * | |
96c271c1 | 12 | * * |
13 | *************************************************************************/ | |
bce887f5 | 14 | |
9b2de807 | 15 | |
16 | #ifndef ALIANALYSISTASKJETCHEM_H | |
17 | #define ALIANALYSISTASKJETCHEM_H | |
18 | ||
bce887f5 | 19 | /* Copyright(c) 1998-2013, ALICE Experiment at CERN, All rights reserved. * |
9b2de807 | 20 | * See cxx source for full Copyright notice */ |
21 | ||
96c271c1 | 22 | /* $Id$ */ |
9b2de807 | 23 | |
bce887f5 | 24 | class AliAODv0; |
25 | class AliAODVertex; | |
26 | class AliAnalysisCentralitySelector; | |
27 | class AliPIDResponse; | |
28 | class TString; | |
29 | class TList; | |
30 | class AliAODMCParticle; | |
31 | class AliAODTrack; | |
32 | ||
96c271c1 | 33 | #include "AliAnalysisTaskFragmentationFunction.h" |
8c483a6b | 34 | #include "AliPID.h" |
bce887f5 | 35 | #include "AliAnalysisTaskSE.h" |
36 | #include "AliAODMCHeader.h" | |
9b2de807 | 37 | |
96c271c1 | 38 | class AliAnalysisTaskJetChem : public AliAnalysisTaskFragmentationFunction { |
9b2de807 | 39 | |
96c271c1 | 40 | public: |
41 | ||
42 | //---------------------------------------- | |
43 | class AliFragFuncHistosInvMass : public TObject | |
44 | { | |
9b2de807 | 45 | |
96c271c1 | 46 | public: |
9b2de807 | 47 | |
96c271c1 | 48 | AliFragFuncHistosInvMass(const char* name = "FFIMhistos", |
49 | Int_t nJetPt = 0, Float_t jetPtMin = 0, Float_t jetPtMax = 0, | |
50 | Int_t nInvMass = 0, Float_t invMassMin=0, Float_t invMassMax=0, | |
51 | Int_t nPt = 0, Float_t ptMin = 0, Float_t ptMax = 0, | |
52 | Int_t nXi = 0, Float_t xiMin = 0, Float_t xiMax = 0, | |
53 | Int_t nZ = 0, Float_t zMin = 0, Float_t zMax = 0); | |
54 | AliFragFuncHistosInvMass(const AliFragFuncHistosInvMass& copy); | |
55 | AliFragFuncHistosInvMass& operator=(const AliFragFuncHistosInvMass &o); | |
56 | virtual ~AliFragFuncHistosInvMass(); | |
9b2de807 | 57 | |
96c271c1 | 58 | virtual void DefineHistos(); |
59 | virtual void FillFF(Float_t trackPt, Float_t invM, Float_t jetPt,Bool_t incrementJetPt); | |
60 | virtual void AddToOutput(TList* list) const; | |
9b2de807 | 61 | |
96c271c1 | 62 | private: |
9b2de807 | 63 | |
96c271c1 | 64 | Int_t fNBinsJetPt; // FF histos bins |
65 | Float_t fJetPtMin; // FF histos limits | |
66 | Float_t fJetPtMax; // FF histos limits | |
67 | Int_t fNBinsInvMass; // FF histos bins | |
68 | Float_t fInvMassMin; // FF histos limits | |
69 | Float_t fInvMassMax; // FF histos limits | |
70 | Int_t fNBinsPt; // FF histos bins | |
71 | Float_t fPtMin; // FF histos limits | |
72 | Float_t fPtMax; // FF histos limits | |
73 | Int_t fNBinsXi; // FF histos bins | |
74 | Float_t fXiMin; // FF histos limits | |
75 | Float_t fXiMax; // FF histos limits | |
76 | Int_t fNBinsZ; // FF histos bins | |
77 | Float_t fZMin; // FF histos limits | |
78 | Float_t fZMax; // FF histos limits | |
79 | ||
bce887f5 | 80 | |
81 | ||
96c271c1 | 82 | TH3F* fh3TrackPt; //! FF: track transverse momentum |
83 | TH3F* fh3Xi; //! FF: xi | |
84 | TH3F* fh3Z; //! FF: z | |
bce887f5 | 85 | TH1F* fh1JetPt; //! jet pt of all jets |
9b2de807 | 86 | |
96c271c1 | 87 | TString fNameFF; // histo names prefix |
88 | ||
89 | ClassDef(AliFragFuncHistosInvMass, 1); | |
90 | }; | |
91 | ||
96c271c1 | 92 | |
93 | //---------------------------------------- | |
9b2de807 | 94 | |
96c271c1 | 95 | AliAnalysisTaskJetChem(); |
96 | AliAnalysisTaskJetChem(const char *name); | |
97 | AliAnalysisTaskJetChem(const AliAnalysisTaskJetChem ©); | |
98 | AliAnalysisTaskJetChem& operator=(const AliAnalysisTaskJetChem &o); | |
99 | virtual ~AliAnalysisTaskJetChem(); | |
100 | ||
101 | virtual void UserCreateOutputObjects(); | |
102 | virtual void UserExec(Option_t *option); | |
103 | ||
6ec586c1 | 104 | enum { kTrackUndef =0, kOnFly, kOnFlyPID, kOnFlydEdx, kOnFlyPrim, kOffl, kOfflPID, kOffldEdx, kOfflPrim }; |
105 | enum { kK0, kLambda, kAntiLambda }; | |
106 | ||
96c271c1 | 107 | static void SetProperties(TH3F* h,const char* x, const char* y,const char* z); |
108 | ||
b6eb08db | 109 | Bool_t IsAccepteddEdx(Double_t mom, Double_t signal, AliPID::EParticleType n, Double_t cutnSig) const;//not used anymore |
110 | Bool_t IsK0InvMass(Double_t mass) const; | |
111 | Int_t GetListOfV0s(TList *list, Int_t type, Int_t particletype, AliAODVertex* primVertex, AliAODEvent* aod); | |
112 | Int_t GetListOfParticles(TList *list, Int_t type, Int_t particletype, AliAODVertex* primVertex); | |
113 | Int_t GetListOfMCParticles(TList *outputlist, Int_t particletype, AliAODEvent* mcaodevent); | |
6ec586c1 | 114 | void GetTracksInCone(TList* inputlist, TList* outputlist, const AliAODJet* jet, Double_t radius, Double_t& sumPt, Double_t minPt, Double_t maxPt, Bool_t& isBadPt); |
b6eb08db | 115 | void GetTracksInPerpCone(TList* inputlist, TList* outputlist, const AliAODJet* jet, Double_t radius, Double_t& sumPerpPt); |
116 | Bool_t MCLabelCheck(AliAODv0* v0, Int_t particletype, const AliAODTrack* trackNeg, const AliAODTrack* trackPos, TList *listmc, Int_t& negDaughterpdg, Int_t& posDaughterpdg, Int_t& motherType, Int_t& v0Label, Double_t& MCPt, Bool_t& fPhysicalPrimary, Int_t& MCv0PDGCode); | |
117 | Bool_t IsParticleMatching(const AliAODMCParticle* mcp0, Int_t v0Label); | |
bce887f5 | 118 | Bool_t DaughterTrackCheck(AliAODv0* v0, Int_t& nnum, Int_t& pnum); |
119 | Int_t IsTrackInjected(AliAODv0 *v0, AliAODMCHeader *header, TClonesArray *arrayMC); | |
120 | TString GetGenerator(Int_t label, AliAODMCHeader* header); | |
121 | Double_t SmearJetPt(Double_t jetPt, Int_t cl, Double_t jetRadius, Double_t ptmintrack, Double_t& jetPtSmear); | |
122 | AliAODJet* GetMedianCluster(); | |
123 | ||
96c271c1 | 124 | virtual void SetK0Type(Int_t i){ fK0Type = i; } |
125 | virtual void SetFilterMaskK0(UInt_t i) {fFilterMaskK0 = i;} | |
126 | ||
b6eb08db | 127 | Bool_t IsLaInvMass(Double_t mass) const; //La and ALa mass check |
bce887f5 | 128 | virtual void SetLaType(Int_t i){ fLaType = i; } |
129 | virtual void SetFilterMaskLa(UInt_t i) {fFilterMaskLa = i;} | |
130 | ||
131 | virtual void SetALaType(Int_t i){ fALaType = i; } | |
132 | virtual void SetFilterMaskALa(UInt_t i) {fFilterMaskALa = i;} | |
133 | ||
134 | virtual void SetSelectArmenteros(Bool_t b) {IsArmenterosSelected = b;} | |
135 | //virtual void SetEventSelectionMask(UInt_t i){fEvtSelectionMask = i;} //already inherited by AliAnalysisFragmentationFunctionTask | |
136 | //virtual void UsePhysicsSelection(Bool_t b) {fUsePhysicsSelection = b;} //already inherited by AliAnalysisFragmentationFunctionTask | |
96c271c1 | 137 | |
b6eb08db | 138 | void CalculateInvMass(AliAODv0* v0vtx, Int_t particletype, Double_t& invM, Double_t& trackPt); |
bce887f5 | 139 | |
b6eb08db | 140 | Bool_t AcceptBetheBloch(AliAODv0 *v0, AliPIDResponse *PIDResponse, Int_t particletype); //don't use this method for MC Analysis |
bce887f5 | 141 | |
142 | ||
143 | Double_t MyRapidity(Double_t rE, Double_t rPz) const; | |
144 | ||
145 | //-- K0s | |
146 | ||
119061b1 | 147 | void SetFFInvMassHistoBins(Int_t nJetPt = 19, Float_t jetPtMin = 5., Float_t jetPtMax = 100., //previous 19, 5.,100. |
dec584ae | 148 | Int_t nInvM = 400, Float_t invMMin = 0.3, Float_t invMMax = 0.7, //previous 0.4 to 0.6 |
bce887f5 | 149 | Int_t nPt = 200, Float_t ptMin = 0., Float_t ptMax = 20., //previous 0. to 10. |
96c271c1 | 150 | Int_t nXi = 35, Float_t xiMin = 0., Float_t xiMax = 7., |
151 | Int_t nZ = 11, Float_t zMin = 0., Float_t zMax = 1.1) | |
152 | { fFFIMNBinsJetPt = nJetPt; fFFIMJetPtMin = jetPtMin; fFFIMJetPtMax = jetPtMax; | |
153 | fFFIMNBinsInvM = nInvM; fFFIMInvMMin = invMMin; fFFIMInvMMax = invMMax; fFFIMNBinsPt = nPt; fFFIMPtMin = ptMin; fFFIMPtMax = ptMax; | |
154 | fFFIMNBinsXi = nXi; fFFIMXiMin = xiMin; fFFIMXiMax = xiMax; fFFIMNBinsZ = nZ; fFFIMZMin = zMin; fFFIMZMax = zMax; } | |
155 | ||
bce887f5 | 156 | |
96c271c1 | 157 | |
bce887f5 | 158 | // --La and ALa |
159 | ||
119061b1 | 160 | void SetFFInvMassLaHistoBins(Int_t nJetPt = 19, Float_t jetPtMin = 5., Float_t jetPtMax = 100., //La |
bce887f5 | 161 | //Int_t nInvM = 140, Float_t invMMin = 1.06, Float_t invMMax = 1.2,//original inv. mass range, now I shifted to Vits slightly larger mass window |
162 | Int_t nInvM = 200, Float_t invMMin = 1.05, Float_t invMMax = 1.25, | |
163 | Int_t nPt = 200, Float_t ptMin = 0., Float_t ptMax = 20., | |
164 | Int_t nXi = 35, Float_t xiMin = 0., Float_t xiMax = 7., | |
165 | Int_t nZ = 11, Float_t zMin = 0., Float_t zMax = 1.1) | |
166 | ||
167 | { fFFIMLaNBinsJetPt = nJetPt; fFFIMLaJetPtMin = jetPtMin; fFFIMLaJetPtMax = jetPtMax; | |
168 | fFFIMLaNBinsInvM = nInvM; fFFIMLaInvMMin = invMMin; fFFIMLaInvMMax = invMMax; fFFIMLaNBinsPt = nPt; fFFIMLaPtMin = ptMin; fFFIMLaPtMax = ptMax; | |
169 | fFFIMLaNBinsXi = nXi; fFFIMLaXiMin = xiMin; fFFIMLaXiMax = xiMax; fFFIMLaNBinsZ = nZ; fFFIMLaZMin = zMin; fFFIMLaZMax = zMax; } | |
170 | ||
171 | ||
bce887f5 | 172 | |
96c271c1 | 173 | // consts |
bce887f5 | 174 | |
bce887f5 | 175 | |
176 | //-- | |
177 | Bool_t fAnalysisMC; | |
178 | Double_t fDeltaVertexZ; | |
a6ffe507 | 179 | Double_t fCutjetEta; |
bce887f5 | 180 | Double_t fCuttrackNegNcls; |
181 | Double_t fCuttrackPosNcls; | |
182 | Double_t fCutPostrackRap; | |
183 | Double_t fCutNegtrackRap; | |
184 | Double_t fCutRap; | |
185 | Double_t fCutPostrackEta; | |
186 | Double_t fCutNegtrackEta; | |
187 | Double_t fCutEta; | |
188 | Double_t fCutV0cosPointAngle; | |
bce887f5 | 189 | Bool_t fKinkDaughters; |
190 | Bool_t fRequireTPCRefit; | |
191 | Double_t fCutArmenteros; | |
192 | Double_t fCutV0DecayMin; | |
193 | Double_t fCutV0DecayMax; | |
194 | Double_t fCutV0totMom; | |
195 | Double_t fCutDcaV0Daughters; | |
196 | Double_t fCutDcaPosToPrimVertex; | |
197 | Double_t fCutDcaNegToPrimVertex; | |
198 | Double_t fCutV0RadiusMin; | |
199 | Double_t fCutV0RadiusMax; | |
200 | Double_t fCutBetheBloch; | |
201 | Double_t fCutRatio; | |
202 | ||
203 | // cuts | |
204 | void SetCuttrackPosNcls(Double_t posNcls){fCuttrackPosNcls=posNcls; Printf("AliAnalysisTaskJetChem:: SetCuttrackPosNcls %f",posNcls);} | |
205 | void SetCuttrackNegNcls(Double_t negNcls){fCuttrackNegNcls=negNcls; Printf("AliAnalysisTaskJetChem:: SetCuttrackNegNcls %f",negNcls);} | |
206 | void SetCuttrackPosRap(Double_t posRap){fCutPostrackRap=posRap; Printf("AliAnalysisTaskJetChem:: SetCuttrackPosRap %f",posRap);} | |
207 | void SetCuttrackNegRap(Double_t negRap){fCutNegtrackRap=negRap; Printf("AliAnalysisTaskJetChem:: SetCuttrackNegRap %f",negRap);} | |
208 | void SetCutV0Rap(Double_t v0Rap){fCutRap=v0Rap; Printf("AliAnalysisTaskJetChem:: SetCutV0Rap %f",v0Rap);} | |
209 | void SetCuttrackPosEta(Double_t posEta){fCutPostrackEta=posEta; Printf("AliAnalysisTaskJetChem:: SetCuttrackPosEta %f",posEta);} | |
210 | void SetCuttrackNegEta(Double_t negEta){fCutNegtrackEta=negEta; Printf("AliAnalysisTaskJetChem:: SetCuttrackNegEta %f",negEta);} | |
211 | void SetCutV0Eta(Double_t v0Eta){fCutEta=v0Eta; Printf("AliAnalysisTaskJetChem:: SetCutV0Eta %f",v0Eta);} | |
a6ffe507 | 212 | void SetCutJetEta(Double_t jetEta){fCutjetEta=jetEta; Printf("AliAnalysisTaskJetChem:: SetCutjetEta %f",jetEta);} |
bce887f5 | 213 | void SetCosOfPointingAngle(Double_t cospointAng){fCutV0cosPointAngle=cospointAng; Printf("AliAnalysisTaskJetChem:: SetCosOfPointingAngle %f",cospointAng);} |
bce887f5 | 214 | void SetAcceptKinkDaughters(Bool_t isKinkDaughtersAccepted){fKinkDaughters=isKinkDaughtersAccepted; Printf("AliAnalysisTaskJetChem:: SetAcceptKinkDaughters %i", isKinkDaughtersAccepted);} |
215 | void SetRequireTPCRefit(Bool_t isTPCRefit){fRequireTPCRefit=isTPCRefit; Printf("AliAnalysisTaskJetChem:: SetRequireTPCRefit %i", isTPCRefit);} | |
216 | void SetCutArmenteros(Double_t armenteros){fCutArmenteros=armenteros; Printf("AliAnalysisTaskJetChem:: SetCutArmenteros %f", armenteros);} | |
217 | void SetCutV0DecayMin(Double_t decayMin){fCutV0DecayMin=decayMin; Printf("AliAnalysisTaskJetChem:: SetCutDecayMin %f", decayMin);} | |
218 | void SetCutV0DecayMax(Double_t decayMax){fCutV0DecayMax=decayMax; Printf("AliAnalysisTaskJetChem:: SetCutDecayMax %f", decayMax);} | |
219 | void SetCutV0totMom(Double_t v0totMom){fCutV0totMom=v0totMom; Printf("AliAnalysisTaskJetChem:: SetCutV0totMom %f", v0totMom);} | |
220 | void SetCutDcaV0Daughters(Double_t dcav0daughters){fCutDcaV0Daughters=dcav0daughters; Printf("AliAnalysisTaskJetChem:: SetCutDcaV0Daughters %f", dcav0daughters);} | |
221 | void SetCutDcaPosToPrimVertex(Double_t dcaPosToPrimVertex){fCutDcaPosToPrimVertex=dcaPosToPrimVertex; Printf("AliAnalysisTaskJetChem:: SetCutDcaPosToPrimVertex %f", dcaPosToPrimVertex);} | |
222 | void SetCutDcaNegToPrimVertex(Double_t dcaNegToPrimVertex){fCutDcaNegToPrimVertex=dcaNegToPrimVertex; Printf("AliAnalysisTaskJetChem:: SetCutDcaNegToPrimVertex %f", dcaNegToPrimVertex);} | |
223 | void SetCutV0RadiusMin(Double_t v0RadiusMin){fCutV0RadiusMin=v0RadiusMin; Printf("AliAnalysisTaskJetChem:: SetCutV0RadiusMin %f", v0RadiusMin);} | |
224 | void SetCutV0RadiusMax(Double_t v0RadiusMax){fCutV0RadiusMax=v0RadiusMax; Printf("AliAnalysisTaskJetChem:: SetCutV0RadiusMax %f", v0RadiusMax);} | |
225 | void SetCutBetheBloch(Double_t cutBetheBloch){fCutBetheBloch=cutBetheBloch; Printf("AliAnalysisTaskJetChem:: SetCutBetheBloch %f", cutBetheBloch);} | |
226 | void SetCutRatioTPC(Double_t cutRatioTPC){fCutRatio=cutRatioTPC; Printf("AliAnalysisTaskJetChem:: SetCutRatioTPC %f", cutRatioTPC);} | |
227 | void SetAnalysisMC(Bool_t analysisMC) {fAnalysisMC = analysisMC;} | |
228 | void SetDeltaZVertexCut(Float_t deltaVtxZ){fDeltaVertexZ = deltaVtxZ;} | |
229 | ||
96c271c1 | 230 | private: |
231 | ||
a6ffe507 | 232 | Int_t fK0Type; // K0 cuts |
96c271c1 | 233 | UInt_t fFilterMaskK0; //! K0 legs cuts |
234 | TList* fListK0s; //! K0 list | |
bce887f5 | 235 | AliPIDResponse *fPIDResponse; // PID |
96c271c1 | 236 | |
237 | AliFragFuncQATrackHistos* fV0QAK0; //! track QA: V0s in K0 inv mass range | |
238 | AliFragFuncHistos* fFFHistosRecCutsK0Evt; //! inclusive FF for K0 evt | |
239 | AliFragFuncHistosInvMass* fFFHistosIMK0AllEvt; //! K0 pt spec for all events | |
240 | AliFragFuncHistosInvMass* fFFHistosIMK0Jet; //! K0 FF all dPhi | |
241 | AliFragFuncHistosInvMass* fFFHistosIMK0Cone; //! K0 FF jet cone | |
bce887f5 | 242 | |
a6ffe507 | 243 | Int_t fLaType; // La cuts |
bce887f5 | 244 | UInt_t fFilterMaskLa; //! La legs cuts |
245 | TList* fListLa; //! La list | |
246 | ||
247 | AliFragFuncHistosInvMass* fFFHistosIMLaAllEvt; //! La pt spec for all events | |
248 | AliFragFuncHistosInvMass* fFFHistosIMLaJet; //! La FF all dPhi | |
249 | AliFragFuncHistosInvMass* fFFHistosIMLaCone; //! La FF jet cone | |
bce887f5 | 250 | |
a6ffe507 | 251 | Int_t fALaType; // ALa cuts |
bce887f5 | 252 | |
253 | UInt_t fFilterMaskALa; //! ALa legs cuts | |
254 | TList* fListALa; //! ALa list | |
255 | TList* fListFeeddownLaCand; //! feeddown from Xi (-,0) | |
256 | TList* fListFeeddownALaCand; //! feeddown from Xibar (+,0) | |
257 | TList* jetConeFDLalist; //! feeddown from Xi (-,0) in jet cone | |
258 | TList* jetConeFDALalist; //! feeddown from Xibar (+,0) in jet cone | |
259 | TList* fListMCgenK0s; //! MC generated K0s | |
260 | TList* fListMCgenLa; //! MC generated La | |
261 | TList* fListMCgenALa; //! MC generated ALa | |
262 | TList* fListMCgenK0sCone; //! MC generated K0s in cone around jet axis, particles are from fragmentation but also from underlying event | |
263 | TList* fListMCgenLaCone; //! MC generated Lambdas in cone around jet axis, particles are from fragmentation but also from underlying event | |
264 | TList* fListMCgenALaCone; //! MC generated Antilambdas in cone around jet axis, particles are from fragmentation but also from underlying event | |
265 | ||
266 | Bool_t IsArmenterosSelected; //Armenteros-Podolanski Cut (is/isn't) applied | |
267 | ||
268 | AliFragFuncHistosInvMass* fFFHistosIMALaAllEvt; //! ALa pt spec for all events | |
269 | AliFragFuncHistosInvMass* fFFHistosIMALaJet; //! ALa FF all dPhi | |
270 | AliFragFuncHistosInvMass* fFFHistosIMALaCone; //! ALa FF jet cone | |
bce887f5 | 271 | |
272 | // histogram bins | |
273 | ||
96c271c1 | 274 | |
96c271c1 | 275 | |
bce887f5 | 276 | //--K0s |
277 | ||
96c271c1 | 278 | Int_t fFFIMNBinsJetPt; // FF histos bins |
279 | Float_t fFFIMJetPtMin; // FF histos limits | |
280 | Float_t fFFIMJetPtMax; // FF histos limits | |
bce887f5 | 281 | |
96c271c1 | 282 | Int_t fFFIMNBinsInvM; // FF histos bins |
283 | Float_t fFFIMInvMMin; // FF histos bins | |
284 | Float_t fFFIMInvMMax; // FF histos bins | |
bce887f5 | 285 | |
96c271c1 | 286 | Int_t fFFIMNBinsPt; // FF histos bins |
287 | Float_t fFFIMPtMin; // FF histos limits | |
288 | Float_t fFFIMPtMax; // FF histos limits | |
bce887f5 | 289 | |
96c271c1 | 290 | Int_t fFFIMNBinsXi; // FF histos bins |
291 | Float_t fFFIMXiMin; // FF histos limits | |
292 | Float_t fFFIMXiMax; // FF histos limits | |
bce887f5 | 293 | |
96c271c1 | 294 | Int_t fFFIMNBinsZ; // FF histos bins |
295 | Float_t fFFIMZMin; // FF histos limits | |
296 | Float_t fFFIMZMax; // FF histos limits | |
bce887f5 | 297 | |
298 | //--La | |
299 | ||
300 | Int_t fFFIMLaNBinsJetPt; // FF histos bins | |
301 | Float_t fFFIMLaJetPtMin; // FF histos limits | |
302 | Float_t fFFIMLaJetPtMax; // FF histos limits | |
303 | ||
304 | Int_t fFFIMLaNBinsInvM; // FF histos bins | |
305 | Float_t fFFIMLaInvMMin; // FF histos bins | |
306 | Float_t fFFIMLaInvMMax; // FF histos bins | |
307 | ||
308 | Int_t fFFIMLaNBinsPt; // FF histos bins | |
309 | Float_t fFFIMLaPtMin; // FF histos limits | |
310 | Float_t fFFIMLaPtMax; // FF histos limits | |
311 | ||
312 | Int_t fFFIMLaNBinsXi; // FF histos bins | |
313 | Float_t fFFIMLaXiMin; // FF histos limits | |
314 | Float_t fFFIMLaXiMax; // FF histos limits | |
315 | ||
316 | Int_t fFFIMLaNBinsZ; // FF histos bins | |
317 | Float_t fFFIMLaZMin; // FF histos limits | |
318 | Float_t fFFIMLaZMax; // FF histos limits | |
319 | ||
96c271c1 | 320 | |
d7fb73a4 | 321 | |
96c271c1 | 322 | // Histograms |
bce887f5 | 323 | |
324 | TH1F* fh1EvtAllCent; | |
325 | TH1F* fh1Evt; | |
326 | TH1F* fh1K0Mult; | |
327 | TH1F* fh1dPhiJetK0; | |
328 | TH1F* fh1LaMult; | |
329 | TH1F* fh1dPhiJetLa; | |
330 | TH1F* fh1ALaMult; | |
331 | TH1F* fh1dPhiJetALa; | |
332 | TH1F* fh1JetEta; | |
333 | TH1F* fh1JetPhi; | |
334 | TH2F* fh2JetEtaPhi; | |
335 | TH1F* fh1V0JetPt; | |
336 | TH2F* fh2FFJetTrackEta; //charged jet track eta distribution | |
337 | TH1F* fh1trackPosNCls; | |
338 | TH1F* fh1trackNegNCls; | |
339 | TH1F* fh1trackPosRap; | |
340 | TH1F* fh1trackNegRap; | |
341 | TH1F* fh1V0Rap; | |
342 | TH1F* fh1trackPosEta; | |
343 | TH1F* fh1trackNegEta; | |
344 | TH1F* fh1V0Eta; | |
345 | TH1F* fh1V0totMom; | |
346 | TH1F* fh1CosPointAngle; | |
bce887f5 | 347 | TH1F* fh1DecayLengthV0; |
348 | TH2F* fh2ProperLifetimeK0sVsPtBeforeCut; | |
349 | TH2F* fh2ProperLifetimeK0sVsPtAfterCut; | |
bce887f5 | 350 | TH1F* fh1V0Radius; |
351 | TH1F* fh1DcaV0Daughters; | |
352 | TH1F* fh1DcaPosToPrimVertex; | |
353 | TH1F* fh1DcaNegToPrimVertex; | |
354 | TH2F* fh2ArmenterosBeforeCuts; | |
355 | TH2F* fh2ArmenterosAfterCuts; | |
bce887f5 | 356 | TH2F* fh2BBLaPos; |
357 | TH2F* fh2BBLaNeg; | |
bce887f5 | 358 | TH1F* fh1PosDaughterCharge; |
359 | TH1F* fh1NegDaughterCharge; | |
360 | TH1F* fh1PtMCK0s; | |
361 | TH1F* fh1PtMCLa; | |
362 | TH1F* fh1PtMCALa; | |
363 | TH1F* fh1EtaK0s; | |
364 | TH1F* fh1EtaLa; | |
365 | TH1F* fh1EtaALa; | |
366 | TH3F* fh3InvMassEtaTrackPtK0s; | |
367 | TH3F* fh3InvMassEtaTrackPtLa; | |
368 | TH3F* fh3InvMassEtaTrackPtALa; | |
bce887f5 | 369 | TH1F* fh1TrackMultCone; |
370 | TH2F* fh2TrackMultCone; | |
8bde4137 | 371 | TH2F* fh2NJK0; |
372 | TH2F* fh2NJLa; | |
373 | TH2F* fh2NJALa; | |
bce887f5 | 374 | TH2F* fh2MCgenK0Cone; |
375 | TH2F* fh2MCgenLaCone; | |
376 | TH2F* fh2MCgenALaCone; | |
377 | TH2F* fh2MCEtagenK0Cone; | |
378 | TH2F* fh2MCEtagenLaCone; | |
379 | TH2F* fh2MCEtagenALaCone; | |
380 | TH1F* fh1FFIMK0ConeSmear; | |
381 | TH1F* fh1FFIMLaConeSmear; | |
382 | TH1F* fh1FFIMALaConeSmear; | |
383 | TH3F* fh3MCrecK0Cone; | |
384 | TH3F* fh3MCrecLaCone; | |
385 | TH3F* fh3MCrecALaCone; | |
386 | TH3F* fh3MCrecK0ConeSmear; | |
387 | TH3F* fh3MCrecLaConeSmear; | |
388 | TH3F* fh3MCrecALaConeSmear; | |
389 | TH3F* fh3SecContinCone; | |
390 | TH3F* fh3StrContinCone; | |
3ed9f288 | 391 | THnSparse* fhnK0sIncl; |
392 | THnSparse* fhnK0sCone; | |
393 | THnSparse* fhnLaIncl; | |
394 | THnSparse* fhnLaCone; | |
395 | THnSparse* fhnALaIncl; | |
396 | THnSparse* fhnALaCone; | |
bce887f5 | 397 | TH3F* fh3IMK0PerpCone; |
398 | TH3F* fh3IMLaPerpCone; | |
399 | TH3F* fh3IMALaPerpCone; | |
400 | TH3F* fh3IMK0MedianCone; | |
401 | TH3F* fh3IMLaMedianCone; | |
402 | TH3F* fh3IMALaMedianCone; | |
a6ffe507 | 403 | TH1F* fh1MedianEta; |
dec584ae | 404 | TH1F* fh1JetPtMedian; //for normalisation by total number of median cluster jets TH3F* fh3IMALaMedianCone; |
bce887f5 | 405 | TH1F* fh1MCMultiplicityPrimary; |
406 | TH1F* fh1MCMultiplicityTracks; | |
bce887f5 | 407 | TH3F* fh3FeedDownLa; |
6ec586c1 | 408 | TH3F* fh3FeedDownALa; |
dec584ae | 409 | TH3F* fh3FeedDownLaCone; |
6ec586c1 | 410 | TH3F* fh3FeedDownALaCone; |
bce887f5 | 411 | TH1F* fh1MCProdRadiusK0s; |
412 | TH1F* fh1MCProdRadiusLambda; | |
413 | TH1F* fh1MCProdRadiusAntiLambda; | |
414 | TH1F* fh1MCPtV0s; | |
415 | TH1F* fh1MCPtK0s; | |
416 | TH1F* fh1MCPtLambda; | |
417 | TH1F* fh1MCPtAntiLambda; | |
418 | TH1F* fh1MCXiPt; | |
419 | TH1F* fh1MCXibarPt; | |
420 | TH2F* fh2MCEtaVsPtK0s; | |
421 | TH2F* fh2MCEtaVsPtLa; | |
422 | TH2F* fh2MCEtaVsPtALa; | |
423 | TH1F* fh1MCRapK0s; | |
424 | TH1F* fh1MCRapLambda; | |
425 | TH1F* fh1MCRapAntiLambda; | |
426 | TH1F* fh1MCEtaAllK0s; | |
427 | TH1F* fh1MCEtaK0s; | |
428 | TH1F* fh1MCEtaLambda; | |
429 | TH1F* fh1MCEtaAntiLambda; | |
9b2de807 | 430 | |
9b2de807 | 431 | |
e2c14c6e | 432 | |
96c271c1 | 433 | ClassDef(AliAnalysisTaskJetChem, 3); |
434 | }; | |
9b2de807 | 435 | |
436 | #endif | |
bce887f5 | 437 | |
438 |