ITS/AliITSRecPoint.cxx

   1 /**************************************************************************
   2  * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
   3  *                                                                        *
   4  * Author: The ALICE Off-line Project.                                    *
   5  * Contributors are mentioned in the code where appropriate.              *
   6  *                                                                        *
   7  * Permission to use, copy, modify and distribute this software and its   *
   8  * documentation strictly for non-commercial purposes is hereby granted   *
   9  * without fee, provided that the above copyright notice appears in all   *
  10  * copies and that both the copyright notice and this permission notice   *
  11  * appear in the supporting documentation. The authors make no claims     *
  12  * about the suitability of this software for any purpose. It is          *
  13  * provided "as is" without express or implied warranty.                  *
  14  **************************************************************************/
  15
  16 /* $Id$ */
  17
  18 ///////////////////////////////////////////////////////////////////////////////
  19 //  Reconstructed space point class for set:ITS
  20 //  Reconstructed points are expressed simultaneously in two different
  21 //  reference frames, both differing from the global system.
  22 //  The first is referred to the sensor (see AliITSsegmentation for the
  23 //  definition) and each point is represented by two coordinates: fXloc and
  24 //  fZloc. This system in the code is referred to as "local"
  25 //  The second is used for tracking (V2, SA and MI versions) and the X axis
  26 //  represents the radial coordinate (this system is, in the bending plane,
  27 //  a rotated system w.r.t. the global reference system).
  28 //  Each reaconstructed point is represented by two coordinates: fY and fZ,
  29 //  inherited from AliCluster. This system in the code is referred to as
  30 //  "trackingV2".
  31 ///////////////////////////////////////////////////////////////////////////////
  32
  33
  34 #include "AliITSRecPoint.h"
  35 #include "AliAlignObj.h"
  36
  37 ClassImp(AliITSRecPoint)
  38
  39 //_____________________________________________________________
  40 AliITSRecPoint::AliITSRecPoint(): AliCluster(),
  41 fXloc(0),
  42 fZloc(0),
  43 fdEdX(0),
  44 fIndex(0),
  45 fQ(0),
  46 fLayer(0),
  47 fNz(0),
  48 fNy(0),
  49 fChargeRatio(0),
  50 fType(0),
  51 fDeltaProb(0)
  52 {
  53     // default constructor
  54 }
  55
  56 //________________________________________________________________________
  57 AliITSRecPoint::AliITSRecPoint(Int_t *lab,Float_t *hit, Int_t *info, Bool_t local):
  58 AliCluster(AliAlignObj::LayerToVolUID((info[2]+AliAlignObj::kSPD1),lab[3]&0x3FF),hit,0,0,lab),
  59 fXloc(0),
  60 fZloc(0),
  61 fdEdX(0),
  62 fIndex(lab[3]),
  63 fQ(hit[4]),
  64 fLayer(info[2]),
  65 fNz(info[1]),
  66 fNy(info[0]),
  67 fChargeRatio(0),
  68 fType(0),
  69 fDeltaProb(0)
  70 {
  71   //standard constructor used in AliITSClusterFinderV2
  72
  73   if (!local) { // Cluster V2
  74     Double_t txyz[3] = {GetX(), GetY(), GetZ()};
  75     Double_t lxyz[3] = {0, 0, 0};
  76     GetTracking2LocalMatrix()->LocalToMaster(txyz,lxyz);
  77     fXloc = lxyz[0]; fZloc = lxyz[2];
  78   }
  79   else {
  80     switch (fLayer) {
  81     case 0:
  82     case 1:
  83       fdEdX = 0;
  84       break;
  85     case 2:
  86     case 3:
  87       fdEdX=fQ*1e-6;
  88       break;
  89     case 4:
  90     case 5:
  91       fdEdX=fQ*2.16;
  92       break;
  93     default:
  94       AliError(Form("Wrong ITS layer %d (0 -> 5)",fLayer));
  95       break;
  96     }
  97     fXloc = hit[0];
  98     fZloc = hit[1];
  99     Double_t lxyz[3] = {fXloc, 0, fZloc};
 100     Double_t txyz[3] = {0, 0, 0};
 101     GetTracking2LocalMatrix()->MasterToLocal(lxyz,txyz);
 102
 103     SetX(0.); SetY(txyz[1]); SetZ(txyz[2]);
 104
 105   }
 106
 107 }
 108
 109 //_______________________________________________________________________
 110 AliITSRecPoint::AliITSRecPoint(const AliITSRecPoint& pt):AliCluster(pt),
 111 fXloc(pt.fXloc),
 112 fZloc(pt.fZloc),
 113 fdEdX(pt.fdEdX),
 114 fIndex(pt.fIndex),
 115 fQ(pt.fQ),
 116 fLayer(pt.fLayer),
 117 fNz(pt.fNz),
 118 fNy(pt.fNy),
 119 fChargeRatio(pt.fChargeRatio),
 120 fType(pt.fType),
 121 fDeltaProb(pt.fDeltaProb)
 122 {
 123   //Copy constructor
 124
 125 }
 126
 127 //______________________________________________________________________
 128 AliITSRecPoint& AliITSRecPoint::operator=(const AliITSRecPoint& source){
 129   // Assignment operator
 130
 131   this->~AliITSRecPoint();
 132   new(this) AliITSRecPoint(source);
 133   return *this;
 134
 135 }
 136
 137 //----------------------------------------------------------------------
 138 void AliITSRecPoint::Print(ostream *os){
 139     ////////////////////////////////////////////////////////////////////////
 140     // Standard output format for this class.
 141     ////////////////////////////////////////////////////////////////////////
 142 #if defined __GNUC__
 143 #if __GNUC__ > 2
 144     ios::fmtflags fmt;
 145 #else
 146     Int_t fmt;
 147 #endif
 148 #else
 149 #if defined __ICC || defined __ECC || defined __xlC__
 150     ios::fmtflags fmt;
 151 #else
 152     Int_t fmt;
 153 #endif
 154 #endif
 155
 156     fmt = os->setf(ios::fixed);  // set fixed floating point output
 157     *os << GetLabel(0) << " " << GetLabel(1) << " " << GetLabel(2) << " ";
 158     *os << fXloc << " " << fZloc << " " << fQ << " ";
 159     fmt = os->setf(ios::scientific); // set scientific for dEdX.
 160     *os << fdEdX << " ";
 161     fmt = os->setf(ios::fixed); // every fixed
 162     *os << GetSigmaY2() << " " << GetSigmaZ2();
 163     os->flags(fmt); // reset back to old formating.
 164     return;
 165 }
 166 //----------------------------------------------------------------------
 167 void AliITSRecPoint::Read(istream *is){
 168 ////////////////////////////////////////////////////////////////////////
 169 // Standard input format for this class.
 170 ////////////////////////////////////////////////////////////////////////
 171
 172   Int_t lab[4];
 173   Float_t hit[5];
 174   lab[3] = 0; // ??
 175   *is >> lab[0] >> lab[1] >> lab[2] >> hit[0] >> hit[1] >> hit[4];
 176   *is >> fdEdX >> hit[2] >> hit[3];
 177   Int_t info[3] = {0,0,0};
 178   AliITSRecPoint rp(lab,hit,info,kTRUE);
 179   *this = rp;
 180
 181   return;
 182 }
 183 //----------------------------------------------------------------------
 184 ostream &operator<<(ostream &os,AliITSRecPoint &p){
 185 ////////////////////////////////////////////////////////////////////////
 186 // Standard output streaming function.
 187 ////////////////////////////////////////////////////////////////////////
 188
 189     p.Print(&os);
 190     return os;
 191 }
 192 //----------------------------------------------------------------------
 193 istream &operator>>(istream &is,AliITSRecPoint &r){
 194 ////////////////////////////////////////////////////////////////////////
 195 // Standard input streaming function.
 196 ////////////////////////////////////////////////////////////////////////
 197
 198     r.Read(&is);
 199     return is;
 200 }
 201 //----------------------------------------------------------------------