TOF/AliTOFSDigit.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 //                                                                         //
  20 //  TOF sdigit: member variables                                           //
  21 //  fSector  : TOF sector                                                  //
  22 //  fPlate   : TOF plate                                                   //
  23 //  fStrip   : strips number                                               //
  24 //  fPadx    : pad number along x                                          //
  25 //  fPadz    : pad number along z                                          //
  26 //  fTdc     : TArrayI of TDC values                                       //
  27 //  fAdc     : TArrayI of ADC values                                       //
  28 //                                                                         //
  29 //  Getters, setters and member functions  defined here                    //
  30 //                                                                         //
  31 // -- Authors: F. Pierella, A. Seganti, D. Vicinanza                       //
  32 //_________________________________________________________________________//
  33
  34 #include "AliLog.h"
  35
  36 #include "AliTOFGeometry.h"
  37 #include "AliTOFSDigit.h"
  38
  39 ClassImp(AliTOFSDigit)
  40
  41 ////////////////////////////////////////////////////////////////////////
  42 AliTOFSDigit::AliTOFSDigit():
  43   fSector(-1),
  44   fPlate(-1),
  45   fStrip(-1),
  46   fPadx(-1),
  47   fPadz(-1),
  48   fNDigits(0),
  49   fTdc(0x0),
  50   fAdc(0x0),
  51   fTracks(0x0)
  52 {
  53   //
  54   // default ctor
  55   //
  56 }
  57
  58 ////////////////////////////////////////////////////////////////////////
  59 AliTOFSDigit::AliTOFSDigit(Int_t tracknum, Int_t * const vol,Int_t * const digit):
  60   TObject(),
  61   fSector(-1),
  62   fPlate(-1),
  63   fStrip(-1),
  64   fPadx(-1),
  65   fPadz(-1),
  66   fNDigits(0),
  67   fTdc(0x0),
  68   fAdc(0x0),
  69   fTracks(0x0)
  70 {
  71   //
  72   // Constructor of digit object
  73   //
  74
  75   fSector = vol[0];
  76   fPlate  = vol[1];
  77   fStrip  = vol[2];
  78   fPadx   = vol[3];
  79   fPadz   = vol[4];
  80   fNDigits = 1;
  81   fTdc = new TArrayI(fNDigits);
  82   (*fTdc)[0] = digit[0];
  83   fAdc = new TArrayI(fNDigits);
  84   (*fAdc)[0] = digit[1];
  85   fTracks = new TArrayI(kMAXDIGITS*fNDigits);
  86   (*fTracks)[0] = tracknum;
  87   for (Int_t i = 1; i <kMAXDIGITS*fNDigits; i++) {
  88     (*fTracks)[i] = -1;
  89   }
  90 }
  91
  92 ////////////////////////////////////////////////////////////////////////
  93 AliTOFSDigit::AliTOFSDigit(const AliTOFSDigit & digit):
  94   TObject(digit),
  95   fSector(digit.fSector),
  96   fPlate(digit.fPlate),
  97   fStrip(digit.fStrip),
  98   fPadx(digit.fPadx),
  99   fPadz(digit.fPadz),
 100   fNDigits(digit.fNDigits),
 101   fTdc(0x0),
 102   fAdc(0x0),
 103   fTracks(0x0)
 104 {
 105   //
 106   // copy ctor for AliTOFSDigit object
 107   //
 108   fTdc = new TArrayI(*digit.fTdc);
 109   fAdc = new TArrayI(*digit.fAdc);
 110   fTracks = new TArrayI(*digit.fTracks);
 111 }
 112
 113 ////////////////////////////////////////////////////////////////////////
 114 AliTOFSDigit& AliTOFSDigit::operator=(const AliTOFSDigit & digit)
 115 {
 116   //
 117   // copy ctor for AliTOFSDigit object
 118   //
 119
 120   if (this == &digit)
 121     return *this;
 122
 123   TObject::operator=(digit);
 124   fSector = digit.fSector;
 125   fPlate  = digit.fPlate;
 126   fStrip  = digit.fStrip;
 127   fPadx   = digit.fPadx;
 128   fPadz   = digit.fPadz;
 129   fNDigits = digit.fNDigits;
 130   fTdc = digit.fTdc;
 131   fAdc = digit.fAdc;
 132   fTracks = digit.fTracks;
 133   return *this;
 134
 135 }
 136
 137 ////////////////////////////////////////////////////////////////////////
 138 AliTOFSDigit::AliTOFSDigit(Int_t sector, Int_t plate, Int_t strip, Int_t padx,
 139                            Int_t padz, Int_t tdc, Int_t adc):
 140   fSector(sector),
 141   fPlate(plate),
 142   fStrip(strip),
 143   fPadx(padx),
 144   fPadz(padz),
 145   fNDigits(1),
 146   fTdc(0x0),
 147   fAdc(0x0),
 148   fTracks(0x0)
 149 {
 150   //
 151   // Constructor for sdigit
 152   //
 153   fTdc = new TArrayI(fNDigits);
 154   (*fTdc)[0] = tdc;
 155   fAdc = new TArrayI(fNDigits);
 156   (*fAdc)[0] = adc;
 157   // no tracks were specified, set them to -1
 158   fTracks = new TArrayI(kMAXDIGITS*fNDigits);
 159   for (Int_t i = 0; i <kMAXDIGITS*fNDigits; i++) {
 160     (*fTracks)[i] = -1;
 161   }
 162 }
 163
 164 ////////////////////////////////////////////////////////////////////////
 165 void AliTOFSDigit::GetLocation(Int_t *Loc) const
 166 {
 167   //
 168   // Get the coordinates of the digit
 169   // in terms of Sector - Plate - Strip - Pad
 170   //
 171
 172   Loc[0]=fSector;
 173   Loc[1]=fPlate;
 174   Loc[2]=fStrip;
 175   Loc[3]=fPadx;
 176   Loc[4]=fPadz;
 177 }
 178
 179 ////////////////////////////////////////////////////////////////////////
 180 void AliTOFSDigit::Update(Float_t tdcbin, Int_t tdc, Int_t adc, Int_t track)
 181 {
 182   //
 183   // Add charge and track
 184   //
 185
 186   Int_t sameTime = -1;
 187   Float_t tdcwindow = AliTOFGeometry::DeadTime()/tdcbin;
 188   for (Int_t i = 0; i < fNDigits; i++) {
 189     if (TMath::Abs(tdc-fTdc->At(i)) < tdcwindow) {
 190       sameTime = i;
 191       break;
 192     }
 193   }
 194
 195   if (sameTime >= 0) { // another time measurement happens during the
 196                        // dead time of the hit pad => it corresponds
 197                        // to the same time measurement
 198     (*fAdc)[sameTime] += adc;
 199
 200     // update track index array in case the current digit track index
 201     // is different from -1
 202     if (track!=-1) {
 203
 204       //Find the first -1 value of the track index array and replace
 205       //it by the current digit track index
 206       for (Int_t iTrack=0; iTrack<kMAXDIGITS; iTrack++) {
 207         if (track==(*fTracks)[sameTime*kMAXDIGITS+iTrack]) break;
 208         if ((*fTracks)[sameTime*kMAXDIGITS+iTrack] == -1) {
 209           (*fTracks)[sameTime*kMAXDIGITS+iTrack] = track;
 210           break;
 211         }
 212         // write warning about many tracks going to this pad at same time
 213         if (iTrack == kMAXDIGITS-1) {
 214           AliDebug(1,Form("Update: Many different tracks in the same TOF pad"
 215                           " (%2d %1d %2d %1d %2d)\n",
 216                           fSector,fPlate,fStrip,fPadz,fPadx));
 217           //ToAliWarning(PrintPad());
 218         }
 219       }
 220     }
 221
 222   } else { // they are two different time measurements
 223
 224     // add new time slot
 225     fNDigits++;
 226     fTdc->Set(fNDigits);
 227     (*fTdc)[fNDigits-1] = tdc;
 228     fAdc->Set(fNDigits);
 229     (*fAdc)[fNDigits-1] = adc;
 230     fTracks->Set(fNDigits*kMAXDIGITS);
 231     (*fTracks)[(fNDigits-1)*kMAXDIGITS] = track;
 232     for (Int_t i = 1; i <kMAXDIGITS; i++)
 233       (*fTracks)[(fNDigits-1)*kMAXDIGITS+i] = -1;
 234
 235   }
 236
 237 }
 238
 239 ////////////////////////////////////////////////////////////////////////
 240 void AliTOFSDigit::Update(AliTOFSDigit * const sdig)
 241 {
 242
 243   //
 244   // Perform the sum with sdig
 245   //
 246
 247   Float_t tdcbin = AliTOFGeometry::TdcBinWidth();// [ps] hardwired for the time being
 248
 249   Int_t track = -1;
 250   Int_t adc = -1;
 251   Int_t tdc = -1;
 252
 253   // start loop on all sdig locations
 254   Int_t nlocations = sdig->GetNDigits();
 255   for (Int_t j = 0; j < nlocations; j++) {
 256     tdc = (Int_t)sdig->GetTdc(j);
 257     adc = (Int_t)sdig->GetAdc(j);
 258
 259     // getting here only the first track number
 260     //Int_t track = GetTrack(j,0);
 261
 262     // getting here all the track numbers
 263     for (Int_t iTrack = 0; iTrack<kMAXDIGITS; iTrack++) {
 264       track = sdig->GetTrack(j,iTrack);
 265       Update(tdcbin, tdc, adc, track);
 266     } // end loop on tracks
 267
 268   } // end loop on sdig locations
 269
 270
 271 }
 272
 273 ////////////////////////////////////////////////////////////////////////
 274 AliTOFSDigit::~AliTOFSDigit()
 275 {
 276   //
 277   // dtor
 278   //
 279   delete fTdc;
 280   delete fAdc;
 281   delete fTracks;
 282 }
 283
 284 ////////////////////////////////////////////////////////////////////////
 285
 286 Int_t AliTOFSDigit::GetTotPad() const
 287 {
 288   //
 289   // Get the "total" index of the pad inside a Sector
 290   // starting from the digits data.
 291   //
 292
 293   Int_t pad = AliTOFGeometry::NpadZ()*fPadx + fPadz;
 294   Int_t before=0;
 295
 296   switch(fPlate){
 297   case 0:
 298     //before = 0;
 299     break;
 300   case 1:
 301     before = AliTOFGeometry::NStripC();
 302     break;
 303   case 2:
 304     before = AliTOFGeometry::NStripB() +   AliTOFGeometry::NStripC();
 305     break;
 306   case 3:
 307     before = AliTOFGeometry::NStripA() +   AliTOFGeometry::NStripB() + AliTOFGeometry::NStripC();
 308     break;
 309   case 4:
 310     before = AliTOFGeometry::NStripA() + 2*AliTOFGeometry::NStripB() + AliTOFGeometry::NStripC();
 311     break;
 312   }
 313
 314   Int_t strip = fStrip + before;
 315   Int_t padTot = AliTOFGeometry::NpadXStrip()*strip + pad;
 316   return padTot;
 317 }