Inheritance diagram for reco::shower::LArPandoraModularShowerCreation:

Public Types
using	ModuleType = EDProducer

template<typename UserConfig , typename KeysToIgnore = void>
using	Table = Modifier::Table< UserConfig, KeysToIgnore >

Public Member Functions
	LArPandoraModularShowerCreation (fhicl::ParameterSet const &pset)

void	doBeginJob (SharedResources const &resources)

void	doEndJob ()

void	doRespondToOpenInputFile (FileBlock const &fb)

void	doRespondToCloseInputFile (FileBlock const &fb)

void	doRespondToOpenOutputFiles (FileBlock const &fb)

void	doRespondToCloseOutputFiles (FileBlock const &fb)

bool	doBeginRun (RunPrincipal &rp, ModuleContext const &mc)

bool	doEndRun (RunPrincipal &rp, ModuleContext const &mc)

bool	doBeginSubRun (SubRunPrincipal &srp, ModuleContext const &mc)

bool	doEndSubRun (SubRunPrincipal &srp, ModuleContext const &mc)

bool	doEvent (EventPrincipal &ep, ModuleContext const &mc, std::atomic< std::size_t > &counts_run, std::atomic< std::size_t > &counts_passed, std::atomic< std::size_t > &counts_failed)

void	fillProductDescriptions ()

void	registerProducts (ProductDescriptions &productsToRegister)

ModuleDescription const &	moduleDescription () const

void	setModuleDescription (ModuleDescription const &)

std::array< std::vector< ProductInfo >, NumBranchTypes > const &	getConsumables () const

void	sortConsumables (std::string const &current_process_name)

std::unique_ptr< Worker >	makeWorker (WorkerParams const &wp)

template<typename T , BranchType BT>
ViewToken< T >	consumesView (InputTag const &tag)

template<typename T , BranchType BT>
ViewToken< T >	mayConsumeView (InputTag const &tag)

Protected Member Functions
ConsumesCollector &	consumesCollector ()

template<typename T , BranchType = InEvent>
ProductToken< T >	consumes (InputTag const &)

template<typename Element , BranchType = InEvent>
ViewToken< Element >	consumesView (InputTag const &)

template<typename T , BranchType = InEvent>
void	consumesMany ()

template<typename T , BranchType = InEvent>
ProductToken< T >	mayConsume (InputTag const &)

template<typename Element , BranchType = InEvent>
ViewToken< Element >	mayConsumeView (InputTag const &)

template<typename T , BranchType = InEvent>
void	mayConsumeMany ()

Private Member Functions
void	produce (art::Event &evt)

template<class T >
art::Ptr< T >	GetProducedElementPtr (const std::string &InstanceName, const reco::shower::ShowerElementHolder &ShowerEleHolder, const int &iter=-1)

Private Attributes
unsigned int	fNumPlanes

const art::InputTag	fPFParticleLabel

const bool	fAllowPartialShowers

const int	fVerbose

const bool	fUseAllParticles

const std::string	fShowerStartPositionLabel

const std::string	fShowerDirectionLabel

const std::string	fShowerEnergyLabel

const std::string	fShowerLengthLabel

const std::string	fShowerOpeningAngleLabel

const std::string	fShowerdEdxLabel

const std::string	fShowerBestPlaneLabel

std::vector< std::unique_ptr< ShowerRecoTools::IShowerTool > >	fShowerTools

std::vector< std::string >	fShowerToolNames

reco::shower::ShowerProducedPtrsHolder	uniqueproducerPtrs

art::ServiceHandle< geo::Geometry >	fGeom

Detailed Description

Definition at line 37 of file LArPandoraModularShowerCreation_module.cc.

Member Typedef Documentation

using art::EDProducer::ModuleType = EDProducer

inherited

Definition at line 17 of file EDProducer.h.

template<typename UserConfig , typename KeysToIgnore = void>

using art::detail::Producer::Table = Modifier::Table<UserConfig, KeysToIgnore>

inherited

Definition at line 26 of file Producer.h.

Constructor & Destructor Documentation

reco::shower::LArPandoraModularShowerCreation::LArPandoraModularShowerCreation ( fhicl::ParameterSet const & pset )

Definition at line 113 of file LArPandoraModularShowerCreation_module.cc.

References fAllowPartialShowers, fGeom, fNumPlanes, fPFParticleLabel, fShowerBestPlaneLabel, fShowerdEdxLabel, fShowerDirectionLabel, fShowerEnergyLabel, fShowerLengthLabel, fShowerOpeningAngleLabel, fShowerStartPositionLabel, fShowerToolNames, fShowerTools, fUseAllParticles, fVerbose, fhicl::ParameterSet::get(), art::InputTag::instance(), art::InputTag::label(), geo::GeometryCore::Nplanes(), reco::shower::ShowerProducedPtrsHolder::PrintPtrs(), art::InputTag::process(), art::ProductRegistryHelper::producesCollector(), fhicl::ParameterSet::put(), fhicl::ParameterSet::put_or_replace(), reco::shower::ShowerProducedPtrsHolder::SetShowerUniqueProduerPtr(), and uniqueproducerPtrs.

   : EDProducer{pset}
   , fPFParticleLabel(pset.get<art::InputTag>("PFParticleLabel"))
   , fAllowPartialShowers(pset.get<bool>("AllowPartialShowers"))
   , fVerbose(pset.get<int>("Verbose", 0))
   , fUseAllParticles(pset.get<bool>("UseAllParticles", false))
   , fShowerStartPositionLabel(pset.get<std::string>("ShowerStartPositionLabel"))
   , fShowerDirectionLabel(pset.get<std::string>("ShowerDirectionLabel"))
   , fShowerEnergyLabel(pset.get<std::string>("ShowerEnergyLabel"))
   , fShowerLengthLabel(pset.get<std::string>("ShowerLengthLabel"))
   , fShowerOpeningAngleLabel(pset.get<std::string>("ShowerOpeningAngleLabel"))
   , fShowerdEdxLabel(pset.get<std::string>("ShowerdEdxLabel"))
   , fShowerBestPlaneLabel(pset.get<std::string>("ShowerBestPlaneLabel"))
 {
   //Intialise the tools
   auto tool_psets = pset.get<std::vector<fhicl::ParameterSet>>("ShowerFinderTools");
   for (auto& tool_pset : tool_psets) {
 
     const std::string tool_name(tool_pset.get<std::string>("tool_type"));
     // If the PFPaticle label is not set for the tool, make it use the one for the module
     // Note we also need to set the label in the Alg, via the base tool
     if (!tool_pset.has_key("PFParticleLabel")) {
 
       // I cannot pass an art::InputTag as it is mangled, so lets make a string instead
       const std::string PFParticleLabelString(fPFParticleLabel.label() + ":" +
                                               fPFParticleLabel.instance() + ":" +
                                               fPFParticleLabel.process());
 
       tool_pset.put<std::string>("PFParticleLabel", PFParticleLabelString);
       fhicl::ParameterSet base_pset = tool_pset.get<fhicl::ParameterSet>("BaseTools");
       fhicl::ParameterSet alg_pset = base_pset.get<fhicl::ParameterSet>("LArPandoraShowerAlg");
       alg_pset.put<std::string>("PFParticleLabel", PFParticleLabelString);
       base_pset.put_or_replace<fhicl::ParameterSet>("LArPandoraShowerAlg", alg_pset);
       tool_pset.put_or_replace<fhicl::ParameterSet>("BaseTools", base_pset);
 
       if (tool_pset.has_key("LArPandoraShowerCheatingAlg")) {
         fhicl::ParameterSet cheat_alg_pset =
           tool_pset.get<fhicl::ParameterSet>("LArPandoraShowerCheatingAlg");
         cheat_alg_pset.put<std::string>("PFParticleLabel", PFParticleLabelString);
         cheat_alg_pset.put_or_replace<fhicl::ParameterSet>("LArPandoraShowerAlg", alg_pset);
         tool_pset.put_or_replace<fhicl::ParameterSet>("LArPandoraShowerCheatingAlg",
                                                       cheat_alg_pset);
       }
     }
 
     // If we have not explicitly set verboseness for a given tool, use global level
     if (!tool_pset.has_key("Verbose")) { tool_pset.put<int>("Verbose", fVerbose); }
 
     fShowerTools.push_back(art::make_tool<ShowerRecoTools::IShowerTool>(tool_pset));
     fShowerToolNames.push_back(tool_name);
 
     fNumPlanes = fGeom->Nplanes();
   }
 
   //  Initialise the EDProducer ptr in the tools
   std::vector<std::string> SetupTools;
   for (unsigned int i = 0; i < fShowerTools.size(); ++i) {
     if (std::find(SetupTools.begin(), SetupTools.end(), fShowerToolNames[i]) != SetupTools.end()) {
       continue;
     }
     fShowerTools[i]->SetPtr(&producesCollector());
     fShowerTools[i]->InitaliseProducerPtr(uniqueproducerPtrs);
     fShowerTools[i]->InitialiseProducers();
   }
 
   //Initialise the other paramters.
 
   produces<std::vector<recob::Shower>>();
   produces<art::Assns<recob::Shower, recob::Hit>>();
   produces<art::Assns<recob::Shower, recob::Cluster>>();
   produces<art::Assns<recob::Shower, recob::SpacePoint>>();
   produces<art::Assns<recob::Shower, recob::PFParticle>>();
 
   // Output -- showers and associations with hits and clusters
   uniqueproducerPtrs.SetShowerUniqueProduerPtr(type<std::vector<recob::Shower>>(), "shower");
   uniqueproducerPtrs.SetShowerUniqueProduerPtr(type<art::Assns<recob::Shower, recob::Cluster>>(),
                                                "clusterAssociationsbase");
   uniqueproducerPtrs.SetShowerUniqueProduerPtr(type<art::Assns<recob::Shower, recob::Hit>>(),
                                                "hitAssociationsbase");
   uniqueproducerPtrs.SetShowerUniqueProduerPtr(type<art::Assns<recob::Shower, recob::SpacePoint>>(),
                                                "spShowerAssociationsbase");
   uniqueproducerPtrs.SetShowerUniqueProduerPtr(type<art::Assns<recob::Shower, recob::PFParticle>>(),
                                                "pfShowerAssociationsbase");
 
   uniqueproducerPtrs.PrintPtrs();
 }

Member Function Documentation

template<typename T , BranchType BT>

ProductToken< T > art::ModuleBase::consumes ( InputTag const & tag )

protectedinherited

Definition at line 61 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::consumes().

   {
     return collector_.consumes<T, BT>(tag);
   }

ConsumesCollector & art::ModuleBase::consumesCollector ( )

protectedinherited

Definition at line 57 of file ModuleBase.cc.

References art::ModuleBase::collector_.

   {
     return collector_;
   }

template<typename T , BranchType BT>

void art::ModuleBase::consumesMany ( )

protectedinherited

Definition at line 75 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::consumesMany().

   {
     collector_.consumesMany<T, BT>();
   }

template<typename Element , BranchType = InEvent>

ViewToken<Element> art::ModuleBase::consumesView ( InputTag const & )

protectedinherited

template<typename T , BranchType BT>

ViewToken<T> art::ModuleBase::consumesView ( InputTag const & tag )

inherited

Definition at line 68 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::consumesView().

   {
     return collector_.consumesView<T, BT>(tag);
   }

void art::detail::Producer::doBeginJob ( SharedResources const & resources )

inherited

Definition at line 22 of file Producer.cc.

References art::detail::Producer::beginJobWithFrame(), and art::detail::Producer::setupQueues().

   {
     setupQueues(resources);
     ProcessingFrame const frame{ScheduleID{}};
     beginJobWithFrame(frame);
   }

bool art::detail::Producer::doBeginRun	(	RunPrincipal &	rp,
		ModuleContext const &	mc
	)

inherited

Definition at line 65 of file Producer.cc.

References art::detail::Producer::beginRunWithFrame(), art::RangeSet::forRun(), art::RunPrincipal::makeRun(), r, art::RunPrincipal::runID(), and art::ModuleContext::scheduleID().

   {
     auto r = rp.makeRun(mc, RangeSet::forRun(rp.runID()));
     ProcessingFrame const frame{mc.scheduleID()};
     beginRunWithFrame(r, frame);
     r.commitProducts();
     return true;
   }

bool art::detail::Producer::doBeginSubRun	(	SubRunPrincipal &	srp,
		ModuleContext const &	mc
	)

inherited

Definition at line 85 of file Producer.cc.

References art::detail::Producer::beginSubRunWithFrame(), art::RangeSet::forSubRun(), art::SubRunPrincipal::makeSubRun(), art::ModuleContext::scheduleID(), and art::SubRunPrincipal::subRunID().

   {
     auto sr = srp.makeSubRun(mc, RangeSet::forSubRun(srp.subRunID()));
     ProcessingFrame const frame{mc.scheduleID()};
     beginSubRunWithFrame(sr, frame);
     sr.commitProducts();
     return true;
   }

void art::detail::Producer::doEndJob ( )

inherited

Definition at line 30 of file Producer.cc.

References art::detail::Producer::endJobWithFrame().

   {
     ProcessingFrame const frame{ScheduleID{}};
     endJobWithFrame(frame);
   }

bool art::detail::Producer::doEndRun	(	RunPrincipal &	rp,
		ModuleContext const &	mc
	)

inherited

Definition at line 75 of file Producer.cc.

References art::detail::Producer::endRunWithFrame(), art::RunPrincipal::makeRun(), r, art::ModuleContext::scheduleID(), and art::Principal::seenRanges().

   {
     auto r = rp.makeRun(mc, rp.seenRanges());
     ProcessingFrame const frame{mc.scheduleID()};
     endRunWithFrame(r, frame);
     r.commitProducts();
     return true;
   }

bool art::detail::Producer::doEndSubRun	(	SubRunPrincipal &	srp,
		ModuleContext const &	mc
	)

inherited

Definition at line 95 of file Producer.cc.

References art::detail::Producer::endSubRunWithFrame(), art::SubRunPrincipal::makeSubRun(), art::ModuleContext::scheduleID(), and art::Principal::seenRanges().

   {
     auto sr = srp.makeSubRun(mc, srp.seenRanges());
     ProcessingFrame const frame{mc.scheduleID()};
     endSubRunWithFrame(sr, frame);
     sr.commitProducts();
     return true;
   }

bool art::detail::Producer::doEvent	(	EventPrincipal &	ep,
		ModuleContext const &	mc,
		std::atomic< std::size_t > &	counts_run,
		std::atomic< std::size_t > &	counts_passed,
		std::atomic< std::size_t > &	counts_failed
	)

inherited

Definition at line 105 of file Producer.cc.

References art::detail::Producer::checkPutProducts_, e, art::EventPrincipal::makeEvent(), art::detail::Producer::produceWithFrame(), and art::ModuleContext::scheduleID().

   {
     auto e = ep.makeEvent(mc);
     ++counts_run;
     ProcessingFrame const frame{mc.scheduleID()};
     produceWithFrame(e, frame);
     e.commitProducts(checkPutProducts_, &expectedProducts<InEvent>());
     ++counts_passed;
     return true;
   }

void art::detail::Producer::doRespondToCloseInputFile ( FileBlock const & fb )

inherited

Definition at line 44 of file Producer.cc.

References art::detail::Producer::respondToCloseInputFileWithFrame().

   {
     ProcessingFrame const frame{ScheduleID{}};
     respondToCloseInputFileWithFrame(fb, frame);
   }

void art::detail::Producer::doRespondToCloseOutputFiles ( FileBlock const & fb )

inherited

Definition at line 58 of file Producer.cc.

References art::detail::Producer::respondToCloseOutputFilesWithFrame().

   {
     ProcessingFrame const frame{ScheduleID{}};
     respondToCloseOutputFilesWithFrame(fb, frame);
   }

void art::detail::Producer::doRespondToOpenInputFile ( FileBlock const & fb )

inherited

Definition at line 37 of file Producer.cc.

References art::detail::Producer::respondToOpenInputFileWithFrame().

   {
     ProcessingFrame const frame{ScheduleID{}};
     respondToOpenInputFileWithFrame(fb, frame);
   }

void art::detail::Producer::doRespondToOpenOutputFiles ( FileBlock const & fb )

inherited

Definition at line 51 of file Producer.cc.

References art::detail::Producer::respondToOpenOutputFilesWithFrame().

   {
     ProcessingFrame const frame{ScheduleID{}};
     respondToOpenOutputFilesWithFrame(fb, frame);
   }

void art::Modifier::fillProductDescriptions ( )

inherited

Definition at line 10 of file Modifier.cc.

References art::ProductRegistryHelper::fillDescriptions(), and art::ModuleBase::moduleDescription().

   {
     ProductRegistryHelper::fillDescriptions(moduleDescription());
   }

std::array< std::vector< ProductInfo >, NumBranchTypes > const & art::ModuleBase::getConsumables ( ) const

inherited

Definition at line 43 of file ModuleBase.cc.

References art::ModuleBase::collector_, and art::ConsumesCollector::getConsumables().

   {
     return collector_.getConsumables();
   }

template<class T >

art::Ptr< T > reco::shower::LArPandoraModularShowerCreation::GetProducedElementPtr	(	const std::string &	InstanceName,
		const reco::shower::ShowerElementHolder &	ShowerEleHolder,
		const int &	iter = `-1`
	)

private

Definition at line 83 of file LArPandoraModularShowerCreation_module.cc.

References reco::shower::ShowerElementHolder::CheckElement(), reco::shower::ShowerProducedPtrsHolder::CheckUniqueProduerPtr(), reco::shower::ShowerProducedPtrsHolder::GetArtPtr(), reco::shower::ShowerElementHolder::GetShowerNumber(), and uniqueproducerPtrs.

 {
 
   bool check_element = ShowerEleHolder.CheckElement(InstanceName);
   if (!check_element) {
     throw cet::exception("LArPandoraModularShowerCreation")
       << "To get a element that does not exist" << std::endl;
   }
 
   bool check_ptr = uniqueproducerPtrs.CheckUniqueProduerPtr(InstanceName);
   if (!check_ptr) {
     throw cet::exception("LArPandoraModularShowerCreation")
       << "Tried to get a ptr that does not exist" << std::endl;
   }
 
   //Get the number of the shower we are on.
   int index;
   if (iter != -1) { index = iter; }
   else {
     index = ShowerEleHolder.GetShowerNumber();
   }
 
   //Make the ptr
   art::Ptr<T> artptr = uniqueproducerPtrs.GetArtPtr<T>(InstanceName, index);
   return artptr;
 }

std::unique_ptr< Worker > art::ModuleBase::makeWorker ( WorkerParams const & wp )

inherited

Definition at line 37 of file ModuleBase.cc.

References art::ModuleBase::doMakeWorker(), and art::NumBranchTypes.

   {
     return doMakeWorker(wp);
   }

template<typename T , BranchType BT>

ProductToken< T > art::ModuleBase::mayConsume ( InputTag const & tag )

protectedinherited

Definition at line 82 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::mayConsume().

   {
     return collector_.mayConsume<T, BT>(tag);
   }

template<typename T , BranchType BT>

void art::ModuleBase::mayConsumeMany ( )

protectedinherited

Definition at line 96 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::mayConsumeMany().

   {
     collector_.mayConsumeMany<T, BT>();
   }

template<typename Element , BranchType = InEvent>

ViewToken<Element> art::ModuleBase::mayConsumeView ( InputTag const & )

protectedinherited

template<typename T , BranchType BT>

ViewToken<T> art::ModuleBase::mayConsumeView ( InputTag const & tag )

inherited

Definition at line 89 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::mayConsumeView().

   {
     return collector_.mayConsumeView<T, BT>(tag);
   }

ModuleDescription const & art::ModuleBase::moduleDescription ( ) const

inherited

Definition at line 13 of file ModuleBase.cc.

References art::errors::LogicError.

Referenced by art::OutputModule::doRespondToOpenInputFile(), art::OutputModule::doWriteEvent(), art::Modifier::fillProductDescriptions(), art::OutputModule::makePlugins_(), art::OutputWorker::OutputWorker(), produce(), art::Modifier::registerProducts(), and art::OutputModule::registerProducts().

   {
     if (md_.has_value()) {
       return *md_;
     }
 
     throw Exception{errors::LogicError,
                     "There was an error while calling moduleDescription().\n"}
       << "The moduleDescription() base-class member function cannot be called\n"
          "during module construction.  To determine which module is "
          "responsible\n"
          "for calling it, find the '<module type>:<module "
          "label>@Construction'\n"
          "tag in the message prefix above.  Please contact artists@fnal.gov\n"
          "for guidance.\n";
   }

void reco::shower::LArPandoraModularShowerCreation::produce ( art::Event & evt )

privatevirtual

Check all the producers

Implements art::EDProducer.

Definition at line 201 of file LArPandoraModularShowerCreation_module.cc.

 {
 
   //Ptr makers for the products
   uniqueproducerPtrs.SetPtrMakers(evt);
   reco::shower::ShowerElementHolder showerEleHolder;
 
   //Get the PFParticles
   auto const pfpHandle = evt.getValidHandle<std::vector<recob::PFParticle>>(fPFParticleLabel);
   std::vector<art::Ptr<recob::PFParticle>> pfps;
   art::fill_ptr_vector(pfps, pfpHandle);
 
   //Handle to access the pandora hits assans
   auto const clusterHandle = evt.getValidHandle<std::vector<recob::Cluster>>(fPFParticleLabel);
 
   //Get the assoications to hits, clusters and spacespoints
   const art::FindManyP<recob::Hit>& fmh =
     showerEleHolder.GetFindManyP<recob::Hit>(clusterHandle, evt, fPFParticleLabel);
   const art::FindManyP<recob::Cluster>& fmcp =
     showerEleHolder.GetFindManyP<recob::Cluster>(pfpHandle, evt, fPFParticleLabel);
   const art::FindManyP<recob::SpacePoint>& fmspp =
     showerEleHolder.GetFindManyP<recob::SpacePoint>(pfpHandle, evt, fPFParticleLabel);
 
   //Holder to pass to the functions, contains the 6 properties of the shower
   // - Start Poistion
   // - Direction
   // - Initial Track
   // - Initial Track Hits
   // - Energy
   // - dEdx
   // - Length
   // - Opening Angle
 
   int shower_iter = 0;
   //Loop of the pf particles
   for (auto const& pfp : pfps) {
 
     //Update the shower iterator
     showerEleHolder.SetShowerNumber(shower_iter);
 
     //loop only over showers unless otherwise specified
     if (!fUseAllParticles && pfp->PdgCode() != 11 && pfp->PdgCode() != 22) continue;
 
     //Get the associated hits,clusters and spacepoints
     const std::vector<art::Ptr<recob::Cluster>> showerClusters = fmcp.at(pfp.key());
     const std::vector<art::Ptr<recob::SpacePoint>> showerSpacePoints = fmspp.at(pfp.key());
 
     // Check the pfp has at least 1 cluster (i.e. not a pfp neutrino)
     if (!showerClusters.size()) continue;
 
     if (fVerbose > 1)
       mf::LogInfo("LArPandoraModularShowerCreation")
         << "Running on shower: " << shower_iter << std::endl;
 
     //Calculate the shower properties
     //Loop over the shower tools
     int err = 0;
     for (unsigned int i = 0; i < fShowerTools.size(); i++) {
 
       //Calculate the metric
       if (fVerbose > 1)
         mf::LogInfo("LArPandoraModularShowerCreation")
           << "Running shower tool: " << fShowerToolNames[i] << std::endl;
       std::string evd_disp_append = fShowerToolNames[i] + "_iteration" + std::to_string(0) + "_" +
                                     this->moduleDescription().moduleLabel();
 
       err = fShowerTools[i]->RunShowerTool(pfp, evt, showerEleHolder, evd_disp_append);
 
       if (err && fVerbose) {
         mf::LogError("LArPandoraModularShowerCreation")
           << "Error in shower tool: " << fShowerToolNames[i] << " with code: " << err << std::endl;
       }
     }
 
     //If we are are not allowing partial shower check all of the things
     if (!fAllowPartialShowers) {
       // If we recieved an error call from a tool return;
 
       // Check everything we need is in the shower element holder
       if (!showerEleHolder.CheckElement(fShowerStartPositionLabel)) {
         if (fVerbose)
           mf::LogError("LArPandoraModularShowerCreation")
             << "The start position is not set in the element holder. bailing" << std::endl;
         continue;
       }
       if (!showerEleHolder.CheckElement(fShowerDirectionLabel)) {
         if (fVerbose)
           mf::LogError("LArPandoraModularShowerCreation")
             << "The direction is not set in the element holder. bailing" << std::endl;
         continue;
       }
       if (!showerEleHolder.CheckElement(fShowerEnergyLabel)) {
         if (fVerbose)
           mf::LogError("LArPandoraModularShowerCreation")
             << "The energy is not set in the element holder. bailing" << std::endl;
         continue;
       }
       if (!showerEleHolder.CheckElement(fShowerdEdxLabel)) {
         if (fVerbose)
           mf::LogError("LArPandoraModularShowerCreation")
             << "The dEdx is not set in the element holder. bailing" << std::endl;
         continue;
       }
       if (!showerEleHolder.CheckElement(fShowerBestPlaneLabel)) {
         if (fVerbose)
           mf::LogError("LArPandoraModularShowerCreation")
             << "The BestPlane is not set in the element holder. bailing" << std::endl;
         continue;
       }
       if (!showerEleHolder.CheckElement(fShowerLengthLabel)) {
         if (fVerbose)
           mf::LogError("LArPandoraModularShowerCreation")
             << "The length is not set in the element holder. bailing" << std::endl;
         continue;
       }
       if (!showerEleHolder.CheckElement(fShowerOpeningAngleLabel)) {
         if (fVerbose)
           mf::LogError("LArPandoraModularShowerCreation")
             << "The opening angle is not set in the element holder. bailing" << std::endl;
         continue;
       }
 
       //Check All of the products that have been asked to be checked.
       bool elements_are_set = showerEleHolder.CheckAllElementTags();
       if (!elements_are_set) {
         if (fVerbose)
           mf::LogError("LArPandoraModularShowerCreation")
             << "Not all the elements in the property holder which should be set are not. Bailing. "
             << std::endl;
         continue;
       }
 
       bool producers_are_set = uniqueproducerPtrs.CheckAllProducedElements(showerEleHolder);
       if (!producers_are_set) {
         if (fVerbose)
           mf::LogError("LArPandoraModularShowerCreation")
             << "Not all the elements in the property holder which are produced are not set. "
                "Bailing. "
             << std::endl;
         continue;
       }
     }
 
     //Get the properties
     geo::Point_t ShowerStartPosition(-999, -999, -999);
     geo::Vector_t ShowerDirection(-999, -999, -999);
     std::vector<double> ShowerEnergy(fNumPlanes, -999);
     std::vector<double> ShowerdEdx(fNumPlanes, -999);
     int BestPlane(-999);
     double ShowerLength(-999);
     double ShowerOpeningAngle(-999);
 
     geo::Point_t ShowerStartPositionErr(-999, -999, -999);
     geo::Vector_t ShowerDirectionErr(-999, -999, -999);
     std::vector<double> ShowerEnergyErr(fNumPlanes, -999);
     std::vector<double> ShowerdEdxErr(fNumPlanes, -999);
 
     err = 0;
     if (showerEleHolder.CheckElement(fShowerStartPositionLabel))
       err += showerEleHolder.GetElementAndError(
         fShowerStartPositionLabel, ShowerStartPosition, ShowerStartPositionErr);
     if (showerEleHolder.CheckElement(fShowerDirectionLabel))
       err += showerEleHolder.GetElementAndError(
         fShowerDirectionLabel, ShowerDirection, ShowerDirectionErr);
     if (showerEleHolder.CheckElement(fShowerEnergyLabel))
       err += showerEleHolder.GetElementAndError(fShowerEnergyLabel, ShowerEnergy, ShowerEnergyErr);
     if (showerEleHolder.CheckElement(fShowerdEdxLabel))
       err += showerEleHolder.GetElementAndError(fShowerdEdxLabel, ShowerdEdx, ShowerdEdxErr);
     if (showerEleHolder.CheckElement(fShowerBestPlaneLabel))
       err += showerEleHolder.GetElement(fShowerBestPlaneLabel, BestPlane);
     if (showerEleHolder.CheckElement(fShowerLengthLabel))
       err += showerEleHolder.GetElement(fShowerLengthLabel, ShowerLength);
     if (showerEleHolder.CheckElement(fShowerOpeningAngleLabel))
       err += showerEleHolder.GetElement(fShowerOpeningAngleLabel, ShowerOpeningAngle);
 
     if (err) {
       throw cet::exception("LArPandoraModularShowerCreation")
         << "Error in LArPandoraModularShowerCreation Module. A Check on a shower property failed "
         << std::endl;
     }
 
     if (fVerbose > 1) {
       //Check the shower
       std::cout << "Shower Vertex: X:" << ShowerStartPosition.X()
                 << " Y: " << ShowerStartPosition.Y() << " Z: " << ShowerStartPosition.Z()
                 << std::endl;
       std::cout << "Shower Direction: X:" << ShowerDirection.X() << " Y: " << ShowerDirection.Y()
                 << " Z: " << ShowerDirection.Z() << std::endl;
       std::cout << "Shower dEdx:";
       for (unsigned int i = 0; i < fNumPlanes; i++) {
         std::cout << " Plane " << i << ": " << ShowerdEdx.at(i);
       }
       std::cout << std::endl;
       std::cout << "Shower Energy:";
       for (unsigned int i = 0; i < fNumPlanes; i++) {
         std::cout << " Plane " << i << ": " << ShowerEnergy.at(i);
       }
       std::cout << std::endl;
       std::cout << "Shower Best Plane: " << BestPlane << std::endl;
       std::cout << "Shower Length: " << ShowerLength << std::endl;
       std::cout << "Shower Opening Angle: " << ShowerOpeningAngle << std::endl;
 
       //Print what has been created in the shower
       showerEleHolder.PrintElements();
     }
 
     if (ShowerdEdx.size() != fNumPlanes) {
       throw cet::exception("LArPandoraModularShowerCreation")
         << "dEdx vector is wrong size: " << ShowerdEdx.size()
         << " compared to Nplanes: " << fNumPlanes << std::endl;
     }
     if (ShowerEnergy.size() != fNumPlanes) {
       throw cet::exception("LArPandoraModularShowerCreation")
         << "Energy vector is wrong size: " << ShowerEnergy.size()
         << " compared to Nplanes: " << fNumPlanes << std::endl;
     }
 
     //Make the shower
     using namespace geo::vect;
     recob::Shower shower(convertTo<TVector3>(ShowerDirection),
                          convertTo<TVector3>(ShowerDirectionErr),
                          convertTo<TVector3>(ShowerStartPosition),
                          convertTo<TVector3>(ShowerDirectionErr),
                          ShowerEnergy,
                          ShowerEnergyErr,
                          ShowerdEdx,
                          ShowerdEdxErr,
                          BestPlane,
                          util::kBogusI,
                          ShowerLength,
                          ShowerOpeningAngle);
     showerEleHolder.SetElement(shower, "shower");
     ++shower_iter;
     art::Ptr<recob::Shower> ShowerPtr =
       this->GetProducedElementPtr<recob::Shower>("shower", showerEleHolder);
 
     //Associate the pfparticle
     uniqueproducerPtrs.AddSingle<art::Assns<recob::Shower, recob::PFParticle>>(
       ShowerPtr, pfp, "pfShowerAssociationsbase");
 
     //Add the hits for each "cluster"
     for (auto const& cluster : showerClusters) {
 
       //Associate the clusters
       std::vector<art::Ptr<recob::Hit>> ClusterHits = fmh.at(cluster.key());
       uniqueproducerPtrs.AddSingle<art::Assns<recob::Shower, recob::Cluster>>(
         ShowerPtr, cluster, "clusterAssociationsbase");
 
       //Associate the hits
       for (auto const& hit : ClusterHits) {
         uniqueproducerPtrs.AddSingle<art::Assns<recob::Shower, recob::Hit>>(
           ShowerPtr, hit, "hitAssociationsbase");
       }
     }
 
     //Associate the spacepoints
     for (auto const& sp : showerSpacePoints) {
       uniqueproducerPtrs.AddSingle<art::Assns<recob::Shower, recob::SpacePoint>>(
         ShowerPtr, sp, "spShowerAssociationsbase");
     }
 
     //Loop over the tool data products and add them.
     uniqueproducerPtrs.AddDataProducts(showerEleHolder);
 
     //AddAssociations
     int assn_err = 0;
     for (auto const& fShowerTool : fShowerTools) {
       //AddAssociations
       assn_err += fShowerTool->AddAssociations(pfp, evt, showerEleHolder);
     }
     if (!fAllowPartialShowers && assn_err > 0) {
       if (fVerbose)
         mf::LogError("LArPandoraModularShowerCreation")
           << "A association failed and not allowing partial showers. The association will not be "
              "added to the event "
           << std::endl;
     }
 
     //Reset the showerproperty holder.
     showerEleHolder.ClearShower();
   }
 
   //Put everything in the event.
   uniqueproducerPtrs.MoveAllToEvent(evt);
 
   //Reset the ptrs to the data products
   uniqueproducerPtrs.reset();
 }

void art::Modifier::registerProducts ( ProductDescriptions & productsToRegister )

inherited

Definition at line 16 of file Modifier.cc.

References art::ModuleBase::moduleDescription(), and art::ProductRegistryHelper::registerProducts().

   {
     ProductRegistryHelper::registerProducts(productsToRegister,
                                             moduleDescription());
   }

void art::ModuleBase::setModuleDescription ( ModuleDescription const & md )

inherited

Definition at line 31 of file ModuleBase.cc.

References art::ModuleBase::md_.

   {
     md_ = md;
   }

void art::ModuleBase::sortConsumables ( std::string const & current_process_name )

inherited

Definition at line 49 of file ModuleBase.cc.

References art::ModuleBase::collector_, and art::ConsumesCollector::sortConsumables().

   {
     // Now that we know we have seen all the consumes declarations,
     // sort the results for fast lookup later.
     collector_.sortConsumables(current_process_name);
   }

Member Data Documentation

const bool reco::shower::LArPandoraModularShowerCreation::fAllowPartialShowers

private

Definition at line 55 of file LArPandoraModularShowerCreation_module.cc.

Referenced by LArPandoraModularShowerCreation(), and produce().

art::ServiceHandle<geo::Geometry> reco::shower::LArPandoraModularShowerCreation::fGeom

private

Definition at line 76 of file LArPandoraModularShowerCreation_module.cc.

Referenced by LArPandoraModularShowerCreation().

unsigned int reco::shower::LArPandoraModularShowerCreation::fNumPlanes

private

Definition at line 53 of file LArPandoraModularShowerCreation_module.cc.

Referenced by LArPandoraModularShowerCreation(), and produce().

const art::InputTag reco::shower::LArPandoraModularShowerCreation::fPFParticleLabel

private

Definition at line 54 of file LArPandoraModularShowerCreation_module.cc.

Referenced by LArPandoraModularShowerCreation(), and produce().

const std::string reco::shower::LArPandoraModularShowerCreation::fShowerBestPlaneLabel

private

Definition at line 66 of file LArPandoraModularShowerCreation_module.cc.

Referenced by LArPandoraModularShowerCreation(), and produce().

const std::string reco::shower::LArPandoraModularShowerCreation::fShowerdEdxLabel

private

Definition at line 65 of file LArPandoraModularShowerCreation_module.cc.

Referenced by LArPandoraModularShowerCreation(), and produce().

const std::string reco::shower::LArPandoraModularShowerCreation::fShowerDirectionLabel

private

Definition at line 61 of file LArPandoraModularShowerCreation_module.cc.

Referenced by LArPandoraModularShowerCreation(), and produce().

const std::string reco::shower::LArPandoraModularShowerCreation::fShowerEnergyLabel

private

Definition at line 62 of file LArPandoraModularShowerCreation_module.cc.

Referenced by LArPandoraModularShowerCreation(), and produce().

const std::string reco::shower::LArPandoraModularShowerCreation::fShowerLengthLabel

private

Definition at line 63 of file LArPandoraModularShowerCreation_module.cc.

Referenced by LArPandoraModularShowerCreation(), and produce().

const std::string reco::shower::LArPandoraModularShowerCreation::fShowerOpeningAngleLabel

private

Definition at line 64 of file LArPandoraModularShowerCreation_module.cc.

Referenced by LArPandoraModularShowerCreation(), and produce().

const std::string reco::shower::LArPandoraModularShowerCreation::fShowerStartPositionLabel

private

Definition at line 60 of file LArPandoraModularShowerCreation_module.cc.

Referenced by LArPandoraModularShowerCreation(), and produce().

std::vector<std::string> reco::shower::LArPandoraModularShowerCreation::fShowerToolNames

private

Definition at line 70 of file LArPandoraModularShowerCreation_module.cc.

Referenced by LArPandoraModularShowerCreation(), and produce().

std::vector<std::unique_ptr<ShowerRecoTools::IShowerTool> > reco::shower::LArPandoraModularShowerCreation::fShowerTools

private

Definition at line 69 of file LArPandoraModularShowerCreation_module.cc.

Referenced by LArPandoraModularShowerCreation(), and produce().

const bool reco::shower::LArPandoraModularShowerCreation::fUseAllParticles

private

Definition at line 57 of file LArPandoraModularShowerCreation_module.cc.

Referenced by LArPandoraModularShowerCreation(), and produce().

const int reco::shower::LArPandoraModularShowerCreation::fVerbose

private

Definition at line 56 of file LArPandoraModularShowerCreation_module.cc.

Referenced by LArPandoraModularShowerCreation(), and produce().

reco::shower::ShowerProducedPtrsHolder reco::shower::LArPandoraModularShowerCreation::uniqueproducerPtrs

private

Definition at line 73 of file LArPandoraModularShowerCreation_module.cc.

Referenced by GetProducedElementPtr(), LArPandoraModularShowerCreation(), and produce().

The documentation for this class was generated from the following file:

larpandora/v09_22_05/source/larpandora/LArPandoraEventBuilding/LArPandoraShower/LArPandoraModularShowerCreation_module.cc

Public Types

Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation