1 #ifndef ALICALOPHOTONCUTS_H
2 #define ALICALOPHOTONCUTS_H
4 // Class handling all kinds of selection cuts for Gamma Conversion analysis
5 // Authors: Svein Lindal, Daniel Lohner *
7 #include "AliConversionPhotonBase.h"
8 #include "AliAODConversionMother.h"
9 #include "AliAODTrack.h"
10 #include "AliESDtrack.h"
11 #include "AliVTrack.h"
12 #include "AliVCluster.h"
13 #include "AliAODTrack.h"
15 #include "AliAnalysisCuts.h"
18 #include "AliAnalysisUtils.h"
19 #include "AliAnalysisManager.h"
23 class AliConversionPhotonBase;
28 class AliAnalysisCuts;
31 class AliAnalysisManager;
32 class AliAODMCParticle;
36 class AliCaloPhotonCuts : public AliAnalysisCuts {
45 kDistanceToBadChannel,
68 //handeling of CutString
69 static const char * fgkCutNames[kNCuts];
70 Bool_t SetCutIds(TString cutString);
72 Bool_t SetCut(cutIds cutID, Int_t cut);
73 Bool_t UpdateCutString();
75 void PrintCutsWithValues();
77 Bool_t InitializeCutsFromCutString(const TString analysisCutSelection);
78 TString GetCutNumber();
79 Int_t GetClusterType() {return fClusterType; };
82 AliCaloPhotonCuts(const char *name="ClusterCuts", const char * title="Cluster Cuts");
83 AliCaloPhotonCuts(const AliCaloPhotonCuts&);
84 AliCaloPhotonCuts& operator=(const AliCaloPhotonCuts&);
87 virtual ~AliCaloPhotonCuts();
89 virtual Bool_t IsSelected(TObject* /*obj*/) {return kTRUE;}
90 virtual Bool_t IsSelected(TList* /*list*/) {return kTRUE;}
92 Bool_t ClusterIsSelected(AliVCluster* cluster, AliVEvent *event, Bool_t isMC);
93 Bool_t ClusterIsSelectedMC(TParticle *particle,AliStack *fMCStack);
94 Bool_t ClusterIsSelectedAODMC(AliAODMCParticle *particle,TClonesArray *aodmcArray);
96 void InitCutHistograms(TString name="");
97 void SetFillCutHistograms(TString name="") {if(!fHistograms){InitCutHistograms(name);}}
98 TList* GetCutHistograms() {return fHistograms;}
99 void FillClusterCutIndex(Int_t photoncut) {if(fHistCutIndex)fHistCutIndex->Fill(photoncut);}
101 void SetExtendedMatching(Bool_t extendedMatching) {fExtendedMatching = extendedMatching;}
104 Bool_t AcceptanceCuts(AliVCluster* cluster, AliVEvent *event);
105 Bool_t ClusterQualityCuts(AliVCluster* cluster,AliVEvent *event, Bool_t isMC);
107 Bool_t MatchConvPhotonToCluster(AliAODConversionPhoton* convPhoton, AliVCluster* cluster, AliVEvent* event );
109 // Set Individual Cuts
110 Bool_t SetClusterTypeCut(Int_t);
111 Bool_t SetMinEtaCut(Int_t);
112 Bool_t SetMaxEtaCut(Int_t);
113 Bool_t SetMinPhiCut(Int_t);
114 Bool_t SetMaxPhiCut(Int_t);
115 Bool_t SetDistanceToBadChannelCut(Int_t);
116 Bool_t SetTimingCut(Int_t);
117 Bool_t SetTrackMatchingCut(Int_t);
118 Bool_t SetExoticCellCut(Int_t);
119 Bool_t SetMinEnergyCut(Int_t);
120 Bool_t SetMinNCellsCut(Int_t);
121 Bool_t SetMaxM02(Int_t);
122 Bool_t SetMinM02(Int_t);
123 Bool_t SetMaxM20(Int_t);
124 Bool_t SetMinM20(Int_t);
125 Bool_t SetDispersion(Int_t);
126 Bool_t SetNLM(Int_t);
132 Int_t fClusterType; // which cluster do we have
133 Double_t fMinEtaCut; // min eta cut
134 Double_t fMaxEtaCut; // max eta cut
135 Bool_t fUseEtaCut; // flag for switching on eta cut
136 Double_t fMinPhiCut; // phi cut
137 Double_t fMaxPhiCut; // phi cut
138 Bool_t fUsePhiCut; // flag for switching on phi cut
139 Double_t fMinDistanceToBadChannel; // minimum distance to bad channel
140 Bool_t fUseDistanceToBadChannel; // flag for switching on distance to bad channel cut
141 Double_t fMaxTimeDiff; // maximum time difference to triggered collision
142 Bool_t fUseTimeDiff; // flag for switching on time difference cut
143 Double_t fMinDistTrackToCluster; // minimum distance between track and cluster
144 Bool_t fUseDistTrackToCluster; // flag for switching on distance between track and cluster cut
145 Bool_t fExtendedMatching; // flag for switching on extended matching histograms
146 Double_t fExoticCell; // exotic cell cut
147 Bool_t fUseExoticCell; // flag for switching on exotic cell cut
148 Double_t fMinEnergy; // minium energy per cluster
149 Bool_t fUseMinEnergy; // flag for switching on minimum energy cut
150 Int_t fMinNCells; // minimum number of cells
151 Bool_t fUseNCells; // flag for switching on minimum N Cells cut
152 Double_t fMaxM02; // maximum M02
153 Double_t fMinM02; // minimum M02
154 Bool_t fUseM02; // flag for switching on M02 cut
155 Double_t fMaxM20; // maximum M20
156 Double_t fMinM20; // minimum M20
157 Bool_t fUseM20; // flag for switching on M20 cut
158 Double_t fMaxDispersion; // maximum dispersion
159 Bool_t fUseDispersion; // flag for switching on dispersion cut
160 Int_t fMinNLM; // minimum number of local maxima in cluster
161 Int_t fMaxNLM; // maximum number of local maxima in cluster
162 Bool_t fUseNLM; // flag for switching on NLM cut
165 TObjString* fCutString; // cut number used for analysis
168 TH1F* fHistCutIndex; // bookkeeping for cuts
169 TH1F* fHistAcceptanceCuts; // bookkeeping for acceptance cuts
170 TH1F* fHistClusterIdentificationCuts; // bookkeeping for cluster identification cuts
172 TH2F* fHistClusterEtavsPhiBeforeAcc; // eta-phi-distribution before acceptance cuts
173 TH2F* fHistClusterEtavsPhiAfterAcc; // eta-phi-distribution of all after acceptance cuts
174 TH2F* fHistClusterEtavsPhiAfterQA; // eta-phi-distribution of all after cluster quality cuts
175 // TH1F* fHistDistanceToBadChannelBeforeAcc; // distance to bad channel before acceptance cuts
176 // TH1F* fHistDistanceToBadChannelAfterAcc; // distance to bad channel after acceptance cuts
177 TH2F* fHistClusterTimevsEBeforeQA; // Cluster time vs E before cluster quality cuts
178 TH2F* fHistClusterTimevsEAfterQA; // Cluster time vs E after cluster quality cuts
179 // TH2F* fHistExoticCellBeforeQA; // Exotic cell: 1-Ecross/E cell vs Ecluster before acceptance cuts
180 // TH2F* fHistExoticCellAfterQA; // Exotic cell: 1-Ecross/E cell vs Ecluster after cluster quality cuts
181 // TH1F* fHistNMatchedTracks; // number of matched tracks
182 TH1F* fHistEnergyOfClusterBeforeQA; // enery per cluster before acceptance cuts
183 TH1F* fHistEnergyOfClusterAfterQA; // enery per cluster after cluster quality cuts
184 TH1F* fHistNCellsBeforeQA; // number of cells per cluster before acceptance cuts
185 TH1F* fHistNCellsAfterQA; // number of cells per cluster after cluster quality cuts
186 TH1F* fHistM02BeforeQA; // M02 before acceptance cuts
187 TH1F* fHistM02AfterQA; // M02 after cluster quality cuts
188 TH1F* fHistM20BeforeQA; // M20 before acceptance cuts
189 TH1F* fHistM20AfterQA; // M20 after cluster quality cuts
190 TH1F* fHistDispersionBeforeQA; // dispersion before acceptance cuts
191 TH1F* fHistDispersionAfterQA; // dispersion after cluster quality cuts
192 // TH1F* fHistNLMBeforeQA; // number of local maxima in cluster before acceptance cuts
193 // TH1F* fHistNLMAfterQA; // number of local maxima in cluster after cluster quality cuts
195 //Track matching histograms
196 TH1F* fHistClusterRBeforeQA; // cluster position in R=SQRT(x^2+y^2) (before QA)
197 TH1F* fHistClusterRAfterQA; // cluster position in R=SQRT(x^2+y^2) for matched tracks (After QA)
198 TH2F* fHistClusterdEtadPhiBeforeQA; // 2-dim plot dEta vs. dPhi
199 TH2F* fHistClusterdEtadPhiAfterQA; // 2-dim plot dEta vs. dPhi for matched tracks (after QA)
200 TH1F* fHistDistanceTrackToClusterBeforeQA; // distance cluster to track before acceptance cuts
201 TH1F* fHistDistanceTrackToClusterAfterQA; // distance cluster to track after cluster quality cuts
203 //Extended track matching histograms
204 TH2F* fHistClusterdEtadPhiPosTracksBeforeQA; // 2-dim plot dEta vs. dPhi
205 TH2F* fHistClusterdEtadPhiNegTracksBeforeQA; // 2-dim plot dEta vs. dPhi
206 TH2F* fHistClusterdEtadPtBeforeQA; // 2-dim plot dEta vs. Pt
207 TH2F* fHistClusterdPhidPtBeforeQA; // 2-dim plot dEta vs. Pt
208 TH2F* fHistClusterM20Pt_dPhiBeforeQA; // 2-dim plot M20 vs. Pt for given dPhi>0.05
209 TH2F* fHistClusterM02Pt_dPhiBeforeQA; // 2-dim plot M02 vs. Pt for given dPhi>0.05
210 TH2F* fHistClusterM20M02BeforeQA; // 2-dim plot M20 vs. M02
211 TH2F* fHistClusterM20M02AfterQA; // 2-dim plot M20 vs. M20
215 ClassDef(AliCaloPhotonCuts,2)