ATLAS Offline Software
Functions
FPGATrackSimLayerStudyTool.cxx File Reference
#include "FPGATrackSimBinning/FPGATrackSimLayerStudyTool.h"
#include "AthenaBaseComps/AthMsgStreamMacros.h"
#include "FPGATrackSimBinning/IFPGATrackSimBinDesc.h"
#include "FPGATrackSimBinning/FPGATrackSimBinUtil.h"
#include "FPGATrackSimBinning/FPGATrackSimBinnedHits.h"
#include "FPGATrackSimObjects/FPGATrackSimTrackPars.h"
#include "TH1D.h"
#include "TH2D.h"
#include "TTree.h"
#include <algorithm>
#include <bit>
#include <numbers>

Go to the source code of this file.

Functions

void FPGATrackSimLayerStudyTool::fillBinLevelOutput ATLAS_NOT_THREAD_SAFE (const FPGATrackSimBinUtil::IdxSet &idx, const FPGATrackSimBinnedHits::BinEntry &data)
 
void FPGATrackSimLayerStudyTool::fillBinningSummary ATLAS_NOT_THREAD_SAFE (const std::vector< std::shared_ptr< const FPGATrackSimHit >> &hits)
 
void FPGATrackSimLayerStudyTool::parseTruthInfo ATLAS_NOT_THREAD_SAFE (std::vector< FPGATrackSimTruthTrack > const &truthtracks)
 

Function Documentation

◆ ATLAS_NOT_THREAD_SAFE() [1/3]

void FPGATrackSimLayerStudyTool::fillBinLevelOutput ATLAS_NOT_THREAD_SAFE ( const FPGATrackSimBinUtil::IdxSet idx,
const FPGATrackSimBinnedHits::BinEntry data 
)

Definition at line 179 of file FPGATrackSimLayerStudyTool.cxx.

181 {
182  setBinPlotsActive(idx);
183 
184  // fill all truth
185  m_ptDist[0]->Fill(std::abs(1/m_truthpars.qOverPt));
186  m_etaDist[0]->Fill(1/m_truthpars.eta);
187  m_phiDist[0]->Fill(1/m_truthpars.phi);
188  m_d0Dist[0]->Fill(1/m_truthpars.d0);
189  m_z0Dist[0]->Fill(1/m_truthpars.z0);
190 
191  if (m_binPlotsActive) {
192  for (auto& hit : data.hits) {
193  m_phiShift_road->Fill(hit.phiShift);
194  m_etaShift_road->Fill(hit.etaShift);
195  m_phiShift2D_road->Fill(hit.phiShift, hit.hitptr->getR());
196  m_etaShift2D_road->Fill(hit.etaShift, hit.hitptr->getR());
197  }
198 
199  // fill param monitoring
200  for (int i = 0; i < 2; i++) {
201  // i=0 all, i=1 no missed layers, i=2 is one missed layer
202  if ((i == 0)|| (data.lyrCnt() >= m_binnedhits->getNLayers() - (i - 1))) {
203  // all layer hit
204  m_ptDist[i]->Fill(std::abs(1 / m_truthpars.qOverPt));
205  m_etaDist[i]->Fill(1 / m_truthpars.eta);
206  m_phiDist[i]->Fill(1 / m_truthpars.phi);
207  m_d0Dist[i]->Fill(1 / m_truthpars.d0);
208  m_z0Dist[i]->Fill(1 / m_truthpars.z0);
209  }
210  }
211 
212  // Module mapping and Layer definition studies
213  // first sort hits by r+z radii
214  std::vector<FPGATrackSimBinUtil::StoredHit> sorted_hits = data.hits;
215  std::sort(sorted_hits.begin(), sorted_hits.end(),
216  [](const auto &hit1, const auto &hit2) {
217  return hit1.rzrad() < hit2.rzrad();
218  });
219 
220  // Fill tree
221  m_bin_tree_bin = std::vector<unsigned>(idx);
222  ClearTreeVectors();
223  for (auto &hit : sorted_hits) {
224 
225  m_bin_tree_r.push_back(hit.hitptr->getR());
226  m_bin_tree_z.push_back(hit.hitptr->getZ());
227  m_bin_tree_id.push_back(hit.hitptr->getIdentifier());
228  m_bin_tree_hash.push_back(hit.hitptr->getIdentifierHash());
229  m_bin_tree_layer.push_back(hit.hitptr->getLayerDisk());
230  m_bin_tree_side.push_back(hit.hitptr->getSide());
231  m_bin_tree_etamod.push_back(hit.hitptr->getEtaModule());
232  m_bin_tree_phimod.push_back(hit.hitptr->getPhiModule());
233  m_bin_tree_dettype.push_back((int)hit.hitptr->getDetType());
234  m_bin_tree_detzone.push_back((int)hit.hitptr->getDetectorZone());
235  }
236  m_bin_tree->Fill();
237  m_binsFilledCnt++;
238  }
239 }

◆ ATLAS_NOT_THREAD_SAFE() [2/3]

void FPGATrackSimLayerStudyTool::fillBinningSummary ATLAS_NOT_THREAD_SAFE ( const std::vector< std::shared_ptr< const FPGATrackSimHit >> &  hits)

Definition at line 241 of file FPGATrackSimLayerStudyTool.cxx.

243 {
244  m_inputHits->Fill(hits.size());
245 
246  for (auto &step : m_binnedhits->getBinTool().steps()) {
247  for (auto bin : m_binnedhits->binnedHits()[step->stepNum()])
248  m_hitsPerStepBin[step->stepNum()]->Fill(bin.data().hitCnt);
249  }
250 
251  for (auto bin :m_binnedhits->binnedHits()[m_binnedhits->getBinTool().lastStep()->stepNum()]) {
252  for (unsigned lyr = 0; lyr < m_binnedhits->getNLayers(); lyr++) {
253  unsigned cnt = bin.data().hitsInLyr(lyr);
254  m_hitsPerLayer->Fill(lyr, cnt);
255  m_hitsPerLayer2D->Fill(lyr, cnt);
256  }
257  }
258 
259  m_binsFilled->Fill(m_binsFilledCnt);
260  m_binsFilledCnt=0;
261 
262 }

◆ ATLAS_NOT_THREAD_SAFE() [3/3]

void FPGATrackSimLayerStudyTool::parseTruthInfo ATLAS_NOT_THREAD_SAFE ( std::vector< FPGATrackSimTruthTrack > const truthtracks)

Definition at line 293 of file FPGATrackSimLayerStudyTool.cxx.

293  {
294  ATH_MSG_DEBUG("In parseTruthInfo, truthtracks size = " << truthtracks.size());
295  m_truthIsValid = false;
296 
297  const IFPGATrackSimBinDesc* bindesc = m_binnedhits->getBinTool().binDesc();
298 
299  if (truthtracks.size() == 0) return;
300 
301  // Convert to binning parameters and find truth bin for each step
302  m_truthpars = (truthtracks)[0].getPars();
303  m_truthpars[FPGATrackSimTrackPars::IHIP] =
304  m_truthpars[FPGATrackSimTrackPars::IHIP] * 1000;
305  m_truthparset = bindesc->trackParsToParSet(m_truthpars);
306  m_truthbin.clear();
307  for (auto &step : m_binnedhits->getBinTool().steps()) {
308  m_truthbin.push_back(step->binIdx(m_truthparset));
309  }
310 
311  // histogram parameters
312  for (unsigned i = 0; i < FPGATrackSimTrackPars::NPARS; i++) {
313  m_truthpars_hists[i]->Fill(m_truthparset[i]);
314  }
315 
316  // a closure test
317  FPGATrackSimTrackPars recovered = bindesc->parSetToTrackPars(m_truthparset);
318  ATH_MSG_DEBUG("parset:" << m_truthparset << " " << m_truthpars
319  << " ?= " << recovered << " closure:"
320  << " " << recovered[FPGATrackSimTrackPars::IHIP] - m_truthpars[FPGATrackSimTrackPars::IHIP]
321  << " " << recovered[FPGATrackSimTrackPars::IPHI] - m_truthpars[FPGATrackSimTrackPars::IPHI]
322  << " " << recovered[FPGATrackSimTrackPars::ID0] - m_truthpars[FPGATrackSimTrackPars::ID0]
323  << " " << recovered[FPGATrackSimTrackPars::IETA] - m_truthpars[FPGATrackSimTrackPars::IETA]
324  << " " << recovered[FPGATrackSimTrackPars::IZ0] - m_truthpars[FPGATrackSimTrackPars::IZ0]);
325 
326  m_truth_tree_phi = m_truthpars[FPGATrackSimTrackPars::IPHI];
327  m_truth_tree_qOverPt = m_truthpars[FPGATrackSimTrackPars::IHIP];
328  m_truth_tree_d0 = m_truthpars[FPGATrackSimTrackPars::ID0];
329  m_truth_tree_z0 = m_truthpars[FPGATrackSimTrackPars::IZ0];
330  m_truth_tree_eta = m_truthpars[FPGATrackSimTrackPars::IETA];
331  m_truth_tree_parset = std::vector<double>(m_truthparset);
332  m_truth_tree->Fill();
333 
334  // print if there are multiple tracks for debugging single track MC
335  if (truthtracks.size() > 1) {
336  for (unsigned i = 0; i < truthtracks.size(); i++) {
337  ATH_MSG_INFO("Multiple truth" << i << " of " << truthtracks.size()
338  << " " << (truthtracks)[i].getPars());
339  }
340  }
341 
342  // find truth bin for later plotting selections
343  ATH_MSG_DEBUG("truthbin " << truthtracks.size()
344  << " " << m_truthpars << " " << m_truthbin);
345 
346  // Check if the truth track falls in the binning range
347  if (!m_binnedhits->getBinTool().inRange(m_truthparset)) {
348  ATH_MSG_INFO("Truth out of range because truth parset = " << m_truthparset << " wrt min = " << m_binnedhits->getBinTool().parMin() << ", max = " << m_binnedhits->getBinTool().parMax());
349  return;
350  }
351 
352  // Check that truth track falls in an actual bin
353  // this should alway pass, except for weird events
354  m_truthIsValid = true;
355  for (auto &step : m_binnedhits->getBinTool().steps()) {
356  if (!step->validBinsFull()[m_truthbin[step->stepNum()]]) {
357  ATH_MSG_INFO("Truth Bin not valid! Step " << step->stepName() << " "
358  << m_truthbin[step->stepNum()]
359  << " : " << m_truthpars);
360  std::vector<FPGATrackSimBinUtil::IdxSet> idxsets =
362  std::vector<unsigned>({0, 1, 2, 3, 4}),
363  m_truthbin[step->stepNum()]);
364  for (FPGATrackSimBinUtil::IdxSet &idxset : idxsets) {
365  ATH_MSG_INFO("Truth Box "
366  << bindesc->parSetToTrackPars(step->binLowEdge(idxset)));
367  }
368  m_truthIsValid = false;
369  }
370  }
371 }
FPGATrackSimTrackPars::IHIP
@ IHIP
Definition: FPGATrackSimTrackPars.h:49
data
char data[hepevt_bytes_allocation_ATLAS]
Definition: HepEvt.cxx:11
FPGATrackSimTrackPars::ID0
@ ID0
Definition: FPGATrackSimTrackPars.h:49
TRTCalib_Extractor.hits
hits
Definition: TRTCalib_Extractor.py:35
ATH_MSG_INFO
#define ATH_MSG_INFO(x)
Definition: AthMsgStreamMacros.h:31
FPGATrackSimTrackPars
Definition: FPGATrackSimTrackPars.h:22
IFPGATrackSimBinDesc::trackParsToParSet
virtual const FPGATrackSimBinUtil::ParSet trackParsToParSet(const FPGATrackSimTrackPars &pars) const =0
FPGATrackSimTrackPars::IZ0
@ IZ0
Definition: FPGATrackSimTrackPars.h:49
bin
Definition: BinsDiffFromStripMedian.h:43
lumiFormat.i
int i
Definition: lumiFormat.py:85
ATH_MSG_DEBUG
#define ATH_MSG_DEBUG(x)
Definition: AthMsgStreamMacros.h:29
IFPGATrackSimBinDesc::parSetToTrackPars
virtual const FPGATrackSimTrackPars parSetToTrackPars(const FPGATrackSimBinUtil::ParSet &parset) const =0
FPGATrackSimBinUtil::makeVariationSet
std::vector< IdxSet > makeVariationSet(const std::vector< unsigned > &scanpars, const IdxSet &idx)
Definition: FPGATrackSimBinUtil.cxx:93
FPGATrackSimTrackPars::IPHI
@ IPHI
Definition: FPGATrackSimTrackPars.h:49
trigbs_pickEvents.cnt
cnt
Definition: trigbs_pickEvents.py:71
std::sort
void sort(typename std::reverse_iterator< DataModel_detail::iterator< DVL > > beg, typename std::reverse_iterator< DataModel_detail::iterator< DVL > > end, const Compare &comp)
Specialization of sort for DataVector/List.
Definition: DVL_algorithms.h:623
FPGATrackSimTrackPars::NPARS
@ NPARS
Definition: FPGATrackSimTrackPars.h:49
IFPGATrackSimBinDesc
Definition: IFPGATrackSimBinDesc.h:47
LArNewCalib_DelayDump_OFC_Cali.idx
idx
Definition: LArNewCalib_DelayDump_OFC_Cali.py:69
LArCellBinning.step
step
Definition: LArCellBinning.py:158
FPGATrackSimBinUtil::IdxSet
Definition: FPGATrackSimBinUtil.h:50
FPGATrackSimTrackPars::IETA
@ IETA
Definition: FPGATrackSimTrackPars.h:49