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 // $Id$
  17
  18 //-----------------------------------------------------------------------------
  19 /// \class AliMUONDigitMaker
  20 /// MUON Digit maker from rawdata.
  21 ///
  22 /// Raw2Digits:
  23 /// Using real mapping  for tracker
  24 /// Indranil Das (Adapted for runloader: Ch. Finck) july 05
  25 ///
  26 /// Implemented non-constant buspatch numbers for tracking
  27 /// with correct DDL id.
  28 /// (Ch. Finck, dec 05)
  29 ///
  30 /// Add reader for scaler trigger events
  31 /// Use memcpy instead of assignment elt by elt
  32 /// (Ch. Finck, Jan 06)
  33 ///
  34 /// Using new interface with AliMUONRawStreamTracker(Trigger)
  35 /// (New interface of AliMUONRawReader class)
  36 /// (further details could be found in Alice-note)
  37 /// (Ch. Finck, March 06)
  38 ///
  39 /// Add (S)Digit maker tracker (for free)
  40 /// and for trigger. Create trigger inverse mapping.
  41 ///
  42 /// \author Ch. Finck, oct 06
  43 //-----------------------------------------------------------------------------
  44
  45 #include "AliMUONDigitMaker.h"
  46
  47 #include "AliMUONDDLTrigger.h"
  48 #include "AliMUONDarcHeader.h"
  49 #include "AliMUONVDigit.h"
  50 #include "AliMUONVDigitStore.h"
  51 #include "AliMUONGlobalTrigger.h"
  52 #include "AliMUONLocalStruct.h"
  53 #include "AliMUONLocalTrigger.h"
  54 #include "AliMUONLogger.h"
  55 #include "AliMUONRawStreamTrackerHP.h"
  56 #include "AliMUONRawStreamTriggerHP.h"
  57 #include "AliMUONRegHeader.h"
  58 #include "AliMUONVTriggerStore.h"
  59 #include "AliMpCDB.h"
  60 #include "AliMpDetElement.h"
  61 #include "AliMpTriggerCrate.h"
  62 #include "AliMpLocalBoard.h"
  63 #include "AliMpCathodType.h"
  64 #include "AliMpDDLStore.h"
  65 #include "AliMpDEManager.h"
  66 #include "AliMpPad.h"
  67 #include "AliMpSegmentation.h"
  68 #include "AliMpVSegmentation.h"
  69 #include "AliCodeTimer.h"
  70 #include "AliLog.h"
  71 #include "AliRawReader.h"
  72 #include <TArrayS.h>
  73
  74 /// \cond CLASSIMP
  75 ClassImp(AliMUONDigitMaker) // Class implementation in ROOT context
  76 /// \endcond
  77
  78 //__________________________________________________________________________
  79 AliMUONDigitMaker::AliMUONDigitMaker(Bool_t enableErrorLogger, Bool_t a, Bool_t b) :
  80 TObject(),
  81 fScalerEvent(kFALSE),
  82 fMakeTriggerDigits(kFALSE),
  83 fRawStreamTracker(new AliMUONRawStreamTrackerHP),
  84 fRawStreamTrigger(new AliMUONRawStreamTriggerHP),
  85 fDigitStore(0x0),
  86 fTriggerStore(0x0),
  87 fLogger(new AliMUONLogger(10000)){
  88   /// ctor
  89
  90   if  ( !a || !b ) AliFatal("no longer supported");
  91
  92   AliDebug(1,"");
  93
  94   // Standard Constructor
  95   if (enableErrorLogger)
  96   {
  97     fRawStreamTracker->EnabbleErrorLogger();
  98     fRawStreamTrigger->EnabbleErrorLogger();
  99   }
 100   else
 101   {
 102     fRawStreamTracker->DisableWarnings();
 103   }
 104
 105   SetMakeTriggerDigits();
 106
 107   // Load mapping
 108   if ( ! AliMpCDB::LoadDDLStore() ) {
 109     AliFatal("Could not access mapping from OCDB !");
 110   }
 111 }
 112
 113 //__________________________________________________________________________
 114 AliMUONDigitMaker::AliMUONDigitMaker(Bool_t enableErrorLogger) :
 115 TObject(),
 116     fScalerEvent(kFALSE),
 117     fMakeTriggerDigits(kFALSE),
 118     fRawStreamTracker(new AliMUONRawStreamTrackerHP),
 119     fRawStreamTrigger(new AliMUONRawStreamTriggerHP),
 120     fDigitStore(0x0),
 121     fTriggerStore(0x0),
 122   fLogger(new AliMUONLogger(10000))
 123 {
 124   /// ctor
 125
 126   AliDebug(1,"");
 127
 128   // Standard Constructor
 129   if (enableErrorLogger)
 130   {
 131     fRawStreamTracker->EnabbleErrorLogger();
 132     fRawStreamTrigger->EnabbleErrorLogger();
 133   }
 134   else
 135   {
 136     fRawStreamTracker->DisableWarnings();
 137   }
 138
 139   SetMakeTriggerDigits();
 140
 141   // Load mapping
 142   if ( ! AliMpCDB::LoadDDLStore() ) {
 143     AliFatal("Could not access mapping from OCDB !");
 144   }
 145 }
 146
 147 //__________________________________________________________________________
 148 AliMUONDigitMaker::~AliMUONDigitMaker()
 149 {
 150   /// clean up
 151   /// and time processing measure
 152
 153   delete fRawStreamTracker;
 154   delete fRawStreamTrigger;
 155   delete fLogger;
 156 }
 157
 158 //____________________________________________________________________
 159 void
 160 AliMUONDigitMaker::Print(Option_t*) const
 161 {
 162   /// Printout
 163
 164   cout << "RawStreamerTracker class=" << fRawStreamTracker->ClassName()
 165        << " MakeTriggerDigits=" << fMakeTriggerDigits
 166        << " ScalerEvent=" << fScalerEvent
 167        << " DigitStore=" << fDigitStore
 168        << " TriggerStore=" << fTriggerStore << endl;
 169
 170   if ( fLogger ) fLogger->Print();
 171 }
 172
 173 //____________________________________________________________________
 174 Int_t
 175 AliMUONDigitMaker::Raw2Digits(AliRawReader* rawReader,
 176                                         AliMUONVDigitStore* digitStore,
 177                                         AliMUONVTriggerStore* triggerStore)
 178 {
 179   /// Main method to creates digit
 180   /// for tracker
 181   /// and trigger
 182
 183   AliDebug(1,Form("rawReader=%p digitStore=%p triggerStore=%p",
 184                   rawReader,digitStore,triggerStore));
 185
 186   fDigitStore = digitStore;
 187   fTriggerStore = triggerStore;
 188
 189   if (!fDigitStore && !fTriggerStore)
 190   {
 191     fLogger->Log("No digit or trigger store given. Nothing to do...");
 192     return kTriggerBAD & kTrackerBAD;
 193   }
 194
 195   Int_t tracker(kOK);
 196   Int_t trigger(kOK);
 197
 198   if ( fDigitStore )
 199   {
 200     fDigitStore->Clear(); // insure we start with an empty container
 201     tracker = ReadTrackerDDL(rawReader);
 202   }
 203
 204   if ( fTriggerStore || fMakeTriggerDigits )
 205   {
 206     if ( fTriggerStore ) fTriggerStore->Clear();
 207     if ( fMakeTriggerDigits && !fDigitStore )
 208     {
 209       fLogger->Log("Asking for trigger digits but digitStore is null");
 210     }
 211     else
 212     {
 213       trigger = ReadTriggerDDL(rawReader);
 214     }
 215   }
 216
 217   return tracker | trigger;
 218 }
 219
 220 //____________________________________________________________________
 221 Int_t
 222 AliMUONDigitMaker::ReadTrackerDDL(AliRawReader* rawReader)
 223 {
 224   /// Reading tracker DDL
 225   /// filling the fDigitStore container, which must not be null
 226
 227   AliDebug(1,"");
 228
 229   AliCodeTimerAuto("",0);
 230
 231   // elex info
 232   Int_t    buspatchId;
 233   UChar_t  channelId;
 234   UShort_t manuId;
 235   UShort_t charge;
 236
 237   fRawStreamTracker->SetReader(rawReader);
 238   fRawStreamTracker->First();
 239
 240   while ( fRawStreamTracker->Next(buspatchId,manuId,channelId,charge,kTRUE) )
 241   {
 242     // getting DE from buspatch
 243     Int_t detElemId = AliMpDDLStore::Instance()->GetDEfromBus(buspatchId);
 244
 245     AliMpDetElement* de = AliMpDDLStore::Instance()->GetDetElement(detElemId);
 246
 247     if (!de)
 248       {
 249         fLogger->Log(Form("DE %04d does not exist !", detElemId));
 250         continue;
 251       }
 252
 253     if (!de->IsConnectedChannel(manuId,channelId))
 254       {
 255         // non connected pad, do nothing (this is not an error !)
 256         continue;
 257       }
 258
 259     const AliMpVSegmentation* seg
 260       = AliMpSegmentation::Instance()->GetMpSegmentationByElectronics(detElemId,
 261                                                                       manuId);
 262
 263     if (!seg)
 264     {
 265       fLogger->Log(Form("(DE,MANUID)=(%04d,%04d) is not valid",detElemId,manuId));
 266       continue;
 267     }
 268
 269     AliMp::CathodType cathodeType = de->GetCathodType(seg->PlaneType());
 270
 271     AliMpPad pad = seg->PadByLocation(manuId,channelId,kFALSE);
 272
 273     if (!pad.IsValid())
 274     {
 275       fLogger->Log(Form("No pad for detElemId: %d, manuId: %d, channelId: %d",
 276                     detElemId, manuId, channelId));
 277       continue;
 278     }
 279
 280     AliMUONVDigit* digit = fDigitStore->Add(detElemId,manuId,channelId,cathodeType,
 281                                             AliMUONVDigitStore::kDeny);
 282
 283     if (!digit)
 284     {
 285       fLogger->Log(Form("Digit DE %04d Manu %04d Channel %02d could not be added",
 286                     detElemId, manuId, channelId));
 287       continue;
 288     }
 289
 290     digit->SetPadXY(pad.GetIx(),pad.GetIy());
 291
 292     digit->SetADC(charge);
 293
 294   }
 295
 296   if ( fRawStreamTracker->IsErrorMessage() )
 297   {
 298     return kTrackerBAD;
 299   }
 300
 301   return kOK;
 302 }
 303
 304 //____________________________________________________________________
 305 Int_t
 306 AliMUONDigitMaker::ReadTriggerDDL(AliRawReader* rawReader)
 307 {
 308   /// reading tracker DDL like ReadTriggerDDL but with fast decoder interface.
 309   /// filling the fTriggerStore container, which must not be null
 310
 311   const AliMUONRawStreamTriggerHP::AliHeader*          darcHeader  = 0x0;
 312   const AliMUONRawStreamTriggerHP::AliRegionalHeader*  regHeader   = 0x0;
 313   const AliMUONRawStreamTriggerHP::AliLocalStruct*     localStruct = 0x0;
 314
 315   Int_t loCircuit;
 316
 317   fRawStreamTrigger->SetReader(rawReader);
 318
 319   while (fRawStreamTrigger->NextDDL())
 320   {
 321     darcHeader = fRawStreamTrigger->GetHeaders();
 322
 323     // fill global trigger information
 324     if (fTriggerStore)
 325     {
 326       if (darcHeader->GetGlobalFlag())
 327       {
 328           AliMUONGlobalTrigger globalTrigger;
 329           globalTrigger.SetFromGlobalResponse(darcHeader->GetGlobalOutput());
 330           globalTrigger.SetFromGlobalInput(darcHeader->GetGlobalHeader()->fInput);
 331           fTriggerStore->SetGlobal(globalTrigger);
 332       }
 333     }
 334
 335     Int_t nReg = fRawStreamTrigger->GetRegionalHeaderCount();
 336
 337     for(Int_t iReg = 0; iReg < nReg ;iReg++)
 338     {   //reg loop
 339
 340
 341       // crate info
 342       AliMpTriggerCrate* crate = AliMpDDLStore::Instance()->
 343                                 GetTriggerCrate(fRawStreamTrigger->GetDDL(), iReg);
 344
 345       if (!crate) {
 346         fLogger->Log(Form("Missing crate number %d in DDL %d\n", iReg, fRawStreamTrigger->GetDDL()));
 347         continue;
 348       }
 349
 350       regHeader =  fRawStreamTrigger->GetRegionalHeader(iReg);
 351
 352       Int_t nLocal = regHeader->GetLocalStructCount();
 353       for(Int_t iLocal = 0; iLocal < nLocal; iLocal++)
 354       {
 355
 356         localStruct = regHeader->GetLocalStruct(iLocal);
 357
 358         // if card exist
 359         if (localStruct) {
 360
 361           loCircuit = crate->GetLocalBoardId(localStruct->GetId());
 362
 363           if ( !loCircuit ) continue; // empty slot
 364
 365           AliMpLocalBoard* localBoard = AliMpDDLStore::Instance()->GetLocalBoard(loCircuit, kTRUE);
 366
 367           // skip copy cards
 368           if( !localBoard->IsNotified())
 369              continue;
 370
 371           if (fTriggerStore)
 372           {
 373             // fill local trigger
 374             AliMUONLocalTrigger localTrigger;
 375             localTrigger.SetLocalStruct(loCircuit, *localStruct);
 376             fTriggerStore->Add(localTrigger);
 377           }
 378
 379           if ( fMakeTriggerDigits )
 380           {
 381             //FIXEME should find something better than a TArray
 382             TArrayS xyPattern[2];
 383
 384             localStruct->GetXPattern(xyPattern[0]);
 385             localStruct->GetYPattern(xyPattern[1]);
 386
 387             TriggerDigits(loCircuit, xyPattern, *fDigitStore);
 388           }
 389         } // if triggerY
 390       } // iLocal
 391     } // iReg
 392   } // NextDDL
 393
 394   return kOK;
 395 }
 396
 397 //____________________________________________________________________
 398 Int_t AliMUONDigitMaker::TriggerDigits(Int_t nBoard,
 399                                        const TArrayS* xyPattern,
 400                                        AliMUONVDigitStore& digitStore) const
 401 {
 402   /// make digits for trigger from pattern, and add them to digitStore
 403
 404   AliCodeTimerAuto("",0);
 405
 406   Int_t detElemId;
 407
 408   AliMpLocalBoard* localBoard = AliMpDDLStore::Instance()->GetLocalBoard(nBoard);
 409
 410   Int_t n,b;
 411
 412   // loop over x1-4 and y1-4
 413   for (Int_t iChamber = 0; iChamber < 4; ++iChamber)
 414   {
 415     for (Int_t iCath = 0; iCath < 2; ++iCath)
 416     {
 417       Int_t pattern = (Int_t)xyPattern[iCath].At(iChamber);
 418       if (!pattern) continue;
 419
 420       // get detElemId
 421       detElemId = AliMpDDLStore::Instance()->GetDEfromLocalBoard(nBoard, iChamber);
 422
 423         const AliMpVSegmentation* seg
 424           = AliMpSegmentation::Instance()
 425           ->GetMpSegmentation(detElemId, AliMp::GetCathodType(iCath));
 426
 427         // loop over the 16 bits of pattern
 428         for (Int_t ibitxy = 0; ibitxy < 16; ++ibitxy)
 429         {
 430           if ((pattern >> ibitxy) & 0x1)
 431           {
 432             // not quite sure about this
 433             Int_t offset = 0;
 434             if (iCath && localBoard->GetSwitch(AliMpLocalBoard::kZeroAllYLSB)) offset = -8;
 435
 436             AliMpPad pad = seg->PadByLocation(nBoard,ibitxy+offset,kTRUE);
 437
 438             if (!pad.IsValid())
 439             {
 440               fLogger->Log(Form("No pad for detElemId: %d, nboard %d, ibitxy: %d\n",
 441                               detElemId, nBoard, ibitxy));
 442               continue;
 443             }
 444
 445             n = pad.GetLocalBoardId(0); // always take first location so that digits are not inserted several times
 446             b = pad.GetLocalBoardChannel(0);
 447
 448             AliDebug(1,Form("Using localBoard %d ixy %d instead of %d,%d",
 449                             n,b,nBoard,ibitxy));
 450
 451             AliMUONVDigit* digit = digitStore.Add(detElemId,n,b,iCath,AliMUONVDigitStore::kDeny);
 452
 453             if (!digit)
 454             {
 455                 AliDebug(1, Form("Digit DE %04d LocalBoard %03d ibitxy %02d cath %d already in store",
 456                                  detElemId,nBoard,ibitxy,iCath));
 457                 continue;
 458             }
 459
 460             Int_t padX = pad.GetIx();
 461             Int_t padY = pad.GetIy();
 462
 463             // fill digit
 464             digit->SetPadXY(padX,padY);
 465             digit->SetCharge(1.);
 466           }// xyPattern
 467         }// ibitxy
 468     }// cath
 469   } // ichamber
 470
 471   return kTRUE;
 472 }
 473
 474 //______________________________________________________________________________
 475 Bool_t
 476 AliMUONDigitMaker::TriggerToDigitsStore(const AliMUONVTriggerStore& triggerStore,
 477                                         AliMUONVDigitStore& digitStore) const
 478 {
 479   //
 480   /// make (S)Digit for trigger
 481   //
 482
 483   digitStore.Clear();
 484
 485   AliMUONLocalTrigger* locTrg;
 486   TIter next(triggerStore.CreateLocalIterator());
 487
 488   while ( ( locTrg = static_cast<AliMUONLocalTrigger*>(next()) ) )
 489   {
 490     if (locTrg->IsNull()) continue;
 491
 492     TArrayS xyPattern[2];
 493     locTrg->GetXPattern(xyPattern[0]);
 494     locTrg->GetYPattern(xyPattern[1]);
 495
 496     Int_t nBoard = locTrg->LoCircuit();
 497     TriggerDigits(nBoard, xyPattern, digitStore);
 498   }
 499   return kTRUE;
 500 }
 501
 502 //______________________________________________________________________________
 503 void
 504 AliMUONDigitMaker::SetTryRecover(Bool_t flag)
 505 {
 506   /// Instruct the decoder to try to recover corrupted raw data.
 507   /// Only use for specific cases for which you know it will work...
 508   fRawStreamTracker->TryRecover(flag);
 509 }