]>
Commit | Line | Data |
---|---|---|
10d100d4 | 1 | #ifndef ALIITSPIDRESPONSE_H |
2 | #define ALIITSPIDRESPONSE_H | |
3 | /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * | |
4 | * See cxx source for full Copyright notice */ | |
5 | ||
6 | /* $Id$ */ | |
7 | ||
8 | //------------------------------------------------------- | |
9 | // ITS PID response class | |
10 | // | |
11 | // | |
12 | //------------------------------------------------------- | |
13 | //#include <Rtypes.h> | |
14 | #include <TObject.h> | |
10d100d4 | 15 | #include "AliPID.h" |
16 | ||
17 | class AliITSPIDResponse : public TObject { | |
18 | ||
19 | public: | |
15e979c9 | 20 | AliITSPIDResponse(Bool_t isMC=kFALSE); |
56576f1e | 21 | //AliITSPIDResponse(Double_t *param); |
10d100d4 | 22 | ~AliITSPIDResponse() {} |
99daa709 | 23 | |
24 | void SetBetheBlochParamsITSTPC(Double_t* param){ | |
25 | for(Int_t iPar=0; iPar<5; iPar++) fBBtpcits[iPar]=param[iPar]; | |
26 | } | |
27 | void SetBetheBlochParamsITSsa(Double_t* param){ | |
28 | for(Int_t iPar=0; iPar<5; iPar++) fBBsa[iPar]=param[iPar]; | |
29 | } | |
88f46717 | 30 | void SetElectronBetheBlochParamsITSsa(Double_t* param){ |
31 | for(Int_t iPar=0; iPar<5; iPar++) fBBsaElectron[iPar]=param[iPar]; | |
32 | } | |
99daa709 | 33 | |
15e979c9 | 34 | Double_t BetheAleph(Double_t p,Double_t mass) const; |
35 | Double_t Bethe(Double_t p, Double_t mass, Bool_t iSA=kFALSE) const; | |
36 | Double_t GetResolution(Double_t bethe, Int_t nPtsForPid=4, Bool_t isSA=kFALSE) const; | |
b52bfc67 | 37 | void GetITSProbabilities(Float_t mom, Double_t qclu[4], Double_t condprobfun[AliPID::kSPECIES],Bool_t isMC=kFALSE) const; |
15e979c9 | 38 | Float_t GetNumberOfSigmas(Float_t mom, Float_t signal, AliPID::EParticleType type, Int_t nPtsForPid=4, Bool_t isSA=kFALSE) const { |
39 | Float_t bethe = Bethe(mom,AliPID::ParticleMass(type),isSA); | |
40 | return (signal - bethe)/GetResolution(bethe,nPtsForPid,isSA); | |
41 | } | |
8abeb05b | 42 | Int_t GetParticleIdFromdEdxVsP(Float_t mom, Float_t signal, Bool_t isSA=kFALSE) const; |
10d100d4 | 43 | |
44 | private: | |
45 | ||
15e979c9 | 46 | |
10d100d4 | 47 | // Data members for truncated mean method |
15e979c9 | 48 | Float_t fRes; // relative dEdx resolution |
10d100d4 | 49 | Double_t fKp1; // ALEPH BB param 1 |
50 | Double_t fKp2; // ALEPH BB param 2 | |
51 | Double_t fKp3; // ALEPH BB param 3 | |
52 | Double_t fKp4; // ALEPH BB param 4 | |
15e979c9 | 53 | Double_t fKp5; // ALEPH BB param |
54 | Double_t fBBsa[5]; // parameters of BB for SA tracks | |
88f46717 | 55 | Double_t fBBsaElectron[5]; // parameters of BB for SA tracks |
15e979c9 | 56 | Double_t fBBtpcits[5]; // parameters of BB for TPC+ITS tracks |
57 | Float_t fResolSA[5]; // resolutions vs. n. of SDD/SSD points | |
58 | Float_t fResolTPCITS[5]; // resolutions vs. n. of SDD/SSD points | |
10d100d4 | 59 | |
88f46717 | 60 | ClassDef(AliITSPIDResponse,3) // ITS PID class |
10d100d4 | 61 | }; |
62 | ||
63 | #endif | |
64 | ||
65 |