9 #include "Pandora/AlgorithmHeaders.h" 24 TestBeamEventValidationAlgorithm::TestBeamEventValidationAlgorithm()
30 TestBeamEventValidationAlgorithm::~TestBeamEventValidationAlgorithm()
36 void TestBeamEventValidationAlgorithm::FillValidationInfo(
const MCParticleList *
const pMCParticleList,
37 const CaloHitList *
const pCaloHitList,
const PfoList *
const pPfoList,
ValidationInfo &validationInfo)
const 39 if (pMCParticleList && pCaloHitList)
42 LArMCParticleHelper::SelectReconstructableMCParticles(
43 pMCParticleList, pCaloHitList, m_primaryParameters, LArMCParticleHelper::IsBeamParticle, targetMCParticleToHitsMap);
44 LArMCParticleHelper::SelectReconstructableMCParticles(
45 pMCParticleList, pCaloHitList, m_primaryParameters, LArMCParticleHelper::IsCosmicRay, targetMCParticleToHitsMap);
52 LArMCParticleHelper::SelectReconstructableMCParticles(
53 pMCParticleList, pCaloHitList, parameters, LArMCParticleHelper::IsBeamParticle, allMCParticleToHitsMap);
54 LArMCParticleHelper::SelectReconstructableMCParticles(
55 pMCParticleList, pCaloHitList, parameters, LArMCParticleHelper::IsCosmicRay, allMCParticleToHitsMap);
63 PfoList allConnectedPfos;
64 LArPfoHelper::GetAllConnectedPfos(*pPfoList, allConnectedPfos);
66 PfoList finalStatePfos;
67 for (
const ParticleFlowObject *
const pPfo : allConnectedPfos)
69 if (pPfo->GetParentPfoList().empty())
70 finalStatePfos.push_back(pPfo);
74 LArMCParticleHelper::GetPfoToReconstructable2DHitsMap(
82 LArMCParticleHelper::GetPfoMCParticleHitSharingMaps(
87 this->InterpretMatching(validationInfo, interpretedMCToPfoHitSharingMap);
93 void TestBeamEventValidationAlgorithm::ProcessOutput(
94 const ValidationInfo &validationInfo,
const bool useInterpretedMatching,
const bool printToScreen,
const bool fillTree)
const 96 static int eventNumber{-1};
98 if (printToScreen && useInterpretedMatching)
99 std::cout <<
"---INTERPRETED-MATCHING-OUTPUT------------------------------------------------------------------" << std::endl;
100 else if (printToScreen)
101 std::cout <<
"---RAW-MATCHING-OUTPUT--------------------------------------------------------------------------" << std::endl;
110 PfoVector primaryPfoVector;
111 LArMonitoringHelper::GetOrderedPfoVector(validationInfo.
GetPfoToHitsMap(), primaryPfoVector);
113 int pfoIndex(0), testBeamPfoIndex(0);
116 for (
const Pfo *
const pPrimaryPfo : primaryPfoVector)
118 pfoToIdMap.insert(PfoToIdMap::value_type(pPrimaryPfo, ++pfoIndex));
119 const Pfo *
const pRecoTestBeam(LArPfoHelper::IsTestBeamFinalState(pPrimaryPfo) ? LArPfoHelper::GetParentPfo(pPrimaryPfo) :
nullptr);
121 if (pRecoTestBeam && !testBeamPfoToIdMap.count(pRecoTestBeam))
122 testBeamPfoToIdMap.insert(PfoToIdMap::value_type(pRecoTestBeam, ++testBeamPfoIndex));
127 MCParticleList associatedMCPrimaries;
129 int nCorrectTB(0), nTotalTB(0), nCorrectCR(0), nTotalCR(0), nFakeTB(0), nFakeCR(0), nSplitTB(0), nSplitCR(0), nLost(0);
130 int mcPrimaryIndex(0), nTargetMatches(0), nTargetTBMatches(0), nTargetCRMatches(0);
131 IntVector mcPrimaryId, mcPrimaryPdg, nMCHitsTotal, nMCHitsU, nMCHitsV, nMCHitsW;
132 FloatVector mcPrimaryE, mcPrimaryPX, mcPrimaryPY, mcPrimaryPZ;
133 FloatVector mcPrimaryVtxX, mcPrimaryVtxY, mcPrimaryVtxZ, mcPrimaryEndX, mcPrimaryEndY, mcPrimaryEndZ;
134 IntVector nPrimaryMatchedPfos, nPrimaryMatchedTBPfos, nPrimaryMatchedCRPfos;
135 IntVector bestMatchPfoId, bestMatchPfoPdg, bestMatchPfoIsTB;
136 IntVector bestMatchPfoNHitsTotal, bestMatchPfoNHitsU, bestMatchPfoNHitsV, bestMatchPfoNHitsW;
137 IntVector bestMatchPfoNSharedHitsTotal, bestMatchPfoNSharedHitsU, bestMatchPfoNSharedHitsV, bestMatchPfoNSharedHitsW;
138 FloatVector bestMatchPfoX0;
140 std::stringstream targetSS;
142 for (
const MCParticle *
const pMCPrimary : mcPrimaryVector)
144 const bool hasMatch(mcToPfoHitSharingMap.count(pMCPrimary) && !mcToPfoHitSharingMap.at(pMCPrimary).empty());
147 if (!isTargetPrimary && !hasMatch)
150 associatedMCPrimaries.push_back(pMCPrimary);
154 const int mcNuanceCode(LArMCParticleHelper::GetNuanceCode(LArMCParticleHelper::GetParentMCParticle(pMCPrimary)));
155 const int isBeamParticle(LArMCParticleHelper::IsBeamParticle(pMCPrimary));
156 const int isCosmicRay(LArMCParticleHelper::IsCosmicRay(pMCPrimary));
158 const int nTargetPrimaries(associatedMCPrimaries.size());
159 const CartesianVector &targetVertex(LArMCParticleHelper::GetParentMCParticle(pMCPrimary)->GetVertex());
160 const float targetVertexX(targetVertex.GetX()), targetVertexY(targetVertex.GetY()), targetVertexZ(targetVertex.GetZ());
163 targetSS << (!isTargetPrimary ?
"(Non target) " :
"") <<
"PrimaryId " << mcPrimaryIndex <<
", TB " << isBeamParticle <<
", CR " 164 << isCosmicRay <<
", MCPDG " << pMCPrimary->GetParticleId() <<
", Energy " << pMCPrimary->GetEnergy() <<
", Dist. " 165 << (pMCPrimary->GetEndpoint() - pMCPrimary->GetVertex()).GetMagnitude() <<
", nMCHits " << mcPrimaryHitList.size() <<
" (" 166 << LArMonitoringHelper::CountHitsByType(TPC_VIEW_U, mcPrimaryHitList) <<
", " 167 << LArMonitoringHelper::CountHitsByType(TPC_VIEW_V, mcPrimaryHitList) <<
", " 168 << LArMonitoringHelper::CountHitsByType(TPC_VIEW_W, mcPrimaryHitList) <<
")" << std::endl;
170 mcPrimaryId.push_back(mcPrimaryIndex);
171 mcPrimaryPdg.push_back(pMCPrimary->GetParticleId());
172 mcPrimaryE.push_back(pMCPrimary->GetEnergy());
173 mcPrimaryPX.push_back(pMCPrimary->GetMomentum().GetX());
174 mcPrimaryPY.push_back(pMCPrimary->GetMomentum().GetY());
175 mcPrimaryPZ.push_back(pMCPrimary->GetMomentum().GetZ());
176 mcPrimaryVtxX.push_back(pMCPrimary->GetVertex().GetX());
177 mcPrimaryVtxY.push_back(pMCPrimary->GetVertex().GetY());
178 mcPrimaryVtxZ.push_back(pMCPrimary->GetVertex().GetZ());
179 mcPrimaryEndX.push_back(pMCPrimary->GetEndpoint().GetX());
180 mcPrimaryEndY.push_back(pMCPrimary->GetEndpoint().GetY());
181 mcPrimaryEndZ.push_back(pMCPrimary->GetEndpoint().GetZ());
182 nMCHitsTotal.push_back(mcPrimaryHitList.size());
183 nMCHitsU.push_back(LArMonitoringHelper::CountHitsByType(TPC_VIEW_U, mcPrimaryHitList));
184 nMCHitsV.push_back(LArMonitoringHelper::CountHitsByType(TPC_VIEW_V, mcPrimaryHitList));
185 nMCHitsW.push_back(LArMonitoringHelper::CountHitsByType(TPC_VIEW_W, mcPrimaryHitList));
187 int matchIndex(0), nPrimaryMatches(0), nPrimaryTBMatches(0), nPrimaryCRMatches(0);
189 float recoVertexX(std::numeric_limits<float>::max()), recoVertexY(std::numeric_limits<float>::max()),
190 recoVertexZ(std::numeric_limits<float>::max());
194 const CaloHitList &sharedHitList(pfoToSharedHits.second);
195 const CaloHitList &pfoHitList(validationInfo.
GetPfoToHitsMap().at(pfoToSharedHits.first));
197 const bool isRecoTestBeam(LArPfoHelper::IsTestBeam(pfoToSharedHits.first));
198 const bool isGoodMatch(this->IsGoodMatch(mcPrimaryHitList, pfoHitList, sharedHitList));
200 const int pfoId(pfoToIdMap.at(pfoToSharedHits.first));
202 if (0 == matchIndex++)
204 bestMatchPfoId.push_back(pfoId);
205 bestMatchPfoPdg.push_back(pfoToSharedHits.first->GetParticleId());
206 bestMatchPfoIsTB.push_back(isRecoTestBeam ? 1 : 0);
207 bestMatchPfoNHitsTotal.push_back(pfoHitList.size());
208 bestMatchPfoNHitsU.push_back(LArMonitoringHelper::CountHitsByType(TPC_VIEW_U, pfoHitList));
209 bestMatchPfoNHitsV.push_back(LArMonitoringHelper::CountHitsByType(TPC_VIEW_V, pfoHitList));
210 bestMatchPfoNHitsW.push_back(LArMonitoringHelper::CountHitsByType(TPC_VIEW_W, pfoHitList));
211 bestMatchPfoNSharedHitsTotal.push_back(sharedHitList.size());
212 bestMatchPfoNSharedHitsU.push_back(LArMonitoringHelper::CountHitsByType(TPC_VIEW_U, sharedHitList));
213 bestMatchPfoNSharedHitsV.push_back(LArMonitoringHelper::CountHitsByType(TPC_VIEW_V, sharedHitList));
214 bestMatchPfoNSharedHitsW.push_back(LArMonitoringHelper::CountHitsByType(TPC_VIEW_W, sharedHitList));
215 bestMatchPfoX0.push_back(pfoToSharedHits.first->GetPropertiesMap().count(
"X0") ? pfoToSharedHits.first->GetPropertiesMap().at(
"X0")
216 : std::numeric_limits<float>::max());
220 const Vertex *
const pRecoVertex(isRecoTestBeam ? LArPfoHelper::GetTestBeamInteractionVertex(pfoToSharedHits.first)
221 : LArPfoHelper::GetVertex(pfoToSharedHits.first));
222 recoVertexX = pRecoVertex->GetPosition().GetX();
223 recoVertexY = pRecoVertex->GetPosition().GetY();
224 recoVertexZ = pRecoVertex->GetPosition().GetZ();
226 catch (
const StatusCodeException &)
235 if (isRecoTestBeam && isGoodMatch)
237 if (!isRecoTestBeam && isGoodMatch)
240 targetSS <<
"-" << (!isGoodMatch ?
"(Below threshold) " :
"") <<
"MatchedPfoId " << pfoId <<
", TB " << isRecoTestBeam
241 <<
", CR " << (!isRecoTestBeam) <<
", PDG " << pfoToSharedHits.first->GetParticleId() <<
", nMatchedHits " 242 << sharedHitList.size() <<
" (" << LArMonitoringHelper::CountHitsByType(TPC_VIEW_U, sharedHitList) <<
", " 243 << LArMonitoringHelper::CountHitsByType(TPC_VIEW_V, sharedHitList) <<
", " 244 << LArMonitoringHelper::CountHitsByType(TPC_VIEW_W, sharedHitList) <<
")" 245 <<
", nPfoHits " << pfoHitList.size() <<
" (" << LArMonitoringHelper::CountHitsByType(TPC_VIEW_U, pfoHitList) <<
", " 246 << LArMonitoringHelper::CountHitsByType(TPC_VIEW_V, pfoHitList) <<
", " 247 << LArMonitoringHelper::CountHitsByType(TPC_VIEW_W, pfoHitList) <<
")" << std::endl;
250 if (mcToPfoHitSharingMap.at(pMCPrimary).empty())
252 targetSS <<
"-No matched Pfo" << std::endl;
253 bestMatchPfoId.push_back(-1);
254 bestMatchPfoPdg.push_back(0);
255 bestMatchPfoIsTB.push_back(0);
256 bestMatchPfoNHitsTotal.push_back(0);
257 bestMatchPfoNHitsU.push_back(0);
258 bestMatchPfoNHitsV.push_back(0);
259 bestMatchPfoNHitsW.push_back(0);
260 bestMatchPfoNSharedHitsTotal.push_back(0);
261 bestMatchPfoNSharedHitsU.push_back(0);
262 bestMatchPfoNSharedHitsV.push_back(0);
263 bestMatchPfoNSharedHitsW.push_back(0);
264 bestMatchPfoX0.push_back(std::numeric_limits<float>::max());
267 nPrimaryMatchedPfos.push_back(nPrimaryMatches);
268 nPrimaryMatchedTBPfos.push_back(nPrimaryTBMatches);
269 nPrimaryMatchedCRPfos.push_back(nPrimaryCRMatches);
270 nTargetMatches += nPrimaryMatches;
271 nTargetTBMatches += nPrimaryTBMatches;
272 nTargetCRMatches += nPrimaryCRMatches;
276 PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), m_treeName.c_str(),
"fileIdentifier", m_fileIdentifier));
277 PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), m_treeName.c_str(),
"eventNumber", eventNumber));
278 PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), m_treeName.c_str(),
"mcNuanceCode", mcNuanceCode));
279 PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), m_treeName.c_str(),
"isBeamParticle", isBeamParticle));
280 PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), m_treeName.c_str(),
"isCosmicRay", isCosmicRay));
281 PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), m_treeName.c_str(),
"nTargetPrimaries", nTargetPrimaries));
282 PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), m_treeName.c_str(),
"targetVertexX", targetVertexX));
283 PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), m_treeName.c_str(),
"targetVertexY", targetVertexY));
284 PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), m_treeName.c_str(),
"targetVertexZ", targetVertexZ));
285 PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), m_treeName.c_str(),
"recoVertexX", recoVertexX));
286 PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), m_treeName.c_str(),
"recoVertexY", recoVertexY));
287 PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), m_treeName.c_str(),
"recoVertexZ", recoVertexZ));
288 PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), m_treeName.c_str(),
"mcPrimaryId", &mcPrimaryId));
289 PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), m_treeName.c_str(),
"mcPrimaryPdg", &mcPrimaryPdg));
290 PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), m_treeName.c_str(),
"mcPrimaryE", &mcPrimaryE));
291 PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), m_treeName.c_str(),
"mcPrimaryPX", &mcPrimaryPX));
292 PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), m_treeName.c_str(),
"mcPrimaryPY", &mcPrimaryPY));
293 PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), m_treeName.c_str(),
"mcPrimaryPZ", &mcPrimaryPZ));
294 PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), m_treeName.c_str(),
"mcPrimaryVtxX", &mcPrimaryVtxX));
295 PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), m_treeName.c_str(),
"mcPrimaryVtxY", &mcPrimaryVtxY));
296 PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), m_treeName.c_str(),
"mcPrimaryVtxZ", &mcPrimaryVtxZ));
297 PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), m_treeName.c_str(),
"mcPrimaryEndX", &mcPrimaryEndX));
298 PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), m_treeName.c_str(),
"mcPrimaryEndY", &mcPrimaryEndY));
299 PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), m_treeName.c_str(),
"mcPrimaryEndZ", &mcPrimaryEndZ));
300 PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), m_treeName.c_str(),
"mcPrimaryNHitsTotal", &nMCHitsTotal));
301 PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), m_treeName.c_str(),
"mcPrimaryNHitsU", &nMCHitsU));
302 PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), m_treeName.c_str(),
"mcPrimaryNHitsV", &nMCHitsV));
303 PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), m_treeName.c_str(),
"mcPrimaryNHitsW", &nMCHitsW));
304 PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), m_treeName.c_str(),
"nPrimaryMatchedPfos", &nPrimaryMatchedPfos));
305 PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), m_treeName.c_str(),
"nPrimaryMatchedTBPfos", &nPrimaryMatchedTBPfos));
306 PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), m_treeName.c_str(),
"nPrimaryMatchedCRPfos", &nPrimaryMatchedCRPfos));
307 PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), m_treeName.c_str(),
"bestMatchPfoId", &bestMatchPfoId));
308 PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), m_treeName.c_str(),
"bestMatchPfoPdg", &bestMatchPfoPdg));
309 PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), m_treeName.c_str(),
"bestMatchPfoNHitsTotal", &bestMatchPfoNHitsTotal));
310 PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), m_treeName.c_str(),
"bestMatchPfoNHitsU", &bestMatchPfoNHitsU));
311 PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), m_treeName.c_str(),
"bestMatchPfoNHitsV", &bestMatchPfoNHitsV));
312 PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), m_treeName.c_str(),
"bestMatchPfoNHitsW", &bestMatchPfoNHitsW));
313 PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), m_treeName.c_str(),
"bestMatchPfoNSharedHitsTotal", &bestMatchPfoNSharedHitsTotal));
314 PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), m_treeName.c_str(),
"bestMatchPfoNSharedHitsU", &bestMatchPfoNSharedHitsU));
315 PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), m_treeName.c_str(),
"bestMatchPfoNSharedHitsV", &bestMatchPfoNSharedHitsV));
316 PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), m_treeName.c_str(),
"bestMatchPfoNSharedHitsW", &bestMatchPfoNSharedHitsW));
317 PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), m_treeName.c_str(),
"bestMatchPfoX0", &bestMatchPfoX0));
318 PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), m_treeName.c_str(),
"nTargetMatches", nTargetMatches));
319 PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), m_treeName.c_str(),
"nTargetTBMatches", nTargetTBMatches));
320 PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), m_treeName.c_str(),
"nTargetCRMatches", nTargetCRMatches));
321 PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), m_treeName.c_str(),
"bestMatchPfoIsTB", &bestMatchPfoIsTB));
324 if (isBeamParticle || isCosmicRay)
328 const int interactionTypeInt(static_cast<int>(interactionType));
331 const int isCorrectTB(isBeamParticle && (nTargetTBMatches == 1) && (nTargetCRMatches == 0));
332 const int isCorrectCR(isCosmicRay && (nTargetTBMatches == 0) && (nTargetCRMatches == 1));
333 const int isFakeTB(isCosmicRay && (nTargetTBMatches > 0));
334 const int isFakeCR(!isCosmicRay && (nTargetCRMatches > 0));
335 const int isSplitTB(!isCosmicRay && (nTargetTBMatches > 1));
336 const int isSplitCR(isCosmicRay && (nTargetCRMatches > 1));
337 const int isLost(nTargetMatches == 0);
339 std::stringstream outcomeSS;
340 outcomeSS << LArInteractionTypeHelper::ToString(interactionType) <<
" (Nuance " << mcNuanceCode <<
", TB " << isBeamParticle
341 <<
", CR " << isCosmicRay <<
")" << std::endl;
363 outcomeSS <<
"IsCorrectTB ";
365 outcomeSS <<
"IsCorrectCR ";
367 outcomeSS <<
"IsFakeTB ";
369 outcomeSS <<
"IsFakeCR ";
371 outcomeSS <<
"isSplitTB ";
373 outcomeSS <<
"IsSplitCR ";
375 outcomeSS <<
"IsLost ";
376 if (nTargetTBMatches > 0)
377 outcomeSS <<
"(NTBMatches: " << nTargetTBMatches <<
") ";
378 if (nTargetCRMatches > 0)
379 outcomeSS <<
"(NCRMatches: " << nTargetCRMatches <<
") ";
381 std::cout << outcomeSS.str() << std::endl << targetSS.str() << std::endl;
385 PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), m_treeName.c_str(),
"interactionType", interactionTypeInt));
386 PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), m_treeName.c_str(),
"isCorrectTB", isCorrectTB));
387 PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), m_treeName.c_str(),
"isCorrectCR", isCorrectCR));
388 PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), m_treeName.c_str(),
"isFakeTB", isFakeTB));
389 PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), m_treeName.c_str(),
"isFakeCR", isFakeCR));
390 PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), m_treeName.c_str(),
"isSplitTB", isSplitTB));
391 PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), m_treeName.c_str(),
"isSplitCR", isSplitCR));
392 PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), m_treeName.c_str(),
"isLost", isLost));
393 PANDORA_MONITORING_API(FillTree(this->GetPandora(), m_treeName.c_str()));
396 targetSS.str(std::string());
398 associatedMCPrimaries.clear();
400 nTargetTBMatches = 0;
401 nTargetCRMatches = 0;
403 mcPrimaryPdg.clear();
404 nMCHitsTotal.clear();
412 mcPrimaryVtxX.clear();
413 mcPrimaryVtxY.clear();
414 mcPrimaryVtxZ.clear();
415 mcPrimaryEndX.clear();
416 mcPrimaryEndY.clear();
417 mcPrimaryEndZ.clear();
418 nPrimaryMatchedPfos.clear();
419 nPrimaryMatchedTBPfos.clear();
420 nPrimaryMatchedCRPfos.clear();
421 bestMatchPfoId.clear();
422 bestMatchPfoPdg.clear();
423 bestMatchPfoIsTB.clear();
424 bestMatchPfoNHitsTotal.clear();
425 bestMatchPfoNHitsU.clear();
426 bestMatchPfoNHitsV.clear();
427 bestMatchPfoNHitsW.clear();
428 bestMatchPfoNSharedHitsTotal.clear();
429 bestMatchPfoNSharedHitsU.clear();
430 bestMatchPfoNSharedHitsV.clear();
431 bestMatchPfoNSharedHitsW.clear();
432 bestMatchPfoX0.clear();
436 if (useInterpretedMatching)
438 std::stringstream summarySS;
439 summarySS <<
"---SUMMARY--------------------------------------------------------------------------------------" << std::endl;
441 summarySS <<
"#CorrectTB: " << nCorrectTB <<
"/" << nTotalTB
442 <<
", Fraction: " << (nTotalTB > 0 ?
static_cast<float>(nCorrectTB) / static_cast<float>(nTotalTB) : 0.f) << std::endl;
444 summarySS <<
"#CorrectCR: " << nCorrectCR <<
"/" << nTotalCR
445 <<
", Fraction: " << (nTotalCR > 0 ?
static_cast<float>(nCorrectCR) / static_cast<float>(nTotalCR) : 0.f) << std::endl;
447 summarySS <<
"#FakeTB: " << nFakeTB <<
" ";
449 summarySS <<
"#FakeCR: " << nFakeCR <<
" ";
451 summarySS <<
"#SplitTB: " << nSplitTB <<
" ";
453 summarySS <<
"#SplitCR: " << nSplitCR <<
" ";
455 summarySS <<
"#Lost: " << nLost <<
" ";
456 if (nFakeTB || nFakeCR || nSplitTB || nSplitCR || nLost)
457 summarySS << std::endl;
459 std::cout << summarySS.str();
463 std::cout <<
"------------------------------------------------------------------------------------------------" << std::endl
469 StatusCode TestBeamEventValidationAlgorithm::ReadSettings(
const TiXmlHandle xmlHandle)
471 return EventValidationBaseAlgorithm::ReadSettings(xmlHandle);
std::unordered_map< const pandora::ParticleFlowObject *, unsigned int > PfoToIdMap
const LArMCParticleHelper::PfoContributionMap & GetPfoToHitsMap() const
Get the pfo to hits map.
Header file for the pfo helper class.
std::unordered_map< const pandora::MCParticle *, pandora::CaloHitList > MCContributionMap
const LArMCParticleHelper::MCContributionMap & GetTargetMCParticleToHitsMap() const
Get the target mc particle to hits map.
Header file for the interaction type helper class.
unsigned int m_minPrimaryGoodHits
the minimum number of primary good Hits
const LArMCParticleHelper::MCContributionMap & GetAllMCParticleToHitsMap() const
Get the all mc particle to hits map.
std::map< const pandora::MCParticle *, PfoToSharedHitsVector > MCParticleToPfoHitSharingMap
const LArMCParticleHelper::MCParticleToPfoHitSharingMap & GetMCToPfoHitSharingMap() const
Get the mc to pfo hit sharing map.
void SetInterpretedMCToPfoHitSharingMap(const LArMCParticleHelper::MCParticleToPfoHitSharingMap &interpretedMCToPfoHitSharingMap)
Set the interpreted mc to pfo hit sharing map.
std::vector< int > IntVector
Header file for the lar monitoring helper helper class.
void SetMCToPfoHitSharingMap(const LArMCParticleHelper::MCParticleToPfoHitSharingMap &mcToPfoHitSharingMap)
Set the mc to pfo hit sharing map.
unsigned int m_minHitsForGoodView
the minimum number of Hits for a good view
void SetTargetMCParticleToHitsMap(const LArMCParticleHelper::MCContributionMap &targetMCParticleToHitsMap)
Set the target mc particle to hits map.
std::unordered_map< const pandora::MCParticle *, int > MCParticleIntMap
std::vector< art::Ptr< simb::MCParticle > > MCParticleVector
const LArMCParticleHelper::MCParticleToPfoHitSharingMap & GetInterpretedMCToPfoHitSharingMap() const
Get the interpreted mc to pfo hit sharing map.
float m_minHitSharingFraction
the minimum Hit sharing fraction
InteractionType
InteractionType enum.
void SetAllMCParticleToHitsMap(const LArMCParticleHelper::MCContributionMap &allMCParticleToHitsMap)
Set the all mc particle to hits map.
std::pair< const pandora::ParticleFlowObject *, pandora::CaloHitList > PfoCaloHitListPair
void SetPfoToHitsMap(const LArMCParticleHelper::PfoContributionMap &pfoToHitsMap)
Set the pfo to hits map.
Header file for the test beam event validation algorithm.
boost::graph_traits< ModuleGraph >::vertex_descriptor Vertex
std::unordered_map< const pandora::ParticleFlowObject *, pandora::CaloHitList > PfoContributionMap
std::map< const pandora::ParticleFlowObject *, MCParticleToSharedHitsVector > PfoToMCParticleHitSharingMap