Inheritance diagram for beam::ReadBeamInfo:

Classes
struct	Config

Public Types
using	Parameters = art::EDAnalyzer::Table< Config >

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

using	ModuleType = EDAnalyzer

Public Member Functions
	ReadBeamInfo (Parameters const &pset)

virtual	~ReadBeamInfo ()

void	analyze (const art::Event &event) override

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
art::InputTag	fBeamLabel
	The name of the beam event producer. More...

Additional Inherited Members
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 34 of file ReadBeamInfo_module.cc.

Member Typedef Documentation

using beam::ReadBeamInfo::Parameters = art::EDAnalyzer::Table<Config>

Definition at line 48 of file ReadBeamInfo_module.cc.

Constructor & Destructor Documentation

beam::ReadBeamInfo::ReadBeamInfo ( ReadBeamInfo::Parameters const & pset )

explicit

Definition at line 66 of file ReadBeamInfo_module.cc.

                                                                 : art::EDAnalyzer(pset){
 
 //  fBeamLabel = pset.get<art::InputTag>("BeamLabel");
   fBeamLabel = pset().BeamLabel();
 
 }

beam::ReadBeamInfo::~ReadBeamInfo ( )

virtual

Definition at line 73 of file ReadBeamInfo_module.cc.

                                {
 
 }

Member Function Documentation

void beam::ReadBeamInfo::analyze ( const art::Event & event )

override

Definition at line 77 of file ReadBeamInfo_module.cc.

                                                  {
 
   // Just for MC!
   if(evt.isRealData()) return;
 
   // ProtoDUNE beam generator information
 
   std::vector<art::Ptr<simb::MCTruth> > beamTruth;
   auto beamTruthListHandle = evt.getHandle< std::vector<simb::MCTruth> >(fBeamLabel);
   if (beamTruthListHandle){
     art::fill_ptr_vector(beamTruth,beamTruthListHandle);
   }
   else{
     // If we requested this info but it doesn't exist, don't use it.
     mf::LogError("ReadBeamInfo") << "Requested protoDUNE beam generator information with name " << fBeamLabel << " but none exists";
     return;
   }
 
   // There is only one MC truth, so use it to get the number of primaries
   Int_t nPrimaries = beamTruth[0]->NParticles();
  
   for(Int_t iPartp = 0; iPartp < nPrimaries; ++iPartp){
     const simb::MCParticle& partp(beamTruth[0]->GetParticle(iPartp));
 
     // Access the beam particle variables
     TVector3 beamVtx(partp.Vx(), partp.Vy(), partp.Vz());
     TVector3 beamMom(partp.Px(), partp.Py(), partp.Pz());
     TVector3 beamDir = beamMom.Unit();
     double beamTime = partp.T();
     double beamMomentum = partp.P();
     double beamEnergy = partp.E();
     int beamPDG = partp.PdgCode();
 
     // There are two types of beam particles...
     // The ones we trigger on, "GoodParticles" have Process() == "primary"
     // The background ones have Process() == "primaryBackground"
     // This is how we identify which ones are the good particle.
     bool beamIsGood = partp.Process()=="primary";
 
     std::string type = "good";
     if(!beamIsGood){
       type = "background";
     }
 
     std::cout << " -- Read " << type << " beam particle with pdg code " << beamPDG << " from event " << evt.event() << std::endl;
 
     if(beamIsGood){
       std::cout << " - Vtx    = (" << beamVtx.X() << ", " << beamVtx.Y() << ", " << beamVtx.Z() << ")" << std::endl; 
       std::cout << " - Time   = " << beamTime << std::endl;
       std::cout << " - Dir    = (" << beamDir.X() << ", " << beamDir.Y() << ", " << beamDir.Z() << ")" << std::endl; 
       std::cout << " - Mom    = (" << beamMom.X() << ", " << beamMom.Y() << ", " << beamMom.Z() << ")" << std::endl; 
       std::cout << " - |Mom|  = " << beamMomentum << std::endl;
       std::cout << " - Energy = " << beamEnergy << std::endl;
       std::cout << " - PDG    = " << beamPDG << std::endl;
     }
 
   } // End loop over the beam particles
 
 }

Member Data Documentation

art::InputTag beam::ReadBeamInfo::fBeamLabel

private

The name of the beam event producer.

Definition at line 59 of file ReadBeamInfo_module.cc.

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

protoduneana/protoduneana/TutorialExamples/ReadBeamInfo_module.cc

Classes

Public Types

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation