Electron neutrino event selection. More...

Inheritance diagram for ana::SBNOsc::CosmicArrivalTimes:

Public Member Functions
	CosmicArrivalTimes ()

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

void	Finalize ()

bool	ProcessEvent (const gallery::Event &event, 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
std::string	fMCParticleTag

unsigned	event_ind

TH1D *	fCosmicEnterTimes

TH2D *	fCosmicEnterTimesPE

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

Electron neutrino event selection.

Definition at line 39 of file CosmicArrivalTimes.cc.

Constructor & Destructor Documentation

ana::SBNOsc::CosmicArrivalTimes::CosmicArrivalTimes ( )

inline

Constructor.

Definition at line 42 of file CosmicArrivalTimes.cc.

42 {}

Member Function Documentation

void ana::SBNOsc::CosmicArrivalTimes::Finalize ( )

inlinevirtual

Finalize and write objects to the output file.

Implements core::ProcessorBase.

Definition at line 57 of file CosmicArrivalTimes.cc.

                   {
     fOutputFile->cd();
     fCosmicEnterTimes->Write();
     fCosmicEnterTimesPE->Write();
   }

void ana::SBNOsc::CosmicArrivalTimes::Initialize ( fhicl::ParameterSet * config = NULL )

inlinevirtual

Initialization.

Parameters

config A configuration, as a FHiCL ParameterSet object

Implements core::ProcessorBase.

Definition at line 49 of file CosmicArrivalTimes.cc.

                                                   {
     fMCParticleTag = config ? config->get<std::string>("MCParticleTag", "largeant") : "largeant";
     fCosmicEnterTimes = new TH1D("cosmic_enter", "cosmic_enter", 700, -60., 10.); 
     fCosmicEnterTimesPE = new TH2D("cosmic_enter_v_pe", "cosmic_enter_v_pe", 700, -60., 10., 100, 0, 10000);
     event_ind = 0;
   }

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

inlinevirtual

Process one event.

Parameters

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

Returns: True to keep event

Implements core::ProcessorBase.

Definition at line 70 of file CosmicArrivalTimes.cc.

                                                                                                                                 {
     std::cout << "New Event!\n";
     auto const &intime = event.getValidHandle<std::vector<simb::MCTruth>>({"GenInTimeSorter", "intime"});
     //auto const &outtime = event.getValidHandle<std::vector<simb::MCTruth>>({"GenInTimeSorter", "outtime"});
 
     art::FindManyP<simb::MCParticle, sim::GeneratedParticleInfo> intime_to_particles(intime, event, fMCParticleTag);
     //art::FindManyP<simb::MCParticle, sim::GeneratedParticleInfo> outtime_to_particles(outtime, event, fMCParticleTag);
 
     double in_time_photon_energy = 0.;
     unsigned in_time_photons = 0;
     const std::vector<sim::SimPhotons> &photon_list = *event.getValidHandle<std::vector<sim::SimPhotons>>("larg4intime");
     for (const sim::SimPhotons &photons: photon_list) {
       for (const sim::OnePhoton &photon: photons) {
         if (photon.Time > -200. && photon.Time < 1800.) {
           in_time_photon_energy += photon.Energy;
           in_time_photons ++;
         }
       }
     }
     gallery::Handle<std::vector<sim::SimPhotons>> reflected;
     event.getByLabel({"larg4intime", "Reflected"}, reflected);
     if (reflected.isValid()) {
     std::cout << "Reflected!\n";
     const std::vector<sim::SimPhotons> &photon_list = *event.getValidHandle<std::vector<sim::SimPhotons>>({"larg4intime", "Reflected"});
     for (const sim::SimPhotons &photons: photon_list) {
       for (const sim::OnePhoton &photon: photons) {
         if (photon.Time > -200. && photon.Time < 1800.) {
           in_time_photon_energy += photon.Energy;
           in_time_photons ++;
         }
       }
     }
     }
 
     for (const art::Ptr<simb::MCParticle> part: intime_to_particles.at(0)) {
       if (abs(part->PdgCode()) != 13) continue;
       unsigned n_traj = part->NumberTrajectoryPoints();
       bool next = false;
       for (unsigned i = 0; i < n_traj; i++) {
         for (const geo::BoxBoundedGeo &vol: fActiveVolumes) {
           if (vol.ContainsPosition(part->Position(i).Vect())) {
             fCosmicEnterTimes->Fill(part->Position(i).T() / 1000.) /* ns -> us*/;
             fCosmicEnterTimesPE->Fill(part->Position(i).T() / 1000., in_time_photons); 
             next = true; 
             break;
           }
         }
         if (next) break;
       }
     }
 
     /*
     for (const art::Ptr<simb::MCParticle> part: outtime_to_particles.at(0)) {
       unsigned n_traj = part->NumberTrajectoryPoints();
       bool next = false;
       for (unsigned i = 0; i < n_traj; i++) {
         for (const geo::BoxBoundedGeo &vol: fActiveVolumes) {
           if (vol.ContainsPosition(part->Position(i).Vect())) {
             fCosmicEnterTimes->Fill(part->Position(i).T() / 1000.);
             next = true; 
             break;
           }
         }
         if (next) break;
       }
     }*/
 
     return false; 
   }

Member Data Documentation

unsigned ana::SBNOsc::CosmicArrivalTimes::event_ind

protected

Definition at line 142 of file CosmicArrivalTimes.cc.

TH1D* ana::SBNOsc::CosmicArrivalTimes::fCosmicEnterTimes

protected

Definition at line 143 of file CosmicArrivalTimes.cc.

TH2D* ana::SBNOsc::CosmicArrivalTimes::fCosmicEnterTimesPE

protected

Definition at line 144 of file CosmicArrivalTimes.cc.

std::string ana::SBNOsc::CosmicArrivalTimes::fMCParticleTag

protected

Definition at line 141 of file CosmicArrivalTimes.cc.

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

CosmicArrivalTimes.cc

Public Member Functions

Protected Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation