]> git.uio.no Git - u/mrichter/AliRoot.git/blame - STEER/STEERBase/AliHMPIDPIDResponse.h
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
* Add full support for signal delta (signal-expected, signal/expected)
[u/mrichter/AliRoot.git] / STEER / STEERBase / AliHMPIDPIDResponse.h
CommitLineData
567624b5 1#ifndef ALIHMPIDPIDRESPONSE_H
2#define ALIHMPIDPIDRESPONSE_H
3
4/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
5 * See cxx source for full Copyright notice */
6
7//***********************************************************
8// Class AliHMPIDPIDResponse
9//
10// HMPID class to perfom particle identification
11//
12// Author: G. Volpe, giacomo.volpe@cern.ch
13//***********************************************************
14
15
16#include <TNamed.h> //base class
17#include <TVector3.h>
18#include <TVector2.h>
19
20#include "AliPID.h"
21
22class AliVTrack;
23class TObjArray;
24class TGeoHMatrix;
25
26class AliHMPIDPIDResponse : public TNamed
27{
28public :
29 AliHMPIDPIDResponse(); //ctor
30 AliHMPIDPIDResponse(const AliHMPIDPIDResponse& c); //copy constructor
31 AliHMPIDPIDResponse &operator=(const AliHMPIDPIDResponse& c); //dummy assignment operator
32 virtual ~AliHMPIDPIDResponse() {;} //dtor
33
34 enum EChamberData{kMinCh=0,kMaxCh=6,kMinPc=0,kMaxPc=5}; //Segmenation
35 enum EPadxData{kPadPcX=80,kMinPx=0,kMaxPx=79,kMaxPcx=159}; //Segmentation structure along x
36 enum EPadyData{kPadPcY=48,kMinPy=0,kMaxPy=47,kMaxPcy=143}; //Segmentation structure along y
37
38 Double_t GetExpectedSignal (const AliVTrack *vTrk, AliPID::EParticleType specie ) const;
39 Double_t GetExpectedSigma (const AliVTrack *vTrk, AliPID::EParticleType specie ) const; //Find the sigma for a given ThetaCerTh
40 Double_t GetNumberOfSigmas (const AliVTrack *vTrk, AliPID::EParticleType specie ) const; //Find the expected Cherenkov angle for a given track
41 void GetProbability (const AliVTrack *vTrk, Int_t nSpecies,Double_t *prob ) const; //Find the PID probability array
1d59271b 42 Double_t GetSignalDelta (const AliVTrack *vTrk, AliPID::EParticleType specie, Bool_t ratio=kFALSE) const;
567624b5 43 void Propagate (const TVector3 dir, TVector3 &pos, Double_t z ) const; //propagate photon alogn the line
44 void Refract (TVector3 &dir, Double_t n1, Double_t n2 ) const; //refract photon on the boundary
45 TVector2 TracePhot (Double_t xRa, Double_t yRa, Double_t thRa, Double_t phRa, Double_t ckovThe,Double_t ckovPhi) const; //trace photon created by track to PC
46 void Trs2Lors (Double_t thRa, Double_t phRa, TVector3 dirCkov, Double_t &thetaCer,Double_t &phiCer) const; //TRS to LORS
47 TVector2 TraceForward (Double_t xRa, Double_t yRa, TVector3 dirCkov ) const; //tracing forward a photon from (x,y) to PC
48 void SetTrack (Double_t xRad, Double_t yRad, Double_t theta,Double_t phi ) {fTrkDir.SetMagThetaPhi(1,theta,phi); fTrkPos.Set(xRad,yRad);} //set track parameter at RAD
49 Double_t RadThick ( ) const {return 1.5;} //Radiator thickness
50 Double_t WinThick ( ) const {return 0.5;} //Window thickness
51 Double_t GapThick ( ) const {return 8.0;} //Proximity gap thicknes
52 Double_t GetRefIdx ( ) const {return fRefIdx;} //running refractive index
53 Double_t WinIdx ( ) const {return 1.5787;} //Mean refractive index of WIN material (SiO2)
54 Double_t GapIdx ( ) const {return 1.0005;} //Mean refractive index of GAP material (CH4)
55 static Bool_t IsInside (Float_t x,Float_t y,Float_t d=0 ) {return x>-d&&y>-d&&x<fgkMaxPcX[kMaxPc]+d&&y<fgkMaxPcY[kMaxPc]+d; } //is point inside chamber boundaries?
56 static Bool_t IsInDead (Float_t x,Float_t y ); //is the point in a dead area?
57 static Float_t SizeAllX ( ) {return fgAllX;} //all PCs size x, [cm]
58 static Float_t SizeAllY ( ) {return fgAllY;} //all PCs size y, [cm]
59 static void IdealPosition(Int_t iCh,TGeoHMatrix *m ); //ideal position of given chamber
60
61 Double_t SigLoc (Double_t trkTheta,Double_t trkPhi,Double_t ckovTh,Double_t ckovPh,Double_t beta) const; //error due to cathode segmetation
62 Double_t SigGeom (Double_t trkTheta,Double_t trkPhi,Double_t ckovTh,Double_t ckovPh,Double_t beta) const; //error due to unknown photon origin
63 Double_t SigCrom (Double_t trkTheta,Double_t ckovTh,Double_t ckovPh,Double_t beta ) const; //error due to unknonw photon energy
64 Double_t Sigma2 (Double_t trkTheta,Double_t trkPhi,Double_t ckovTh,Double_t ckovPh ) const; //photon candidate sigma^2
65 Double_t GetNMean (const AliVTrack *vTrk ) const;
66 static Double_t SigmaCorrFact(Int_t iPart, Double_t occupancy ) ; //correction factor for theoretical resolution
67
68 void SetRefIndexArray (TObjArray *array ) {fRefIndexArray = array;}
69 TObjArray* GetRefIndexArray ( ) const {return fRefIndexArray;}
70
71//
72private:
73
74 Double_t ExpectedSignal (const AliVTrack *vTrk, Double_t nmean, AliPID::EParticleType specie ) const;
75 Double_t ExpectedSigma (const AliVTrack *vTrk, Double_t nmean, AliPID::EParticleType specie ) const; //Find the sigma for a given ThetaCerTh
76
77protected:
78
79 static /*const*/ Float_t fgkMinPcX[6]; //limits PC
80 static /*const*/ Float_t fgkMinPcY[6]; //limits PC
81 static /*const*/ Float_t fgkMaxPcX[6]; //limits PC
82 static /*const*/ Float_t fgkMaxPcY[6];
83
84 static Float_t fgCellX, fgCellY, fgPcX, fgPcY, fgAllX, fgAllY; //definition of HMPID geometric parameters
85
86 TGeoHMatrix *fM[7]; //pointers to matrices defining HMPID chambers rotations-translations
87
88 Double_t fRefIdx; //running refractive index of C6F14
89 TVector3 fTrkDir; //track direction in LORS at RAD
90 TVector2 fTrkPos; //track positon in LORS at RAD
91 TObjArray *fRefIndexArray; //array of refracive index funxtion;
92
93 ClassDef(AliHMPIDPIDResponse,1)
94};
95#endif // #ifdef AliHMPIDPIDResponse_cxx
96