fPmax(0), // zero at 0.5 GeV/c for pp
fTime0(0)
{
+ // Reset T0 info
+ ResetT0info();
+ SetMomBoundary();
}
//_________________________________________________________________________
AliTOFPIDResponse::AliTOFPIDResponse(Double_t *param):
//
//fPmax=TMath::Exp(-0.5*3*3)/fSigma; // ~3 sigma at 0.5 GeV/c for PbPb
+
+ // Reset T0 info
+ ResetT0info();
+ SetMomBoundary();
}
//_________________________________________________________________________
Double_t
// If the operation is not possible, return a negative value.
//
- Double_t dpp = 0.01; //mean relative pt resolution;
- if (mom>0.5) dpp = 0.01*mom;
+ Double_t dpp=0.018*mass/mom; //mean relative pt resolution;
Double_t sigma = dpp*time/(1.+ mom*mom/(mass*mass));
+
+ Int_t index = GetMomBin(mom);
+
+ Double_t t0res = fT0resolution[index];
+
+ return TMath::Sqrt(sigma*sigma + 50.0*50.0/mom/mom + fSigma*fSigma + t0res*t0res);
+
+}
+//_________________________________________________________________________
+Int_t AliTOFPIDResponse::GetMomBin(Float_t p) const{
+ //
+ // Returns the momentum bin index
+ //
- return TMath::Sqrt(sigma*sigma + fSigma*fSigma);
+ Int_t i=0;
+ while(p > fPCutMin[i] && i < fNmomBins) i++;
+ if(i > 0) i--;
+
+ return i;
}
+//_________________________________________________________________________
+void AliTOFPIDResponse::SetMomBoundary(){
+ //
+ // Set boundaries for momentum bins
+ //
+ fPCutMin[0] = 0.3;
+ fPCutMin[1] = 0.5;
+ fPCutMin[2] = 0.6;
+ fPCutMin[3] = 0.7;
+ fPCutMin[4] = 0.8;
+ fPCutMin[5] = 0.9;
+ fPCutMin[6] = 1;
+ fPCutMin[7] = 1.2;
+ fPCutMin[8] = 1.5;
+ fPCutMin[9] = 2;
+ fPCutMin[10] = 3;
+}
+//_________________________________________________________________________
+Float_t AliTOFPIDResponse::GetStartTime(Float_t mom) const {
+ //
+ // Returns event_time value as estimated by TOF combinatorial algorithm
+ //
+
+ Int_t ibin = GetMomBin(mom);
+ return GetT0bin(ibin);
+
+}
+//_________________________________________________________________________
+Float_t AliTOFPIDResponse::GetStartTimeRes(Float_t mom) const {
+ //
+ // Returns event_time resolution as estimated by TOF combinatorial algorithm
+ //
+
+ Int_t ibin = GetMomBin(mom);
+ return GetT0binRes(ibin);
+
+}
+//_________________________________________________________________________
+Int_t AliTOFPIDResponse::GetStartTimeMask(Float_t mom) const {
+ //
+ // Returns event_time mask
+ //
+
+ Int_t ibin = GetMomBin(mom);
+ return GetT0binMask(ibin);
+
+}