ATLAS Offline Software
Loading...
Searching...
No Matches
TileTriggerDefaultCalibTool.cxx
Go to the documentation of this file.
1/*
2 Copyright (C) 2002-2025 CERN for the benefit of the ATLAS collaboration
3*/
4
5#include "TileTriggerDefaultCalibTool.h"
6
7#include "xAODTrigL1Calo/TriggerTowerContainer.h"
8#include "TrigT1CaloCalibToolInterfaces/IL1CaloTTIdTools.h"
9#include "TrigT1CaloCalibToolInterfaces/IL1CaloOfflineTriggerTowerTools.h"
10
11#include "Identifier/HWIdentifier.h"
12#include "StoreGate/ReadHandle.h"
13
14#include "CaloIdentifier/CaloLVL1_ID.h"
15#include "CaloIdentifier/TileID.h"
16
17#include "TileCalibBlobObjs/TileCalibUtils.h"
18#include "TileIdentifier/TileHWID.h"
19#include "TileEvent/TileRawChannelContainer.h"
20#include "TileConditions/TileCablingService.h"
21
22#include "TFile.h"
23#include "TTree.h"
24#include "TString.h"
25
26#include <vector>
27#include <cmath>
28
29
30namespace {
31template <class T>
32inline
33T square(T x) { return x*x; }
34}
35
36
37TileTriggerDefaultCalibTool::TileTriggerDefaultCalibTool(const std::string& type, const std::string& name,const IInterface* pParent):
38 AthAlgTool(type, name, pParent)
39 , m_TT_ID(nullptr)
40 , m_tileHWID(nullptr)
41 , m_tileID(nullptr)
42 , m_tileCablingService(nullptr)
43 , m_charge(0)
44 , m_ipmt(0)
45 , m_ipmtCount(0)
46 , m_ipmtOld(0)
47 , m_DACvalue(0)
48 , m_nEvtGlobal(0)
49{
50 declareInterface<ITileCalibTool>( this );
51 declareProperty("MaxNTriggerTowers", m_maxNTT=7200);
52 declareProperty("NtupleID", m_ntupleID="h3000");
53 declareProperty("NumEventPerPMT", m_nevpmt=195); // changed from 200 to 195
54
55 m_meanTile = new float[Tile::MAX_ROS][Tile::MAX_DRAWER][Tile::MAX_CHAN]();
56 m_rmsTile = new float[Tile::MAX_ROS][Tile::MAX_DRAWER][Tile::MAX_CHAN]();
57 m_meanTileDAC = new float[Tile::MAX_ROS][Tile::MAX_DRAWER][Tile::MAX_CHAN]();
58 m_rmsTileDAC = new float[Tile::MAX_ROS][Tile::MAX_DRAWER][Tile::MAX_CHAN]();
59 m_ietaTile = new int[Tile::MAX_ROS][Tile::MAX_DRAWER][Tile::MAX_CHAN]();
60 m_iphiTile = new int[Tile::MAX_ROS][Tile::MAX_DRAWER][Tile::MAX_CHAN]();
61 m_ipmtTile = new int[Tile::MAX_ROS][Tile::MAX_DRAWER][Tile::MAX_CHAN]();
62 m_nEvtTile = new int[Tile::MAX_ROS][Tile::MAX_DRAWER][Tile::MAX_CHAN]();
63
64 m_meanL1Calo = new float[Tile::MAX_ROS][Tile::MAX_DRAWER][Tile::MAX_CHAN]();
65 m_rmsL1Calo = new float[Tile::MAX_ROS][Tile::MAX_DRAWER][Tile::MAX_CHAN]();
66 m_meanL1CaloDAC = new float[Tile::MAX_ROS][Tile::MAX_DRAWER][Tile::MAX_CHAN]();
67 m_rmsL1CaloDAC = new float[Tile::MAX_ROS][Tile::MAX_DRAWER][Tile::MAX_CHAN]();
68 m_ietaL1Calo = new int[Tile::MAX_ROS][Tile::MAX_DRAWER][Tile::MAX_CHAN]();
69 m_iphiL1Calo = new int[Tile::MAX_ROS][Tile::MAX_DRAWER][Tile::MAX_CHAN]();
70 m_ipmtL1Calo = new int[Tile::MAX_ROS][Tile::MAX_DRAWER][Tile::MAX_CHAN]();
71 m_nEvtL1Calo = new int[Tile::MAX_ROS][Tile::MAX_DRAWER][Tile::MAX_CHAN]();
72
73 m_meanTileL1Calo = new float[Tile::MAX_ROS][Tile::MAX_DRAWER][Tile::MAX_CHAN]();
74 m_rmsTileL1Calo = new float[Tile::MAX_ROS][Tile::MAX_DRAWER][Tile::MAX_CHAN]();
75}
76
77TileTriggerDefaultCalibTool::~TileTriggerDefaultCalibTool()
78{
79 delete[] m_meanTile;
80 delete[] m_rmsTile;
81 delete[] m_meanTileDAC;
82 delete[] m_rmsTileDAC;
83 delete[] m_ietaTile;
84 delete[] m_iphiTile;
85 delete[] m_ipmtTile;
86 delete[] m_nEvtTile;
87 delete[] m_meanL1Calo;
88 delete[] m_rmsL1Calo;
89 delete[] m_meanL1CaloDAC;
90 delete[] m_rmsL1CaloDAC;
91 delete[] m_ietaL1Calo;
92 delete[] m_iphiL1Calo;
93 delete[] m_ipmtL1Calo;
94 delete[] m_nEvtL1Calo;
95 delete[] m_meanTileL1Calo;
96 delete[] m_rmsTileL1Calo;
97}
98
99StatusCode TileTriggerDefaultCalibTool::initialize()
100{
101 ATH_MSG_INFO ( "initialize()" );
102
103 // Initialize arrays for results
104 m_charge = 0;
105 m_ipmt = 0;
106 m_ipmtOld = 0;
107 m_ipmtCount = 0;
108
109 ATH_CHECK( m_rawChannelContainerKey.initialize() );
110 ATH_CHECK( m_triggerTowerContainerKey.initialize() );
111
112 ATH_CHECK( m_l1CaloTTIdTools.retrieve() );
113 ATH_MSG_DEBUG("L1CaloTTIdTools retrieved");
114
115 ATH_CHECK( detStore()->retrieve(m_TT_ID) );
116 ATH_CHECK( detStore()->retrieve(m_tileHWID) );
117 ATH_CHECK( detStore()->retrieve(m_tileID) );
118
119 m_tileCablingService = TileCablingService::getInstance();
120
121 //=== get TileCondToolEmscale
122 CHECK( m_tileToolEmscale.retrieve() );
123
124 CHECK( m_dqStatusKey.initialize() );
125
126 return StatusCode::SUCCESS;
127}
128
129StatusCode TileTriggerDefaultCalibTool::initNtuple(int runNumber, int runType, TFile * rootFile )
130{
131 ATH_MSG_INFO ( "initialize(" << runNumber << "," << runType << "," << rootFile << ")" );
132 return StatusCode::SUCCESS;
133}
134
135StatusCode TileTriggerDefaultCalibTool::execute()
136{
137 ATH_MSG_DEBUG ( "executeTrigger()" );
138
139 const EventContext& ctx = Gaudi::Hive::currentContext();
140 const TileDQstatus* dqStatus = SG::makeHandle (m_dqStatusKey, ctx).get();
141
142 // Get TileRawChannelContainer
143 SG::ReadHandle<TileRawChannelContainer> container(m_rawChannelContainerKey);
144 ATH_CHECK( container.isValid() );
145
146 ATH_MSG_DEBUG ( "second executeTrigger()" );
147 // declare an array of pmt id for the trigger tower loop
148 //coverity[STACK_USE]
149 HWIdentifier chanIds[2][16][64][6];
150 memset(chanIds, 0, sizeof(chanIds));
151
152 // Get event's ext CIS parameters
153 const uint32_t *cispar = dqStatus->cispar();
154
155 // Mapping pmt2chan
156 int chan_bar[54] = {4 ,1 ,2 ,3 ,0 ,-1 ,
157 8 ,5 ,6 ,7 ,13,-1 ,
158 10,9 ,12,11,14,-1 ,
159 18,15,16,17,25,-1 ,
160 20,19,22,21,24,-1 ,
161 26,23,28,27,41,-1 ,
162 32,29,34,33,44,-1 ,
163 38,35,40,39,-1,-1 ,
164 36,37,42,47,46,45 };
165
166 int chan_ext[36] = {4, 5 ,2 ,3 ,16,-1 ,
167 8, 9 ,17,-1,-1,-1 ,
168 6, 7 ,14,15,38,-1 ,
169 10,11,22,23,37,-1 ,
170 20,21,30,35,-1,-1 ,
171 31,32,39,36,41,40 };
172
173 ATH_MSG_DEBUG ( "cispar[16] " << cispar[16] << ", cispas[17] " << cispar[17] << ", cispar[18] " << cispar[18] );
174 if (cispar[16] == 0x07){
175 if (cispar[18]>5) {
176 ATH_MSG_WARNING ( "Bad CISpar detected, using pmt index and charge from previous event: " << m_ipmt << " " << m_charge );
177 } else {
178 m_charge = cispar[17];
179 m_ipmt = cispar[18];
180 //m_itower = cispar[19];
181 //m_idrawer = cispar[20];
182 }
183 }
184 else if (cispar[16] == 0x107){ // bad CISpar
185 ATH_MSG_WARNING ( "Bad CISpar detected, using pmt index and charge from previous event: " << m_ipmt << " " << m_charge );
186 }
187 else
188 return StatusCode::SUCCESS;
189
190 ATH_MSG_DEBUG ( "Pattern/ipmt/charge: " << cispar[16] << " " << m_ipmt << " " << m_charge );
191 if (m_charge<1.0e-20) {
192 ATH_MSG_WARNING ( "Bad charge " << m_charge << " - skipping event" );
193 return StatusCode::SUCCESS;
194 }
195
196 if (m_ipmtOld != m_ipmt){
197 if (m_ipmt < m_ipmtOld) return StatusCode::SUCCESS;
198 m_ipmtCount = 1;
199 m_ipmtOld = m_ipmt;
200 }
201 else
202 m_ipmtCount += 1;
203
204 if (m_ipmtCount > m_nevpmt) return StatusCode::SUCCESS; //takes into account only the first m_nevpmt events for each pmt to avoid problem with wrong CISpar
205
206 m_DACvalue = m_charge;
207
208 // Create iterator over RawChannelContainer for Tile
209 TileRawChannelContainer::const_iterator itColl = (*container).begin();
210 TileRawChannelContainer::const_iterator itCollEnd = (*container).end();
211 TileRawChannelCollection::const_iterator it, itEnd;
212 TileRawChannelUnit::UNIT RChUnit = container->get_unit();
213
214 // Go through all TileRawChannelCollections
215 for(; itColl != itCollEnd; ++itColl) {
216
217 //go through all TileRawChannels in collection (loop over all channels)
218 it = (*itColl)->begin();
219 itEnd = (*itColl)->end();
220
221 for(; it != itEnd; ++it) {
222
223 // get hardware id to identify adc channel
224 HWIdentifier hwid=(*it)->adc_HWID();
225 //ATH_MSG_DEBUG("HWID1: "<< m_tileHWID->to_string(hwid));
226 int ros = m_tileHWID->ros(hwid); // LBA=1 LBC=2 EBA=3 EBC=4
227 int drawer = m_tileHWID->drawer(hwid); // 0 to 63
228 int chan = m_tileHWID->channel(hwid); // 0 to 47 channel not PMT
229 int gain = m_tileHWID->adc(hwid); // low=0 high=1
230 unsigned int drawerIdx = TileCalibUtils::getDrawerIdx(ros, drawer);
231
232 // Get the PMT Id, create on the fly a mapping to the trigger tower id
233 // this mapping is to be used below in the loop over all trigger towers
234 Identifier tt_id = (*it)->tt_ID();
235 if (!tt_id.is_valid()) continue;
236 //Identifier pmt_id = (*it)->pmt_ID();
237
238 int pos_neg_z = m_TT_ID->pos_neg_z(tt_id);
239 int ieta = m_TT_ID->eta(tt_id);
240 int iphi = m_TT_ID->phi(tt_id);
241
242 if (ros < 3) {
243 if (ieta < 0 || ieta > 8) continue;
244 } else {
245 if (ieta < 9 || ieta > 14) continue;
246 }
247
248 if((ros==1) || (ros==2)){
249 // (Manually fold constants to avoid cppcheck warning.)
250 if (chan != chan_bar[(ieta*6)/*-1*/ + m_ipmt/*+1*/]) continue; // check if the chan is firing
251 }
252 else if ((ros==3) || (ros==4)){
253 if ((chan != chan_ext[((ieta-9)*6)/*-1*/ + m_ipmt/*+1*/]) || (chan == 0)) continue; // check if the chan is firing
254 }
255 else
256 continue;
257
258 float amp = (*it)->amplitude();
259 float amp_pC = m_tileToolEmscale->channelCalib(drawerIdx, chan, gain, amp, RChUnit, TileRawChannelUnit::PicoCoulombs);
260 //float amp_GeV = m_tileToolEmscale->channelCalib(drawerIdx, chan, gain, amp, RChUnit, TileRawChannelUnit::MegaElectronVolts) * 0.001;
261 //float pC2GeV = m_tileToolEmscale->channelCalib(drawerIdx, chan, 0, 1.0, TileRawChannelUnit::PicoCoulombs, TileRawChannelUnit::MegaElectronVolts) * 0.001;
262 //log << MSG::DEBUG << "ros " << ros << ", pos_neg_z " << pos_neg_z << ", drawer " << drawer <<" ieta " << ieta << ", chan " << chan << ", amplitude " << amp << endreq;
263
264 m_meanTile[ros][drawer][chan] += amp_pC;
265 m_rmsTile[ros][drawer][chan] += square(amp_pC);
266 m_meanTileDAC[ros][drawer][chan] += (amp_pC-m_charge);
267 m_rmsTileDAC[ros][drawer][chan] += square(amp_pC-m_charge);
268 m_ietaTile[ros][drawer][chan] = ieta;
269 m_iphiTile[ros][drawer][chan] = iphi;
270 m_ipmtTile[ros][drawer][chan] = m_ipmt;
271 ++m_nEvtTile[ros][drawer][chan];
272
273 if (pos_neg_z < 0) pos_neg_z = 0;
274 chanIds[pos_neg_z][ieta][iphi][m_ipmt] = hwid;
275 }
276
277 } // end of loop over raw channels for Tile
278
279 // loop over all L1Calo trigger channels, calculate the average and RMS
280 SG::ReadHandle<xAOD::TriggerTowerContainer> triggerTowers(m_triggerTowerContainerKey);
281 ATH_CHECK( triggerTowers.isValid() );
282
283 int ntt = 0;
284
285 for (const xAOD::TriggerTower* triggerTower : *triggerTowers) {
286
287 if(ntt>m_maxNTT) break;
288
289 if (!triggerTower->sampling()) continue; // is not hadronic
290
291 double eta = triggerTower->eta();
292 double phi = triggerTower->phi();
293 double tt_ene = 0.; //triggerTower->e();
294
295 // temporary, compute id out of the tower position
296 int pos_neg_z = m_l1CaloTTIdTools->pos_neg_z(eta);
297 int region = m_l1CaloTTIdTools->regionIndex(eta);
298 int ieta = m_l1CaloTTIdTools->etaIndex(eta);
299 int iphi = m_l1CaloTTIdTools->phiIndex(eta, phi);
300 Identifier tt_id = m_TT_ID->tower_id(pos_neg_z, 1, region, ieta, iphi);
301 //Not used, but could be useful, eg. coolId is used to get proper pedestal value
302 //Identifier id = m_ttTool->identifier(triggerTower->eta(), triggerTower->phi(), 1); //alternative way to get tt_id
303 //unsigned int coolId = m_ttTool->channelID(tt_id).id();
304
305 // check if it is a Tile tower
306 if (tt_id.is_valid() && m_TT_ID->is_tile(tt_id))
307 ++ntt;
308 else continue;
309
310 //compute L1Calo energy from ADC:
311 // - take max ADC value
312 // - subtract pedestal (about 32 ADC counts), apply 0.25 conversion factor to GeV
313 const std::vector<uint16_t>& adc = triggerTower->adc();
314 uint16_t max = *(std::max_element(adc.begin(), adc.end()));
315 uint16_t ped = adc.front();
316 tt_ene = (max - ped) * 0.25;
317 if (tt_ene < 0.) tt_ene = 0.;
318
319 // check boundaries
320 if ((ieta < 0) || (ieta > 14)) continue;
321
322 if (pos_neg_z < 0) pos_neg_z = 0;
323 HWIdentifier hwid = chanIds[pos_neg_z][ieta][iphi][m_ipmt];
324 int ros = m_tileHWID->ros(hwid); // LBA=1 LBC=2 EBA=3 EBC=4
325 if(ros >= 5 || ros <=0){continue;} // ignoring special channels, like those connected to inner MBTS in EBA/EBC 39-42,55-58
326 int drawer = m_tileHWID->drawer(hwid); // 0 to 63
327 int chan = m_tileHWID->channel(hwid); // 0 to 47 channel not PMT
328 unsigned int drawerIdx = TileCalibUtils::getDrawerIdx(ros, drawer);
329 float pC2GeV = m_tileToolEmscale->channelCalib(drawerIdx, chan, 0, 1.0, TileRawChannelUnit::PicoCoulombs, TileRawChannelUnit::MegaElectronVolts) * 0.001;
330 double tt_ene_pC = tt_ene / pC2GeV;
331
332 //comparing results from triggerTower->e() and energy computed from ADC (default method)
333 if (triggerTower->e() == 0)
334 ATH_MSG_DEBUG ( "Tower partition " << ros
335 << ", drawer "<< drawer
336 << ", chan: " << chan
337 << ", eta "<< eta
338 << ", phi: " << phi
339 << ", amplitude : " << tt_ene
340 << ", old " << triggerTower->e()
341 << ", ratio old/new " << (tt_ene == 0 ? 0 : triggerTower->e() / tt_ene ));
342
343 if ((ros==1) || (ros==2)){
344 if (chan != chan_bar[(ieta*6)-1 + m_ipmt+1]) continue; // check if the chan is firing
345 }
346 else if ((ros==3) || (ros==4)) {
347 if ((chan != chan_ext[((ieta-9)*6)-1 + m_ipmt+1]) || (chan == 0)) continue; // check if the chan is firing
348 }
349 else continue;
350
351 m_meanL1Calo[ros][drawer][chan] += tt_ene_pC;
352 m_rmsL1Calo[ros][drawer][chan] += square(tt_ene_pC);
353 m_meanL1CaloDAC[ros][drawer][chan] += tt_ene_pC - m_charge;
354 m_rmsL1CaloDAC[ros][drawer][chan] += square(tt_ene_pC - m_charge);
355 m_ietaL1Calo[ros][drawer][chan] = ieta;
356 m_iphiL1Calo[ros][drawer][chan] = iphi;
357 m_ipmtL1Calo[ros][drawer][chan] = m_ipmt;
358 ++m_nEvtL1Calo[ros][drawer][chan];
359
360 }// end of trigger tower loop for L1Calo
361
362 ++m_nEvtGlobal; // Number of MaxEvt from Athena
363
364 for (int ros=0;ros<5;ros++) {
365 for (int drawer=0;drawer<64;drawer++) {
366 for (int chan=0;chan<48;chan++) {
367 // ???
368 m_rmsTileL1Calo[ros][drawer][chan] += m_meanTile[ros][drawer][chan]*m_meanL1Calo[ros][drawer][chan];
369 }
370 }
371 }
372
373 return StatusCode::SUCCESS;
374
375}
376
377StatusCode TileTriggerDefaultCalibTool::finalizeCalculations()
378{
379 ATH_MSG_INFO ( "finalizeCalculations()" );
380
381 for (int ros=0;ros<5;ros++) {
382 for (int drawer=0;drawer<64;drawer++) {
383 for (int chan=0;chan<48;chan++) {
384 if (m_nEvtTile[ros][drawer][chan] > 0){
385 m_meanTile[ros][drawer][chan] /= m_nEvtTile[ros][drawer][chan];
386 m_rmsTile[ros][drawer][chan] /= m_nEvtTile[ros][drawer][chan];
387 m_meanTileDAC[ros][drawer][chan] /= m_nEvtTile[ros][drawer][chan];
388 m_rmsTileDAC[ros][drawer][chan] /= m_nEvtTile[ros][drawer][chan];
389 }
390
391 if(m_rmsTile[ros][drawer][chan] <= square(m_meanTile[ros][drawer][chan]))
392 m_rmsTile[ros][drawer][chan] = 0;
393 else
394 m_rmsTile[ros][drawer][chan] = sqrt(m_rmsTile[ros][drawer][chan] - square(m_meanTile[ros][drawer][chan]));
395
396 if(m_rmsTileDAC[ros][drawer][chan] <= square(m_meanTileDAC[ros][drawer][chan]))
397 m_rmsTileDAC[ros][drawer][chan] = 0;
398 else
399 m_rmsTileDAC[ros][drawer][chan] = sqrt(m_rmsTileDAC[ros][drawer][chan] - square(m_meanTileDAC[ros][drawer][chan]));
400
401 if (m_nEvtL1Calo[ros][drawer][chan] > 0){
402 m_meanL1Calo[ros][drawer][chan] /= m_nEvtL1Calo[ros][drawer][chan];
403 m_rmsL1Calo[ros][drawer][chan] /= m_nEvtL1Calo[ros][drawer][chan];
404 m_meanL1CaloDAC[ros][drawer][chan] /= m_nEvtL1Calo[ros][drawer][chan];
405 m_rmsL1CaloDAC[ros][drawer][chan] /= m_nEvtL1Calo[ros][drawer][chan];
406 }
407
408 if(m_rmsL1Calo[ros][drawer][chan] <= square(m_meanL1Calo[ros][drawer][chan]))
409 m_rmsL1Calo[ros][drawer][chan] = 0;
410 else
411 m_rmsL1Calo[ros][drawer][chan] = sqrt(m_rmsL1Calo[ros][drawer][chan] - square(m_meanL1Calo[ros][drawer][chan]));
412
413 if(m_rmsL1CaloDAC[ros][drawer][chan] <= square(m_meanL1CaloDAC[ros][drawer][chan]))
414 m_rmsL1CaloDAC[ros][drawer][chan] = 0;
415 else
416 m_rmsL1CaloDAC[ros][drawer][chan] = sqrt(m_rmsL1CaloDAC[ros][drawer][chan] - square(m_meanL1CaloDAC[ros][drawer][chan]));
417
418 m_meanTileL1Calo[ros][drawer][chan] = m_meanTile[ros][drawer][chan] - m_meanL1Calo[ros][drawer][chan];
419 if(m_rmsTileL1Calo[ros][drawer][chan] <= square(m_meanTileL1Calo[ros][drawer][chan]))
420 m_rmsTileL1Calo[ros][drawer][chan] = 0;
421 else
422 m_rmsTileL1Calo[ros][drawer][chan] = sqrt(m_rmsTileL1Calo[ros][drawer][chan] - square(m_meanTileL1Calo[ros][drawer][chan]));
423 }
424 }
425 }
426
427 return StatusCode::SUCCESS;
428}
429
430StatusCode TileTriggerDefaultCalibTool::writeNtuple(int runNumber, int runType, TFile * rootFile)
431{
432 ATH_MSG_INFO ( "writeNtuple(" << runNumber << "," << runType << "," << rootFile << ")" );
433
434 TTree *t = new TTree(m_ntupleID.c_str(), "TileCalib-Ntuple");
435 TString ind = TString::Format("[%d][%d][%d]",Tile::MAX_ROS,Tile::MAX_DRAWER,Tile::MAX_CHAN);
436
437 t->Branch("meanTile",m_meanTile,ind.Format("meanTile%s/F",(const char *)ind));
438 t->Branch("rmsTile",m_rmsTile,ind.Format("rmsTile%s/F",(const char *)ind));
439 t->Branch("meanTileDAC",m_meanTileDAC,ind.Format("meanTileDAC%s/F",(const char *)ind));
440 t->Branch("rmsTileDAC",m_rmsTileDAC,ind.Format("rmsTileDAC%s/F",(const char *)ind));
441 t->Branch("ietaTile",m_ietaTile,ind.Format("ietaTile%s/I",(const char *)ind));
442 t->Branch("iphiTile",m_iphiTile,ind.Format("iphiTile%s/I",(const char *)ind));
443 t->Branch("ipmtTile",m_ipmtTile,ind.Format("ipmtTile%s/I",(const char *)ind));
444 t->Branch("nEvtTile",m_nEvtTile,ind.Format("nEvtTile%s/I",(const char *)ind));
445 t->Branch("meanL1Calo",m_meanL1Calo,ind.Format("meanL1Calo%s/F",(const char *)ind));
446 t->Branch("rmsL1Calo",m_rmsL1Calo,ind.Format("rmsL1Calo%s/F",(const char *)ind));
447 t->Branch("meanL1CaloDAC",m_meanL1CaloDAC,ind.Format("meanL1CaloDAC%s/F",(const char *)ind));
448 t->Branch("rmsL1CaloDAC",m_rmsL1CaloDAC,ind.Format("rmsL1CaloDAC%s/F",(const char *)ind));
449 t->Branch("ietaL1Calo",m_ietaL1Calo,ind.Format("ietaL1Calo%s/I",(const char *)ind));
450 t->Branch("iphiL1Calo",m_iphiL1Calo,ind.Format("iphiL1Calo%s/I",(const char *)ind));
451 t->Branch("ipmtL1Calo",m_ipmtL1Calo,ind.Format("ipmtL1Calo%s/I",(const char *)ind));
452 t->Branch("nEvtL1Calo",m_nEvtL1Calo,ind.Format("nEvtL1Calo%s/I",(const char *)ind));
453 t->Branch("meanTileL1Calo",m_meanTileL1Calo,ind.Format("meanTileL1Calo%s/F",(const char *)ind));
454 t->Branch("rmsTileL1Calo",m_rmsTileL1Calo,ind.Format("rmsTileL1Calo%s/F",(const char *)ind));
455 t->Branch("nEvtGlobal",&m_nEvtGlobal,"nEvtGlobal/I");
456 t->Branch("DACvalue",&m_DACvalue,"DACvalue/F");
457
458 // Fill with current values (i.e. tree will have only one entry for this whole run)
459 t->Fill();
460 t->Write();
461
462 return StatusCode::SUCCESS;
463}
464
465StatusCode TileTriggerDefaultCalibTool::finalize()
466{
467 ATH_MSG_DEBUG ( "finalizeTrigger()" );
468 return StatusCode::SUCCESS;
469}
eta
Scalar eta() const
pseudorapidity method
Definition AmgMatrixBasePlugin.h:83
phi
Scalar phi() const
phi method
Definition AmgMatrixBasePlugin.h:67
ATH_CHECK
#define ATH_CHECK
Evaluate an expression and check for errors.
Definition AthCheckMacros.h:40
ATH_MSG_INFO
#define ATH_MSG_INFO(x)
Definition AthMsgStreamMacros.h:31
ATH_MSG_WARNING
#define ATH_MSG_WARNING(x)
Definition AthMsgStreamMacros.h:32
ATH_MSG_DEBUG
#define ATH_MSG_DEBUG(x)
Definition AthMsgStreamMacros.h:29
CaloLVL1_ID.h
CHECK
#define CHECK(...)
Evaluate an expression and check for errors.
Definition Control/AthenaKernel/AthenaKernel/errorcheck.h:422
ReadHandle.h
Handle class for reading from StoreGate.
x
#define x
max
#define max(a, b)
Definition cfImp.cxx:41
AthAlgTool::AthAlgTool
AthAlgTool(const std::string &type, const std::string &name, const IInterface *parent)
Constructor with parameters:
Definition AthAlgTool.cxx:16
AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T, V, H > &t)
Definition AthCommonDataStore.h:145
AthCommonDataStore< AthCommonMsg< AlgTool > >::detStore
const ServiceHandle< StoreGateSvc > & detStore() const
Definition AthCommonDataStore.h:95
DataVector::const_iterator
DataModel_detail::const_iterator< DataVector > const_iterator
Standard const_iterator.
Definition DataVector.h:838
IdentifiableContainerMT::const_iterator
Definition IdentifiableContainerMT.h:82
Identifier::is_valid
bool is_valid() const
Check if id is in a valid state.
SG::ReadHandle::isValid
virtual bool isValid() override final
Can the handle be successfully dereferenced?
TileCablingService::getInstance
static const TileCablingService * getInstance()
get pointer to service instance
Definition TileCablingService.cxx:24
TileCalibUtils::MAX_ROS
static const unsigned int MAX_ROS
Number of ROSs.
Definition TileCalibUtils.h:138
TileCalibUtils::MAX_DRAWER
static const unsigned int MAX_DRAWER
Number of drawers in ROS 1-4.
Definition TileCalibUtils.h:139
TileCalibUtils::getDrawerIdx
static unsigned int getDrawerIdx(unsigned int ros, unsigned int drawer)
Returns a drawer hash.
Definition TileCalibUtils.cxx:60
TileCalibUtils::MAX_CHAN
static const unsigned int MAX_CHAN
Number of channels in drawer.
Definition TileCalibUtils.h:141
TileDQstatus
Class that holds Data Quality fragment information and provides functions to extract the data quality...
Definition TileDQstatus.h:49
TileDQstatus::cispar
const uint32_t * cispar() const
CIS parameters.
Definition TileDQstatus.h:152
TileTriggerDefaultCalibTool::initNtuple
virtual StatusCode initNtuple(int runNumber, int runType, TFile *rootfile) override
Definition TileTriggerDefaultCalibTool.cxx:129
TileTriggerDefaultCalibTool::m_l1CaloTTIdTools
ToolHandle< LVL1::IL1CaloTTIdTools > m_l1CaloTTIdTools
Definition TileTriggerDefaultCalibTool.h:68
TileTriggerDefaultCalibTool::m_meanL1Calo
float(* m_meanL1Calo)[Tile::MAX_DRAWER][Tile::MAX_CHAN]
Definition TileTriggerDefaultCalibTool.h:84
TileTriggerDefaultCalibTool::m_rmsTileDAC
float(* m_rmsTileDAC)[Tile::MAX_DRAWER][Tile::MAX_CHAN]
Definition TileTriggerDefaultCalibTool.h:77
TileTriggerDefaultCalibTool::m_rmsTileL1Calo
float(* m_rmsTileL1Calo)[Tile::MAX_DRAWER][Tile::MAX_CHAN]
Definition TileTriggerDefaultCalibTool.h:94
TileTriggerDefaultCalibTool::finalizeCalculations
virtual StatusCode finalizeCalculations() override
Definition TileTriggerDefaultCalibTool.cxx:377
TileTriggerDefaultCalibTool::m_nEvtL1Calo
int(* m_nEvtL1Calo)[Tile::MAX_DRAWER][Tile::MAX_CHAN]
Definition TileTriggerDefaultCalibTool.h:91
TileTriggerDefaultCalibTool::m_rawChannelContainerKey
SG::ReadHandleKey< TileRawChannelContainer > m_rawChannelContainerKey
Definition TileTriggerDefaultCalibTool.h:63
TileTriggerDefaultCalibTool::m_ietaTile
int(* m_ietaTile)[Tile::MAX_DRAWER][Tile::MAX_CHAN]
Definition TileTriggerDefaultCalibTool.h:78
TileTriggerDefaultCalibTool::finalize
virtual StatusCode finalize() override
Definition TileTriggerDefaultCalibTool.cxx:465
TileTriggerDefaultCalibTool::m_meanTileDAC
float(* m_meanTileDAC)[Tile::MAX_DRAWER][Tile::MAX_CHAN]
Definition TileTriggerDefaultCalibTool.h:76
TileTriggerDefaultCalibTool::execute
virtual StatusCode execute() override
Definition TileTriggerDefaultCalibTool.cxx:135
TileTriggerDefaultCalibTool::writeNtuple
virtual StatusCode writeNtuple(int runNumber, int runType, TFile *rootfile) override
Definition TileTriggerDefaultCalibTool.cxx:430
TileTriggerDefaultCalibTool::m_ietaL1Calo
int(* m_ietaL1Calo)[Tile::MAX_DRAWER][Tile::MAX_CHAN]
Definition TileTriggerDefaultCalibTool.h:88
TileTriggerDefaultCalibTool::m_meanTileL1Calo
float(* m_meanTileL1Calo)[Tile::MAX_DRAWER][Tile::MAX_CHAN]
Definition TileTriggerDefaultCalibTool.h:93
TileTriggerDefaultCalibTool::initialize
virtual StatusCode initialize() override
Definition TileTriggerDefaultCalibTool.cxx:99
TileTriggerDefaultCalibTool::m_meanL1CaloDAC
float(* m_meanL1CaloDAC)[Tile::MAX_DRAWER][Tile::MAX_CHAN]
Definition TileTriggerDefaultCalibTool.h:86
TileTriggerDefaultCalibTool::m_ipmtL1Calo
int(* m_ipmtL1Calo)[Tile::MAX_DRAWER][Tile::MAX_CHAN]
Definition TileTriggerDefaultCalibTool.h:90
TileTriggerDefaultCalibTool::m_triggerTowerContainerKey
SG::ReadHandleKey< xAOD::TriggerTowerContainer > m_triggerTowerContainerKey
Definition TileTriggerDefaultCalibTool.h:65
TileTriggerDefaultCalibTool::m_rmsL1CaloDAC
float(* m_rmsL1CaloDAC)[Tile::MAX_DRAWER][Tile::MAX_CHAN]
Definition TileTriggerDefaultCalibTool.h:87
TileTriggerDefaultCalibTool::m_meanTile
float(* m_meanTile)[Tile::MAX_DRAWER][Tile::MAX_CHAN]
Definition TileTriggerDefaultCalibTool.h:74
TileTriggerDefaultCalibTool::m_rmsL1Calo
float(* m_rmsL1Calo)[Tile::MAX_DRAWER][Tile::MAX_CHAN]
Definition TileTriggerDefaultCalibTool.h:85
TileTriggerDefaultCalibTool::m_rmsTile
float(* m_rmsTile)[Tile::MAX_DRAWER][Tile::MAX_CHAN]
Definition TileTriggerDefaultCalibTool.h:75
TileTriggerDefaultCalibTool::m_tileCablingService
const TileCablingService * m_tileCablingService
Definition TileTriggerDefaultCalibTool.h:58
TileTriggerDefaultCalibTool::m_iphiTile
int(* m_iphiTile)[Tile::MAX_DRAWER][Tile::MAX_CHAN]
Definition TileTriggerDefaultCalibTool.h:79
TileTriggerDefaultCalibTool::m_nEvtTile
int(* m_nEvtTile)[Tile::MAX_DRAWER][Tile::MAX_CHAN]
Definition TileTriggerDefaultCalibTool.h:81
TileTriggerDefaultCalibTool::m_dqStatusKey
SG::ReadHandleKey< TileDQstatus > m_dqStatusKey
Definition TileTriggerDefaultCalibTool.h:61
TileTriggerDefaultCalibTool::m_ipmtTile
int(* m_ipmtTile)[Tile::MAX_DRAWER][Tile::MAX_CHAN]
Definition TileTriggerDefaultCalibTool.h:80
TileTriggerDefaultCalibTool::m_tileToolEmscale
ToolHandle< TileCondToolEmscale > m_tileToolEmscale
main Tile Calibration tool
Definition TileTriggerDefaultCalibTool.h:59
TileTriggerDefaultCalibTool::m_iphiL1Calo
int(* m_iphiL1Calo)[Tile::MAX_DRAWER][Tile::MAX_CHAN]
Definition TileTriggerDefaultCalibTool.h:89
TileTriggerDefaultCalibTool::TileTriggerDefaultCalibTool
TileTriggerDefaultCalibTool(const std::string &type, const std::string &name, const IInterface *pParent)
Definition TileTriggerDefaultCalibTool.cxx:37
SG::makeHandle
SG::ReadCondHandle< T > makeHandle(const SG::ReadCondHandleKey< T > &key, const EventContext &ctx=Gaudi::Hive::currentContext())
Definition ReadCondHandle.h:269
xAOD::TriggerTower
TriggerTower_v2 TriggerTower
Define the latest version of the TriggerTower class.
Definition Event/xAOD/xAODTrigL1Calo/xAODTrigL1Calo/TriggerTower.h:16
type
Generated on for ATLAS Offline Software by doxygen 1.14.0