AliPhysics  608b256 (608b256)
AliAnalysisTaskGammaConvCalo.h
Go to the documentation of this file.
1 #ifndef ALIANLYSISTASKGAMMACONVCALO_cxx
2 #define ALIANLYSISTASKGAMMACONVCALO_cxx
3 
4 #include "AliAnalysisTaskSE.h"
5 #include "AliESDtrack.h"
6 #include "AliV0ReaderV1.h"
10 #include "AliCaloPhotonCuts.h"
11 #include "AliConvEventCuts.h"
13 #include "AliConversionMesonCuts.h"
14 #include "AliAnalysisManager.h"
15 #include "AliAnalysisTaskConvJet.h"
16 #include "TProfile2D.h"
17 #include "TH3.h"
18 #include "TH3F.h"
19 #include "THnSparse.h"
20 #include <vector>
21 #include <map>
22 
24  public:
25 
27  AliAnalysisTaskGammaConvCalo(const char *name);
29 
30  virtual void UserCreateOutputObjects();
31  virtual Bool_t Notify();
32  virtual void UserExec(Option_t *);
33  virtual void Terminate(const Option_t*);
34  void InitBack();
35 
36  void SetV0ReaderName(TString name){fV0ReaderName=name; return;}
37  void SetIsHeavyIon(Int_t flag){
38  fIsHeavyIon = flag;
39  }
40  // Function to set correction task setting
42 
43  // base functions for selecting photon and meson candidates in reconstructed data
44  void ProcessClusters();
47 
48  // MC functions
49  void SetIsMC ( Int_t isMC) { fIsMC = isMC ;}
50  void ProcessMCParticles ();
51  void ProcessAODMCParticles ();
53  void ProcessTruePhotonCandidates ( AliAODConversionPhoton* TruePhotonCandidate);
54  void ProcessTrueClusterCandidates ( AliAODConversionPhoton* TruePhotonCandidate,
55  Float_t clusM02);
57  Float_t clusM02);
58  void ProcessTruePhotonCandidatesAOD ( AliAODConversionPhoton* TruePhotonCandidate);
60  AliAODConversionPhoton *TrueGammaCandidate0,
61  AliAODConversionPhoton *TrueGammaCandidate1,
62  Bool_t matched);
64  AliAODConversionPhoton *TrueGammaCandidate0,
65  AliAODConversionPhoton *TrueGammaCandidate1,
66  Bool_t matched);
69  void ProcessJets();
70 
71  // switches for additional analysis streams or outputs
73  void SetLightOutput ( Bool_t flag ) { fDoLightOutput = flag ;}
74  void SetDoMesonAnalysis ( Bool_t flag ) { fDoMesonAnalysis = flag ;}
75  void SetDoMesonQA ( Int_t flag ) { fDoMesonQA = flag ;}
76  void SetDoPhotonQA ( Int_t flag ) { fDoPhotonQA = flag ;}
77  void SetDoClusterQA ( Int_t flag ) { fDoClusterQA = flag ;}
78  void SetUseTHnSparse ( Bool_t flag ) { fDoTHnSparse = flag ;}
79  void SetPlotHistsExtQA ( Bool_t flag ) { fSetPlotHistsExtQA = flag ;}
82  void SetAllowOverlapHeaders ( Bool_t allowOverlapHeader ) { fAllowOverlapHeaders = allowOverlapHeader ;}
83 
84  // Setting the cut lists for the conversion photons
85  void SetEventCutList ( Int_t nCuts,
86  TList *CutArray) {
87  fnCuts = nCuts ;
88  fEventCutArray = CutArray ;
89  }
90 
91  // Setting the cut lists for the conversion photons
93  TList *CutArray) {
94  fnCuts = nCuts ;
95  fCutArray = CutArray ;
96  }
97 
98  // Setting the cut lists for the calo photons
99  void SetCaloCutList ( Int_t nCuts,
100  TList *CutArray) {
101  fnCuts = nCuts ;
102  fClusterCutArray = CutArray ;
103  }
104 
105  // Setting the cut lists for the meson
106  void SetMesonCutList ( Int_t nCuts,
107  TList *CutArray) {
108  fnCuts = nCuts ;
109  fMesonCutArray = CutArray ;
110  }
111 
112  // BG HandlerSettings
113  void CalculateBackground ();
114  void CalculateBackgroundRP ();
115  void RotateParticle ( AliAODConversionPhoton *gamma );
117  Double_t previousEventEP,
118  Double_t thisEventEP );
121  Int_t pdgCode[] );
124  void UpdateEventByEventData ();
125 
126  // Additional functions for convenience
127  void SetLogBinningXTH2 ( TH2* histoRebin );
129  Int_t pdgCode );
130  Bool_t CheckVectorOnly ( vector<Int_t> &vec,
131  Int_t tobechecked );
132  Bool_t CheckVectorForDoubleCount ( vector<Int_t> &vec,
133  Int_t tobechecked );
134 
135  void FillMultipleCountMap ( map<Int_t,Int_t> &ma,
136  Int_t tobechecked );
138  Int_t> &ma,
139  TH1F* hist );
140 
141  // Function to enable MC label sorting
142  void SetEnableSortingOfMCClusLabels (Bool_t enableSort) { fEnableSortForClusMC = enableSort;}
143 
145 
146  protected:
147  AliV0ReaderV1* fV0Reader; // basic photon Selection Task
149  TString fCorrTaskSetting; // Correction Task Special Name
150  AliGammaConversionAODBGHandler** fBGHandler; // BG handler for Conversion
151  AliConversionAODBGHandlerRP** fBGHandlerRP; // BG handler for Conversion (possibility to mix with respect to RP)
152  AliGammaConversionAODBGHandler** fBGClusHandler; // BG handler for Cluster
153  AliConversionAODBGHandlerRP** fBGClusHandlerRP; // BG handler for Cluster (possibility to mix with respect to RP)
154  AliVEvent* fInputEvent; // current event
155  AliMCEvent* fMCEvent; // corresponding MC event
156  TList** fCutFolder; // Array of lists for containers belonging to cut
157  TList** fESDList; // Array of lists with histograms with reconstructed properties
158  TList** fBackList; // Array of lists with BG THnSparseF
159  TList** fMotherList; // Array of lists with Signal THnSparseF
160  TList** fGammaERM02; // Array of lists with conv photon shower shape trees
161  TList** fInvMassShowerShape; // Array of lists with invMass shower shape trees
162  TList** fTrueList; // Array of lists with histograms with MC validated reconstructed properties
163  TList** fMCList; // Array of lists with histograms with pure MC information
165  TList* fOutputContainer; // Output container
166  TClonesArray* fReaderGammas; // Array with conversion photons selected by V0Reader Cut
167  TList* fGammaCandidates; // current list of photon candidates
169  TList* fEventCutArray; // List with Event Cuts
170  AliConvEventCuts* fEventCuts; // EventCutObject
171  TList* fCutArray; // List with Conversion Cuts
172  AliConversionPhotonCuts* fConversionCuts; // ConversionCutObject
173  TList* fClusterCutArray; // List with Cluster Cuts
174  AliCaloPhotonCuts* fCaloPhotonCuts; // CaloPhotonCutObject
175  TList* fMesonCutArray; // List with Meson Cuts
176  AliConversionMesonCuts* fMesonCuts; // MesonCutObject
178  Bool_t fDoJetAnalysis; // Bool to produce Jet Plots
179  Bool_t fDoJetQA; // Bool to produce Jet QA Plots
180  TList** fJetHistograms; // Jet Histograms
181  TList** fTrueJetHistograms; // True Jet Histograms
182  Int_t fMaxPtNearEMCalPlace; // Place in jet vector of highest pt jet that is near the EMCal
183  Bool_t fJetNearEMCal; // If a jet is near the EMCal in the current event
184 
185  //histograms for Conversions reconstructed quantities
190  UChar_t fCharCatPhoton;
192  // 0: garbage,
193  // 1: background
194  // 2: secondary photon not from eta or k0s,
195  // 3: secondary photon from eta,
196  // 4: secondary photon from k0s,
197  // 5: dalitz
198  // 6: primary gamma
209 
224 
225  //histograms for mesons reconstructed quantities
242 
243  // histograms for rec photons tagged by Calo
246  // histograms for rec photon clusters
252 
253  //histograms for pure MC quantities
255  //histograms for pure MC quantities for gamma part
268 
269  //histograms for pure MC quantities for mesaon analysis
270  TH1F** fHistoMCPi0Pt;
273  TH1F** fHistoMCEtaPt;
294 
295  // histograms for pure MC quantites for dir gamma analysis part
304 
305  // MC validated reconstructed quantities mesons
358 
359  // MC validated reconstructed quantities photons
369 
392 
393  // MC validated reconstructed quantities for direct photon reco
404 
405  TH1F** fHistoPtJet; // Histogram of Jet Pt
406  TH1F** fHistoJetEta; // Histogram of Jet Eta
407  TH1F** fHistoJetPhi; // Histogram of Jet Phi
408  TH1F** fHistoJetArea; // Histogram of Jet Area
409  TH1F** fHistoNJets; // Histogram of number of jets
410  TH1F** fHistoEventwJets; // Histogram of number of events with jets > 0
411  TH1F** fHistoJetPi0PtRatio; // Histogram of PtPi0/PtJet
412  TH1F** fHistoDoubleCounting; // Histogram if NM candidates are defined within multiple jets
413  TH2F** fHistoJetMotherInvMassPt; // Histogram of NM candidates with a jet in the event
414  TH2F** fHistoPi0InJetMotherInvMassPt; // Histogram of NM candidates that are inside a jet
415  TH2F** fHistoMotherBackJetInvMassPt; // Histogram of Backgrouns candidates that are involved with jets
416  TH2F** fHistoRJetPi0Cand; // Histogram of RJetPi0Cand vs Pt
417  TH2F** fHistoEtaPhiJetPi0Cand; // Histogram of delta eta and delta phi distr between jet and NM candidates
418  TH2F** fHistoEtaPhiJetWithPi0Cand; // Histogram of delta eta and delta phi distr when pi0 is inside a jet
419  TH2F** fHistoJetFragmFunc; // Histogram to determine fragmentation function
420  TH2F** fHistoJetFragmFuncChargPart; // Histogram with fragmentation function as function of charged particles
421  TH2F** fHistoTruevsRecJetPt; // Histogram of true jet pt vs reconstructed jet pt
422  TH2F** fHistoTruePi0JetMotherInvMassPt; // Histogram of true pi0s in an event with a jet
423  TH2F** fHistoTruePi0InJetMotherInvMassPt; // Histogram of true pi0s in a jet
424  TH2F** fHistoTruePrimaryPi0JetInvMassPt; // Histogram of true primary pi0s in an event with a jet
425  TH2F** fHistoTruePrimaryPi0inJetInvMassPt; // Histogram of true primary pi0s in a jet
426  TH2F** fHistoTruePrimaryPi0InJetInvMassTruePt; // Histogram of true primary pi0s in a jet with their true pt
427  TH1F** fHistoTrueDoubleCountingPi0Jet; // Histogram of when a true pi0 is defined to be in multiple jets
428  TH2F** fHistoTrueEtaJetMotherInvMassPt; // Histogram of true etas in an event with a jet
429  TH2F** fHistoTrueEtaInJetMotherInvMassPt; // Histogram of true etas in a jet
430  TH2F** fHistoTruePrimaryEtaJetInvMassPt; // Histogram of true primary etas in an event with a jet
431  TH2F** fHistoTruePrimaryEtainJetInvMassPt; // Histogram of true primary etas in a jet
432  TH1F** fHistoTrueDoubleCountingEtaJet; // Histogram of when a true eta is defined to be in multiple jets
433  TH2F** fHistoTruePi0JetFragmFunc; // Histogram to determine true pi0 fragmentation function
434  TH2F** fHistoTrueEtaJetFragmFunc; // Histogram to determine true eta fragmentation function
435  TH2F** fHistoTruePi0JetFragmFuncChargPart; // Histogram with fragmentation function as function of charged particles
436  TH2F** fHistoTrueEtaJetFragmFuncChargPart; // Histogram with fragmentation function as function of charged particles
437  TH1F** fHistoMCPi0JetInAccPt; // Histogram with weighted pi0 in a jet event in acceptance, pT
438  TH1F** fHistoMCPi0inJetInAccPt; // Histogram with weighted pi0 in a jet in acceptance, pT
439  TH1F** fHistoMCEtaJetInAccPt; // Histogram with weighted eta in a jet event in acceptance, pT
440  TH1F** fHistoMCEtainJetInAccPt; // Histogram with weighted eta in a jet in acceptance, pT
441  TH1F** fHistoMCPi0JetEventGenerated; // Histogram with mesons in a jet event generated, pT
442  TH1F** fHistoMCPi0inJetGenerated; // Histogram with mesons in a jet generated, pT
443  TH1F** fHistoMCEtaJetEventGenerated; // Histogram with mesons in a jet event generated, pT
444  TH1F** fHistoMCEtainJetGenerated; // Histogram with mesons in a jet generated, pT
445  TH2F** fHistoTrueSecondaryPi0FromK0sJetInvMassPt; // Histogram with validated secondary mothers from K0s in an event with a jet, invMass, pt
446  TH2F** fHistoTrueSecondaryPi0FromK0sinJetInvMassPt; // Histogram with validated secondary mothers from K0s in a jet, invMass, pt
447  TH2F** fHistoTrueSecondaryPi0FromLambdaJetInvMassPt; // Histogram with validated secondary mothers from lambda in an event with a jet, invMass, pt
448  TH2F** fHistoTrueSecondaryPi0FromLambdainJetInvMassPt; // Histogram with validated secondary mothers from lambda in a jet, invMass, pt
449  TH2F** fHistoTrueSecondaryPi0FromK0lJetInvMassPt; // Histogram with validated secondary mothers from K0l in an event with a jet, invMass, pt
450  TH2F** fHistoTrueSecondaryPi0FromK0linJetInvMassPt; // Histogram with validated secondary mothers from K0l in a jet, invMass, pt
451  TH2F** fHistoTrueSecondaryPi0InvJetMassPt; // Histogram with validated secondary mothers in an event with a jet, invMass, pt
452  TH2F** fHistoTrueSecondaryPi0InvinJetMassPt; // Histogram with validated secondary mothers in a jet, invMass, pt
453  TH2F** fHistoMotherPi0inJetPtY; // Histogram with the rapidity of the validated pi0s in jets
454  TH2F** fHistoMotherEtainJetPtY; // Histogram with the rapidity of the validated etas in jets
455  TH2F** fHistoMotherPi0inJetPtPhi; // Histogram with the phi of the validated pi0s in jets
456  TH2F** fHistoMotherEtainJetPtPhi; // Histogram with the phi of the validated etas in jets
457  TH2F** fHistoUnfoldingAsData; // Histogram to use for jet pi0 unfolding
458  TH2F** fHistoUnfoldingMissed; // Histogram to use for jet pi0 unfolding
459  TH2F** fHistoUnfoldingReject; // Histogram to use for jet pi0 unfolding
460 
461  vector<Double_t> fVectorJetPt; // Vector of JetPt
462  vector<Double_t> fVectorJetPx; // Vector of JetPx
463  vector<Double_t> fVectorJetPy; // Vector of JetPy
464  vector<Double_t> fVectorJetPz; // Vector of JetPz
465  vector<Double_t> fVectorJetEta; // Vector of JetEta
466  vector<Double_t> fVectorJetPhi; // Vector of JetPhi
467  vector<Double_t> fVectorJetArea; // Vector of JetArea
468  vector<Double_t> fVectorJetNChargPart; // vector of Number of charged particles in Jet
469  vector<Double_t> fTrueVectorJetPt; // Vector of True JetPt
470  vector<Double_t> fTrueVectorJetPx; // Vector of True JetPx
471  vector<Double_t> fTrueVectorJetPy; // Vector of True JetPy
472  vector<Double_t> fTrueVectorJetPz; // Vector of True JetPz
473  vector<Double_t> fTrueVectorJetEta; // Vector of True JetEta
474  vector<Double_t> fTrueVectorJetPhi; // Vector of True JetPhi
475  vector<Double_t> fTrueVectorJetNChargPart; // Vector of number of charged particles in Jet
476 
477  // variable to keep track of multiple & missing reco
478  vector<Int_t> fVectorRecTruePi0s;
479  vector<Int_t> fVectorRecTrueEtas;
492  map<Int_t,Int_t> fMapMultipleCountTruePi0s;
493  map<Int_t,Int_t> fMapMultipleCountTrueEtas;
496 
497  // further QA and calibration hist for Calo clusters
502 
503  // event histograms
504  TH1F** fHistoNEvents;
507  TH1F** fHistoVertexZ;
508  TH1F** fHistoVertexX;
509  TH1F** fHistoVertexY;
514  TProfile** fProfileEtaShift;
517 
518  // hists for nonlineartiy calibration
519 // TH2F** fHistoTruePi0NonLinearity; //! E_truth/E_rec vs E_rec for TruePi0s
520 // TH2F** fHistoTrueEtaNonLinearity; //! E_truth/E_rec vs E_rec for TrueEtas
521 
522  // additional variables
523  Double_t fEventPlaneAngle; // EventPlaneAngle
524  TRandom3 fRandom; // random
525  Int_t fNGammaCandidates; // number of gamma candidates in event
526  Double_t* fUnsmearedPx; //[fNGammaCandidates]
527  Double_t* fUnsmearedPy; //[fNGammaCandidates]
528  Double_t* fUnsmearedPz; //[fNGammaCandidates]
529  Double_t* fUnsmearedE; //[fNGammaCandidates]
530  Int_t* fMCEventPos; //[fNGammaCandidates]
531  Int_t* fMCEventNeg; //[fNGammaCandidates]
532  Int_t* fESDArrayPos; //[fNGammaCandidates]
533  Int_t* fESDArrayNeg; //[fNGammaCandidates]
534  Int_t fnCuts; // number of cuts to be analysed in parallel
535  Int_t fiCut; // current cut
536  Bool_t fMoveParticleAccordingToVertex; // boolean for BG calculation
537  Int_t fIsHeavyIon; // switch for pp = 0, PbPb = 1, pPb = 2
538  Bool_t fDoLightOutput; // switch for running light output, kFALSE -> normal mode, kTRUE -> light mode
539  Bool_t fDoMesonAnalysis; // flag for meson analysis
540  Int_t fDoMesonQA; // flag for meson QA
541  Int_t fDoPhotonQA; // flag for photon QA
542  Int_t fDoClusterQA; // flag for cluster QA
543  Bool_t fIsFromDesiredHeader; // flag for MC headers
544  Bool_t fIsOverlappingWithOtherHeader; // flag for particles in MC overlapping between headers
545  Int_t fIsMC; // flag for MC information
546  Bool_t fDoTHnSparse; // flag for using THnSparses for background estimation
547  Bool_t fSetPlotHistsExtQA; // flag for extended QA hists
548  Double_t fWeightJetJetMC; // weight for Jet-Jet MC
549  Bool_t fDoConvGammaShowerShapeTree; // flag for tree with conv gamma R vs energy vs shower shape
550  Bool_t fEnableSortForClusMC; // switch on sorting for MC labels in cluster
551  Bool_t fDoPrimaryTrackMatching; // switch for basic track matching for primaries
552  Bool_t fDoInvMassShowerShapeTree; // flag for producing tree tESDInvMassShowerShape
553  TTree* tBrokenFiles; // tree for keeping track of broken files
554  TObjString* fFileNameBroken; // string object for broken file name
555  Bool_t fAllowOverlapHeaders; // enable overlapping headers for cluster selection
556  Int_t fTrackMatcherRunningMode; // CaloTrackMatcher running mode
557 
558  private:
559  AliAnalysisTaskGammaConvCalo(const AliAnalysisTaskGammaConvCalo&); // Prevent copy-construction
561 
562  ClassDef(AliAnalysisTaskGammaConvCalo, 49);
563 };
564 
565 #endif
Float_t tESDIMClusterM02
energy of clus gamma for tree
TH1F ** fHistoTruePrimaryConvGammaPt
array of histos with combinatorial BG, pt, source
TH2F ** fHistoMotherInvMassECalib
array of histograms with invariant mass cut of 0.45 && pi0cand->M() < 0.65 ,eta/phi of conversion pho...
TH2F ** fHistoDoubleCountTrueEtaInvMassPt
array of histos with double counted pi0s, invMass, pT
TH2F ** fHistoTruePi0CaloConvertedPhotonInvMassPt
array of histos with validated eta, photon leading, invMass, pt
TH1F ** fHistoClusGammaE
array of histos with cluster, pt
TH1F ** fHistoTrueK0sWithPi0DaughterMCPt
array of histos with validated secondary mothers from K0s, invMass, pt
TH2F ** fHistoTruePi0InvMassPt
array of histos with secondary pi0, pT gamma 2, source in acceptance
void ProcessTrueMesonCandidatesAOD(AliAODConversionMother *Pi0Candidate, AliAODConversionPhoton *TrueGammaCandidate0, AliAODConversionPhoton *TrueGammaCandidate1, Bool_t matched)
TH2F ** fHistoMCSecPi0PtGamma1vsSource
array of histos with unweighted pi0, pT gamma, different legs
void SetEventCutList(Int_t nCuts, TList *CutArray)
void ProcessTrueClusterCandidatesAOD(AliAODConversionPhoton *TruePhotonCandidate, Float_t clusM02)
double Double_t
Definition: External.C:58
TH2F ** fHistoMotherPi0PtAlpha
array of histograms with invariant mass cut of 0.45 && pi0cand->M() < 0.65, pt, Y ...
Bool_t CheckVectorOnly(vector< Int_t > &vec, Int_t tobechecked)
TH2F ** fHistoTruePrimaryPi0W0WeightingInvMassPt
array of histos with validated weighted primary mothers, invMass, pt
TProfile2D ** fProfileTruePrimaryEtaWeightsInvMassPt
array of profiles with weights for validated primary mothers, invMass, pt
void FillMultipleCountHistoAndClear(map< Int_t, Int_t > &ma, TH1F *hist)
TH2F ** fHistoTruePi0PureGammaInvMassECalib
array of histogram with pure pi0 signal inv Mass, energy of cluster
Definition: External.C:236
TH1F ** fHistoTrueClusDalitzMergedPt
array of histos with validated Dalitz decay, pt
TH1F ** fHistoNGoodESDTracks
array of histos with event information without event weights
TH2F ** fHistoTruePi0CaloConvPhotonConvRPt
array of histos with validated pi0 matched with conv photon, converted photon leading, invMass, pt
Float_t tESDGammaConvR
energy of conv gamma for tree
TH2F ** fHistoTrueEtaCaloConvPhotonConvRPt
array of histos with validated eta matched with conv photon, converted photon leading, invMass, pt
TH1F ** fHistoMCDecayGammaRhoPt
array of histos with decay gamma from pi0, pT
TH2F ** fHistoMCEtaPtJetPt
array of histos with weighted pi0, pT, hardest jet pt
TProfile ** fProfileEtaShift
array of histos with V0 counts
TH1F ** fHistoClusRejectedHeadersGammaPt
array of histos with cluster, pt all headers
AliGammaConversionAODBGHandler ** fBGClusHandler
Int_t tESDIMClusterClassification
leading cellID of cluster
map< Long_t, Int_t > tESDmapIsClusterAcceptedWithoutTrackMatch
sum of Et of tracks within R<0.2
Float_t tESDClusterMaxECell
NCells of matched cluster for tree.
TH1F ** fHistoMCDecayGammaOmegaPt
array of histos with decay gamma from eta, pT
TH1F ** fHistoMCSecEtaPt
array of histos with secondary pi0, pT, source in acceptance
TH1F ** fHistoMultipleCountTruePi0
vector containing labels of validated cluster photons
TList * fOutputContainer
Array of lists of output histograms for cluster photons.
TH2F ** fHistoTrueEtaCaloMergedClusterPartConvInvMassPt
array of histos with validated mothers, merged cluster part conv, invMass, pt
TH2F ** fHistoTrueSecondaryConvGammaPt
array of histos with validated primary conversion photon, rec pt, mc pt conversion photon...
TH2F ** fHistoNGoodESDTracksVsNGammaCandidates
array of histos with number of gamma candidates per event
TH2F ** fHistoTruePrimaryEtaW0WeightingInvMassPt
array of histos with validated unweighted primary mothers, invMass, pt
vector< Int_t > fVectorDoubleCountTrueEtas
vector containing labels of validated pi0
map< Int_t, Int_t > fMapMultipleCountTruePi0s
array of histos how often TrueClusterGammas are counted
TH2F ** fHistoTruePi0InvMassECalib
array of histos with TruePi0s: cluster E vs M02
TH2F ** fHistoTruePrimaryEtaInvMassPt
array of histos with validated weighted primary mothers, invMass, pt
Float_t tESDClusterNCells
Phi of matched cluster for tree.
TH2F ** fHistoTruePrimaryConvGammaESDPtMCPt
array of histos with validated primary
TH2F ** fHistoTrueBckGGInvMassPt
array of histos with lambda with reconstructed pi0 as daughter, pt
TH2F ** fHistoMCPi0WOWeightInAccPtGammaLeg
array of histos with weighted pi0, pT gamma, different legs
TH2F ** fHistoTrueMotherCaloEMNonLeadingInvMassPt
array of histos with validated mothers, merged cluster invMass, pt
TH1F ** fHistoMultipleCountTrueEta
array of histos how often TruePi0s are counted
TH2F ** fHistoTrueClusGammaEM02
map containing cluster photon labels that are counted at least twice
TH1F ** fHistoVertexZ
array of histos with number of good tracks (2010 Standard track cuts)
TH1F ** fHistoMCDecayGammaEtapPt
array of histos with decay gamma from omega, pT
Float_t tESDClusterEta
M20 of matched cluster for tree.
virtual void Terminate(const Option_t *)
TH1F ** fHistoTruePrimaryClusConvGammaPt
array of histos with validated primary cluster, photons, rec Pt, MC pt
TH2F ** fHistoTrueEtaPtAlpha
array of histos with validated pi0, pt, alpha
TH2F ** fHistoMotherPi0ConvPhotonEtaPhi
array of histograms with invariant mass cut of 0.45 && pi0cand->M() < 0.65, pt, openAngle ...
TH1F ** fHistoTrueClusGammaPt
array of histos with validated secondary conversion photon from Lambda, rec pt, mc pt conversion phot...
TH2F ** fHistoTrueMotherEtaConvPhotonEtaPhi
array of histograms with invariant mass cut of 0.05 && pi0cand->M() < 0.17 ,eta/phi of conversion pho...
Float_t tESDClusterNLM
MaxEnergyCell of matched cluster for tree.
TH2F ** fHistoPhotonPairMixedEventPtconv
array of histo for pairs vs. pt of converted photon
TH2F ** fHistoTruePi0CaloPhotonInvMassPt
array of histos with rejected eta, invMass, pt
TH2F ** fHistoTrueSecondaryConvGammaFromXFromLambdaMCPtESDPt
array of histos with validated secondary conversion photon from K0l, rec pt, mc pt conversion photon...
TH2F ** fHistoMotherEtaPtY
array of histograms with invariant mass cut of 0.05 && pi0cand->M() < 0.17, pt, Y ...
TTree ** fTreeConvGammaPtDcazCat
histogram conversion photon pT
TH2F ** fHistoTruePi0PtOpenAngle
array of histos with validated eta, pt, alpha
TH2F ** fHistoMCPi0InAccPtGammaLeg
array of histos with unweighted pi0, pT gamma, different legs
TH2F ** fHistoDoubleCountTruePi0InvMassPt
array of strings containing the stack position of the reconstructed validated eta ...
TH2F ** fHistoMCEtaPtAlpha
array of histos with weighted pi0, pT, alpha
void RotateParticleAccordingToEP(AliAODConversionPhoton *gamma, Double_t previousEventEP, Double_t thisEventEP)
UChar_t fCharCatPhoton
dcaz of conversion for tree
TH1F ** fHistoTrueClusDalitzPt
array of histos with validated merged partially converted photons, pt
TH1F ** fHistoMCPi0InAccPt
array of histos without evt weight eta, pT
TH2F ** fHistoTruePi0MatchedInvMassPt
array of histos with validated eta, invMass, pt
TH2F ** fHistoMCPi0PtJetPt
array of histos with secondary eta, source
TH2F ** fHistoMotherPi0PtOpenAngle
array of histograms with invariant mass cut of 0.45 && pi0cand->M() < 0.65, pt, alpha ...
TH1F ** fHistoNEvents
array of histogram with pure pi0 signal (only pure gammas) inv Mass, energy of cluster ...
TH2F ** fHistoTruePrimaryPi0W0WeightsPhotonPairPtconv
array of histos with validated primary pi0&#39;s vs conversion photon pT
TH2F ** fHistoMCPi0WOWeightPtGammaLeg
array of histos with weighted pi0, pT gamma, different legs
TProfile ** fProfileJetJetXSection
array of profiles with eta shift
TH2F ** fHistoMCSecPi0PtvsSource
array of histos with weighted primary particles, pT vs source
TH1F ** fHistoMCConvGammaPt
array of histos with decay gamma from Sigma0, pT
Short_t fDCAzPhoton
pt of conversion for tree
void SetConversionCutList(Int_t nCuts, TList *CutArray)
TH2F ** fHistoTrueMotherPi0ConvPhotonEtaPhi
array of histos with validated weighted primary eta, MCpt, resol pt
TH2F ** fHistoTrueSecondaryPi0FromLambdaPhotonPairPtconv
array of histos with validated secondary pi0&#39;s from K0l vs conversion photon pT
TH1F ** fHistoVertexX
array of histos with vertex z distribution for selected events
AliConversionAODBGHandlerRP ** fBGClusHandlerRP
TH2F ** fHistoTrueSecondaryPi0FromK0lPhotonPairPtconv
array of histos with validated secondary pi0&#39;s from K0s vs conversion photon pT
TH2F ** fHistoMCSecPi0InAccPtGamma2vsSource
array of histos with secondary pi0, pT gamma 2, source
Double_t fEventPlaneAngle
array of histos with ntrials for jetjet
TH2F ** fHistoSPDClusterTrackletBackground
array of histos with number of good tracks vs gamma candidates
THnSparseF ** fSparseMotherInvMassPtZM
array of histogram with signal + BG for same event photon pairs, inv Mass, pt
TH1F ** fHistoMCPi0Pt
array of histos wiht secondary converted gamma, pt
int Int_t
Definition: External.C:63
TH1F ** fHistoTruePrimaryPi0MissingPtconv
array of histos with validated primary pi0&#39;s vs conversion photon pT, double counting ...
TH1F ** fHistoTrueSecondaryClusGammaFromK0sPt
array of histos with validated secondary cluster, photons, pt
TH2F ** fHistoPhotonPairPtconv
array of histogram with BG for mixed event photon pairs, inv Mass, energy of cluster ...
Class handling all kinds of selection cuts for Gamma Calo analysis.
TH1F ** fHistoTrueClusShowerPt
array of histos with validated photon from electron, pt
Definition: External.C:204
TH1F ** fHistoNGammaCandidates
array of histos with vertex y distribution for selected events
float Float_t
Definition: External.C:68
TH1F ** fHistoClusOverlapHeadersGammaPt
array of histos with cluster, E
THnSparseF ** fSparseMotherBackInvMassPtZM
array of histogram with BG for mixed event photon pairs, inv Mass, pt
AliGammaConversionAODBGHandler ** fBGHandler
void SetEnableSortingOfMCClusLabels(Bool_t enableSort)
Float_t tESDIMClusterIsoSumClusterEt
matching residual closest track <-> cluster
TH2F ** fHistoMotherEtaPtAlpha
array of histograms with invariant mass cut of 0.05 && pi0cand->M() < 0.17, pt, alpha ...
TH1F ** fHistoClusGammaPt
array of histo for pairs vs. pt of converted photon
TH2F ** fHistoTrueSecondaryPi0InvMassPt
array of histograms with invariant mass cut of 0.45 && pi0cand->M() < 0.65 ,eta/phi of conversion pho...
TH2F ** fHistoCombinatorialPt
array of histos with validated conversion photon, pt
Float_t tESDIMClusterIsoSumTrackEt
sum of Et of clusters within R<0.2
Float_t tESDClusterPhi
Eta of matched cluster for tree.
TH2F ** fHistoTrueEtaPtY
array of histos with validated pi0, pt, Y
TH2F ** fHistoTruePrimaryClusConvGammaESDPtMCPt
array of histos with validated primary cluster, converted photons, pt
TH2F ** fHistoMCSecondaryConvGammaPt
array of histos with converted gamma, pT
TH2F ** fHistoTrueSecondaryConvGammaMCPt
array of histos with validated secondary conversion photon, rec pt for different sources ...
TH2F ** fHistoTruePi0CaloMergedClusterPartConvInvMassPt
array of histos with validated mothers, EM non leading, invMass, pt
TH2F ** fHistoMCPi0PtGammaLeg
array of histos with weighted eta, pT, hardest jet pt
UShort_t fPtGamma
tree with dca for conversions
UChar_t fCharPhotonMCInfo
category of conversion for tree
TH1F ** fHistoMCPi0WOWeightPt
array of histos with weighted pi0, pT
TH2F ** fHistoTruePrimaryPi0PhotonPairPtconv
array of histos with validated secondary cluster from Lambda, photons, pt
TH1F ** fHistoTrueClusPhotonFromElecMotherPt
array of histos with validated Dalitz decay, more than one decay product in cluster, pt
TH2F ** fHistoTrueBckFullMesonContainedInOneClusterInvMassPt
array of histos with pure gamma gamma combinatorial BG, invMass, pt
TH1F ** fHistoMCDecayGammaPi0Pt
array of histos with all secondary gamma, pT
TH2F ** fHistoTrueSecondaryConvGammaFromXFromK0sMCPtESDPt
array of histos with validated secondary conversion photon, MC pt for different sources ...
TH2F ** fHistoMotherInvMassPtAlpha
array of THnSparseF with BG for same event photon pairs, inv Mass, pt
TH2F ** fHistoMCSecPi0InAccPtGamma1vsSource
array of histos with secondary pi0, pT gamma 1, source
TH1F ** fHistoNV0Tracks
array of histos with SPD tracklets vs SPD clusters for background rejection
Int_t tESDIMClusterLeadCellID
M20 of cluster for tree.
TH2F ** fHistoMCPi0PtAlpha
array of histos with weighted eta, pT, Y
Class handling all kinds of selection cuts for Gamma Conversion analysis.
TH2F ** fHistoTruePrimaryEtaMCPtResolPt
array of histos with validated weighted primary pi0, MCpt, resol pt
TH1F ** fHistoMCPi0WOEvtWeightPt
array of histos with unweighted pi0, pT
TH1F ** fHistoMultipleCountTrueConvGamma
array of histos how often TrueEtas are counted
TH1F ** fHistoTrueEtaWithPi0DaughterMCPt
array of histos with validated secondary mothers from eta, invMass, pt
AliConversionPhotonCuts * fConversionCuts
TH2F ** fHistoDoubleCountTrueClusterGammaPt
array of histos with double counted photons, R, pT
TH1F ** fHistoMCEtaInAccPt
array of histos without evt weight pi0 in acceptance, pT
Int_t mode
Definition: anaM.C:41
TH2F ** fHistoTrueSecondaryPi0FromLambdaInvMassPt
array of histos with eta with reconstructed pi0 as daughter, pt
TH2F ** fHistoTruePi0CaloConvertedPhotonMatchedInvMassPt
array of histos with validated pi0, converted photon leading, invMass, pt
TH2F ** fHistoTrueSecondaryPi0FromK0sPhotonPairPtconv
array of histos with validated secondary pi0&#39;s vs conversion photon pT
TH1F ** fHistoMCSecPi0Source
array of histos with secondary pi0, pT, source
AliConversionAODBGHandlerRP ** fBGHandlerRP
TH2F ** fHistoMCPi0PtY
array of histos without evt weight eta in acceptance, pT
TH2F ** fHistoTrueEtaPtOpenAngle
array of histos with validated pi0, pt, openAngle
TH2F ** fHistoMCSecPi0PtGamma2vsSource
array of histos with secondary pi0, pT gamma 1, source in acceptance
TH2F ** fHistoTrueEtaCaloConvertedPhotonMatchedInvMassPt
array of histos with validated eta, converted photon leading, invMass, pt
Float_t tESDIMMesonInvMass
tree with cluster shower shape for meson cand
TH2F ** fHistoMCAllSecondaryGammaPt
array of histos with all gamma in EMCAL acceptance, pT
void SetMesonCutList(Int_t nCuts, TList *CutArray)
short Short_t
Definition: External.C:23
TH2F ** fHistoMotherInvMassPt
needed for filling of tree -> compute cluster isolation
Float_t tESDIMClusTrackDeltaEta
pT of closest matched track
TH2F ** fHistoMotherEtaConvPhotonEtaPhi
array of histograms with invariant mass cut of 0.05 && pi0cand->M() < 0.17 ,eta/phi of conversion pho...
TH2F ** fHistoMotherPi0PtY
array of histograms with alpha cut of 0.1 for inv mass vs pt
TH1F ** fHistoPtJet
array of histos with validated secondary pi0&#39;s vs conversion photon pT, missing vs source ...
Float_t tESDIMClusterM20
M02 of cluster for tree.
TH2F ** fHistoTrueSecondaryPi0DCPtconvSource
array of histos with validated secondary pi0&#39;s from Lambda vs conversion photon pT ...
TH1F ** fHistoTruePrimaryClusGammaPt
array of histos with number of labels in cluster
TH2F ** fHistoTrueSecondaryConvGammaFromXFromK0lMCPtESDPt
array of histos with validated secondary conversion photon from K0s, rec pt, mc pt conversion photon...
TH2F ** fHistoTruePi0PtY
array of histos with contamination BG, invMass, pt
void FillPhotonCombinatorialBackgroundHist(AliAODConversionPhoton *TruePhotonCandidate, Int_t pdgCode[])
Bool_t isMC
TH1F ** fHistoTrueClusMergedPartConvGammaPt
array of histos with validated merged photons, electrons, dalitz, pt
TTree ** tESDGammaERM02
MC info of conversion for tree.
TH1F ** fHistoMCPi0WOEvtWeightInAccPt
array of histos without weight pi0 in acceptance, pT
TH1F ** fHistoTrueConvGammaPt
array of histos with validated eta, pt, openAngle
TH1F ** fHistoMCDecayGammaPhiPt
array of histos with decay gamma from eta&#39;, pT
TH1F ** fHistoTrueClusElectronPt
array of histos with validated cluster (electron or photon), pt
Float_t tESDIMClusTrackDeltaPhi
matching residual closest track <-> cluster
void RotateParticle(AliAODConversionPhoton *gamma)
TH1F ** fHistoClusAllHeadersGammaPt
array of histos with cluster, pt overlapping with other headers
TH1F ** fHistoTrueClusMergedGammaPt
array of histos with validated converted photon, fully contained, pt
TH2F ** fHistoMCPrimaryPtvsSource
array of histos with weighted eta, pT, alpha
void ProcessTrueClusterCandidates(AliAODConversionPhoton *TruePhotonCandidate, Float_t clusM02)
TList * fEventCutArray
current list of cluster candidates
TH1F ** fHistoTrueClusNMothers
array of histos with pi0/eta/eta_prime in subleading contribution
TH2F ** fHistoTrueBckAsymEClustersInvMassPt
array of histos with pi0 fully contained in one cluster, invMass, pt
TH1F ** fHistoMCEtaWOEvtWeightPt
array of histos with unweighted eta, pT
Definition: External.C:220
TH1F ** fHistoVertexY
array of histos with vertex x distribution for selected events
TH2F ** fHistoTrueNLabelsInClusPt
array of histos with cluster with largest energy by hadron
Class handling all kinds of selection cuts for Gamma Conversion analysis.
TH2F ** fHistoMCEtaPtY
array of histos with weighted pi0, pT, Y
TH2F ** fHistoTruePi0PtAlpha
array of histos with validated eta, pt, Y
TH2F ** fHistoMotherBackInvMassECalib
array of histogram with signal + BG for same event photon pairs, inv Mass, energy of cluster ...
Float_t tESDClusE
tree with E_gamma vs R_conv vs M02 vs M20
TH2F ** fHistoTruePi0CaloMergedClusterInvMassPt
array of histos with validated mothers, electron leading, invMass, pt
TH1F ** fHistoTrueClusEMNonLeadingPt
array of histos with number of different particles (pi0/eta/eta_prime) contributing to cluster ...
Float_t tESDIMMesonPt
invMass of meson cand
TH1F ** fHistoTruePrimaryPi0DCPtconv
array of histos with validated primary pi0&#39;s vs conversion photon pT
vector< Int_t > fVectorDoubleCountTrueConvGammas
vector containing labels of validated eta
TProfile2D ** fProfileTruePrimaryPi0WeightsInvMassPt
array of histos with validated unweighted primary mothers, invMass, pt
Int_t GetSourceClassification(Int_t daughter, Int_t pdgCode)
Float_t tESDClusterM02
conv R of conv gamma for tree
TH1F ** fHistoTrueClusSubLeadingPt
array of histos with validated shower, pt
map< Int_t, Int_t > fMapMultipleCountTrueConvGammas
map containing eta labels that are counted at least twice
TH2F ** fHistoTrueSecondaryPi0FromK0lInvMassPt
array of histos with K0s with reconstructed pi0 as daughter, pt
TH1F ** fHistoMultipleCountTrueClusterGamma
array of histos how often TrueConvGammas are counted
TH2F ** fHistoMCSecPi0InAccPtvsSource
array of histos with secondary pi0, source
TH1F ** fHistoMCEtaPt
array of histos without event weights pi0, pT
vector< Int_t > fVectorDoubleCountTruePi0s
array of histos with double counted cluster photons
TH1F ** fHistoMCEtaWOWeightPt
array of histos with weighted eta, pT
Float_t tESDIMClusMatchedTrackPt
classification of cluster in MC
void ProcessTruePhotonCandidates(AliAODConversionPhoton *TruePhotonCandidate)
TH2F ** fHistoTrueEtaMatchedInvMassPt
array of histos with rejected pi0, invMass, pt
TH2F ** fHistoTrueSecondaryPi0MissingPtconvSource
array of histos with validated secondary pi0&#39;s vs conversion photon pT, double counting vs source ...
Class handling all kinds of selection cuts for Gamma Conversion analysis.
TH2F ** fHistoTruePrimaryPi0InvMassPt
array of histos with validated mothers, merged cluster part conv, invMass, pt
TH1F ** fHistoMCEtaWOWeightInAccPt
array of histos with weighted eta in acceptance, pT
TH1F ** fHistoMCPi0WOWeightInAccPt
array of histos with weighted pi0 in acceptance, pT
void ProcessTrueMesonCandidates(AliAODConversionMother *Pi0Candidate, AliAODConversionPhoton *TrueGammaCandidate0, AliAODConversionPhoton *TrueGammaCandidate1, Bool_t matched)
unsigned short UShort_t
Definition: External.C:28
TH1F ** fHistoMCAllGammaPt
array of histos for header names
TH2F ** fHistoTrueEtaCaloMergedClusterInvMassPt
array of histos with validated mothers, merged cluster invMass, pt
const char Option_t
Definition: External.C:48
void ProcessTruePhotonCandidatesAOD(AliAODConversionPhoton *TruePhotonCandidate)
map< Int_t, Int_t > fMapMultipleCountTrueClusterGammas
map containing photon labels that are counted at least twice
TH1F ** fHistoTrueK0lWithPi0DaughterMCPt
array of histos with validated secondary mothers from K0l, invMass, pt
TH1F ** fHistoTrueSecondaryClusGammaFromLambdaPt
array of histos with validated secondary cluster from K0l, photons, pt
TH2F ** fHistoTruePrimaryClusGammaESDPtMCPt
array of histos with validated primary cluster, photons, pt
void FillMultipleCountMap(map< Int_t, Int_t > &ma, Int_t tobechecked)
TH2F ** fHistoMotherMatchedInvMassPt
array of histogram with signal + BG for same event photon pairs, inv Mass, pt
void SetCorrectionTaskSetting(TString setting)
TH2F ** fHistoTrueSecondaryPi0FromEtaInvMassPt
array of histos with K0l with reconstructed pi0 as daughter, pt
TH2F ** fHistoTrueSecondaryPi0FromK0sInvMassPt
array of histos with validated secondary mothers, invMass, pt
TH1F ** fHistoMCAllGammaEMCALAccPt
array of histos with all gamma, pT
void MoveParticleAccordingToVertex(AliAODConversionPhoton *particle, const AliGammaConversionAODBGHandler::GammaConversionVertex *vertex)
TH2F ** fHistoTrueBckContInvMassPt
array of histos with asymmetry energy distributions of clusters, invMass, pt
bool Bool_t
Definition: External.C:53
Bool_t CheckVectorForDoubleCount(vector< Int_t > &vec, Int_t tobechecked)
TH1I ** fHistoMCHeaders
array of histos with cluster, pt rejected headers
Implementation of a sample jet analysis task.
TH1F ** fHistoTrueSecondaryClusGammaPt
array of histos with validated primary cluster, converted photons, rec Pt, MC pt
TH1F ** fHistoTrueLambdaWithPi0DaughterMCPt
array of histos with validated secondary mothers from Lambda, invMass, pt
TH1F ** fHistoTrueSecondaryClusGammaFromK0lPt
array of histos with validated secondary cluster from K0s, photons, pt
TH2F ** fHistoTrueSecondaryPi0PhotonPairPtconv
array of histos with validated primary pi0&#39;s vs conversion photon pT, missing
void SetCaloCutList(Int_t nCuts, TList *CutArray)
vector< Int_t > fVectorDoubleCountTrueClusterGammas
vector containing labels of validated photons
TH2F ** fHistoTrueEtaCaloPhotonInvMassPt
array of histos with validated pi0, photon leading, invMass, pt
TH1F ** fHistoMCEtaWOEvtWeightInAccPt
array of histos without weight eta in acceptance, pT
TTree ** tESDInvMassShowerShape
NLM of matched cluster for tree.
Float_t tESDClusterM20
M02 of matched cluster for tree.
TH1F ** fHistoMCDecayGammaSigmaPt
array of histos with decay gamma from phi, pT
TH1F ** fHistoNEventsWOWeight
array of histos with event information
void SetAllowOverlapHeaders(Bool_t allowOverlapHeader)
TH2F ** fHistoMotherEtaPtOpenAngle
array of histograms with invariant mass cut of 0.05 && pi0cand->M() < 0.17, pt, openAngle ...
TH2F ** fHistoTruePrimaryPi0MCPtResolPt
array of profiles with weights for validated primary mothers, invMass, pt
TH2F ** fHistoDoubleCountTrueConvGammaRPt
array of histos with double counted etas, invMass, pT
TH2F ** fHistoTrueEtaCaloElectronInvMassPt
array of histos with validated mothers, electron leading, invMass, pt
TH1F ** fHistoTrueClusConvGammaPt
array of histos with validated electron, pt
TH1F ** fHistoMCDecayGammaEtaPt
array of histos with decay gamma from rho, pT
TH2F ** fHistoMotherBackInvMassPt
array of THnSparseF with signal + BG for same event photon pairs, inv Mass, pt
map< Int_t, Int_t > fMapMultipleCountTrueEtas
map containing pi0 labels that are counted at least twice
TH1F ** fHistoMCSecEtaSource
array of histos with secondary eta, pT
vector< Int_t > fVectorRecTrueEtas
array of strings containing the stack position of the reconstructed validated pi0 ...
TH1F ** fHistoJetJetNTrials
array of profiles with xsection for jetjet
TH2F ** fHistoTrueEtaInvMassPt
array of histos with validated pi0, invMass, pt
TH1F ** fHistoTrueClusConvGammaFullyPt
array of histos with validated converted photon, pt
TH2F ** fHistoTrueClusPi0EM02
array of histos with TruePhotons: cluster E vs M02
AliAnalysisTaskGammaConvCalo & operator=(const AliAnalysisTaskGammaConvCalo &)