Additional protection to avoid rare problem related to the simulation of

[u/mrichter/AliRoot.git] / ITS / AliITSRecPoint.cxx

diff --git a/ITS/AliITSRecPoint.cxx b/ITS/AliITSRecPoint.cxx
index 701dae0e02172c04ec702c914229f6281ee00285..462353f09285085234ab3c45183ea304da3cf91f 100644 (file)
--- a/ITS/AliITSRecPoint.cxx
+++ b/ITS/AliITSRecPoint.cxx
@@ -15,20 +15,105 @@
 
 /* $Id$ */
 
-////////////////////////////////////////////////
-//  Reconstructed point class for set:ITS     //
-////////////////////////////////////////////////
+///////////////////////////////////////////////////////////////////////////////
+//  Reconstructed space point class for set:ITS   
+//  Reconstructed points are expressed simultaneously in two different 
+//  reference frames, both differing from the global system.
+//  The first is referred to the sensor (see AliITSsegmentation for the
+//  definition) and each point is represented by two coordinates: fXloc and
+//  fZloc. This system in the code is referred to as "local"
+//  The second is used for tracking (V2, SA and MI versions) and the X axis 
+//  represents the radial coordinate (this system is, in the bending plane, 
+//  a rotated system w.r.t. the global reference system). 
+//  Each reaconstructed point is represented by two coordinates: fY and fZ, 
+//  inherited from AliCluster. This system in the code is referred to as 
+//  "trackingV2".
+///////////////////////////////////////////////////////////////////////////////
 
 
 #include "AliITSRecPoint.h"
+#include "AliITSgeom.h"
 ClassImp(AliITSRecPoint)
 
-AliITSRecPoint::AliITSRecPoint() {
+//_____________________________________________________________
+AliITSRecPoint::AliITSRecPoint(): AliCluster() {
     // default creator
-    fTracks[0]=fTracks[1]=fTracks[2]=-3; 
-    fX=fZ=fQ=fdEdX=0.;
-    fSigmaX2=fSigmaZ2=0.;
+    fXloc=fZloc=fdEdX=0.;
+    fQ=0; fLayer=0; fNz=fNy=1; fType=0; fDeltaProb=0;fChargeRatio=0;
+    fGeom=0;
 }
+
+//_____________________________________________________________
+AliITSRecPoint::AliITSRecPoint(AliITSgeom* geom): AliCluster() {
+    // default creator
+    fXloc=fZloc=fdEdX=0.;
+    fQ=0; fLayer=0; fNz=fNy=1; fType=0; fDeltaProb=0;fChargeRatio=0;
+    fGeom=geom;
+}
+
+//________________________________________________________________________
+AliITSRecPoint::AliITSRecPoint(Int_t module,AliITSgeom* geom,Int_t *lab,Float_t *hit, Int_t *info):AliCluster(lab,hit){
+  //standard constructor used in AliITSClusterFinderV2
+  fIndex=lab[3];
+  fQ=hit[4];
+  fNy    = info[0];
+  fNz    = info[1];
+  fLayer = info[2];
+  fChargeRatio = 0;
+  fType=0;
+  fDeltaProb=0.;
+  
+  fGeom = geom;
+  fGeom->TrackingV2ToDetL(module,fY,fZ,fXloc,fZloc);
+  if(module<fGeom->GetStartSDD()) fdEdX=0.;
+  if(module>=fGeom->GetStartSDD() && module<fGeom->GetStartSSD()){
+    fdEdX=fQ*1e-6;
+  }
+  if(module>=fGeom->GetStartSSD()) fdEdX=fQ*2.16;
+  
+  
+}
+//_______________________________________________________________________
+AliITSRecPoint::AliITSRecPoint(const AliITSRecPoint& pt):AliCluster(pt){
+  //Copy constructor
+  fXloc = pt.fXloc;
+  fZloc = pt.fZloc;
+  fdEdX = pt.fdEdX;
+  fIndex= pt.fIndex;
+  fQ = pt.fQ;
+  fNy = pt.fNy;
+  fNz = pt.fNz;
+  fLayer = pt.fLayer;
+  fChargeRatio = pt.fChargeRatio;
+  fType = pt.fType;
+  fDeltaProb = pt.fDeltaProb;
+  fGeom = pt.fGeom;
+
+}
+
+//______________________________________________________________________
+AliITSRecPoint& AliITSRecPoint::operator=(const AliITSRecPoint& /* source */){
+  // Assignment operator
+  // Assignment is not allowed. The method is protected to avoid misuse.
+  Fatal("= operator","Assignment operator not allowed\n");
+  return *this;
+}
+
+//________________________________________________________________________
+AliITSRecPoint::AliITSRecPoint(Int_t *lab,Float_t *hit, Int_t *info):AliCluster(lab,hit){
+  //standard constructor used in AliITSClusterFinderV2
+  fXloc=fZloc=fdEdX=0.;
+  fIndex=lab[3];
+  fQ=hit[4];
+  fNy    = info[0];
+  fNz    = info[1];
+  fLayer = info[2];
+  fChargeRatio = 0;
+  fType=0;
+  fDeltaProb=0.;  
+  fGeom = 0;
+}
+
 //----------------------------------------------------------------------
 void AliITSRecPoint::Print(ostream *os){
     ////////////////////////////////////////////////////////////////////////
@@ -41,7 +126,7 @@ void AliITSRecPoint::Print(ostream *os){
     Int_t fmt;
 #endif
 #else
-#if defined __ICC || defined __ECC || defined _xlc_
+#if defined __ICC || defined __ECC || defined __xlC__
     ios::fmtflags fmt;
 #else
     Int_t fmt;
@@ -50,11 +135,11 @@ void AliITSRecPoint::Print(ostream *os){
  
     fmt = os->setf(ios::fixed);  // set fixed floating point output
     *os << fTracks[0]<< " " << fTracks[1] << " " << fTracks[2] << " ";
-    *os << fX << " " << fZ << " " << fQ << " ";
+    *os << fXloc << " " << fZloc << " " << fQ << " ";
     fmt = os->setf(ios::scientific); // set scientific for dEdX.
     *os << fdEdX << " ";
     fmt = os->setf(ios::fixed); // every fixed
-    *os << fSigmaX2 << " " << fSigmaZ2;
+    *os << fSigmaY2 << " " << fSigmaZ2;
     os->flags(fmt); // reset back to old formating.
     return;
 }
@@ -65,8 +150,8 @@ void AliITSRecPoint::Read(istream *is){
 ////////////////////////////////////////////////////////////////////////
  
 
-    *is >> fTracks[0] >> fTracks[1] >> fTracks[2] >> fX >> fZ >> fQ;
-    *is >> fdEdX >> fSigmaX2 >> fSigmaZ2;
+    *is >> fTracks[0] >> fTracks[1] >> fTracks[2] >> fXloc >> fZloc >> fQ;
+    *is >> fdEdX >> fSigmaY2 >> fSigmaZ2;
     return;
 }
 //----------------------------------------------------------------------