Inheritance diagram for filt::NucleonDecayFilter:

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

virtual	~NucleonDecayFilter ()

virtual bool	filter (art::Event &e)

void	reconfigure (fhicl::ParameterSet const &pset)

void	beginJob ()

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
void	FillVars (std::vector< int > &TrackIDVec, int &numParts, int ThisID)

Private Attributes
double	ActiveBounds [6]

std::map< int, const simb::MCParticle * >	truthmap

art::ServiceHandle< geo::Geometry >	geom

bool	fWantMuon

bool	fWantPion

bool	fWantPi0

bool	fWantKaon

bool	fWantElec

TTree *	MyOutTree

int	Run

int	SubRun

int	Event

int	nMuon

int	nPion

int	nPi0

int	nKaon

int	nElec

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 50 of file NucleonDecayFilter_module.cc.

Constructor & Destructor Documentation

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

explicit

virtual filt::NucleonDecayFilter::~NucleonDecayFilter ( )

inlinevirtual

Definition at line 53 of file NucleonDecayFilter_module.cc.

53 {};

Member Function Documentation

void filt::NucleonDecayFilter::beginJob ( )

virtual

Reimplemented from art::EDFilter.

Definition at line 195 of file NucleonDecayFilter_module.cc.

                                     {
     ActiveBounds[0] = ActiveBounds[2] = ActiveBounds[4] = DBL_MAX;
     ActiveBounds[1] = ActiveBounds[3] = ActiveBounds[5] = -DBL_MAX;
     // ----- FixMe: Assume single cryostats ------
     auto const* geom = lar::providerFrom<geo::Geometry>();
     for (geo::TPCGeo const& TPC: geom->IterateTPCs()) {
       // get center in world coordinates
       const double origin[3] = {0.};
       double center[3] = {0.};
       TPC.LocalToWorld(origin, center);
       double tpcDim[3] = {TPC.HalfWidth(), TPC.HalfHeight(), 0.5*TPC.Length() };
       if( center[0] - tpcDim[0] < ActiveBounds[0] ) ActiveBounds[0] = center[0] - tpcDim[0];
       if( center[0] + tpcDim[0] > ActiveBounds[1] ) ActiveBounds[1] = center[0] + tpcDim[0];
       if( center[1] - tpcDim[1] < ActiveBounds[2] ) ActiveBounds[2] = center[1] - tpcDim[1];
       if( center[1] + tpcDim[1] > ActiveBounds[3] ) ActiveBounds[3] = center[1] + tpcDim[1];
       if( center[2] - tpcDim[2] < ActiveBounds[4] ) ActiveBounds[4] = center[2] - tpcDim[2];
       if( center[2] + tpcDim[2] > ActiveBounds[5] ) ActiveBounds[5] = center[2] + tpcDim[2];
     } // for all TPC
     std::cout << "Active Boundaries: "
               << "\n\tx: " << ActiveBounds[0] << " to " << ActiveBounds[1]
               << "\n\ty: " << ActiveBounds[2] << " to " << ActiveBounds[3]
               << "\n\tz: " << ActiveBounds[4] << " to " << ActiveBounds[5]
               << std::endl;
   
     art::ServiceHandle<art::TFileService> tfs;
     MyOutTree = tfs->make<TTree>("FilterTree","FilterTree");
     MyOutTree->Branch( "Run",   &Run   , "Run/I"    );
     MyOutTree->Branch( "SubRun",&SubRun, "SubRun/I" );
     MyOutTree->Branch( "Event", &Event , "Event/I"  );
     MyOutTree->Branch( "nMuon", &nMuon , "nMuon/I"  );
     MyOutTree->Branch( "nPion", &nPion , "nPion/I"  );
     MyOutTree->Branch( "nPi0" , &nPi0  , "nPi0/I"   );
     MyOutTree->Branch( "nKaon", &nKaon , "nKaon/I"  );
     MyOutTree->Branch( "nElec", &nElec , "nElec/I"  );
   }

void filt::NucleonDecayFilter::FillVars	(	std::vector< int > &	TrackIDVec,
		int &	numParts,
		int	ThisID
	)

private

Definition at line 231 of file NucleonDecayFilter_module.cc.

                                                                                            {
     std::vector<int>::iterator it=std::find( TrackIDVec.begin(), TrackIDVec.end(), abs(ThisID) );
     if ( it==TrackIDVec.end() ) {
       TrackIDVec.push_back ( abs(ThisID) ); // Push back this particle type IDVec
       ++numParts;
     }
   }

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

virtual

Implements art::EDFilter.

Definition at line 104 of file NucleonDecayFilter_module.cc.

                                               {
     Run    = e.run();
     SubRun = e.subRun();
     Event  = e.event();
     nMuon  = nPion = nPi0 = nKaon = nElec = 0;
 
     // Any providers I need.
     auto const* geo = lar::providerFrom<geo::Geometry>();
     
     // Make a map of MCParticles which I can access later.
     auto truth = e.getHandle<std::vector<simb::MCParticle> >("largeant");
     truthmap.clear();
     for (size_t i=0; i<truth->size(); i++)
       truthmap[truth->at(i).TrackId()]=&((*truth)[i]);
     //std::cout << "The primary muon in this event had an energy of " << truthmap[1]->E() << std::endl;
 
     // Get a vector of sim channels.
     auto simchannels = e.getHandle<std::vector<sim::SimChannel> >("largeant");
 
     // Make vectors to hold all of my particle TrackIDs
     std::vector<int> MuonVec, PionVec, Pi0Vec, KaonVec, ElecVec;
     
     // ------ Now loop through all of my hits ------
     for (auto const& simchannel:*simchannels) {
     // ------ Only want to look at collection plane hits ------
       if (geo->SignalType(simchannel.Channel()) != geo::kCollection) continue;
       // ------ Loop through all the IDEs for this channel ------
       for (auto const& tdcide:simchannel.TDCIDEMap()) {
         auto const& idevec=tdcide.second;
         // ------ Look at each individual IDE ------
         for (auto const& ide:idevec) {
           int ideTrackID = ide.trackID;
           double idePos[3];
           idePos[0] = ide.x;
           idePos[1] = ide.y;
           idePos[2] = ide.z;
           // ------ Is the hit in a TPC? ------
           geo::TPCID tpcid=geo->FindTPCAtPosition(idePos);
           if (!(geo->HasTPC(tpcid)) ) {
             //std::cout << "Outside the Active volume I found at the top!" << std::endl;
             continue;
           }
           // ------ Now to work out which particles in particular I want to save more information about... ------
           // ----------------- I want to write out the IDE to the relevant parent particle type -----------------
           // ------- This means looping back through the IDE's parents until I find an interesting TrackID ------
           bool WrittenOut = false;
           while ( ideTrackID != 0 && !WrittenOut ) {
             const simb::MCParticle& part=*( truthmap[ abs(ideTrackID) ] );
             int PdgCode=part.PdgCode();
             WrittenOut = true; // If one of my particles don't want to go further back.
             // ========== Muons ==========
             if      ( (PdgCode == -13  || PdgCode == 13)  ) 
               FillVars( MuonVec, nMuon, ideTrackID );
             // ========== Pions ==========
             else if ( (PdgCode == -211 || PdgCode == 211) && part.Process() != "pi+Inelastic" && part.Process() != "pi-Inelastic" ) 
               FillVars( PionVec, nPion, ideTrackID );
             // ========== Pi0s  ==========
             else if ( PdgCode == 111  )
               FillVars( Pi0Vec , nPi0 , ideTrackID );
             // ========== Kaons ===========
             else if ( (PdgCode == 321 || PdgCode == -321) && part.Process() != "kaon+Inelastic" && part.Process() != "kaon-Inelastic" ) 
               FillVars( KaonVec, nKaon, ideTrackID );
             // ========== Elecs ===========
             else if ( (PdgCode == -11 || PdgCode == 11)  )
               FillVars( ElecVec, nElec, ideTrackID );
             // ========== If still not one of my intersting particles I need to find this particles parent ==========
             else {
               ideTrackID = part.Mother();
               PdgCode = truthmap[ abs(ideTrackID) ]->PdgCode();
               WrittenOut = false; // If not one of my particles change back to false.
             } // If not one of chosen particles.
           } // While loop
         } // Each IDE ( ide:idevec )
       } // IDE vector for SimChannel ( tdcide:simcahnnel.TPCIDEMap() )
     } // Loop through simchannels
     /*
     std::cout << "Looking at Run " << Run << " SubRun " << SubRun << " Event " << Event
               << ". There were " << nMuon << " muons, " << nPion << " pions, " << nPi0 << " pi0s, " << nKaon << " Kaons, " << nElec << " Electrons."
               << std::endl;
     */
     // Fill Tree and decide whether to keep this event...
     MyOutTree->Fill();
     bool EventPasses = true;
     if ( (fWantMuon && nMuon == 0) || nMuon > MaxPart) EventPasses = false;
     if ( (fWantPion && nPion == 0) || nPi0  > MaxPart) EventPasses = false;
     if ( (fWantPi0  && nPi0  == 0) || nPion > MaxPart) EventPasses = false;
     if ( (fWantKaon && nKaon == 0) || nKaon > MaxPart) EventPasses = false;
     if ( (fWantElec && nElec == 0) || nElec > MaxPart) EventPasses = false;
     return EventPasses;
   }

void filt::NucleonDecayFilter::reconfigure ( fhicl::ParameterSet const & pset )

Definition at line 89 of file NucleonDecayFilter_module.cc.

                                                                     {
     fWantMuon = pset.get<bool>("WantMuon");
     fWantPion = pset.get<bool>("WantPion");
     fWantPi0  = pset.get<bool>("WantPi0" );
     fWantKaon = pset.get<bool>("WantKaon");
     fWantElec = pset.get<bool>("WantElec");
     std::cout << "Your filter is selecting only events with: "
               << "Muons " << fWantMuon << ", "
               << "Pions " << fWantPion << ", "
               << "Pi0s "  << fWantPi0  << ", "
               << "Kaons " << fWantKaon << ", "
               << "Elecs " << fWantElec << "."
               << std::endl;
   }