Number of sigma pedestal cut increased to 4
[u/mrichter/AliRoot.git] / STEER / AliEMCALPIDResponse.h
1 #ifndef AliEMCALPIDResponse_h\r
2 #define AliEMCALPIDResponse_h\r
3 \r
4 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *\r
5  * See cxx source for full Copyright notice                               */\r
6 \r
7 //////////////////////////////////////////////////////////////////////////\r
8 //                                                                      //\r
9 // AliEMCALPIDResponse                                                  //\r
10 //                                                                      //\r
11 // EMCAL class to perfom PID                                            //\r
12 // This is a prototype and still under development                      //\r
13 //                                                                      //\r
14 //                                                                      //\r
15 //////////////////////////////////////////////////////////////////////////\r
16 \r
17 #include "AliPID.h"\r
18 \r
19 class AliEMCALPIDResponse \r
20 {\r
21 public : \r
22              AliEMCALPIDResponse();    //ctor\r
23     virtual ~AliEMCALPIDResponse() {;} //dtor\r
24 \r
25 \r
26     // Getters\r
27     Int_t     GetPtBin(Float_t pt) const;\r
28     Int_t     GetNPtBins() const {return fNptBins;};\r
29 \r
30     Float_t   GetLowEoP() const {return fLowEoP;};\r
31     Float_t   GetHighEoP() const {return fHighEoP;};\r
32 \r
33     Double_t  GetExpectedSignal( Float_t pt, AliPID::EParticleType n,  Int_t charge) const;\r
34     Double_t  GetExpectedSigma ( Float_t pt, AliPID::EParticleType n,  Int_t charge) const;\r
35     Double_t  GetNumberOfSigmas( Float_t pt,  Float_t eop, AliPID::EParticleType n,  Int_t charge) const;\r
36     Double_t  GetExpectedNorm  ( Float_t pt, AliPID::EParticleType n,  Int_t charge) const;\r
37     Double_t  GetLowProb       ( Float_t pt, AliPID::EParticleType n,  Int_t charge) const;\r
38     Double_t  GetHighProb      ( Float_t pt, AliPID::EParticleType n,  Int_t charge) const;\r
39 \r
40     //Setters\r
41     void   SetPtBoundary();\r
42     void   SetParametrizations();\r
43 \r
44     // EMCAL probability -> should go to another place?\r
45     Double_t ComputeEMCALProbability( Float_t pt, Float_t eop, Int_t charge, Double_t *pEMCAL) const;\r
46 \r
47 protected:\r
48   \r
49 private:\r
50   AliEMCALPIDResponse( AliEMCALPIDResponse& r);                //dummy copy ructor\r
51   AliEMCALPIDResponse &operator=( AliEMCALPIDResponse& r);     //dummy assignment operator\r
52 \r
53   TF1 *fNorm;                            // Gauss function for normalizing NON electron probabilities \r
54 \r
55   static const Int_t fNptBins   = 6;     // number of momentum bins\r
56   static const Float_t fLowEoP  = 0.5;   // lower E/p threshold for NON electrons\r
57   static const Float_t fHighEoP = 1.5;   // upper E/p threshold for NON electrons\r
58 \r
59   Float_t fPtCutMin[fNptBins+1];                       // min values for pt bins\r
60   Float_t fMeanEoP[2*AliPID::kSPECIES][fNptBins];      // mean value of E/p distribution (charge dependent)\r
61   Float_t fSigmaEoP[2*AliPID::kSPECIES][fNptBins];     // mean value of E/p distribution (charge dependent)\r
62   Float_t fProbLow[2*AliPID::kSPECIES][fNptBins];      // probability below E/p threshold for NON electrons (charge dependent)\r
63   Float_t fProbHigh[2*AliPID::kSPECIES][fNptBins];     // probability above E/p threshold for NON electrons (charge dependent)\r
64 \r
65 \r
66   ClassDef(AliEMCALPIDResponse,0)\r
67 };\r
68 \r
69 #endif // #ifdef AliEMCALPIDResponse_cxx\r
70 \r