]> git.uio.no Git - u/mrichter/AliRoot.git/blame - PWGLF/STRANGENESS/Cascades/AliAnalysisTaskExtractCascade.h
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
Stupid bug fix in new superlight mode (from Zurich airport)
[u/mrichter/AliRoot.git] / PWGLF / STRANGENESS / Cascades / AliAnalysisTaskExtractCascade.h
CommitLineData
76029adc 1/**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
3 * *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
6 * *
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 **************************************************************************/
15
16// +-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
17//
18// Modified version of AliAnalysisTaskCheckCascade.h
19// Used bits of code from AliAnalysisTaskCheckPerformanceStrange
20//
21// --- David Dobrigkeit Chinellato
22//
23// +-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
24
25#ifndef ALIANALYSISTASKEXTRACTCASCADE_H
26#define ALIANALYSISTASKEXTRACTCASCADE_H
27
28class TList;
29class TH1F;
30class TH2F;
31class TH3F;
32class TVector3;
33class THnSparse;
34
35class AliESDpid;
36class AliESDtrackCuts;
37class AliESDEvent;
38class AliPhysicsSelection;
39class AliCFContainer;
2add50c1 40class AliAnalysisUtils;
76029adc 41
42//#include "TString.h"
43//#include "AliESDtrackCuts.h"
44#include "AliAnalysisTaskSE.h"
45
46class AliAnalysisTaskExtractCascade : public AliAnalysisTaskSE {
df14ae99 47public:
48 AliAnalysisTaskExtractCascade();
49 AliAnalysisTaskExtractCascade(const char *name);
50 virtual ~AliAnalysisTaskExtractCascade();
51
52 virtual void UserCreateOutputObjects();
53 virtual void UserExec(Option_t *option);
54 virtual void Terminate(Option_t *);
55 Double_t MyRapidity(Double_t rE, Double_t rPz) const;
56
57 void SetIsNuclear (Bool_t lIsNuclear = kTRUE ) { fkIsNuclear = lIsNuclear; }
58 void SetINT7Trigger (Bool_t lSwitchINT7 = kTRUE ) { fkSwitchINT7 = lSwitchINT7; }
59 void SetCentralityEstimator (TString lCentralityEstimator = "V0M" ) { fCentralityEstimator = lCentralityEstimator; }
60 void SetpAVertexSelection (Bool_t lpAVertexSelection = kTRUE) {fkpAVertexSelection = lpAVertexSelection; }
61 void SetEtaRefMult ( Double_t lEtaRefMult = 0.5 ) { fEtaRefMult = lEtaRefMult; }
62
63 //---------------------------------------------------------------------------------------
64 //Task Configuration: Meant to enable quick re-execution of vertexer if needed
65 void SetRunVertexers ( Bool_t lRunVertexers = kTRUE) { fkRunVertexers = lRunVertexers; }
66 //---------------------------------------------------------------------------------------
67 //Set Debug Mode
68 void SetDebugMode ( Bool_t lDebugMode = kTRUE) { fkDebugMode = lDebugMode; }
47d2266a 69 //---------------------------------------------------------------------------------------
70 //Set Peripheral event debugging mode (Pb-Pb X-check)
b21a3912 71 void SetSelectCentrality ( Bool_t lSelectCentrality = kTRUE, Double_t lCentSelLow = 0.0, Double_t lCentSelHigh = 10.0) {
72 fkSelectCentrality = lSelectCentrality;
73 fCentSel_Low = lCentSelLow;
74 fCentSel_High = lCentSelHigh;
75 }
47d2266a 76 //---------------------------------------------------------------------------------------
165ff019 77 void SetLowPtCutoff ( Double_t lLowPtCutoff = 1.0) {
78 fLowPtCutoff = lLowPtCutoff;
79 }
80 //---------------------------------------------------------------------------------------
df14ae99 81 void SetSuperLightMode ( Bool_t lSuperLight = kTRUE) {
82 fkSuperLight = lSuperLight;
83 }
84 //---------------------------------------------------------------------------------------
35709701 85 void SetLightMode ( Bool_t lLight = kTRUE) {
86 fkLight = lLight;
87 }
88 //---------------------------------------------------------------------------------------
165ff019 89 //Setters for the V0 Vertexer Parameters
df14ae99 90 void SetV0VertexerMaxChisquare ( Double_t lParameter ){ fV0VertexerSels[0] = lParameter; }
91 void SetV0VertexerDCAFirstToPV ( Double_t lParameter ){ fV0VertexerSels[1] = lParameter; }
92 void SetV0VertexerDCASecondtoPV ( Double_t lParameter ){ fV0VertexerSels[2] = lParameter; }
93 void SetV0VertexerDCAV0Daughters ( Double_t lParameter ){ fV0VertexerSels[3] = lParameter; }
94 void SetV0VertexerCosinePA ( Double_t lParameter ){ fV0VertexerSels[4] = lParameter; }
95 void SetV0VertexerMinRadius ( Double_t lParameter ){ fV0VertexerSels[5] = lParameter; }
96 void SetV0VertexerMaxRadius ( Double_t lParameter ){ fV0VertexerSels[6] = lParameter; }
97 //---------------------------------------------------------------------------------------
98 //Setters for the Cascade Vertexer Parameters
99 void SetCascVertexerMaxChisquare ( Double_t lParameter ){ fCascadeVertexerSels[0] = lParameter; }
100 void SetCascVertexerMinV0ImpactParameter ( Double_t lParameter ){ fCascadeVertexerSels[1] = lParameter; }
101 void SetCascVertexerV0MassWindow ( Double_t lParameter ){ fCascadeVertexerSels[2] = lParameter; }
102 void SetCascVertexerDCABachToPV ( Double_t lParameter ){ fCascadeVertexerSels[3] = lParameter; }
103 void SetCascVertexerDCACascadeDaughters ( Double_t lParameter ){ fCascadeVertexerSels[4] = lParameter; }
104 void SetCascVertexerCascadeCosinePA ( Double_t lParameter ){ fCascadeVertexerSels[5] = lParameter; }
105 void SetCascVertexerCascadeMinRadius ( Double_t lParameter ){ fCascadeVertexerSels[6] = lParameter; }
106 void SetCascVertexerCascadeMaxRadius ( Double_t lParameter ){ fCascadeVertexerSels[7] = lParameter; }
107 //---------------------------------------------------------------------------------------
108 //Setters for SuperLight analysis
109 void SetSuperLightV0Radius ( Double_t lParameter ){ fCut_V0Radius = lParameter; }
110 void SetSuperLightCascRadius ( Double_t lParameter ){ fCut_CascRadius = lParameter; }
111 void SetSuperLightV0Mass ( Double_t lParameter ){ fCut_V0Mass = lParameter; }
112 void SetSuperLightV0CosPA ( Double_t lParameter ){ fCut_V0CosPA = lParameter; }
113 void SetSuperLightCascCosPA ( Double_t lParameter ){ fCut_CascCosPA = lParameter; }
114 void SetSuperLightDCANegToPV ( Double_t lParameter ){ fCut_DCANegToPV = lParameter; }
115 void SetSuperLightDCAPosToPV ( Double_t lParameter ){ fCut_DCAPosToPV = lParameter; }
116 void SetSuperLightDCABachToPV ( Double_t lParameter ){ fCut_DCABachToPV = lParameter; }
117 void SetSuperLightDCAV0Daughters ( Double_t lParameter ){ fCut_DCAV0Daughters = lParameter; }
118 void SetSuperLightDCACascDaughters ( Double_t lParameter ){ fCut_DCACascDaughters = lParameter; }
119 void SetSuperLightDCAV0ToPV ( Double_t lParameter ){ fCut_DCAV0ToPV = lParameter; }
120 void SetSuperLightCTau ( Double_t lParameter ){ fCut_CTau = lParameter; }
121 //---------------------------------------------------------------------------------------
122
123 //---------------------------------------------------------------------------------------
124private:
125 // Note : In ROOT, "//!" means "do not stream the data from Master node to Worker node" ...
126 // your data member object is created on the worker nodes and streaming is not needed.
127 // http://root.cern.ch/download/doc/11InputOutput.pdf, page 14
128 TList *fListHist; //! List of Cascade histograms
129 TTree *fTreeCascade; //! Output Tree, Cascades
130
131 //Objects that have to be streamed:
132 AliPIDResponse *fPIDResponse; // PID response object
133 AliESDtrackCuts *fESDtrackCuts; // ESD track cuts used for primary track definition
134 AliAnalysisUtils *fUtils; // analysis utils (for pA vertex selection)
135
136 //Objects Controlling Task Behaviour
137 // (have to be streamed too or configuration is lost)
138 Bool_t fkIsNuclear; //if true, replace multiplicity est. by centrality (default FALSE)
139 Bool_t fkSwitchINT7; //if true, skip FASTOnly (default FALSE)
140 TString fCentralityEstimator; //Centrality Estimator String value (default V0M)
141 Bool_t fkpAVertexSelection; //if true, select vertex with pPb Methods
142 Double_t fEtaRefMult; //Reference multiplicity eta
143 //Objects Controlling Task Behaviour: has to be streamed!
144 Bool_t fkRunVertexers; // if true, re-run vertexer with loose cuts. CARE MUST BE TAKEN in PbPb!
145 Bool_t fkDebugMode; //Store extra debugging information in tree
b21a3912 146 Bool_t fkSelectCentrality; //Switch to skip anything other than 60-80% V0M
df14ae99 147 Double_t fV0VertexerSels[7]; // Array to store the 7 values for the different selections V0 related
148 Double_t fCascadeVertexerSels[8]; // Array to store the 8 values for the different selections Casc. related
b21a3912 149
150 Double_t fCentSel_Low;
151 Double_t fCentSel_High;
165ff019 152
153 Double_t fLowPtCutoff; //Reduction of data volume
df14ae99 154
b68d5ca6 155 //Double_t fV0Sels[7]; // Array to store the 7 values for the different selections V0 related
156 //Double_t fCascSels[8]; // Array to store the 8 values for the different selections Casc. related
df14ae99 157
158 //===========================================================================================
159 // Variables for tree, cascades
160 //===========================================================================================
161
162 Int_t fTreeCascVarCharge; //!
163 Float_t fTreeCascVarMassAsXi; //!
164 Float_t fTreeCascVarMassAsOmega; //!
165 Float_t fTreeCascVarPt; //!
166 Float_t fTreeCascVarPtMC; //!
167 Float_t fTreeCascVarRapMC; //!
168 Float_t fTreeCascVarRapXi; //!
169 Float_t fTreeCascVarRapOmega; //!
170 Float_t fTreeCascVarNegEta; //!
171 Float_t fTreeCascVarPosEta; //!
172 Float_t fTreeCascVarBachEta; //!
173 Float_t fTreeCascVarDCACascDaughters; //!
174 Float_t fTreeCascVarDCABachToPrimVtx; //!
175 Float_t fTreeCascVarDCAV0Daughters; //!
176 Float_t fTreeCascVarDCAV0ToPrimVtx; //!
177 Float_t fTreeCascVarDCAPosToPrimVtx; //!
178 Float_t fTreeCascVarDCANegToPrimVtx; //!
179 Float_t fTreeCascVarCascCosPointingAngle; //!
180 Float_t fTreeCascVarCascRadius; //!
181 Float_t fTreeCascVarV0Mass; //!
182 Float_t fTreeCascVarV0CosPointingAngle; //!
183 Float_t fTreeCascVarV0CosPointingAngleSpecial; //!
184 Float_t fTreeCascVarV0Radius; //!
185 Int_t fTreeCascVarLeastNbrClusters; //!
186 Int_t fTreeCascVarMultiplicity; //!
187 Int_t fTreeCascVarMultiplicityV0A; //!
188 Int_t fTreeCascVarMultiplicityZNA; //!
189 Int_t fTreeCascVarMultiplicityTRK; //!
190 Int_t fTreeCascVarMultiplicitySPD; //!
191 Float_t fTreeCascVarDistOverTotMom; //!
192 Int_t fTreeCascVarPID; //!
193 Int_t fTreeCascVarPIDBachelor; //!
194 Int_t fTreeCascVarPIDNegative; //!
195 Int_t fTreeCascVarPIDPositive; //!
196 Float_t fTreeCascVarBachTransMom; //!
197 Float_t fTreeCascVarPosTransMom; //!
198 Float_t fTreeCascVarNegTransMom; //!
199 Float_t fTreeCascVarPosTransMomMC; //!
200 Float_t fTreeCascVarNegTransMomMC; //!
201
202 Float_t fTreeCascVarNegNSigmaPion; //!
203 Float_t fTreeCascVarNegNSigmaProton; //!
204 Float_t fTreeCascVarPosNSigmaPion; //!
205 Float_t fTreeCascVarPosNSigmaProton; //!
206 Float_t fTreeCascVarBachNSigmaPion; //!
207 Float_t fTreeCascVarBachNSigmaKaon; //!
208
209 Bool_t fTreeCascVarkITSRefitBachelor; //!
210 Bool_t fTreeCascVarkITSRefitNegative; //!
211 Bool_t fTreeCascVarkITSRefitPositive; //!
d55fa19c 212
213 //Debugging information, if requested
214 //Part A: EbyE info, Run number
215 Int_t fTreeCascVarRunNumber; //!
216 ULong64_t fTreeCascVarEventNumber; //!
217
218 //Part B: Shared Clusters
80521867 219 Int_t fTreeCascVarNegClusters; //!
220 Int_t fTreeCascVarPosClusters; //!
221 Int_t fTreeCascVarBachClusters; //!
d55fa19c 222 Int_t fTreeCascVarNegSharedClusters; //!
223 Int_t fTreeCascVarPosSharedClusters; //!
224 Int_t fTreeCascVarBachSharedClusters; //!
225
226 //Part C: All momenta
227 Float_t fTreeCascVarNegPx; //!
228 Float_t fTreeCascVarNegPy; //!
229 Float_t fTreeCascVarNegPz; //!
230 Float_t fTreeCascVarPosPx; //!
231 Float_t fTreeCascVarPosPy; //!
232 Float_t fTreeCascVarPosPz; //!
233 Float_t fTreeCascVarBachPx; //!
234 Float_t fTreeCascVarBachPy; //!
235 Float_t fTreeCascVarBachPz; //!
df14ae99 236
3eb670b7 237 Float_t fTreeCascVarV0DecayX; //!
238 Float_t fTreeCascVarV0DecayY; //!
239 Float_t fTreeCascVarV0DecayZ; //!
240 Float_t fTreeCascVarCascadeDecayX; //!
241 Float_t fTreeCascVarCascadeDecayY; //!
242 Float_t fTreeCascVarCascadeDecayZ; //!
d55fa19c 243
df14ae99 244 Bool_t fTreeCascVarBadCascadeJai; //!
45ecf7a7 245 Float_t fTreeCascVarDeltaDCA; //!
246
df14ae99 247 //===========================================================================================
248 // Histograms
249 //===========================================================================================
250
251 TH1F *fHistV0MultiplicityBeforeTrigSel; //! V0 multiplicity distribution
252 TH1F *fHistV0MultiplicityForTrigEvt; //! V0 multiplicity distribution
253 TH1F *fHistV0MultiplicityForSelEvt; //! V0 multiplicity distribution
254 TH1F *fHistV0MultiplicityForSelEvtNoTPCOnly; //! V0 multiplicity distribution
255 TH1F *fHistV0MultiplicityForSelEvtNoTPCOnlyNoPileup; //! V0 multiplicity distribution
256
257 //V0M Centrality (default)
258 TH1F *fHistMultiplicityBeforeTrigSel; //! multiplicity distribution
259 TH1F *fHistMultiplicityForTrigEvt; //! multiplicity distribution
260 TH1F *fHistMultiplicity; //! multiplicity distribution
261 TH1F *fHistMultiplicityNoTPCOnly; //! multiplicity distribution
262 TH1F *fHistMultiplicityNoTPCOnlyNoPileup; //! multiplicity distribution
263
264 //V0A Centrality
265 TH1F *fHistMultiplicityV0ABeforeTrigSel; //! multiplicity distribution
5d43de53 266 TH1F *fHistMultiplicityV0AForTrigEvt; //! multiplicity distribution
267 TH1F *fHistMultiplicityV0A; //! multiplicity distribution
268 TH1F *fHistMultiplicityV0ANoTPCOnly; //! multiplicity distribution
269 TH1F *fHistMultiplicityV0ANoTPCOnlyNoPileup; //! multiplicity distribution
df14ae99 270
271 //ZNA Centrality
272 TH1F *fHistMultiplicityZNABeforeTrigSel; //! multiplicity distribution
5d43de53 273 TH1F *fHistMultiplicityZNAForTrigEvt; //! multiplicity distribution
274 TH1F *fHistMultiplicityZNA; //! multiplicity distribution
275 TH1F *fHistMultiplicityZNANoTPCOnly; //! multiplicity distribution
276 TH1F *fHistMultiplicityZNANoTPCOnlyNoPileup; //! multiplicity distribution
df14ae99 277
278 //TRK Centrality
279 TH1F *fHistMultiplicityTRKBeforeTrigSel; //! multiplicity distribution
5d43de53 280 TH1F *fHistMultiplicityTRKForTrigEvt; //! multiplicity distribution
281 TH1F *fHistMultiplicityTRK; //! multiplicity distribution
282 TH1F *fHistMultiplicityTRKNoTPCOnly; //! multiplicity distribution
283 TH1F *fHistMultiplicityTRKNoTPCOnlyNoPileup; //! multiplicity distribution
df14ae99 284
285 //SPD Centrality
286 TH1F *fHistMultiplicitySPDBeforeTrigSel; //! multiplicity distribution
5d43de53 287 TH1F *fHistMultiplicitySPDForTrigEvt; //! multiplicity distribution
288 TH1F *fHistMultiplicitySPD; //! multiplicity distribution
289 TH1F *fHistMultiplicitySPDNoTPCOnly; //! multiplicity distribution
290 TH1F *fHistMultiplicitySPDNoTPCOnlyNoPileup; //! multiplicity distribution
df14ae99 291
292 //---> Generated Histograms
293
294 TH1F *fHistPVx; //! PVx distrib
295 TH1F *fHistPVy; //! PVy distrib
296 TH1F *fHistPVz; //! PVz distrib
297 TH1F *fHistPVxAnalysis; //! PVx distrib
298 TH1F *fHistPVyAnalysis; //! PVy distrib
299 TH1F *fHistPVzAnalysis; //! PVz distrib
300
301 //=======================================================
302 // --- Superlight Output Mode - Experimental ---
303 //=======================================================
304
35709701 305 Bool_t fkLight; //Switch for intermediate mode
df14ae99 306 Bool_t fkSuperLight; //Switch for super light output mode
307
308 //Selection criteria for superlight analysis
309 Double_t fCut_V0Radius;
310 Double_t fCut_CascRadius;
311 Double_t fCut_V0Mass;
312 Double_t fCut_V0CosPA;
313 Double_t fCut_CascCosPA;
314 Double_t fCut_DCANegToPV;
315 Double_t fCut_DCAPosToPV;
316 Double_t fCut_DCABachToPV;
317 Double_t fCut_DCAV0Daughters;
318 Double_t fCut_DCACascDaughters;
319 Double_t fCut_DCAV0ToPV;
320 Double_t fCut_CTau;
321
322 //---> Super-lightweight histogram output (rebinnable)
323
324 TH2F *f2dHist_MassVsPt_XiMinus;
325 TH2F *f2dHist_MassVsPt_XiPlus;
326 TH2F *f2dHist_MassVsPt_OmegaMinus;
327 TH2F *f2dHist_MassVsPt_OmegaPlus;
328
329 AliAnalysisTaskExtractCascade(const AliAnalysisTaskExtractCascade&); // not implemented
330 AliAnalysisTaskExtractCascade& operator=(const AliAnalysisTaskExtractCascade&); // not implemented
331
332 ClassDef(AliAnalysisTaskExtractCascade, 11);
76029adc 333};
334
335#endif