]> git.uio.no Git - u/mrichter/AliRoot.git/blame - PWGLF/STRANGENESS/Cascades/AliAnalysisTaskStrangenessVsMultiplicity.h
Merge branch 'TPCdev' of https://git.cern.ch/reps/AliRoot into TPCdev
[u/mrichter/AliRoot.git] / PWGLF / STRANGENESS / Cascades / AliAnalysisTaskStrangenessVsMultiplicity.h
CommitLineData
af0d4791 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 AliAnalysisTaskStrangenessVsMultiplicity_H
26#define AliAnalysisTaskStrangenessVsMultiplicity_H
27
28class TList;
29class TH1F;
30class TH2F;
31class TH3F;
32class TVector3;
33class THnSparse;
34
35class AliESDpid;
36class AliESDtrackCuts;
37class AliAnalysisUtils;
b40a01ae 38class AliPPVsMultUtils;
af0d4791 39class AliESDEvent;
40class AliPhysicsSelection;
41class AliCFContainer;
42
43//#include "TString.h"
44//#include "AliESDtrackCuts.h"
45//#include "AliAnalysisTaskSE.h"
46
47class AliAnalysisTaskStrangenessVsMultiplicity : public AliAnalysisTaskSE {
48 public:
49 AliAnalysisTaskStrangenessVsMultiplicity();
50 AliAnalysisTaskStrangenessVsMultiplicity(const char *name);
51 virtual ~AliAnalysisTaskStrangenessVsMultiplicity();
52
53 virtual void UserCreateOutputObjects();
54 virtual void UserExec(Option_t *option);
55 virtual void Terminate(Option_t *);
56 Double_t MyRapidity(Double_t rE, Double_t rPz) const;
57
58 void SetSaveV0s (Bool_t lSaveV0s = kTRUE ) { fkSaveV0Tree = lSaveV0s; }
59 void SetSaveCascades (Bool_t lSaveCascades = kTRUE ) { fkSaveCascadeTree = lSaveCascades; }
60
61//---------------------------------------------------------------------------------------
62 //Task Configuration: Meant to enable quick re-execution of vertexer if needed
63 void SetRunVertexers ( Bool_t lRunVertexers = kTRUE) { fkRunVertexers = lRunVertexers; }
758b26b0 64//---------------------------------------------------------------------------------------
65 //Task Configuration: Skip Event Selections after trigger (VZERO test)
66 void SetSkipEventSelection ( Bool_t lSkipEventSelection = kTRUE) { fkSkipEventSelection = lSkipEventSelection; }
af0d4791 67//---------------------------------------------------------------------------------------
68//Setters for the V0 Vertexer Parameters
69 void SetV0VertexerMaxChisquare ( Double_t lParameter ){ fV0VertexerSels[0] = lParameter; }
70 void SetV0VertexerDCAFirstToPV ( Double_t lParameter ){ fV0VertexerSels[1] = lParameter; }
71 void SetV0VertexerDCASecondtoPV ( Double_t lParameter ){ fV0VertexerSels[2] = lParameter; }
72 void SetV0VertexerDCAV0Daughters ( Double_t lParameter ){ fV0VertexerSels[3] = lParameter; }
73 void SetV0VertexerCosinePA ( Double_t lParameter ){ fV0VertexerSels[4] = lParameter; }
74 void SetV0VertexerMinRadius ( Double_t lParameter ){ fV0VertexerSels[5] = lParameter; }
75 void SetV0VertexerMaxRadius ( Double_t lParameter ){ fV0VertexerSels[6] = lParameter; }
76//---------------------------------------------------------------------------------------
77//Setters for the Cascade Vertexer Parameters
78 void SetCascVertexerMaxChisquare ( Double_t lParameter ){ fCascadeVertexerSels[0] = lParameter; }
79 void SetCascVertexerMinV0ImpactParameter ( Double_t lParameter ){ fCascadeVertexerSels[1] = lParameter; }
80 void SetCascVertexerV0MassWindow ( Double_t lParameter ){ fCascadeVertexerSels[2] = lParameter; }
81 void SetCascVertexerDCABachToPV ( Double_t lParameter ){ fCascadeVertexerSels[3] = lParameter; }
82 void SetCascVertexerDCACascadeDaughters ( Double_t lParameter ){ fCascadeVertexerSels[4] = lParameter; }
83 void SetCascVertexerCascadeCosinePA ( Double_t lParameter ){ fCascadeVertexerSels[5] = lParameter; }
84 void SetCascVertexerCascadeMinRadius ( Double_t lParameter ){ fCascadeVertexerSels[6] = lParameter; }
85 void SetCascVertexerCascadeMaxRadius ( Double_t lParameter ){ fCascadeVertexerSels[7] = lParameter; }
86//---------------------------------------------------------------------------------------
87
88 private:
89 // Note : In ROOT, "//!" means "do not stream the data from Master node to Worker node" ...
90 // your data member object is created on the worker nodes and streaming is not needed.
91 // http://root.cern.ch/download/doc/11InputOutput.pdf, page 14
92 TList *fListHist; //! List of Cascade histograms
93 TTree *fTreeEvent; //! Output Tree, Events
94 TTree *fTreeV0; //! Output Tree, V0s
95 TTree *fTreeCascade; //! Output Tree, Cascades
96
97 AliPIDResponse *fPIDResponse; // PID response object
98 AliESDtrackCuts *fESDtrackCuts; // ESD track cuts used for primary track definition
15863da5 99 AliPPVsMultUtils *fPPVsMultUtils; //
100 AliAnalysisUtils *fUtils; // analysis utils (for MV pileup selection)
af0d4791 101
102 //Objects Controlling Task Behaviour
103 Bool_t fkSaveV0Tree; //if true, save TTree
104 Bool_t fkSaveCascadeTree; //if true, save TTree
105
106 //Objects Controlling Task Behaviour: has to be streamed!
107 Bool_t fkRunVertexers; // if true, re-run vertexer with loose cuts *** only for CASCADES! ***
d8c18e93 108 Bool_t fkSkipEventSelection; // if true, will only perform TRIGGER selection (currently kMB, to change)
109 Bool_t fkApplyTrackletsVsClustersCut; //if true, applies Tracklet vs clusters cut together with PS
af0d4791 110 Double_t fV0VertexerSels[7]; // Array to store the 7 values for the different selections V0 related
111 Double_t fCascadeVertexerSels[8]; // Array to store the 8 values for the different selections Casc. related
112
113//===========================================================================================
114// Variables for Event Tree
115//===========================================================================================
116 Float_t fAmplitude_V0A; //!
117 Float_t fAmplitude_V0C; //!
118 Float_t fAmplitude_V0AEq; //!
119 Float_t fAmplitude_V0CEq; //!
120 Float_t fCentrality_V0A; //!
121 Float_t fCentrality_V0C; //!
122 Float_t fCentrality_V0M; //!
123 Float_t fCentrality_V0AEq; //!
124 Float_t fCentrality_V0CEq; //!
b40a01ae 125 Float_t fCentrality_V0MEq; //!
126 Float_t fCustomCentrality_V0M; //!
127 Float_t fCustomCentrality_V0MEq; //!
af0d4791 128 Int_t fRefMultEta5; //!
129 Int_t fRefMultEta8; //!
130 Int_t fRunNumber; //!
47529a6d 131
132 //Event Characterization Variables - optional
133 Bool_t fEvSel_HasAtLeastSPDVertex; //!
134 Bool_t fEvSel_VtxZCut; //!
135 Bool_t fEvSel_IsNotPileup; //!
15863da5 136 Bool_t fEvSel_IsNotPileupMV; //!
47529a6d 137 Bool_t fEvSel_IsNotPileupInMultBins; //!
7f83b30f 138 Bool_t fEvSel_Triggered; //!
47529a6d 139
140 //Other Selections: more dedicated filtering to be studied!
141 Int_t fEvSel_nTracklets; //!
142 Int_t fEvSel_nSPDClusters; //!
7f83b30f 143 Float_t fEvSel_VtxZ; //! the actual value
1b914674 144 Int_t fEvSel_nSPDPrimVertices; //! pileup vertices
145 Float_t fEvSel_distZ; //! distance between largest vertices
212fd60c 146 Int_t fEvSel_nContributors; //!
147 Int_t fEvSel_nContributorsPileup; //!
47529a6d 148
af0d4791 149
161c57ac 150//===========================================================================================
151// Variables for V0 Tree
152//===========================================================================================
153 Float_t fTreeVariableChi2V0; //!
154 Float_t fTreeVariableDcaV0Daughters; //!
155 Float_t fTreeVariableDcaV0ToPrimVertex; //!
156 Float_t fTreeVariableDcaPosToPrimVertex; //!
157 Float_t fTreeVariableDcaNegToPrimVertex; //!
158 Float_t fTreeVariableV0CosineOfPointingAngle; //!
159 Float_t fTreeVariableV0Radius; //!
160 Float_t fTreeVariablePt; //!
161 Float_t fTreeVariableRapK0Short; //!
162 Float_t fTreeVariableRapLambda; //!
163 Float_t fTreeVariableInvMassK0s; //!
164 Float_t fTreeVariableInvMassLambda; //!
165 Float_t fTreeVariableInvMassAntiLambda; //!
166 Float_t fTreeVariableAlphaV0; //!
167 Float_t fTreeVariablePtArmV0;//!
168 Float_t fTreeVariableNegEta; //!
169 Float_t fTreeVariablePosEta; //!
170
171 Float_t fTreeVariableNSigmasPosProton; //!
172 Float_t fTreeVariableNSigmasPosPion; //!
173 Float_t fTreeVariableNSigmasNegProton; //!
174 Float_t fTreeVariableNSigmasNegPion; //!
175
176 Float_t fTreeVariableDistOverTotMom;//!
177 Int_t fTreeVariableLeastNbrCrossedRows;//!
178 Float_t fTreeVariableLeastRatioCrossedRowsOverFindable;//!
179
180 //Event Multiplicity Variables
181 Float_t fTreeVariableCentV0A; //!
182 Float_t fTreeVariableCentV0C; //!
183 Float_t fTreeVariableCentV0M; //!
184 Float_t fTreeVariableCentV0AEq; //!
185 Float_t fTreeVariableCentV0CEq; //!
186 Float_t fTreeVariableCentV0MEq; //!
6cc116f9 187 Float_t fTreeVariableCustomCentV0M; //!
161c57ac 188 Float_t fTreeVariableAmpV0A; //!
189 Float_t fTreeVariableAmpV0C; //!
190 Float_t fTreeVariableAmpV0AEq; //!
191 Float_t fTreeVariableAmpV0CEq; //!
192 Int_t fTreeVariableRefMultEta8; //!
193 Int_t fTreeVariableRefMultEta5; //!
194 Int_t fTreeVariableRunNumber; //! //want to re-quantile per run? here's your ticket
195
af0d4791 196//===========================================================================================
197// Variables for Cascade Candidate Tree
198//===========================================================================================
199 Int_t fTreeCascVarCharge; //!
200 Float_t fTreeCascVarMassAsXi; //!
201 Float_t fTreeCascVarMassAsOmega; //!
202 Float_t fTreeCascVarPt; //!
203 Float_t fTreeCascVarRapXi; //!
204 Float_t fTreeCascVarRapOmega; //!
205 Float_t fTreeCascVarNegEta; //!
206 Float_t fTreeCascVarPosEta; //!
207 Float_t fTreeCascVarBachEta; //!
208 Float_t fTreeCascVarDCACascDaughters; //!
209 Float_t fTreeCascVarDCABachToPrimVtx; //!
210 Float_t fTreeCascVarDCAV0Daughters; //!
211 Float_t fTreeCascVarDCAV0ToPrimVtx; //!
212 Float_t fTreeCascVarDCAPosToPrimVtx; //!
213 Float_t fTreeCascVarDCANegToPrimVtx; //!
214 Float_t fTreeCascVarCascCosPointingAngle; //!
215 Float_t fTreeCascVarCascRadius; //!
216 Float_t fTreeCascVarV0Mass; //!
217 Float_t fTreeCascVarV0CosPointingAngle; //!
218 Float_t fTreeCascVarV0CosPointingAngleSpecial; //!
219 Float_t fTreeCascVarV0Radius; //!
220 Int_t fTreeCascVarLeastNbrClusters; //!
221 Float_t fTreeCascVarDistOverTotMom; //!
222
223 //TPC dEdx
224 Float_t fTreeCascVarNegNSigmaPion; //!
225 Float_t fTreeCascVarNegNSigmaProton; //!
226 Float_t fTreeCascVarPosNSigmaPion; //!
227 Float_t fTreeCascVarPosNSigmaProton; //!
228 Float_t fTreeCascVarBachNSigmaPion; //!
229 Float_t fTreeCascVarBachNSigmaKaon; //!
230
231 //Event Multiplicity Variables
232 Float_t fTreeCascVarCentV0A; //!
233 Float_t fTreeCascVarCentV0C; //!
234 Float_t fTreeCascVarCentV0M; //!
235 Float_t fTreeCascVarCentV0AEq; //!
236 Float_t fTreeCascVarCentV0CEq; //!
237 Float_t fTreeCascVarCentV0MEq; //!
6cc116f9 238 Float_t fTreeCascVarCustomCentV0M; //!
af0d4791 239 Float_t fTreeCascVarAmpV0A; //!
240 Float_t fTreeCascVarAmpV0C; //!
241 Float_t fTreeCascVarAmpV0AEq; //!
242 Float_t fTreeCascVarAmpV0CEq; //!
243 Int_t fTreeCascVarRefMultEta8; //!
244 Int_t fTreeCascVarRefMultEta5; //!
26d39641 245 Int_t fTreeCascVarRunNumber; //! //want to re-quantile per run? here's your ticket
af0d4791 246
247//===========================================================================================
248// Histograms
249//===========================================================================================
250
251 TH1D *fHistEventCounter; //!
252
253 AliAnalysisTaskStrangenessVsMultiplicity(const AliAnalysisTaskStrangenessVsMultiplicity&); // not implemented
254 AliAnalysisTaskStrangenessVsMultiplicity& operator=(const AliAnalysisTaskStrangenessVsMultiplicity&); // not implemented
255
256 ClassDef(AliAnalysisTaskStrangenessVsMultiplicity, 11);
257};
258
259#endif