another memory leak fix (incorrect push)
[u/mrichter/AliRoot.git] / AD / AliADReconstructor.cxx
CommitLineData
5e319bd5 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: AliADReconstructor.cxx 20956 2007-09-26 14:22:18Z mrodrigu $ */
17//////////////////////////////////////////////////////////////////////////////
18// //
19// Class for AD reconstruction //
20//////////////////////////////////////////////////////////////////////////////
fcb6ff8b 21#include <TParameter.h>
5e319bd5 22
23#include "AliRawReader.h"
fcb6ff8b 24#include "AliGRPObject.h"
25#include "AliCDBManager.h"
26#include "AliCDBStorage.h"
27#include "AliCDBEntry.h"
28#include "AliESDEvent.h"
5e319bd5 29
30#include "AliADReconstructor.h"
5e319bd5 31#include "AliADdigit.h"
7b752923 32#include "AliESDAD.h"
fcb6ff8b 33#include "AliADConst.h"
34#include "AliADCalibData.h"
35#include "AliADRawStream.h"
5e319bd5 36
37ClassImp(AliADReconstructor)
fcb6ff8b 38//_____________________________________________________________________________
5e319bd5 39AliADReconstructor:: AliADReconstructor():
40 AliReconstructor(),
41 fESDAD(0x0),
fcb6ff8b 42 fCalibData(NULL),
5e319bd5 43 fDigitsArray(0)
fcb6ff8b 44
5e319bd5 45{
fcb6ff8b 46 fCalibData = GetCalibData();
5e319bd5 47 // Default constructor
48 // Get calibration data
49
50}
51
52//_____________________________________________________________________________
53AliADReconstructor& AliADReconstructor::operator =
54 (const AliADReconstructor& /*reconstructor*/)
55{
56// assignment operator
57
58 Fatal("operator =", "assignment operator not implemented");
59 return *this;
60}
61
62//_____________________________________________________________________________
63AliADReconstructor::~AliADReconstructor()
64{
65// destructor
66 delete fESDAD;
67 delete fDigitsArray;
68}
69
70//_____________________________________________________________________________
71void AliADReconstructor::Init()
72{
73// initializer
74 fESDAD = new AliESDAD;
75}
76
fcb6ff8b 77//_____________________________________________________________________________
78void AliADReconstructor::ConvertDigits(AliRawReader* rawReader, TTree* digitsTree) const
5e319bd5 79{
fcb6ff8b 80// converts RAW to digits
81
82 if (!digitsTree) {
83 AliError("No digits tree!");
84 return;
85 }
86
87 if (!fDigitsArray){
88 fDigitsArray = new TClonesArray("AliADdigit", 16);
89 digitsTree->Branch("ADDigit", &fDigitsArray);
90 }
91
92 rawReader->Reset();
93 AliADRawStream rawStream(rawReader);
94 if (rawStream.Next()) {
95
96 Int_t aBBflagsADA = 0;
97 Int_t aBBflagsADC = 0;
98 Int_t aBGflagsADA = 0;
99 Int_t aBGflagsADC = 0;
100
101 for(Int_t iChannel=0; iChannel < 16; ++iChannel) {
102 Int_t offlineCh = rawStream.GetOfflineChannel(iChannel);
103 // ADC charge samples
104 Short_t chargeADC[kNClocks];
105 for(Int_t iClock=0; iClock < kNClocks; ++iClock) {
106 chargeADC[iClock] = rawStream.GetPedestal(iChannel,iClock);
107 }
108 // Integrator flag
109 Bool_t integrator = rawStream.GetIntegratorFlag(iChannel,kNClocks/2);
110 // Beam-beam and beam-gas flags
111 if(offlineCh<8) {
112 if (rawStream.GetBBFlag(iChannel,kNClocks/2)) aBBflagsADC |= (1 << offlineCh);
113 if (rawStream.GetBGFlag(iChannel,kNClocks/2)) aBGflagsADC |= (1 << offlineCh);
114 } else {
115 if (rawStream.GetBBFlag(iChannel,kNClocks/2)) aBBflagsADA |= (1 << (offlineCh-32));
116 if (rawStream.GetBGFlag(iChannel,kNClocks/2)) aBGflagsADA |= (1 << (offlineCh-32));
117 }
118 // HPTDC data (leading time and width)
119 Int_t board = AliADCalibData::GetBoardNumber(offlineCh);
120 Float_t time = rawStream.GetTime(iChannel)*fCalibData->GetTimeResolution(board);
121 Float_t width = rawStream.GetWidth(iChannel)*fCalibData->GetWidthResolution(board);
122 // Add a digit
123 if(!fCalibData->IsChannelDead(iChannel)){
124 new ((*fDigitsArray)[fDigitsArray->GetEntriesFast()]) AliADdigit(offlineCh, time, width,integrator, chargeADC);
125 }
126 }
127 // Store the BB and BG flags in the digits tree (user info)
128 digitsTree->GetUserInfo()->Add(new TParameter<int>("BBflagsADA",aBBflagsADA));
129 digitsTree->GetUserInfo()->Add(new TParameter<int>("BBflagsADC",aBBflagsADC));
130 digitsTree->GetUserInfo()->Add(new TParameter<int>("BGflagsADA",aBGflagsADA));
131 digitsTree->GetUserInfo()->Add(new TParameter<int>("BGflagsADC",aBGflagsADC));
132
133
134 digitsTree->Fill();
135 }
136
137 fDigitsArray->Clear();
5e319bd5 138
5e319bd5 139}
140
fcb6ff8b 141//_____________________________________________________________________________
5e319bd5 142void AliADReconstructor::FillESD(TTree* digitsTree, TTree* /*clustersTree*/,AliESDEvent* esd) const
143{
144
145 printf("Running AD Reconstruction \n");
146
147 // fills ESD with AD Digits
148
149 if (!digitsTree)
150 {
151 AliError("No digits tree!");
152 return;
153 }
154
155 TBranch* digitBranch = digitsTree->GetBranch("ADdigit");
156 if (!digitBranch) {
157 AliError("No AD digits branch found!");
158 return;
159 }
160 digitBranch->SetAddress(&fDigitsArray);
161
162 digitsTree->GetEvent(0);
163
164 Bool_t ADHits[16];
165 for(Int_t i = 0; i < 16; i++) { ADHits[i] = kFALSE; }
166
167 Int_t nDigits = fDigitsArray->GetEntriesFast();
168
169 for (Int_t d=0; d<nDigits; d++) {
170 AliADdigit* digit = (AliADdigit*) fDigitsArray->At(d);
aa8120bb 171 Int_t module = digit->PMNumber();
5e319bd5 172 // printf("AD Module: %d\n",module);
173 ADHits[module] = kTRUE;
174 }
175 if (!esd) {
176 AliError("NO AD ESD branch found!");
177 return;
178}
179 fESDAD->SetADBitCell(ADHits);
180
181 if (esd)
182 {
183 AliDebug(1, Form("Writing AD data to ESD Tree"));
184 esd->SetADData(fESDAD);
185 }
186
187 fDigitsArray->Clear();
188}
189
fcb6ff8b 190//_____________________________________________________________________________
191AliADCalibData* AliADReconstructor::GetCalibData() const
192{
193 AliCDBManager *man = AliCDBManager::Instance();
194
195 AliCDBEntry *entry=0;
196
197 //entry = man->Get("AD/Calib/Data");
198 //if(!entry){
199 //AliWarning("Load of calibration data from default storage failed!");
200 //AliWarning("Calibration data will be loaded from local storage ($ALICE_ROOT)");
201
202 man->SetDefaultStorage("local://$ALICE_ROOT/OCDB");
203 man->SetRun(1);
204 entry = man->Get("AD/Calib/Data");
205 //}
206 // Retrieval of data in directory AD/Calib/Data:
207
208 AliADCalibData *calibdata = 0;
209
210 if (entry) calibdata = (AliADCalibData*) entry->GetObject();
211 if (!calibdata) AliFatal("No calibration data from calibration database !");
212
213 return calibdata;
214}
215
5e319bd5 216