ITS/AliITSRawStreamSPD.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 /// This class provides access to ITS SPD digits in raw data.
  21 ///
  22 ///////////////////////////////////////////////////////////////////////////////
  23
  24 #include "AliITSRawStreamSPD.h"
  25 #include "AliRawReader.h"
  26 #include "AliLog.h"
  27
  28 ClassImp(AliITSRawStreamSPD)
  29
  30
  31   // this map has to change, waiting for the new geometry
  32 const Int_t AliITSRawStreamSPD::fgkDDLModuleMap[kDDLsNumber][kModulesPerDDL] = {
  33   { 0, 1, 4, 5, 80, 81, 84, 85, 88, 89, 92, 93},
  34   { 8, 9,12,13, 96, 97,100,101,104,105,108,109},
  35   {16,17,20,21,112,113,116,117,120,121,124,125},
  36   {24,25,28,29,128,129,132,133,136,137,140,141},
  37   {32,33,36,37,144,145,148,149,152,153,156,157},
  38   {40,41,44,45,160,161,164,165,168,169,172,173},
  39   {48,49,52,53,176,177,180,181,184,185,188,189},
  40   {56,57,60,61,192,193,196,197,200,201,204,205},
  41   {64,65,68,69,208,209,212,213,216,217,220,221},
  42   {72,73,76,77,224,225,228,229,232,233,236,237},
  43   { 2, 3, 6, 7, 82, 83, 86, 87, 90, 91, 94, 95},
  44   {10,11,14,15, 98, 99,102,103,106,107,110,111},
  45   {18,19,22,23,114,115,118,119,122,123,126,127},
  46   {26,27,30,31,130,131,134,135,138,139,142,143},
  47   {34,35,38,39,146,147,150,151,154,155,158,159},
  48   {42,43,46,47,162,163,166,167,170,171,174,175},
  49   {50,51,54,55,178,179,182,183,186,187,190,191},
  50   {58,59,62,63,194,195,198,199,202,203,206,207},
  51   {66,67,70,71,210,211,214,215,218,219,222,223},
  52   {74,75,78,79,226,227,230,231,234,235,238,239}
  53 };
  54
  55
  56 AliITSRawStreamSPD::AliITSRawStreamSPD(AliRawReader* rawReader) :
  57   AliITSRawStream(rawReader),
  58   fEventNumber(-1),fChipAddr(0),fHalfStaveNr(0),fCol(0),fRow(0),
  59   fData(0),fOffset(0),fHitCount(0),
  60   fDataChar1(0),fDataChar2(0),fDataChar3(0),fDataChar4(0),
  61   fFirstWord(kTRUE),fPrevEventId(0xffffffff)
  62 {
  63   // create an object to read ITS SPD raw digits
  64   fRawReader->Select("ITSSPD");
  65   // reset calib header words
  66   for (UInt_t iword=0; iword<kCalHeadLenMax; iword++) {
  67     fCalHeadWord[iword]=0xffffffff;
  68   }
  69   NewEvent();
  70 }
  71
  72 Bool_t AliITSRawStreamSPD::ReadNextShort()
  73 {
  74   // read next 16 bit word into fData
  75   if (fFirstWord) {
  76     Bool_t b1 = fRawReader->ReadNextChar(fDataChar1);
  77     if (!b1) return kFALSE;
  78     Bool_t  b2, b3, b4;
  79     b2 = fRawReader->ReadNextChar(fDataChar2);
  80     b3 = fRawReader->ReadNextChar(fDataChar3);
  81     b4 = fRawReader->ReadNextChar(fDataChar4);
  82     if (!(b2 && b3 && b4)) {
  83       return kFALSE;
  84     }
  85     fData = fDataChar3+(fDataChar4<<8);
  86     if ((*fRawReader->GetEventId())!=fPrevEventId) { // if new event...
  87       NewEvent();
  88       fPrevEventId=(*fRawReader->GetEventId());
  89     }
  90     fFirstWord=kFALSE;
  91   }
  92   else {
  93     fFirstWord=kTRUE;
  94     fData = fDataChar1+(fDataChar2<<8);
  95   }
  96
  97   return kTRUE;
  98 }
  99
 100 Bool_t AliITSRawStreamSPD::ReadNextInt()
 101 {
 102   // reads next 32 bit into fDataChar1..4
 103   // (if first 16 bits read already, just completes the present word)
 104   if (fFirstWord) {
 105     if (ReadNextShort() && ReadNextShort()) {
 106       return kTRUE;
 107     }
 108   }
 109   else {
 110     if (ReadNextShort()) {
 111       return kTRUE;
 112     }
 113   }
 114   return kFALSE;
 115 }
 116
 117 void AliITSRawStreamSPD::NewEvent()
 118 {
 119   // call this to reset flags for a new event
 120   for (UInt_t eqId=0; eqId<20; eqId++) {
 121     fCalHeadRead[eqId]=kFALSE;
 122   }
 123   fEventNumber=-1;
 124 }
 125
 126 Bool_t AliITSRawStreamSPD::ReadCalibHeader()
 127 {
 128   // read the extra calibration header
 129   // returns kTRUE if the header is present and has length > 0
 130
 131   Int_t ddlID = fRawReader->GetDDLID();
 132   if (ddlID==-1) { // we may need to read one word to get the blockAttr
 133     if (!ReadNextShort()) return kFALSE;
 134     ddlID = fRawReader->GetDDLID();
 135   }
 136   UChar_t attr = fRawReader->GetBlockAttributes();
 137   if (ddlID>=0 && ddlID<20) fCalHeadRead[ddlID]=kTRUE;
 138   if ((attr & 0x40) == 0x40) { // is the header present?
 139     if (ReadNextInt()) {
 140       // length of cal header:
 141       UInt_t calLen = fDataChar1+(fDataChar2<<8)+(fDataChar3<<16)+(fDataChar4<<24);
 142       if (calLen>kCalHeadLenMax) {
 143         fRawReader->AddMajorErrorLog(kCalHeaderLengthErr,Form("Header length %d > max = %d",calLen,kCalHeadLenMax));
 144         AliWarning(Form("Header length problem. %d > %d (max)",calLen,kCalHeadLenMax));
 145         return kFALSE;
 146       }
 147       else if (calLen>0) {
 148         for (UInt_t iword=0; iword<calLen; iword++) {
 149           if (ReadNextInt()) {
 150             fCalHeadWord[iword] = fDataChar1+(fDataChar2<<8)+(fDataChar3<<16)+(fDataChar4<<24);
 151           }
 152           else {
 153             fRawReader->AddMajorErrorLog(kCalHeaderLengthErr,"header length problem");
 154             AliWarning("header length problem");
 155             return kFALSE;
 156           }
 157         }
 158         return kTRUE;
 159       }
 160     }
 161   }
 162
 163   return kFALSE;
 164 }
 165
 166 Bool_t AliITSRawStreamSPD::Next()
 167 {
 168 // read the next raw digit
 169 // returns kFALSE if there is no digit left
 170
 171   Int_t ddlID=-1;
 172   fPrevModuleID = fModuleID;
 173
 174   while (ReadNextShort()) {
 175
 176     ddlID = fRawReader->GetDDLID();
 177     if (ddlID>=0 && ddlID<20) {
 178       if (!fCalHeadRead[ddlID]) {
 179         ReadCalibHeader();
 180       }
 181     }
 182     else {
 183       fRawReader->AddMajorErrorLog(kDDLNumberErr,Form("Wrong DDL number %d",ddlID));
 184       AliWarning(Form("DDL number error (= %d) , setting it to 19", ddlID));
 185       ddlID=19;
 186     }
 187
 188     if ((fData & 0xC000) == 0x4000) {         // header
 189       fHitCount = 0;
 190       UShort_t eventNumber = (fData >> 4) & 0x007F;
 191       if (fEventNumber < 0) {
 192         fEventNumber = eventNumber;
 193       }
 194       else if (eventNumber != fEventNumber) {
 195         fRawReader->AddMajorErrorLog(kEventNumberErr,Form("Reading event number %d instead of %d",eventNumber,fEventNumber));
 196         AliWarning(Form("mismatching event numbers: %d != %d",eventNumber, fEventNumber));
 197       }
 198       fChipAddr = fData & 0x000F;
 199       if (fChipAddr>9) {
 200         fRawReader->AddMajorErrorLog(kChipAddrErr,Form("Overflow chip address %d - set to 9",fChipAddr));
 201         AliWarning(Form("overflow chip addr (= %d) , setting it to 9", fChipAddr));
 202         fChipAddr=9;
 203       }
 204       fHalfStaveNr = (fData & 0x3800)>>11;
 205       if (fHalfStaveNr>5 || fRawReader->TestBlockAttribute(fHalfStaveNr)) {
 206         fRawReader->AddMajorErrorLog(kStaveNumberErr,Form("Half stave number error %d - set to 5",fHalfStaveNr));
 207         AliWarning(Form("half stave number error(=%d) , setting it to 5", fHalfStaveNr));
 208         fHalfStaveNr=5;
 209       }
 210       // translate  ("online") ddl, hs, chip nr  to  ("offline") module id :
 211       fModuleID = fgkDDLModuleMap[ddlID][fHalfStaveNr*2+fChipAddr/5];
 212       fOffset = 32 * (fChipAddr % 5);
 213     }
 214     else if ((fData & 0xC000) == 0x0000) {    // trailer
 215       UShort_t hitCount = fData & 0x1FFF;
 216       if (hitCount != fHitCount){
 217         fRawReader->AddMajorErrorLog(kNumbHitsErr,Form("Number of hits %d instead of %d",hitCount,fHitCount));
 218         AliWarning(Form("wrong number of hits: %d != %d", fHitCount, hitCount));
 219       }
 220     }
 221     else if ((fData & 0xC000) == 0x8000) {    // pixel hit
 222       fHitCount++;
 223       fCol = (fData & 0x001F);
 224       fRow = (fData >> 5) & 0x00FF;
 225
 226       // translate  ("online") chipcol, chiprow  to  ("offline") col (coord1), row (coord2):
 227       // This will change, waiting for new geometry!!!
 228       fCoord1 = fCol;
 229       //      if      (fModuleID < 80 && ddlID < 10) fCoord1=31-fCoord1;
 230       //      else if (fModuleID >=80 && ddlID >=10) fCoord1=31-fCoord1;
 231       fCoord1 += fOffset;
 232       //      if (ddlID>=10) fCoord1=159-fCoord1;
 233       fCoord2 = fRow;
 234       //      if (fModuleID<80) fCoord2=255-fCoord2;
 235
 236       return kTRUE;
 237     }
 238     else {                                    // fill word
 239       if ((fData & 0xC000) != 0xC000) {
 240         fRawReader->AddMajorErrorLog(kWrongWordErr,"Wrong fill word");
 241         AliWarning("wrong fill word!");
 242       }
 243     }
 244
 245   }
 246
 247   return kFALSE;
 248 }
 249 Bool_t AliITSRawStreamSPD::GetHalfStavePresent(UInt_t hs) {
 250   // Reads the half stave present status from the block attributes
 251   Int_t ddlID = fRawReader->GetDDLID();
 252   if (ddlID==-1) {
 253     fRawReader->AddMinorErrorLog(kHalfStaveStatusErr,"DDL ID = -1. Cannot read block attributes.");
 254     AliWarning("DDL ID = -1. Cannot read block attributes. Return kFALSE.");
 255     return kFALSE;
 256   }
 257   else {
 258     if (hs>=6) {
 259       fRawReader->AddMinorErrorLog(kHalfStaveStatusErr,Form( "HS >= 6 requested (%d). Return kFALSE.",hs));
 260       AliWarning(Form("HS >= 6 requested (%d). Return kFALSE.",hs));
 261       return kFALSE;
 262     }
 263     UChar_t attr = fRawReader->GetBlockAttributes();
 264     if (((attr>>hs) & 0x01) == 0x01) { // bit set means not present
 265       return kFALSE;
 266     }
 267     else {
 268       return kTRUE;
 269     }
 270   }
 271 }
 272
 273 Bool_t AliITSRawStreamSPD::GetHhalfStaveScanned(UInt_t hs) const {
 274   if (hs<6) return (Bool_t)((fCalHeadWord[0]>>(6+hs)) & (0x00000001));
 275   else return kFALSE;
 276 }
 277 Bool_t AliITSRawStreamSPD::GetHchipPresent(UInt_t hs, UInt_t chip) const {
 278   if (hs<6 && chip<10) return ((( fCalHeadWord[hs/3+3]>>((hs%3)*10+chip)) & 0x00000001) == 1);
 279   else return kFALSE;
 280 }
 281 UInt_t AliITSRawStreamSPD::GetHdacHigh(UInt_t hs) const {
 282   if (hs<6) return (fCalHeadWord[hs/2+7]>>(24-16*(hs%2)) & 0x000000ff);
 283   else return 0;
 284 }
 285 UInt_t AliITSRawStreamSPD::GetHdacLow(UInt_t hs) const {
 286   if (hs<6) return (fCalHeadWord[hs/2+7]>>(16-16*(hs%2)) & 0x000000ff);
 287   else return 0;
 288 }
 289 UInt_t AliITSRawStreamSPD::GetHTPAmp(UInt_t hs) const {
 290   if (hs<6) return fCalHeadWord[hs+10];
 291   else return 0;
 292 }