MUON/MUONRawStreamTrigger.C

   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 //
  17 // Macro for reading tracker raw data
  18 // Ch. Finck, Subatech Febuary
  19 //
  20
  21 #if !defined(__CINT__) || defined(__MAKECINT__)
  22
  23 // RAW includes
  24 #include "AliRawReaderDate.h"
  25 #include "AliRawReaderFile.h"
  26 #include "AliRawReaderRoot.h"
  27
  28 // MUON includes
  29 #include "AliMUONRawStreamTrigger.h"
  30 #include "AliMUONDarcHeader.h"
  31 #include "AliMUONRegHeader.h"
  32 #include "AliMUONLocalStruct.h"
  33 #include "AliMUONDDLTrigger.h"
  34 #include "AliMpTriggerCrate.h"
  35 #include "AliMpDDLStore.h"
  36 #include "AliMpCDB.h"
  37
  38 #include "TStopwatch.h"
  39
  40
  41 #endif
  42
  43 // Macro to read rawdata for trigger
  44 // Ch. Finck, Subatech, April. 2006
  45 // Implement "digits" iterator
  46 // This macro is interface with AliRawReader for RAW
  47 // The different stucture of the patload are readout and stored in TClonesArray
  48 // with AliMUONRawStreamTrigger classe.
  49 // The macro just simpy read again the TClonesArray contents.
  50 // The parameter of each structure could be seen in the container classes
  51 // AliMUONDarcHeader, AliMUONRegHeader, AliMUONLocalStruct.
  52 // The class AliMUONDDLTrigger manages the structure containers.
  53 // The number of structures in the rawdata file could be set.
  54 // The DATE format reading is no more supported please use the MUONTRGda code
  55
  56 void MUONRawStreamTrigger(Int_t maxEvent = 1, Int_t minDDL = 0, Int_t maxDDL = 1, TString fileName = "./")
  57 {
  58
  59    TStopwatch timer;
  60    timer.Start(kTRUE);
  61
  62    AliRawReader* rawReader = 0x0;
  63
  64    if (fileName.EndsWith("/")) {
  65      rawReader = new AliRawReaderFile(fileName);// DDL files
  66    } else if (fileName.EndsWith(".root")) {
  67      rawReader = new AliRawReaderRoot(fileName);
  68    }
  69
  70
  71    // Load mapping
  72      if ( ! AliMpCDB::LoadDDLStore() ) {
  73        printf("Could not access mapping from OCDB !\n");
  74      }
  75
  76    // raw stream
  77    AliMUONRawStreamTrigger* rawStream   = new AliMUONRawStreamTrigger(rawReader);
  78
  79    // set the number ofreg & local that are PRESENT in the rawdata file
  80    // it's NOT the number to be read.
  81    // default wise set to 8, 16 respectively.
  82    //    rawStream->SetMaxReg(2);
  83    //    rawStream->SetMaxLoc(xx);
  84
  85    // containers
  86    AliMUONDDLTrigger*       ddlTrigger  = 0x0;
  87    AliMUONDarcHeader*       darcHeader  = 0x0;
  88    AliMUONRegHeader*        regHeader   = 0x0;
  89    AliMUONLocalStruct*      localStruct = 0x0;
  90
  91
  92    // Loop over events
  93    Int_t iEvent = 0;
  94
  95    while (rawReader->NextEvent()) {
  96
  97      if (iEvent == maxEvent)
  98        break;
  99
 100      printf("Event %d\n",iEvent++);
 101
 102      // read DDL while < 2 DDL
 103      while(rawStream->NextDDL()) {
 104
 105       if (rawStream->GetDDL() < minDDL || rawStream->GetDDL() > maxDDL)
 106          continue;
 107
 108        printf("\niDDL %d\n", rawStream->GetDDL());
 109
 110        ddlTrigger = rawStream->GetDDLTrigger();
 111        darcHeader = ddlTrigger->GetDarcHeader();
 112
 113        printf("Global output %x\n", (Int_t)darcHeader->GetGlobalOutput());
 114
 115        // loop over regional structures
 116        Int_t nReg = darcHeader->GetRegHeaderEntries();
 117        for(Int_t iReg = 0; iReg < nReg ;iReg++){   //REG loop
 118
 119 //       printf("RegionalId %d\n", iReg);
 120
 121          regHeader =  darcHeader->GetRegHeaderEntry(iReg);
 122          //  printf("Reg length %d\n",regHeader->GetHeaderLength());
 123
 124          // crate info
 125          AliMpTriggerCrate* crate = AliMpDDLStore::Instance()->
 126                                     GetTriggerCrate(rawStream->GetDDL(), iReg);
 127
 128          // loop over local structures
 129          Int_t nLocal = regHeader->GetLocalEntries();
 130          for(Int_t iLocal = 0; iLocal < nLocal; iLocal++) {
 131
 132            localStruct = regHeader->GetLocalEntry(iLocal);
 133
 134            Int_t iLocCard = crate->GetLocalBoardId(localStruct->GetId());
 135
 136            if ( !iLocCard ) continue; // empty slot
 137
 138            // check if trigger
 139            if (localStruct->GetTriggerX()
 140                || localStruct->GetTriggerY()) { // no empty data
 141
 142              printf("LocalId %d\n", localStruct->GetId());
 143
 144              Int_t loStripX  = (Int_t)localStruct->GetXPos();
 145              Int_t loStripY  = (Int_t)localStruct->GetYPos();
 146              Int_t loDev     = (Int_t)localStruct->GetXDev();
 147
 148              printf("iLocCard: %d, XPos: %d, YPos: %d Dev: %d\n", iLocCard, loStripX, loStripY, loDev);
 149
 150            }
 151          } // iLocal
 152        } // iReg
 153      } // NextDDL
 154    }// NextEvent
 155
 156    delete rawReader;
 157    delete rawStream;
 158
 159    timer.Print();
 160
 161 }
 162
 163
 164 void MUONRawStreamTriggerSimple(Int_t maxEvent = 1, TString fileName = "./")
 165 {
 166   /// Reads the raw data in fileName, using a simplified interface (iterator
 167   /// over local structure response).
 168
 169   TStopwatch timer;
 170   timer.Start(kTRUE);
 171
 172    AliRawReader* rawReader = 0x0;
 173
 174    if (fileName.EndsWith("/")) {
 175      rawReader = new AliRawReaderFile(fileName);// DDL files
 176    } else if (fileName.EndsWith(".root")) {
 177      rawReader = new AliRawReaderRoot(fileName);
 178    }
 179
 180    // raw stream
 181    AliMUONRawStreamTrigger* rawStream   = new AliMUONRawStreamTrigger(rawReader);
 182
 183    // set the number of reg & local that are PRESENT in the rawdata file
 184    // it's NOT the number to be read.
 185    // default wise set to 8, 16 respectively.
 186    //    rawStream->SetMaxReg(2);
 187    //    rawStream->SetMaxLoc(xx);
 188
 189    UChar_t id;
 190    UChar_t dec;
 191    Bool_t trigY;
 192    UChar_t yPos;
 193    UChar_t sXDev;
 194    UChar_t xDev;
 195    UChar_t xPos;
 196
 197    Bool_t triggerX;
 198    Bool_t triggerY;
 199
 200
 201    TArrayS xPattern;
 202    TArrayS yPattern;
 203
 204    // Loop over events
 205    Int_t iEvent = 0;
 206
 207    while (rawReader->NextEvent()) {
 208
 209      if (iEvent == maxEvent)
 210          break;
 211
 212      printf("Event %d\n",iEvent++);
 213
 214      rawStream->First();
 215
 216      // read while there are digits
 217      while( rawStream->Next(id, dec, trigY, yPos, sXDev, xDev, xPos,
 218                             triggerX, triggerY, xPattern, yPattern) )
 219      {
 220        if ( triggerX || triggerY )  // no empty data
 221            printf("iLocCard: %d, XPos: %d, YPos: %d Dev: %d\n", id, xPos, yPos, xDev);
 222
 223      }// Next
 224
 225    }// NextEvent
 226
 227    delete rawReader;
 228    delete rawStream;
 229
 230    timer.Print();
 231
 232 }