9 #include "Helpers/MCParticleHelper.h"
10 #include "Pandora/AlgorithmHeaders.h"
21 StatusCode CheatingCosmicRayRemovalAlgorithm::Run()
23 const MCParticleList *pMCParticleList(
nullptr);
24 PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::GetList(*
this, m_mcParticleListName, pMCParticleList));
26 const CaloHitList *pCaloHitList(
nullptr);
27 PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::GetList(*
this, m_inputCaloHitListName, pCaloHitList));
29 CaloHitList outputCaloHitList;
31 for (
const CaloHit *pCaloHit : *pCaloHitList)
35 const MCParticle *
const pMCParticle(MCParticleHelper::GetMainMCParticle(pCaloHit));
36 if (!LArMCParticleHelper::IsCosmicRay(LArMCParticleHelper::GetParentMCParticle(pMCParticle)))
37 outputCaloHitList.push_back(pCaloHit);
39 catch (
const StatusCodeException &)
41 std::cout <<
"CheatingCosmicRayRemovalAlgorithm::Run - Unable to determine MCParticle origin for an input CaloHit, which will be skipped."
47 PANDORA_THROW_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::SaveList(*
this, outputCaloHitList, m_outputCaloHitListName));
48 PANDORA_THROW_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::ReplaceCurrentList<CaloHit>(*
this, m_outputCaloHitListName));
50 return STATUS_CODE_SUCCESS;
55 StatusCode CheatingCosmicRayRemovalAlgorithm::ReadSettings(
const TiXmlHandle xmlHandle)
57 PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, XmlHelper::ReadValue(xmlHandle,
"InputCaloHitListName", m_inputCaloHitListName));
58 PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, XmlHelper::ReadValue(xmlHandle,
"MCParticleListName", m_mcParticleListName));
59 PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, XmlHelper::ReadValue(xmlHandle,
"OutputCaloHitListName", m_outputCaloHitListName));
61 return STATUS_CODE_SUCCESS;
Header file for the cheating cosmic ray removal algorithm class.
Header file for the lar monte carlo particle helper helper class.
required by fuzzyCluster table::sbnd_g4_services gaushitTruthMatch pandora
BEGIN_PROLOG could also be cout