PWGGA/CaloTrackCorrelations/AliAnaPhoton.h

   1 #ifndef ALIANAPHOTON_H
   2 #define ALIANAPHOTON_H
   3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
   4  * See cxx source for full Copyright notice     */
   5
   6 //_________________________________________________________________________
   7 //
   8 // Class for the photon identification.
   9 // Clusters from calorimeters are identified as photons
  10 // and kept in the AOD. Few histograms produced.
  11 // Produces input for other analysis classes like AliAnaPi0,
  12 // AliAnaParticleHadronCorrelation ...
  13 //
  14
  15 //-- Author: Gustavo Conesa (INFN-LNF)
  16
  17 // --- ROOT system ---
  18 class TH2F ;
  19 class TH1F;
  20 class TString ;
  21 class TObjString;
  22 class TList ;
  23
  24 // --- ANALYSIS system ---
  25 #include "AliAnaCaloTrackCorrBaseClass.h"
  26
  27 class AliAnaPhoton : public AliAnaCaloTrackCorrBaseClass {
  28
  29  public:
  30   AliAnaPhoton() ;              // default ctor
  31   virtual ~AliAnaPhoton() { ; } // virtual dtor
  32
  33   //---------------------------------------
  34   // General analysis frame methods
  35   //---------------------------------------
  36
  37   TObjString * GetAnalysisCuts();
  38
  39   TList      * GetCreateOutputObjects();
  40
  41   void         Init();
  42
  43   void         InitParameters();
  44
  45   void         MakeAnalysisFillAOD()  ;
  46
  47   void         MakeAnalysisFillHistograms() ;
  48
  49   void         Print(const Option_t * opt)const;
  50
  51
  52   // Analysis methods
  53
  54   Bool_t       ClusterSelected(AliVCluster* cl, TLorentzVector mom, Int_t nlm) ;
  55
  56   void         FillAcceptanceHistograms();
  57
  58   void         FillEMCALTriggerClusterBCHistograms(const Int_t idcalo,       const Float_t ecluster,   const Float_t tofcluster,
  59                                                    const Float_t etacluster, const Float_t phicluster);
  60
  61   void         FillShowerShapeHistograms( AliVCluster* cluster, Int_t mcTag) ;
  62
  63   void         SwitchOnFillShowerShapeHistograms()    { fFillSSHistograms      = kTRUE  ; }
  64   void         SwitchOffFillShowerShapeHistograms()   { fFillSSHistograms      = kFALSE ; }
  65
  66   void         SwitchOnOnlySimpleSSHistoFill()        { fFillOnlySimpleSSHisto = kTRUE  ; }
  67   void         SwitchOffOnlySimpleHistoFill()         { fFillOnlySimpleSSHisto = kFALSE ; }
  68
  69   void         FillTrackMatchingResidualHistograms(AliVCluster* calo, Int_t cut);
  70
  71   void         SwitchOnTMHistoFill()                  { fFillTMHisto           = kTRUE  ; }
  72   void         SwitchOffTMHistoFill()                 { fFillTMHisto           = kFALSE ; }
  73
  74   void         FillClusterPileUpHistograms(AliVCluster * calo,       const Bool_t matched,
  75                                            const Float_t ecluster,   const Float_t ptcluster,
  76                                            const Float_t etacluster, const Float_t phicluster,
  77                                            const Float_t l0cluster);
  78
  79   void         FillPileUpHistograms(Float_t energy, Float_t pt, Float_t time) ;
  80   void         FillPileUpHistogramsPerEvent() ;
  81
  82   void         SwitchOnFillPileUpHistograms()         { fFillPileUpHistograms  = kTRUE  ; }
  83   void         SwitchOffFillPileUpHistograms()        { fFillPileUpHistograms  = kFALSE ; }
  84
  85   void         SwitchOnFillEMCALBCHistograms()        { fFillEMCALBCHistograms = kTRUE  ; }
  86   void         SwitchOffFillEMCALBCHistograms()       { fFillEMCALBCHistograms = kFALSE ; }
  87
  88
  89   // Analysis parameters setters getters
  90
  91   TString      GetCalorimeter()                 const { return fCalorimeter        ; }
  92   void         SetCalorimeter(TString  & det)         { fCalorimeter = det         ; }
  93
  94   // ** Cluster selection methods **
  95
  96   void         SetMinDistanceToBadChannel(Float_t m1, Float_t m2, Float_t m3) {
  97                 fMinDist = m1; fMinDist2 = m2; fMinDist3 = m3; }
  98
  99   void         SetTimeCut(Double_t min, Double_t max) { fTimeCutMin = min;
 100                                                         fTimeCutMax = max          ; }
 101   Double_t     GetTimeCutMin()                  const { return fTimeCutMin         ; }
 102   Double_t     GetTimeCutMax()                  const { return fTimeCutMax         ; }
 103
 104   void         SetNCellCut(Int_t n)                   { fNCellsCut = n             ; }
 105   Double_t     GetNCellCut()                    const { return fNCellsCut          ; }
 106
 107   void         SetNLMCut(Int_t min, Int_t max)        { fNLMCutMin = min;
 108     fNLMCutMax = max                ; }
 109   Int_t        GetNLMCutMin()                   const { return fNLMCutMin          ; }
 110   Int_t        GetNLMCutMax()                   const { return fNLMCutMax          ; }
 111
 112
 113   Bool_t       IsTrackMatchRejectionOn()        const { return fRejectTrackMatch   ; }
 114   void         SwitchOnTrackMatchRejection()          { fRejectTrackMatch = kTRUE  ; }
 115   void         SwitchOffTrackMatchRejection()         { fRejectTrackMatch = kFALSE ; }
 116
 117   void         FillNOriginHistograms(Int_t n)         { fNOriginHistograms = n ;
 118     if(n > 14) fNOriginHistograms = 14; }
 119   void         FillNPrimaryHistograms(Int_t n)        { fNPrimaryHistograms= n ;
 120     if(n > 7)  fNPrimaryHistograms = 7; }
 121
 122   // For histograms in arrays, index in the array, corresponding to a particle
 123   enum mcTypes    { kmcPhoton = 0,        kmcPi0Decay = 1,       kmcOtherDecay = 2,
 124                     kmcPi0 = 3,           kmcEta = 4,            kmcElectron = 5,
 125                     kmcConversion = 6,    kmcOther = 7,          kmcAntiNeutron = 8,
 126                     kmcAntiProton = 9,    kmcPrompt = 10,        kmcFragmentation = 11,
 127                     kmcISR = 12,          kmcString = 13                               };
 128
 129   enum mcPTypes   { kmcPPhoton = 0,       kmcPPi0Decay = 1,       kmcPOtherDecay = 2,  kmcPOther = 3,
 130                     kmcPPrompt = 4,       kmcPFragmentation = 5,  kmcPISR = 6           };
 131
 132   enum mcssTypes  { kmcssPhoton = 0,      kmcssOther = 1,       kmcssPi0 = 2,
 133                     kmcssEta = 3,         kmcssConversion = 4,  kmcssElectron = 5       };
 134
 135   private:
 136
 137   TString  fCalorimeter ;                // Calorimeter where the gamma is searched;
 138   Float_t  fMinDist ;                    // Minimal distance to bad channel to accept cluster
 139   Float_t  fMinDist2;                    // Cuts on Minimal distance to study acceptance evaluation
 140   Float_t  fMinDist3;                    // One more cut on distance used for acceptance-efficiency study
 141   Bool_t   fRejectTrackMatch ;           // If PID on, reject clusters which have an associated TPC track
 142   Bool_t   fFillTMHisto;                 // Fill track matching plots
 143   Double_t fTimeCutMin  ;                // Remove clusters/cells with time smaller than this value, in ns
 144   Double_t fTimeCutMax  ;                // Remove clusters/cells with time larger than this value, in ns
 145   Int_t    fNCellsCut ;                  // Accept for the analysis clusters with more than fNCellsCut cells
 146   Int_t    fNLMCutMin  ;                 // Remove clusters/cells with number of local maxima smaller than this value
 147   Int_t    fNLMCutMax  ;                 // Remove clusters/cells with number of local maxima larger than this value
 148   Bool_t   fFillSSHistograms ;           // Fill shower shape histograms
 149   Bool_t   fFillOnlySimpleSSHisto;       // Fill selected cluster histograms, selected SS histograms
 150   Int_t    fNOriginHistograms;           // Fill only NOriginHistograms of the 14 defined types
 151   Int_t    fNPrimaryHistograms;          // Fill only NPrimaryHistograms of the 7 defined types
 152   Bool_t   fFillPileUpHistograms;        // Fill pile-up related histograms
 153   Bool_t   fFillEMCALBCHistograms;       // Fill eta-phi BC dependent histograms
 154
 155   //Histograms
 156   TH1F * fhClusterCuts[10];              //! control histogram on the different photon selection cuts
 157   TH2F * fhNCellsE;                      //! number of cells in cluster vs E
 158   TH2F * fhCellsE;                       //! energy of cells in cluster vs E of cluster
 159   TH2F * fhMaxCellDiffClusterE;          //! Fraction of energy carried by cell with maximum energy
 160   TH2F * fhTimePt;                       //! time of photon cluster vs pt
 161
 162   TH2F * fhEtaPhi  ;                     //! Pseudorapidity vs Phi of clusters for E > 0.5
 163   TH2F * fhEtaPhiEMCALBC0  ;             //! Pseudorapidity vs Phi of clusters for E > 0.5
 164   TH2F * fhEtaPhiEMCALBC1  ;             //! Pseudorapidity vs Phi of clusters for E > 0.5
 165   TH2F * fhEtaPhiEMCALBCN  ;             //! Pseudorapidity vs Phi of clusters for E > 0.5
 166
 167   TH2F * fhEtaPhiTriggerEMCALBC[11] ;     //! Pseudorapidity vs Phi of clusters for E > 2
 168   TH2F * fhTimeTriggerEMCALBC  [11] ;     //! Time distribution of clusters, when trigger is in a given BC
 169   TH2F * fhTimeTriggerEMCALBCPileUpSPD[11]; //! Time distribution of clusters, when trigger is in a given BC, tagged as pile-up SPD
 170
 171   TH2F * fhEtaPhiTriggerEMCALBCUM[11] ;     //! Pseudorapidity vs Phi of clusters for E > 2, not matched to trigger
 172   TH2F * fhTimeTriggerEMCALBCUM  [11] ;     //! Time distribution of clusters, when trigger is in a given BC, not matched to trigger
 173
 174   TH2F * fhEtaPhiTriggerEMCALBCCluster  [11] ; //! Pseudorapidity vs Phi of trigger clusters
 175   TH2F * fhTimeTriggerEMCALBCCluster    [11] ; //! Time distribution of clusters, when trigger cluster is in a given BC
 176   TH2F * fhEtaPhiTriggerEMCALBCUMCluster[11] ; //! Pseudorapidity vs Phi of highest E cluster  in event, not matched to trigger
 177   TH2F * fhTimeTriggerEMCALBCUMCluster  [11] ; //! Time distribution of highest energy cluster in event, when trigger is in a given BC, not
 178
 179   TH2F * fhEtaPhiTriggerEMCALBCClusterOverTh     ; //! Pseudorapidity vs Phi of trigger clusters, over nominal threshold
 180   TH2F * fhEtaPhiTriggerEMCALBCUMClusterOverTh   ; //! Pseudorapidity vs Phi of highest E cluster  in event, not matched to trigger, over nominal threshold
 181   TH2F * fhEtaPhiTriggerEMCALBCClusterBelowTh1   ; //! Pseudorapidity vs Phi of trigger clusters, 1 GeV below nominal threshold
 182   TH2F * fhEtaPhiTriggerEMCALBCUMClusterBelowTh1 ; //! Pseudorapidity vs Phi of highest E cluster  in event, not matched to trigger, 2 GeV below nominal threshold
 183   TH2F * fhEtaPhiTriggerEMCALBCClusterBelowTh2   ; //! Pseudorapidity vs Phi of trigger clusters, 1 GeV below nominal threshold
 184   TH2F * fhEtaPhiTriggerEMCALBCUMClusterBelowTh2 ; //! Pseudorapidity vs Phi of highest E cluster  in event, not matched to trigger, 2 GeV below nominal threshold
 185
 186   TH2F * fhEtaPhiTriggerEMCALBCExotic            ; //! Pseudorapidity vs Phi of trigger exotic clusters
 187   TH2F * fhTimeTriggerEMCALBCExotic              ; //! Time distribution of clusters, when trigger exotic cluster
 188   TH2F * fhEtaPhiTriggerEMCALBCUMExotic          ; //! Pseudorapidity vs Phi of highest E exotic cluster  in event, not matched to trigger
 189   TH2F * fhTimeTriggerEMCALBCUMExotic            ; //! Time distribution of highest energy exotic cluster in event, not matched to trigger
 190
 191   TH2F * fhEtaPhiTriggerEMCALBCBad               ; //! Pseudorapidity vs Phi of trigger exotic clusters
 192   TH2F * fhTimeTriggerEMCALBCBad                 ; //! Time distribution of clusters, when trigger exotic
 193   TH2F * fhEtaPhiTriggerEMCALBCUMBad             ; //! Pseudorapidity vs Phi of highest E exotic cluster  in event, not matched to trigger
 194   TH2F * fhTimeTriggerEMCALBCUMBad               ; //! Time distribution of highest energy exotic cluster in event, not matched to trigger
 195
 196   TH2F * fhEtaPhiTriggerEMCALBCBadExotic         ; //! Pseudorapidity vs Phi of trigger exotic and bad clusters
 197   TH2F * fhTimeTriggerEMCALBCBadExotic           ; //! Time distribution of clusters, when trigger exotic and bad cluster
 198   TH2F * fhEtaPhiTriggerEMCALBCUMBadExotic       ; //! Pseudorapidity vs Phi of highest E exotic cluster  in event, not matched to trigger
 199   TH2F * fhTimeTriggerEMCALBCUMBadExotic         ; //! Time distribution of highest energy exotic cluster in event, not matched to trigger
 200
 201   TH2F * fhEtaPhiTriggerEMCALBCExoticCluster     ; //! Pseudorapidity vs Phi of trigger exotic clusters
 202   TH2F * fhTimeTriggerEMCALBCExoticCluster       ; //! Time distribution of clusters, when trigger exotic cluster
 203   TH2F * fhEtaPhiTriggerEMCALBCUMExoticCluster   ; //! Pseudorapidity vs Phi of highest E exotic cluster  in event, not matched to trigger
 204   TH2F * fhTimeTriggerEMCALBCUMExoticCluster     ; //! Time distribution of highest energy exotic cluster in event, not matched to trigger
 205
 206   TH2F * fhEtaPhiTriggerEMCALBCBadCluster        ; //! Pseudorapidity vs Phi of trigger bad clusters
 207   TH2F * fhTimeTriggerEMCALBCBadCluster          ; //! Time distribution of clusters, when trigger bad cluster is in a given BC
 208   TH2F * fhEtaPhiTriggerEMCALBCUMBadCluster      ; //! Pseudorapidity vs Phi of highest E bad cluster  in event, not matched to trigger
 209   TH2F * fhTimeTriggerEMCALBCUMBadCluster        ; //! Time distribution of highest energy bad cluster in event, when trigger is in a given BC, not
 210
 211   TH2F * fhEtaPhiTriggerEMCALBCBadExoticCluster  ; //! Pseudorapidity vs Phi of trigger exotic and bad clusters
 212   TH2F * fhTimeTriggerEMCALBCBadExoticCluster    ; //! Time distribution of clusters, when trigger exotic and bad cluster
 213   TH2F * fhEtaPhiTriggerEMCALBCUMBadExoticCluster; //! Pseudorapidity vs Phi of highest E exotic and bad cluster in event, not matched to trigger
 214   TH2F * fhTimeTriggerEMCALBCUMBadExoticCluster  ; //! Time distribution of highest energy exotic and bad cluster in event, not matched to trigger
 215
 216   TH2F * fhTimeTriggerEMCALBCBadMaxCell          ; //! Time distribution of trigger clusters, when trigger bad max cell
 217   TH2F * fhTimeTriggerEMCALBCUMBadMaxCell        ; //! Time distribution of highest energy bad max cell cluster in event, when trigger is not found
 218   TH2F * fhTimeTriggerEMCALBCBadMaxCellExotic    ; //! Time distribution of trigger clusters, when trigger exotic cluster with bad max cell
 219   TH2F * fhTimeTriggerEMCALBCUMBadMaxCellExotic  ; //! Time distribution of highest energy exotic with bad max cell cluster in event, when trigger is not found
 220
 221   TH2F * fhEtaPhiTriggerEMCALBCUMReMatchOpenTime ;  //! Pseudorapidity vs Phi of highest E bad cluster  in event, not matched to trigger
 222   TH2F * fhTimeTriggerEMCALBCUMReMatchOpenTime   ;  //! Time distribution of highest energy bad max cell cluster in event, when trigger is not found
 223   TH2F * fhEtaPhiTriggerEMCALBCUMReMatchCheckNeigh; //! Pseudorapidity vs Phi of highest E bad cluster  in event, not matched to trigger
 224   TH2F * fhTimeTriggerEMCALBCUMReMatchCheckNeigh ;  //! Time distribution of highest energy bad max cell cluster in event, when trigger is not found
 225   TH2F * fhEtaPhiTriggerEMCALBCUMReMatchBoth     ;  //! Pseudorapidity vs Phi of highest E bad cluster  in event, not matched to trigger
 226   TH2F * fhTimeTriggerEMCALBCUMReMatchBoth       ;  //! Time distribution of highest energy bad max cell cluster in event, when trigger is not found
 227
 228   TH2F * fhEtaPhiNoTrigger ;             //! Pseudorapidity vs Phi of highest E exotic cluster  in event, no trigger at all
 229   TH2F * fhTimeNoTrigger   ;             //! Time distribution of highest energy exotic cluster in event, no trigger at all
 230
 231   TH1F * fhEPhoton    ;                  //! Number of identified photon vs energy
 232   TH1F * fhPtPhoton   ;                  //! Number of identified photon vs transerse momentum
 233   TH2F * fhPhiPhoton  ;                  //! Azimuthal angle of identified  photon vs transerse momentum
 234   TH2F * fhEtaPhoton  ;                  //! Pseudorapidity of identified  photon vs transerse momentum
 235   TH2F * fhEtaPhiPhoton  ;               //! Pseudorapidity vs Phi of identified  photon for E > 0.5
 236   TH2F * fhEtaPhi05Photon  ;             //! Pseudorapidity vs Phi of identified  photon for E < 0.5
 237   TH2F * fhEtaPhiPhotonEMCALBC0  ;       //! Pseudorapidity vs Phi of identified  photon for E > 0.5
 238   TH2F * fhEtaPhiPhotonEMCALBC1  ;       //! Pseudorapidity vs Phi of identified  photon for E > 0.5
 239   TH2F * fhEtaPhiPhotonEMCALBCN  ;       //! Pseudorapidity vs Phi of identified  photon for E > 0.5
 240   TH2F * fhEtaPhiPhotonTriggerEMCALBC[11]; //! Pseudorapidity vs Phi of photons for E > 0.5
 241   TH2F * fhTimePhotonTriggerEMCALBC  [11]; //! Time distribution of photons, when trigger is in a given BC
 242   TH2F * fhTimePhotonTriggerEMCALBCPileUpSPD[11] ; //! Time distribution of photons, when trigger is in a given BC, tagged as pile-up SPD
 243   TH2F * fhEtaPhiPhotonTriggerEMCALBCUM[11]; //! Pseudorapidity vs Phi of photons for E > 2, not matched to trigger
 244   TH2F * fhTimePhotonTriggerEMCALBCUM  [11]; //! Time distribution of photons, when trigger is in a given BC, not matched to trigger
 245
 246   TH2F * fhPtCentralityPhoton    ;       //! centrality  vs photon pT
 247   TH2F * fhPtEventPlanePhoton    ;       //! event plane vs photon pT
 248
 249   //Shower shape
 250   TH2F * fhNLocMax;                       //! number of maxima in selected clusters
 251
 252   TH2F * fhDispE;                         //! cluster dispersion vs E
 253   TH2F * fhLam0E;                         //! cluster lambda0 vs  E
 254   TH2F * fhLam1E;                         //! cluster lambda1 vs  E
 255
 256   TH2F * fhDispETRD;                      //! cluster dispersion vs E, SM covered by TRD
 257   TH2F * fhLam0ETRD;                      //! cluster lambda0 vs  E, SM covered by TRD
 258   TH2F * fhLam1ETRD;                      //! cluster lambda1 vs  E, SM covered by TRD
 259
 260   TH2F * fhDispETM;                       //! cluster dispersion vs E, cut on Track Matching residual
 261   TH2F * fhLam0ETM;                       //! cluster lambda0 vs  E, cut on Track Matching residual
 262   TH2F * fhLam1ETM;                       //! cluster lambda1 vs  E, cut on Track Matching residual
 263
 264   TH2F * fhDispETMTRD;                    //! cluster dispersion vs E, SM covered by TRD, cut on Track Matching residual
 265   TH2F * fhLam0ETMTRD;                    //! cluster lambda0 vs  E, SM covered by TRD, cut on Track Matching residual
 266   TH2F * fhLam1ETMTRD;                    //! cluster lambda1 vs  E, SM covered by TRD, cut on Track Matching residual
 267
 268   TH2F * fhNCellsLam0LowE;                //! number of cells in cluster vs lambda0
 269   TH2F * fhNCellsLam1LowE;                //! number of cells in cluster vs lambda1
 270   TH2F * fhNCellsDispLowE;                //! number of cells in cluster vs dispersion
 271   TH2F * fhNCellsLam0HighE;               //! number of cells in cluster vs lambda0, E>2
 272   TH2F * fhNCellsLam1HighE;               //! number of cells in cluster vs lambda1, E>2
 273   TH2F * fhNCellsDispHighE;               //! number of cells in cluster vs dispersion, E>2
 274
 275   TH2F * fhEtaLam0LowE;                   //! cluster eta vs lambda0, E<2
 276   TH2F * fhPhiLam0LowE;                   //! cluster phi vs lambda0, E<2
 277   TH2F * fhEtaLam0HighE;                  //! cluster eta vs lambda0, E>2
 278   TH2F * fhPhiLam0HighE;                  //! cluster phi vs lambda0, E>2
 279   TH2F * fhLam0DispLowE;                  //! cluster lambda0 vs dispersion, E<2
 280   TH2F * fhLam0DispHighE;                 //! cluster lambda0 vs dispersion, E>2
 281   TH2F * fhLam1Lam0LowE;                  //! cluster lambda1 vs lambda0, E<2
 282   TH2F * fhLam1Lam0HighE;                 //! cluster lambda1 vs lambda0, E>2
 283   TH2F * fhDispLam1LowE;                  //! cluster disp vs lambda1, E<2
 284   TH2F * fhDispLam1HighE;                 //! cluster disp vs lambda1, E>2
 285
 286   TH2F * fhDispEtaE ;                     //! shower dispersion in eta direction
 287   TH2F * fhDispPhiE ;                     //! shower dispersion in phi direction
 288   TH2F * fhSumEtaE ;                      //! shower dispersion in eta direction
 289   TH2F * fhSumPhiE ;                      //! shower dispersion in phi direction
 290   TH2F * fhSumEtaPhiE ;                   //! shower dispersion in eta and phi direction
 291   TH2F * fhDispEtaPhiDiffE ;              //! shower dispersion eta - phi
 292   TH2F * fhSphericityE ;                  //! shower sphericity in eta vs phi
 293   TH2F * fhDispSumEtaDiffE ;              //! difference of 2 eta dispersions
 294   TH2F * fhDispSumPhiDiffE ;              //! difference of 2 phi dispersions
 295   TH2F * fhDispEtaDispPhi[7] ;            //! shower dispersion in eta direction vs phi direction for 5 E bins [0-2],[2-4],[4-6],[6-10],[> 10]
 296   TH2F * fhLambda0DispEta[7] ;            //! shower shape correlation l0 vs disp eta
 297   TH2F * fhLambda0DispPhi[7] ;            //! shower shape correlation l0 vs disp phi
 298
 299   //Fill MC dependent histograms, Origin of this cluster is ...
 300
 301   TH2F * fhMCDeltaE[14]  ;                      //! MC-Reco E distribution coming from MC particle
 302   TH2F * fhMCDeltaPt[14] ;                      //! MC-Reco pT distribution coming from MC particle
 303   TH2F * fhMC2E[14]  ;                          //! E distribution, Reco vs MC coming from MC particle
 304   TH2F * fhMC2Pt[14] ;                          //! pT distribution, Reco vs MC coming from MC particle
 305
 306   TH1F * fhMCE[14];                             //! Number of identified photon vs cluster energy coming from MC particle
 307   TH1F * fhMCPt[14];                            //! Number of identified photon vs cluster pT     coming from MC particle
 308   TH2F * fhMCPhi[14];                           //! Phi of identified photon coming from MC particle
 309   TH2F * fhMCEta[14];                           //! eta of identified photon coming from MC particle
 310
 311   TH1F * fhEPrimMC[7];                          //! Number of generated photon vs energy
 312   TH1F * fhPtPrimMC[7];                         //! Number of generated photon vs pT
 313   TH2F * fhPhiPrimMC[7];                        //! Phi of generted photon
 314   TH2F * fhYPrimMC[7];                          //! Rapidity of generated photon
 315
 316   TH1F * fhEPrimMCAcc[7];                       //! Number of generated photon vs energy, in calorimeter acceptance
 317   TH1F * fhPtPrimMCAcc[7];                      //! Number of generated photon vs pT, in calorimeter acceptance
 318   TH2F * fhPhiPrimMCAcc[7];                     //! Phi of generted photon, in calorimeter acceptance
 319   TH2F * fhYPrimMCAcc[7];                       //! Rapidity of generated photon, in calorimeter acceptance
 320
 321   // Shower Shape MC
 322
 323   TH2F * fhMCELambda0[6] ;                      //! E vs Lambda0     from MC particle
 324   TH2F * fhMCELambda1[6] ;                      //! E vs Lambda1     from MC particle
 325   TH2F * fhMCEDispersion[6] ;                   //! E vs Dispersion  from MC particle
 326
 327   TH2F * fhMCPhotonELambda0NoOverlap ;          //! E vs Lambda0     from MC photons, no overlap
 328   TH2F * fhMCPhotonELambda0TwoOverlap ;         //! E vs Lambda0     from MC photons, 2 particles overlap
 329   TH2F * fhMCPhotonELambda0NOverlap ;           //! E vs Lambda0     from MC photons, N particles overlap
 330
 331   TH2F * fhMCLambda0vsClusterMaxCellDiffE0[6];  //! Lambda0 vs fraction of energy of max cell for E < 2 GeV
 332   TH2F * fhMCLambda0vsClusterMaxCellDiffE2[6];  //! Lambda0 vs fraction of energy of max cell for 2< E < 6 GeV
 333   TH2F * fhMCLambda0vsClusterMaxCellDiffE6[6];  //! Lambda0 vs fraction of energy of max cell for E > 6 GeV
 334   TH2F * fhMCNCellsvsClusterMaxCellDiffE0[6];   //! NCells  vs fraction of energy of max cell for E < 2
 335   TH2F * fhMCNCellsvsClusterMaxCellDiffE2[6];   //! NCells  vs fraction of energy of max cell for 2 < E < 6 GeV
 336   TH2F * fhMCNCellsvsClusterMaxCellDiffE6[6];   //! NCells  vs fraction of energy of max cell for E > 6
 337   TH2F * fhMCNCellsE[6];                        //! NCells per cluster vs energy
 338   TH2F * fhMCMaxCellDiffClusterE[6];            //! Fraction of energy carried by cell with maximum energy
 339
 340   TH2F * fhMCEDispEta[6] ;                      //! shower dispersion in eta direction
 341   TH2F * fhMCEDispPhi[6] ;                      //! shower dispersion in phi direction
 342   TH2F * fhMCESumEtaPhi[6] ;                    //! shower dispersion in eta vs phi direction
 343   TH2F * fhMCEDispEtaPhiDiff[6] ;               //! shower dispersion in eta -phi direction
 344   TH2F * fhMCESphericity[6] ;                   //! shower sphericity, eta vs phi
 345   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]
 346   TH2F * fhMCLambda0DispEta[7][6] ;             //! shower shape correlation l0 vs disp eta
 347   TH2F * fhMCLambda0DispPhi[7][6] ;             //! shower shape correlation l0 vs disp phi
 348
 349   //Embedding
 350   TH2F * fhEmbeddedSignalFractionEnergy ;       //! Fraction of photon energy of embedded signal vs cluster energy
 351
 352   TH2F * fhEmbedPhotonELambda0FullSignal ;      //!  Lambda0 vs E for embedded photons with more than 90% of the cluster energy
 353   TH2F * fhEmbedPhotonELambda0MostlySignal ;    //!  Lambda0 vs E for embedded photons with 90%<fraction<50%
 354   TH2F * fhEmbedPhotonELambda0MostlyBkg ;       //!  Lambda0 vs E for embedded photons with 50%<fraction<10%
 355   TH2F * fhEmbedPhotonELambda0FullBkg ;         //!  Lambda0 vs E for embedded photons with less than 10% of the cluster energy
 356
 357   TH2F * fhEmbedPi0ELambda0FullSignal ;         //!  Lambda0 vs E for embedded photons with more than 90% of the cluster energy
 358   TH2F * fhEmbedPi0ELambda0MostlySignal ;       //!  Lambda0 vs E for embedded photons with 90%<fraction<50%
 359   TH2F * fhEmbedPi0ELambda0MostlyBkg ;          //!  Lambda0 vs E for embedded photons with 50%<fraction<10%
 360   TH2F * fhEmbedPi0ELambda0FullBkg ;            //!  Lambda0 vs E for embedded photons with less than 10% of the cluster energy
 361
 362   // Track Matching
 363   TH2F * fhTrackMatchedDEta[2]           ;      //! Eta distance between track and cluster vs cluster E, after and before photon cuts
 364   TH2F * fhTrackMatchedDPhi[2]           ;      //! Phi distance between track and cluster vs cluster E, after and before photon cuts
 365   TH2F * fhTrackMatchedDEtaDPhi[2]       ;      //! Eta vs Phi distance between track and cluster, E cluster > 0.5 GeV, after and before
 366
 367   TH2F * fhTrackMatchedDEtaPos[2]        ;      //! Eta distance between track and cluster vs cluster E, after and before photon cuts
 368   TH2F * fhTrackMatchedDPhiPos[2]        ;      //! Phi distance between track and cluster vs cluster E, after and before photon cuts
 369   TH2F * fhTrackMatchedDEtaDPhiPos[2]    ;      //! Eta vs Phi distance between track and cluster, E cluster > 0.5 GeV, after and before
 370
 371   TH2F * fhTrackMatchedDEtaNeg[2]        ;      //! Eta distance between track and cluster vs cluster E, after and before photon cuts
 372   TH2F * fhTrackMatchedDPhiNeg[2]        ;      //! Phi distance between track and cluster vs cluster E, after and before photon cuts
 373   TH2F * fhTrackMatchedDEtaDPhiNeg[2]    ;      //! Eta vs Phi distance between track and cluster, E cluster > 0.5 GeV, after and before photon cuts
 374
 375   TH2F * fhTrackMatchedDEtaTRD[2]        ;      //! Eta distance between track and cluster vs cluster E, after and before photon cuts, behind TRD
 376   TH2F * fhTrackMatchedDPhiTRD[2]        ;      //! Phi distance between track and cluster vs cluster E, after and before photon cuts, behind TRD
 377
 378   TH2F * fhTrackMatchedDEtaMCOverlap[2]  ;      //! Eta distance between track and cluster vs cluster E, several particle overlap, after and before photon cuts
 379   TH2F * fhTrackMatchedDPhiMCOverlap[2]  ;      //! Phi distance between track and cluster vs cluster E, several particle overlap, after and before photon cuts
 380   TH2F * fhTrackMatchedDEtaMCNoOverlap[2];      //! Eta distance between track and cluster vs cluster E, not other particle overlap, after and before photon cuts
 381   TH2F * fhTrackMatchedDPhiMCNoOverlap[2];      //! Phi distance between track and cluster vs cluster E, not other particle overlap, after and before photon cuts
 382   TH2F * fhTrackMatchedDEtaMCConversion[2];     //! Eta distance between track and cluster vs cluster E, originated in conversion, after and before photon cuts
 383   TH2F * fhTrackMatchedDPhiMCConversion[2];     //! Phi distance between track and cluster vs cluster E, originated in conversion, after and before photon cuts
 384
 385   TH2F * fhTrackMatchedMCParticle[2];           //! Trace origin of matched particle
 386   TH2F * fhdEdx[2];                             //! matched track dEdx vs cluster E, after and before photon cuts
 387   TH2F * fhEOverP[2];                           //! matched track E cluster over P track vs cluster E, after dEdx cut, after and before photon cuts
 388   TH2F * fhEOverPTRD[2];                        //! matched track E cluster over P track vs cluster E, after dEdx cut, after and before photon cuts, behind TRD
 389
 390   // Pile-up
 391   TH1F * fhPtPileUp[7];                         //! pT distribution of clusters before any selection
 392   TH1F * fhPtChargedPileUp[7];                  //! pT distribution of track matched clusters
 393   TH1F * fhPtPhotonPileUp[7];                   //! pT distribution of selected photons
 394   TH2F * fhLambda0PileUp[7];                    //! E vs M02 distribution of clusters, before any selection
 395   TH2F * fhLambda0ChargedPileUp[7];             //! E vs M02 distribution of clusters, track matched clusters
 396   TH2F * fhClusterTimeDiffPileUp[7];            //! E vs Time difference inside cluster, before any selection
 397   TH2F * fhClusterTimeDiffChargedPileUp[7];     //! E vs Time difference inside cluster for track matched clusters
 398   TH2F * fhClusterTimeDiffPhotonPileUp[7];      //! E vs Time difference inside cluster for selected photons
 399   TH2F * fhClusterEFracLongTimePileUp[7];       //! E vs fraction of cluster energy from cells with large time
 400   TH2F * fhTimePtNoCut;                         //! time of cluster vs Pt, no cut
 401   TH2F * fhTimePtSPD;                           //! time of cluster vs Pt, IsSPDPileUp
 402   TH2F * fhTimePtPhotonNoCut;                   //! time of photon cluster vs Pt, no cut
 403   TH2F * fhTimePtPhotonSPD;                     //! time of photon cluster vs Pt, IsSPDPileUp
 404   TH2F * fhTimeNPileUpVertSPD;                  //! time of cluster vs n pile-up vertices from SPD
 405   TH2F * fhTimeNPileUpVertTrack;                //! time of cluster vs n pile-up vertices from Tracks
 406   TH2F * fhTimeNPileUpVertContributors;         //! time of cluster vs n pile-up vertex from SPD contributors
 407   TH2F * fhTimePileUpMainVertexZDistance;       //! time of cluster vs difference of z main vertex and pile-up vertex
 408   TH2F * fhTimePileUpMainVertexZDiamond;        //! time of cluster vs difference of z diamond and pile-up vertex
 409   TH2F * fhClusterMultSPDPileUp[4];             //! E max cluster vs event cluster multiplicity, for tmax-tdiff cuts, pile up event
 410   TH2F * fhClusterMultNoPileUp[4];              //! E max cluster vs event cluster multiplicity, for tmax-tdiff cuts, not pile up event
 411   TH2F * fhEtaPhiBC0;                           //! eta/phi of clusters in BC=0
 412   TH2F * fhEtaPhiBCPlus;                        //! eta/phi of clusters in BC>0
 413   TH2F * fhEtaPhiBCMinus;                       //! eta/phi of clusters in BC<0
 414   TH2F * fhEtaPhiBC0PileUpSPD;                  //! eta/phi of clusters in BC=0, SPD pile-up
 415   TH2F * fhEtaPhiBCPlusPileUpSPD;               //! eta/phi of clusters in BC>0, SPD pile-up
 416   TH2F * fhEtaPhiBCMinusPileUpSPD;              //! eta/phi of clusters in BC<0, SPD pile-up
 417
 418   TH2F * fhPtNPileUpSPDVtx;                         //! cluster pt vs number of spd pile-up vertices
 419   TH2F * fhPtNPileUpTrkVtx;                     //! cluster pt vs number of track pile-up vertices
 420   TH2F * fhPtNPileUpSPDVtxTimeCut;                  //! cluster pt vs number of spd pile-up vertices, time cut +-25 ns
 421   TH2F * fhPtNPileUpTrkVtxTimeCut;              //! cluster pt vs number of track pile-up vertices, time cut +- 25 ns
 422   TH2F * fhPtNPileUpSPDVtxTimeCut2;                 //! cluster pt vs number of spd pile-up vertices, time cut +-75 ns
 423   TH2F * fhPtNPileUpTrkVtxTimeCut2;             //! cluster pt vs number of track pile-up vertices, time cut +- 75 ns
 424
 425   TH2F * fhPtPhotonNPileUpSPDVtx;                     //! photon pt vs number of spd pile-up vertices
 426   TH2F * fhPtPhotonNPileUpTrkVtx;               //! photon pt vs number of track pile-up vertices
 427   TH2F * fhPtPhotonNPileUpSPDVtxTimeCut;        //! photon pt vs number of spd pile-up vertices, time cut +-25 ns
 428   TH2F * fhPtPhotonNPileUpTrkVtxTimeCut;        //! photon pt vs number of track pile-up vertices, time cut +- 25 ns
 429   TH2F * fhPtPhotonNPileUpSPDVtxTimeCut2;       //! photon pt vs number of spd pile-up vertices, time cut +-75 ns
 430   TH2F * fhPtPhotonNPileUpTrkVtxTimeCut2;       //! photon pt vs number of track pile-up vertices, time cut +- 75 ns
 431
 432   AliAnaPhoton(              const AliAnaPhoton & g) ; // cpy ctor
 433   AliAnaPhoton & operator = (const AliAnaPhoton & g) ; // cpy assignment
 434
 435   ClassDef(AliAnaPhoton,34)
 436
 437 } ;
 438
 439 #endif//ALIANAPHOTON_H
 440
 441
 442