1 /**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
7 * Permission to use, copy, modify and distribute this software and its *
8 * documentation strictly for non-commercial purposes is hereby granted *
9 * without fee, provided that the above copyright notice appears in all *
10 * copies and that both the copyright notice and this permission notice *
11 * appear in the supporting documentation. The authors make no claims *
12 * about the suitability of this software for any purpose. It is *
13 * provided "as is" without express or implied warranty. *
14 **************************************************************************/
16 // +-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
18 // Modified version of AliAnalysisTaskCheckCascade.h
19 // Used bits of code from AliAnalysisTaskCheckPerformanceStrange
21 // --- David Dobrigkeit Chinellato
23 // +-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
25 #ifndef ALIANALYSISTASKEXTRACTPERFORMANCEV0_H
26 #define ALIANALYSISTASKEXTRACTPERFORMANCEV0_H
36 class AliESDtrackCuts;
38 class AliPhysicsSelection;
41 //#include "TString.h"
42 //#include "AliESDtrackCuts.h"
43 #include "AliAnalysisTaskSE.h"
45 class AliAnalysisTaskExtractPerformanceV0 : public AliAnalysisTaskSE {
47 AliAnalysisTaskExtractPerformanceV0();
48 AliAnalysisTaskExtractPerformanceV0(const char *name);
49 virtual ~AliAnalysisTaskExtractPerformanceV0();
51 virtual void UserCreateOutputObjects();
52 virtual void UserExec(Option_t *option);
53 virtual void Terminate(Option_t *);
54 Double_t MyRapidity(Double_t rE, Double_t rPz) const;
55 void CheckChargeV0(AliESDv0 *thisv0);
57 void SetIsNuclear (Bool_t lIsNuclear = kTRUE ) { fkIsNuclear = lIsNuclear; }
58 void SetIsLowEnergyPP (Bool_t lLowEnergyPP = kTRUE ) { fkLowEnergyPP = lLowEnergyPP; }
59 void SetUseOnTheFly (Bool_t lUseOnTheFly = kTRUE ) { fkUseOnTheFly = lUseOnTheFly; }
62 // Note : In ROOT, "//!" means "do not stream the data from Master node to Worker node" ...
63 // your data member object is created on the worker nodes and streaming is not needed.
64 // http://root.cern.ch/download/doc/11InputOutput.pdf, page 14
65 TList *fListHistV0; //! List of Cascade histograms
66 TTree *fTree; //! Output Tree, V0
68 AliPIDResponse *fPIDResponse; // PID response object
69 AliESDtrackCuts *fESDtrackCuts; // ESD track cuts used for primary track definition
71 //Objects Controlling Task Behaviour
73 Bool_t fkIsNuclear; //if true, replace multiplicity est. by centrality (default FALSE)
74 Bool_t fkLowEnergyPP; //if true, skip FASTOnly (default FALSE)
75 Bool_t fkUseOnTheFly; //if true, will use On-the-fly V0s instead of Offline V0s (default FALSE)
77 //===========================================================================================
78 // Variables for Tree, V0s
79 //===========================================================================================
80 Int_t fTreeVariablePrimaryStatus; //!
81 Int_t fTreeVariablePrimaryStatusMother; //!
82 Float_t fTreeVariableChi2V0; //!
83 Float_t fTreeVariableDcaV0Daughters; //!
84 Float_t fTreeVariableDcaV0ToPrimVertex; //!
85 Float_t fTreeVariableDcaPosToPrimVertex; //!
86 Float_t fTreeVariableDcaNegToPrimVertex; //!
87 Float_t fTreeVariableV0CosineOfPointingAngle; //!
88 Float_t fTreeVariableV0Radius; //!
89 Float_t fTreeVariablePt; //!
90 Float_t fTreeVariablePtMC; //!
91 Float_t fTreeVariableRapK0Short; //!
92 Float_t fTreeVariableRapLambda; //!
93 Float_t fTreeVariableRapMC; //!
94 Float_t fTreeVariableInvMassK0s; //!
95 Float_t fTreeVariableInvMassLambda; //!
96 Float_t fTreeVariableInvMassAntiLambda; //!
97 Float_t fTreeVariableAlphaV0; //!
98 Float_t fTreeVariablePtArmV0;//!
99 Float_t fTreeVariableNegTotMomentum; //!
100 Float_t fTreeVariablePosTotMomentum; //!
101 Float_t fTreeVariableNegTransvMomentum; //!
102 Float_t fTreeVariablePosTransvMomentum; //!
103 Float_t fTreeVariableNegTransvMomentumMC; //!
104 Float_t fTreeVariablePosTransvMomentumMC; //!
106 Float_t fTreeVariableNSigmasPosProton; //!
107 Float_t fTreeVariableNSigmasPosPion; //!
108 Float_t fTreeVariableNSigmasNegProton; //!
109 Float_t fTreeVariableNSigmasNegPion; //!
111 Float_t fTreeVariablePtMother; //!
112 Float_t fTreeVariableV0CreationRadius; //!
113 Int_t fTreeVariablePID; //!
114 Int_t fTreeVariablePIDPositive; //!
115 Int_t fTreeVariablePIDNegative; //!
116 Int_t fTreeVariablePIDMother; //!
117 Int_t fTreeVariableIndexStatus; //!
118 Int_t fTreeVariableIndexStatusMother; //!
120 Int_t fTreeVariableRunNumber; //!
121 ULong64_t fTreeVariableEventNumber; //!
123 //Note: TDistOverTotMom needs a mass hypothesis to be converted to proper decaylength.
124 Float_t fTreeVariableDistOverTotMom;//!
126 Float_t fTreeVariablePosEta; //!
127 Float_t fTreeVariableNegEta; //!
129 Float_t fTreeVariableVertexZ; //!
131 Int_t fTreeVariableLeastNbrCrossedRows;//!
132 Float_t fTreeVariableLeastRatioCrossedRowsOverFindable;//!
133 Int_t fTreeVariableMultiplicity;//!
134 Int_t fTreeVariableMultiplicityMC;//!
136 //===========================================================================================
138 //===========================================================================================
140 TH1F *fHistV0MultiplicityBeforeTrigSel; //! V0 multiplicity distribution
141 TH1F *fHistV0MultiplicityForTrigEvt; //! V0 multiplicity distribution
142 TH1F *fHistV0MultiplicityForSelEvt; //! V0 multiplicity distribution
143 TH1F *fHistV0MultiplicityForSelEvtNoTPCOnly; //! V0 multiplicity distribution
144 TH1F *fHistV0MultiplicityForSelEvtNoTPCOnlyNoPileup; //! V0 multiplicity distribution
146 TH1F *fHistMultiplicityBeforeTrigSel; //! multiplicity distribution
147 TH1F *fHistMultiplicityForTrigEvt; //! multiplicity distribution
148 TH1F *fHistMultiplicity; //! multiplicity distribution
149 TH1F *fHistMultiplicityNoTPCOnly; //! multiplicity distribution
150 TH1F *fHistMultiplicityNoTPCOnlyNoPileup; //! multiplicity distribution
152 //Raw Data for J/Psi paper Technique
153 TH2F *f2dHistMultiplicityVsTrueBeforeTrigSel; //! multiplicity distribution
154 TH2F *f2dHistMultiplicityVsTrueForTrigEvt; //! multiplicity distribution
155 TH2F *f2dHistMultiplicityVsTrue; //! multiplicity distribution
156 TH2F *f2dHistMultiplicityVsTrueNoTPCOnly; //! multiplicity distribution
157 TH2F *f2dHistMultiplicityVsTrueNoTPCOnlyNoPileup; //! multiplicity distribution
159 //Raw Data for Vertex Z position estimator change
160 TH2F *f2dHistMultiplicityVsVertexZBeforeTrigSel; //! multiplicity distribution
161 TH2F *f2dHistMultiplicityVsVertexZForTrigEvt; //! multiplicity distribution
162 TH2F *f2dHistMultiplicityVsVertexZ; //! multiplicity distribution
163 TH2F *f2dHistMultiplicityVsVertexZNoTPCOnly; //! multiplicity distribution
164 TH2F *f2dHistMultiplicityVsVertexZNoTPCOnlyNoPileup; //! multiplicity distribution
166 TH1F *fHistGenVertexZBeforeTrigSel; //! multiplicity distribution
167 TH1F *fHistGenVertexZForTrigEvt; //! multiplicity distribution
168 TH1F *fHistGenVertexZ; //! multiplicity distribution
169 TH1F *fHistGenVertexZNoTPCOnly; //! multiplicity distribution
170 TH1F *fHistGenVertexZNoTPCOnlyNoPileup; //! multiplicity distribution
172 //---> Filled At Analysis Scope
174 TH3F *f3dHistPrimAnalysisPtVsYVsMultLambda; //! Lambda
175 TH3F *f3dHistPrimAnalysisPtVsYVsMultAntiLambda; //! AntiLambda
176 TH3F *f3dHistPrimAnalysisPtVsYVsMultK0Short; //! K0Short
178 //---> Containers for monte carlo information for calculating efficiency!
180 TH3F *f3dHistPrimRawPtVsYVsMultLambda; //! Lambda
181 TH3F *f3dHistPrimRawPtVsYVsMultAntiLambda; //! AntiLambda
182 TH3F *f3dHistPrimRawPtVsYVsMultK0Short; //! K0Short
184 TH3F *f3dHistPrimRawPtVsYVsMultMCLambda; //! Lambda
185 TH3F *f3dHistPrimRawPtVsYVsMultMCAntiLambda; //! AntiLambda
186 TH3F *f3dHistPrimRawPtVsYVsMultMCK0Short; //! K0Short
188 TH3F *f3dHistPrimRawPtVsYVsVertexZLambda; //! Lambda
189 TH3F *f3dHistPrimRawPtVsYVsVertexZAntiLambda; //! AntiLambda
190 TH3F *f3dHistPrimRawPtVsYVsVertexZK0Short; //! K0Short
192 TH3F *f3dHistPrimCloseToPVPtVsYVsMultLambda; //! Lambda
193 TH3F *f3dHistPrimCloseToPVPtVsYVsMultAntiLambda; //! AntiLambda
194 TH3F *f3dHistPrimCloseToPVPtVsYVsMultK0Short; //! K0Short
196 //---> Filled vs Decay Length
198 TH3F *f3dHistPrimRawPtVsYVsDecayLengthLambda; //! Lambda
199 TH3F *f3dHistPrimRawPtVsYVsDecayLengthAntiLambda; //! AntiLambda
200 TH3F *f3dHistPrimRawPtVsYVsDecayLengthK0Short; //! K0Short
202 //---> Needed for FeedDown Corrections
204 TH3F *f3dHistGenPtVsYVsMultXiMinus; //! Generated Xi- Distrib
205 TH3F *f3dHistGenPtVsYVsMultXiPlus; //! Generated Xi+ Distrib
206 TH3F *f3dHistGenPtVsYVsMultOmegaMinus; //! Generated Omega- Distrib
207 TH3F *f3dHistGenPtVsYVsMultOmegaPlus; //! Generated Omega+ Distrib
209 TH3F *f3dHistGenSelectedPtVsYVsMultXiMinus; //! Generated Xi- Distrib, at event selection level
210 TH3F *f3dHistGenSelectedPtVsYVsMultXiPlus; //! Generated Xi+ Distrib, at event selection level
211 TH3F *f3dHistGenSelectedPtVsYVsMultOmegaMinus; //! Generated Omega- Distrib, at event selection level
212 TH3F *f3dHistGenSelectedPtVsYVsMultOmegaPlus; //! Generated Omega+ Distrib, at event selection level
214 TH1F *fHistPVx; //! PVx distrib
215 TH1F *fHistPVy; //! PVy distrib
216 TH1F *fHistPVz; //! PVz distrib
217 TH1F *fHistPVxAnalysis; //! PVx distrib
218 TH1F *fHistPVyAnalysis; //! PVy distrib
219 TH1F *fHistPVzAnalysis; //! PVz distrib
220 TH1F *fHistPVxAnalysisHasHighPtLambda; //! PVx distrib
221 TH1F *fHistPVyAnalysisHasHighPtLambda; //! PVy distrib
222 TH1F *fHistPVzAnalysisHasHighPtLambda; //! PVz distrib
224 TH1F *fHistSwappedV0Counter; //! Swapped v0 counter
226 AliAnalysisTaskExtractPerformanceV0(const AliAnalysisTaskExtractPerformanceV0&); // not implemented
227 AliAnalysisTaskExtractPerformanceV0& operator=(const AliAnalysisTaskExtractPerformanceV0&); // not implemented
229 ClassDef(AliAnalysisTaskExtractPerformanceV0, 11);