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
void SwitchOnFillOriginHisto() { fFillOriginHisto = kTRUE ; }
void SwitchOffFillOriginHisto() { fFillOriginHisto = kFALSE ; }
+
+ void SwitchOnFillArmenterosThetaStarHisto() { fFillArmenterosThetaStar = kTRUE ; }
+ void SwitchOffFillArmenterosThetaStarHisto() { fFillArmenterosThetaStar = kFALSE ; }
//MC analysis related methods
void SwitchOnMultipleCutAnalysisInSimulation() { fMultiCutAnaSim = kTRUE ; }
void SwitchOffMultipleCutAnalysisInSimulation() { fMultiCutAnaSim = kFALSE ; }
+ void SwitchOnCheckAcceptanceInSector() { fCheckAccInSector = kTRUE ; }
+ void SwitchOffCheckAcceptanceInSector(){ fCheckAccInSector = kFALSE ; }
+
void FillAcceptanceHistograms();
- void FillMCVersusRecDataHistograms(const Int_t index1, const Int_t index2,
- const Float_t pt1, const Float_t pt2,
- const Int_t ncells1, const Int_t ncells2,
- const Double_t mass, const Double_t pt, const Double_t asym,
- const Double_t deta, const Double_t dphi);
+ void FillMCVersusRecDataHistograms(Int_t index1, Int_t index2,
+ Float_t pt1, Float_t pt2,
+ Int_t ncells1, Int_t ncells2,
+ Double_t mass, Double_t pt, Double_t asym,
+ Double_t deta, Double_t dphi);
+
+ 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
Bool_t fFillAngleHisto; // Fill histograms with pair opening angle
Bool_t fFillAsymmetryHisto; // Fill histograms with asymmetry vs pt
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
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
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
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 ...
+
AliAnaPi0( const AliAnaPi0 & api0) ; // cpy ctor
AliAnaPi0 & operator = (const AliAnaPi0 & api0) ; // cpy assignment
- ClassDef(AliAnaPi0,24)
+ ClassDef(AliAnaPi0,29)
} ;