PHOS/AliPHOSRawDigiProducer.cxx

   1 /**************************************************************************
   2  * Copyright(c) 2007, 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 //This class produces PHOS digits of one event
  19 //using AliPHOSRawDecoder.
  20 //
  21 //   For example:
  22 //   TClonesArray *digits = new TClonesArray("AliPHOSDigit",100);
  23 //   AliRawReader* rawReader = new AliRawReaderDate("2006run2211.raw");
  24 //   AliPHOSRawDecoder dc(rawReader);
  25 //   while (rawReader->NextEvent()) {
  26 //     AliPHOSRawDigiProducer producer;
  27 //     producer.MakeDigits(digits,&dc);
  28 //   }
  29
  30 // Author: Boris Polichtchouk
  31
  32 // --- ROOT system ---
  33 #include "TClonesArray.h"
  34
  35 // --- AliRoot header files ---
  36 #include "AliPHOSRawDigiProducer.h"
  37 #include "AliPHOSRawDecoder.h"
  38 #include "AliPHOSGeometry.h"
  39 #include "AliPHOSDigit.h"
  40
  41 ClassImp(AliPHOSRawDigiProducer)
  42
  43 //--------------------------------------------------------------------------------------
  44 void AliPHOSRawDigiProducer::MakeDigits(TClonesArray *digits, AliPHOSRawDecoder* decoder)
  45 {
  46   //Makes the job.
  47   //TClonesArray *digits and raw data decoder should be provided by calling function.
  48
  49   digits->Clear();
  50
  51   AliPHOSGeometry* geo = AliPHOSGeometry::GetInstance();
  52   if(!geo) geo = AliPHOSGeometry::GetInstance("IHEP");
  53
  54   Int_t    iDigit   = 0 ;
  55   Double_t time     = 0. ;
  56   Int_t    iOldDigit;
  57   Bool_t   seen,lowGainFlag;
  58   Int_t relId[4], absId =0;
  59
  60   while (decoder->NextDigit()) {
  61
  62     lowGainFlag = decoder->IsLowGain();
  63     time = decoder->GetTime();
  64
  65     relId[0] = decoder->GetModule();
  66     relId[1] = 0;
  67     relId[2] = decoder->GetRow();
  68     relId[3] = decoder->GetColumn();
  69     geo->RelToAbsNumbering(relId, absId);
  70
  71     // Add low gain digit only
  72     //if the high gain digit does not exist in the digits array
  73
  74     seen = kFALSE;
  75
  76     if(lowGainFlag) {
  77       for (iOldDigit=iDigit-1; iOldDigit>=0; iOldDigit--) {
  78         if ((dynamic_cast<AliPHOSDigit*>(digits->At(iOldDigit)))->GetId() == absId) {
  79           seen = kTRUE;
  80           break;
  81         }
  82       }
  83       if (!seen) {
  84         new((*digits)[iDigit]) AliPHOSDigit(-1,absId,(Float_t)decoder->GetEnergy(),time);
  85         iDigit++;
  86       }
  87     }
  88
  89     // Add high gain digit only if it is not saturated;
  90     // replace low gain digit by a high gain one
  91     else {
  92       if (decoder->GetEnergy() >= 1023) continue;
  93       for (iOldDigit=iDigit-1; iOldDigit>=0; iOldDigit--) {
  94         if ((dynamic_cast<AliPHOSDigit*>(digits->At(iOldDigit)))->GetId() == absId) {
  95           digits->RemoveAt(iOldDigit);
  96           new((*digits)[iOldDigit]) AliPHOSDigit(-1,absId,(Float_t)decoder->GetEnergy(),time);
  97           seen = kTRUE;
  98           break;
  99         }
 100       }
 101       if (!seen) {
 102         new((*digits)[iDigit]) AliPHOSDigit(-1,absId,(Float_t)decoder->GetEnergy(),time);
 103         iDigit++;
 104       }
 105     }
 106
 107   }
 108
 109   digits->Compress();
 110   digits->Sort();
 111 }