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 #include <TGeoMatrix.h>
  34 #include "AliITSRecPoint.h"
  35 #include "AliAlignObj.h"
  36
  37 using std::ios;
  38 ClassImp(AliITSRecPoint)
  39
  40 //_____________________________________________________________
  41 AliITSRecPoint::AliITSRecPoint():
  42 AliCluster(),
  43 fXloc(0),
  44 fZloc(0),
  45 fdEdX(0),
  46 fIndex(0),
  47 fQ(0),
  48 fLayer(0),
  49 fNz(0),
  50 fNy(0),
  51 fChargeRatio(0),
  52 fType(0),
  53 fDeltaProb(0),
  54 fDriftTime(0.),
  55 fDriftSide(0)
  56 {
  57     // default constructor
  58 }
  59
  60 //________________________________________________________________________
  61 AliITSRecPoint::AliITSRecPoint(Int_t *lab,Float_t *hit, Int_t *info, Bool_t local):
  62 AliCluster(AliGeomManager::LayerToVolUID((info[2]+AliGeomManager::kSPD1),lab[3]&0x3FF),hit,0,0,lab),
  63 fXloc(0),
  64 fZloc(0),
  65 fdEdX(0),
  66 fIndex(lab[3]),
  67 fQ(hit[4]),
  68 fLayer(info[2]),
  69 fNz(info[1]),
  70 fNy(info[0]),
  71 fChargeRatio(0),
  72 fType(0),
  73 fDeltaProb(0),
  74 fDriftTime(0.),
  75 fDriftSide(0)
  76 {
  77   //standard constructor used in AliITSClusterFinderV2
  78
  79   if (!local) { // Cluster V2
  80     Double_t txyz[3] = {GetX(), GetY(), GetZ()};
  81     Double_t lxyz[3] = {0, 0, 0};
  82     GetTracking2LocalMatrix()->LocalToMaster(txyz,lxyz);
  83     fXloc = lxyz[0]; fZloc = lxyz[2];
  84     if(fLayer==4) hit[5]=-hit[5];
  85     if( (fLayer==4) || (fLayer==5) ) SetSigmaYZ(hit[5]);
  86   }
  87   else {
  88     switch (fLayer) {
  89     case 0:
  90     case 1:
  91       fdEdX = 0;
  92       break;
  93     case 2:
  94     case 3:
  95       fdEdX=fQ*1e-6;
  96       break;
  97     case 4:
  98       fdEdX=fQ*2.16;
  99       SetSigmaYZ(hit[5]);
 100       break;
 101     case 5:
 102       fdEdX=fQ*2.16;
 103       hit[5]=-hit[5];
 104       SetSigmaYZ(hit[5]);
 105       break;
 106     default:
 107       AliError(Form("Wrong ITS layer %d (0 -> 5)",fLayer));
 108       break;
 109     }
 110     fXloc = hit[0];
 111     fZloc = hit[1];
 112     Double_t lxyz[3] = {fXloc, 0, fZloc};
 113     Double_t txyz[3] = {0, 0, 0};
 114     GetTracking2LocalMatrix()->MasterToLocal(lxyz,txyz);
 115
 116     SetX(0.); SetY(txyz[1]); SetZ(txyz[2]);
 117
 118   }
 119
 120 }
 121
 122 //_______________________________________________________________________
 123 AliITSRecPoint::AliITSRecPoint(const AliITSRecPoint& pt):
 124   AliCluster(pt),
 125   fXloc(pt.fXloc),
 126   fZloc(pt.fZloc),
 127   fdEdX(pt.fdEdX),
 128   fIndex(pt.fIndex),
 129   fQ(pt.fQ),
 130   fLayer(pt.fLayer),
 131   fNz(pt.fNz),
 132   fNy(pt.fNy),
 133   fChargeRatio(pt.fChargeRatio),
 134   fType(pt.fType),
 135   fDeltaProb(pt.fDeltaProb),
 136   fDriftTime(pt.fDriftTime),
 137   fDriftSide(pt.fDriftSide)
 138 {
 139   //Copy constructor
 140
 141 }
 142
 143 //______________________________________________________________________
 144 AliITSRecPoint& AliITSRecPoint::operator=(const AliITSRecPoint& source){
 145   // Assignment operator
 146
 147   this->~AliITSRecPoint();
 148   new(this) AliITSRecPoint(source);
 149   return *this;
 150
 151 }
 152
 153 //----------------------------------------------------------------------
 154 void AliITSRecPoint::Print(ostream *os){
 155     ////////////////////////////////////////////////////////////////////////
 156     // Standard output format for this class.
 157     ////////////////////////////////////////////////////////////////////////
 158 #if defined __GNUC__
 159 #if __GNUC__ > 2
 160     ios::fmtflags fmt;
 161 #else
 162     Int_t fmt;
 163 #endif
 164 #else
 165 #if defined __ICC || defined __ECC || defined __xlC__
 166     ios::fmtflags fmt;
 167 #else
 168     Int_t fmt;
 169 #endif
 170 #endif
 171
 172     fmt = os->setf(ios::fixed);  // set fixed floating point output
 173     *os << GetLabel(0) << " " << GetLabel(1) << " " << GetLabel(2) << " ";
 174     fmt = os->setf(ios::scientific); // set scientific for dEdX.
 175     *os << GetX() <<" " << GetY() << " " << GetZ() << " " ;
 176     *os << GetSigmaY2() << " " << GetSigmaZ2() << " " << GetSigmaYZ() << " ";
 177     fmt = os->setf(ios::fixed);
 178     *os << GetVolumeId() << " "<< Misalign() /*fIsMisaligned*/ << " ";
 179     fmt = os->setf(ios::scientific); // set scientific for dEdX.
 180     *os << fXloc << " " << fZloc << " " << fdEdX << " ";
 181     fmt = os->setf(ios::fixed); // every fixed
 182     *os << fIndex <<" " << fQ << " "<<fLayer <<" "<<fNz<<" "<<fNy<<" ";
 183     *os << fChargeRatio<<" " << fType << " " << fDeltaProb << " " << fDriftTime<< " " << fDriftSide;
 184     os->flags(fmt); // reset back to old formating.
 185     return;
 186 }
 187
 188 //----------------------------------------------------------------------
 189 Int_t AliITSRecPoint::GetNpixels() const {
 190 //
 191 // returns the number of pixels used for the SPD clusters
 192 //
 193
 194  if(fLayer > 1) return -1;
 195  else return fType;
 196
 197 }
 198
 199 //----------------------------------------------------------------------
 200 Int_t AliITSRecPoint::GetSPDclusterType() const {
 201 //
 202 // returns an Int_t with encoded information on cluster size
 203 //    type <= 16: cluster type identifier according to conventional numbering
 204 //    type >  16: Npixels+1000*Ny+1000000*Nz
 205 //
 206
 207  Int_t type = -1;
 208  if(fLayer > 1) return type;
 209  else {
 210
 211    switch (fType) {
 212      case 1 : type = 1 ;break;
 213      case 2 : if(fNy == 2) type = 2;
 214               else         type = 3;
 215               break;
 216      case 3 : if(fNy == 3)      type = 4;
 217               else if(fNz == 3) type = 6;
 218               else              type = 5;
 219               break;
 220      case 4 : if(fNz == 1)                  type = 7;
 221               else if(fNz == 2 && fNy == 2) type = 8;
 222               else if(fNy == 2 && fNz == 3) type = 11;
 223               else if(fNy == 3 && fNz == 2) type = 9;
 224               else                          type = 15;
 225               break;
 226      case 5 : if(fNy == 3 && fNz == 2)  type = 10;
 227               if(fNy == 2 && fNz == 3 ) type = 12;
 228               if(fNy == 5)              type = 16;
 229               else                      type = fType+1000*fNy+1000000*fNz;
 230               break;
 231      case 6 : if(fNy ==3 && fNz == 2) type = 13;
 232               if(fNy ==2 && fNz == 3) type = 14;
 233               else                    type = fType+1000*fNy+1000000*fNz;
 234               break;
 235      default: type = fType+1000*fNy+1000000*fNz;
 236               break;
 237    }
 238
 239    return type;
 240  }
 241 }
 242
 243 //----------------------------------------------------------------------
 244 Int_t AliITSRecPoint::GetSDDclusterType() const {
 245 // returns an Int_t with encoded information on cluster size
 246 // Byte1 = fNz Byte0=fNy, other two bytes empty for extra information
 247 // max. allowed cluster size = 255
 248   Int_t typ=(fNz&0xFF)<<8;
 249   typ+=fNy&0xFF;
 250   if(fDriftSide==1) typ+=1<<16;
 251   return typ;
 252 }
 253 //----------------------------------------------------------------------
 254 void  AliITSRecPoint::DecodeSDDclusterType(Int_t cluType, Int_t &cluSizAn, Int_t& cluSizTb, Int_t &drSide){
 255 // Extract cluster sizes and drift side from cluster type
 256   cluSizTb=cluType&0xFF;
 257   cluSizAn=(cluType>>8)&0xFF;
 258   drSide=(cluType>>16);
 259   return;
 260 }
 261 //----------------------------------------------------------------------
 262 Int_t AliITSRecPoint::GetSSDclusterType() const {
 263 // returns an Int_t with encoded information on cluster size
 264 // Byte1 = fNz Byte0=fNy, other two bytes empty for extra information
 265 // max. allowed cluster size = 255
 266   Int_t typ=(fNz&0xFF)<<8;
 267   typ+=fNy&0xFF;
 268   return typ;
 269 }
 270
 271 //----------------------------------------------------------------------
 272 void AliITSRecPoint::Read(istream *is){
 273 ////////////////////////////////////////////////////////////////////////
 274 // Standard input format for this class.
 275 ////////////////////////////////////////////////////////////////////////
 276     Bool_t mis;
 277     Int_t lab[4];
 278     Float_t hit[6];
 279     lab[3] = 0; // ??
 280     *is >> lab[0] >> lab[1] >> lab[2];
 281     SetLabel(lab[0],0); SetLabel(lab[1],1); SetLabel(lab[2],2);
 282     *is >> hit[0] >> hit[1] >> hit[2] >> hit[3] >> hit[4] >> hit[5];
 283     SetX(hit[0]);SetY(hit[1]);SetZ(hit[2]);SetSigmaY2(hit[3]);
 284     SetSigmaZ2(hit[4]);//fSigmaYZ=hit[5];
 285     *is >> lab[0] >> mis;
 286     SetVolumeId(lab[0]);// fIsMisalinged = mis;
 287     *is >> fXloc >> fZloc >> fdEdX;
 288     *is >> fIndex >> fQ >> fLayer >> fNz >> fNy >> fChargeRatio >> fType;
 289     *is >> fDeltaProb >> fDriftTime >> fDriftSide;
 290
 291     return;
 292 }
 293 //----------------------------------------------------------------------
 294 ostream &operator<<(ostream &os,AliITSRecPoint &p){
 295 ////////////////////////////////////////////////////////////////////////
 296 // Standard output streaming function.
 297 ////////////////////////////////////////////////////////////////////////
 298
 299     p.Print(&os);
 300     return os;
 301 }
 302 //----------------------------------------------------------------------
 303 istream &operator>>(istream &is,AliITSRecPoint &r){
 304 ////////////////////////////////////////////////////////////////////////
 305 // Standard input streaming function.
 306 ////////////////////////////////////////////////////////////////////////
 307
 308     r.Read(&is);
 309     return is;
 310 }
 311 //----------------------------------------------------------------------