From: marian <marian@f7af4fe6-9843-0410-8265-dc069ae4e863>
Date: Wed, 17 Oct 2007 06:40:51 +0000 (+0000)
Subject: Adding things to the "debug" streamer (Marian)
X-Git-Url: http://git.uio.no/git/?a=commitdiff_plain;h=9b42742ad944d8493258d97fab5b78f2c0f5cc20;p=u%2Fmrichter%2FAliRoot.git

Adding things to the "debug" streamer (Marian)
---

diff --git a/TPC/TPCcalib/AliTPCcalibTracksGain.cxx b/TPC/TPCcalib/AliTPCcalibTracksGain.cxx
index e55a4d49bee..f219692f1e3 100644
--- a/TPC/TPCcalib/AliTPCcalibTracksGain.cxx
+++ b/TPC/TPCcalib/AliTPCcalibTracksGain.cxx
@@ -515,6 +515,37 @@ void  AliTPCcalibTracksGain::DumpTrack(AliTPCseed * track){
   
 }
 
+
+
+//_____________________________________________________________________________
+Float_t AliTPCcalibTracksGain::TPCBetheBloch(Float_t p, Float_t mass)
+{
+  //
+  // Bethe-Bloch energy loss formula
+  //	
+  Float_t bg = p /mass;
+  const Double_t kp1=0.76176e-1;
+  const Double_t kp2=10.632;
+  const Double_t kp3=0.13279e-4;
+  const Double_t kp4=1.8631;
+  const Double_t kp5=1.9479;
+
+  Double_t dbg = (Double_t) bg;
+
+  Double_t beta = dbg/TMath::Sqrt(1.+dbg*dbg);
+
+  Double_t aa = TMath::Power(beta,kp4);
+  Double_t bb = TMath::Power(1./dbg,kp5);
+
+  bb=TMath::Log(kp3+bb);
+  
+  return ((Float_t)((kp2-aa-bb)*kp1/aa));
+}
+
+
+
+
+
 Bool_t   AliTPCcalibTracksGain::GetDedx(AliTPCseed * track, Int_t padType, Int_t *rows){
   //
   // GetDedx for given sector for given track
@@ -647,6 +678,9 @@ Bool_t   AliTPCcalibTracksGain::GetDedx(AliTPCseed * track, Int_t padType, Int_t
   }
 
   inonEdge=0;
+  Float_t momenta = track->GetP();
+  Float_t mdedx   = track->GetdEdx();
+
   for (Int_t i=0; i<nclusters; i++){
     Int_t rowSorted = rowIn[index[i]]; 
     AliTPCclusterMI* cluster = track->GetClusterPointer(rowSorted);
@@ -658,9 +692,7 @@ Bool_t   AliTPCcalibTracksGain::GetDedx(AliTPCseed * track, Int_t padType, Int_t
     Float_t dedge    = cluster->GetX()*ktany-TMath::Abs(cluster->GetY());
     Float_t fraction = Float_t(i)/Float_t(nclusters);
     Float_t fraction2= Float_t(inonEdge)/Float_t(nclustersNE);
-    Float_t momenta = track->GetP();
-    Float_t mdedx   = track->GetdEdx();
-    (*fDebugStream)<<"dEdx"<<
+    (*fDebugStream)<<"dEdxCl"<<
       "Cl.="<<cluster<<    //cluster of interest
       "P="<<momenta<<     // track momenta
       "dedx="<<mdedx<<    // mean dedx - corrected for angle
@@ -678,6 +710,24 @@ Bool_t   AliTPCcalibTracksGain::GetDedx(AliTPCseed * track, Int_t padType, Int_t
 
   }
 
+  Float_t prim =  TPCBetheBloch(momenta, 0.1057);   //dedx under muon assumption
+  (*fDebugStream)<<"dEdxTr"<<
+    "P="<<momenta<<     // track momenta
+    "sector="<<lastSector<<// sector to investigate
+    "ncl0="<<nclusters<<    // number of clusters
+    "ncl="<<inonEdge<<      // number of clusters used for dEdx
+    "dedx="<<mdedx<<    // mean dedx - corrected for angle
+    "prim="<<prim<<     // primary ionization parameters
+    "IPad="<<padType<<   // pad type 0..2
+    "xc="<<xcenter<<     // x center of chamber
+    "dedxQ.="<<&dedxQ<<  // dedxQ  - total charge
+    "dedxM.="<<&dedxM<<  // dedxM  - maximal charge
+    "parY.="<<&parY<<      // line fit
+    "parZ.="<<&parZ<<      // line fit
+    "meanPos.="<<&meanPos<< // mean position (dx, dx^2, y,y^2, z, z^2)
+    "\n";
+
+  return 0;
 
 }
 
diff --git a/TPC/TPCcalib/AliTPCcalibTracksGain.h b/TPC/TPCcalib/AliTPCcalibTracksGain.h
index 0e648ebb555..aa6f99270b2 100644
--- a/TPC/TPCcalib/AliTPCcalibTracksGain.h
+++ b/TPC/TPCcalib/AliTPCcalibTracksGain.h
@@ -55,6 +55,9 @@ public :
    TLinearFitter*  GetFitter(UInt_t segment, UInt_t padType, UInt_t fitType);
    Double_t        GetPadLength(Double_t lx);
    Int_t           GetPadType(Double_t lx);
+  //
+  //
+  static Float_t   TPCBetheBloch(Float_t p, Float_t mass=0.1057);