Inheritance diagram for protoana::ProtoDUNETruthBeamParticle:

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

virtual	~ProtoDUNETruthBeamParticle ()

void	beginJob ()

void	analyze (const art::Event &evt)

void	reconfigure (fhicl::ParameterSet const &pset)

void	endJob ()

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

fhicl::ParameterSetID	selectorConfig () 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)

Private Attributes
int	nbeamparticles

int	run [NMAXBEAMPARTICLES]

int	subrun [NMAXBEAMPARTICLES]

int	event [NMAXBEAMPARTICLES]

int	pdg [NMAXBEAMPARTICLES]

int	goodparticle [NMAXBEAMPARTICLES]

float	mom [NMAXBEAMPARTICLES]

float	ene [NMAXBEAMPARTICLES]

float	startpos [NMAXBEAMPARTICLES][4]

float	p [NMAXBEAMPARTICLES][3]

TTree *	truth

std::vector< int >	p_pdg

bool	p_savegoodparticle

bool	selectall

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

using	ModuleType = EDAnalyzer

Protected Member Functions inherited from art::Observer
std::string const &	processName () const

bool	wantAllEvents () const noexcept

bool	wantEvent (ScheduleID id, Event const &e) const

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

	Observer (fhicl::ParameterSet const &config)

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

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 41 of file ProtoDUNETruthBeamParticle_module.cc.

Constructor & Destructor Documentation

protoana::ProtoDUNETruthBeamParticle::ProtoDUNETruthBeamParticle ( fhicl::ParameterSet const & pset )

explicit

Definition at line 75 of file ProtoDUNETruthBeamParticle_module.cc.

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

protoana::ProtoDUNETruthBeamParticle::~ProtoDUNETruthBeamParticle ( )

virtual

Definition at line 82 of file ProtoDUNETruthBeamParticle_module.cc.

82 {}

Member Function Documentation

void protoana::ProtoDUNETruthBeamParticle::analyze ( const art::Event & evt )

Definition at line 117 of file ProtoDUNETruthBeamParticle_module.cc.

                                                                    {
     
   // Access truth info
   std::vector< art::Ptr<simb::MCTruth> > mclist;
   auto mctruthhandle = evt.getHandle< std::vector<simb::MCTruth> >("generator");
   if (mctruthhandle){ 
     art::fill_ptr_vector(mclist,mctruthhandle);
   }
   else{
     mf::LogError("protoana::ProtoDUNETruthBeamParticle::analyze") << "Requested protoDUNE beam generator information does not exist!";
     return;
   }
   
   art::Ptr<simb::MCTruth> mctruth;
   mctruth = mclist[0];
   
   nbeamparticles = 0;
   for(int i = 0; i < mctruth->NParticles(); ++i){
     const simb::MCParticle& part(mctruth->GetParticle(i));
 
     // Select partciles with the chosen pdg
     if(!selectall){
       bool found = false;
       for(unsigned int j = 0; j < p_pdg.size(); j++){
         if(part.PdgCode() == p_pdg[j]){
           found = true;
           break;
         }
       }
       if(!found) continue;
     }   
     
     // Option to select only good beam particles
     bool isgoodparticle   = (part.Process() == "primary");
     if(p_savegoodparticle && !isgoodparticle) continue;
 
     nbeamparticles++;
 
     run[i]            = evt.run();
     subrun[i]         = evt.subRun();
     event[i]          = evt.id().event();
     goodparticle[i]   = (part.Process() == "primary");
     startpos[i][0]    = part.Vx();
     startpos[i][1]    = part.Vy();
     startpos[i][2]    = part.Vz();
     startpos[i][3]    = part.T();
     p[i][0]           = part.Px();
     p[i][1]           = part.Py();
     p[i][2]           = part.Pz();
     mom[i]            = part.P();
     ene[i]            = part.E();
     pdg[i]            = part.PdgCode();
   }
   
   if(nbeamparticles > 0)
     truth->Fill();
   
 }

void protoana::ProtoDUNETruthBeamParticle::beginJob ( )

virtual

Reimplemented from art::EDAnalyzer.

Definition at line 85 of file ProtoDUNETruthBeamParticle_module.cc.

                                                   {
   
   art::ServiceHandle<art::TFileService> tfs;
   
   // Define output tree
   truth = tfs->make<TTree>("truth", "Beam Particle Truth Tree");
   truth->Branch("nbeamparticles",   &nbeamparticles,       "nbeamparticles/I");
   truth->Branch("run",              run,                   "run[nbeamparticles]/I");
   truth->Branch("subrun",           subrun,                "subrun[nbeamparticles]/I");
   truth->Branch("event",            event,                 "event[nbeamparticles]/I");
   truth->Branch("pdg",              pdg,                   "pdg[nbeamparticles]/I");
   truth->Branch("goodparticle",     goodparticle,          "goodparticle[nbeamparticles]/I");
   truth->Branch("mom",              mom,                   "mom[nbeamparticles]/F");
   truth->Branch("ene",              ene,                   "ene[nbeamparticles]/F");
   truth->Branch("startpos",         startpos,              "startpos[nbeamparticles][4]/F");
   truth->Branch("p",                p,                     "p[nbeamparticles][3]/F");
   
 }

void protoana::ProtoDUNETruthBeamParticle::endJob ( )

virtual

Reimplemented from art::EDAnalyzer.

Definition at line 176 of file ProtoDUNETruthBeamParticle_module.cc.

176 {}

void protoana::ProtoDUNETruthBeamParticle::reconfigure ( fhicl::ParameterSet const & pset )

Definition at line 105 of file ProtoDUNETruthBeamParticle_module.cc.

                                                                                  {
 
   p_savegoodparticle  = pset.get<bool>("SaveOnlyGoodParticle");
   p_pdg               = pset.get< std::vector<int> >("Pdg");
   if(p_pdg[0] == 0)
     selectall = true;
   else
     selectall = false;
 
 }