* provided "as is" without express or implied warranty. *
**************************************************************************/
+/* $Id$ */
+
//-----------------------------------------------------------------
// Implementation of the ITS PID class
// Very naive one... Should be made better by the detector experts...
#include "AliITSpidESD.h"
#include "AliESDtrack.h"
+#include "AliExternalTrackParam.h"
ClassImp(AliITSpidESD)
//Default constructor
}
-Double_t AliITSpidESD::Bethe(Double_t bg) {
- //
- // This is the Bethe-Bloch function normalised to 1 at the minimum
- //
- Double_t bg2=bg*bg;
- Double_t bethe;
- /*if (bg<3.5)
- bethe=(1.+ bg2)/bg2*(log(5940*bg2) - bg2/(1.+ bg2));
- else*/ // not 100% clear why...
- bethe=(1.+ bg2)/bg2*(TMath::Log(3.5*5940*bg) - bg2/(1.+ bg2));
- return bethe/11.091;
-}
+Double_t AliITSpidESD::Bethe(Double_t p,Double_t mass) {
+ // returns AliExternalTrackParam::BetheBloch normalized to 1 at the minimum
+ Double_t density=2.33; // g/cm3
+ Double_t thickness=0.03; // cm
+ Double_t meanMIPSi=116.24; // keV in 300 microns of Si
+ Double_t conv=density*1E6*thickness/meanMIPSi;
+ Float_t betagamma=p/mass;
+ return conv*AliExternalTrackParam::BetheBlochSolid(betagamma);
+
+}