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 */
7 /* History of cvs commits:
10 * Revision 1.2 2007/08/17 12:40:04 schutz
11 * New analysis classes by Gustavo Conesa
13 * Revision 1.1.2.1 2007/07/26 10:32:09 schutz
14 * new analysis classes in the the new analysis framework
18 //_________________________________________________________________________
19 // Class that contains methods to select candidate pairs to neutral meson
20 // 2 main selections, invariant mass around pi0 (also any other mass),
21 // apperture angle to distinguish from combinatorial.
22 // There is a 3rd cut based on the gamma correlation on phi or pt.
23 //-- Author: Gustavo Conesa (INFN-LNF)
25 // --- ROOT system ---
29 class TLorentzVector ;
34 //--- ANALYSIS system ---
37 class AliNeutralMesonSelection : public TObject {
41 AliNeutralMesonSelection() ; // default ctor
42 AliNeutralMesonSelection(const AliNeutralMesonSelection & g) ; // cpy ctor
43 AliNeutralMesonSelection & operator = (const AliNeutralMesonSelection & g) ;//cpy assignment
44 virtual ~AliNeutralMesonSelection() ; //virtual dtor
46 TList * GetCreateOutputObjects();
48 Double_t GetAngleMaxParam(Int_t i) const {return fAngleMaxParam.At(i) ; }
49 void SetAngleMaxParam(Int_t i, Double_t par){fAngleMaxParam.AddAt(par,i) ; }
51 Double_t GetInvMassMaxCut() const {return fInvMassMaxCut ; }
52 Double_t GetInvMassMinCut() const {return fInvMassMinCut ; }
53 void SetInvMassCutRange(Double_t invmassmin, Double_t invmassmax)
54 {fInvMassMaxCut =invmassmax; fInvMassMinCut =invmassmin;}
56 Double_t GetMass() const {return fM ; }
57 void SetMass(Double_t m) { fM =m ; }
59 Bool_t AreNeutralMesonSelectionHistosKept() const { return fKeepNeutralMesonHistos ; }
60 void KeepNeutralMesonSelectionHistos(Bool_t keep) { fKeepNeutralMesonHistos = keep ; }
62 void InitParameters();
63 Bool_t IsAngleInWindow(const Float_t angle, const Float_t e) const ;
64 void Print(const Option_t * opt) const;
66 Bool_t SelectPair(TLorentzVector particlei, TLorentzVector particlej) ;
69 Double_t fM ; //mass of the neutral meson
70 Double_t fInvMassMaxCut ; // Invariant Mass cut maximum
71 Double_t fInvMassMinCut ; // Invariant Masscut minimun
72 TArrayD fAngleMaxParam ; //Max opening angle selection parameters
73 Double_t fMinPt; // Minimum pt
74 Bool_t fKeepNeutralMesonHistos ; // Keep neutral meson selection histograms
77 TH2F * fhAnglePairNoCut ; //Aperture angle of decay photons, no cuts
78 TH2F * fhAnglePairOpeningAngleCut ; //Aperture angle of decay photons, cut on opening angle
79 TH2F * fhAnglePairAllCut ; //Aperture angle of decay photons, all cuts
80 TH2F * fhInvMassPairNoCut ; //Invariant mass of decay photons, no cuts
81 TH2F * fhInvMassPairOpeningAngleCut ; //Invariant mass of decay photons, cut on opening angle
82 TH2F * fhInvMassPairAllCut ; //Invariant mass of decay photons, all cuts
84 ClassDef(AliNeutralMesonSelection,1)
89 #endif //ALINEUTRALMESONSELECTION_H