Inheritance diagram for filt::LArG4ParticleFilter:

Public Types
using	ModuleType = EDFilter

template<typename UserConfig >
using	Table = Modifier::Table< UserConfig >

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

virtual bool	filter (art::Event &e)

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 &)

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)

Static Public Attributes
static constexpr bool	Pass {true}

static constexpr bool	Fail {false}

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
bool	IsInterestingParticle (const art::Ptr< simb::MCParticle > &particle, const unsigned int index) const

bool	PDGCheck (const art::Ptr< simb::MCParticle > &particle, const unsigned int index) const

bool	IsPrimaryCheck (const art::Ptr< simb::MCParticle > &particle, const unsigned int index) const

bool	MinMomentumCheck (const art::Ptr< simb::MCParticle > &particle, const unsigned int index) const

bool	MaxMomentumCheck (const art::Ptr< simb::MCParticle > &particle, const unsigned int index) const

bool	StartInTPCCheck (const art::Ptr< simb::MCParticle > &particle, const unsigned int index) const

bool	StopInTPCCheck (const art::Ptr< simb::MCParticle > &particle, const unsigned int index) const

bool	TPCTrajLengthCheck (const art::Ptr< simb::MCParticle > &particle, const unsigned int index) const

double	CalculateLength (const std::vector< TVector3 > &position_segment) const

void	VerifyDataMembers () const

void	VerifyElementRequest (const unsigned int index) const

Private Attributes
art::ServiceHandle< geo::Geometry const >	fGeom

const std::string	fLArG4ModuleLabel

const std::vector< int >	fInterestingPDGs

const std::vector< int >	fIsPrimary

const std::vector< double >	fParticleMinMomentum

const std::vector< double >	fParticleMaxMomentum

const std::vector< int >	fStartInTPC

const std::vector< int >	fStopInTPC

const std::vector< double >	fParticleMinTPCLength

const bool	fRequireAllInterestingParticles

std::vector< bool >	fFoundInterestingParticles

Detailed Description

Definition at line 24 of file LArG4ParticleFilter_module.cc.

Member Typedef Documentation

using art::EDFilter::ModuleType = EDFilter

inherited

Definition at line 18 of file EDFilter.h.

template<typename UserConfig >

using art::detail::Filter::Table = Modifier::Table<UserConfig>

inherited

Definition at line 32 of file Filter.h.

Constructor & Destructor Documentation

filt::LArG4ParticleFilter::LArG4ParticleFilter ( fhicl::ParameterSet const & pset )

explicit

Member Function Documentation

double filt::LArG4ParticleFilter::CalculateLength ( const std::vector< TVector3 > & position_segment ) const

private

Definition at line 270 of file LArG4ParticleFilter_module.cc.

Referenced by TPCTrajLengthCheck().

   {
     double length = 0;
     //Check how many points we have in the segment.  If it is one or less, there is nothing to calculate so return 0
     if (position_segment.size() <= 1) return length;
 
     //Now we need to compare every adjacent pair of positions to work out the length of this segment
     for (unsigned int i = 1; i < position_segment.size(); i++) {
       length += (position_segment[i] - position_segment[i - 1]).Mag();
     }
 
     return length;
   }

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::Filter::doBeginJob ( SharedResources const & resources )

inherited

Definition at line 29 of file Filter.cc.

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

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

inherited

Definition at line 72 of file Filter.cc.

References 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()};
     bool const rc = beginRunWithFrame(r, frame);
     r.commitProducts();
     return rc;
   }

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

inherited

Definition at line 92 of file Filter.cc.

References art::SubRunPrincipal::makeSubRun(), art::ModuleContext::scheduleID(), and art::SubRunPrincipal::subRunID().

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

void art::detail::Filter::doEndJob ( )

inherited

Definition at line 37 of file Filter.cc.

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

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

inherited

Definition at line 82 of file Filter.cc.

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

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

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

inherited

Definition at line 102 of file Filter.cc.

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

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

bool art::detail::Filter::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 112 of file Filter.cc.

References e, art::EventPrincipal::makeEvent(), and art::ModuleContext::scheduleID().

   {
     auto e = ep.makeEvent(mc);
     ++counts_run;
     ProcessingFrame const frame{mc.scheduleID()};
     bool const rc = filterWithFrame(e, frame);
     e.commitProducts(checkPutProducts_, &expectedProducts<InEvent>());
     if (rc) {
       ++counts_passed;
     } else {
       ++counts_failed;
     }
     return rc;
   }

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

inherited

Definition at line 51 of file Filter.cc.

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

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

inherited

Definition at line 65 of file Filter.cc.

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

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

inherited

Definition at line 44 of file Filter.cc.

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

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

inherited

Definition at line 58 of file Filter.cc.

   {
     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());
   }

bool filt::LArG4ParticleFilter::filter ( art::Event & e )

virtual

Implements art::EDFilter.

Definition at line 80 of file LArG4ParticleFilter_module.cc.

References fFoundInterestingParticles, fInterestingPDGs, fLArG4ModuleLabel, fRequireAllInterestingParticles, art::ProductRetriever::getByLabel(), and IsInterestingParticle().

   {
     //Reset the found vector
     for (unsigned int i = 0; i < fFoundInterestingParticles.size(); i++) {
       fFoundInterestingParticles[i] = false;
     }
 
     //Get the vector of particles
     art::Handle<std::vector<simb::MCParticle>> particles;
     e.getByLabel(fLArG4ModuleLabel, particles);
     //Loop through the particles
     for (unsigned int part_i = 0; part_i < particles->size(); part_i++) {
       //Get the particle
       const art::Ptr<simb::MCParticle> particle(particles, part_i);
       //Loop over the list of particles we want and compare it with the particle we are looking it
       for (unsigned int interest_i = 0; interest_i < fInterestingPDGs.size(); interest_i++) {
         if (IsInterestingParticle(particle, interest_i)) {
           if (!fRequireAllInterestingParticles)
             return true; //If we only require at least one of the particles be found then we have already done our job
           fFoundInterestingParticles[interest_i] = true;
           bool foundThemAll = true;
           for (unsigned int found_i = 0; found_i < fFoundInterestingParticles.size(); found_i++) {
             if (fFoundInterestingParticles[found_i] == false) {
               foundThemAll = false;
               break;
             }
           }
           if (foundThemAll) { return true; }
         }
       }
     }
     //Assume that the event is not worth saving
     return false;
   }

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();
   }

bool filt::LArG4ParticleFilter::IsInterestingParticle	(	const art::Ptr< simb::MCParticle > &	particle,
		const unsigned int	index
	)		const

private

Definition at line 115 of file LArG4ParticleFilter_module.cc.

References IsPrimaryCheck(), MaxMomentumCheck(), MinMomentumCheck(), PDGCheck(), StartInTPCCheck(), StopInTPCCheck(), TPCTrajLengthCheck(), and VerifyElementRequest().

Referenced by filter().

   {
     VerifyElementRequest(index);
     //Run the checks
     if (!PDGCheck(particle, index)) return false;
     if (!IsPrimaryCheck(particle, index)) return false;
     if (!MinMomentumCheck(particle, index)) return false;
     if (!MaxMomentumCheck(particle, index)) return false;
     if (!StartInTPCCheck(particle, index)) return false;
     if (!StopInTPCCheck(particle, index)) return false;
     if (!TPCTrajLengthCheck(particle, index)) return false;
 
     return true;
   }

bool filt::LArG4ParticleFilter::IsPrimaryCheck	(	const art::Ptr< simb::MCParticle > &	particle,
		const unsigned int	index
	)		const

private

Definition at line 140 of file LArG4ParticleFilter_module.cc.

References fIsPrimary, and simb::MCParticle::Mother().

Referenced by IsInterestingParticle().

   {
     const int isPrimaryCondition(fIsPrimary[index]);
     if (isPrimaryCondition == -1) return true;
     bool particlePrimaryStatus(particle->Mother() > 0 ? false : true);
     if (particlePrimaryStatus != isPrimaryCondition) return false;
     return true;
   }

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);
   }

bool filt::LArG4ParticleFilter::MaxMomentumCheck	(	const art::Ptr< simb::MCParticle > &	particle,
		const unsigned int	index
	)		const

private

Definition at line 159 of file LArG4ParticleFilter_module.cc.

References fParticleMaxMomentum, and simb::MCParticle::Momentum().

Referenced by IsInterestingParticle().

   {
     if (fParticleMaxMomentum[index] > 0 &&
         particle->Momentum(0).Vect().Mag() > fParticleMaxMomentum[index])
       return false;
     return true;
   }

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);
   }

bool filt::LArG4ParticleFilter::MinMomentumCheck	(	const art::Ptr< simb::MCParticle > &	particle,
		const unsigned int	index
	)		const

private

Definition at line 150 of file LArG4ParticleFilter_module.cc.

References fParticleMinMomentum, and simb::MCParticle::Momentum().

Referenced by IsInterestingParticle().

   {
     if (fParticleMinMomentum[index] > 0 &&
         particle->Momentum(0).Vect().Mag() < fParticleMinMomentum[index])
       return false;
     return true;
   }

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(), reco::shower::LArPandoraModularShowerCreation::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";
   }

bool filt::LArG4ParticleFilter::PDGCheck	(	const art::Ptr< simb::MCParticle > &	particle,
		const unsigned int	index
	)		const

private

Definition at line 131 of file LArG4ParticleFilter_module.cc.

References fInterestingPDGs, and simb::MCParticle::PdgCode().

Referenced by IsInterestingParticle().

   {
     int pdg = fInterestingPDGs[index];
     if (pdg == 0) return true; //User does not care what the PDG is
     if (particle->PdgCode() != pdg) return false;
     return true;
   }

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);
   }

bool filt::LArG4ParticleFilter::StartInTPCCheck	(	const art::Ptr< simb::MCParticle > &	particle,
		const unsigned int	index
	)		const

private

Definition at line 168 of file LArG4ParticleFilter_module.cc.

References fGeom, geo::GeometryCore::FindTPCAtPosition(), fStartInTPC, geo::CryostatID::isValid, simb::MCParticle::Position(), and geo::vect::toPoint().

Referenced by IsInterestingParticle().

   {
     //Firstly check if we even care if the particle starts in the TPC or not
     int demand = fStartInTPC[index];
     if (demand == 0)
       return true; //We don't care if the particle starts in the TPC or not so pass the check
     //Get TPC corresponding to starting position of particle
     geo::TPCID tpcid = fGeom->FindTPCAtPosition(geo::vect::toPoint(particle->Position(0).Vect()));
     //Now we need to compare if we have a TPC that we started in with whether we wanted to start in a TPC at all
     if (tpcid.isValid) {
       //The particle DID start in a TPC.  Now, did we WANT this to happen
       return demand == 1;
     }
     else {
       //The particle did NOT start in a TPC.  Did we WANT this to happen?
       return demand == 2;
     }
 
     //Assume true by default
     return true;
   }

bool filt::LArG4ParticleFilter::StopInTPCCheck	(	const art::Ptr< simb::MCParticle > &	particle,
		const unsigned int	index
	)		const

private

Definition at line 191 of file LArG4ParticleFilter_module.cc.

References fGeom, geo::GeometryCore::FindTPCAtPosition(), fStopInTPC, geo::CryostatID::isValid, simb::MCParticle::NumberTrajectoryPoints(), simb::MCParticle::Position(), and geo::vect::toPoint().

Referenced by IsInterestingParticle().

   {
     //Firstly check if we even care if the particle stops in the TPC or not
     int demand = fStopInTPC[index];
     if (demand == 0)
       return true; //We don't care if the particle stops in the TPC or not so pass the check
     //Get final position of particle
     auto const final_position = particle->Position(particle->NumberTrajectoryPoints() - 1).Vect();
 
     geo::TPCID tpcid = fGeom->FindTPCAtPosition(geo::vect::toPoint(final_position));
     //Now we need to compare if we have a TPC that we stopped in with whether we wanted to stop in a TPC at all
     if (tpcid.isValid) {
       //The particle DID stop in a TPC.  Now, did we WANT this to happen
       return demand == 1;
     }
     else {
       //The particle did NOT stop in a TPC.  Did we WANT this to happen?
       return demand == 2;
     }
 
     //Assume true by default
     return true;
   }

bool filt::LArG4ParticleFilter::TPCTrajLengthCheck	(	const art::Ptr< simb::MCParticle > &	particle,
		const unsigned int	index
	)		const

private

Definition at line 216 of file LArG4ParticleFilter_module.cc.

References CalculateLength(), fGeom, geo::GeometryCore::FindTPCAtPosition(), fParticleMinTPCLength, geo::CryostatID::isValid, simb::MCParticle::NumberTrajectoryPoints(), simb::MCParticle::Position(), and geo::vect::toPoint().

Referenced by IsInterestingParticle().

   {
     double min_traj_length = fParticleMinTPCLength[index];
 
     //Firstly, if we don't care about the TPC trajectory length then pass the check
     if (min_traj_length < 0) return true;
 
     //Now the hard bit
     //To do this, we need to collect the sequential particle positions which are contained in the TPC into segments.  The reason for doing this is that the particle may enter a TPC, leave it and enter it again and if this isn't taken into account, the length might be grossly overestimated
     //It is easiest to store the positions in a vector of vectors
     bool OK = false;
     std::vector<std::vector<TVector3>> position_segments;
     //We are also going to need an empty vector to store in the above vector
     std::vector<TVector3> position_segment;
     //Loop through the trajectory points
     for (unsigned int i = 0; i < particle->NumberTrajectoryPoints(); i++) {
       //Extract the current position of the particle
       geo::TPCID curr_tpcid =
         fGeom->FindTPCAtPosition(geo::vect::toPoint(particle->Position(i).Vect()));
       //There are a couple of things to check here.  If the particle is currently in the TPC, then we need to store that particular position.  If it is NOT in the TPC, then its either exited the TPC or has not yet entered.  If it has just exited, then the position_segment should have some positions stored in it, it which case we now need to store this segment.  If it has not yet entered the TPC, then we don't need to do anything
       //If it is currently in the TPC
       if (curr_tpcid.isValid) position_segment.push_back(particle->Position(i).Vect());
       //It has just exited the TPC
       else if (position_segment.size() > 0) {
         //Store the segment
         position_segments.push_back(position_segment);
         //Now reset the segment
         position_segment.clear();
       }
       //There is nothing to do because the particle has remained outside of the TPC
     }
     //We need to check once more if the position_segment vector has been filled
     if (position_segment.size() > 0) {
       position_segments.push_back(position_segment);
       position_segment.clear();
     }
     //Now lets check the length of each segment
     //Firstly, if we didn't store a segment then the particle fails the check
     if (position_segments.size() == 0) return false;
     //Now loop through the segments and check if they are above threshold
     for (unsigned int i = 0; i < position_segments.size(); i++) {
       double segment_length = CalculateLength(position_segments[i]);
       if (segment_length > min_traj_length) {
         //We found a track segment in the TPC which passes the length threshold so don't flag as bad
         OK = true;
         break;
       }
     }
     if (!OK) return false;
 
     return true;
   }

void filt::LArG4ParticleFilter::VerifyDataMembers ( ) const

private

Definition at line 284 of file LArG4ParticleFilter_module.cc.

References art::errors::Configuration, fInterestingPDGs, fIsPrimary, fParticleMaxMomentum, fParticleMinMomentum, fParticleMinTPCLength, fStartInTPC, and fStopInTPC.

   {
     if (fInterestingPDGs.size() != fIsPrimary.size())
       throw art::Exception(art::errors::Configuration)
         << "Config error: InterestingPDGs and IsPrimary FCL vectors are different sizes\n";
     if (fInterestingPDGs.size() != fParticleMinMomentum.size())
       throw art::Exception(art::errors::Configuration)
         << "Config error: InterestingPDGs and ParticleMinMomentum FCL vectors are different "
            "sizes\n";
     if (fInterestingPDGs.size() != fParticleMaxMomentum.size())
       throw art::Exception(art::errors::Configuration)
         << "Config error: InterestingPDGs and ParticleMaxMomentum FCL vectors are different "
            "sizes\n";
     if (fInterestingPDGs.size() != fStartInTPC.size())
       throw art::Exception(art::errors::Configuration)
         << "Config error: InterestingPDGs and StartInTPC FCL vectors are different sizes\n";
     if (fInterestingPDGs.size() != fStopInTPC.size())
       throw art::Exception(art::errors::Configuration)
         << "Config error: InterestingPDGs and StopInTPC FCL vectors are different sizes\n";
     if (fInterestingPDGs.size() != fParticleMinTPCLength.size())
       throw art::Exception(art::errors::Configuration)
         << "Config error: InterestingPDGs and ParticleMinTPCLength FCL vectors are different "
            "sizes\n";
 
     for (unsigned int iIsPrimary = 0; iIsPrimary < fIsPrimary.size(); iIsPrimary++)
       if (!(-1 == fIsPrimary[iIsPrimary] || 0 == fIsPrimary[iIsPrimary] ||
             1 == fIsPrimary[iIsPrimary]))
         throw art::Exception(art::errors::Configuration)
           << "Config error: Element " << iIsPrimary << " of the IsPrimary vector equals "
           << fIsPrimary[iIsPrimary]
           << ".  Valid options are -1 (not set), 0 (not primary) or 1 (is primary)" << std::endl;
 
     for (unsigned int iStartInTPC = 0; iStartInTPC < fStartInTPC.size(); iStartInTPC++)
       if (!(0 == fStartInTPC[iStartInTPC] || 1 == fStartInTPC[iStartInTPC] ||
             2 == fStartInTPC[iStartInTPC]))
         throw art::Exception(art::errors::Configuration)
           << "Config error: Element " << iStartInTPC << " of the StartInTPC vector equals "
           << fStartInTPC[iStartInTPC]
           << ".  Valid options are 0 (not set), 1 (start in TPC) or 2 (do not start in TPC)"
           << std::endl;
 
     for (unsigned int iStopInTPC = 0; iStopInTPC < fStopInTPC.size(); iStopInTPC++)
       if (!(0 == fStopInTPC[iStopInTPC] || 1 == fStopInTPC[iStopInTPC] ||
             2 == fStopInTPC[iStopInTPC]))
         throw art::Exception(art::errors::Configuration)
           << "Config error: Element " << iStopInTPC << " of the StopInTPC vector equals "
           << fStopInTPC[iStopInTPC]
           << ".  Valid options are 0 (not set), 1 (stop in TPC) or 2 (do not stop in TPC)"
           << std::endl;
 
     return;
   }

void filt::LArG4ParticleFilter::VerifyElementRequest ( const unsigned int index ) const

private

Definition at line 337 of file LArG4ParticleFilter_module.cc.

References DEFINE_ART_MODULE, fInterestingPDGs, fIsPrimary, fParticleMaxMomentum, fParticleMinMomentum, fParticleMinTPCLength, fStartInTPC, fStopInTPC, and art::errors::InvalidNumber.

Referenced by IsInterestingParticle().

   {
     if (index > fInterestingPDGs.size())
       throw art::Exception(art::errors::InvalidNumber)
         << "Bounds error:  Requested element " << index
         << " from InterestingPDGs vector which is size " << fInterestingPDGs.size() << "\n";
     if (index > fIsPrimary.size())
       throw art::Exception(art::errors::InvalidNumber)
         << "Bounds error:  Requested element " << index << " from IsPrimary vector which is size "
         << fIsPrimary.size() << "\n";
     if (index > fParticleMinMomentum.size())
       throw art::Exception(art::errors::InvalidNumber)
         << "Bounds error:  Requested element " << index
         << " from ParticleMinMomentum vector which is size " << fParticleMinMomentum.size() << "\n";
     if (index > fParticleMaxMomentum.size())
       throw art::Exception(art::errors::InvalidNumber)
         << "Bounds error:  Requested element " << index
         << " from ParticleMaxMomentum vector which is size " << fParticleMaxMomentum.size() << "\n";
     if (index > fStartInTPC.size())
       throw art::Exception(art::errors::InvalidNumber)
         << "Bounds error:  Requested element " << index << " from StartInTPC vector which is size "
         << fStartInTPC.size() << "\n";
     if (index > fStopInTPC.size())
       throw art::Exception(art::errors::InvalidNumber)
         << "Bounds error:  Requested element " << index << " from StopInTPC vector which is size "
         << fStopInTPC.size() << "\n";
     if (index > fParticleMinTPCLength.size())
       throw art::Exception(art::errors::InvalidNumber)
         << "Bounds error:  Requested element " << index
         << " from ParticleMinTPCLength vector which is size " << fParticleMinTPCLength.size()
         << "\n";
     return;
   }