Inheritance diagram for simfilter::FilterNoDirtNeutrinos:

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

Public Member Functions inherited from art::EDFilter
	EDFilter (fhicl::ParameterSet const &pset)

template<typename Config >
	EDFilter (Table< Config > const &config)

std::string	workerType () const

Public Member Functions inherited from art::detail::Filter
virtual	~Filter () noexcept

	Filter (fhicl::ParameterSet const &)

	Filter (Filter const &)=delete

	Filter (Filter &&)=delete

Filter &	operator= (Filter const &)=delete

Filter &	operator= (Filter &&)=delete

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)

Public Member Functions inherited from art::Modifier
	~Modifier () noexcept

	Modifier ()

	Modifier (Modifier const &)=delete

	Modifier (Modifier &&)=delete

Modifier &	operator= (Modifier const &)=delete

Modifier &	operator= (Modifier &&)=delete

Public Member Functions inherited from art::ModuleBase
virtual	~ModuleBase () noexcept

	ModuleBase ()

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)

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

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

Private Member Functions
bool	filter (art::Event &) override

Private Attributes
std::string	fLArG4ModuleLabel

std::string	fGenModuleLabel

bool	fKeepCryostatNeutrinos

Additional Inherited Members
Public Types inherited from art::EDFilter
using	ModuleType = EDFilter

using	WorkerType = WorkerT< EDFilter >

Public Types inherited from art::detail::Filter
template<typename UserConfig >
using	Table = Modifier::Table< UserConfig >

Public Types inherited from art::Modifier
template<typename UserConfig , typename UserKeysToIgnore = void>
using	Table = ProducerTable< UserConfig, detail::ModuleConfig, UserKeysToIgnore >

Static Public Attributes inherited from art::detail::Filter
static constexpr bool	Pass {true}

static constexpr bool	Fail {false}

Protected Member Functions inherited from art::ModuleBase
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 ()

Detailed Description

Definition at line 33 of file FilterNoDirtNeutrinos_module.cc.

Constructor & Destructor Documentation

simfilter::FilterNoDirtNeutrinos::FilterNoDirtNeutrinos ( fhicl::ParameterSet const & pset )

explicit

Definition at line 54 of file FilterNoDirtNeutrinos_module.cc.

                                                                             : EDFilter{pset},
     fLArG4ModuleLabel    (pset.get< std::string > ("LArG4ModuleLabel"   , "NoLabel")       )
     , fGenModuleLabel    (pset.get< std::string > ("GenModuleLabel"  , "NoLabel")       )
     , fKeepCryostatNeutrinos    (pset.get< bool > ("KeepCryostatNeutrinos", false)      )
   {
   }

Member Function Documentation

bool simfilter::FilterNoDirtNeutrinos::filter ( art::Event & evt )

overrideprivatevirtual

Implements art::EDFilter.

Definition at line 62 of file FilterNoDirtNeutrinos_module.cc.

   {
     bool interactionDesired(false);
     //get the list of particles from this event
     art::ServiceHandle<geo::Geometry const> geom;
 
     // * MC truth information
     art::Handle< std::vector<simb::MCTruth> > mctruthListHandle;
     std::vector<art::Ptr<simb::MCTruth> > mclist;
     art::Handle<std::vector<simb::MCParticle> > mcpHandle;
 
 
 
     if (evt.getByLabel(fGenModuleLabel,mctruthListHandle))
       art::fill_ptr_vector(mclist, mctruthListHandle);
     evt.getByLabel(fLArG4ModuleLabel,mcpHandle);
 
     //    std::cout << "FilterNoDirtNeutrinos: mclist.size() is " << mclist.size()<< std::endl ;
 
   std::set<art::Ptr<simb::MCTruth> > mctSetGENIE;
   for(size_t i=0; i<mctruthListHandle->size(); ++i)
     {
       art::Ptr<simb::MCTruth> mct_ptr(mctruthListHandle,i);
       if( mctSetGENIE.find(mct_ptr) == mctSetGENIE.end() ) mctSetGENIE.insert(mct_ptr);
     }
 
   // Get the MCTruths from associations to our particles
   art::FindOneP<simb::MCTruth> assMCT(mcpHandle, evt, "largeant");
 
   double xmin;
   double xmax;
   double ymin;
   double ymax;
   double zmin;
   double zmax;
 
   if(fKeepCryostatNeutrinos)
   {
     // Get cryostat (box) volume boundary.
     xmin = geom->DetHalfWidth() - geom->CryostatHalfWidth();
     xmax = geom->DetHalfWidth() + geom->CryostatHalfWidth();
     ymin = -geom->CryostatHalfHeight();
     ymax = geom->CryostatHalfHeight();
     zmin = geom->DetLength()/2. - geom->DetLength()/2.;
     zmax = geom->DetLength()/2. + geom->DetLength()/2.;
   }
   else
   {
     // Get fiducial volume boundary.
     xmin = 0.;
     xmax = 2.*geom->DetHalfWidth();
     ymin = -geom->DetHalfHeight();
     ymax = geom->DetHalfHeight();
     zmin = 0.;
     zmax = geom->DetLength();
   }
 
   //  std::cout << "FilterNoDirtNeutrinos: mcpHandle->size() is " << mcpHandle->size()<< std::endl ;
   // Now let's loop over G4 MCParticle list and track back MCTruth
   bool inTPC (false);
   for(size_t i=0; i < mcpHandle->size() && !inTPC; ++i)
     {
       const art::Ptr<simb::MCParticle> mcp_ptr(mcpHandle,i);
       const art::Ptr<simb::MCTruth> &mct = assMCT.at(i);
       if( mctSetGENIE.find(mct) == mctSetGENIE.end() )
         {
           // This is non-genie
           continue;
         }
       else
         {
         // This is genie
 
           const simb::MCParticle* part(&mcpHandle->at(i));
           // for c2: pdg and trackID no longer used
           //int pdg = part->PdgCode();
           //int trackID = part->TrackId();
 
           //    std::cout << "FilterNoDirtNeutrinos: i is " << i << std::endl ;
           // Now walk through trajectory and see if it enters the TPC
           int n = part->NumberTrajectoryPoints();
           for(int j = 0; j < n && !inTPC; ++j)
             {
               //            std::cout << "FilterNoDirtNeutrinos: Loop  counter on NumTrajPt j is " << j << std::endl ;
 
               TVector3 pos = part->Position(j).Vect();
               if(pos.X() >= xmin &&
                  pos.X() <= xmax &&
                  pos.Y() >= ymin &&
                  pos.Y() <= ymax &&
                  pos.Z() >= zmin &&
                  pos.Z() <= zmax)
                 {
                   interactionDesired = true;
                   //              std::cout << "FilterNoDirtNeutrinos: Genie daughter found in TPC. G4Particle " << i << " , TrackID/pdg " << trackID << "/ " << pdg << " is discovered." << std::endl ;
                   std::cout << "FilterNoDirtNeutrinos: Genie daughter found in TPC. G4Particle " << std::endl ;
                   inTPC=true;
                 }
             } // trajectory loop
         } // end Genie particle
     } // loop on MCPHandle
 
   return interactionDesired;
 
   } // end FilterNoDirtNeutrinos()function