LArSoft  v09_90_00
Liquid Argon Software toolkit - https://larsoft.org/
CheatingBeamParticleIdTool.cc
Go to the documentation of this file.
1 
9 #include "Pandora/AlgorithmHeaders.h"
10 
12 
15 
16 using namespace pandora;
17 
18 namespace lar_content
19 {
20 
21 CheatingBeamParticleIdTool::CheatingBeamParticleIdTool() :
22  m_minWeightFraction(0.5f)
23 {
24 }
25 
26 //------------------------------------------------------------------------------------------------------------------------------------------
27 
28 void CheatingBeamParticleIdTool::SelectOutputPfos(const pandora::Algorithm *const /*pAlgorithm*/,
29  const SliceHypotheses &testBeamSliceHypotheses, const SliceHypotheses &crSliceHypotheses, PfoList &selectedPfos)
30 {
31  if (testBeamSliceHypotheses.size() != crSliceHypotheses.size())
32  throw StatusCodeException(STATUS_CODE_INVALID_PARAMETER);
33 
34  for (unsigned int sliceIndex = 0, nSlices = testBeamSliceHypotheses.size(); sliceIndex < nSlices; ++sliceIndex)
35  {
36  float beamParticleWeight(0.f), totalWeight(0.f);
37  const PfoList &testBeamPfoList(testBeamSliceHypotheses.at(sliceIndex));
38 
39  for (const Pfo *const pTestBeamPfo : testBeamPfoList)
40  {
41  if (!LArPfoHelper::IsTestBeam(pTestBeamPfo))
42  throw StatusCodeException(STATUS_CODE_INVALID_PARAMETER);
43 
44  PfoList downstreamPfos;
45  LArPfoHelper::GetAllDownstreamPfos(pTestBeamPfo, downstreamPfos);
46 
47  float thisBeamParticleWeight(0.f), thisTotalWeight(0.f);
48  CheatingSliceIdBaseTool::GetTargetParticleWeight(&downstreamPfos, thisBeamParticleWeight, thisTotalWeight, LArMCParticleHelper::IsBeamParticle);
49 
50  beamParticleWeight += thisBeamParticleWeight;
51  totalWeight += thisTotalWeight;
52  }
53 
54  const float beamWeightFraction(totalWeight < std::numeric_limits<float>::epsilon() ? 0.f : beamParticleWeight / totalWeight);
55 
56  if (beamWeightFraction > m_minWeightFraction)
57  {
58  const PfoList &sliceOutput(testBeamSliceHypotheses.at(sliceIndex));
59  selectedPfos.insert(selectedPfos.end(), sliceOutput.begin(), sliceOutput.end());
60  }
61  else
62  {
63  const PfoList &sliceOutput(crSliceHypotheses.at(sliceIndex));
64  selectedPfos.insert(selectedPfos.end(), sliceOutput.begin(), sliceOutput.end());
65  }
66  }
67 }
68 
69 //------------------------------------------------------------------------------------------------------------------------------------------
70 
71 StatusCode CheatingBeamParticleIdTool::ReadSettings(const TiXmlHandle xmlHandle)
72 {
73  PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, XmlHelper::ReadValue(xmlHandle, "MinimumWeightFraction", m_minWeightFraction));
74 
75  return STATUS_CODE_SUCCESS;
76 }
77 
78 } // namespace lar_content
pandora::StatusCode ReadSettings(const pandora::TiXmlHandle xmlHandle)
std::vector< pandora::PfoList > SliceHypotheses
Header file for the pfo helper class.
TFile f
Definition: plotHisto.C:6
void SelectOutputPfos(const pandora::Algorithm *const pAlgorithm, const SliceHypotheses &nuSliceHypotheses, const SliceHypotheses &crSliceHypotheses, pandora::PfoList &selectedPfos)
Select which reconstruction hypotheses to use; neutrino outcomes or cosmic-ray muon outcomes for each...
Header file for the cheating beam particle id tool class.
static bool IsTestBeam(const pandora::ParticleFlowObject *const pPfo)
Whether a pfo is a test beam particle.
Header file for the lar monte carlo particle helper helper class.
static bool IsBeamParticle(const pandora::MCParticle *const pMCParticle)
Returns true if passed a primary beam MCParticle.
float m_minWeightFraction
The minimum weight fraction for identifying a slice as a beam particle.
static void GetTargetParticleWeight(const pandora::PfoList *const pPfoList, float &targetParticleWeight, float &totalWeight, std::function< bool(const pandora::MCParticle *const)> fCriteria)
Get the target particle weight in a list of pfos.
static void GetAllDownstreamPfos(const pandora::PfoList &inputPfoList, pandora::PfoList &outputPfoList)
Get a flat list of all pfos, recursively, of all daughters associated with those pfos in an input lis...