Bethe-Bloch formula normalized to one at the minimum (Yuri)
authorhristov <hristov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 18 Oct 2007 07:28:08 +0000 (07:28 +0000)
committerhristov <hristov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 18 Oct 2007 07:28:08 +0000 (07:28 +0000)
ITS/AliITSReconstructor.cxx
ITS/AliITSpidESD.cxx
ITS/AliITSpidESD1.cxx
ITS/AliITSpidESD1.h

index 13dcdb9ac5cebc6dd4ada5f8c7afcb9a567ed2b6..3d65b7cfba8d5afa7b4ccb5a4a35d3734dacf639 100644 (file)
@@ -159,10 +159,10 @@ AliTracker* AliITSReconstructor::CreateTracker() const
   }
   else{
     Info("FillESD","ITS default PID\n");
-    Double_t parITS[] = {0.15, 10.}; //PH positions of the MIP peak
+    Double_t parITS[] = {79.,0.13, 10.}; //IB: this is  "pp tuning"
     nc->fItsPID = new AliITSpidESD1(parITS);
   }
-
   return tracker;
   
 }
index 2270234b86a2ddc98c90149de5631801ae5b96cd..2d464c08aa9b8de332a4ac084c0516e869e92e06 100644 (file)
@@ -49,7 +49,7 @@ Double_t AliITSpidESD::Bethe(Double_t p,Double_t mass) {
   Double_t beta=mom/En;
   Double_t Tmax=2*me*beta*beta*gamma*gamma/(1+2*gamma*me/Mass+(me/Mass)*(me/Mass));
   Double_t deltaover2=28.816*1e-6*TMath::Sqrt(dens*ZovA)+TMath::Log(beta*gamma)-0.5;
-  Double_t FNor=80./105.48;
+  Double_t FNor=0.009164; //normalizing to 1 at the minimum of ionization
 
   return K*ZovA*1/(beta*beta)*(0.5*TMath::Log(2*me*beta*beta*gamma*gamma*Tmax/(I*I))-beta*beta-deltaover2)*2.33*1000*0.03*FNor;
 
index 9683ae3023b8752860bf2b94f25293783f3d467c..b57d8459de335ae7bc2c09c6f411c95b52d1c418 100755 (executable)
@@ -29,6 +29,7 @@
 ClassImp(AliITSpidESD1)
 
 AliITSpidESD1::AliITSpidESD1(): AliITSpidESD(),
+fMIP(0),
 fRes(0),
 fRange(0) 
 {
@@ -36,8 +37,9 @@ fRange(0)
 }
 //_________________________________________________________________________
 AliITSpidESD1::AliITSpidESD1(Double_t *param): AliITSpidESD(),
-fRes(param[0]),
-fRange(param[1])
+fMIP(param[0]),
+fRes(param[1]),
+fRange(param[2])
 {
   //
   //  The main constructor
@@ -57,7 +59,7 @@ Int_t AliITSpidESD1::MakePID(AliESDEvent *event)
     if ((t->GetStatus()&AliESDtrack::kITSin )==0)
       if ((t->GetStatus()&AliESDtrack::kITSout)==0) continue;
     Double_t mom=t->GetP();
-    Double_t dedx=t->GetITSsignal();
+    Double_t dedx=t->GetITSsignal()/fMIP;
     Int_t ns=AliPID::kSPECIES;
     Double_t p[10];
     for (Int_t j=0; j<ns; j++) {
index f18eb04533e1de7edd2df353500f29b4b8252ef6..2f5e032b0b15571a8ed00d9024c39d322a953820 100755 (executable)
@@ -25,9 +25,10 @@ public:
   {AliError("Method should not be used!"); return 0;}
 
 private:
+  Double_t fMIP;          // dE/dx for MIP (relative units)
   Double_t fRes;          // relative dEdx resolution
   Double_t fRange;        // one particle type PID range (in sigmas)
-  ClassDef(AliITSpidESD1,2)   // ITS PID class
+  ClassDef(AliITSpidESD1,3)   // ITS PID class
 };
 
 #endif