AliPhysics  2b88e80 (2b88e80)
AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero.h
Go to the documentation of this file.
1 
2 #ifndef ALIANALYSISTASKNEUTRALMESONTOPIPLPIMIPIZERO_H
3 #define ALIANALYSISTASKNEUTRALMESONTOPIPLPIMIPIZERO_H
4 
5 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
6  * See cxx source for full Copyright notice */
7 
8 
9 #include "AliAnalysisTaskSE.h"
10 #include "AliV0ReaderV1.h"
11 #include "AliKFConversionPhoton.h"
12 #include "AliPrimaryPionSelector.h"
13 #include "AliConversionMesonCuts.h"
14 #include "AliConvEventCuts.h"
15 #include "AliCaloPhotonCuts.h"
16 #include "AliGammaConversionAODBGHandler.h"
17 #include "TProfile2D.h"
18 #include <vector>
19 
20 class AliESDInputHandler;
21 class AliMCEventHandler;
22 class AliESDEvent;
23 class AliESDtrack;
24 class AliESDtrackCuts;
25 class AliESDpidCuts;
26 class AliTriggerAnalysis;
27 
28 class AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero: public AliAnalysisTaskSE
29 {
30  public:
31 
32  AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero();
33  AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero( const char* name );
34  virtual ~AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero();
35 
36  virtual void UserExec(Option_t *);
37  virtual void UserCreateOutputObjects();
38  virtual Bool_t Notify();
39  virtual void Terminate(const Option_t *);
40 
41  void SetV0ReaderName(TString name){fV0ReaderName=name; return;}
42 
43  void SetMoveParticleAccordingToVertex(Bool_t flag){fMoveParticleAccordingToVertex = flag;}
44  void SetIsHeavyIon(Int_t flag){
45  if (flag == 1 || flag ==2 ){
46  fIsHeavyIon = 1;
47  } else {
48  fIsHeavyIon = 0;
49  }
50  }
51 
52  void SetIsMC(Bool_t isMC){fIsMC=isMC;}
53  void SetLightOutput(Bool_t flag){fDoLightOutput = flag;}
54  void SetEventCutList(Int_t nCuts, TList *CutArray){
55  fnCuts= nCuts;
56  fEventCutArray = CutArray;
57  }
58  void SetConversionCutList(TList *CutArray){ fGammaCutArray = CutArray;}
59  void SetClusterCutList(TList *CutArray){ fClusterCutArray = CutArray;}
60  void SetPionCutList(TList *CutArray){ fPionCutArray = CutArray;}
61  void SetNeutralPionCutList(TList *CutArray){ fNeutralPionMesonCutArray = CutArray; }
62  void SetMesonCutList(TList *CutArray){ fMesonCutArray = CutArray; }
63  void SetDoMesonQA(Int_t flag){ fDoMesonQA = flag; }
64  void SetNeutralPionMode(Int_t mode){fNeutralPionMode = mode; }
65  void SetTolerance(Double_t tol){fTolerance=tol;}
66 
67  private:
68 
69  void InitBack();
70 
71  // routines for photon selection from conversions
72  void ProcessConversionPhotonCandidates();
73  void ProcessTrueConversionPhotonCandidates(AliAODConversionPhoton*);
74 
75  // routines for photon selection from clusters
76  void ProcessCaloPhotonCandidates();
77  void ProcessTrueCaloPhotonCandidates(AliAODConversionPhoton *TruePhotonCandidate);
78 
79  void ProcessTrueMesonCandidates(AliAODConversionMother *Pi0Candidate, AliAODConversionMother *TrueNeutralPionCandidate, AliAODConversionPhoton *TrueVirtualGammaCandidate);
80  void MoveParticleAccordingToVertex(AliAODConversionMother* particle,const AliGammaConversionAODBGHandler::GammaConversionVertex *vertex);
81 
82  void FixPzToMatchPDGInvMassPi0(AliAODConversionMother* particle);
83 
84  // routines for neutral pion candidates from pure conversion
85  void ProcessNeutralPionCandidatesPureConversions();
86  void ProcessTrueNeutralPionCandidatesPureConversions(AliAODConversionMother *Pi0Candidate, AliAODConversionPhoton *TrueGammaCandidate0, AliAODConversionPhoton *TrueGammaCandidate1);
87  void ProcessTrueNeutralPionCandidatesPureConversionsAOD(AliAODConversionMother *Pi0Candidate, AliAODConversionPhoton *TrueGammaCandidate0, AliAODConversionPhoton *TrueGammaCandidate1);
88 
89  // routines for neutral pion candidates from pure calo
90  void ProcessNeutralPionCandidatesPureCalo();
91  void ProcessTrueNeutralPionCandidatesPureCalo(AliAODConversionMother *Pi0Candidate, AliAODConversionPhoton *TrueGammaCandidate0, AliAODConversionPhoton *TrueGammaCandidate1);
92 
93  // routines for neutral pion candidates from mixed conv + calo
94  void ProcessNeutralPionCandidatesMixedConvCalo();
95  void ProcessTrueNeutralPionCandidatesMixedConvCalo( AliAODConversionMother *Pi0Candidate, AliAODConversionPhoton *TrueGammaCandidate0, AliAODConversionPhoton *TrueGammaCandidate1);
96 
97  void ProcessPionCandidates();
98  void ProcessMCParticles();
99  void CalculateMesonCandidates();
100  void CalculateBackground();
101  void UpdateEventByEventData();
102 
103  Bool_t IsPiPlPiMiPiZeroDecay(TParticle *fMCMother) const;
104  Bool_t IsEtaPiPlPiMiPiZeroDaughter( Int_t label ) const;
105  Bool_t IsOmegaPiPlPiMiPiZeroDaughter( Int_t label ) const;
106  Bool_t GammaIsNeutralMesonPiPlPiMiPiZeroDaughter( Int_t label ) const;
107 
108  Bool_t CheckVectorForDoubleCount(vector<Int_t> &vec, Int_t tobechecked);
109 
110  Bool_t KinematicCut(AliAODConversionMother *negpion, AliAODConversionMother *pospion, AliAODConversionMother *neutpion, AliAODConversionMother *omega);
111 
112 
113  AliV0ReaderV1* fV0Reader; // V0Reader for basic conversion photon selection
114  TString fV0ReaderName;
115  AliPrimaryPionSelector* fPionSelector; // primary charged pion selector, basic selection of pi+,pi-
116  AliGammaConversionAODBGHandler** fBGHandlerPiPl; // BG handler Pos Pion
117  AliGammaConversionAODBGHandler** fBGHandlerPiMi; // BG handler Neg Pion
118  AliESDEvent* fESDEvent; // current event
119  AliMCEvent* fMCEvent; // current MC event
120  TList** fCutFolder; // list of output folders with main cut name
121  TList** fESDList; // list with main output histograms for data
122  TList** fTrueList; // list with validated reconstructed MC histograms
123  TList** fTrueTreeList; // list containing TTree's for MC True
124  TList** fMCList; // list with pure MC histograms
125  TList* fOutputContainer; // output container
126  TClonesArray* fReaderGammas; // array with photon from fV0Reader
127  vector<Int_t> fSelectorNegPionIndex; // array with pion indices for negative pions from fPionSelector
128  vector<Int_t> fSelectorPosPionIndex; // array with pion indices for positive pions from fPionSelector
129  TList* fGoodConvGammas; // good conv gammas after selection
130  TList* fClusterCandidates;
131  TList* fNeutralPionCandidates; // good neutral pion candidates
132  TList* fNeutralPionSidebandCandidates; // good neutral pion candidates from sideband
133  TList* fPosPionCandidates; // good positive pion candidates
134  TList* fNegPionCandidates; // good negative pion candidates
135  TList* fGoodVirtualParticles; // combination of pi+pi- candidates
136  TList* fEventCutArray; // array with event cuts
137  TList* fGammaCutArray; // array with Conversion Cuts
138  TList* fClusterCutArray; // array with Cluster Cuts
139  TList* fPionCutArray; // array with charged pion cuts
140  TList* fNeutralPionMesonCutArray; // array with neutral pion cuts
141  TList* fMesonCutArray; // array with neutral meson cuts
142  AliConvEventCuts* fEventCuts; // current event cuts
143  AliConversionPhotonCuts* fConversionCuts; // current conversion cuts
144  AliCaloPhotonCuts* fClusterCuts; // current cluster cuts
145 
146  // TTrees
147  TTree** fTreePiPiSameMother; // Tree containing info about the mother of two pions who have the same mother,
148  // if ID isn't covered by current implementations
149  TTree** fTreePiPiPiSameMother; // Tree containing info about the mother of three pions who have the same mother,
150  TTree** fTreeEventInfoOmega; // Tree containing information about an event where omega->pi+pi-pi0 was found
151  TTree** fTreeEventInfoEta; // Tree containing information about an event where eta->pi+pi-pi0 was found
152  // if ID isn't covered by current implementations
153  Short_t fCasePiPi; // 0:PiPlPiMi 1:PiMiPiZero 1:PiPlPiZero
154  Float_t fSamePiPiMotherID; // ID of mother of two pions
155  Float_t fSamePiPiMotherInvMass; // Invariant mass of mother of two pions
156  Float_t fSamePiPiMotherPt; // pT of mother of two pions
157  Float_t fSamePiPiPiMotherID; // ID of mother of two pions
158  Float_t fSamePiPiPiMotherInvMass; // Invariant mass of mother of two pions
159  Float_t fSamePiPiPiMotherPt; // pT of mother of two pions
160  Float_t fV0MultiplicityOmegaEvent; // V0 multiplicity of an event where a true omega was found
161  Float_t fTrackMultiplicityOmegaEvent; // track multiplicity of an event where a true omega was found
162  Float_t fZVertexOmegaEvent; // z position of primary vertex of an event where a true omega was found
163  Float_t fPtOmega; // pT of a true omega
164  Float_t fV0MultiplicityEtaEvent; // V0 multiplicity of an event where a true Eta was found
165  Float_t fTrackMultiplicityEtaEvent; // track multiplicity of an event where a true Eta was found
166  Float_t fZVertexEtaEvent; // z position of primary vertex of an event where a true Eta was found
167  Float_t fPtEta; // pT of a true Eta
168  Float_t fPDGMassPi0; // PDG mass of pi0
169  // reconstructed particles
170  TH1F** fHistoConvGammaPt; // array of histos of conversion photon, pt
171  TH1F** fHistoConvGammaEta; // array of histos of conversion photon, eta
172  TH1F** fHistoClusterGammaPt; // array of histos of Cluster photon, pt
173  TH1F** fHistoClusterGammaEta; // array of histos of Cluster photon, eta
174  TH1F** fHistoNegPionPt; // array of histos of negative pion, pt
175  TH1F** fHistoPosPionPt; // array of histos of positive pion, pt
176  TH1F** fHistoNegPionPhi; // array of histos of negative pion, phi
177  TH1F** fHistoPosPionPhi; // array of histos of positive pion, phi
178  TH1F** fHistoNegPionEta; // array of histos of negative pion, eta
179  TH1F** fHistoPosPionEta; // array of histos of positive pion, eta
180  TH2F** fHistoNegPionClsTPC; // array of histos of negative pion, findable tpc cluster, pT
181  TH2F** fHistoPosPionClsTPC; // array of histos of positive pion, findable tpc cluster, pT
182  TH2F** fHistoPionDCAxy; // array of histos of pion, dca_xy, pT
183  TH2F** fHistoPionDCAz; // array of histos of pion, dca_z, pT
184  TH2F** fHistoPionTPCdEdxNSigma; // array of histos of pion, p, TPC nSigma dEdx pion
185  TH2F** fHistoPionTPCdEdx; // array of histos of pion, p, TPC dEdx
186  TH2F** fHistoPionPionInvMassPt; // array of histos of pion pion, invMass, pT_{pi+pi-}
187  TH2F** fHistoGammaGammaInvMassPt; // array of histos of gamma-gamma, invMass, pT_{gamma gamma}
188  TH2F** fHistoMotherInvMassPt; // array of histos of pi+pi-pi0 same event, invMass, pT_{pi+pi-pi0}
189  TH2F** fHistoMotherInvMassPtRejectedKinematic; // array of histos of rejected pi+pi-pi0 same event, invMass, pT_{pi+pi-pi0}
190  TH2F** fHistoBackInvMassPtGroup1; // Event mixing background group 1 (pi+ and pi- from same event)
191  TH2F** fHistoBackInvMassPtGroup2; // Event mixing background group 2 (pi+ and pi0 from same event)
192  TH2F** fHistoBackInvMassPtGroup3; // Event mixing background group 3 (pi- and pi0 from same event)
193  TH2F** fHistoBackInvMassPtGroup4; // Event mixing background group 4 (no pion from same event)
194  TH2F** fHistoMotherLikeSignBackInvMassPt; // array of histos of pi+pi+pi0 likesign mixed event, invMass, pT_{pi+pi-pi0}
195 
196  // angle distributions
197  TH2F** fHistoAngleOmegaPiPlPiMi; // angle between combined Pi+ and Pi- and omega
198  TH2F** fHistoAngleOmegaPiZero; // angle between Pi0 and omega
199  TH2F** fHistoAngleOmegaPiPl; // angle between Pi+ and omega
200  TH2F** fHistoAngleOmegaPiMi; // angle between Pi- and omega
201  TH2F** fHistoAnglePiPlPiMi; // angle between Pi+ and Pi-
202  TH2F** fHistoAnglePiZeroPiMi; // angle between Pi0 and Pi-
203  TH2F** fHistoAnglePiPlPiZero; // angle between Pi+ and Pi0
204  TH2F** fHistoAngleSum; // angle between omega and Pi0 + angle between Pi+ and Pi- + angle between Pi0 and Pi-
205  TH2F** fHistoTrueAngleSum;
206 
207  TH2F** fHistoMotherInvMassSubPi0; // invariant mass of (pi+,pi-,pi0) - invariant mass of pi0
208  TH2F** fHistoBackInvMassPtGroup1SubPi0; // background group 1, invMass-invMass(pi0), pT_{pi+pi-pi0} (pi+ and pi- from same event)
209  TH2F** fHistoBackInvMassPtGroup2SubPi0; // background group 2, invMass-invMass(pi0), pT_{pi+pi-pi0} (pi+ and pi0 from same event)
210  TH2F** fHistoBackInvMassPtGroup3SubPi0; // background group 3, invMass-invMass(pi0), pT_{pi+pi-pi0} (pi+ and pi0 from same event)
211  TH2F** fHistoBackInvMassPtGroup4SubPi0; // background group 4, invMass-invMass(pi0), pT_{pi+pi-pi0} (no pion from same event)
212  TH2F** fHistoMotherLikeSignBackInvMassSubPi0Pt; // array of histos of pi+pi+pi0 likesign mixed event, invMass-invMass(pi0), pT_{pi+pi-pi0}
213 
214  TH2F** fHistoMotherInvMassFixedPzPi0; // invariant mass of (pi+,pi-,pi0) - invariant mass of pi0
215  TH2F** fHistoBackInvMassPtGroup1FixedPzPi0; // background group 1 mixed event, invMass, pT_{pi+pi-pi0}, the Pz of the pi0 was fixed such that
216  // its invMass matches the PDG value
217  TH2F** fHistoBackInvMassPtGroup2FixedPzPi0; // background group 2 mixed event, invMass, pT_{pi+pi-pi0}, the Pz of the pi0 was fixed such that
218  // its invMass matches the PDG value
219  TH2F** fHistoBackInvMassPtGroup3FixedPzPi0; // background group 3 mixed event, invMass, pT_{pi+pi-pi0}, the Pz of the pi0 was fixed such that
220  // its invMass matches the PDG value
221  TH2F** fHistoBackInvMassPtGroup4FixedPzPi0; // background group 4 mixed event, invMass, pT_{pi+pi-pi0}, the Pz of the pi0 was fixed such that
222  // its invMass matches the PDG value
223  TH2F** fHistoMotherLikeSignBackInvMassFixedPzPi0Pt; // array of histos of pi+pi+pi0 likesign mixed event, invMass, pT_{pi+pi+pi0}, the Pz of the pi0 was fixed such that
224  // its invMass matches the PDG value
225  // pure MC properties
226  TH1F** fHistoMCAllGammaPt; // array of histos of all produced gammas in the specified y range
227  TH1F** fHistoMCConvGammaPt; // array of histos of all converted gammas in the specified y range
228  TH1F** fHistoMCAllPosPionsPt; // array of histos with all produced primary positive pions in the specified y range
229  TH1F** fHistoMCAllNegPionsPt; // array of histos with all produced primary negative pions in the specified y range
230  TH1F** fHistoMCGammaFromNeutralMesonPt; // array of histos of all produced gammas from omega or eta via pi+pi-pi0 in the specified y range/
231  TH1F** fHistoMCPosPionsFromNeutralMesonPt; // array of histos of all produced positive pions from omega or eta via pi+pi-pi0 in the specified y range/
232  TH1F** fHistoMCNegPionsFromNeutralMesonPt; // array of histos of all produced negative pions from omega or eta via pi+pi-pi0 in the specified y range/
233  TH1F** fHistoMCEtaPiPlPiMiPiZeroPt; // array of histos of produced etas via pi+pi-pi0 in the specified y range
234  TH1F** fHistoMCEtaPiPlPiMiPiZeroInAccPt; // array of histos of produced etas via pi+pi-pi0 in the specified y range,
235  // with decay products in respective y, eta ranges
236  TH1F** fHistoMCOmegaPiPlPiMiPiZeroPt; // array of histos of produced omegas via pi+pi-pi0 in the specified y range
237  TH1F** fHistoMCOmegaPiPlPiMiPiZeroInAccPt; // array of histos of produced omegas via pi+pi-pi0 in the specified y range,
238  // with decay products in respective y, eta ranges
239 
240  // reconstructed particles MC validated
241  TH2F** fHistoTrueMotherPiPlPiMiPiZeroInvMassPt; // histos with reconstructed validated eta or omega, inv mass, pT
242  TH2F** fHistoTrueMotherOmegaPiPlPiMiPiZeroInvMassPt; // histos with reconstructed validated omega, inv mass, pT
243  TH2F** fHistoTrueMotherEtaPiPlPiMiPiZeroInvMassPt; // histos with reconstructed validated eta, inv mass, pT
244  TH2F** fHistoTrueMotherGammaGammaInvMassPt; // histos with reconstructed validated pi0, inv mass, pT
245  TH2F** fHistoTrueMotherGammaGammaFromEtaInvMassPt; // histos with reconstructed validated pi0, inv mass, pT
246  TH2F** fHistoTrueMotherGammaGammaFromOmegaInvMassPt; // histos with reconstructed validated pi0, inv mass, pT
247  TH1F** fHistoTrueConvGammaPt; // histos with reconstructed validated conv gamma, pT
248  TH1F** fHistoTrueConvGammaFromNeutralMesonPt; // histos with reconstructed validated conv gamma from eta or omega via pi0, pT
249  TH1F** fHistoTrueClusterGammaPt; // histos with reconstructed validated cluster gamma, pT
250  TH1F** fHistoTrueClusterGammaFromNeutralMesonPt; // histos with reconstructed validated cluster gamma from eta or omega via pi0, pT
251  TH1F** fHistoTruePosPionPt; // histos with reconstructed validated positive pion, pT
252  TH1F** fHistoTruePosPionFromNeutralMesonPt; // histos with reconstructed validated positive pion from eta or omega, pT
253  TH1F** fHistoTrueNegPionPt; // histos with reconstructed validated negative pion, pT
254  TH1F** fHistoTrueNegPionFromNeutralMesonPt; // histos with reconstructed validated negative pion from eta or omega, pT
255  TH2F** fHistoTruePionPionInvMassPt; // histos with reconstructed validated two pion, invariant mass, pT
256  TH2F** fHistoTruePionPionFromSameMotherInvMassPt; // histos with reconstructed validated two pion from same mother, invariant mass, pT
257  TH2F** fHistoTruePionPionFromEtaInvMassPt; // histos with reconstructed validated two pion from eta , invariant mass, pT
258  TH2F** fHistoTruePionPionFromOmegaInvMassPt; // histos with reconstructed validated two pion from omega, invariant mass, pT
259 
260  TH2F** fHistoTruePiPlPiMiSameMotherFromEtaInvMassPt; // histos with reconstructed validated pi+ pi- from omega, invariant mass, pT
261  TH2F** fHistoTruePiPlPiMiSameMotherFromOmegaInvMassPt; // histos with reconstructed validated pi+ pi- from eta, invariant mass, pT
262  TH2F** fHistoTruePiPlPiMiSameMotherFromRhoInvMassPt; // histos with reconstructed validated pi+ pi- from rho0, invariant mass, pT
263  TH2F** fHistoTruePiPlPiMiSameMotherFromEtaPrimeInvMassPt; // histos with reconstructed validated pi+ pi- from etaprime, invariant mass, pT
264  TH2F** fHistoTruePiPlPiMiSameMotherFromK0sInvMassPt; // histos with reconstructed validated pi+ pi- from K0s, invariant mass, pT
265  TH2F** fHistoTruePiPlPiMiSameMotherFromK0lInvMassPt; // histos with reconstructed validated pi+ pi- from K0s, invariant mass, pT
266 
267  TH2F** fHistoTruePiMiPiZeroSameMotherFromEtaInvMassPt; // histos with reconstructed validated pi0 pi- from omega, invariant mass, pT
268  TH2F** fHistoTruePiMiPiZeroSameMotherFromOmegaInvMassPt; // histos with reconstructed validated pi0 pi- from eta, invariant mass, pT
269  TH2F** fHistoTruePiMiPiZeroSameMotherFromRhoInvMassPt; // histos with reconstructed validated pi0 pi- from rho0, invariant mass, pT
270  TH2F** fHistoTruePiMiPiZeroSameMotherFromK0lInvMassPt; // histos with reconstructed validated pi0 pi- from rho0, invariant mass, pT
271 
272  TH2F** fHistoTruePiPlPiZeroSameMotherFromEtaInvMassPt; // histos with reconstructed validated pi0 pi+ from omega, invariant mass, pT
273  TH2F** fHistoTruePiPlPiZeroSameMotherFromOmegaInvMassPt; // histos with reconstructed validated pi0 pi+ from eta, invariant mass, pT
274  TH2F** fHistoTruePiPlPiZeroSameMotherFromRhoInvMassPt; // histos with reconstructed validated pi0 pi+ from rho0, invariant mass, pT
275  TH2F** fHistoTruePiPlPiZeroSameMotherFromK0lInvMassPt; // histos with reconstructed validated pi0 pi+ from K0l, invariant mass, pT
276  TH2F** fHistoTruePiPlPiMiPiZeroPureCombinatoricalInvMassPt; // histos with reconstructed validated pi+pi-pi0 that are pure combinatorical (do not share a mother)
277  TH2F** fHistoTruePiPlPiMiPiZeroContaminationInvMassPt; // histos with reconstructed pi+pi-pi0 that are not actually pions
278 
279  TH2F** fHistoDoubleCountTruePi0InvMassPt;
280  TH2F** fHistoDoubleCountTrueEtaInvMassPt;
281  TH2F** fHistoDoubleCountTrueOmegaInvMassPt;
282  TH2F** fHistoDoubleCountTrueConvGammaRPt;
283  vector<Int_t> fVectorDoubleCountTruePi0s;
284  vector<Int_t> fVectorDoubleCountTrueEtas;
285  vector<Int_t> fVectorDoubleCountTrueOmegas;
286  vector<Int_t> fVectorDoubleCountTrueConvGammas;
287  // Event properties
288  TH1I** fHistoNEvents; // histo for event counting
289  TH1I** fHistoNGoodESDTracks; // histo number of reconstructed primary tracks
290  TProfile** fProfileEtaShift; // profile for eta shift bookkeeping
291  TH2F** fHistoSPDClusterTrackletBackground;
292 
293 
294  TRandom3 fRandom; // random number
295  Int_t fnCuts; // number of cuts to be run in parallel
296  Int_t fiCut; // current cut
297  Int_t fNumberOfESDTracks; // integer with number of primary tracks in this event
298  Bool_t fMoveParticleAccordingToVertex; // Flag to move parice to the vertex
299  Int_t fIsHeavyIon; // Flag for collision system 0: pp, 1: PbPb, 2: pPb
300  Bool_t fDoMesonAnalysis; // Flag for switching on meson analysis
301  Int_t fDoMesonQA; // Switching for meson QA 0: no QA 1: small QA 2: big QA
302  Bool_t fIsFromMBHeader; // Flag for particle whether it belongs to accepted header
303  Bool_t fIsMC; // Flag for MC
304  Bool_t fDoLightOutput; // Flag to turn on light output
305  Int_t fNeutralPionMode; // Flag how neutral pion is reconstructed 0=PCM-PCM, 1=PCM-Calo, 2=Calo-Calo
306  Double_t fTolerance; // tolerance in rad for angle cuts
307 private:
308  AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero( const AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero& ); // Not implemented
309  AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero& operator=( const AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero& ); // Not implemented
310 
311  ClassDef(AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero, 22);
312 };
313 
314 #endif // ALIANALYSISTASKNEUTRALMESONTOPIPLPIMIPIZERO_H
315 
AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero::operator=
AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero & operator=(const AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero &)
AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero::ProcessTrueNeutralPionCandidatesPureCalo
void ProcessTrueNeutralPionCandidatesPureCalo(AliAODConversionMother *Pi0Candidate, AliAODConversionPhoton *TrueGammaCandidate0, AliAODConversionPhoton *TrueGammaCandidate1)
Definition: AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero.cxx:1927
Double_t
double Double_t
Definition: External.C:58
AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero::fVectorDoubleCountTruePi0s
vector< Int_t > fVectorDoubleCountTruePi0s
array of histos with double counted photons, R, pT
Definition: AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero.h:283
TH2F
Definition: External.C:236
AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero::fVectorDoubleCountTrueOmegas
vector< Int_t > fVectorDoubleCountTrueOmegas
vector containing labels of validated eta
Definition: AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero.h:285
AliConversionMesonCuts.h
AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero::fVectorDoubleCountTrueEtas
vector< Int_t > fVectorDoubleCountTrueEtas
vector containing labels of validated pi0
Definition: AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero.h:284
AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero::ProcessTrueMesonCandidates
void ProcessTrueMesonCandidates(AliAODConversionMother *Pi0Candidate, AliAODConversionMother *TrueNeutralPionCandidate, AliAODConversionPhoton *TrueVirtualGammaCandidate)
Definition: AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero.cxx:3250
AliPrimaryPionSelector.h
AliCaloPhotonCuts.h
Int_t
int Int_t
Definition: External.C:63
AliCaloPhotonCuts
Class handling all kinds of selection cuts for Gamma Calo analysis.
Definition: AliCaloPhotonCuts.h:76
TH1I
Definition: External.C:204
AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero::ProcessTrueCaloPhotonCandidates
void ProcessTrueCaloPhotonCandidates(AliAODConversionPhoton *TruePhotonCandidate)
Definition: AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero.cxx:1616
Float_t
float Float_t
Definition: External.C:68
AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero::fRandom
TRandom3 fRandom
array of histos with SPD tracklets vs SPD clusters for background rejection
Definition: AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero.h:294
AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero::fHistoNEvents
TH1I ** fHistoNEvents
vector containing labels of validated photons
Definition: AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero.h:288
AliGammaConversionAODBGHandler.h
AliConversionPhotonCuts
Class handling all kinds of selection cuts for Gamma Conversion analysis.
Definition: AliConversionPhotonCuts.h:76
AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero::fHistoDoubleCountTrueConvGammaRPt
TH2F ** fHistoDoubleCountTrueConvGammaRPt
array of histos with double counted omegas, invMass, pT
Definition: AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero.h:282
mode
Int_t mode
Definition: anaM.C:41
Short_t
short Short_t
Definition: External.C:23
AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero::ProcessTrueNeutralPionCandidatesPureConversions
void ProcessTrueNeutralPionCandidatesPureConversions(AliAODConversionMother *Pi0Candidate, AliAODConversionPhoton *TrueGammaCandidate0, AliAODConversionPhoton *TrueGammaCandidate1)
Definition: AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero.cxx:1999
isMC
Bool_t isMC
Definition: ReadDvsMultiplicity.C:60
AliConvEventCuts.h
AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero::fVectorDoubleCountTrueConvGammas
vector< Int_t > fVectorDoubleCountTrueConvGammas
vector containing labels of validated omega
Definition: AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero.h:286
AliV0ReaderV1.h
AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero::fHistoDoubleCountTrueEtaInvMassPt
TH2F ** fHistoDoubleCountTrueEtaInvMassPt
array of histos with double counted pi0s, invMass, pT
Definition: AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero.h:280
AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero::fHistoDoubleCountTrueOmegaInvMassPt
TH2F ** fHistoDoubleCountTrueOmegaInvMassPt
array of histos with double counted etas, invMass, pT
Definition: AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero.h:281
AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero::KinematicCut
Bool_t KinematicCut(AliAODConversionMother *negpion, AliAODConversionMother *pospion, AliAODConversionMother *neutpion, AliAODConversionMother *omega)
Definition: AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero.cxx:3221
AliConvEventCuts
Class handling all kinds of selection cuts for Gamma Conversion analysis.
Definition: AliConvEventCuts.h:55
Option_t
const char Option_t
Definition: External.C:48
AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero::MoveParticleAccordingToVertex
void MoveParticleAccordingToVertex(AliAODConversionMother *particle, const AliGammaConversionAODBGHandler::GammaConversionVertex *vertex)
Definition: AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero.cxx:3497
Bool_t
bool Bool_t
Definition: External.C:53
AliKFConversionPhoton.h
AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero::ProcessTrueNeutralPionCandidatesMixedConvCalo
void ProcessTrueNeutralPionCandidatesMixedConvCalo(AliAODConversionMother *Pi0Candidate, AliAODConversionPhoton *TrueGammaCandidate0, AliAODConversionPhoton *TrueGammaCandidate1)
Definition: AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero.cxx:2252
AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero::CheckVectorForDoubleCount
Bool_t CheckVectorForDoubleCount(vector< Int_t > &vec, Int_t tobechecked)
Definition: AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero.cxx:3612
AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero::ProcessTrueNeutralPionCandidatesPureConversionsAOD
void ProcessTrueNeutralPionCandidatesPureConversionsAOD(AliAODConversionMother *Pi0Candidate, AliAODConversionPhoton *TrueGammaCandidate0, AliAODConversionPhoton *TrueGammaCandidate1)
Definition: AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero.cxx:2085