Inheritance diagram for protoana::BeamExample:

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

	BeamExample (BeamExample const &)=delete

	BeamExample (BeamExample &&)=delete

BeamExample &	operator= (BeamExample const &)=delete

BeamExample &	operator= (BeamExample &&)=delete

virtual void	beginJob () override

virtual void	endJob () override

void	analyze (art::Event const &e) 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
std::string	fCalorimetryTag

std::string	fTrackerTag

std::string	fShowerTag

std::string	fPFParticleTag

std::string	fGeneratorTag

bool	fVerbose

protoana::ProtoDUNEDataUtils	dataUtil

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

Constructor & Destructor Documentation

protoana::BeamExample::BeamExample ( fhicl::ParameterSet const & p )

explicit

Definition at line 84 of file BeamExample_module.cc.

   :
   EDAnalyzer(p),
   fCalorimetryTag(p.get<std::string>("CalorimetryTag")),
   fTrackerTag(p.get<std::string>("TrackerTag")),
   fShowerTag(p.get<std::string>("ShowerTag")),
   fPFParticleTag(p.get<std::string>("PFParticleTag")),
   fGeneratorTag(p.get<std::string>("GeneratorTag")),
   fVerbose(p.get<bool>("Verbose")),
   dataUtil(p.get<fhicl::ParameterSet>("DataUtils"))
 {
 
 }

protoana::BeamExample::BeamExample ( BeamExample const & )

delete

protoana::BeamExample::BeamExample ( BeamExample && )

delete

Member Function Documentation

void protoana::BeamExample::analyze ( art::Event const & e )

overridevirtual

Implements art::EDAnalyzer.

Definition at line 103 of file BeamExample_module.cc.

 {
 
   bool beamTriggerEvent = false;
   // If this event is MC then we can check what the true beam particle is
   if(!evt.isRealData()){
     // Get the truth utility to help us out
     protoana::ProtoDUNETruthUtils truthUtil;
     // Firstly we need to get the list of MCTruth objects from the generator. The standard protoDUNE
     // simulation has fGeneratorTag = "generator"
     auto mcTruths = evt.getValidHandle<std::vector<simb::MCTruth>>(fGeneratorTag);
     // mcTruths is basically a pointer to an std::vector of simb::MCTruth objects. There should only be one
     // of these, so we pass the first element into the function to get the good particle
     const simb::MCParticle* geantGoodParticle = truthUtil.GetGeantGoodParticle((*mcTruths)[0],evt);
     if(geantGoodParticle != 0x0){
       std::cout << "Found GEANT particle corresponding to the good particle with pdg = " << geantGoodParticle->PdgCode() << std::endl;
     }
   }
   else{
     // For data we can see if this event comes from a beam trigger
     beamTriggerEvent = dataUtil.IsBeamTrigger(evt);
     if(beamTriggerEvent){
       std::cout << "This data event has a beam trigger" << std::endl;
     }
   }
 
   /*
   // Now we want to access the output from Pandora. This comes in the form of particle flow objects (recob::PFParticle).
   // The primary PFParticles are those we want to consider and these PFParticles then have a hierarchy of daughters that
   // describe the whole interaction of a given primary particle
   //
   //                     / daughter track
   //                    /
   //  primary track    /   
   //  ---------------- ---- daughter track
   //                   \
   //                   /\-
   //                   /\\-- daughter shower
   //
   // The above primary PFParticle will have links to three daughter particles, two track-like and one shower-like
   */
 
   // Get the PFParticle utility
   protoana::ProtoDUNEPFParticleUtils pfpUtil;
 
   // Get all of the PFParticles, by default from the "pandora" product
   auto recoParticles = evt.getValidHandle<std::vector<recob::PFParticle>>(fPFParticleTag);
 
   // We'd like to find the beam particle. Pandora tries to do this for us, so let's use the PFParticle utility 
   // to look for it. Pandora reconstructs slices containing one (or sometimes more) primary PFParticles. These
   // are tagged as either beam or cosmic for ProtoDUNE. This function automatically considers only those
   // PFParticles considered as primary
   std::vector<const recob::PFParticle*> beamParticles = pfpUtil.GetPFParticlesFromBeamSlice(evt,fPFParticleTag);
 
   if(beamParticles.size() == 0){
     std::cerr << "We found no beam particles for this event... moving on" << std::endl;
     return;
   }
 
   // We can now look at these particles
   for(const recob::PFParticle* particle : beamParticles){
 
     // "particle" is the pointer to our beam particle. The recob::Track or recob::Shower object
     // of this particle might be more helpful. These return null pointers if not track-like / shower-like
     const recob::Track* thisTrack = pfpUtil.GetPFParticleTrack(*particle,evt,fPFParticleTag,fTrackerTag);
     const recob::Shower* thisShower = pfpUtil.GetPFParticleShower(*particle,evt,fPFParticleTag,fShowerTag);
     if(thisTrack != 0x0) std::cout << "Beam particle is track-like" << std::endl;
     if(thisShower != 0x0) std::cout << "Beam particle is shower-like" << std::endl;
 
     // Find the particle vertex. We need the tracker tag here because we need to do a bit of
     // additional work if the PFParticle is track-like to find the vertex. 
     const TVector3 vtx = pfpUtil.GetPFParticleVertex(*particle,evt,fPFParticleTag,fTrackerTag);
 
     // Now we can look for the interaction point of the particle if one exists, i.e where the particle
     // scatters off an argon nucleus. Shower-like objects won't have an interaction point, so we can
     // check this by making sure we get a sensible position
     const TVector3 interactionVtx = pfpUtil.GetPFParticleSecondaryVertex(*particle,evt,fPFParticleTag,fTrackerTag);
 
     // Let's get the daughter PFParticles... we can do this simply without the utility
     for(const int daughterID : particle->Daughters()){
       // Daughter ID is the element of the original recoParticle vector
       const recob::PFParticle *daughterParticle = &(recoParticles->at(daughterID));
       std::cout << "Daughter " << daughterID << " has " << daughterParticle->NumDaughters() << " daughters" << std::endl;
     }
  
     // For actually studying the objects, it is easier to have the daughters in their track and shower forms.
     // We can use the utility to get a vector of track-like and a vector of shower-like daughters
     const std::vector<const recob::Track*> trackDaughters = pfpUtil.GetPFParticleDaughterTracks(*particle,evt,fPFParticleTag,fTrackerTag);  
     const std::vector<const recob::Shower*> showerDaughters = pfpUtil.GetPFParticleDaughterShowers(*particle,evt,fPFParticleTag,fShowerTag);  
     std::cout << "Beam particle has " << trackDaughters.size() << " track-like daughters and " << showerDaughters.size() << " shower-like daughters." << std::endl;
 
     // At this point we have access to the primary particle track or shower, plus the track and shower daughters of this event.
     // For other information that can be extracted from the PFParticle objects please see the ProtoDUNEPFParticleUtils header file
     // /dunetpc/dune/ProtoDUNE/Analysis/ProtoDUNEPFParticleUtils.h
 
   } 
 
 }

void protoana::BeamExample::beginJob ( )

overridevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 98 of file BeamExample_module.cc.

 {
 
 }

void protoana::BeamExample::endJob ( )

overridevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 202 of file BeamExample_module.cc.

 {
 
 }

BeamExample& protoana::BeamExample::operator= ( BeamExample const & )

delete

BeamExample& protoana::BeamExample::operator= ( BeamExample && )

delete

Member Data Documentation

protoana::ProtoDUNEDataUtils protoana::BeamExample::dataUtil

private

Definition at line 79 of file BeamExample_module.cc.

std::string protoana::BeamExample::fCalorimetryTag

private

Definition at line 72 of file BeamExample_module.cc.

std::string protoana::BeamExample::fGeneratorTag

private

Definition at line 76 of file BeamExample_module.cc.

std::string protoana::BeamExample::fPFParticleTag

private

Definition at line 75 of file BeamExample_module.cc.

std::string protoana::BeamExample::fShowerTag

private

Definition at line 74 of file BeamExample_module.cc.

std::string protoana::BeamExample::fTrackerTag

private

Definition at line 73 of file BeamExample_module.cc.

bool protoana::BeamExample::fVerbose

private

Definition at line 77 of file BeamExample_module.cc.

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

protoduneana/protoduneana/Utilities/BeamExample/BeamExample_module.cc

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation