Inheritance diagram for MyGapFilter::GapFilter:

Public Member Functions
	GapFilter (fhicl::ParameterSet const &p)

virtual	~GapFilter ()

void	beginJob ()

void	reconfigure (fhicl::ParameterSet const &pset)

void	analyze (const art::Event &evt)

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

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

std::string	workerType () const

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

	Analyzer (fhicl::ParameterSet const &pset)

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

void	doBeginJob ()

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

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::Observer
	~Observer () noexcept

	Observer (Observer const &)=delete

	Observer (Observer &&)=delete

Observer &	operator= (Observer const &)=delete

Observer &	operator= (Observer &&)=delete

void	registerProducts (ProductDescriptions &, ModuleDescription const &)

void	fillDescriptions (ModuleDescription const &)

std::string const &	processName () const

bool	wantAllEvents () const

bool	wantEvent (Event const &e)

fhicl::ParameterSetID	selectorConfig () const

Handle< TriggerResults >	getTriggerResults (Event const &e) const

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)

Public Attributes
TH2F *	fPrimaryMom_vs_TrackLength

Private Member Functions
void	ResetVars ()

Private Attributes
TTree *	fTree

int	run

int	subrun

int	event

double	evttime

double	efield [3]

int	t0

int	nhits

int	hit_cryostat [kMaxHits]

int	hit_tpc [kMaxHits]

int	hit_plane [kMaxHits]

int	hit_wire [kMaxHits]

int	hit_channel [kMaxHits]

double	hit_peakT [kMaxHits]

double	hit_charge [kMaxHits]

double	hit_ph [kMaxHits]

std::vector< int >	gap1

std::vector< int >	gap2

std::vector< int >	gap3

std::vector< int >	gap4

std::vector< int >	gap5

std::vector< int >	cross12

std::vector< int >	cross34

std::vector< int >	cross14

std::vector< int >	cross23

std::vector< int >	cross1and2or4

std::vector< int >	cross2and1or3

std::vector< int >	cross1or3and2or4

bool	TPC1 = false

bool	TPC7 = false

bool	TPC5GAP1 = false

bool	GAP1 = false

bool	TPC5GAP2 = false

bool	GAP2 = false

bool	TPC3GAP3 = false

bool	GAP3 = false

bool	TPC3GAP4 = false

bool	GAP4 = false

bool	TPC3 = false

bool	TPC5 = false

bool	GAP5 = false

std::string	fHitsModuleLabel

std::string	fSimulationProducerLabel

Additional Inherited Members
Public Types inherited from art::EDAnalyzer
using	WorkerType = WorkerT< EDAnalyzer >

using	ModuleType = EDAnalyzer

Protected Member Functions inherited from art::Observer
	Observer (fhicl::ParameterSet const &config)

	Observer (std::vector< std::string > const &paths, fhicl::ParameterSet const &config)

detail::ProcessAndEventSelectors &	processAndEventSelectors ()

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 60 of file GapFilter_module.cc.

Constructor & Destructor Documentation

MyGapFilter::GapFilter::GapFilter ( fhicl::ParameterSet const & p )

explicit

Definition at line 145 of file GapFilter_module.cc.

     : EDAnalyzer(pset)
   {
     // Read in the parameters from the .fcl file.
     this->reconfigure(pset);
   }

MyGapFilter::GapFilter::~GapFilter ( )

virtual

Definition at line 162 of file GapFilter_module.cc.

   {
     // Clean up dynamic memory and other resources here.
   }

Member Function Documentation

void MyGapFilter::GapFilter::analyze ( const art::Event & evt )

Definition at line 167 of file GapFilter_module.cc.

   {
     // Implementation of required member function here.
     ResetVars();
 
     art::ServiceHandle<geo::Geometry> geom;
     
     run = evt.run();
     subrun = evt.subRun();
     event = evt.id().event();
 
     //Reset vectors per event
     gap1.clear();
     gap2.clear();
     gap3.clear();
     gap4.clear();
     gap5.clear();
     cross12.clear();
     cross34.clear();
     cross14.clear();
     cross23.clear();
     cross1and2or4.clear();
     cross2and1or3.clear();
     cross1or3and2or4.clear();
 
     //Art handles
 
     art::Handle< std::vector<simb::MCParticle> > particleHandle;
     evt.getByLabel(fSimulationProducerLabel, particleHandle);
     art::Handle< std::vector<sim::SimChannel> > simChannelHandle;
     evt.getByLabel(fSimulationProducerLabel, simChannelHandle);
     std::map< int, const simb::MCParticle* > particleMap;
 
     // Get hit list
     art::Handle< std::vector<recob::Hit> > hitListHandle;
     std::vector<art::Ptr<recob::Hit> > hitlist;
     if (evt.getByLabel(fHitsModuleLabel,hitListHandle))
       art::fill_ptr_vector(hitlist, hitListHandle);
     nhits = hitlist.size();
   
     //loop to determine hits and assign bool metrics for deciding on crossers
     for (size_t k = 0; k < hitlist.size(); ++k){          
       unsigned int channel = hitlist[k]->Channel();
       geo::WireID wireid   = hitlist[k]->WireID();
       hit_cryostat[k]      = wireid.Cryostat;
       hit_tpc[k]           = wireid.TPC;
       hit_plane[k]         = wireid.Plane;
       hit_wire[k]          = wireid.Wire;
       hit_channel[k]       = channel;
       hit_peakT[k]         = hitlist[k]->PeakTime();
       hit_charge[k]        = hitlist[k]->Integral();
       hit_ph[k]            = hitlist[k]->PeakAmplitude();
    
       if (hit_channel[k] == 511  && hit_charge[k] != 0) TPC1                           = true;
       if (hit_channel[k] == 1424 && hit_charge[k] != 0) TPC5GAP1                       = true;
       if (hit_channel[k] == 1535 && hit_charge[k] != 0) TPC5GAP2                       = true;
       if (hit_channel[k] == 912  && hit_charge[k] != 0) TPC3GAP3                       = true;
       if (hit_channel[k] == 1023 && hit_charge[k] != 0) TPC3GAP4                       = true;
       if (hit_channel[k] == 1936 && hit_charge[k] != 0) TPC7                           = true;
       if (hit_channel[k] >= 912  && hit_channel[k] <= 1023 && hit_charge[k] != 0) TPC3 = true;  
       if (hit_channel[k] >= 1424 && hit_channel[k] <= 1535 && hit_charge[k] != 0) TPC5 = true;  
     }  
 
     //Use bools to decide whether an event crosses a gap
     
     if (TPC1 == true && TPC5GAP1 == true) {
       GAP1 = true;
       gap1.push_back(event);
     }
     else GAP1 = false;
         
     if (TPC5GAP2 == true && TPC7 == true) {
       GAP2 = true;
       gap2.push_back(event);
     }
     else GAP2 = false;
         
     if (TPC1 == true && TPC3GAP3 == true) {
       GAP3 = true;
       gap3.push_back(event);
     }
     else GAP3 = false;  
 
     if (TPC3GAP4 == true && TPC7 == true) {
       GAP4 = true;
       gap4.push_back(event);
     }
     else GAP4 = false;
 
     if (TPC3 == true && TPC5 == true) {
       GAP5 = true;
       gap5.push_back(event);
     }
     else GAP5 = false;
     
     //Use gap crossers to determine multiple gap crossing events
     if (GAP1 == true && GAP2 == true) cross12.push_back(event);
     if (GAP3 == true && GAP4 == true) cross34.push_back(event);
     if (GAP1 == true && GAP4 == true) cross14.push_back(event);
     if (GAP2 == true && GAP3 == true) cross23.push_back(event);
     if (GAP1 == true && (GAP2 == true || GAP4 == true)) cross1and2or4.push_back(event);
     if (GAP2 == true && (GAP1 == true || GAP3 == true)) cross2and1or3.push_back(event);
     if ((GAP1 == true || GAP3 == true) && (GAP2 == true || GAP4 == true)) cross1or3and2or4.push_back(event);
 
     //Reset bools for next event
 
     TPC1     = false; 
     TPC7     = false;
     TPC5GAP1 = false;
     GAP1     = false;    
     TPC5GAP2 = false;
     GAP2     = false;    
     TPC3GAP3 = false;
     GAP3     = false; 
     TPC3GAP4 = false;
     GAP4     = false;
     TPC3     = false;
     TPC5     = false;
     GAP5     = false;
   
     fTree->Fill();
   }

void MyGapFilter::GapFilter::beginJob ( )

virtual

Reimplemented from art::EDAnalyzer.

Definition at line 290 of file GapFilter_module.cc.

   {    
     // Implementation of optional member function here.
     art::ServiceHandle<art::TFileService> tfs;
   
     fTree = tfs->make<TTree>("hitdumper","analysis tree");
     fTree->Branch("run",&run,"run/I");
     fTree->Branch("subrun",&subrun,"subrun/I");
     fTree->Branch("event",&event,"event/I");
     fTree->Branch("evttime",&evttime,"evttime/D");
     fTree->Branch("efield",efield,"efield[3]/D");
     fTree->Branch("t0",&t0,"t0/I");
 
     fTree->Branch("nhits",&nhits,"nhits/I");
     fTree->Branch("hit_cryostat",hit_cryostat,"hit_cryostat[nhits]/I");
     fTree->Branch("hit_tpc",hit_tpc,"hit_tpc[nhits]/I");
     fTree->Branch("hit_plane",hit_plane,"hit_plane[nhits]/I");
     fTree->Branch("hit_wire",hit_wire,"hit_wire[nhits]/I");
     fTree->Branch("hit_channel",hit_channel,"hit_channel[nhits]/I");
     fTree->Branch("hit_peakT",hit_peakT,"hit_peakT[nhits]/D");
     fTree->Branch("hit_charge",hit_charge,"hit_charge[nhits]/D");
     fTree->Branch("hit_ph",hit_ph,"hit_ph[nhits]/D");
 
     fTree->Branch("gap1", &gap1);
     fTree->Branch("gap2", &gap2);
     fTree->Branch("gap3", &gap3);
     fTree->Branch("gap4", &gap4);
     fTree->Branch("gap5", &gap5);
 
     fTree->Branch("cross12", &cross12);
     fTree->Branch("cross34", &cross34);
 
     fTree->Branch("cross14", &cross14);
     fTree->Branch("cross23", &cross23);
 
     fTree->Branch("cross1and2or4", &cross1and2or4);
     fTree->Branch("cross2and1or3", &cross2and1or3);
     fTree->Branch("cross1or3and2or4", &cross1or3and2or4);
   }

void MyGapFilter::GapFilter::reconfigure ( fhicl::ParameterSet const & pset )

Definition at line 152 of file GapFilter_module.cc.

   {
     // Read parameters from the .fcl file. The names in the arguments
     // to p.get<TYPE> must match names in the .fcl file.
 
     fHitsModuleLabel         = p.get< std::string >("HitsModuleLabel");
     fSimulationProducerLabel = p.get< std::string >("SimulationProducerLabel");
     return;
   }

void MyGapFilter::GapFilter::ResetVars ( void )

private

Definition at line 330 of file GapFilter_module.cc.

                            {
 
     run = -99999;
     subrun = -99999;
     event = -99999;
     evttime = -99999;
     for (int i = 0; i<3; ++i){
       efield[i] = -99999;
     }
     t0 = -99999;
     nhits = -99999;
     for (int i = 0; i<kMaxHits; ++i){
       hit_cryostat[i] = -99999;
       hit_tpc[i] = -99999;
       hit_plane[i] = -99999;
       hit_wire[i] = -99999;
       hit_channel[i] = -99999;
       hit_peakT[i] = -99999;
       hit_charge[i] = -99999;
       hit_ph[i] = -99999;
     }
   }

Member Data Documentation

std::vector<int> MyGapFilter::GapFilter::cross12

private

Definition at line 106 of file GapFilter_module.cc.

std::vector<int> MyGapFilter::GapFilter::cross14

private

Definition at line 110 of file GapFilter_module.cc.

std::vector<int> MyGapFilter::GapFilter::cross1and2or4

private

Definition at line 114 of file GapFilter_module.cc.

std::vector<int> MyGapFilter::GapFilter::cross1or3and2or4

private

Definition at line 116 of file GapFilter_module.cc.

std::vector<int> MyGapFilter::GapFilter::cross23

private

Definition at line 111 of file GapFilter_module.cc.

std::vector<int> MyGapFilter::GapFilter::cross2and1or3

private

Definition at line 115 of file GapFilter_module.cc.

std::vector<int> MyGapFilter::GapFilter::cross34

private

Definition at line 107 of file GapFilter_module.cc.

double MyGapFilter::GapFilter::efield[3]

private

Definition at line 84 of file GapFilter_module.cc.

int MyGapFilter::GapFilter::event

private

Definition at line 82 of file GapFilter_module.cc.

double MyGapFilter::GapFilter::evttime

private

Definition at line 83 of file GapFilter_module.cc.

std::string MyGapFilter::GapFilter::fHitsModuleLabel

private

Definition at line 139 of file GapFilter_module.cc.

TH2F* MyGapFilter::GapFilter::fPrimaryMom_vs_TrackLength

Definition at line 71 of file GapFilter_module.cc.

std::string MyGapFilter::GapFilter::fSimulationProducerLabel

private

Definition at line 140 of file GapFilter_module.cc.

TTree* MyGapFilter::GapFilter::fTree

private

Definition at line 77 of file GapFilter_module.cc.

std::vector<int> MyGapFilter::GapFilter::gap1

private

Definition at line 99 of file GapFilter_module.cc.

bool MyGapFilter::GapFilter::GAP1 = false

private

Definition at line 123 of file GapFilter_module.cc.

std::vector<int> MyGapFilter::GapFilter::gap2

private

Definition at line 100 of file GapFilter_module.cc.

bool MyGapFilter::GapFilter::GAP2 = false

private

Definition at line 126 of file GapFilter_module.cc.

std::vector<int> MyGapFilter::GapFilter::gap3

private

Definition at line 101 of file GapFilter_module.cc.

bool MyGapFilter::GapFilter::GAP3 = false

private

Definition at line 129 of file GapFilter_module.cc.

std::vector<int> MyGapFilter::GapFilter::gap4

private

Definition at line 102 of file GapFilter_module.cc.

bool MyGapFilter::GapFilter::GAP4 = false

private

Definition at line 132 of file GapFilter_module.cc.

std::vector<int> MyGapFilter::GapFilter::gap5

private

Definition at line 103 of file GapFilter_module.cc.

bool MyGapFilter::GapFilter::GAP5 = false

private

Definition at line 136 of file GapFilter_module.cc.

int MyGapFilter::GapFilter::hit_channel[kMaxHits]

private

Definition at line 93 of file GapFilter_module.cc.

double MyGapFilter::GapFilter::hit_charge[kMaxHits]

private

Definition at line 95 of file GapFilter_module.cc.

int MyGapFilter::GapFilter::hit_cryostat[kMaxHits]

private

Definition at line 89 of file GapFilter_module.cc.

double MyGapFilter::GapFilter::hit_peakT[kMaxHits]

private

Definition at line 94 of file GapFilter_module.cc.

double MyGapFilter::GapFilter::hit_ph[kMaxHits]

private

Definition at line 96 of file GapFilter_module.cc.

int MyGapFilter::GapFilter::hit_plane[kMaxHits]

private

Definition at line 91 of file GapFilter_module.cc.

int MyGapFilter::GapFilter::hit_tpc[kMaxHits]

private

Definition at line 90 of file GapFilter_module.cc.

int MyGapFilter::GapFilter::hit_wire[kMaxHits]

private

Definition at line 92 of file GapFilter_module.cc.

int MyGapFilter::GapFilter::nhits

private

Definition at line 88 of file GapFilter_module.cc.

int MyGapFilter::GapFilter::run

private

Definition at line 80 of file GapFilter_module.cc.

int MyGapFilter::GapFilter::subrun

private

Definition at line 81 of file GapFilter_module.cc.

int MyGapFilter::GapFilter::t0

private

Definition at line 85 of file GapFilter_module.cc.

bool MyGapFilter::GapFilter::TPC1 = false

private

Definition at line 119 of file GapFilter_module.cc.

bool MyGapFilter::GapFilter::TPC3 = false

private

Definition at line 134 of file GapFilter_module.cc.

bool MyGapFilter::GapFilter::TPC3GAP3 = false

private

Definition at line 128 of file GapFilter_module.cc.

bool MyGapFilter::GapFilter::TPC3GAP4 = false

private

Definition at line 131 of file GapFilter_module.cc.

bool MyGapFilter::GapFilter::TPC5 = false

private

Definition at line 135 of file GapFilter_module.cc.

bool MyGapFilter::GapFilter::TPC5GAP1 = false

private

Definition at line 122 of file GapFilter_module.cc.

bool MyGapFilter::GapFilter::TPC5GAP2 = false

private

Definition at line 125 of file GapFilter_module.cc.

bool MyGapFilter::GapFilter::TPC7 = false

private

Definition at line 120 of file GapFilter_module.cc.

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

dunetpc/dune/Gaps/GapFilter_module.cc

Public Member Functions

Public Attributes

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation