SDD drift time (needed for alignment) added to ITS rec points (F. Prino)

diff --git a/ITS/AliITSClusterFinderV2SDD.cxx b/ITS/AliITSClusterFinderV2SDD.cxx
index 28b20a366f4b8ad2beafc5dec13327272e2c9d5e..2547aeda4b7b577f633f7b5f8936069d41c600d9 100644 (file)
--- a/ITS/AliITSClusterFinderV2SDD.cxx
+++ b/ITS/AliITSClusterFinderV2SDD.cxx
@@ -226,7 +226,8 @@ FindClustersSDD(AliBin* bins[2], Int_t nMaxBin, Int_t nzBins,
        Float_t timebin=y-0.5;  // to have time bin in range 0.-255. amd centered on the mid of the bin
        if(s==1) zAnode += GetSeg()->NpzHalf();  // right side has anodes from 256. to 511.
        Float_t zdet = GetSeg()->GetLocalZFromAnode(zAnode);
-       Float_t driftTime = GetSeg()->GetDriftTimeFromTb(timebin) - rsdd->GetTimeOffset();
+       Float_t driftTimeUncorr = GetSeg()->GetDriftTimeFromTb(timebin);
+       Float_t driftTime=driftTimeUncorr-rsdd->GetTimeOffset();
        Float_t driftPathMicron = cal->GetDriftPath(driftTime,zAnode);
        const Double_t kMicronTocm = 1.0e-4; 
        Float_t xdet=(driftPathMicron-GetSeg()->Dx())*kMicronTocm; // xdet is negative
@@ -261,6 +262,7 @@ FindClustersSDD(AliBin* bins[2], Int_t nMaxBin, Int_t nzBins,
 
        AliITSRecPoint cc(milab,hit,info);
        cc.SetType(npeaks);
+       cc.SetDriftTime(driftTimeUncorr);
        if(clusters) new (cl[ncl]) AliITSRecPoint(cc); 
        else {
          fDetTypeRec->AddRecPoint(cc);

diff --git a/ITS/AliITSRecPoint.cxx b/ITS/AliITSRecPoint.cxx
index 3ef645174e633b5f960198d1e1a287476b123dac..b429fd26a3dfcdb54bf5f1c0f06880fd2a2a3405 100644 (file)
--- a/ITS/AliITSRecPoint.cxx
+++ b/ITS/AliITSRecPoint.cxx
@@ -48,7 +48,8 @@ fNz(0),
 fNy(0),
 fChargeRatio(0),
 fType(0),
-fDeltaProb(0)
+fDeltaProb(0),
+fDriftTime(0.)
 {
     // default constructor
 }
@@ -66,7 +67,8 @@ fNz(info[1]),
 fNy(info[0]),
 fChargeRatio(0),
 fType(0),
-fDeltaProb(0)
+fDeltaProb(0),
+fDriftTime(0.)
 {
   //standard constructor used in AliITSClusterFinderV2
 
@@ -118,7 +120,8 @@ fNz(pt.fNz),
 fNy(pt.fNy),
 fChargeRatio(pt.fChargeRatio),
 fType(pt.fType),
-fDeltaProb(pt.fDeltaProb)
+fDeltaProb(pt.fDeltaProb),
+fDriftTime(pt.fDriftTime)
 {
   //Copy constructor
 
@@ -164,7 +167,7 @@ void AliITSRecPoint::Print(ostream *os){
     *os << fXloc << " " << fZloc << " " << fdEdX << " ";
     fmt = os->setf(ios::fixed); // every fixed
     *os << fIndex <<" " << fQ << " "<<fLayer <<" "<<fNz<<" "<<fNy<<" ";
-    *os << fChargeRatio<<" " << fType << " " << fDeltaProb;
+    *os << fChargeRatio<<" " << fType << " " << fDeltaProb << " " << fDriftTime;
     os->flags(fmt); // reset back to old formating.
     return;
 }
@@ -186,7 +189,7 @@ void AliITSRecPoint::Read(istream *is){
     SetVolumeId(lab[0]);// fIsMisalinged = mis;
     *is >> fXloc >> fZloc >> fdEdX;
     *is >> fIndex >> fQ >> fLayer >> fNz >> fNy >> fChargeRatio >> fType;
-    *is >> fDeltaProb;
+    *is >> fDeltaProb >> fDriftTime;
 
     return;
 }

diff --git a/ITS/AliITSRecPoint.h b/ITS/AliITSRecPoint.h
index f88fcb7bcba8ef78cb3d8372b04777c6ba1ef5d8..0a37c8a2e3215e8790650dbde2c12f913d933ddb 100644 (file)
--- a/ITS/AliITSRecPoint.h
+++ b/ITS/AliITSRecPoint.h
@@ -58,7 +58,8 @@ class AliITSRecPoint : public AliCluster {
   void SetPhiR(Float_t y) { fChargeRatio=y; }
   void SetType(Int_t type){ fType=type;}
   void SetDeltaProbability(Float_t prob){fDeltaProb = prob;}
-  
+  void SetDriftTime(Float_t tim) {fDriftTime=tim;}
+ 
   Int_t IsUsed() const {return (fQ<0)?1:0;}
   Float_t GetQ() const {return TMath::Abs(fQ);}
   Int_t GetDetectorIndex() const { return 0x3FF&fIndex; }
@@ -71,7 +72,8 @@ class AliITSRecPoint : public AliCluster {
   Int_t GetNindex() const { return 0xFFC00&fIndex; }  //SSD clusters only
   Int_t GetType() const {return fType;}  // type of the cluster
   Float_t GetDeltaProbability() const{return fDeltaProb;} //probability to belong to the delta ray
-  
+  Float_t GetDriftTime() const{return  fDriftTime;}
+
  protected:
 
   Float_t   fXloc ;        //X of cluster (local coordinates)
@@ -85,9 +87,10 @@ class AliITSRecPoint : public AliCluster {
   Short_t  fNy;       //number of digits in y direction 
   Float_t  fChargeRatio; //charge ratio
   Int_t    fType;         //quality factor of the cluster
-  Float_t  fDeltaProb;    // probability to be deleta electron
+  Float_t  fDeltaProb;    // probability to be delta electron
+  Float_t  fDriftTime;    // drift time in SDD
     
-  ClassDef(AliITSRecPoint,5)  // AliITSRecPoint class
+  ClassDef(AliITSRecPoint,6)  // AliITSRecPoint class
 };
 // Input and output function for standard C++ input/output.
 ostream& operator<<(ostream &os,AliITSRecPoint &source);