minor changes in the normalization (Marian)
[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
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}