An example selection analysis. More...

#include <ExampleSelection.h>

Inheritance diagram for ana::ExampleAnalysis::ExampleSelection:

Public Member Functions
	ExampleSelection ()

void	Initialize (fhicl::ParameterSet *config=NULL)

void	Finalize ()

bool	ProcessEvent (const gallery::Event &ev, const std::vector< event::Interaction > &truth, std::vector< event::RecoInteraction > &reco)

Public Member Functions inherited from core::SelectionBase
	SelectionBase ()

virtual	~SelectionBase ()

Public Member Functions inherited from core::ProcessorBase
	ProcessorBase ()

virtual	~ProcessorBase ()

virtual void	FillTree ()

virtual void	FillRecoTree ()

virtual void	EventCleanup ()

template<class T >
TBranch *	AddBranch (std::string name, T *obj)

template<class T >
TBranch *	AddRecoBranch (std::string name, T *obj)

Protected Attributes
unsigned	fEventCounter
	Count processed events. More...

art::InputTag	fTruthTag
	art tag for MCTruth information More...

int	fMyParam
	A parameter from the configuration file. More...

int	fNuCount
	Number of neutrino interactions in the event. More...

int	fMyVar
	Another variable of interest. More...

TH2D *	fNuVertexXZHist
	Neutrino vertex XZ projection. More...

Protected Attributes inherited from core::ProcessorBase
unsigned long	fEventIndex
	An incrementing index. More...

Experiment	fExperimentID
	Experiment identifier. More...

ProviderManager *	fProviderManager
	Interface for provider access. More...

std::string	fOutputFilename
	The output filename. More...

std::string	fProviderConfig
	A custom provider config fcl file. More...

std::vector< geo::BoxBoundedGeo >	fActiveVolumes
	List of active volumes in configured detector. More...

bool	fWriteTree
	Enable writing of the main tree. More...

TFile *	fOutputFile
	The output ROOT file. More...

TTree *	fTree
	The output ROOT tree. More...

event::Event *	fEvent
	The standard output event data structure. More...

bool	fWriteRecoTree
	Enable writing of the reco tree. More...

TTree *	fRecoTree
	The output reco ROOT tree. More...

event::RecoEvent *	fRecoEvent
	The standard output reco event data structure. More...

TTree *	fSubRunTree
	Subrun output tree. More...

SubRun *	fSubRun
	Standard output subrun structure. More...

TTree *	fFileMetaTree
	File metadata output tree. More...

FileMeta *	fFileMeta
	standard output file metadata structure More...

TParameter< int > *	fExperimentParameter
	Saves value of experiment enum. More...

std::set< std::pair< int, int > >	fSubRunCache
	Cache stored subruns. More...

art::InputTag	fTruthTag
	art tag for MCTruth information More...

art::InputTag	fFluxTag
	art tag for MCFlux information More...

std::vector< art::InputTag >	fWeightTags
	art tag(s) for MCEventWeight information More...

art::InputTag	fMCTrackTag
	art tag for MCTrack More...

art::InputTag	fMCShowerTag
	art tag for MCShower More...

art::InputTag	fMCParticleTag
	art tag for MCParticle More...

std::string	fGeneratorProcess
	process_name of process used to run genie. Used to extract subrun/POT information More...

Additional Inherited Members
Public Attributes inherited from core::ProcessorBase
std::vector < event::RecoInteraction > *	fReco
	Reco interaction list. More...

Protected Member Functions inherited from core::ProcessorBase
virtual void	Initialize (char *config=NULL)

virtual void	Setup (char *config=NULL)

virtual void	Setup (fhicl::ParameterSet *config=NULL)

virtual void	Teardown ()

void	BuildEventTree (gallery::Event &ev)

void	SetupServices (gallery::Event &ev)

void	UpdateSubRuns (gallery::Event &ev)

void	UpdateFileMeta (gallery::Event &ev)

Detailed Description

An example selection analysis.

This selection analysis doesn't actually select events, it just demonstrates the framework!

Definition at line 38 of file ExampleSelection.h.

Constructor & Destructor Documentation

ana::ExampleAnalysis::ExampleSelection::ExampleSelection ( )

Constructor.

Definition at line 19 of file ExampleSelection.cxx.

20 : SelectionBase(), fNuCount(0), fEventCounter(0) {}

core::SelectionBase::SelectionBase

SelectionBase()

Definition: SelectionBase.cxx:5

ana::ExampleAnalysis::ExampleSelection::fEventCounter

unsigned fEventCounter

Count processed events.

Definition: ExampleSelection.h:68

ana::ExampleAnalysis::ExampleSelection::fNuCount

int fNuCount

Number of neutrino interactions in the event.

Definition: ExampleSelection.h:75

Member Function Documentation

void ana::ExampleAnalysis::ExampleSelection::Finalize ( )

virtual

Finalize and write objects to the output file.

Implements core::ProcessorBase.

Definition at line 58 of file ExampleSelection.cxx.

                                 {
   // Output our histograms to the ROOT file
   fOutputFile->cd();
   fNuVertexXZHist->Write();
 }

void ana::ExampleAnalysis::ExampleSelection::Initialize ( fhicl::ParameterSet * config = NULL )

virtual

Initialization.

Here we load configuration parameters, set up histograms for output, and add our own branches to the output tree.

Parameters

config A configuration, as FHiCL ParameterSet object

Implements core::ProcessorBase.

Definition at line 23 of file ExampleSelection.cxx.

                                                            {
   // Make a histogram
   fNuVertexXZHist = new TH2D("nu_vtx_XZ", "",
                              100, -1000, 1000, 100, -1000, 1000);
 
   // Load configuration parameters
   fMyParam = 0;
   fTruthTag = { "generator" };
 
   if (config) {
     fhicl::ParameterSet pconfig = \
       config->get<fhicl::ParameterSet>("ExampleAnalysis");
 
     fMyParam = pconfig.get<int>("parameter", 0);
 
     fTruthTag = \
       { pconfig.get<std::string>("MCTruthTag", "generator") };
   }
 
   // Add custom branches
   AddBranch("nucount", &fNuCount);
   AddBranch("myvar", &fMyVar);
 
   // Use LArSoft service functions via the ProviderManager
   if (fProviderManager) {
     std::cout << "Detector: "
               << fProviderManager->GetGeometryProvider()->DetectorName()
               << std::endl;
   }
 
   // Use some library code
   hello();
 }

bool ana::ExampleAnalysis::ExampleSelection::ProcessEvent	(	const gallery::Event &	ev,
		const std::vector< event::Interaction > &	truth,
		std::vector< event::RecoInteraction > &	reco
	)

virtual

Process one event.

Parameters

ev	A single event, as a gallery::Event
reco	Reconstructed interactions

Returns: True to keep event

Implements core::ProcessorBase.

Definition at line 65 of file ExampleSelection.cxx.

                                            {
   if (fEventCounter % 10 == 0) {
     std::cout << "ExampleSelection: Processing event "
               << fEventCounter << std::endl;
   }
   fEventCounter++;
 
   // Grab a data product from the event
   auto const& mctruths = \
     *ev.getValidHandle<std::vector<simb::MCTruth>>(fTruthTag);
 
   // Fill in the custom branches
   fNuCount = mctruths.size();  // Number of neutrinos in this event
   fMyVar = fMyParam;
 
   // Iterate through the neutrinos
   for (size_t i=0; i<mctruths.size(); i++) {
     auto const& mctruth = mctruths.at(i);
 
     // Fill neutrino vertex position histogram
     if (mctruth.NeutrinoSet()) {
       fNuVertexXZHist->Fill(mctruth.GetNeutrino().Nu().Vx(),
                             mctruth.GetNeutrino().Nu().Vz());
     }
     // Add in the "reconstructed" interaction
     //
     // Contruct truth information from the provided vector
     event::RecoInteraction interaction(i);
     // get "reconstructed" energy
     interaction.reco_energy = util::ECCQE(mctruth.GetNeutrino().Nu().Momentum().Vect(), mctruth.GetNeutrino().Lepton().E());
     reco.push_back(interaction);
   }
 
   return true;
 }

Member Data Documentation

unsigned ana::ExampleAnalysis::ExampleSelection::fEventCounter

protected

Count processed events.

Definition at line 68 of file ExampleSelection.h.

int ana::ExampleAnalysis::ExampleSelection::fMyParam

protected

A parameter from the configuration file.

Definition at line 72 of file ExampleSelection.h.

int ana::ExampleAnalysis::ExampleSelection::fMyVar

protected

Another variable of interest.

Definition at line 76 of file ExampleSelection.h.

int ana::ExampleAnalysis::ExampleSelection::fNuCount

protected

Number of neutrino interactions in the event.

Custom data branches

Definition at line 75 of file ExampleSelection.h.

TH2D* ana::ExampleAnalysis::ExampleSelection::fNuVertexXZHist

protected

Neutrino vertex XZ projection.

Histograms

Definition at line 79 of file ExampleSelection.h.

art::InputTag ana::ExampleAnalysis::ExampleSelection::fTruthTag

protected

art tag for MCTruth information

Configuration parameters

Definition at line 71 of file ExampleSelection.h.

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

Public Member Functions

Protected Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation