1 #ifndef ALIDIELECTRONSIGNALBASE_H
2 #define ALIDIELECTRONSIGNALBASE_H
4 /* Copyright(c) 1998-2009, ALICE Experiment at CERN, All rights reserved. *
5 * See cxx source for full Copyright notice */
7 //#############################################################
9 //# Class AliDielectronSignalBase #
10 //# Manage Cuts on the legs of the pair #
13 //# Anton Andronic, GSI / A.Andronic@gsi.de #
14 //# Ionut C. Arsene, GSI / I.C.Arsene@gsi.de #
15 //# Julian Book, Uni Ffm / Julian.Book@cern.ch #
16 //# Frederick Kramer, Uni Ffm, / Frederick.Kramer@cern.ch #
17 //# Magnus Mager, CERN / Magnus.Mager@cern.ch #
18 //# WooJin J. Park, GSI / W.J.Park@gsi.de #
19 //# Jens Wiechula, Uni HD / Jens.Wiechula@cern.ch #
21 //#############################################################
32 class AliDielectronSignalBase : public TNamed {
34 enum EBackgroundMethod {
40 AliDielectronSignalBase();
41 AliDielectronSignalBase(const char*name, const char* title);
43 virtual ~AliDielectronSignalBase();
45 void SetIntegralRange(Double_t min, Double_t max) {fIntMin=min;fIntMax=max;}
46 void SetFitRange(Double_t min, Double_t max) {fFitMin=min; fFitMax=max;}
47 void SetRebin(Int_t factor) {fRebin = factor;}
48 void SetMethod(EBackgroundMethod method) {fMethod = method;}
50 const TVectorD& GetValues() const {return fValues;}
51 const TVectorD& GetErrors() const {return fErrors;}
53 Double_t GetIntegralMin() const { return fIntMin; }
54 Double_t GetIntegralMax() const { return fIntMax; }
55 Double_t GetSignal() const { return fValues(0);}
56 Double_t GetSignalError() const { return fErrors(0);}
57 Double_t GetBackground() const { return fValues(1);}
58 Double_t GetBackgroundError() const { return fErrors(1);}
59 Double_t GetSignificance() const { return fValues(2);}
60 Double_t GetSignificanceError() const { return fErrors(2);}
61 Double_t GetSB() const { return fValues(3);}
62 Double_t GetSBError() const { return fErrors(3);}
63 Double_t GetMass() const { return fValues(4);}
64 Double_t GetMassError() const { return fErrors(4);}
65 Double_t GetMassWidth() const { return fValues(5);}
66 Double_t GetMassWidthError() const { return fErrors(5);}
68 TH1F* GetSignalHistogram() {return fHistSignal;}
69 TH1F* GetBackgroundHistogram() {return fHistBackground;}
71 virtual void Print(Option_t *option="") const;
74 This function needs to be implemented by the signal extraction classes.
75 Here all the work should be done.
77 The signal extraction is done on the mass spectra.
78 The TObjArray should contain the Inv. Mass spectra of the 10 possible combinations
79 for single and mixed events defined in AliDielectron.cxx
81 virtual void Process(TObjArray * const /*arrhist*/) = 0;
85 TH1F *fHistSignal; // histogram of pure signal
86 TH1F *fHistBackground; // histogram of background (fitted=0, like-sign=1, event mixing=2)
87 TH1F *fHistDataPM; // histogram of selected +- pair candidates
88 TH1F *fHistDataPP; // histogram of selected ++ pair candidates
89 TH1F *fHistDataMM; // histogram of selected -- pair candidates
91 TVectorD fValues; // values
92 TVectorD fErrors; // value errors
94 Double_t fIntMin; // signal extraction range min
95 Double_t fIntMax; // signal extraction range max
96 Double_t fFitMin; // fit range lowest inv. mass
97 Double_t fFitMax; // fit range highest inv. mass
99 Int_t fRebin; // histogram rebin factor
100 EBackgroundMethod fMethod; // method for background substraction
102 Bool_t fProcessed; // flag
104 void SetSignificanceAndSOB(); // calculate the significance and S/B
105 TPaveText* DrawStats(Double_t x1=0., Double_t y1=0., Double_t x2=0., Double_t y2=0.);
107 AliDielectronSignalBase(const AliDielectronSignalBase &c);
108 AliDielectronSignalBase &operator=(const AliDielectronSignalBase &c);
110 ClassDef(AliDielectronSignalBase,3) // Dielectron SignalBase
113 inline void AliDielectronSignalBase::SetSignificanceAndSOB()
116 // Calculate S/B and significance
119 fValues(3) = (fValues(1)>0 ? fValues(0)/fValues(1) : 0);
120 Float_t epsSig = (fValues(0)>0 ? fErrors(0)/fValues(0) : 0);
121 Float_t epsBknd = (fValues(1)>0 ? fErrors(1)/fValues(1) : 0);
122 fErrors(3) = fValues(3)*TMath::Sqrt(epsSig*epsSig + epsBknd*epsBknd);
124 fValues(2) = ((fValues(0)+fValues(1))>0 ? fValues(0)/TMath::Sqrt(fValues(0)+fValues(1)) : 0);
125 Float_t s = fValues(0); Float_t b = fValues(1);
126 fErrors(2) = ((s+b)>0 ? TMath::Sqrt((s*(s+2*b)*(s+2*b)+b*s*s)/(4*TMath::Power(s+b,3))) : 0);