Simplify AliAnaPhoton, move conversion utilities to other analysis class
[u/mrichter/AliRoot.git] / PWG4 / PartCorrDep / AliAnaPhoton.h
index caee5f47bcb92603ee4006cfeecd0f7e1e4c266e..369b7f6ca6d33e11d3fb44139c82882889ed0d1b 100755 (executable)
 // --- ROOT system ---
 class TH2F ;
 class TH1F;
-class TH3D;
 class TString ;
 class TObjString;
+class TList ;
 
 // --- ANALYSIS system ---
 #include "AliAnaPartCorrBaseClass.h"
-class AliStack;
-class TParticle;
-
-class TList ;
 
 class AliAnaPhoton : public AliAnaPartCorrBaseClass {
 
  public: 
-  AliAnaPhoton() ; // default ctor
-  virtual ~AliAnaPhoton() ; //virtual dtor
+  AliAnaPhoton() ;              // default ctor
+  virtual ~AliAnaPhoton() { ; } // virtual dtor
  private:
-  AliAnaPhoton(const AliAnaPhoton & g) ; // cpy ctor
-  AliAnaPhoton & operator = (const AliAnaPhoton & g) ;//cpy assignment
+  AliAnaPhoton(const AliAnaPhoton & g) ;               // cpy ctor
+  AliAnaPhoton & operator = (const AliAnaPhoton & g) ; // cpy assignment
 
  public:
        
@@ -65,16 +61,13 @@ class AliAnaPhoton : public AliAnaPartCorrBaseClass {
   
   void         FillAcceptanceHistograms();
 
-  //           Fill Shower Shape histograms
   void         FillShowerShapeHistograms( AliVCluster* cluster, const Int_t mcTag) ;
   
   void         SwitchOnFillShowerShapeHistograms()    { fFillSSHistograms = kTRUE  ; }
   void         SwitchOffFillShowerShapeHistograms()   { fFillSSHistograms = kFALSE ; }  
   
   
-  //---------------------------------------
   // Analysis parameters setters getters
-  //---------------------------------------
   
   TString      GetCalorimeter()                 const { return fCalorimeter        ; }
   void         SetCalorimeter(TString  & det)         { fCalorimeter = det         ; }
@@ -95,37 +88,7 @@ class AliAnaPhoton : public AliAnaPartCorrBaseClass {
   Bool_t       IsTrackMatchRejectionOn()        const { return fRejectTrackMatch   ; }
   void         SwitchOnTrackMatchRejection()          { fRejectTrackMatch = kTRUE  ; }
   void         SwitchOffTrackMatchRejection()         { fRejectTrackMatch = kFALSE ; }  
-  
-  // ** Conversion pair analysis **
-  
-  Float_t      GetMassCut()                     const { return fMassCut            ; }
-  void         SetMassCut(Float_t m)                  { fMassCut    = m            ; }
-  
-  Bool_t       IsCheckConversionOn()            const { return fCheckConversion    ; }
-  void         SwitchOnConversionChecker()            { fCheckConversion = kTRUE   ; }
-  void         SwitchOffConversionChecker()           { fCheckConversion = kFALSE  ; }  
-       
-  Bool_t       AreConvertedPairsInAOD()         const { return fAddConvertedPairsToAOD   ; }
-  void         SwitchOnAdditionConvertedPairsToAOD()  { fAddConvertedPairsToAOD = kTRUE  ; 
-                                                        fCheckConversion        = kTRUE  ; }
-  void         SwitchOffAdditionConvertedPairsToAOD() { fAddConvertedPairsToAOD = kFALSE ; }  
-       
-  Bool_t       AreConvertedPairsRemoved()       const { return fRemoveConvertedPair      ; }
-  void         SwitchOnConvertedPairsRemoval()        { fRemoveConvertedPair  = kTRUE    ; 
-                                                        fCheckConversion      = kTRUE    ; }
-  void         SwitchOffConvertedPairsRemoval()       { fRemoveConvertedPair  = kFALSE   ; }    
-  
-  void         SetConvAsymCut(Float_t c)              { fConvAsymCut = c           ; }
-  Float_t      GetConvAsymCut()                 const { return fConvAsymCut        ; }
-  
-  void         SetConvDEtaCut(Float_t c)              { fConvDEtaCut = c           ; }
-  Float_t      GetConvDEtaCut()                 const { return fConvDEtaCut        ; }
-  
-  void         SetConvDPhiCut(Float_t min, Float_t max)  { fConvDPhiMinCut = min   ;  
-                                                           fConvDPhiMaxCut = max   ; }
-  Float_t      GetConvDPhiMinCut()              const { return fConvDPhiMinCut     ; }
-  Float_t      GetConvDPhiMaxCut()              const { return fConvDPhiMaxCut     ; }
-  
+         
   void         FillNOriginHistograms(Int_t n)         { fNOriginHistograms = n ; 
     if(n > 14) fNOriginHistograms = 14; }
   void         FillNPrimaryHistograms(Int_t n)        { fNPrimaryHistograms= n ;
@@ -158,16 +121,6 @@ class AliAnaPhoton : public AliAnaPartCorrBaseClass {
   Int_t    fNOriginHistograms;           // Fill only NOriginHistograms of the 14 defined types
   Int_t    fNPrimaryHistograms;          // Fill only NPrimaryHistograms of the 7 defined types
 
-  //Conversion pairs selection cuts
-  Bool_t   fCheckConversion;             // Combine pairs of clusters with mass close to 0
-  Bool_t   fRemoveConvertedPair;         // Remove conversion pairs
-  Bool_t   fAddConvertedPairsToAOD;      // Put Converted pairs in AOD
-  Float_t  fMassCut;                     // Mass cut for the conversion pairs selection  
-  Float_t  fConvAsymCut;                 // Select conversion pairs when asymmetry is smaller than cut
-  Float_t  fConvDEtaCut;                 // Select conversion pairs when deta of pair smaller than cut
-  Float_t  fConvDPhiMinCut;              // Select conversion pairs when dphi of pair lager than cut
-  Float_t  fConvDPhiMaxCut;              // Select conversion pairs when dphi of pair smaller than cut
-
   //Histograms 
   TH2F * fhNCellsE;                      //! number of cells in cluster vs E 
   TH2F * fhMaxCellDiffClusterE;          //! Fraction of energy carried by cell with maximum energy
@@ -179,28 +132,6 @@ class AliAnaPhoton : public AliAnaPartCorrBaseClass {
   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
   
-  //Conversion pairs
-  TH1F * fhPtPhotonConv   ;              //! Number of identified photon vs transerse momentum 
-  TH2F * fhEtaPhiPhotonConv  ;           //! Pseudorapidity vs Phi of identified  photon for transerse momentum > 0.5, for converted
-  TH2F * fhEtaPhi05PhotonConv  ;         //! Pseudorapidity vs Phi of identified  photon for transerse momentum < 0.5, for converted
-  TH2F * fhConvDeltaEta;                 //! Small mass photons, correlation in eta
-  TH2F * fhConvDeltaPhi;                 //! Small mass photons, correlation in phi
-  TH2F * fhConvDeltaEtaPhi;              //! Small mass photons, correlation in phi and eta
-  TH2F * fhConvAsym;                     //! Small mass photons, correlation in energy asymmetry
-  TH2F * fhConvPt;                       //! Small mass photons, pT of pair
-  
-  //Vertex distance
-  TH2F * fhConvDistEta;                   //! Approx distance to vertex vs cluster Eta 
-  TH2F * fhConvDistEn;                    //! Approx distance to vertex vs Energy
-  TH2F * fhConvDistMass;                  //! Approx distance to vertex vs Mass
-  TH2F * fhConvDistEtaCutEta;             //! Approx distance to vertex vs cluster Eta, dEta < 0.05 
-  TH2F * fhConvDistEnCutEta;              //! Approx distance to vertex vs Energy, dEta < 0.05
-  TH2F * fhConvDistMassCutEta;            //! Approx distance to vertex vs Mass, dEta < 0.05
-  TH2F * fhConvDistEtaCutMass;            //! Approx distance to vertex vs cluster Eta, dEta < 0.05, m < 10 MeV 
-  TH2F * fhConvDistEnCutMass;             //! Approx distance to vertex vs Energy, dEta < 0.05, m < 10 MeV
-  TH2F * fhConvDistEtaCutAsy;             //! Approx distance to vertex vs cluster Eta, dEta < 0.05, m < 10 MeV, A < 0.1
-  TH2F * fhConvDistEnCutAsy;              //! Approx distance to vertex vs energy, dEta < 0.05, m < 10 MeV, A < 0.1
-
   //Shower shape
   
   TH2F * fhDispE;                         //! cluster dispersion vs E
@@ -231,78 +162,35 @@ class AliAnaPhoton : public AliAnaPartCorrBaseClass {
     
   //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
+  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 * 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 
   
-  //Conversion pairs analysis histograms
-  TH1F * fhPtConversionTagged;                //! Number of identified gamma from Conversion , tagged as conversion 
-  TH1F * fhPtAntiNeutronTagged;               //! Number of identified gamma from AntiNeutrons gamma, tagged as conversion 
-  TH1F * fhPtAntiProtonTagged;                //! Number of identified gamma from AntiProtons gamma, tagged as conversion 
-  TH1F * fhPtUnknownTagged;                   //! Number of identified gamma from unknown, tagged as conversion 
+  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   
   
-  TH2F * fhEtaPhiConversion  ;                //! Pseudorapidity vs Phi for transerse momentum > 0.5, for MC converted
-  TH2F * fhEtaPhi05Conversion  ;              //! Pseudorapidity vs Phi for transerse momentum < 0.5, for MC converted
-  
-  TH2F * fhConvDeltaEtaMCConversion;          //! Small mass cluster pairs, correlation in eta, origin of both clusters is conversion
-  TH2F * fhConvDeltaPhiMCConversion;          //! Small mass cluster pairs, correlation in phi, origin of both clusters is conversion
-  TH2F * fhConvDeltaEtaPhiMCConversion;       //! Small mass cluster pairs, correlation in eta-phi, origin of both clusters is conversion
-  TH2F * fhConvAsymMCConversion;              //! Small mass cluster pairs, correlation in energy asymmetry, origin of both clusters is conversion
-  TH2F * fhConvPtMCConversion;                //! Small mass cluster pairs, pt of pair, origin of both clusters is conversion
-  TH2F * fhConvDispersionMCConversion;        //! Small mass cluster pairs, dispersion of cluster 1 vs cluster 2
-  TH2F * fhConvM02MCConversion;               //! Small mass cluster pairs, m02 of cluster 1 vs cluster 2 
-
-  TH2F * fhConvDeltaEtaMCAntiNeutron;         //! Small mass cluster pairs, correlation in eta, origin of both clusters is anti neutron
-  TH2F * fhConvDeltaPhiMCAntiNeutron;         //! Small mass cluster pairs, correlation in phi, origin of both clusters is anti neutron
-  TH2F * fhConvDeltaEtaPhiMCAntiNeutron;      //! Small mass cluster pairs, correlation in eta-phi, origin of both clusters is anti neutron
-  TH2F * fhConvAsymMCAntiNeutron;             //! Small mass cluster pairs, correlation in energy asymmetry, origin of both clusters is anti neutron
-  TH2F * fhConvPtMCAntiNeutron;               //! Small mass cluster pairs, pt of pair, origin of both clusters is anti neutron
-  TH2F * fhConvDispersionMCAntiNeutron;       //! Small mass cluster pairs, dispersion of cluster 1 vs cluster 2, origin of both clusters is anti neutron
-  TH2F * fhConvM02MCAntiNeutron;              //! Small mass cluster pairs, m02 of cluster 1 vs cluster 2, origin of both clusters is anti neutron
-
-  TH2F * fhConvDeltaEtaMCAntiProton;          //! Small mass cluster pairs, correlation in eta, origin of both clusters is anti proton
-  TH2F * fhConvDeltaPhiMCAntiProton;          //! Small mass cluster pairs, correlation in phi, origin of both clusters is anti proton
-  TH2F * fhConvDeltaEtaPhiMCAntiProton;       //! Small mass cluster pairs, correlation in eta-phi, origin of both clusters is anti proton
-  TH2F * fhConvAsymMCAntiProton;              //! Small mass cluster pairs, correlation in energy asymmetry, origin of both clusters is anti proton
-  TH2F * fhConvPtMCAntiProton;                //! Small mass cluster pairs, pt of pairs, origin of both clusters is anti proton
-  TH2F * fhConvDispersionMCAntiProton;        //! Small mass cluster pairs, dispersion of cluster 1 vs cluster 2, origin of both clusters is anti proton
-  TH2F * fhConvM02MCAntiProton;               //! Small mass cluster pairs, m02 of cluster 1 vs cluster 2, origin of both clusters is anti proton
-
-  TH2F * fhConvDeltaEtaMCString;              //! Small mass cluster pairs, correlation in eta, origin of both clusters is string
-  TH2F * fhConvDeltaPhiMCString;              //! Small mass cluster pairs, correlation in phi, origin of both clusters is string
-  TH2F * fhConvDeltaEtaPhiMCString;           //! Small mass cluster pairs, correlation in eta-phi, origin of both clusters is string
-  TH2F * fhConvAsymMCString;                  //! Small mass cluster pairs, correlation in energy asymmetry, origin of both clusters is string
-  TH2F * fhConvPtMCString;                    //! Small mass cluster pairs, pt of pairs, origin of both clusters is string
-  TH2F * fhConvDispersionMCString;            //! Small mass cluster pairs, dispersion of cluster 1 vs cluster 2, origin of both clusters is string
-  TH2F * fhConvM02MCString;                   //! Small mass cluster pairs, m02 of cluster 1 vs cluster 2, origin of both clusters is string
-  TH2F * fhConvDistMCConversion;              //! Calculated conversion distance vs real distance to vertex       
-  TH2F * fhConvDistMCConversionCuts;          //! Calculated conversion distance vs real distance to vertex        
-
   // 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 * 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 * 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
@@ -314,23 +202,22 @@ class AliAnaPhoton : public AliAnaPartCorrBaseClass {
   TH2F * fhMCMaxCellDiffClusterE[6];            //! Fraction of energy carried by cell with maximum energy
 
   //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
   
-   ClassDef(AliAnaPhoton,17)
+   ClassDef(AliAnaPhoton,18)
 
 } ;
  
-
 #endif//ALIANAPHOTON_H