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 //_________________________________________________________________________
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 ;
23 //--- ANALYSIS system ---
26 class AliNeutralMesonSelection : public TObject {
30 AliNeutralMesonSelection() ; // default ctor
31 AliNeutralMesonSelection(const AliNeutralMesonSelection & g) ; // cpy ctor
32 AliNeutralMesonSelection & operator = (const AliNeutralMesonSelection & g) ;//cpy assignment
33 virtual ~AliNeutralMesonSelection() ; //virtual dtor
35 TList * GetCreateOutputObjects();
37 Double_t GetAngleMaxParam(Int_t i) const {return fAngleMaxParam.At(i) ; }
38 void SetAngleMaxParam(Int_t i, Double_t par){fAngleMaxParam.AddAt(par,i) ; }
40 Double_t GetInvMassMaxCut() const {return fInvMassMaxCut ; }
41 Double_t GetInvMassMinCut() const {return fInvMassMinCut ; }
42 void SetInvMassCutRange(Double_t invmassmin, Double_t invmassmax)
43 {fInvMassMaxCut =invmassmax; fInvMassMinCut =invmassmin;}
45 Double_t GetMass() const {return fM ; }
46 void SetMass(Double_t m) { fM =m ; }
48 Bool_t AreNeutralMesonSelectionHistosKept() const { return fKeepNeutralMesonHistos ; }
49 void KeepNeutralMesonSelectionHistos(Bool_t keep) { fKeepNeutralMesonHistos = keep ; }
51 void InitParameters();
52 Bool_t IsAngleInWindow(const Float_t angle, const Float_t e) const ;
53 void Print(const Option_t * opt) const;
55 Bool_t SelectPair(TLorentzVector particlei, TLorentzVector particlej) ;
58 Double_t fM ; //mass of the neutral meson
59 Double_t fInvMassMaxCut ; // Invariant Mass cut maximum
60 Double_t fInvMassMinCut ; // Invariant Masscut minimun
61 TArrayD fAngleMaxParam ; //Max opening angle selection parameters
63 Bool_t fKeepNeutralMesonHistos ; // Keep neutral meson selection histograms
66 TH2F * fhAnglePairNoCut ; //Aperture angle of decay photons, no cuts
67 TH2F * fhAnglePairOpeningAngleCut ; //Aperture angle of decay photons, cut on opening angle
68 TH2F * fhAnglePairAllCut ; //Aperture angle of decay photons, all cuts
69 TH2F * fhInvMassPairNoCut ; //Invariant mass of decay photons, no cuts
70 TH2F * fhInvMassPairOpeningAngleCut ; //Invariant mass of decay photons, cut on opening angle
71 TH2F * fhInvMassPairAllCut ; //Invariant mass of decay photons, all cuts
73 ClassDef(AliNeutralMesonSelection,1)
78 #endif //ALINEUTRALMESONSELECTION_H