Inheritance diagram for simfilter::FilterNoDirtNeutrinos:

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

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

Private Attributes
std::string	fLArG4ModuleLabel

std::string	fGenModuleLabel

bool	fKeepCryostatNeutrinos

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 )

overrideprivate

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