1 #ifndef ALINEUTRALMESONSELECTION_H
2 #define ALINEUTRALMESONSELECTION_H
3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * See cxx source for full Copyright notice */
5 /* $Id: AliNeutralMesonSelection.h 27413 2008-07-18 13:28:12Z gconesab $ */
7 //_________________________________________________________________________
8 // Class that contains methods to select candidate pairs to neutral meson
9 // 2 main selections, invariant mass around pi0 (also any other mass),
10 // apperture angle to distinguish from combinatorial.
11 // There is a 3rd cut based on the gamma correlation on phi or pt.
12 //-- Author: Gustavo Conesa (INFN-LNF)
14 // --- ROOT system ---
18 class TLorentzVector ;
22 //--- ANALYSIS system ---
24 class AliNeutralMesonSelection : public TObject {
27 AliNeutralMesonSelection() ; // default ctor
28 virtual ~AliNeutralMesonSelection() { ; } //virtual dtor
30 AliNeutralMesonSelection(const AliNeutralMesonSelection & g) ; // cpy ctor
31 AliNeutralMesonSelection & operator = (const AliNeutralMesonSelection & g) ;//cpy assignment
36 TList * GetCreateOutputObjects();
38 void InitParameters();
40 void Print(const Option_t * opt) const;
42 Bool_t AreNeutralMesonSelectionHistosKept() const { return fKeepNeutralMesonHistos ; }
43 void KeepNeutralMesonSelectionHistos(Bool_t keep) { fKeepNeutralMesonHistos = keep ; }
45 Bool_t SelectPair(TLorentzVector particlei, TLorentzVector particlej, TString calo) ;
47 void SetParticle(TString particleName) ; // Do some default settings for "Pi0" or "Eta"
49 // Asymmetry selection
51 Float_t GetAsymmetryCut() const { return fAsymmetryCut ; }
52 void SetAsymmetryCut(Float_t asy) { fAsymmetryCut = asy ; }
54 void SwitchOnAsymmetryCut() { fUseAsymmetryCut = kTRUE ; }
55 void SwitchOffAsymmetryCut() { fUseAsymmetryCut = kFALSE ; }
57 //Opening angle selection
59 Double_t GetAngleMaxParam(Int_t i) const { return fAngleMaxParam.At(i) ; }
60 void SetAngleMaxParam(Int_t i, Double_t par) { fAngleMaxParam.AddAt(par,i) ; }
62 void SetShiftMinAngleCut(Float_t a, Float_t b) { fShiftMinAngle[0] = a ;
63 fShiftMinAngle[1] = b ; }
65 void SwitchOnAngleSelection() { fUseAngleCut = kTRUE ; }
66 void SwitchOffAngleSelection() { fUseAngleCut = kFALSE ; }
68 Bool_t IsAngleInWindow(const Float_t angle, const Float_t e) const ;
70 //Invariant mass selection
72 Double_t GetInvMassMaxCut() const { return fInvMassMaxCut ; }
73 Double_t GetInvMassMinCut() const { return fInvMassMinCut ; }
75 void SetInvMassCutRange(Double_t invmassmin, Double_t invmassmax)
76 { fInvMassMaxCut =invmassmax; fInvMassMinCut =invmassmin ; }
78 void SetInvMassCutMaxParameters(Float_t a, Float_t b, Float_t c)
79 { fInvMassMaxCutParam[0] = a ;
80 fInvMassMaxCutParam[1] = b ;
81 fInvMassMaxCutParam[2] = c ; }
83 Double_t GetMass() const { return fM ; }
84 void SetMass(Double_t m) { fM = m ; }
86 //Histogrammes setters and getters
88 virtual void SetHistoERangeAndNBins(Float_t min, Float_t max, Int_t n) {
94 Int_t GetHistoNEBins() const { return fHistoNEBins ; }
95 Float_t GetHistoEMin() const { return fHistoEMin ; }
96 Float_t GetHistoEMax() const { return fHistoEMax ; }
98 virtual void SetHistoAngleRangeAndNBins(Float_t min, Float_t max, Int_t n) {
99 fHistoNAngleBins = n ;
100 fHistoAngleMax = max ;
101 fHistoAngleMin = min ;
104 Int_t GetHistoNAngleBins() const { return fHistoNAngleBins ; }
105 Float_t GetHistoAngleMin() const { return fHistoAngleMin ; }
106 Float_t GetHistoAngleMax() const { return fHistoAngleMax ; }
108 virtual void SetHistoIMRangeAndNBins(Float_t min, Float_t max, Int_t n) {
114 Int_t GetHistoNIMBins() const { return fHistoNIMBins ; }
115 Float_t GetHistoIMMin() const { return fHistoIMMin ; }
116 Float_t GetHistoIMMax() const { return fHistoIMMax ; }
121 Float_t fAsymmetryCut ; // Asymmetry cut
122 Bool_t fUseAsymmetryCut; // Use the asymmetry cut
124 Double_t fM ; // Mass of the neutral meson
125 Double_t fInvMassMaxCut ; // Invariant Mass cut maximum
126 Double_t fInvMassMinCut ; // Invariant Masscut minimun
127 Double_t fInvMassMaxCutParam[3]; // Variable invariant mass max cut, for pi0 in EMCAL
129 TArrayD fAngleMaxParam ; // Max opening angle selection parameters
130 Bool_t fUseAngleCut ; // Select pairs depending on their opening angle
131 Float_t fShiftMinAngle[2] ; // Correction shift for min angle from true kinematic limit, resolution effects
133 Bool_t fKeepNeutralMesonHistos ; // Keep neutral meson selection histograms
136 TH2F * fhAnglePairNoCut ; //! Aperture angle of decay photons, no cuts
137 TH2F * fhAnglePairOpeningAngleCut ; //! Aperture angle of decay photons, cut on opening angle
138 TH2F * fhAnglePairAsymmetryCut ; //! Aperture angle of decay photons, asymmetry cut
139 TH2F * fhAnglePairAllCut ; //! Aperture angle of decay photons, all cuts
141 TH2F * fhInvMassPairNoCut ; //! Invariant mass of decay photons, no cuts
142 TH2F * fhInvMassPairOpeningAngleCut ; //! Invariant mass of decay photons, cut on opening angle
143 TH2F * fhInvMassPairAsymmetryCut ; //! Invariant mass of decay photons, asymmetry cut
144 TH2F * fhInvMassPairAllCut ; //! Invariant mass of decay photons, all cuts
146 TH2F * fhAsymmetryNoCut ; //! Asymmetry of decay photons, no cuts
147 TH2F * fhAsymmetryOpeningAngleCut ; //! Asymmetry of decay photons, cut on opening angle
148 TH2F * fhAsymmetryAllCut ; //! Asymmetry of decay photons, all cuts
150 //Histograms binning and range
151 Int_t fHistoNEBins ; // Number of bins in pi0 E axis
152 Float_t fHistoEMax ; // Maximum value of pi0 E histogram range
153 Float_t fHistoEMin ; // Minimum value of pi0 E histogram range
155 Int_t fHistoNAngleBins ; // Number of bins in angle axis
156 Float_t fHistoAngleMax ; // Maximum value of angle histogram range
157 Float_t fHistoAngleMin ; // Minimum value of angle histogram range
159 Int_t fHistoNIMBins ; // Number of bins in Invariant Mass axis
160 Float_t fHistoIMMax ; // Maximum value of Invariant Mass histogram range
161 Float_t fHistoIMMin ; // Minimum value of Invariant Mass histogram range
163 ClassDef(AliNeutralMesonSelection,6)
167 #endif //ALINEUTRALMESONSELECTION_H