]>
Commit | Line | Data |
---|---|---|
d3da6dc4 | 1 | #ifndef AliHMPIDParam_h |
2 | #define AliHMPIDParam_h | |
3010c308 | 3 | /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * |
4 | * See cxx source for full Copyright notice */ | |
d3da6dc4 | 5 | |
3010c308 | 6 | /* $Id$ */ |
7 | ||
8 | #include <TMath.h> | |
d3da6dc4 | 9 | #include <TNamed.h> //base class |
10 | #include <TGeoManager.h> //Instance() | |
11 | #include <TVector3.h> //Lors2Mars() Mars2Lors() | |
12 | ||
d3da6dc4 | 13 | // Class providing all the needed parametrised information |
14 | // to construct the geometry, to define segmentation and to provide response model | |
15 | // In future will also provide all the staff needed for alignment and calibration | |
16 | ||
17 | class AliHMPIDParam :public TNamed | |
18 | { | |
19 | public: | |
20 | //ctor&dtor | |
21 | virtual ~AliHMPIDParam() {for(Int_t i=0;i<7;i++) delete fM[i]; delete fgInstance; fgInstance=0;} | |
22 | void Print(Option_t *opt="") const; //print current parametrization | |
23 | static inline AliHMPIDParam* Instance(); //pointer to AliHMPIDParam singleton | |
24 | ||
25 | Double_t MeanIdxRad () {return 1.29204;}//??????????? | |
26 | Double_t MeanIdxWin () {return 1.57819;}//??????????? | |
27 | static Int_t Stack(Int_t evt=-1,Int_t tid=-1); //Print stack info for event and tid | |
28 | static Int_t StackCount(Int_t pid,Int_t evt); //Counts stack particles of given sort in given event | |
1d4857c5 | 29 | static void IdealPosition(Int_t iCh,TGeoHMatrix *m); //ideal position of given chamber |
30 | //trasformation methodes | |
d3da6dc4 | 31 | void Lors2Mars (Int_t c,Float_t x,Float_t y,Double_t *m,Int_t pl=kPc)const{Double_t z=0; switch(pl){case kPc:z=8.0;break; case kAnod:z=7.806;break; case kRad:z=-1.25; break;} Double_t l[3]={x-fX,y-fY,z}; fM[c]->LocalToMaster(l,m); } |
32 | TVector3 Lors2Mars (Int_t c,Float_t x,Float_t y, Int_t pl=kPc)const{Double_t m[3];Lors2Mars(c,x,y,m,pl); return TVector3(m); }//MRS->LRS | |
33 | void Mars2Lors (Int_t c,Double_t *m,Float_t &x,Float_t &y )const{Double_t l[3];fM[c]->MasterToLocal(m,l);x=l[0]+fX;y=l[1]+fY;}//MRS->LRS | |
32675400 | 34 | void Mars2LorsVec(Int_t c,Double_t *m,Float_t &th,Float_t &ph )const{Double_t l[3]; fM[c]->MasterToLocalVect(m,l); Float_t pt=TMath::Sqrt(l[0]*l[0]+l[1]*l[1]); th=TMath::ATan(l[2]/pt); ph=TMath::ATan(l[0]/pt);} |
d3da6dc4 | 35 | TVector3 Norm (Int_t c )const{Double_t n[3]; Norm(c,n); return TVector3(n); }//norm |
36 | void Norm (Int_t c,Double_t *n )const{Double_t l[3]={0,0,1};fM[c]->LocalToMasterVect(l,n); }//norm | |
f3bae3e2 | 37 | |
d3da6dc4 | 38 | enum EPlaneId {kPc,kRad,kAnod}; //3 planes in chamber |
39 | protected: | |
40 | AliHMPIDParam(); //default ctor is protected to enforce it to be singleton | |
41 | static AliHMPIDParam *fgInstance; //static pointer to instance of AliHMPIDParam singleton | |
42 | TGeoHMatrix *fM[7]; //poiners to matrices defining HMPID chambers rotations-translations | |
43 | Float_t fX; //x shift of LORS with respect to rotated MARS | |
44 | Float_t fY; //y shift of LORS with respect to rotated MARS | |
45 | ClassDef(AliHMPIDParam,0) //HMPID main parameters class | |
46 | }; | |
cf7e313e | 47 | |
d3da6dc4 | 48 | //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
49 | AliHMPIDParam* AliHMPIDParam::Instance() | |
50 | { | |
51 | // Return pointer to the AliHMPIDParam singleton. | |
52 | // Arguments: none | |
53 | // Returns: pointer to the instance of AliHMPIDParam or 0 if no geometry | |
1d4857c5 | 54 | if(!fgInstance) new AliHMPIDParam; |
d3da6dc4 | 55 | return fgInstance; |
56 | }//Instance() | |
57 | //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ | |
58 | ||
1d4857c5 | 59 | |
d3da6dc4 | 60 | #endif |