Merge branch 'master' of https://git.cern.ch/reps/AliRoot

[u/mrichter/AliRoot.git] / PWGGA / CaloTrackCorrelations / AliAnaPi0.h

diff --git a/PWGGA/CaloTrackCorrelations/AliAnaPi0.h b/PWGGA/CaloTrackCorrelations/AliAnaPi0.h
index 1d06afd9d80d64f7967c5d657933cce3e3ae1b70..9e17c54dce021a629201451efe1a3d57f8e6f521 100755 (executable)
--- a/PWGGA/CaloTrackCorrelations/AliAnaPi0.h
+++ b/PWGGA/CaloTrackCorrelations/AliAnaPi0.h
@@ -44,11 +44,6 @@ class AliAnaPi0 : public AliAnaCaloTrackCorrBaseClass {
   void         MakeAnalysisFillHistograms();
   
   void         InitParameters();
-
-  //Calorimeter options
-  TString      GetCalorimeter()         const   { return fCalorimeter           ; }
-  void         SetCalorimeter(TString & det)    { fCalorimeter         = det    ; }
-  void         SetNumberOfModules(Int_t nmod)   { fNModules            = nmod   ; }
   
   //-------------------------------
   // EVENT Bin Methods
@@ -125,6 +120,9 @@ class AliAnaPi0 : public AliAnaCaloTrackCorrBaseClass {
   void         SwitchOnMultipleCutAnalysisInSimulation()  { fMultiCutAnaSim = kTRUE  ; }
   void         SwitchOffMultipleCutAnalysisInSimulation() { fMultiCutAnaSim = kFALSE ; }
   
+  void         SwitchOnCheckAcceptanceInSector() { fCheckAccInSector   = kTRUE  ; }
+  void         SwitchOffCheckAcceptanceInSector(){ fCheckAccInSector   = kFALSE ; }
+  
   void         FillAcceptanceHistograms();
   void         FillMCVersusRecDataHistograms(Int_t    index1,  Int_t    index2,
                                              Float_t  pt1,     Float_t  pt2,
@@ -132,15 +130,13 @@ class AliAnaPi0 : public AliAnaCaloTrackCorrBaseClass {
                                              Double_t mass,    Double_t pt,     Double_t asym,
                                              Double_t deta,    Double_t dphi);
   
-  void         FillArmenterosThetaStar(Int_t pdg,             TLorentzVector meson,
-                                       TLorentzVector daugh1, TLorentzVector daugh2);
+  void         FillArmenterosThetaStar(Int_t pdg);
 
   
   private:
 
   TList ** fEventsList ;               //![GetNCentrBin()*GetNZvertBin()*GetNRPBin()] Containers for photons in stored events
 
-  TString  fCalorimeter ;              // Select Calorimeter for IM
   Int_t    fNModules ;                 // Number of EMCAL/PHOS modules, set as many histogras as modules 
   
   Bool_t   fUseAngleCut ;              // Select pairs depending on their opening angle
@@ -172,6 +168,14 @@ class AliAnaPi0 : public AliAnaCaloTrackCorrBaseClass {
   Bool_t   fFillOriginHisto;           // Fill histograms depending on their origin
   Bool_t   fFillArmenterosThetaStar;   // Fill armenteros histograms
   
+  Bool_t   fCheckAccInSector;          // Check that the decay pi0 falls in the same SM or sector
+  
+  TLorentzVector fPhotonMom1;          //! photon cluster momentum
+  TLorentzVector fPhotonMom1Boost;     //! photon cluster momentum
+  TLorentzVector fPhotonMom2;          //! photon cluster momentum
+  TLorentzVector fPi0Mom;              //! pi0 cluster momentum
+  TVector3       fProdVertex;          //! production vertex
+  
   //Histograms
   
   //Event characterization
@@ -248,6 +252,9 @@ class AliAnaPi0 : public AliAnaCaloTrackCorrBaseClass {
   TH1F *   fhPrimPi0AccPt ;            //! Spectrum of primary with accepted daughters
   TH2F *   fhPrimPi0Y ;                //! Rapidity distribution of primary particles  vs pT
   TH2F *   fhPrimPi0AccY ;             //! Rapidity distribution of primary with accepted daughters  vs pT
+  TH2F *   fhPrimPi0Yeta ;             //! PseudoRapidity distribution of primary particles  vs pT
+  TH2F *   fhPrimPi0YetaYcut ;         //! PseudoRapidity distribution of primary particles  vs pT, Y<1
+  TH2F *   fhPrimPi0AccYeta ;          //! PseudoRapidity distribution of primary with accepted daughters  vs pT
   TH2F *   fhPrimPi0Phi ;              //! Azimutal distribution of primary particles  vs pT
   TH2F *   fhPrimPi0AccPhi;            //! Azimutal distribution of primary with accepted daughters  vs pT
   TH2F *   fhPrimPi0OpeningAngle ;     //! Opening angle of pair versus pair energy, primaries
@@ -265,6 +272,9 @@ class AliAnaPi0 : public AliAnaCaloTrackCorrBaseClass {
   TH1F *   fhPrimEtaAccPt ;            //! Spectrum of primary with accepted daughters
   TH2F *   fhPrimEtaY ;                //! Rapidity distribution of primary particles vs pT
   TH2F *   fhPrimEtaAccY ;             //! Rapidity distribution of primary with accepted daughters  vs pT
+  TH2F *   fhPrimEtaYeta ;             //! PseudoRapidity distribution of primary particles vs pT
+  TH2F *   fhPrimEtaYetaYcut ;         //! PseudoRapidity distribution of primary particles vs pT, Y<1
+  TH2F *   fhPrimEtaAccYeta ;          //! PseudoRapidity distribution of primary with accepted daughters  vs pT
   TH2F *   fhPrimEtaPhi ;              //! Azimutal distribution of primary particles  vs pT
   TH2F *   fhPrimEtaAccPhi;            //! Azimutal distribution of primary with accepted daughters     vs pT
   TH2F *   fhPrimEtaOpeningAngle ;     //! Opening angle of pair versus pair energy, eta primaries
@@ -298,19 +308,24 @@ class AliAnaPi0 : public AliAnaCaloTrackCorrBaseClass {
   TH2F *   fhMCPi0PtOrigin ;           //! Mass of reoconstructed pi0 pairs  in calorimeter vs mother
   TH2F *   fhMCEtaPtOrigin ;           //! Mass of reoconstructed pi0 pairs  in calorimeter vs mother
 
+  TH2F *   fhMCPi0ProdVertex;          //! Spectrum of selected pi0 vs production vertex
+  TH2F *   fhMCEtaProdVertex;          //! Spectrum of selected eta vs production vertex
+  TH2F *   fhPrimPi0ProdVertex;        //! Spectrum of primary pi0 vs production vertex
+  TH2F *   fhPrimEtaProdVertex;        //! Spectrum of primary eta vs production vertex
+  
   TH2F *   fhReMCFromConversion ;      //! Invariant mass of 2 clusters originated in conversions
   TH2F *   fhReMCFromNotConversion ;   //! Invariant mass of 2 clusters not originated in conversions
   TH2F *   fhReMCFromMixConversion ;   //! Invariant mass of 2 clusters one from conversion and the other not
 
-  TH2F *    fhArmPrimPi0[4];           //! Armenteros plots for primary pi0 in 6 energy bins
-  TH2F *    fhArmPrimEta[4];           //! Armenteros plots for primary eta in 6 energy bins
-  TH2F *    fhCosThStarPrimPi0;        //! cos(theta*) plots vs E for primary pi0, same as asymmetry ...
-  TH2F *    fhCosThStarPrimEta;        //! cos(theta*) plots vs E for primary eta, same as asymmetry ...
+  TH2F *   fhArmPrimPi0[4];            //! Armenteros plots for primary pi0 in 6 energy bins
+  TH2F *   fhArmPrimEta[4];            //! Armenteros plots for primary eta in 6 energy bins
+  TH2F *   fhCosThStarPrimPi0;         //! cos(theta*) plots vs E for primary pi0, same as asymmetry ...
+  TH2F *   fhCosThStarPrimEta;         //! cos(theta*) plots vs E for primary eta, same as asymmetry ...
   
   AliAnaPi0(              const AliAnaPi0 & api0) ; // cpy ctor
   AliAnaPi0 & operator = (const AliAnaPi0 & api0) ; // cpy assignment
   
-  ClassDef(AliAnaPi0,26)
+  ClassDef(AliAnaPi0,29)
 } ;