Inheritance diagram for filt::CFilter:

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

virtual	~CFilter ()

virtual bool	filter (art::Event &e)

void	reconfigure (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 Attributes
std::string	fSimulationProducerLabel

std::string	fRawDigitLabel

int	fCounterType

int	fSelectedPDG

int	fTrigger

int	fnumTracks

std::string	fHitsModuleLabel

std::string	fTrackModuleLabel

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 38 of file CFilter_module.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 64 of file CFilter_module.cc.

                                                  : EDFilter{pset}
   {
     this->reconfigure(pset);
 
   }

virtual filt::CFilter::~CFilter ( )

inlinevirtual

Definition at line 41 of file CFilter_module.cc.

41 { }

Member Function Documentation

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

virtual

Implements art::EDFilter.

Definition at line 84 of file CFilter_module.cc.

   {
     //muon conter geometry
     int counters_loaded=-1;
     std::vector<std::vector<double> > countergeometry;
     //load the muon counter positions from a text file
 
     char  counterfile[]= "/afs/fnal.gov/files/home/room1/melnimr/lr_dev_7/srcs/dunetpc/dune/Geometry/muoncounters.txt";
     counters_loaded=geo::MuonCounter35Alg::loadMuonCounterGeometry(counterfile,countergeometry);
     bool keepFlag=0;
 
     //    int icount=0; int trkind;
     // define the handle as an MCParticle vector and fill it with events from the simulation
     auto particleHandle = event.getHandle< std::vector<simb::MCParticle> >(fSimulationProducerLabel);
 
     // define a sorted map in which to put the particles
     std::map< int, const simb::MCParticle* > particleMap;
     // loop over all the particles, find the primary muon, and get the initial/final positions
     for ( auto const& particle : (*particleHandle) )
       {
         // For the methods you can call to get particle information,
         // see $NUTOOLS_INC/SimulationBase/MCParticle.h.
         int fTrackID = particle.TrackId();
         // Add the address of the MCParticle to the map, with the track ID as the key.
         particleMap[fTrackID] = &particle; 
         // PDG code of every particle in the event.
         //        int fPDG = particle.PdgCode();
         // locate the primary muon
         //        if ( particle.Process() == "primary"  &&  fPDG == fSelectedPDG )
         // {
 
             // A particle has a trajectory, consisting of a set of
             // 4-positions and 4-mommenta.
             size_t numberTrajectoryPoints = particle.NumberTrajectoryPoints();
             // For trajectories, as for vectors and arrays, the
             // first point is #0, not #1.
             int last = numberTrajectoryPoints - 1;
             const TLorentzVector& positionStart = particle.Position(0);
             const TLorentzVector& positionEnd   = particle.Position(last);
             const TLorentzVector& momentumStart = particle.Momentum(0);
             const TLorentzVector& momentumEnd   = particle.Momentum(last);
             // move the initial position (for muon counter studies)
             double x_increment = 0.; // 349.666
             double z_increment = 0.; // 231.065
             TVector3 trackStart(positionStart.X()+x_increment,positionStart.Y(),positionStart.Z()+z_increment);
             //}
 
             std::vector< std::vector<double> > hitcounters;         
             if(counters_loaded){
               //unsigned int counters_hit=0;
               
               geo::MuonCounter35Alg::testTrackInAllCounters(fTrackID,trackStart,momentumStart.Vect(),
                                                             countergeometry,hitcounters);
             }
             
 
             // condition flags for each layer
             bool Layer_1_2 = false;
             bool Layer_3_4_5 = false;
             bool Layer_E = false;
             bool Layer_W = false;
             bool Layer_N_U = false;
             bool Layer_N_L = false;
             bool Layer_S_U = false;
             bool Layer_S_L = false;
             int Trigger= 0;
 
             for(unsigned int ii=0;ii<hitcounters.size();ii++)
               {
                 //      for(unsigned int jj=0;jj<hitcounters[ii].size();jj++)
                 //      {
                 std::cout<< "::Filter, counters triggered are: "<<hitcounters[ii][0] <<std::endl;
                 /*        if (hitcounters[ii][0]==fCounterType) 
                   {
                   keepFlag=1;
                   return keepFlag;
                   }*/
                 // check which layer is hit
                 if( 40 <= hitcounters[ii][0] && hitcounters[ii][0] <= 61){
                   Layer_1_2 = true;
                 }
                 if( hitcounters[ii][0] > 61){
                   Layer_3_4_5 = true;
                 }
                 if (14 <= hitcounters[ii][0] && hitcounters[ii][0] <=19){
                   Layer_N_U = true;
                 }
                 if ( 34 <= hitcounters[ii][0] && hitcounters[ii][0] <= 39){
                   Layer_S_L = true;
                 }
                 if (8 <= hitcounters[ii][0] && hitcounters[ii][0] <= 13){
                   Layer_N_L = true;
                 }
                 if (28 <= hitcounters[ii][0] && hitcounters[ii][0] <= 33){
                   Layer_S_U = true;
                 }
                 if (hitcounters[ii][0] <= 7){
                   Layer_E = true;
                 }
                 if (20 <= hitcounters[ii][0] && hitcounters[ii][0] <= 27){
                   Layer_W = true;
                 }
                 //      }
               }
             // check for a satisfied trigger condition
             if (Layer_1_2 && Layer_3_4_5){
               Trigger = 1;
             }
             if (Layer_N_U && Layer_S_L){
               Trigger = 2;
             }
             if (Layer_N_L && Layer_S_U){
                 Trigger = 3;
             }
             if (Layer_E && Layer_W){
               Trigger = 4;
             }
             std::cout << "Trigger is ....." << Trigger << std::endl;
             if(Trigger==fTrigger)
               {
                 keepFlag=1;
                 return keepFlag;
               }
       }
     return keepFlag;
   }

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

Definition at line 70 of file CFilter_module.cc.

   {
     
     fSimulationProducerLabel = p.get< std::string >("SimulationLabel");
     //    fCounterFile = p.get< std::string >("CounterFile");
     fSelectedPDG = p.get< int >("PDGcode");
     fCounterType=p.get< int >("CounterType");
     fTrigger=p.get< int >("Trigger");
     fnumTracks=1;
 
     //      fnumTracks = pset.get<int>("numTracks");
       
   }

Member Data Documentation

int filt::CFilter::fCounterType

private

Definition at line 51 of file CFilter_module.cc.

std::string filt::CFilter::fHitsModuleLabel

private

Definition at line 59 of file CFilter_module.cc.

int filt::CFilter::fnumTracks

private

Definition at line 58 of file CFilter_module.cc.

std::string filt::CFilter::fRawDigitLabel

private

Definition at line 49 of file CFilter_module.cc.

int filt::CFilter::fSelectedPDG

private

Definition at line 52 of file CFilter_module.cc.

std::string filt::CFilter::fSimulationProducerLabel

private

Definition at line 48 of file CFilter_module.cc.

std::string filt::CFilter::fTrackModuleLabel

private

Definition at line 60 of file CFilter_module.cc.

int filt::CFilter::fTrigger

private

Definition at line 53 of file CFilter_module.cc.

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

dunetpc/dune/CFilter/CFilter_module.cc

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation