ATLAS Offline Software
Loading...
Searching...
No Matches
TRT_RDOAnalysis.cxx
Go to the documentation of this file.
1/*
2 Copyright (C) 2002-2026 CERN for the benefit of the ATLAS collaboration
3*/
4
5
6#include "TRT_RDOAnalysis.h"
7#include "StoreGate/ReadHandle.h"
8
9#include "TTree.h"
10#include "TString.h"
11
12#include <cmath>
13
14TRT_RDOAnalysis::TRT_RDOAnalysis(const std::string& name, ISvcLocator *pSvcLocator)
15 : AthAlgorithm(name, pSvcLocator)
16 , m_inputKey("TRT_RDOs")
17 , m_inputTruthKey("TRT_SDO_Map")
18 , m_trtID(nullptr)
19 , m_rdoID(0)
20 , m_rdoWord(0)
21 , m_barrelEndcap(0)
22 , m_phiModule(0)
23 , m_layerWheel(0)
24 , m_strawLayer(0)
25 , m_straw(0)
26 , m_highLevel(0)
27 , m_timeOverThreshold(0)
28 , m_driftTimeBin(0)
29 , m_trailEdge(0)
30 , m_firstBin(0)
31 , m_lastBin(0)
32 , m_sdoID(0)
33 , m_sdoWord(0)
34 , m_barrelEndcap_sdo(0)
35 , m_phiModule_sdo(0)
36 , m_layerWheel_sdo(0)
37 , m_strawLayer_sdo(0)
38 , m_straw_sdo(0)
39 , m_aboveThresh(0)
40 , m_deadChan(0)
41 , m_RODdata(0)
42 , m_validStraw(0)
43 , m_hit(0)
44 , m_barcode(0)
45 , m_eventIndex(0)
46 , m_charge(0)
47 , m_barcode_vec(0)
48 , m_eventIndex_vec(0)
49 , m_charge_vec(0)
50
51 , m_h_rdoID(0)
52 , m_h_rdoWord(0)
53 , m_h_barrelEndcap(0)
54 , m_h_phiModule(0)
55 , m_h_layerWheel(0)
56 , m_h_strawLayer(0)
57 , m_h_straw(0)
58 , m_h_ToT(0)
59 , m_h_ToT_HL(0)
60 , m_h_driftTimeBin(0)
61 , m_h_trailEdge(0)
62 , m_h_brlPhiMod(0)
63 , m_h_brlLayer(0)
64 , m_h_brlStrawLayer(0)
65 , m_h_brlStraw(0)
66 , m_h_brlToT(0)
67 , m_h_brlToT_HL(0)
68 , m_h_brlDriftTimeBin(0)
69 , m_h_brlTrailEdge(0)
70 , m_h_ecPhiMod(0)
71 , m_h_ecWheel(0)
72 , m_h_ecStrawLayer(0)
73 , m_h_ecStraw(0)
74 , m_h_ecToT(0)
75 , m_h_ecToT_HL(0)
76 , m_h_ecDriftTimeBin(0)
77 , m_h_ecTrailEdge(0)
78 , m_h_sdoID(0)
79 , m_h_sdoWord(0)
80 , m_h_barrelEndcap_sdo(0)
81 , m_h_phiModule_sdo(0)
82 , m_h_layerWheel_sdo(0)
83 , m_h_strawLayer_sdo(0)
84 , m_h_straw_sdo(0)
85 , m_h_barcode(0)
86 , m_h_eventIndex(0)
87 , m_h_charge(0)
88
89 , m_tree(0)
90 , m_ntupleFileName("/ntuples/file1")
91 , m_ntupleDirName("/TRT_RDOAnalysis/")
92 , m_ntupleTreeName("TRT_RDOAna")
93 , m_path("/TRT_RDOAnalysis/")
94 , m_thistSvc("THistSvc", name)
95{
96 declareProperty("InputKey", m_inputKey);
97 declareProperty("InputTruthKey", m_inputTruthKey);
98 declareProperty("NtupleFileName", m_ntupleFileName);
99 declareProperty("NtupleDirectoryName", m_ntupleDirName);
100 declareProperty("NtupleTreeName", m_ntupleTreeName);
101 declareProperty("HistPath", m_path);
102}
103
104StatusCode TRT_RDOAnalysis::initialize() {
105 ATH_MSG_DEBUG( "Initializing TRT_RDOAnalysis" );
106
107 // This will check that the properties were initialized
108 // properly by job configuration.
109 ATH_CHECK( m_inputKey.initialize() );
110 ATH_CHECK( m_inputTruthKey.initialize() );
111
112 // Grab TRT_ID helper
113 ATH_CHECK(detStore()->retrieve(m_trtID, "TRT_ID"));
114
115 // Grab Ntuple and histogramming service for tree
116 ATH_CHECK(m_thistSvc.retrieve());
117
118 // ntuple branches
119 m_tree = new TTree(TString(m_ntupleTreeName), "TRT_RDOAna");
120 std::string fullNtupleName = m_ntupleFileName + m_ntupleDirName + m_ntupleTreeName;
121 ATH_CHECK(m_thistSvc->regTree(fullNtupleName,m_tree));
122 if (m_tree) {
123 // TRT RDO
124 m_tree->Branch("rdoID", &m_rdoID);
125 m_tree->Branch("rdoWord", &m_rdoWord);
126 m_tree->Branch("barrelEndcap", &m_barrelEndcap);
127 m_tree->Branch("phiModule", &m_phiModule);
128 m_tree->Branch("layerWheel", &m_layerWheel);
129 m_tree->Branch("strawLayer", &m_strawLayer);
130 m_tree->Branch("straw", &m_straw);
131 m_tree->Branch("highLevel", &m_highLevel);
132 m_tree->Branch("timeOverThreshold", &m_timeOverThreshold);
133 m_tree->Branch("driftTimeBin", &m_driftTimeBin);
134 m_tree->Branch("trailEdge", &m_trailEdge);
135 m_tree->Branch("firstBin", &m_firstBin);
136 m_tree->Branch("lastBin", &m_lastBin);
137 // TRT SDO deposits
138 m_tree->Branch("sdoID", &m_sdoID);
139 m_tree->Branch("sdoWord", &m_sdoWord);
140 m_tree->Branch("barrelEndcap_sdo", &m_barrelEndcap_sdo);
141 m_tree->Branch("phiModule_sdo", &m_phiModule_sdo);
142 m_tree->Branch("layerWheel_sdo", &m_layerWheel_sdo);
143 m_tree->Branch("strawLayer_sdo", &m_strawLayer_sdo);
144 m_tree->Branch("straw_sdo", &m_straw_sdo);
145 m_tree->Branch("aboveThresh", &m_aboveThresh);
146 m_tree->Branch("deadChan", &m_deadChan);
147 m_tree->Branch("RODdata", &m_RODdata);
148 m_tree->Branch("validStraw", &m_validStraw);
149 m_tree->Branch("hit", &m_hit);
150 m_tree->Branch("barcode", &m_barcode);
151 m_tree->Branch("eventIndex", &m_eventIndex);
152 m_tree->Branch("charge", &m_charge);
153 m_tree->Branch("barcode_vec", &m_barcode_vec);
154 m_tree->Branch("eventIndex_vec", &m_eventIndex_vec);
155 m_tree->Branch("charge_vec", &m_charge_vec);
156 }
157 else {
158 ATH_MSG_ERROR("No tree found!");
159 }
160
161 // HISTOGRAMS
162 m_h_rdoID = new TH1F("h_rdoID", "rdoID", 100, 0, 2e18);
163 m_h_rdoID->StatOverflows();
164 ATH_CHECK(m_thistSvc->regHist(m_path + m_h_rdoID->GetName(), m_h_rdoID));
165
166 m_h_rdoWord = new TH1F("h_rdoWord", "rdoWord", 100, 0, 1.5e8);
167 m_h_rdoWord->StatOverflows();
168 ATH_CHECK(m_thistSvc->regHist(m_path + m_h_rdoWord->GetName(), m_h_rdoWord));
169
170 m_h_barrelEndcap = new TH1F("h_barrelEndcap", "Barrel or Endcap", 100, -3.5, 3.5);
171 m_h_barrelEndcap->StatOverflows();
172 ATH_CHECK(m_thistSvc->regHist(m_path + m_h_barrelEndcap->GetName(), m_h_barrelEndcap));
173
174 m_h_phiModule = new TH1F("h_phiModule", "Phi module", 100, 0, 35);
175 m_h_phiModule->StatOverflows();
176 ATH_CHECK(m_thistSvc->regHist(m_path + m_h_phiModule->GetName(), m_h_phiModule));
177
178 m_h_layerWheel = new TH1F("h_layerWheel", "Barrel layer or Endcap wheel", 100, 0, 20);
179 m_h_layerWheel->StatOverflows();
180 ATH_CHECK(m_thistSvc->regHist(m_path + m_h_layerWheel->GetName(), m_h_layerWheel));
181
182 m_h_strawLayer = new TH1F("h_strawLayer", "Straw layer", 100, 0, 35);
183 m_h_strawLayer->StatOverflows();
184 ATH_CHECK(m_thistSvc->regHist(m_path + m_h_strawLayer->GetName(), m_h_strawLayer));
185
186 m_h_straw = new TH1F("h_straw", "Straw", 100, 0, 1000);
187 m_h_straw->StatOverflows();
188 ATH_CHECK(m_thistSvc->regHist(m_path + m_h_straw->GetName(), m_h_straw));
189
190 m_h_ToT = new TH1F("h_ToT", "Time over Threshold", 100, 0, 100);
191 m_h_ToT->StatOverflows();
192 ATH_CHECK(m_thistSvc->regHist(m_path + m_h_ToT->GetName(), m_h_ToT));
193
194 m_h_ToT_HL = new TH1F("h_ToT_HL", "Time over Threshold (highLevel)", 100, 0, 100);
195 m_h_ToT_HL->StatOverflows();
196 ATH_CHECK(m_thistSvc->regHist(m_path + m_h_ToT_HL->GetName(), m_h_ToT_HL));
197
198 m_h_driftTimeBin = new TH1F("h_driftTimeBin", "Drift Time Bin", 100, 0, 30);
199 m_h_driftTimeBin->StatOverflows();
200 ATH_CHECK(m_thistSvc->regHist(m_path + m_h_driftTimeBin->GetName(), m_h_driftTimeBin));
201
202 m_h_trailEdge = new TH1F("h_trailEdge", "Trailing Edge", 100, 0, 30);
203 m_h_trailEdge->StatOverflows();
204 ATH_CHECK(m_thistSvc->regHist(m_path + m_h_trailEdge->GetName(), m_h_trailEdge));
205
206 m_h_brlPhiMod = new TH1F("h_brlPhiMod", "Barrel phi module", 100, 0, 35);
207 m_h_brlPhiMod->StatOverflows();
208 ATH_CHECK(m_thistSvc->regHist(m_path + m_h_brlPhiMod->GetName(), m_h_brlPhiMod));
209
210 m_h_brlLayer = new TH1F("h_brlLayer", "Barrel layer", 100, 0, 20);
211 m_h_brlLayer->StatOverflows();
212 ATH_CHECK(m_thistSvc->regHist(m_path + m_h_brlLayer->GetName(), m_h_brlLayer));
213
214 m_h_brlStrawLayer = new TH1F("h_brlStrawLayer", "Barrel straw layer", 100, 0, 35);
215 m_h_brlStrawLayer->StatOverflows();
216 ATH_CHECK(m_thistSvc->regHist(m_path + m_h_brlStrawLayer->GetName(), m_h_brlStrawLayer));
217
218 m_h_brlStraw = new TH1F("h_brlStraw", "Barrel straw", 100, 0, 1000);
219 m_h_brlStraw->StatOverflows();
220 ATH_CHECK(m_thistSvc->regHist(m_path + m_h_brlStraw->GetName(), m_h_brlStraw));
221
222 m_h_brlToT = new TH1F("h_brlToT", "Barrel Time over Threshold", 100, 0, 100);
223 m_h_brlToT->StatOverflows();
224 ATH_CHECK(m_thistSvc->regHist(m_path + m_h_brlToT->GetName(), m_h_brlToT));
225
226 m_h_brlToT_HL = new TH1F("h_brlToT_HL", "Barrel Time over Threshold (highLevel)", 100, 0, 100);
227 m_h_brlToT_HL->StatOverflows();
228 ATH_CHECK(m_thistSvc->regHist(m_path + m_h_brlToT_HL->GetName(), m_h_brlToT_HL));
229
230 m_h_brlDriftTimeBin = new TH1F("h_brlDriftTimeBin", "Barrel Drift Time Bin", 100, 0, 30);
231 m_h_brlDriftTimeBin->StatOverflows();
232 ATH_CHECK(m_thistSvc->regHist(m_path + m_h_brlDriftTimeBin->GetName(), m_h_brlDriftTimeBin));
233
234 m_h_brlTrailEdge = new TH1F("h_brlTrailEdge", "Barrel Trailing Edge", 100, 0, 30);
235 m_h_brlTrailEdge->StatOverflows();
236 ATH_CHECK(m_thistSvc->regHist(m_path + m_h_brlTrailEdge->GetName(), m_h_brlTrailEdge));
237
238 m_h_ecPhiMod = new TH1F("h_ecPhiMod", "Endcap phi module", 100, 0, 35);
239 m_h_ecPhiMod->StatOverflows();
240 ATH_CHECK(m_thistSvc->regHist(m_path + m_h_ecPhiMod->GetName(), m_h_ecPhiMod));
241
242 m_h_ecWheel = new TH1F("h_ecWheel", "Endcap wheel", 100, 0, 20);
243 m_h_ecWheel->StatOverflows();
244 ATH_CHECK(m_thistSvc->regHist(m_path + m_h_ecWheel->GetName(), m_h_ecWheel));
245
246 m_h_ecStrawLayer = new TH1F("h_ecStrawLayer", "Endcap straw layer", 100, 0, 35);
247 m_h_ecStrawLayer->StatOverflows();
248 ATH_CHECK(m_thistSvc->regHist(m_path + m_h_ecStrawLayer->GetName(), m_h_ecStrawLayer));
249
250 m_h_ecStraw = new TH1F("h_ecStraw", "Endcap straw", 100, 0, 1000);
251 m_h_ecStraw->StatOverflows();
252 ATH_CHECK(m_thistSvc->regHist(m_path + m_h_ecStraw->GetName(), m_h_ecStraw));
253
254 m_h_ecToT = new TH1F("h_ecToT", "Endcap Time over Threshold", 100, 0, 100);
255 m_h_ecToT->StatOverflows();
256 ATH_CHECK(m_thistSvc->regHist(m_path + m_h_ecToT->GetName(), m_h_ecToT));
257
258 m_h_ecToT_HL = new TH1F("h_ecToT_HL", "Endcap Time over Threshold (highLevel)", 100, 0, 100);
259 m_h_ecToT_HL->StatOverflows();
260 ATH_CHECK(m_thistSvc->regHist(m_path + m_h_ecToT_HL->GetName(), m_h_ecToT_HL));
261
262 m_h_ecDriftTimeBin = new TH1F("h_ecDriftTimeBin", "Endcap Drift Time Bin", 100, 0, 30);
263 m_h_ecDriftTimeBin->StatOverflows();
264 ATH_CHECK(m_thistSvc->regHist(m_path + m_h_ecDriftTimeBin->GetName(), m_h_ecDriftTimeBin));
265
266 m_h_ecTrailEdge = new TH1F("h_ecTrailEdge", "Endcap Trailing Edge", 100, 0, 30);
267 m_h_ecTrailEdge->StatOverflows();
268 ATH_CHECK(m_thistSvc->regHist(m_path + m_h_ecTrailEdge->GetName(), m_h_ecTrailEdge));
269
270 m_h_sdoID = new TH1F("h_sdoID", "sdoID", 100, 0, 2.5e18);
271 m_h_sdoID->StatOverflows();
272 ATH_CHECK(m_thistSvc->regHist(m_path + m_h_sdoID->GetName(), m_h_sdoID));
273
274 m_h_sdoWord = new TH1F("h_sdoWord", "sdoWord", 100, 0, 10);
275 m_h_sdoWord->StatOverflows();
276 ATH_CHECK(m_thistSvc->regHist(m_path + m_h_sdoWord->GetName(), m_h_sdoWord));
277
278 m_h_barrelEndcap_sdo = new TH1F("h_barrelEndcap_sdo", "Barrel or Endcap (sdo)", 100, -3.5, 3.5);
279 m_h_barrelEndcap_sdo->StatOverflows();
280 ATH_CHECK(m_thistSvc->regHist(m_path + m_h_barrelEndcap_sdo->GetName(), m_h_barrelEndcap_sdo));
281
282 m_h_phiModule_sdo = new TH1F("h_phiModule_sdo", "Phi module (sdo)", 100, 0, 35);
283 m_h_phiModule_sdo->StatOverflows();
284 ATH_CHECK(m_thistSvc->regHist(m_path + m_h_phiModule_sdo->GetName(), m_h_phiModule_sdo));
285
286 m_h_layerWheel_sdo = new TH1F("h_layerWheel_sdo", "Barrel layer or Endcap wheel (sdo)", 100, 0, 20);
287 m_h_layerWheel_sdo->StatOverflows();
288 ATH_CHECK(m_thistSvc->regHist(m_path + m_h_layerWheel_sdo->GetName(), m_h_layerWheel_sdo));
289
290 m_h_strawLayer_sdo = new TH1F("h_strawLayer_sdo", "Straw layer (sdo)", 100, 0, 35);
291 m_h_strawLayer_sdo->StatOverflows();
292 ATH_CHECK(m_thistSvc->regHist(m_path + m_h_strawLayer_sdo->GetName(), m_h_strawLayer_sdo));
293
294 m_h_straw_sdo = new TH1F("h_straw_sdo", "Straw (sdo)", 100, 0, 1000);
295 m_h_straw_sdo->StatOverflows();
296 ATH_CHECK(m_thistSvc->regHist(m_path + m_h_straw_sdo->GetName(), m_h_straw_sdo));
297
298 m_h_barcode = new TH1F("h_barcode", "Barcode (sdo)", 100, 0, 2.2e5);
299 m_h_barcode->StatOverflows();
300 ATH_CHECK(m_thistSvc->regHist(m_path + m_h_barcode->GetName(), m_h_barcode));
301
302 m_h_eventIndex = new TH1F("h_eventIndex", "Event index (sdo)", 100, 0, 10);
303 m_h_eventIndex->StatOverflows();
304 ATH_CHECK(m_thistSvc->regHist(m_path + m_h_eventIndex->GetName(), m_h_eventIndex));
305
306 m_h_charge = new TH1F("h_charge", "Charge (sdo)", 100, 0, 500);
307 m_h_charge->StatOverflows();
308 ATH_CHECK(m_thistSvc->regHist(m_path + m_h_charge->GetName(), m_h_charge));
309
310 return StatusCode::SUCCESS;
311}
312
313StatusCode TRT_RDOAnalysis::execute() {
314 ATH_MSG_DEBUG( "In TRT_RDOAnalysis::execute()" );
315
316 m_rdoID->clear();
317 m_rdoWord->clear();
318 m_barrelEndcap->clear();
319 m_phiModule->clear();
320 m_layerWheel->clear();
321 m_strawLayer->clear();
322 m_straw->clear();
323 m_highLevel->clear();
324 m_timeOverThreshold->clear();
325 m_driftTimeBin->clear();
326 m_trailEdge->clear();
327 m_firstBin->clear();
328 m_lastBin->clear();
329 m_sdoID->clear();
330 m_sdoWord->clear();
331 m_barrelEndcap_sdo->clear();
332 m_phiModule_sdo->clear();
333 m_layerWheel_sdo->clear();
334 m_strawLayer_sdo->clear();
335 m_straw_sdo->clear();
336 m_aboveThresh->clear();
337 m_deadChan->clear();
338 m_RODdata->clear();
339 m_validStraw->clear();
340 m_hit->clear();
341 m_barcode->clear();
342 m_eventIndex->clear();
343 m_charge->clear();
344 m_barcode_vec->clear();
345 m_eventIndex_vec->clear();
346 m_charge_vec->clear();
347
348 // RawData
349 SG::ReadHandle<TRT_RDO_Container> p_TRT_RDO_cont (m_inputKey);
350 if(p_TRT_RDO_cont.isValid()) {
351 // loop over RDO container
352 TRT_RDO_Container::const_iterator rdoCont_itr(p_TRT_RDO_cont->begin());
353 const TRT_RDO_Container::const_iterator rdoCont_end(p_TRT_RDO_cont->end());
354
355 for ( ; rdoCont_itr != rdoCont_end; ++rdoCont_itr ) {
356 const TRT_RDO_Collection* p_TRT_RDO_coll(*rdoCont_itr);
357 TRT_RDO_Collection::const_iterator rdo_itr(p_TRT_RDO_coll->begin());
358 const TRT_RDO_Collection::const_iterator rdo_end(p_TRT_RDO_coll->end());
359
360 for ( ; rdo_itr != rdo_end; ++rdo_itr ) {
361 const Identifier rdoID((*rdo_itr)->identify());
362 const unsigned int rdoWord((*rdo_itr)->getWord());
363 const int trtBrlEc(m_trtID->barrel_ec(rdoID));
364 const int trtPhiMod(m_trtID->phi_module(rdoID));
365 const int trtLayerWheel(m_trtID->layer_or_wheel(rdoID));
366 const int trtStrawLayer(m_trtID->straw_layer(rdoID));
367 const int trtStraw(m_trtID->straw(rdoID));
368 const bool trtHL((*rdo_itr)->highLevel());
369 const double trtToT((*rdo_itr)->timeOverThreshold());
370 const int trtDriftTimeBin((*rdo_itr)->driftTimeBin());
371 const int trtTrailEdge = static_cast<const TRT_LoLumRawData*>(*rdo_itr)->trailingEdge();
372 const bool trtFirstBin = static_cast<const TRT_LoLumRawData*>(*rdo_itr)->firstBinHigh();
373 const bool trtLastBin = static_cast<const TRT_LoLumRawData*>(*rdo_itr)->lastBinHigh();
374
375 const unsigned long long rdoID_int = rdoID.get_compact();
376 m_rdoID->push_back(rdoID_int);
377 m_rdoWord->push_back(rdoWord);
378 m_barrelEndcap->push_back(trtBrlEc);
379 m_phiModule->push_back(trtPhiMod);
380 m_layerWheel->push_back(trtLayerWheel);
381 m_strawLayer->push_back(trtStrawLayer);
382 m_straw->push_back(trtStraw);
383 m_highLevel->push_back(trtHL);
384 m_timeOverThreshold->push_back(trtToT);
385 m_driftTimeBin->push_back(trtDriftTimeBin);
386 m_trailEdge->push_back(trtTrailEdge);
387 m_firstBin->push_back(trtFirstBin);
388 m_lastBin->push_back(trtLastBin);
389
390 m_h_rdoID->Fill(rdoID_int);
391 m_h_rdoWord->Fill(rdoWord);
392 m_h_barrelEndcap->Fill(trtBrlEc);
393 m_h_phiModule->Fill(trtPhiMod);
394 m_h_layerWheel->Fill(trtLayerWheel);
395 m_h_strawLayer->Fill(trtStrawLayer);
396 m_h_straw->Fill(trtStraw);
397 m_h_ToT->Fill(trtToT);
398 if (trtHL) {
399 m_h_ToT_HL->Fill(trtToT);
400 }
401 m_h_driftTimeBin->Fill(trtDriftTimeBin);
402 m_h_trailEdge->Fill(trtTrailEdge);
403
404 if (std::abs(trtBrlEc) == 1) {
405 m_h_brlPhiMod->Fill(trtPhiMod);
406 m_h_brlLayer->Fill(trtLayerWheel);
407 m_h_brlStrawLayer->Fill(trtStrawLayer);
408 m_h_brlStraw->Fill(trtStraw);
409 m_h_brlToT->Fill(trtToT);
410 if (trtHL) {
411 m_h_brlToT_HL->Fill(trtToT);
412 }
413 m_h_brlDriftTimeBin->Fill(trtDriftTimeBin);
414 m_h_brlTrailEdge->Fill(trtTrailEdge);
415 }
416 else if (std::abs(trtBrlEc) == 2) {
417 m_h_ecPhiMod->Fill(trtPhiMod);
418 m_h_ecWheel->Fill(trtLayerWheel);
419 m_h_ecStrawLayer->Fill(trtStrawLayer);
420 m_h_ecStraw->Fill(trtStraw);
421 m_h_ecToT->Fill(trtToT);
422 if (trtHL) {
423 m_h_ecToT_HL->Fill(trtToT);
424 }
425 m_h_ecDriftTimeBin->Fill(trtDriftTimeBin);
426 m_h_ecTrailEdge->Fill(trtTrailEdge);
427 }
428 }
429 }
430 }
431
432 // SimData
433 SG::ReadHandle<InDetSimDataCollection> simDataMapTRT (m_inputTruthKey);
434 if(simDataMapTRT.isValid()) {
435 // loop over SDO container
436 InDetSimDataCollection::const_iterator sdo_itr(simDataMapTRT->begin());
437 const InDetSimDataCollection::const_iterator sdo_end(simDataMapTRT->end());
438
439 std::vector<int> barcode_vec;
440 std::vector<int> eventIndex_vec;
441 std::vector<float> charge_vec;
442 for ( ; sdo_itr != sdo_end; ++sdo_itr ) {
443 const Identifier sdoID((*sdo_itr).first);
444 const InDetSimData& sdo((*sdo_itr).second);
445
446 const unsigned long long sdoID_int = sdoID.get_compact();
447 const int sdoWord(sdo.word());
448 const int trtBrlEc_sdo(m_trtID->barrel_ec(sdoID));
449 const int trtPhiMod_sdo(m_trtID->phi_module(sdoID));
450 const int trtLayerWheel_sdo(m_trtID->layer_or_wheel(sdoID));
451 const int trtStrawLayer_sdo(m_trtID->straw_layer(sdoID));
452 const int trtStraw_sdo(m_trtID->straw(sdoID));
453 const bool aboveThresh(TRT_SimHelper::isAboveTR_Thresh(sdo));
454 const bool deadChan(TRT_SimHelper::isDeadChannel(sdo));
455 const bool rodData(TRT_SimHelper::isRodData(sdo));
456 const bool validStraw(TRT_SimHelper::isValidStraw(sdo));
457 const bool hit(TRT_SimHelper::HasHit(sdo));
458
459 m_sdoID->push_back(sdoID_int);
460 m_sdoWord->push_back(sdoWord);
461 m_barrelEndcap_sdo->push_back(trtBrlEc_sdo);
462 m_phiModule_sdo->push_back(trtPhiMod_sdo);
463 m_layerWheel_sdo->push_back(trtLayerWheel_sdo);
464 m_strawLayer_sdo->push_back(trtStrawLayer_sdo);
465 m_straw_sdo->push_back(trtStraw_sdo);
466 m_aboveThresh->push_back(aboveThresh);
467 m_deadChan->push_back(deadChan);
468 m_RODdata->push_back(rodData);
469 m_validStraw->push_back(validStraw);
470 m_hit->push_back(hit);
471
472 m_h_sdoID->Fill(sdoID_int);
473 m_h_sdoWord->Fill(sdoWord);
474 m_h_barrelEndcap_sdo->Fill(trtBrlEc_sdo);
475 m_h_phiModule_sdo->Fill(trtPhiMod_sdo);
476 m_h_layerWheel_sdo->Fill(trtLayerWheel_sdo);
477 m_h_strawLayer_sdo->Fill(trtStrawLayer_sdo);
478 m_h_straw_sdo->Fill(trtStraw_sdo);
479
480 // loop over deposits
481 // InDetSimData::Deposit typedef for std::pair<HepMcParticleLink, float>
482 const std::vector<InDetSimData::Deposit>& deposits = sdo.getdeposits();
483 std::vector<InDetSimData::Deposit>::const_iterator dep_itr(deposits.begin());
484 const std::vector<InDetSimData::Deposit>::const_iterator dep_end(deposits.end());
485
486 for ( ; dep_itr != dep_end; ++dep_itr ) {
487 const HepMcParticleLink& particleLink = (*dep_itr).first;
488 const int bar(HepMC::barcode(particleLink)); // FIXME barcode-based
489 const int eventIx(particleLink.eventIndex());
490 const float charge((*dep_itr).second);
491
492 m_barcode->push_back(bar);
493 m_eventIndex->push_back(eventIx);
494 m_charge->push_back(charge);
495
496 m_h_barcode->Fill(bar);
497 m_h_eventIndex->Fill(eventIx);
498 m_h_charge->Fill(charge);
499
500 barcode_vec.push_back(bar);
501 eventIndex_vec.push_back(eventIx);
502 charge_vec.push_back(charge);
503 }
504 m_barcode_vec->push_back(barcode_vec);
505 m_eventIndex_vec->push_back(eventIndex_vec);
506 m_charge_vec->push_back(charge_vec);
507 barcode_vec.clear();
508 eventIndex_vec.clear();
509 charge_vec.clear();
510 }
511 }
512
513 if (m_tree) {
514 m_tree->Fill();
515 }
516
517 return StatusCode::SUCCESS;
518}
519
520StatusCode TRT_RDOAnalysis::finalize() {
521
522 return StatusCode::SUCCESS;
523
524}
ATH_CHECK
#define ATH_CHECK
Evaluate an expression and check for errors.
Definition AthCheckMacros.h:40
ATH_MSG_ERROR
#define ATH_MSG_ERROR(x)
Definition AthMsgStreamMacros.h:33
ATH_MSG_DEBUG
#define ATH_MSG_DEBUG(x)
Definition AthMsgStreamMacros.h:29
charge
double charge(const T &p)
Definition AtlasPID.h:997
ReadHandle.h
Handle class for reading from StoreGate.
TRT_RDOAnalysis.h
TRT_RDO_Collection
InDetRawDataCollection< TRT_RDORawData > TRT_RDO_Collection
Definition TRT_RDO_Collection.h:23
AthAlgorithm::AthAlgorithm
AthAlgorithm(const std::string &name, ISvcLocator *pSvcLocator)
Constructor with parameters:
Definition AthAlgorithm.cxx:24
AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T, V, H > &t)
Definition AthCommonDataStore.h:145
AthCommonDataStore< AthCommonMsg< Algorithm > >::detStore
const ServiceHandle< StoreGateSvc > & detStore() const
Definition AthCommonDataStore.h:95
InDetRawDataCollection< TRT_RDORawData >::const_iterator
DataModel_detail::const_iterator< DataVector > const_iterator
Definition DataVector.h:838
DataVector< RawDataT >::end
const_iterator end() const noexcept
DataVector< RawDataT >::begin
const_iterator begin() const noexcept
IdentifiableContainerMT::const_iterator
Definition IdentifiableContainerMT.h:82
Identifier::get_compact
value_type get_compact() const
Get the compact id.
InDetSimData::word
int word() const
Definition InDetSimData.h:69
InDetSimData::getdeposits
const std::vector< Deposit > & getdeposits() const
Definition InDetSimData.h:74
SG::ReadHandle::isValid
virtual bool isValid() override final
Can the handle be successfully dereferenced?
TRT_LoLumRawData
Definition TRT_LoLumRawData.h:25
TRT_LoLumRawData::firstBinHigh
bool firstBinHigh() const
Definition TRT_LoLumRawData.h:154
TRT_LoLumRawData::lastBinHigh
bool lastBinHigh() const
Definition TRT_LoLumRawData.h:165
TRT_LoLumRawData::trailingEdge
int trailingEdge() const
Definition TRT_LoLumRawData.h:66
TRT_RDOAnalysis::m_sdoID
std::vector< unsigned long long > * m_sdoID
Definition TRT_RDOAnalysis.h:62
TRT_RDOAnalysis::m_strawLayer
std::vector< int > * m_strawLayer
Definition TRT_RDOAnalysis.h:52
TRT_RDOAnalysis::m_h_eventIndex
TH1 * m_h_eventIndex
Definition TRT_RDOAnalysis.h:123
TRT_RDOAnalysis::m_RODdata
std::vector< bool > * m_RODdata
Definition TRT_RDOAnalysis.h:73
TRT_RDOAnalysis::m_h_ecToT
TH1 * m_h_ecToT
Definition TRT_RDOAnalysis.h:110
TRT_RDOAnalysis::m_h_brlPhiMod
TH1 * m_h_brlPhiMod
Definition TRT_RDOAnalysis.h:97
TRT_RDOAnalysis::finalize
virtual StatusCode finalize() override final
Definition TRT_RDOAnalysis.cxx:520
TRT_RDOAnalysis::m_ntupleFileName
std::string m_ntupleFileName
Definition TRT_RDOAnalysis.h:127
TRT_RDOAnalysis::m_ntupleTreeName
std::string m_ntupleTreeName
Definition TRT_RDOAnalysis.h:129
TRT_RDOAnalysis::m_h_trailEdge
TH1 * m_h_trailEdge
Definition TRT_RDOAnalysis.h:95
TRT_RDOAnalysis::m_h_ecWheel
TH1 * m_h_ecWheel
Definition TRT_RDOAnalysis.h:107
TRT_RDOAnalysis::m_h_straw
TH1 * m_h_straw
Definition TRT_RDOAnalysis.h:91
TRT_RDOAnalysis::m_h_brlToT
TH1 * m_h_brlToT
Definition TRT_RDOAnalysis.h:101
TRT_RDOAnalysis::m_h_brlStraw
TH1 * m_h_brlStraw
Definition TRT_RDOAnalysis.h:100
TRT_RDOAnalysis::m_aboveThresh
std::vector< bool > * m_aboveThresh
Definition TRT_RDOAnalysis.h:71
TRT_RDOAnalysis::m_highLevel
std::vector< bool > * m_highLevel
Definition TRT_RDOAnalysis.h:55
TRT_RDOAnalysis::m_phiModule_sdo
std::vector< int > * m_phiModule_sdo
Definition TRT_RDOAnalysis.h:66
TRT_RDOAnalysis::m_h_ecToT_HL
TH1 * m_h_ecToT_HL
Definition TRT_RDOAnalysis.h:111
TRT_RDOAnalysis::m_ntupleDirName
std::string m_ntupleDirName
Definition TRT_RDOAnalysis.h:128
TRT_RDOAnalysis::m_h_ecTrailEdge
TH1 * m_h_ecTrailEdge
Definition TRT_RDOAnalysis.h:113
TRT_RDOAnalysis::m_sdoWord
std::vector< int > * m_sdoWord
Definition TRT_RDOAnalysis.h:63
TRT_RDOAnalysis::m_h_ecPhiMod
TH1 * m_h_ecPhiMod
Definition TRT_RDOAnalysis.h:106
TRT_RDOAnalysis::m_trtID
const TRT_ID * m_trtID
Definition TRT_RDOAnalysis.h:42
TRT_RDOAnalysis::m_tree
TTree * m_tree
Definition TRT_RDOAnalysis.h:126
TRT_RDOAnalysis::m_rdoID
std::vector< unsigned long long > * m_rdoID
Definition TRT_RDOAnalysis.h:46
TRT_RDOAnalysis::m_h_phiModule_sdo
TH1 * m_h_phiModule_sdo
Definition TRT_RDOAnalysis.h:118
TRT_RDOAnalysis::m_h_layerWheel
TH1 * m_h_layerWheel
Definition TRT_RDOAnalysis.h:89
TRT_RDOAnalysis::m_inputTruthKey
SG::ReadHandleKey< InDetSimDataCollection > m_inputTruthKey
Definition TRT_RDOAnalysis.h:41
TRT_RDOAnalysis::m_eventIndex_vec
std::vector< std::vector< int > > * m_eventIndex_vec
Definition TRT_RDOAnalysis.h:81
TRT_RDOAnalysis::m_inputKey
SG::ReadHandleKey< TRT_RDO_Container > m_inputKey
Definition TRT_RDOAnalysis.h:40
TRT_RDOAnalysis::m_h_ToT_HL
TH1 * m_h_ToT_HL
Definition TRT_RDOAnalysis.h:93
TRT_RDOAnalysis::m_charge_vec
std::vector< std::vector< float > > * m_charge_vec
Definition TRT_RDOAnalysis.h:82
TRT_RDOAnalysis::m_hit
std::vector< bool > * m_hit
Definition TRT_RDOAnalysis.h:75
TRT_RDOAnalysis::m_driftTimeBin
std::vector< int > * m_driftTimeBin
Definition TRT_RDOAnalysis.h:57
TRT_RDOAnalysis::m_barrelEndcap
std::vector< int > * m_barrelEndcap
Definition TRT_RDOAnalysis.h:49
TRT_RDOAnalysis::m_timeOverThreshold
std::vector< double > * m_timeOverThreshold
Definition TRT_RDOAnalysis.h:56
TRT_RDOAnalysis::m_eventIndex
std::vector< int > * m_eventIndex
Definition TRT_RDOAnalysis.h:78
TRT_RDOAnalysis::m_straw_sdo
std::vector< int > * m_straw_sdo
Definition TRT_RDOAnalysis.h:69
TRT_RDOAnalysis::m_h_brlTrailEdge
TH1 * m_h_brlTrailEdge
Definition TRT_RDOAnalysis.h:104
TRT_RDOAnalysis::m_h_ecDriftTimeBin
TH1 * m_h_ecDriftTimeBin
Definition TRT_RDOAnalysis.h:112
TRT_RDOAnalysis::m_h_sdoID
TH1 * m_h_sdoID
Definition TRT_RDOAnalysis.h:115
TRT_RDOAnalysis::m_h_charge
TH1 * m_h_charge
Definition TRT_RDOAnalysis.h:124
TRT_RDOAnalysis::m_h_barcode
TH1 * m_h_barcode
Definition TRT_RDOAnalysis.h:122
TRT_RDOAnalysis::m_h_rdoWord
TH1 * m_h_rdoWord
Definition TRT_RDOAnalysis.h:86
TRT_RDOAnalysis::m_deadChan
std::vector< bool > * m_deadChan
Definition TRT_RDOAnalysis.h:72
TRT_RDOAnalysis::m_validStraw
std::vector< bool > * m_validStraw
Definition TRT_RDOAnalysis.h:74
TRT_RDOAnalysis::m_barcode_vec
std::vector< std::vector< int > > * m_barcode_vec
Definition TRT_RDOAnalysis.h:80
TRT_RDOAnalysis::m_h_barrelEndcap
TH1 * m_h_barrelEndcap
Definition TRT_RDOAnalysis.h:87
TRT_RDOAnalysis::m_h_strawLayer
TH1 * m_h_strawLayer
Definition TRT_RDOAnalysis.h:90
TRT_RDOAnalysis::m_h_ecStraw
TH1 * m_h_ecStraw
Definition TRT_RDOAnalysis.h:109
TRT_RDOAnalysis::m_straw
std::vector< int > * m_straw
Definition TRT_RDOAnalysis.h:53
TRT_RDOAnalysis::m_path
std::string m_path
Definition TRT_RDOAnalysis.h:130
TRT_RDOAnalysis::m_rdoWord
std::vector< unsigned int > * m_rdoWord
Definition TRT_RDOAnalysis.h:47
TRT_RDOAnalysis::m_h_brlStrawLayer
TH1 * m_h_brlStrawLayer
Definition TRT_RDOAnalysis.h:99
TRT_RDOAnalysis::m_barrelEndcap_sdo
std::vector< int > * m_barrelEndcap_sdo
Definition TRT_RDOAnalysis.h:65
TRT_RDOAnalysis::m_h_driftTimeBin
TH1 * m_h_driftTimeBin
Definition TRT_RDOAnalysis.h:94
TRT_RDOAnalysis::m_h_brlDriftTimeBin
TH1 * m_h_brlDriftTimeBin
Definition TRT_RDOAnalysis.h:103
TRT_RDOAnalysis::m_barcode
std::vector< int > * m_barcode
Definition TRT_RDOAnalysis.h:77
TRT_RDOAnalysis::m_phiModule
std::vector< int > * m_phiModule
Definition TRT_RDOAnalysis.h:50
TRT_RDOAnalysis::m_h_sdoWord
TH1 * m_h_sdoWord
Definition TRT_RDOAnalysis.h:116
TRT_RDOAnalysis::m_strawLayer_sdo
std::vector< int > * m_strawLayer_sdo
Definition TRT_RDOAnalysis.h:68
TRT_RDOAnalysis::m_h_ecStrawLayer
TH1 * m_h_ecStrawLayer
Definition TRT_RDOAnalysis.h:108
TRT_RDOAnalysis::TRT_RDOAnalysis
TRT_RDOAnalysis(const std::string &name, ISvcLocator *pSvcLocator)
Definition TRT_RDOAnalysis.cxx:14
TRT_RDOAnalysis::m_h_ToT
TH1 * m_h_ToT
Definition TRT_RDOAnalysis.h:92
TRT_RDOAnalysis::m_layerWheel
std::vector< int > * m_layerWheel
Definition TRT_RDOAnalysis.h:51
TRT_RDOAnalysis::execute
virtual StatusCode execute() override final
Definition TRT_RDOAnalysis.cxx:313
TRT_RDOAnalysis::m_h_brlToT_HL
TH1 * m_h_brlToT_HL
Definition TRT_RDOAnalysis.h:102
TRT_RDOAnalysis::m_thistSvc
ServiceHandle< ITHistSvc > m_thistSvc
Definition TRT_RDOAnalysis.h:131
TRT_RDOAnalysis::m_layerWheel_sdo
std::vector< int > * m_layerWheel_sdo
Definition TRT_RDOAnalysis.h:67
TRT_RDOAnalysis::m_h_barrelEndcap_sdo
TH1 * m_h_barrelEndcap_sdo
Definition TRT_RDOAnalysis.h:117
TRT_RDOAnalysis::m_charge
std::vector< float > * m_charge
Definition TRT_RDOAnalysis.h:79
TRT_RDOAnalysis::m_h_phiModule
TH1 * m_h_phiModule
Definition TRT_RDOAnalysis.h:88
TRT_RDOAnalysis::m_h_straw_sdo
TH1 * m_h_straw_sdo
Definition TRT_RDOAnalysis.h:121
TRT_RDOAnalysis::m_firstBin
std::vector< bool > * m_firstBin
Definition TRT_RDOAnalysis.h:59
TRT_RDOAnalysis::m_h_strawLayer_sdo
TH1 * m_h_strawLayer_sdo
Definition TRT_RDOAnalysis.h:120
TRT_RDOAnalysis::m_h_brlLayer
TH1 * m_h_brlLayer
Definition TRT_RDOAnalysis.h:98
TRT_RDOAnalysis::m_trailEdge
std::vector< int > * m_trailEdge
Definition TRT_RDOAnalysis.h:58
TRT_RDOAnalysis::initialize
virtual StatusCode initialize() override final
Definition TRT_RDOAnalysis.cxx:104
TRT_RDOAnalysis::m_h_rdoID
TH1 * m_h_rdoID
Definition TRT_RDOAnalysis.h:85
TRT_RDOAnalysis::m_lastBin
std::vector< bool > * m_lastBin
Definition TRT_RDOAnalysis.h:60
TRT_RDOAnalysis::m_h_layerWheel_sdo
TH1 * m_h_layerWheel_sdo
Definition TRT_RDOAnalysis.h:119
TRT_SimHelper::HasHit
static bool HasHit(const InDetSimData &sdo)
Definition TRT_SimHelper.h:152
TRT_SimHelper::isDeadChannel
static bool isDeadChannel(const InDetSimData &sdo)
Definition TRT_SimHelper.h:138
TRT_SimHelper::isAboveTR_Thresh
static bool isAboveTR_Thresh(const InDetSimData &sdo)
Definition TRT_SimHelper.h:133
TRT_SimHelper::isValidStraw
static bool isValidStraw(const InDetSimData &sdo)
Definition TRT_SimHelper.h:148
TRT_SimHelper::isRodData
static bool isRodData(const InDetSimData &sdo)
Definition TRT_SimHelper.h:143
HepMC::barcode
int barcode(const T *p)
Definition Barcode.h:16
Generated on for ATLAS Offline Software by doxygen 1.14.0