Moved calibration and cleaning to RawDigiProducer
[u/mrichter/AliRoot.git] / PHOS / AliPHOSRawDigiProducer.cxx
CommitLineData
7bc140a5 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 ---
88fb5e50 33#include "TClonesArray.h"
34
7bc140a5 35// --- AliRoot header files ---
88fb5e50 36#include "AliPHOSRawDigiProducer.h"
37#include "AliPHOSRawDecoder.h"
38#include "AliPHOSGeometry.h"
39#include "AliPHOSDigit.h"
40
41ClassImp(AliPHOSRawDigiProducer)
42
7bc140a5 43//--------------------------------------------------------------------------------------
88fb5e50 44void AliPHOSRawDigiProducer::MakeDigits(TClonesArray *digits, AliPHOSRawDecoder* decoder)
45{
7bc140a5 46 //Makes the job.
47 //TClonesArray *digits and raw data decoder should be provided by calling function.
48
88fb5e50 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
34827459 62 if (decoder->GetEnergy() <= 0.)
63 continue;
64
65 if(decoder->IsOverflow())
66 continue ;
1267d56d 67
88fb5e50 68 lowGainFlag = decoder->IsLowGain();
69 time = decoder->GetTime();
70
71 relId[0] = decoder->GetModule();
72 relId[1] = 0;
73 relId[2] = decoder->GetRow();
74 relId[3] = decoder->GetColumn();
75 geo->RelToAbsNumbering(relId, absId);
76
77 // Add low gain digit only
78 //if the high gain digit does not exist in the digits array
34827459 79 //or has negative energy (overflow)
88fb5e50 80
81 seen = kFALSE;
82
83 if(lowGainFlag) {
84 for (iOldDigit=iDigit-1; iOldDigit>=0; iOldDigit--) {
34827459 85 AliPHOSDigit * dig = dynamic_cast<AliPHOSDigit*>(digits->At(iOldDigit)) ;
86 if (dig->GetId() == absId) {
88fb5e50 87 seen = kTRUE;
34827459 88// //if High Gain overflowed - replace energy
89// if(dig->GetEnergy()<0)
90// dig->SetEnergy((Float_t)decoder->GetEnergy()) ;
88fb5e50 91 break;
92 }
93 }
94 if (!seen) {
95 new((*digits)[iDigit]) AliPHOSDigit(-1,absId,(Float_t)decoder->GetEnergy(),time);
96 iDigit++;
97 }
98 }
88fb5e50 99 // Add high gain digit only if it is not saturated;
100 // replace low gain digit by a high gain one
101 else {
88fb5e50 102 for (iOldDigit=iDigit-1; iOldDigit>=0; iOldDigit--) {
77ea1c6f 103 AliPHOSDigit * dig = dynamic_cast<AliPHOSDigit*>(digits->At(iOldDigit)) ;
104 if (dig->GetId() == absId) {
34827459 105 dig->SetTime(time) ; //Replace time in any case
106 dig->SetEnergy((Float_t)decoder->GetEnergy()) ;
107// if (decoder->GetEnergy() > 0){ //replace Energy only if not saturated
108// dig->SetEnergy((Float_t)decoder->GetEnergy()) ;
109// }
88fb5e50 110 seen = kTRUE;
111 break;
112 }
113 }
114 if (!seen) {
115 new((*digits)[iDigit]) AliPHOSDigit(-1,absId,(Float_t)decoder->GetEnergy(),time);
116 iDigit++;
117 }
118 }
77ea1c6f 119
88fb5e50 120 }
121
122 digits->Compress();
123 digits->Sort();
124}