void Print(Option_t *opt="") const; //print current parametrization
static inline AliHMPIDParam* Instance(); //pointer to AliHMPIDParam singleton
- Double_t MeanIdxRad () {return 1.29204;} //<--TEMPORAR--> to be removed in future Mean ref index C6F14
- Double_t MeanIdxWin () {return 1.57819;} //<--TEMPORAR--> to be removed in future. Mean ref index quartz
- Float_t DistCut () {return 1.0;} //<--TEMPORAR--> to be removed in future. Cut for MIP-TRACK residual
- Float_t QCut () {return 100;} //<--TEMPORAR--> to be removed in future. Separation PHOTON-MIP charge
- Float_t MultCut () {return 200;} //<--TEMPORAR--> to be removed in future. Multiplicity cut to activate WEIGHT procedure
+ Double_t MeanIdxRad ()const {return 1.29204;} //<--TEMPORAR--> to be removed in future Mean ref index C6F14
+ Double_t MeanIdxWin ()const {return 1.57819;} //<--TEMPORAR--> to be removed in future. Mean ref index quartz
+ Float_t DistCut ()const {return 1.0;} //<--TEMPORAR--> to be removed in future. Cut for MIP-TRACK residual
+ Float_t QCut ()const {return 100;} //<--TEMPORAR--> to be removed in future. Separation PHOTON-MIP charge
+ Float_t MultCut ()const {return 200;} //<--TEMPORAR--> to be removed in future. Multiplicity cut to activate WEIGHT procedure
static Int_t Stack(Int_t evt=-1,Int_t tid=-1); //Print stack info for event and tid
static Int_t StackCount(Int_t pid,Int_t evt); //Counts stack particles of given sort in given event
static void IdealPosition(Int_t iCh,TGeoHMatrix *m); //ideal position of given chamber
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); }
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
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
- 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(pt/l[2]); ph=TMath::ATan(l[1]/l[0]);}
+ 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(pt/l[2]);
+ ph=TMath::ATan2(l[1],l[0]);}
TVector3 Norm (Int_t c )const{Double_t n[3]; Norm(c,n); return TVector3(n); }//norm
void Norm (Int_t c,Double_t *n )const{Double_t l[3]={0,0,1};fM[c]->LocalToMasterVect(l,n); }//norm
protected:
AliHMPIDParam(); //default ctor is protected to enforce it to be singleton
static AliHMPIDParam *fgInstance; //static pointer to instance of AliHMPIDParam singleton
- TGeoHMatrix *fM[7]; //poiners to matrices defining HMPID chambers rotations-translations
- Float_t fX; //x shift of LORS with respect to rotated MARS
- Float_t fY; //y shift of LORS with respect to rotated MARS
+ TGeoHMatrix *fM[7]; //pointers to matrices defining HMPID chambers rotations-translations
+ Float_t fX; //x shift of LORS with respect to rotated MARS
+ Float_t fY; //y shift of LORS with respect to rotated MARS
ClassDef(AliHMPIDParam,0) //HMPID main parameters class
};