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]);}
+ void Lors2MarsVec(Int_t c,Double_t *m,Double_t *l )const{fM[c]->LocalToMasterVect(m,l); }//LRS->MRS
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
void Point (Int_t c,Double_t *p,Int_t plane )const{Lors2Mars(c,0,0,p,plane);} //point of given chamber plane
// return gradT*y+t1;
Double_t halfPadSize = 0.5*SizePadY();
Double_t gradT = (TMath::Log(SizePcY()) - TMath::Log(halfPadSize))/(TMath::Log(tHigh)-TMath::Log(tLow));
- return tLow*TMath::Power(y/halfPadSize,1./gradT);
+ if(y<0) y = 0;
+ return tLow + TMath::Power(y/halfPadSize,1./gradT);
}
//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
#endif