MUON/AliMUONDigitMaker.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 /// \class AliMUONDigitMaker
  17 /// MUON Digit maker from rawdata.
  18 ///
  19 /// Raw2Digits:
  20 /// Using real mapping  for tracker
  21 /// Indranil Das (Adapted for runloader: Ch. Finck) july 05
  22 ///
  23 /// Implemented non-constant buspatch numbers for tracking
  24 /// with correct DDL id.
  25 /// (Ch. Finck, dec 05)
  26 ///
  27 /// Add reader for scaler trigger events
  28 /// Use memcpy instead of assignment elt by elt
  29 /// (Ch. Finck, Jan 06)
  30 ///
  31 /// Using new interface with AliMUONRawStreamTracker(Trigger)
  32 /// (New interface of AliMUONRawReader class)
  33 /// (further details could be found in Alice-note)
  34 /// (Ch. Finck, March 06)
  35 ///
  36 /// Add (S)Digit maker tracker (for free)
  37 /// and for trigger. Create trigger inverse mapping.
  38 ///
  39 /// \author Ch. Finck, oct 06
  40
  41 #include "AliMUONDigitMaker.h"
  42
  43 #include "AliLog.h"
  44 #include "AliMUONDDLTrigger.h"
  45 #include "AliMUONDarcHeader.h"
  46 #include "AliMUONVDigit.h"
  47 #include "AliMUONVDigitStore.h"
  48 #include "AliMUONGlobalTrigger.h"
  49 #include "AliMUONLocalStruct.h"
  50 #include "AliMUONLocalTrigger.h"
  51 #include "AliMUONRawStreamTracker.h"
  52 #include "AliMUONRawStreamTrigger.h"
  53 #include "AliMUONRegHeader.h"
  54 #include "AliMUONTriggerCircuit.h"
  55 #include "AliMpTriggerCrate.h"
  56 #include "AliMpLocalBoard.h"
  57 #include "AliMUONVTriggerStore.h"
  58 #include "AliMpCathodType.h"
  59 #include "AliMpDDLStore.h"
  60 #include "AliMpDEManager.h"
  61 #include "AliMpPad.h"
  62 #include "AliMpSegmentation.h"
  63 #include "AliMpVSegmentation.h"
  64 #include "AliRawReader.h"
  65 #include <TArrayS.h>
  66
  67 /// \cond CLASSIMP
  68 ClassImp(AliMUONDigitMaker) // Class implementation in ROOT context
  69 /// \endcond
  70
  71 //__________________________________________________________________________
  72 AliMUONDigitMaker::AliMUONDigitMaker(Bool_t enableErrorLogger)
  73   : TObject(),
  74     fScalerEvent(kFALSE),
  75     fMakeTriggerDigits(kFALSE),
  76     fRawStreamTracker(new AliMUONRawStreamTracker()),
  77     fRawStreamTrigger(new AliMUONRawStreamTrigger()),
  78     fTrackerTimer(),
  79     fTriggerTimer(),
  80     fMappingTimer(),
  81     fDigitStore(0x0),
  82     fTriggerStore(0x0)
  83 {
  84   /// ctor
  85
  86   AliDebug(1,"");
  87
  88   // Standard Constructor
  89   if (enableErrorLogger) {
  90     fRawStreamTracker->EnabbleErrorLogger();
  91     fRawStreamTrigger->EnabbleErrorLogger();
  92   }
  93
  94   fTrackerTimer.Start(kTRUE); fTrackerTimer.Stop();
  95   fTriggerTimer.Start(kTRUE); fTriggerTimer.Stop();
  96   fMappingTimer.Start(kTRUE); fMappingTimer.Stop();
  97
  98   SetMakeTriggerDigits();
  99
 100 }
 101
 102 //__________________________________________________________________________
 103 AliMUONDigitMaker::~AliMUONDigitMaker()
 104 {
 105   /// clean up
 106   /// and time processing measure
 107
 108   delete fRawStreamTracker;
 109   delete fRawStreamTrigger;
 110
 111   AliDebug(1, Form("Execution time for MUON tracker : R:%.2fs C:%.2fs",
 112                fTrackerTimer.RealTime(),fTrackerTimer.CpuTime()));
 113   AliDebug(1, Form("   Execution time for MUON tracker (mapping calls part) "
 114                ": R:%.2fs C:%.2fs",
 115                fMappingTimer.RealTime(),fMappingTimer.CpuTime()));
 116   AliDebug(1, Form("Execution time for MUON trigger : R:%.2fs C:%.2fs",
 117                fTriggerTimer.RealTime(),fTriggerTimer.CpuTime()));
 118
 119 }
 120
 121 //____________________________________________________________________
 122 Int_t AliMUONDigitMaker::Raw2Digits(AliRawReader* rawReader,
 123                                     AliMUONVDigitStore* digitStore,
 124                                     AliMUONVTriggerStore* triggerStore)
 125 {
 126   /// Main method to creates digit
 127   /// for tracker
 128   /// and trigger
 129
 130   AliDebug(1,Form("rawReader=%p digitStore=%p triggerStore=%p",
 131                   rawReader,digitStore,triggerStore));
 132
 133   fDigitStore = digitStore;
 134   fTriggerStore = triggerStore;
 135
 136   if (!fDigitStore && !fTriggerStore)
 137   {
 138     AliError("No digit or trigger store given. Nothing to do...");
 139     return kFALSE;
 140   }
 141
 142   if ( fDigitStore )
 143   {
 144     fDigitStore->Clear(); // insure we start with an empty container
 145     ReadTrackerDDL(rawReader);
 146   }
 147
 148   if ( fTriggerStore || fMakeTriggerDigits )
 149   {
 150     if ( fTriggerStore ) fTriggerStore->Clear();
 151     if ( fMakeTriggerDigits && !fDigitStore )
 152     {
 153       AliError("Asking for trigger digits but digitStore is null");
 154     }
 155     else
 156     {
 157       ReadTriggerDDL(rawReader);
 158     }
 159   }
 160
 161   return kTRUE;
 162 }
 163
 164 //____________________________________________________________________
 165 Int_t AliMUONDigitMaker::ReadTrackerDDL(AliRawReader* rawReader)
 166 {
 167   /// Reading tracker DDL
 168   /// filling the fDigitStore container, which must not be null
 169
 170   AliDebug(1,"");
 171
 172   fTrackerTimer.Start(kFALSE);
 173
 174   // elex info
 175   Int_t    buspatchId;
 176   UChar_t  channelId;
 177   UShort_t manuId;
 178   UShort_t charge;
 179
 180   fRawStreamTracker->SetReader(rawReader);
 181   fRawStreamTracker->First();
 182
 183   while ( fRawStreamTracker->Next(buspatchId,manuId,channelId,charge) )
 184   {
 185     // getting DE from buspatch
 186     Int_t detElemId = AliMpDDLStore::Instance()->GetDEfromBus(buspatchId);
 187
 188     const AliMpVSegmentation* seg
 189       = AliMpSegmentation::Instance()->GetMpSegmentationByElectronics(detElemId,
 190                                                                       manuId);
 191
 192     AliMp::CathodType cathodeType = AliMpDEManager::GetCathod(detElemId,
 193                                                               seg->PlaneType());
 194
 195     AliMpPad pad = seg->PadByLocation(AliMpIntPair(manuId,channelId),kFALSE);
 196
 197     if (!pad.IsValid())
 198     {
 199       AliError(Form("No pad for detElemId: %d, manuId: %d, channelId: %d",
 200                     detElemId, manuId, channelId));
 201       continue;
 202     }
 203
 204     AliMUONVDigit* digit = fDigitStore->Add(detElemId,manuId,channelId,cathodeType,
 205                                             AliMUONVDigitStore::kDeny);
 206     if (!digit)
 207     {
 208       AliError(Form("Digit DE %04d Manu %04d Channel %02d could not be added",
 209                     detElemId, manuId, channelId));
 210       continue;
 211     }
 212
 213     digit->SetPadXY(pad.GetIndices().GetFirst(),
 214                    pad.GetIndices().GetSecond());
 215
 216           digit->SetADC(charge);
 217
 218   }
 219
 220   fTrackerTimer.Stop();
 221
 222   return kTRUE;
 223 }
 224
 225 //____________________________________________________________________
 226 Int_t AliMUONDigitMaker::ReadTriggerDDL(AliRawReader* rawReader)
 227 {
 228   /// reading tracker DDL
 229   /// filling the fTriggerStore container, which must not be null
 230
 231   AliDebug(1,"");
 232
 233   AliMUONDDLTrigger*       ddlTrigger      = 0x0;
 234   AliMUONDarcHeader*       darcHeader      = 0x0;
 235   AliMUONRegHeader*        regHeader       = 0x0;
 236   AliMUONLocalStruct*      localStruct     = 0x0;
 237
 238   Int_t loCircuit;
 239
 240   fTriggerTimer.Start(kFALSE);
 241
 242   fRawStreamTrigger->SetReader(rawReader);
 243
 244   while (fRawStreamTrigger->NextDDL())
 245   {
 246     ddlTrigger = fRawStreamTrigger->GetDDLTrigger();
 247     darcHeader = ddlTrigger->GetDarcHeader();
 248
 249     // fill global trigger information
 250     if (fTriggerStore)
 251     {
 252       if (darcHeader->GetGlobalFlag())
 253       {
 254           AliMUONGlobalTrigger globalTrigger;
 255           globalTrigger.SetFromGlobalResponse(darcHeader->GetGlobalOutput());
 256           fTriggerStore->SetGlobal(globalTrigger);
 257       }
 258     }
 259
 260     Int_t nReg = darcHeader->GetRegHeaderEntries();
 261
 262     for(Int_t iReg = 0; iReg < nReg ;iReg++)
 263     {   //reg loop
 264
 265
 266       // crate info
 267       AliMpTriggerCrate* crate = AliMpDDLStore::Instance()->
 268                                 GetTriggerCrate(fRawStreamTrigger->GetDDL(), iReg);
 269
 270       if (!crate)
 271         AliWarning(Form("Missing crate number %d in DDL %d\n", iReg, fRawStreamTrigger->GetDDL()));
 272
 273
 274       regHeader =  darcHeader->GetRegHeaderEntry(iReg);
 275
 276       Int_t nLocal = regHeader->GetLocalEntries();
 277       for(Int_t iLocal = 0; iLocal < nLocal; iLocal++)
 278       {
 279
 280         localStruct = regHeader->GetLocalEntry(iLocal);
 281
 282         // if card exist
 283         if (localStruct) {
 284
 285           loCircuit = crate->GetLocalBoardId(localStruct->GetId());
 286
 287           if ( !loCircuit ) continue; // empty slot
 288
 289           AliMpLocalBoard* localBoard = AliMpDDLStore::Instance()->GetLocalBoard(loCircuit, false);
 290
 291           // skip copy cards
 292           if( !localBoard->IsNotified())
 293              continue;
 294
 295           if (fTriggerStore)
 296           {
 297             // fill local trigger
 298             AliMUONLocalTrigger localTrigger;
 299             localTrigger.SetLocalStruct(loCircuit, *localStruct);
 300             fTriggerStore->Add(localTrigger);
 301           }
 302
 303           if ( fMakeTriggerDigits )
 304           {
 305             //FIXEME should find something better than a TArray
 306             TArrayS xyPattern[2];
 307
 308             localStruct->GetXPattern(xyPattern[0]);
 309             localStruct->GetYPattern(xyPattern[1]);
 310
 311             TriggerDigits(loCircuit, xyPattern, *fDigitStore);
 312           }
 313         } // if triggerY
 314       } // iLocal
 315     } // iReg
 316   } // NextDDL
 317
 318   fTriggerTimer.Stop();
 319
 320   return kTRUE;
 321
 322 }
 323
 324 //____________________________________________________________________
 325 Int_t AliMUONDigitMaker::TriggerDigits(Int_t nBoard,
 326                                        TArrayS* xyPattern,
 327                                        AliMUONVDigitStore& digitStore) const
 328 {
 329   /// make digits for trigger from pattern, and add them to digitStore
 330
 331   Int_t detElemId;
 332
 333   AliMpLocalBoard* localBoard = AliMpDDLStore::Instance()->GetLocalBoard(nBoard);
 334
 335   Int_t n,b;
 336
 337   // loop over x1-4 and y1-4
 338   for (Int_t iChamber = 0; iChamber < 4; ++iChamber)
 339   {
 340     for (Int_t iCath = 0; iCath < 2; ++iCath)
 341     {
 342       Int_t pattern = (Int_t)xyPattern[iCath].At(iChamber);
 343       if (!pattern) continue;
 344
 345       // get detElemId
 346       detElemId = AliMpDDLStore::Instance()->GetDEfromLocalBoard(nBoard, iChamber);
 347
 348         const AliMpVSegmentation* seg
 349           = AliMpSegmentation::Instance()
 350           ->GetMpSegmentation(detElemId, AliMp::GetCathodType(iCath));
 351
 352         // loop over the 16 bits of pattern
 353         for (Int_t ibitxy = 0; ibitxy < 16; ++ibitxy)
 354         {
 355           if ((pattern >> ibitxy) & 0x1)
 356           {
 357             // not quite sure about this
 358             Int_t offset = 0;
 359             if (iCath && localBoard->GetSwitch(6)) offset = -8;
 360
 361             AliMpPad pad = seg->PadByLocation(AliMpIntPair(nBoard,ibitxy+offset),kTRUE);
 362
 363             if (!pad.IsValid())
 364             {
 365               AliWarning(Form("No pad for detElemId: %d, nboard %d, ibitxy: %d\n",
 366                               detElemId, nBoard, ibitxy));
 367               continue;
 368             }
 369
 370             n = pad.GetLocation(0).GetFirst(); // always take first location so that digits are not inserted several times
 371             b = pad.GetLocation(0).GetSecond();
 372
 373             AliDebug(1,Form("Using localBoard %d ixy %d instead of %d,%d",
 374                             n,b,nBoard,ibitxy));
 375
 376             AliMUONVDigit* digit = digitStore.Add(detElemId,n,b,iCath,AliMUONVDigitStore::kDeny);
 377
 378             if (!digit)
 379             {
 380                 AliDebug(1, Form("Digit DE %04d LocalBoard %03d ibitxy %02d cath %d already in store",
 381                                  detElemId,nBoard,ibitxy,iCath));
 382                 continue;
 383             }
 384
 385             Int_t padX = pad.GetIndices().GetFirst();
 386             Int_t padY = pad.GetIndices().GetSecond();
 387
 388             // fill digit
 389             digit->SetPadXY(padX,padY);
 390             digit->SetCharge(1.);
 391           }// xyPattern
 392         }// ibitxy
 393     }// cath
 394   } // ichamber
 395
 396   return kTRUE;
 397 }