1 #ifndef ALIANALYSISTASKGAMMACONVERSION_H
2 #define ALIANALYSISTASKGAMMACONVERSION_H
4 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
5 * See cxx source for full Copyright notice */
7 ////////////////////////////////////////////////
8 //---------------------------------------------
9 // Class used to do analysis on conversion pairs
10 //---------------------------------------------
11 ////////////////////////////////////////////////
13 #include "AliAnalysisTaskSE.h"
15 #include "AliV0Reader.h"
16 //#include "AliCFManager.h" // for CF
17 //#include "AliCFContainer.h" // for CF
19 class AliGammaConversionAODObject;
21 class AliGammaConversionHistograms;
26 class AliESDInputHandler;
32 class AliESDtrackCuts;
33 class AliCFManager; // for CF
34 class AliCFContainer; // for CF
36 class AliAnalysisTaskGammaConversion : public AliAnalysisTaskSE
40 AliAnalysisTaskGammaConversion();
41 AliAnalysisTaskGammaConversion(const char* name);
42 virtual ~AliAnalysisTaskGammaConversion() ;// virtual destructor
44 // Implementation of interface methods
45 virtual void UserCreateOutputObjects();
47 virtual void LocalInit() {Init();}
48 virtual void UserExec(Option_t *option);
49 virtual void Terminate(Option_t *option);
50 //virtual void ConnectInputData(Option_t * option);
53 void ProcessV0sNoCut();
55 void ProcessGammasForNeutralMesonAnalysis();
56 void ProcessGammasForOmegaMesonAnalysis();
57 void ProcessConvPHOSGammasForNeutralMesonAnalysis();
58 void RecalculateV0ForGamma();
60 void SetCFManager(AliCFManager * const io) {fCFManager = io;};
61 AliCFManager *GetCFManager() const {return fCFManager;}
65 TString GetAODBranchName() const {return fAODBranchName;}
66 void SetAODBranchName(TString name) {fAODBranchName = name ;}
67 void SetForceAOD(Bool_t forceAOD ) { fKFForceAOD = forceAOD; }
68 void FillAODWithConversionGammas();
71 static Bool_t IsGoodImpPar(AliESDtrack *const track);
73 // for GammaJetAnalysis
74 void ProcessGammasForGammaJetAnalysis();
75 void CreateListOfChargedParticles();
76 Double_t GetMinimumDistanceToCharge(Int_t indexHighestPtGamma);
77 void CalculateJetCone(Int_t gammaIndex);
78 Int_t GetIndexHighestPtGamma();
79 void SetESDtrackCuts();
82 void SetMinPtForGammaJet(Double_t minPtForGammaJet){fMinPtForGammaJet=minPtForGammaJet;}
83 void SetMinIsoConeSize(Double_t minIsoConeSize){fMinIsoConeSize=minIsoConeSize;}
84 void SetMinPtIsoCone(Double_t minPtIsoCone){fMinPtIsoCone=minPtIsoCone;}
85 void SetMinPtGamChargedCorr(Double_t minPtGamChargedCorr){fMinPtGamChargedCorr=minPtGamChargedCorr;}
86 void SetMinPtJetCone(Double_t minPtJetCone){fMinPtJetCone=minPtJetCone;}
88 void SetLowPtMapping(Double_t lowPtMapping){fLowPtMapping=lowPtMapping;}
89 void SetHighPtMapping(Double_t highPtMapping){fHighPtMapping=highPtMapping;}
92 void SetHistograms(AliGammaConversionHistograms *const histograms){fHistograms=histograms;}
93 void SetTriggerCINT1B(Bool_t flag){fTriggerCINT1B=flag;}
94 void SetDoMCTruth(Bool_t flag){fDoMCTruth=flag;}
95 void SetDoNeutralMeson(Bool_t flag){fDoNeutralMeson=flag;}
96 void SetDoOmegaMeson(Bool_t flag){fDoOmegaMeson=flag;}
97 void SetDoNeutralMesonV0MCCheck(Bool_t flag){fDoNeutralMesonV0MCCheck=flag;}
98 void SetDoJet(Bool_t flag){fDoJet=flag;}
99 void SetDoChic(Bool_t flag){fDoChic=flag;}
100 void SetRecalculateV0ForGamma(Bool_t flag){fRecalculateV0ForGamma=flag;}
102 void SetElectronMass(Double_t electronMass){fElectronMass = electronMass;}
103 void SetGammaMass(Double_t gammaMass){fGammaMass = gammaMass;}
104 void SetGammaWidth(Double_t gammaWidth){fGammaWidth = gammaWidth;}
105 void SetPi0Mass(Double_t pi0Mass){fPi0Mass = pi0Mass;}
106 void SetPi0Width(Double_t pi0Width){fPi0Width = pi0Width;}
107 void SetEtaMass(Double_t etaMass){fEtaMass = etaMass;}
108 void SetEtaWidth(Double_t etaWidth){fEtaWidth = etaWidth;}
109 void SetMinOpeningAngleGhostCut(Double_t ghostCut){fMinOpeningAngleGhostCut = ghostCut;}
110 void SetV0Reader(AliV0Reader* const reader){fV0Reader=reader; fV0Reader->SetESDtrackCuts(fEsdTrackCuts);}
111 void SetCalculateBackground(Bool_t bg){fCalculateBackground=bg;}
112 void CalculateBackground();
113 void SetWriteNtuple(Bool_t writeNtuple){fWriteNtuple = writeNtuple;}
115 Double_t GetMCOpeningAngle(TParticle* const daughter0, TParticle* const daughter1) const;
116 void CheckV0Efficiency();
117 void SetDeltaAODFileName(TString fn) { fKFDeltaAODFileName = fn; };
119 //////////////////Chi_c Analysis////////////////////////////
120 void GetPID(AliESDtrack *track, Stat_t &pid, Stat_t &weight);
121 double GetSigmaToVertex(AliESDtrack* t);
122 void ElectronBackground(TString hBg, TClonesArray e);
123 void FillAngle(TString histoName,TClonesArray const tlVeNeg, TClonesArray const tlVePos);
124 void FillElectronInvMass(TString histoName, TClonesArray const negativeElectron, TClonesArray const positiveElectron);
125 void FillGammaElectronInvMass(TString histoMass,TString histoDiff, TClonesArray const fKFGammas, TClonesArray const tlVeNeg,TClonesArray const tlVePos);
126 void CleanWithAngleCuts(TClonesArray const negativeElectrons, TClonesArray const positiveElectrons, TClonesArray const gammas);
127 TClonesArray GetTLorentzVector(TClonesArray* esdTrack);
128 void ProcessGammaElectronsForChicAnalysis();
129 ///////////////////////////////////////////////////////////////
131 void SetDoCF(Bool_t flag){fDoCF = flag;}
134 AliAnalysisTaskGammaConversion(const AliAnalysisTaskGammaConversion&); // Not implemented
135 AliAnalysisTaskGammaConversion& operator=(const AliAnalysisTaskGammaConversion&); // Not implemented
137 /// Add reconstructed pions to aod
138 void AddPionToAOD(AliKFParticle * pionkf, Double_t mass, Int_t daughter1, Int_t daughter2);
139 void AddOmegaToAOD(AliKFParticle * pionkf, Double_t mass, Int_t daughter1, Int_t daughter2);
144 kStepReconstructable = 1,
149 kStepdEdx_electronselection = 6,
150 kStepdEdx_pionrejection = 7,
151 kStepNContributors = 8,
165 AliV0Reader* fV0Reader; // The V0 reader object
167 AliStack * fStack; // pointer to the MC particle stack
168 AliMCEventHandler *fMCTruth; // for CF pointer to MCTruth
169 AliMCEvent *fGCMCEvent; // for CF pointer to the MC Event
170 AliESDEvent* fESDEvent; //pointer to the ESDEvent
171 TList * fOutputContainer; // Histogram container
172 AliCFManager *fCFManager; // for CF
173 // AliCFContainer *container; // for CF
177 AliGammaConversionHistograms *fHistograms; // Pointer to the histogram handling class
178 Bool_t fTriggerCINT1B; //Flag to select trigger CINT1B
179 Bool_t fDoMCTruth; // Flag to switch on/off MC truth
180 Bool_t fDoNeutralMeson; // flag
181 Bool_t fDoOmegaMeson; // flag
182 Bool_t fDoJet; // flag
183 Bool_t fDoChic; // flag
184 Bool_t fRecalculateV0ForGamma;//flag
186 TClonesArray * fKFReconstructedGammasTClone; //! transient
187 TClonesArray * fKFReconstructedPi0sTClone; //! transient
188 TClonesArray * fKFRecalculatedGammasTClone; //! transient
189 TClonesArray * fCurrentEventPosElectronTClone; //! transient
190 TClonesArray * fCurrentEventNegElectronTClone; //! transient
191 TClonesArray * fKFReconstructedGammasCutTClone; //! transient
192 TClonesArray * fPreviousEventTLVNegElectronTClone; //! transient
193 TClonesArray * fPreviousEventTLVPosElectronTClone; //! transient
195 // vector<AliKFParticle> fKFReconstructedGammas; // vector containing all reconstructed gammas
196 vector<Int_t> fElectronv1; // vector containing index of electron 1
197 vector<Int_t> fElectronv2; // vector containing index of electron 2
199 vector<Int_t> fGammav1; // vector containing index of gamma 1
200 vector<Int_t> fGammav2; // vector containing index of gamma 2
202 vector<Int_t> fElectronRecalculatedv1; // vector containing index of electron 1
203 vector<Int_t> fElectronRecalculatedv2; // vector containing index of electron 2
205 // AliESDpid * fESDpid; // esd pid
208 ///////Chi_c Analysis///////////////////////////
209 // vector<AliESDtrack*> fCurrentEventPosElectron; // comment here
210 // vector<AliESDtrack*> fCurrentEventNegElectron; // comment here
211 // vector<AliKFParticle> fKFReconstructedGammasCut; // comment here
212 // vector<TLorentzVector> fPreviousEventTLVNegElectron; // comment here
213 // vector<TLorentzVector> fPreviousEventTLVPosElectron; // comment here
214 //////////////////////////////////////////////////
217 Double_t fElectronMass; //electron mass
218 Double_t fGammaMass; //gamma mass
219 Double_t fPi0Mass; //pi0mass
220 Double_t fEtaMass; //eta mass
223 Double_t fGammaWidth; //gamma width cut
224 Double_t fPi0Width; // pi0 width cut
225 Double_t fEtaWidth; // eta width cut
227 Double_t fMinOpeningAngleGhostCut; // minimum angle cut
229 AliESDtrackCuts* fEsdTrackCuts; // Object containing the parameters of the esd track cuts
231 Bool_t fCalculateBackground; //flag to set backgrount calculation on/off
232 Bool_t fWriteNtuple; // flag to set if writing to ntuple on/off
233 TNtuple *fGammaNtuple; // Ntuple for gamma values
234 TNtuple *fNeutralMesonNtuple;// NTuple for mesons
236 Int_t fTotalNumberOfAddedNtupleEntries; // number of added ntuple entries
238 TClonesArray* fChargedParticles; //! transient
239 vector<Int_t> fChargedParticlesId; //! transient
241 Double_t fGammaPtHighest; //! transient
242 Double_t fMinPtForGammaJet; //! transient
243 Double_t fMinIsoConeSize; //! transient
244 Double_t fMinPtIsoCone; //! transient
245 Double_t fMinPtGamChargedCorr; //! transient
246 Double_t fMinPtJetCone; //! transient
247 Int_t fLeadingChargedIndex; //! transient
248 Double_t fLowPtMapping; //! transient
249 Double_t fHighPtMapping; //! transient
250 Bool_t fDoCF; //! transient
252 TClonesArray * fAODGamma; //TClonesArray for gammas to put in AOD
253 TClonesArray * fAODPi0; //TTClonesArray for Pi0s to put in AOD
254 TClonesArray * fAODOmega; //TTClonesArray for omegas to put in AOD
255 TString fAODBranchName; // New AOD branch name
256 Bool_t fKFForceAOD; //Set the Analysis Manager FillAOD variable to true every event
257 TString fKFDeltaAODFileName; //! File name for delta AOD (if any)
258 Bool_t fDoNeutralMesonV0MCCheck; //flag
260 vector<Int_t>fKFReconstructedGammasV0Index; // index of the reconstructed v0s
262 ClassDef(AliAnalysisTaskGammaConversion, 9); // Analysis task for gamma conversions
265 #endif //ALIANALYSISTASKGAMMA_H