// Analysis methods
- Bool_t ClusterSelected(AliVCluster* cl, TLorentzVector mom) ;
+ Bool_t ClusterSelected(AliVCluster* cl, TLorentzVector mom, Int_t nlm) ;
void FillAcceptanceHistograms();
- void FillShowerShapeHistograms( AliVCluster* cluster, const Int_t mcTag) ;
+ void FillEMCALTriggerClusterBCHistograms(Int_t idcalo, Float_t ecluster, Float_t tofcluster,
+ Float_t etacluster, Float_t phicluster);
- void SwitchOnFillShowerShapeHistograms() { fFillSSHistograms = kTRUE ; }
- void SwitchOffFillShowerShapeHistograms() { fFillSSHistograms = kFALSE ; }
+ void FillShowerShapeHistograms( AliVCluster* cluster, Int_t mcTag) ;
- void FillTrackMatchingResidualHistograms(AliVCluster* calo, const Int_t cut);
+ void SwitchOnFillShowerShapeHistograms() { fFillSSHistograms = kTRUE ; }
+ void SwitchOffFillShowerShapeHistograms() { fFillSSHistograms = kFALSE ; }
- void SwitchOnTMHistoFill() { fFillTMHisto = kTRUE ; }
- void SwitchOffTMHistoFill() { fFillTMHisto = kFALSE ; }
+ void SwitchOnOnlySimpleSSHistoFill() { fFillOnlySimpleSSHisto = kTRUE ; }
+ void SwitchOffOnlySimpleHistoFill() { fFillOnlySimpleSSHisto = kFALSE ; }
+
+ void FillTrackMatchingResidualHistograms(AliVCluster* calo, Int_t cut);
+
+ void SwitchOnTMHistoFill() { fFillTMHisto = kTRUE ; }
+ void SwitchOffTMHistoFill() { fFillTMHisto = kFALSE ; }
+
+ void FillClusterPileUpHistograms(AliVCluster * calo, Bool_t matched,
+ Float_t ptcluster, Float_t etacluster,
+ Float_t phicluster, Float_t l0cluster);
+
+ void FillPileUpHistograms(Float_t energy, Float_t pt, Float_t time) ;
+ void FillPileUpHistogramsPerEvent() ;
+
+ void SwitchOnFillPileUpHistograms() { fFillPileUpHistograms = kTRUE ; }
+ void SwitchOffFillPileUpHistograms() { fFillPileUpHistograms = kFALSE ; }
+
+ void SwitchOnFillEMCALBCHistograms() { fFillEMCALBCHistograms = kTRUE ; }
+ void SwitchOffFillEMCALBCHistograms() { fFillEMCALBCHistograms = kFALSE ; }
// Analysis parameters setters getters
void SetNCellCut(Int_t n) { fNCellsCut = n ; }
Double_t GetNCellCut() const { return fNCellsCut ; }
+ void SetNLMCut(Int_t min, Int_t max) { fNLMCutMin = min;
+ fNLMCutMax = max ; }
+ Int_t GetNLMCutMin() const { return fNLMCutMin ; }
+ Int_t GetNLMCutMax() const { return fNLMCutMax ; }
+
+
Bool_t IsTrackMatchRejectionOn() const { return fRejectTrackMatch ; }
void SwitchOnTrackMatchRejection() { fRejectTrackMatch = kTRUE ; }
void SwitchOffTrackMatchRejection() { fRejectTrackMatch = kFALSE ; }
kmcISR = 12, kmcString = 13 };
enum mcPTypes { kmcPPhoton = 0, kmcPPi0Decay = 1, kmcPOtherDecay = 2, kmcPOther = 3,
- kmcPPrompt = 4, kmcPFragmentation = 5, kmcPISR = 6 };
+ kmcPPrompt = 4, kmcPFragmentation = 5, kmcPISR = 6 };
enum mcssTypes { kmcssPhoton = 0, kmcssOther = 1, kmcssPi0 = 2,
kmcssEta = 3, kmcssConversion = 4, kmcssElectron = 5 };
private:
- TString fCalorimeter ; // Calorimeter where the gamma is searched;
- Float_t fMinDist ; // Minimal distance to bad channel to accept cluster
- Float_t fMinDist2; // Cuts on Minimal distance to study acceptance evaluation
- Float_t fMinDist3; // One more cut on distance used for acceptance-efficiency study
- Bool_t fRejectTrackMatch ; // If PID on, reject clusters which have an associated TPC track
- Bool_t fFillTMHisto; // Fill track matching plots
- Double_t fTimeCutMin ; // Remove clusters/cells with time smaller than this value, in ns
- Double_t fTimeCutMax ; // Remove clusters/cells with time larger than this value, in ns
- Int_t fNCellsCut ; // Accept for the analysis clusters with more than fNCellsCut cells
- Bool_t fFillSSHistograms ; // Fill shower shape histograms
- Int_t fNOriginHistograms; // Fill only NOriginHistograms of the 14 defined types
- Int_t fNPrimaryHistograms; // Fill only NPrimaryHistograms of the 7 defined types
-
+ TString fCalorimeter ; // Calorimeter where the gamma is searched;
+ Float_t fMinDist ; // Minimal distance to bad channel to accept cluster
+ Float_t fMinDist2; // Cuts on Minimal distance to study acceptance evaluation
+ Float_t fMinDist3; // One more cut on distance used for acceptance-efficiency study
+ Bool_t fRejectTrackMatch ; // If PID on, reject clusters which have an associated TPC track
+ Bool_t fFillTMHisto; // Fill track matching plots
+ Double_t fTimeCutMin ; // Remove clusters/cells with time smaller than this value, in ns
+ Double_t fTimeCutMax ; // Remove clusters/cells with time larger than this value, in ns
+ Int_t fNCellsCut ; // Accept for the analysis clusters with more than fNCellsCut cells
+ Int_t fNLMCutMin ; // Remove clusters/cells with number of local maxima smaller than this value
+ Int_t fNLMCutMax ; // Remove clusters/cells with number of local maxima larger than this value
+ Bool_t fFillSSHistograms ; // Fill shower shape histograms
+ Bool_t fFillOnlySimpleSSHisto; // Fill selected cluster histograms, selected SS histograms
+ Int_t fNOriginHistograms; // Fill only NOriginHistograms of the 14 defined types
+ Int_t fNPrimaryHistograms; // Fill only NPrimaryHistograms of the 7 defined types
+ Bool_t fFillPileUpHistograms; // Fill pile-up related histograms
+ Bool_t fFillEMCALBCHistograms; // Fill eta-phi BC dependent histograms
+
//Histograms
- TH1F * fhClusterCuts[9]; //! control histogram on the different photon selection cuts
- TH2F * fhNCellsE; //! number of cells in cluster vs E
- TH2F * fhMaxCellDiffClusterE; //! Fraction of energy carried by cell with maximum energy
- TH2F * fhTimeE; //! time of cluster vs E
-
- TH1F * fhEPhoton ; //! Number of identified photon vs energy
- TH1F * fhPtPhoton ; //! Number of identified photon vs transerse momentum
- TH2F * fhPhiPhoton ; //! Azimuthal angle of identified photon vs transerse momentum
- TH2F * fhEtaPhoton ; //! Pseudorapidity of identified photon vs transerse momentum
- TH2F * fhEtaPhiPhoton ; //! Pseudorapidity vs Phi of identified photon for transerse momentum > 0.5
- TH2F * fhEtaPhi05Photon ; //! Pseudorapidity vs Phi of identified photon for transerse momentum < 0.5
+ TH1F * fhClusterCuts[10]; //! control histogram on the different photon selection cuts
+ TH2F * fhNCellsE; //! number of cells in cluster vs E
+ TH2F * fhCellsE; //! energy of cells in cluster vs E of cluster
+ TH2F * fhMaxCellDiffClusterE; //! Fraction of energy carried by cell with maximum energy
+ TH2F * fhTimePt; //! time of photon cluster vs pt
+
+ TH2F * fhEtaPhi ; //! Pseudorapidity vs Phi of clusters for E > 0.5
+ TH2F * fhEtaPhiEMCALBC0 ; //! Pseudorapidity vs Phi of clusters for E > 0.5
+ TH2F * fhEtaPhiEMCALBC1 ; //! Pseudorapidity vs Phi of clusters for E > 0.5
+ TH2F * fhEtaPhiEMCALBCN ; //! Pseudorapidity vs Phi of clusters for E > 0.5
+
+ TH2F * fhEtaPhiTriggerEMCALBC[11] ; //! Pseudorapidity vs Phi of clusters for E > 2
+ TH2F * fhTimeTriggerEMCALBC [11] ; //! Time distribution of clusters, when trigger is in a given BC
+ TH2F * fhTimeTriggerEMCALBCPileUpSPD[11]; //! Time distribution of clusters, when trigger is in a given BC, tagged as pile-up SPD
+
+ TH2F * fhEtaPhiTriggerEMCALBCUM[11] ; //! Pseudorapidity vs Phi of clusters for E > 2, not matched to trigger
+ TH2F * fhTimeTriggerEMCALBCUM [11] ; //! Time distribution of clusters, when trigger is in a given BC, not matched to trigger
+
+ TH2F * fhEtaPhiTriggerEMCALBCCluster [11] ; //! Pseudorapidity vs Phi of trigger clusters
+ TH2F * fhTimeTriggerEMCALBCCluster ; //! Time distribution of clusters, when trigger cluster is in a given BC
+ TH2F * fhEtaPhiTriggerEMCALBCUMCluster[11] ; //! Pseudorapidity vs Phi of highest E cluster in event, not matched to trigger
+ TH2F * fhTimeTriggerEMCALBCUMCluster ; //! Time distribution of highest energy cluster in event, when trigger is in a given BC, not
+
+ TH2F * fhEtaPhiTriggerEMCALBCClusterOverTh ; //! Pseudorapidity vs Phi of trigger clusters, over nominal threshold
+ TH2F * fhEtaPhiTriggerEMCALBCUMClusterOverTh ; //! Pseudorapidity vs Phi of highest E cluster in event, not matched to trigger, over nominal threshold
+ TH2F * fhEtaPhiTriggerEMCALBCClusterBelowTh1 ; //! Pseudorapidity vs Phi of trigger clusters, 1 GeV below nominal threshold
+ TH2F * fhEtaPhiTriggerEMCALBCUMClusterBelowTh1 ; //! Pseudorapidity vs Phi of highest E cluster in event, not matched to trigger, 2 GeV below nominal threshold
+ TH2F * fhEtaPhiTriggerEMCALBCClusterBelowTh2 ; //! Pseudorapidity vs Phi of trigger clusters, 1 GeV below nominal threshold
+ TH2F * fhEtaPhiTriggerEMCALBCUMClusterBelowTh2 ; //! Pseudorapidity vs Phi of highest E cluster in event, not matched to trigger, 2 GeV below nominal threshold
+
+ TH2F * fhEtaPhiTriggerEMCALBCExotic ; //! Pseudorapidity vs Phi of trigger exotic clusters
+ TH2F * fhTimeTriggerEMCALBCExotic ; //! Time distribution of clusters, when trigger exotic cluster
+ TH2F * fhEtaPhiTriggerEMCALBCUMExotic ; //! Pseudorapidity vs Phi of highest E exotic cluster in event, not matched to trigger
+ TH2F * fhTimeTriggerEMCALBCUMExotic ; //! Time distribution of highest energy exotic cluster in event, not matched to trigger
+
+ TH2F * fhEtaPhiTriggerEMCALBCBad ; //! Pseudorapidity vs Phi of trigger exotic clusters
+ TH2F * fhTimeTriggerEMCALBCBad ; //! Time distribution of clusters, when trigger exotic
+ TH2F * fhEtaPhiTriggerEMCALBCUMBad ; //! Pseudorapidity vs Phi of highest E exotic cluster in event, not matched to trigger
+ TH2F * fhTimeTriggerEMCALBCUMBad ; //! Time distribution of highest energy exotic cluster in event, not matched to trigger
+
+ TH2F * fhEtaPhiTriggerEMCALBCBadExotic ; //! Pseudorapidity vs Phi of trigger exotic and bad clusters
+ TH2F * fhTimeTriggerEMCALBCBadExotic ; //! Time distribution of clusters, when trigger exotic and bad cluster
+ TH2F * fhEtaPhiTriggerEMCALBCUMBadExotic ; //! Pseudorapidity vs Phi of highest E exotic cluster in event, not matched to trigger
+ TH2F * fhTimeTriggerEMCALBCUMBadExotic ; //! Time distribution of highest energy exotic cluster in event, not matched to trigger
+
+ TH2F * fhEtaPhiTriggerEMCALBCExoticCluster ; //! Pseudorapidity vs Phi of trigger exotic clusters
+ TH2F * fhTimeTriggerEMCALBCExoticCluster ; //! Time distribution of clusters, when trigger exotic cluster
+ TH2F * fhEtaPhiTriggerEMCALBCUMExoticCluster ; //! Pseudorapidity vs Phi of highest E exotic cluster in event, not matched to trigger
+ TH2F * fhTimeTriggerEMCALBCUMExoticCluster ; //! Time distribution of highest energy exotic cluster in event, not matched to trigger
+
+ TH2F * fhEtaPhiTriggerEMCALBCBadCluster ; //! Pseudorapidity vs Phi of trigger bad clusters
+ TH2F * fhTimeTriggerEMCALBCBadCluster ; //! Time distribution of clusters, when trigger bad cluster is in a given BC
+ TH2F * fhEtaPhiTriggerEMCALBCUMBadCluster ; //! Pseudorapidity vs Phi of highest E bad cluster in event, not matched to trigger
+ TH2F * fhTimeTriggerEMCALBCUMBadCluster ; //! Time distribution of highest energy bad cluster in event, when trigger is in a given BC, not
+
+ TH2F * fhEtaPhiTriggerEMCALBCBadExoticCluster ; //! Pseudorapidity vs Phi of trigger exotic and bad clusters
+ TH2F * fhTimeTriggerEMCALBCBadExoticCluster ; //! Time distribution of clusters, when trigger exotic and bad cluster
+ TH2F * fhEtaPhiTriggerEMCALBCUMBadExoticCluster; //! Pseudorapidity vs Phi of highest E exotic and bad cluster in event, not matched to trigger
+ TH2F * fhTimeTriggerEMCALBCUMBadExoticCluster ; //! Time distribution of highest energy exotic and bad cluster in event, not matched to trigger
+
+ TH2F * fhTimeTriggerEMCALBCBadMaxCell ; //! Time distribution of trigger clusters, when trigger bad max cell
+ TH2F * fhTimeTriggerEMCALBCUMBadMaxCell ; //! Time distribution of highest energy bad max cell cluster in event, when trigger is not found
+ TH2F * fhTimeTriggerEMCALBCBadMaxCellExotic ; //! Time distribution of trigger clusters, when trigger exotic cluster with bad max cell
+ TH2F * fhTimeTriggerEMCALBCUMBadMaxCellExotic ; //! Time distribution of highest energy exotic with bad max cell cluster in event, when trigger is not found
+
+ TH2F * fhEtaPhiTriggerEMCALBCUMReMatchOpenTimeCluster ; //! Pseudorapidity vs Phi of highest E bad cluster in event, not matched to trigger, rematched open time trigger
+ TH2F * fhTimeTriggerEMCALBCUMReMatchOpenTimeCluster ; //! Time distribution of highest energy bad max cell cluster in event, when trigger is not found, rematched open time trigger
+ TH2F * fhEtaPhiTriggerEMCALBCUMReMatchCheckNeighCluster; //! Pseudorapidity vs Phi of highest E bad cluster in event, not matched to trigger, rematched with neigbour patchs
+ TH2F * fhTimeTriggerEMCALBCUMReMatchCheckNeighCluster ; //! Time distribution of highest energy bad max cell cluster in event, when trigger is not found, rematched with neigbour patchs
+ TH2F * fhEtaPhiTriggerEMCALBCUMReMatchBothCluster;//! Pseudorapidity vs Phi of highest E bad cluster in event, not matched to trigger, rematched open both
+ TH2F * fhTimeTriggerEMCALBCUMReMatchBothCluster ; //! Time distribution of highest energy bad max cell cluster in event, when trigger is not found, rematched open both
+
+ TH2F * fhTimeTriggerEMCALBC0UMReMatchOpenTime ; //! Time distribution of clusters, not matched to trigger, rematched open time trigger
+ TH2F * fhTimeTriggerEMCALBC0UMReMatchCheckNeigh ; //! Time distribution of clusters, not matched to trigger, rematched with neighbour patchs
+ TH2F * fhTimeTriggerEMCALBC0UMReMatchBoth ; //! Time distribution of clusters, not matched to trigger, rematched open both
+
+ TH2F * fhEtaPhiNoTrigger ; //! Pseudorapidity vs Phi of highest E exotic cluster in event, no trigger at all
+ TH2F * fhTimeNoTrigger ; //! Time distribution of highest energy exotic cluster in event, no trigger at all
+
+ TH1F * fhEPhoton ; //! Number of identified photon vs energy
+ TH1F * fhPtPhoton ; //! Number of identified photon vs transerse momentum
+ TH2F * fhPhiPhoton ; //! Azimuthal angle of identified photon vs transerse momentum
+ TH2F * fhEtaPhoton ; //! Pseudorapidity of identified photon vs transerse momentum
+ TH2F * fhEtaPhiPhoton ; //! Pseudorapidity vs Phi of identified photon for E > 0.5
+ TH2F * fhEtaPhi05Photon ; //! Pseudorapidity vs Phi of identified photon for E < 0.5
+ TH2F * fhEtaPhiPhotonEMCALBC0 ; //! Pseudorapidity vs Phi of identified photon for E > 0.5
+ TH2F * fhEtaPhiPhotonEMCALBC1 ; //! Pseudorapidity vs Phi of identified photon for E > 0.5
+ TH2F * fhEtaPhiPhotonEMCALBCN ; //! Pseudorapidity vs Phi of identified photon for E > 0.5
+ TH2F * fhEtaPhiPhotonTriggerEMCALBC[11]; //! Pseudorapidity vs Phi of photons for E > 0.5
+ TH2F * fhTimePhotonTriggerEMCALBC [11]; //! Time distribution of photons, when trigger is in a given BC
+ TH2F * fhTimePhotonTriggerEMCALBCPileUpSPD[11] ; //! Time distribution of photons, when trigger is in a given BC, tagged as pile-up SPD
+ TH2F * fhEtaPhiPhotonTriggerEMCALBCUM[11]; //! Pseudorapidity vs Phi of photons for E > 2, not matched to trigger
+ TH2F * fhTimePhotonTriggerEMCALBCUM [11]; //! Time distribution of photons, when trigger is in a given BC, not matched to trigger
+
+ TH2F * fhTimePhotonTriggerEMCALBC0UMReMatchOpenTime ; //! Time distribution of photons in event, when trigger is not found, rematched open time trigger
+ TH2F * fhTimePhotonTriggerEMCALBC0UMReMatchCheckNeigh ; //! Time distribution of photons in event, when trigger is not found, rematched with neigbour patchs
+ TH2F * fhTimePhotonTriggerEMCALBC0UMReMatchBoth ; //! Time distribution of photons in event, when trigger is not found, rematched open both
+
+ TH2F * fhPtCentralityPhoton ; //! centrality vs photon pT
+ TH2F * fhPtEventPlanePhoton ; //! event plane vs photon pT
//Shower shape
-
- TH2F * fhDispE; //! cluster dispersion vs E
- TH2F * fhLam0E; //! cluster lambda0 vs E
- TH2F * fhLam1E; //! cluster lambda1 vs E
-
- TH2F * fhDispETRD; //! cluster dispersion vs E, SM covered by TRD
- TH2F * fhLam0ETRD; //! cluster lambda0 vs E, SM covered by TRD
- TH2F * fhLam1ETRD; //! cluster lambda1 vs E, SM covered by TRD
-
- TH2F * fhDispETM; //! cluster dispersion vs E, cut on Track Matching residual
- TH2F * fhLam0ETM; //! cluster lambda0 vs E, cut on Track Matching residual
- TH2F * fhLam1ETM; //! cluster lambda1 vs E, cut on Track Matching residual
-
- TH2F * fhDispETMTRD; //! cluster dispersion vs E, SM covered by TRD, cut on Track Matching residual
- TH2F * fhLam0ETMTRD; //! cluster lambda0 vs E, SM covered by TRD, cut on Track Matching residual
- TH2F * fhLam1ETMTRD; //! cluster lambda1 vs E, SM covered by TRD, cut on Track Matching residual
-
-
- TH2F * fhNCellsLam0LowE; //! number of cells in cluster vs lambda0
- TH2F * fhNCellsLam1LowE; //! number of cells in cluster vs lambda1
- TH2F * fhNCellsDispLowE; //! number of cells in cluster vs dispersion
- TH2F * fhNCellsLam0HighE; //! number of cells in cluster vs lambda0, E>2
- TH2F * fhNCellsLam1HighE; //! number of cells in cluster vs lambda1, E>2
- TH2F * fhNCellsDispHighE; //! number of cells in cluster vs dispersion, E>2
-
- TH2F * fhEtaLam0LowE; //! cluster eta vs lambda0, E<2
- TH2F * fhPhiLam0LowE; //! cluster phi vs lambda0, E<2
- TH2F * fhEtaLam0HighE; //! cluster eta vs lambda0, E>2
- TH2F * fhPhiLam0HighE; //! cluster phi vs lambda0, E>2
- TH2F * fhLam0DispLowE; //! cluster lambda0 vs dispersion, E<2
- TH2F * fhLam0DispHighE; //! cluster lambda0 vs dispersion, E>2
- TH2F * fhLam1Lam0LowE; //! cluster lambda1 vs lambda0, E<2
- TH2F * fhLam1Lam0HighE; //! cluster lambda1 vs lambda0, E>2
- TH2F * fhDispLam1LowE; //! cluster disp vs lambda1, E<2
- TH2F * fhDispLam1HighE; //! cluster disp vs lambda1, E>2
+ TH2F * fhNLocMax; //! number of maxima in selected clusters
+
+ TH2F * fhDispE; //! cluster dispersion vs E
+ TH2F * fhLam0E; //! cluster lambda0 vs E
+ TH2F * fhLam1E; //! cluster lambda1 vs E
+
+ TH2F * fhDispETRD; //! cluster dispersion vs E, SM covered by TRD
+ TH2F * fhLam0ETRD; //! cluster lambda0 vs E, SM covered by TRD
+ TH2F * fhLam1ETRD; //! cluster lambda1 vs E, SM covered by TRD
+
+ TH2F * fhDispETM; //! cluster dispersion vs E, cut on Track Matching residual
+ TH2F * fhLam0ETM; //! cluster lambda0 vs E, cut on Track Matching residual
+ TH2F * fhLam1ETM; //! cluster lambda1 vs E, cut on Track Matching residual
+
+ TH2F * fhDispETMTRD; //! cluster dispersion vs E, SM covered by TRD, cut on Track Matching residual
+ TH2F * fhLam0ETMTRD; //! cluster lambda0 vs E, SM covered by TRD, cut on Track Matching residual
+ TH2F * fhLam1ETMTRD; //! cluster lambda1 vs E, SM covered by TRD, cut on Track Matching residual
+
+ TH2F * fhNCellsLam0LowE; //! number of cells in cluster vs lambda0
+ TH2F * fhNCellsLam1LowE; //! number of cells in cluster vs lambda1
+ TH2F * fhNCellsDispLowE; //! number of cells in cluster vs dispersion
+ TH2F * fhNCellsLam0HighE; //! number of cells in cluster vs lambda0, E>2
+ TH2F * fhNCellsLam1HighE; //! number of cells in cluster vs lambda1, E>2
+ TH2F * fhNCellsDispHighE; //! number of cells in cluster vs dispersion, E>2
+
+ TH2F * fhEtaLam0LowE; //! cluster eta vs lambda0, E<2
+ TH2F * fhPhiLam0LowE; //! cluster phi vs lambda0, E<2
+ TH2F * fhEtaLam0HighE; //! cluster eta vs lambda0, E>2
+ TH2F * fhPhiLam0HighE; //! cluster phi vs lambda0, E>2
+ TH2F * fhLam0DispLowE; //! cluster lambda0 vs dispersion, E<2
+ TH2F * fhLam0DispHighE; //! cluster lambda0 vs dispersion, E>2
+ TH2F * fhLam1Lam0LowE; //! cluster lambda1 vs lambda0, E<2
+ TH2F * fhLam1Lam0HighE; //! cluster lambda1 vs lambda0, E>2
+ TH2F * fhDispLam1LowE; //! cluster disp vs lambda1, E<2
+ TH2F * fhDispLam1HighE; //! cluster disp vs lambda1, E>2
+ TH2F * fhDispEtaE ; //! shower dispersion in eta direction
+ TH2F * fhDispPhiE ; //! shower dispersion in phi direction
+ TH2F * fhSumEtaE ; //! shower dispersion in eta direction
+ TH2F * fhSumPhiE ; //! shower dispersion in phi direction
+ TH2F * fhSumEtaPhiE ; //! shower dispersion in eta and phi direction
+ TH2F * fhDispEtaPhiDiffE ; //! shower dispersion eta - phi
+ TH2F * fhSphericityE ; //! shower sphericity in eta vs phi
+ TH2F * fhDispSumEtaDiffE ; //! difference of 2 eta dispersions
+ TH2F * fhDispSumPhiDiffE ; //! difference of 2 phi dispersions
+ TH2F * fhDispEtaDispPhi[7] ; //! shower dispersion in eta direction vs phi direction for 5 E bins [0-2],[2-4],[4-6],[6-10],[> 10]
+ TH2F * fhLambda0DispEta[7] ; //! shower shape correlation l0 vs disp eta
+ TH2F * fhLambda0DispPhi[7] ; //! shower shape correlation l0 vs disp phi
+
//Fill MC dependent histograms, Origin of this cluster is ...
- TH2F * fhMCDeltaE[14] ; //! MC-Reco E distribution coming from MC particle
- TH2F * fhMCDeltaPt[14] ; //! MC-Reco pT distribution coming from MC particle
- TH2F * fhMC2E[14] ; //! E distribution, Reco vs MC coming from MC particle
- TH2F * fhMC2Pt[14] ; //! pT distribution, Reco vs MC coming from MC particle
-
- TH1F * fhMCE[14]; //! Number of identified photon vs cluster energy coming from MC particle
- TH1F * fhMCPt[14]; //! Number of identified photon vs cluster pT coming from MC particle
- TH2F * fhMCPhi[14]; //! Phi of identified photon coming from MC particle
- TH2F * fhMCEta[14]; //! eta of identified photon coming from MC particle
+ TH2F * fhMCDeltaE[14] ; //! MC-Reco E distribution coming from MC particle
+ TH2F * fhMCDeltaPt[14] ; //! MC-Reco pT distribution coming from MC particle
+ TH2F * fhMC2E[14] ; //! E distribution, Reco vs MC coming from MC particle
+ TH2F * fhMC2Pt[14] ; //! pT distribution, Reco vs MC coming from MC particle
+
+ TH1F * fhMCE[14]; //! Number of identified photon vs cluster energy coming from MC particle
+ TH1F * fhMCPt[14]; //! Number of identified photon vs cluster pT coming from MC particle
+ TH2F * fhMCPhi[14]; //! Phi of identified photon coming from MC particle
+ TH2F * fhMCEta[14]; //! eta of identified photon coming from MC particle
- TH1F * fhEPrimMC[7]; //! Number of generated photon vs energy
- TH1F * fhPtPrimMC[7]; //! Number of generated photon vs pT
- TH2F * fhPhiPrimMC[7]; //! Phi of generted photon
- TH2F * fhYPrimMC[7]; //! Rapidity of generated photon
-
- TH1F * fhEPrimMCAcc[7]; //! Number of generated photon vs energy, in calorimeter acceptance
- TH1F * fhPtPrimMCAcc[7]; //! Number of generated photon vs pT, in calorimeter acceptance
- TH2F * fhPhiPrimMCAcc[7]; //! Phi of generted photon, in calorimeter acceptance
- TH2F * fhYPrimMCAcc[7]; //! Rapidity of generated photon, in calorimeter acceptance
+ TH1F * fhEPrimMC[7]; //! Number of generated photon vs energy
+ TH1F * fhPtPrimMC[7]; //! Number of generated photon vs pT
+ TH2F * fhPhiPrimMC[7]; //! Phi of generted photon
+ TH2F * fhYPrimMC[7]; //! Rapidity of generated photon
+ TH2F * fhEtaPrimMC[7]; //! Eta of generated photon
+
+ TH1F * fhEPrimMCAcc[7]; //! Number of generated photon vs energy, in calorimeter acceptance
+ TH1F * fhPtPrimMCAcc[7]; //! Number of generated photon vs pT, in calorimeter acceptance
+ TH2F * fhPhiPrimMCAcc[7]; //! Phi of generted photon, in calorimeter acceptance
+ TH2F * fhEtaPrimMCAcc[7]; //! Phi of generted photon, in calorimeter acceptance
+ TH2F * fhYPrimMCAcc[7]; //! Rapidity of generated photon, in calorimeter acceptance
// Shower Shape MC
- TH2F * fhMCELambda0[6] ; //! E vs Lambda0 from MC particle
- TH2F * fhMCELambda1[6] ; //! E vs Lambda1 from MC particle
- TH2F * fhMCEDispersion[6] ; //! E vs Dispersion from MC particle
-
- TH2F * fhMCPhotonELambda0NoOverlap ; //! E vs Lambda0 from MC photons, no overlap
- TH2F * fhMCPhotonELambda0TwoOverlap ; //! E vs Lambda0 from MC photons, 2 particles overlap
- TH2F * fhMCPhotonELambda0NOverlap ; //! E vs Lambda0 from MC photons, N particles overlap
-
- TH2F * fhMCLambda0vsClusterMaxCellDiffE0[6]; //! Lambda0 vs fraction of energy of max cell for E < 2 GeV
- TH2F * fhMCLambda0vsClusterMaxCellDiffE2[6]; //! Lambda0 vs fraction of energy of max cell for 2< E < 6 GeV
- TH2F * fhMCLambda0vsClusterMaxCellDiffE6[6]; //! Lambda0 vs fraction of energy of max cell for E > 6 GeV
- TH2F * fhMCNCellsvsClusterMaxCellDiffE0[6]; //! NCells vs fraction of energy of max cell for E < 2
- TH2F * fhMCNCellsvsClusterMaxCellDiffE2[6]; //! NCells vs fraction of energy of max cell for 2 < E < 6 GeV
- TH2F * fhMCNCellsvsClusterMaxCellDiffE6[6]; //! NCells vs fraction of energy of max cell for E > 6
- TH2F * fhMCNCellsE[6]; //! NCells per cluster vs energy
- TH2F * fhMCMaxCellDiffClusterE[6]; //! Fraction of energy carried by cell with maximum energy
+ TH2F * fhMCELambda0[6] ; //! E vs Lambda0 from MC particle
+ TH2F * fhMCELambda1[6] ; //! E vs Lambda1 from MC particle
+ TH2F * fhMCEDispersion[6] ; //! E vs Dispersion from MC particle
+
+ TH2F * fhMCPhotonELambda0NoOverlap ; //! E vs Lambda0 from MC photons, no overlap
+ TH2F * fhMCPhotonELambda0TwoOverlap ; //! E vs Lambda0 from MC photons, 2 particles overlap
+ TH2F * fhMCPhotonELambda0NOverlap ; //! E vs Lambda0 from MC photons, N particles overlap
+
+ TH2F * fhMCLambda0vsClusterMaxCellDiffE0[6]; //! Lambda0 vs fraction of energy of max cell for E < 2 GeV
+ TH2F * fhMCLambda0vsClusterMaxCellDiffE2[6]; //! Lambda0 vs fraction of energy of max cell for 2< E < 6 GeV
+ TH2F * fhMCLambda0vsClusterMaxCellDiffE6[6]; //! Lambda0 vs fraction of energy of max cell for E > 6 GeV
+ TH2F * fhMCNCellsvsClusterMaxCellDiffE0[6]; //! NCells vs fraction of energy of max cell for E < 2
+ TH2F * fhMCNCellsvsClusterMaxCellDiffE2[6]; //! NCells vs fraction of energy of max cell for 2 < E < 6 GeV
+ TH2F * fhMCNCellsvsClusterMaxCellDiffE6[6]; //! NCells vs fraction of energy of max cell for E > 6
+ TH2F * fhMCNCellsE[6]; //! NCells per cluster vs energy
+ TH2F * fhMCMaxCellDiffClusterE[6]; //! Fraction of energy carried by cell with maximum energy
+
+ TH2F * fhMCEDispEta[6] ; //! shower dispersion in eta direction
+ TH2F * fhMCEDispPhi[6] ; //! shower dispersion in phi direction
+ TH2F * fhMCESumEtaPhi[6] ; //! shower dispersion in eta vs phi direction
+ TH2F * fhMCEDispEtaPhiDiff[6] ; //! shower dispersion in eta -phi direction
+ TH2F * fhMCESphericity[6] ; //! shower sphericity, eta vs phi
+ TH2F * fhMCDispEtaDispPhi[7][6] ; //! shower dispersion in eta direction vs phi direction for 5 E bins [0-2],[2-4],[4-6],[6-10],[> 10]
+ TH2F * fhMCLambda0DispEta[7][6] ; //! shower shape correlation l0 vs disp eta
+ TH2F * fhMCLambda0DispPhi[7][6] ; //! shower shape correlation l0 vs disp phi
//Embedding
- TH2F * fhEmbeddedSignalFractionEnergy ; //! Fraction of photon energy of embedded signal vs cluster energy
+ TH2F * fhEmbeddedSignalFractionEnergy ; //! Fraction of photon energy of embedded signal vs cluster energy
- TH2F * fhEmbedPhotonELambda0FullSignal ; //! Lambda0 vs E for embedded photons with more than 90% of the cluster energy
- TH2F * fhEmbedPhotonELambda0MostlySignal ; //! Lambda0 vs E for embedded photons with 90%<fraction<50%
- TH2F * fhEmbedPhotonELambda0MostlyBkg ; //! Lambda0 vs E for embedded photons with 50%<fraction<10%
- TH2F * fhEmbedPhotonELambda0FullBkg ; //! Lambda0 vs E for embedded photons with less than 10% of the cluster energy
+ TH2F * fhEmbedPhotonELambda0FullSignal ; //! Lambda0 vs E for embedded photons with more than 90% of the cluster energy
+ TH2F * fhEmbedPhotonELambda0MostlySignal ; //! Lambda0 vs E for embedded photons with 90%<fraction<50%
+ TH2F * fhEmbedPhotonELambda0MostlyBkg ; //! Lambda0 vs E for embedded photons with 50%<fraction<10%
+ TH2F * fhEmbedPhotonELambda0FullBkg ; //! Lambda0 vs E for embedded photons with less than 10% of the cluster energy
- TH2F * fhEmbedPi0ELambda0FullSignal ; //! Lambda0 vs E for embedded photons with more than 90% of the cluster energy
- TH2F * fhEmbedPi0ELambda0MostlySignal ; //! Lambda0 vs E for embedded photons with 90%<fraction<50%
- TH2F * fhEmbedPi0ELambda0MostlyBkg ; //! Lambda0 vs E for embedded photons with 50%<fraction<10%
- TH2F * fhEmbedPi0ELambda0FullBkg ; //! Lambda0 vs E for embedded photons with less than 10% of the cluster energy
+ TH2F * fhEmbedPi0ELambda0FullSignal ; //! Lambda0 vs E for embedded photons with more than 90% of the cluster energy
+ TH2F * fhEmbedPi0ELambda0MostlySignal ; //! Lambda0 vs E for embedded photons with 90%<fraction<50%
+ TH2F * fhEmbedPi0ELambda0MostlyBkg ; //! Lambda0 vs E for embedded photons with 50%<fraction<10%
+ TH2F * fhEmbedPi0ELambda0FullBkg ; //! Lambda0 vs E for embedded photons with less than 10% of the cluster energy
// Track Matching
- TH2F * fhTrackMatchedDEta[2] ; //! Eta distance between track and cluster vs cluster E, after and before photon cuts
- TH2F * fhTrackMatchedDPhi[2] ; //! Phi distance between track and cluster vs cluster E, after and before photon cuts
- TH2F * fhTrackMatchedDEtaDPhi[2] ; //! Eta vs Phi distance between track and cluster, E cluster > 0.5 GeV, after and before photon cuts
-
- TH2F * fhTrackMatchedDEtaTRD[2] ; //! Eta distance between track and cluster vs cluster E, after and before photon cuts, behind TRD
- TH2F * fhTrackMatchedDPhiTRD[2] ; //! Phi distance between track and cluster vs cluster E, after and before photon cuts, behind TRD
-
- TH2F * fhTrackMatchedDEtaMCOverlap[2] ; //! Eta distance between track and cluster vs cluster E, several particle overlap, after and before photon cuts
- TH2F * fhTrackMatchedDPhiMCOverlap[2] ; //! Phi distance between track and cluster vs cluster E, several particle overlap, after and before photon cuts
- TH2F * fhTrackMatchedDEtaMCNoOverlap[2]; //! Eta distance between track and cluster vs cluster E, not other particle overlap, after and before photon cuts
- TH2F * fhTrackMatchedDPhiMCNoOverlap[2]; //! Phi distance between track and cluster vs cluster E, not other particle overlap, after and before photon cuts
- TH2F * fhTrackMatchedDEtaMCConversion[2]; //! Eta distance between track and cluster vs cluster E, originated in conversion, after and before photon cuts
- TH2F * fhTrackMatchedDPhiMCConversion[2]; //! Phi distance between track and cluster vs cluster E, originated in conversion, after and before photon cuts
-
- TH2F * fhTrackMatchedMCParticle[2]; //! Trace origin of matched particle
- TH2F * fhdEdx[2]; //! matched track dEdx vs cluster E, after and before photon cuts
- TH2F * fhEOverP[2]; //! matched track E cluster over P track vs cluster E, after dEdx cut, after and before photon cuts
- TH2F * fhEOverPTRD[2]; //! matched track E cluster over P track vs cluster E, after dEdx cut, after and before photon cuts, behind TRD
+ TH2F * fhTrackMatchedDEta[2] ; //! Eta distance between track and cluster vs cluster E, after and before photon cuts
+ TH2F * fhTrackMatchedDPhi[2] ; //! Phi distance between track and cluster vs cluster E, after and before photon cuts
+ TH2F * fhTrackMatchedDEtaDPhi[2] ; //! Eta vs Phi distance between track and cluster, E cluster > 0.5 GeV, after and before
+
+ TH2F * fhTrackMatchedDEtaPos[2] ; //! Eta distance between track and cluster vs cluster E, after and before photon cuts
+ TH2F * fhTrackMatchedDPhiPos[2] ; //! Phi distance between track and cluster vs cluster E, after and before photon cuts
+ TH2F * fhTrackMatchedDEtaDPhiPos[2] ; //! Eta vs Phi distance between track and cluster, E cluster > 0.5 GeV, after and before
+
+ TH2F * fhTrackMatchedDEtaNeg[2] ; //! Eta distance between track and cluster vs cluster E, after and before photon cuts
+ TH2F * fhTrackMatchedDPhiNeg[2] ; //! Phi distance between track and cluster vs cluster E, after and before photon cuts
+ TH2F * fhTrackMatchedDEtaDPhiNeg[2] ; //! Eta vs Phi distance between track and cluster, E cluster > 0.5 GeV, after and before photon cuts
+
+ TH2F * fhTrackMatchedDEtaTRD[2] ; //! Eta distance between track and cluster vs cluster E, after and before photon cuts, behind TRD
+ TH2F * fhTrackMatchedDPhiTRD[2] ; //! Phi distance between track and cluster vs cluster E, after and before photon cuts, behind TRD
+
+ TH2F * fhTrackMatchedDEtaMCOverlap[2] ; //! Eta distance between track and cluster vs cluster E, several particle overlap, after and before photon cuts
+ TH2F * fhTrackMatchedDPhiMCOverlap[2] ; //! Phi distance between track and cluster vs cluster E, several particle overlap, after and before photon cuts
+ TH2F * fhTrackMatchedDEtaMCNoOverlap[2]; //! Eta distance between track and cluster vs cluster E, not other particle overlap, after and before photon cuts
+ TH2F * fhTrackMatchedDPhiMCNoOverlap[2]; //! Phi distance between track and cluster vs cluster E, not other particle overlap, after and before photon cuts
+ TH2F * fhTrackMatchedDEtaMCConversion[2]; //! Eta distance between track and cluster vs cluster E, originated in conversion, after and before photon cuts
+ TH2F * fhTrackMatchedDPhiMCConversion[2]; //! Phi distance between track and cluster vs cluster E, originated in conversion, after and before photon cuts
+
+ TH2F * fhTrackMatchedMCParticle[2]; //! Trace origin of matched particle
+ TH2F * fhdEdx[2]; //! matched track dEdx vs cluster E, after and before photon cuts
+ TH2F * fhEOverP[2]; //! matched track E cluster over P track vs cluster E, after dEdx cut, after and before photon cuts
+ TH2F * fhEOverPTRD[2]; //! matched track E cluster over P track vs cluster E, after dEdx cut, after and before photon cuts, behind TRD
+ // Pile-up
+ TH1F * fhPtPileUp[7]; //! pT distribution of clusters before any selection
+ TH1F * fhPtChargedPileUp[7]; //! pT distribution of track matched clusters
+ TH1F * fhPtPhotonPileUp[7]; //! pT distribution of selected photons
+ TH2F * fhLambda0PileUp[7]; //! E vs M02 distribution of clusters, before any selection
+ TH2F * fhLambda0ChargedPileUp[7]; //! E vs M02 distribution of clusters, track matched clusters
+ TH2F * fhClusterCellTimePileUp[7]; //! E vs Time inside cluster, before any selection, not max cell
+ TH2F * fhClusterTimeDiffPileUp[7]; //! E vs Time difference inside cluster, before any selection
+ TH2F * fhClusterTimeDiffChargedPileUp[7]; //! E vs Time difference inside cluster for track matched clusters
+ TH2F * fhClusterTimeDiffPhotonPileUp[7]; //! E vs Time difference inside cluster for selected photons
+ TH2F * fhClusterEFracLongTimePileUp[7]; //! E vs fraction of cluster energy from cells with large time
+ TH2F * fhTimePtNoCut; //! time of cluster vs Pt, no cut
+ TH2F * fhTimePtSPD; //! time of cluster vs Pt, IsSPDPileUp
+ TH2F * fhTimePtPhotonNoCut; //! time of photon cluster vs Pt, no cut
+ TH2F * fhTimePtPhotonSPD; //! time of photon cluster vs Pt, IsSPDPileUp
+ TH2F * fhTimeNPileUpVertSPD; //! time of cluster vs n pile-up vertices from SPD
+ TH2F * fhTimeNPileUpVertTrack; //! time of cluster vs n pile-up vertices from Tracks
+ TH2F * fhTimeNPileUpVertContributors; //! time of cluster vs n pile-up vertex from SPD contributors
+ 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
+ TH2F * fhClusterMultSPDPileUp[4]; //! E max cluster vs event cluster multiplicity, for tmax-tdiff cuts, pile up event
+ TH2F * fhClusterMultNoPileUp[4]; //! E max cluster vs event cluster multiplicity, for tmax-tdiff cuts, not pile up event
+ TH2F * fhEtaPhiBC0; //! eta/phi of clusters in BC=0
+ TH2F * fhEtaPhiBCPlus; //! eta/phi of clusters in BC>0
+ TH2F * fhEtaPhiBCMinus; //! eta/phi of clusters in BC<0
+ TH2F * fhEtaPhiBC0PileUpSPD; //! eta/phi of clusters in BC=0, SPD pile-up
+ TH2F * fhEtaPhiBCPlusPileUpSPD; //! eta/phi of clusters in BC>0, SPD pile-up
+ TH2F * fhEtaPhiBCMinusPileUpSPD; //! eta/phi of clusters in BC<0, SPD pile-up
+
+ TH2F * fhPtNPileUpSPDVtx; //! cluster pt vs number of spd pile-up vertices
+ TH2F * fhPtNPileUpTrkVtx; //! cluster pt vs number of track pile-up vertices
+ TH2F * fhPtNPileUpSPDVtxTimeCut; //! cluster pt vs number of spd pile-up vertices, time cut +-25 ns
+ TH2F * fhPtNPileUpTrkVtxTimeCut; //! cluster pt vs number of track pile-up vertices, time cut +- 25 ns
+ TH2F * fhPtNPileUpSPDVtxTimeCut2; //! cluster pt vs number of spd pile-up vertices, time cut +-75 ns
+ TH2F * fhPtNPileUpTrkVtxTimeCut2; //! cluster pt vs number of track pile-up vertices, time cut +- 75 ns
+
+ TH2F * fhPtPhotonNPileUpSPDVtx; //! photon pt vs number of spd pile-up vertices
+ TH2F * fhPtPhotonNPileUpTrkVtx; //! photon pt vs number of track pile-up vertices
+ TH2F * fhPtPhotonNPileUpSPDVtxTimeCut; //! photon pt vs number of spd pile-up vertices, time cut +-25 ns
+ TH2F * fhPtPhotonNPileUpTrkVtxTimeCut; //! photon pt vs number of track pile-up vertices, time cut +- 25 ns
+ TH2F * fhPtPhotonNPileUpSPDVtxTimeCut2; //! photon pt vs number of spd pile-up vertices, time cut +-75 ns
+ TH2F * fhPtPhotonNPileUpTrkVtxTimeCut2; //! photon pt vs number of track pile-up vertices, time cut +- 75 ns
+
+ TH2F * fhEClusterSM ; //! cluster E distribution per SM, before any selection, after reader
+ TH2F * fhEPhotonSM ; //! photon-like cluster E distribution per SM
+ TH2F * fhPtClusterSM; //! cluster E distribution per SM, before any selection, after reader
+ TH2F * fhPtPhotonSM ; //! photon-like cluster E distribution per SM
+
AliAnaPhoton( const AliAnaPhoton & g) ; // cpy ctor
AliAnaPhoton & operator = (const AliAnaPhoton & g) ; // cpy assignment
- ClassDef(AliAnaPhoton,23)
+ ClassDef(AliAnaPhoton,36)
} ;
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff