Float_t GetPtThresholds(Int_t i) const { return fPtThresholds[i] ; }
Float_t GetSumPtThresholds(Int_t i) const { return fSumPtThresholds[i]; }
Float_t GetPtFractions(Int_t i) const { return fPtFractions[i] ; }
- Int_t GetNumberOfSMCoveredByTRD() const { return fTRDSMCovered ; }
Int_t GetMCIndex(Int_t mcTag);
void SetPtThresholds(Int_t i, Float_t pt) { fPtThresholds[i] = pt ; }
void SetPtFractions(Int_t i, Float_t pt) { fPtFractions[i] = pt ; }
void SetSumPtThresholds(Int_t i, Float_t pt){ fSumPtThresholds[i] = pt ; }
- void SetNumberOfSMCoveredByTRD(Int_t n) { fTRDSMCovered = n ; }
Bool_t IsReIsolationOn() const { return fReMakeIC ; }
void SwitchOnReIsolation() { fReMakeIC = kTRUE ; }
void SwitchOnDecayTaggedHistoFill() { fFillTaggedDecayHistograms = kTRUE ; }
void SwitchOffDecayTaggedHistoFill() { fFillTaggedDecayHistograms = kFALSE; }
- //Histogrammes setters and getters
-
- void SetHistoPtSumRangeAndNBins(Float_t min, Float_t max, Int_t n) {
- fHistoNPtSumBins = n ; fHistoPtSumMax = max ; fHistoPtSumMin = min ; }
-
- Int_t GetHistoNPtSumBins() const { return fHistoNPtSumBins ; }
- Float_t GetHistoPtSumMin() const { return fHistoPtSumMin ; }
- Float_t GetHistoPtSumMax() const { return fHistoPtSumMax ; }
-
- void SetHistoPtInConeRangeAndNBins(Float_t min, Float_t max, Int_t n) {
- fHistoNPtInConeBins = n ; fHistoPtInConeMax = max ; fHistoPtInConeMin = min ; }
-
- Int_t GetHistoNPtInConeBins() const { return fHistoNPtInConeBins; }
- Float_t GetHistoPtInConeMin() const { return fHistoPtInConeMin ; }
- Float_t GetHistoPtInConeMax() const { return fHistoPtInConeMax ; }
-
// For primary histograms in arrays, index in the array, corresponding to a photon origin
- enum mcPrimTypes { kmcPrimPhoton = 0, kmcPrimPi0Decay = 1, kmcPrimOtherDecay = 2,
- kmcPrimOther = 3, kmcPrimPrompt = 4,
- kmcPrimFrag = 5, kmcPrimISR = 6 } ;
+ enum mcPrimTypes { kmcPrimPhoton = 0, kmcPrimPi0Decay = 1, kmcPrimOtherDecay = 2,
+ kmcPrimPrompt = 3, kmcPrimFrag = 4, kmcPrimISR = 5 } ;
// For histograms in arrays, index in the array, corresponding to any particle origin
enum mcTypes { kmcPhoton = 0, kmcPrompt = 1, kmcFragment = 2,
Bool_t fFillHighMultHistograms; // Fill high multiplicity histograms
Bool_t fFillTaggedDecayHistograms; // Fill histograms for clusters tagged as decay
Bool_t fFillNLMHistograms; // Fill NLM histograms
- Int_t fTRDSMCovered; // From which SM EMCal is covered by TRD
// Analysis data members for multiple cones and pt thresholds
Int_t fNCones ; //! Number of cone sizes to test
//MC
//
- TH2F * fhEtaPrimMC [7]; //! Pt vs Eta of generated photon
- TH2F * fhPhiPrimMC [7]; //! Pt vs Phi of generated photon
- TH1F * fhEPrimMC [7]; //! Number of generated photon vs E
- TH1F * fhPtPrimMCiso[7]; //! Number of generated isolated photon vs pT
+ TH2F * fhEtaPrimMC [6]; //! Pt vs Eta of generated photon
+ TH2F * fhPhiPrimMC [6]; //! Pt vs Phi of generated photon
+ TH1F * fhEPrimMC [6]; //! Number of generated photon vs E
+ TH1F * fhPtPrimMCiso[6]; //! Number of generated isolated photon vs pT
TH1F * fhPtNoIsoMC[9]; //! Number of not isolated mcTypes particle
TH1F * fhPtIsoMC [9]; //! Number of isolated mcTypes particle
TH2F * fhPhiIsoMC [9]; //! Phi of isolated mcTypes particle
TH2F * fhEtaIsoMC [9]; //! eta of isolated mcTypes particle
- TH1F * fhPtThresIsolatedMC[9][5][5]; //! Isolated mcTypes particle with pt threshold
- TH1F * fhPtFracIsolatedMC [9][5][5]; //! Isolated mcTypes particle with pt frac
- TH2F * fhPtSumIsolatedMC [9][5]; //! Isolated mcTypes particle with threshold on cone pt sum
+ TH1F * fhPtDecayIsoMC[9] ; //! Number of isolated Pi0 decay particles (invariant mass tag) for a mcTypes particle
+ TH1F * fhPtDecayNoIsoMC[9] ; //! Number of not isolated Pi0 decay particles (invariant mass tag) for a mcTypes particle
+
TH2F * fhPtLambda0MC[9][2]; //! Shower shape of (non) isolated candidates originated by mcTypes particle (do not apply SS cut previously)
- //
-
+
+ // Multiple cut analysis
TH2F * fhSumPtLeadingPt[5] ; //! Sum Pt in the cone
TH2F * fhPtLeadingPt[5] ; //! Particle Pt in the cone
TH2F * fhPerpSumPtLeadingPt[5] ; //! Sum Pt in the cone at the perpendicular phi region to trigger axis (phi +90)
TH1F * fhPtThresIsolated[5][5] ; //! Isolated particle with pt threshold
TH1F * fhPtFracIsolated[5][5] ; //! Isolated particle with pt threshold frac
- TH1F * fhPtSumIsolated[5][5] ; //! Isolated particle with threshold on cone pt sum
+ TH1F * fhSumPtIsolated[5][5] ; //! Isolated particle with threshold on cone pt sum
TH2F * fhEtaPhiPtThresIso[5][5] ; //! eta vs phi of isolated particles with pt threshold
TH2F * fhEtaPhiPtThresDecayIso[5][5] ; //! eta vs phi of isolated particles with pt threshold
TH2F * fhEtaPhiFracPtSumIso[5][5]; //! Isolated particle with threshold on cone sum density
TH2F * fhEtaPhiFracPtSumDecayIso[5][5]; //! Isolated particle with threshold on cone sum density
+ // Multiple cut MC
+ TH1F * fhPtThresIsolatedMC[9][5][5]; //! Isolated mcTypes particle with pt threshold
+ TH1F * fhPtFracIsolatedMC [9][5][5]; //! Isolated mcTypes particle with pt frac
+ TH1F * fhSumPtIsolatedMC [9][5][5]; //! Isolated mcTypes particle with threshold on cone pt sum
+ TH2F * fhSumPtLeadingPtMC [9][5]; //! mcTypes particle for sum Pt, different cone
+
+
// Track matching studies
TH2F * fhTrackMatchedDEta[2] ; //! Eta distance between track and cluster vs cluster E
TH2F * fhTrackMatchedDPhi[2] ; //! Phi distance between track and cluster vs cluster E
TH2F * fhTimePileUpMainVertexZDistance; //! time of cluster vs difference of z main vertex and pile-up vertex
TH2F * fhTimePileUpMainVertexZDiamond; //! time of cluster vs difference of z diamond and pile-up vertex
- //Histograms settings
- Int_t fHistoNPtSumBins; // Number of bins in PtSum histograms
- Float_t fHistoPtSumMax; // PtSum maximum in histogram
- Float_t fHistoPtSumMin; // PtSum minimum in histogram
- Int_t fHistoNPtInConeBins; // Number of bins in PtInCone histogram
- Float_t fHistoPtInConeMax; // PtInCone maximum in histogram
- Float_t fHistoPtInConeMin; // PtInCone maximum in histogram
-
AliAnaParticleIsolation( const AliAnaParticleIsolation & iso) ; // cpy ctor
AliAnaParticleIsolation & operator = (const AliAnaParticleIsolation & iso) ; // cpy assignment
- ClassDef(AliAnaParticleIsolation,26)
+ ClassDef(AliAnaParticleIsolation,27)
} ;