29 auto props = this->getProperties();
30 for( Gaudi::Details::PropertyBase* prop : props ) {
31 if (prop->ownerTypeName()==this->type()) {
32 ATH_MSG_DEBUG(
"Property:\t" << prop->name() <<
"\t : \t" << prop->toString());
40 return StatusCode::SUCCESS;
68 "RZ_allhits",
"; Z [mm] ; R [mm] ", 500, -2000, 2000, 500, 0, 500));
71 "phiResidual",
"phi residual [mm]", 1000, -10, 10));
73 "etaResidual",
"eta residual [mm]", 1000, -10, 10));
76 "phiTrueBinShift",
"phi TrueBinShift [mm]", 1000, -10, 10));
78 "etaTrueBinShift",
"eta TrueBinShift [mm]", 1000, -10, 10));
85 "; Hits per bin in step", 1000, 0, 1000));
93 ";#Hit Layers per Bin", nLyrs + 1, 0, nLyrs + 1));
97 "; Layer ; Hits ", nLyrs, 0, nLyrs, 20, 0,
111 "; Layer; Hits in Bin Passing Threshold",
112 nLyrs, 0, nLyrs, 20, 0, 20));
134 "PairSetMatchPhi",
";PairSet Match Phi [mm]", 500,
138 "PairSetMatchEta",
";PairSet Match Eta [mm]", 500,
150 ";Phi Curvature ", 5000,
155 ";Eta Curvature ", 5000,
160 ";Delta Phi Curvature ", 5000,
165 ";Delta Eta Curvature ", 500,
170 ";Phi Out Extrap Curved", 2000, -1.0*
m_phiScale,
174 ";Phi In Extrap Curved ", 2000, -1.0*
m_phiScale,
177 ";PairSet Match Phi vs dRin [mm]", 100, 0, 100, 500,
180 ";PairSet Match Eta vs dRin [mm]", 100, 0, 100, 500,
189 return StatusCode::SUCCESS;
209 return StatusCode::SUCCESS;
239 std::vector<FPGATrackSimGenScanMonitoring::eventDispSet *> allsets(
245 for (
auto&
set : allsets) {
249 return StatusCode::SUCCESS;
261 for (
auto& hit :
data.hits) {
270 std::vector<FPGATrackSimBinUtil::StoredHit> sorted_hits =
data.hits;
271 std::sort(sorted_hits.begin(), sorted_hits.end(),
272 [](
const auto &hit1,
const auto &hit2) {
273 return hit1.rzrad() < hit2.rzrad();
279 for (
auto &hit : sorted_hits) {
280 m_tree_r.push_back(hit.hitptr->getR());
281 m_tree_z.push_back(hit.hitptr->getZ());
282 m_tree_id.push_back(hit.hitptr->getIdentifier());
283 m_tree_hash.push_back(hit.hitptr->getIdentifierHash());
297 const std::vector<std::vector<const FPGATrackSimBinUtil::StoredHit *> > & hitsByLayer)
331 bool passedPairFilter) {
337 if (passedPairFilter) {
346 << filteredpairs.
lyrCnt() <<
" "
347 << std::bitset<16>(filteredpairs.
hitLayers));
354 const std::vector<std::shared_ptr<const FPGATrackSimHit>> &
hits)
366 unsigned cnt =
bin.data().hitsInLyr(lyr);
375 std::vector<FPGATrackSimTruthTrack>
const * truthtracks,
bool isSingleParticle)
378 <<
" size = " << (truthtracks ? truthtracks->size() : 0));
406 <<
" ?= " << recovered <<
" closure:"
439 std::vector<FPGATrackSimBinUtil::IdxSet> idxsets =
441 std::vector<unsigned>({0, 1, 2, 3, 4}),
453 std::vector<const FPGATrackSimBinUtil::StoredHit *>
const *lastlyr,
454 std::vector<const FPGATrackSimBinUtil::StoredHit *>
const *lastlastlyr) {
456 auto size_if_nonzero_ptr =
457 [](std::vector<const FPGATrackSimBinUtil::StoredHit *>
const *
ptr) {
466 size_if_nonzero_ptr(lastlastlyr));
473 int lyr = pair.first->hitptr->getLayer();
495 unsigned minlyr =
std::min(lastpair.first->layer, lastpair.second->layer);
496 minlyr =
std::min(minlyr, pair.first->layer);
497 minlyr =
std::min(minlyr, pair.second->layer);
500 ((1 << lastpair.first->layer) | (1 << lastpair.second->layer) |
501 (1 << pair.first->layer) | (1 << pair.second->layer));
502 hitbits = hitbits >> minlyr;
504 int hitlyrs = std::popcount(hitbits);
506 assert(hitbits & 0
x1);
512 }
else if (hitlyrs == 3) {
513 if (hitbits == 0b111)
515 if ((lastpair.first->layer == pair.second->layer) ||
516 (pair.first->layer == lastpair.second->layer)) {
521 }
else if ((hitbits == 0b1101) || (hitbits == 0b1011))
531 <<
" " << lastpair.second->layer <<
" "
532 << pair.first->layer <<
" " << pair.second->layer <<
" : "
533 << minlyr <<
" " << std::bitset<16>(hitbits));
539 std::vector<TH1D *> & histset,
double val,
545 histset[0]->Fill(
val);
546 histset[ppcat]->Fill(
val);
554 const std::vector<FPGATrackSimGenScanTool::HitPairSet> &pairsets,
558 for (
const auto &pairset : pairsets) {
570 const std::vector<FPGATrackSimGenScanTool::IntermediateState> &
states,
571 unsigned allowed_misses) {
573 for (
unsigned i = 0;
i <
states.size();
i++) {
576 for (
auto &pair :
states[
i].pairsets) {
579 unsigned totalInput =
589 const std::string &
name) {
590 TGraph *retv =
new TGraph();
591 retv->SetName(
name.c_str());
596 TGraph *
g,
double x,
double y) {
597 g->SetPoint(
g->GetN(),
x,
y);
602 unsigned int count = m_rZ.size();
605 if (
count > m_maxEvts)
611 for (
auto &hit :
hits) {
612 AddPoint(m_rZ.back(), hit->getZ(), hit->getR());
613 AddPoint(m_xY.back(), hit->getX(), hit->getY());
619 if (m_rZ.size() > m_maxEvts)
return;
622 std::vector<const FPGATrackSimBinUtil::StoredHit *> ptrvec;
623 for (
auto &hit :
hits) {
624 ptrvec.push_back(&hit);
628 return addEvent(ptrvec);
632 unsigned int count = m_rZ.size();
635 if (
count > m_maxEvts)
return;
643 for (
auto &hit :
hits) {
644 AddPoint(m_rZ.back(), hit->hitptr->getZ(), hit->hitptr->getR());
645 AddPoint(m_xY.back(), hit->hitptr->getX(), hit->hitptr->getY());
646 AddPoint(m_rEta.back(), hit->etaShift, hit->hitptr->getR());
647 AddPoint(m_rPhi.back(), hit->phiShift, hit->hitptr->getR());
654 for (
auto &
graphs : std::vector<std::vector<TGraph *> *>({&m_rZ, &m_xY, &m_rPhi, &m_rEta})) {
659 return StatusCode::SUCCESS;