]>
Commit | Line | Data |
---|---|---|
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 | ||
28 | class TList; | |
29 | class TH1F; | |
30 | class TH2F; | |
31 | class TH3F; | |
32 | class TVector3; | |
33 | class THnSparse; | |
34 | ||
35 | class AliESDpid; | |
36 | class AliESDtrackCuts; | |
37 | class AliAnalysisUtils; | |
b40a01ae | 38 | class AliPPVsMultUtils; |
af0d4791 | 39 | class AliESDEvent; |
40 | class AliPhysicsSelection; | |
41 | class AliCFContainer; | |
42 | ||
43 | //#include "TString.h" | |
44 | //#include "AliESDtrackCuts.h" | |
45 | //#include "AliAnalysisTaskSE.h" | |
46 | ||
47 | class 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 |