#include "SimPhotonCounter.h"

Inheritance diagram for opdet::SimPhotonCounter:

Public Types
using	ModuleType = EDAnalyzer

Public Member Functions
	SimPhotonCounter ()

	SimPhotonCounter (size_t s, float t_p1, float t_p2, float t_l1, float t_l2, float min_w=0, float max_w=1e6, float e=1.0)

	SimPhotonCounter (float t_p1, float t_p2, float t_l1, float t_l2, float min_w, float max_w, const std::vector< float > &eV)

void	SetVectorSize (size_t s)

size_t	GetVectorSize () const

void	SetWavelengthRanges (float min_w, float max_w)

float	MinWavelength () const

float	MaxWavelength () const

void	SetTimeRanges (float t_p1, float t_p2, float t_l1, float t_l2)

float	MinPromptTime () const

float	MaxPromptTime () const

float	MinLateTime () const

float	MaxLateTime () const

void	SetQE (size_t i, float e)

float	QE (size_t i) const

void	SetQEVector (const std::vector< float > &eV)

std::vector< float > const &	QEVector () const

void	AddOnePhoton (size_t i_opdet, const sim::OnePhoton &photon)

void	AddSimPhotons (const sim::SimPhotons &photons)

void	ClearVectors ()

const std::vector< float > &	PromptPhotonVector () const

const std::vector< float > &	LatePhotonVector () const

float	PromptPhotonVector (size_t i) const

float	LatePhotonVector (size_t i) const

std::vector< float >	TotalPhotonVector () const

float	TotalPhotonVector (size_t i) const

float	PromptPhotonTotal () const

float	LatePhotonTotal () const

float	PhotonTotal () const

void	Print ()

	SimPhotonCounter (const fhicl::ParameterSet &)

void	analyze (art::Event const &)

void	beginJob ()

void	endJob ()

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)

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)

std::unique_ptr< Worker >	makeWorker (WorkerParams const &wp)

template<typename T , BranchType BT>
ViewToken< T >	consumesView (InputTag const &tag)

template<typename T , BranchType BT>
ViewToken< T >	mayConsumeView (InputTag const &tag)

Protected Member Functions
std::string const &	processName () const

bool	wantAllEvents () const noexcept

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

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

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

Private Member Functions
float	Wavelength (const sim::OnePhoton &ph)

void	storeVisibility (int channel, int nDirectPhotons, int nReflectedPhotons, double reflectedT0=0.0) const

bool	isVisible (double wavelength) const
	Returns if we label as "visibile" a photon with specified wavelength [nm]. More...

Private Attributes
std::vector< float >	_photonVector_prompt

std::vector< float >	_photonVector_late

float	_min_prompt_time

float	_max_prompt_time

float	_min_late_time

float	_max_late_time

std::vector< float >	_qeVector

float	_min_wavelength

float	_max_wavelength

TTree *	fThePhotonTreeAll

TTree *	fThePhotonTreeDetected

TTree *	fTheOpDetTree

TTree *	fTheEventTree

std::vector< std::string >	fInputModule

int	fVerbosity

bool	fMakeDetectedPhotonsTree

bool	fMakeAllPhotonsTree

bool	fMakeOpDetsTree

bool	fMakeOpDetEventsTree

bool	fSavePhotonOrigins

TVector3	initialPhotonPosition

TVector3	finalPhotonPosition

Float_t	fWavelength

Float_t	fTime

Float_t	foriginX

Float_t	foriginY

Float_t	foriginZ

Int_t	fCountOpDetAll

Int_t	fCountOpDetDetected

Int_t	fCountOpDetReflDetected

Float_t	fT0_vis

Int_t	fCountEventAll

Int_t	fCountEventDetected

Int_t	fEventID

Int_t	fOpChannel

bool	fMakeLightAnalysisTree

std::vector< std::vector< std::vector< double > > >	fSignals_vuv

std::vector< std::vector< std::vector< double > > >	fSignals_vis

TTree *	fLightAnalysisTree = nullptr

int	fRun

int	fTrackID

int	fpdg

int	fmotherTrackID

double	fEnergy

double	fdEdx

std::vector< double >	fPosition0

std::vector< std::vector< double > >	fstepPositions

std::vector< double >	fstepTimes

std::vector< std::vector< double > >	fSignalsvuv

std::vector< std::vector< double > >	fSignalsvis

std::string	fProcess

cheat::ParticleInventoryService const *	pi_serv = nullptr

phot::PhotonVisibilityService const *	fPVS = nullptr

bool const	fUseLitePhotons

Static Private Attributes
static constexpr double	kVisibleThreshold = 200.0
	Threshold used to resolve between visible and ultraviolet light. More...

static constexpr double	kVisibleWavelength = 450.0
	Value used when a typical visible light wavelength is needed. More...

static constexpr double	kVUVWavelength = 128.0
	Value used when a typical ultraviolet light wavelength is needed. More...

Detailed Description

Definition at line 23 of file SimPhotonCounter.h.

Member Typedef Documentation

using art::EDAnalyzer::ModuleType = EDAnalyzer

inherited

Definition at line 22 of file EDAnalyzer.h.

Constructor & Destructor Documentation

opdet::SimPhotonCounter::SimPhotonCounter ( )

inline

Definition at line 26 of file SimPhotonCounter.h.

References e.

Referenced by isVisible().

26 {}

opdet::SimPhotonCounter::SimPhotonCounter	(	size_t	s,
		float	t_p1,
		float	t_p2,
		float	t_l1,
		float	t_l2,
		float	min_w = `0`,
		float	max_w = `1e6`,
		float	e = `1.0`
	)

Definition at line 10 of file SimPhotonCounter.cxx.

References _photonVector_late, _photonVector_prompt, _qeVector, e, SetTimeRanges(), and SetWavelengthRanges().

 {
 
   SetWavelengthRanges(min_w, max_w);
   SetTimeRanges(t_p1, t_p2, t_l1, t_l2);
 
   _photonVector_prompt = std::vector<float>(s);
   _photonVector_late = std::vector<float>(s);
   _qeVector = std::vector<float>(s, e);
 }

opdet::SimPhotonCounter::SimPhotonCounter	(	float	t_p1,
		float	t_p2,
		float	t_l1,
		float	t_l2,
		float	min_w,
		float	max_w,
		const std::vector< float > &	eV
	)

Definition at line 28 of file SimPhotonCounter.cxx.

References _photonVector_late, _photonVector_prompt, _qeVector, SetTimeRanges(), and SetWavelengthRanges().

 {
   SetWavelengthRanges(min_w, max_w);
   SetTimeRanges(t_p1, t_p2, t_l1, t_l2);
 
   _photonVector_prompt = std::vector<float>(eV.size());
   _photonVector_late = std::vector<float>(eV.size());
   _qeVector = eV;
 }

opdet::SimPhotonCounter::SimPhotonCounter ( const fhicl::ParameterSet & )

Member Function Documentation

void opdet::SimPhotonCounter::AddOnePhoton	(	size_t	i_opdet,
		const sim::OnePhoton &	photon
	)

Definition at line 71 of file SimPhotonCounter.cxx.

References _max_late_time, _max_prompt_time, _max_wavelength, _min_late_time, _min_prompt_time, _min_wavelength, _photonVector_late, _photonVector_prompt, _qeVector, GetVectorSize(), sim::OnePhoton::Time, and Wavelength().

Referenced by AddSimPhotons().

 {
   if (i_opdet > GetVectorSize())
     throw std::runtime_error("ERROR in SimPhotonCounter: Opdet requested out of range!");
 
   if (Wavelength(photon) < _min_wavelength || Wavelength(photon) > _max_wavelength) return;
 
   if (photon.Time > _min_prompt_time && photon.Time <= _max_prompt_time)
     _photonVector_prompt[i_opdet] += _qeVector[i_opdet];
   else if (photon.Time > _min_late_time && photon.Time < _max_late_time)
     _photonVector_late[i_opdet] += _qeVector[i_opdet];
 }

void opdet::SimPhotonCounter::AddSimPhotons ( const sim::SimPhotons & photons )

Definition at line 84 of file SimPhotonCounter.cxx.

References AddOnePhoton(), and sim::SimPhotons::OpChannel().

 {
   for (size_t i_ph = 0; i_ph < photons.size(); i_ph++)
     AddOnePhoton(photons.OpChannel(), photons[i_ph]);
 }

void opdet::SimPhotonCounter::analyze ( art::Event const & evt )

virtual

Implements art::EDAnalyzer.

Definition at line 298 of file SimPhotonCounter_module.cc.

References opdet::OpDetResponseInterface::detected(), opdet::OpDetResponseInterface::detectedLite(), sim::ParticleList::end(), art::EventID::event(), fCountEventAll, fCountEventDetected, fCountOpDetAll, fCountOpDetDetected, fCountOpDetReflDetected, fdEdx, fEnergy, fEventID, sim::ParticleList::find(), fInputModule, fLightAnalysisTree, fMakeAllPhotonsTree, fMakeDetectedPhotonsTree, fMakeLightAnalysisTree, fMakeOpDetEventsTree, fMakeOpDetsTree, fmotherTrackID, fOpChannel, foriginX, foriginY, foriginZ, fpdg, fProcess, fPVS, fRun, fSignals_vis, fSignals_vuv, fSignalsvis, fSignalsvuv, fstepPositions, fstepTimes, fT0_vis, fTheEventTree, fTheOpDetTree, fThePhotonTreeAll, fThePhotonTreeDetected, fTime, fTrackID, fUseLitePhotons, fVerbosity, fWavelength, Get, art::ProductRetriever::getHandle(), art::ProductRetriever::getMany(), art::ProductRetriever::getView(), art::Event::id(), phot::PhotonVisibilityService::IsBuildJob(), isVisible(), kVisibleThreshold, kVUVWavelength, sim::NoParticleId, cheat::ParticleInventoryService::ParticleList(), pi_serv, art::errors::ProductNotFound, art::Event::run(), sc, phot::PhotonVisibilityService::StoreReflected(), phot::PhotonVisibilityService::StoreReflT0(), storeVisibility(), and opdet::OpDetResponseInterface::wavelength().

   {
     // Lookup event ID from event
     art::EventNumber_t event = evt.id().event();
     fEventID = Int_t(event);
 
     // Service for determining opdet responses
     art::ServiceHandle<opdet::OpDetResponseInterface const> odresponse;
 
     // get the geometry to be able to figure out signal types and chan -> plane mappings
     auto const& wireReadoutGeom = art::ServiceHandle<geo::WireReadout const>()->Get();
 
     // GEANT4 info on the particles (only used if making light analysis tree)
     std::vector<simb::MCParticle> const* mcpartVec = nullptr;
 
     //-------------------------initializing light tree vectors------------------------
     std::vector<double> totalEnergy_track;
     fstepPositions.clear();
     fstepTimes.clear();
     if (fMakeLightAnalysisTree) {
       mcpartVec = evt.getHandle<std::vector<simb::MCParticle>>("largeant").product();
 
       size_t maxNtracks = 1000U; // mcpartVec->size(); --- { to be fixed soon! ]
       fSignals_vuv.clear();
       fSignals_vuv.resize(maxNtracks);
       fSignals_vis.clear();
       fSignals_vis.resize(maxNtracks);
       for (size_t itrack = 0; itrack != maxNtracks; itrack++) {
         fSignals_vuv[itrack].resize(wireReadoutGeom.NOpChannels());
         fSignals_vis[itrack].resize(wireReadoutGeom.NOpChannels());
       }
       totalEnergy_track.resize(maxNtracks, 0.);
       //-------------------------stimation of dedx per trackID----------------------
       //get the list of particles from this event
       const sim::ParticleList* plist = pi_serv ? &(pi_serv->ParticleList()) : nullptr;
 
       // loop over all sim::SimChannels in the event and make sure there are no
       // sim::IDEs with trackID values that are not in the sim::ParticleList
       std::vector<const sim::SimChannel*> sccol;
       for (auto const& mod : fInputModule) {
         evt.getView(mod, sccol);
         //loop over the sim channels collection
         for (size_t sc = 0; sc < sccol.size(); ++sc) {
           const auto& tdcidemap = sccol[sc]->TDCIDEMap();
           //loop over all of the tdc IDE map objects
           for (auto mapitr = tdcidemap.begin(); mapitr != tdcidemap.end(); mapitr++) {
             const std::vector<sim::IDE> idevec = (*mapitr).second;
             //go over all of the IDEs in a given simchannel
             for (size_t iv = 0; iv < idevec.size(); ++iv) {
               if (plist) {
                 if (plist->find(idevec[iv].trackID) == plist->end() &&
                     idevec[iv].trackID != sim::NoParticleId) {
                   mf::LogWarning("LArG4Ana") << idevec[iv].trackID << " is not in particle list";
                 }
               }
               if (idevec[iv].trackID < 0) continue;
 
               totalEnergy_track[idevec[iv].trackID] += idevec[iv].energy / 3.;
             }
           }
         }
       }
     } //End of if(fMakeLightAnalysisTree)
 
     if (!fUseLitePhotons) {
 
       //Reset counters
       fCountEventAll = 0;
       fCountEventDetected = 0;
 
       //Get *ALL* SimPhotonsCollection from Event
       auto photon_handles = evt.getMany<std::vector<sim::SimPhotons>>();
       if (photon_handles.size() == 0)
         throw art::Exception(art::errors::ProductNotFound)
           << "sim SimPhotons retrieved and you requested them.";
 
       for (auto const& mod : fInputModule) {
         //switching off to add reading in of labelled collections: Andrzej, 02/26/19
 
         for (auto const& ph_handle : photon_handles) {
           // Do some checking before we proceed
           if (!ph_handle.isValid()) continue;
           if (ph_handle.provenance()->moduleLabel() != mod)
             continue; //not the most efficient way of doing this, but preserves the logic of the module. Andrzej
 
           bool Reflected = (ph_handle.provenance()->productInstanceName() == "Reflected");
 
           if ((*ph_handle).size() > 0) {
             if (fMakeLightAnalysisTree) {
               //resetting the signalt to save in the analysis tree per event
               const int maxNtracks = 1000;
               for (size_t itrack = 0; itrack != maxNtracks; itrack++) {
                 for (size_t pmt_i = 0; pmt_i != wireReadoutGeom.NOpChannels(); pmt_i++) {
                   fSignals_vuv[itrack][pmt_i].clear();
                   fSignals_vis[itrack][pmt_i].clear();
                 }
               }
             }
           }
 
           if (fVerbosity > 0)
             std::cout << "Found OpDet hit collection of size " << (*ph_handle).size() << std::endl;
 
           if ((*ph_handle).size() > 0) {
             for (auto const& itOpDet : (*ph_handle)) {
               //Reset Counters
               fCountOpDetAll = 0;
               fCountOpDetDetected = 0;
               fCountOpDetReflDetected = 0;
               //Reset t0 for visible light
               fT0_vis = 999.;
 
               //Get data from HitCollection entry
               fOpChannel = itOpDet.OpChannel();
               const sim::SimPhotons& TheHit = itOpDet;
 
               // Loop through OpDet phots.
               //   Note we make the screen output decision outside the loop
               //   in order to avoid evaluating large numbers of unnecessary
               //   if conditions.
 
               for (const sim::OnePhoton& Phot : TheHit) {
                 // Calculate wavelength in nm
                 fWavelength = odresponse->wavelength(Phot.Energy);
 
                 //Get arrival time from phot
                 fTime = Phot.Time;
                 foriginX = Phot.InitialPosition.X();
                 foriginY = Phot.InitialPosition.Y();
                 foriginZ = Phot.InitialPosition.Z();
 
                 // special case for LibraryBuildJob: no working "Reflected" handle and all photons stored in single object - must sort using wavelength instead
                 if (fPVS->IsBuildJob() && !Reflected) {
                   // all photons contained in object with Reflected = false flag
                   // Increment per OpDet counters and fill per phot trees
                   fCountOpDetAll++;
                   if (fMakeAllPhotonsTree) {
                     if (!isVisible(fWavelength) || fPVS->StoreReflected()) {
                       fThePhotonTreeAll->Fill();
                     }
                   }
 
                   if (odresponse->detected(fOpChannel, Phot)) {
                     if (fMakeDetectedPhotonsTree) fThePhotonTreeDetected->Fill();
                     //only store direct direct light
                     if (!isVisible(fWavelength)) fCountOpDetDetected++;
                     // reflected and shifted light is in visible range
                     else if (fPVS->StoreReflected()) {
                       fCountOpDetReflDetected++;
                       // find the first visible arrival time
                       if (fPVS->StoreReflT0() && fTime < fT0_vis) fT0_vis = fTime;
                     }
                     if (fVerbosity > 3)
                       std::cout << "OpDetResponseInterface PerPhoton : Event " << fEventID
                                 << " OpChannel " << fOpChannel << " Wavelength " << fWavelength
                                 << " Detected 1 " << std::endl;
                   }
                   else {
                     if (fVerbosity > 3)
                       std::cout << "OpDetResponseInterface PerPhoton : Event " << fEventID
                                 << " OpChannel " << fOpChannel << " Wavelength " << fWavelength
                                 << " Detected 0 " << std::endl;
                   }
 
                 } // if build library and not reflected
 
                 else {
                   // store in appropriate trees using "Reflected" handle and fPVS->StoreReflected() flag
                   // Increment per OpDet counters and fill per phot trees
                   fCountOpDetAll++;
                   if (fMakeAllPhotonsTree) {
                     if (!Reflected || (fPVS->StoreReflected() && Reflected)) {
                       fThePhotonTreeAll->Fill();
                     }
                   }
 
                   if (odresponse->detected(fOpChannel, Phot)) {
                     if (fMakeDetectedPhotonsTree) fThePhotonTreeDetected->Fill();
                     //only store direct direct light
                     if (!Reflected) fCountOpDetDetected++;
                     // reflected and shifted light is in visible range
                     else if (fPVS->StoreReflected() && Reflected) {
                       fCountOpDetReflDetected++;
                       // find the first visible arrival time
                       if (fPVS->StoreReflT0() && fTime < fT0_vis) fT0_vis = fTime;
                     }
                     if (fVerbosity > 3)
                       std::cout << "OpDetResponseInterface PerPhoton : Event " << fEventID
                                 << " OpChannel " << fOpChannel << " Wavelength " << fWavelength
                                 << " Detected 1 " << std::endl;
                   }
                   else {
                     if (fVerbosity > 3)
                       std::cout << "OpDetResponseInterface PerPhoton : Event " << fEventID
                                 << " OpChannel " << fOpChannel << " Wavelength " << fWavelength
                                 << " Detected 0 " << std::endl;
                   }
                 }
               } // for each photon in collection
 
               // If this is a library building job, fill relevant entry
               if (
                 fPVS->IsBuildJob() &&
                 !Reflected) // for library build job, both componenents stored in first object with Reflected = false
               {
                 storeVisibility(fOpChannel, fCountOpDetDetected, fCountOpDetReflDetected, fT0_vis);
               }
 
               // Incremenent per event and fill Per OpDet trees
               if (fMakeOpDetsTree) fTheOpDetTree->Fill();
               fCountEventAll += fCountOpDetAll;
               fCountEventDetected += fCountOpDetDetected;
 
               // Give per OpDet output
               if (fVerbosity > 2)
                 std::cout << "OpDetResponseInterface PerOpDet : Event " << fEventID << " OpDet "
                           << fOpChannel << " All " << fCountOpDetAll << " Det "
                           << fCountOpDetDetected << std::endl;
             }
 
             // Fill per event tree
             if (fMakeOpDetEventsTree) fTheEventTree->Fill();
 
             // Give per event output
             if (fVerbosity > 1)
               std::cout << "OpDetResponseInterface PerEvent : Event " << fEventID << " All "
                         << fCountOpDetAll << " Det " << fCountOpDetDetected << std::endl;
           }
           else {
             // if empty OpDet hit collection,
             // add an empty record to the per event tree
             if (fMakeOpDetEventsTree) fTheEventTree->Fill();
           }
           if (fMakeLightAnalysisTree) {
             assert(mcpartVec);
             assert(fLightAnalysisTree);
 
             std::cout << "Filling the analysis tree" << std::endl;
             //---------------Filling the analysis tree-----------:
             fRun = evt.run();
             std::vector<double> this_xyz;
 
             //loop over the particles
             for (simb::MCParticle const& pPart : *mcpartVec) {
 
               if (pPart.Process() == "primary") fEnergy = pPart.E();
 
               //resetting the vectors
               fstepPositions.clear();
               fstepTimes.clear();
               fSignalsvuv.clear();
               fSignalsvis.clear();
               fdEdx = -1.;
               //filling the tree fields
               fTrackID = pPart.TrackId();
               fpdg = pPart.PdgCode();
               fmotherTrackID = pPart.Mother();
               fdEdx = totalEnergy_track[fTrackID];
               fSignalsvuv = fSignals_vuv[fTrackID];
               fSignalsvis = fSignals_vis[fTrackID];
               fProcess = pPart.Process();
               //filling the center positions of each step
               for (size_t i_s = 1; i_s < pPart.NumberTrajectoryPoints(); i_s++) {
                 this_xyz.clear();
                 this_xyz.resize(3);
                 this_xyz[0] = pPart.Position(i_s).X();
                 this_xyz[1] = pPart.Position(i_s).Y();
                 this_xyz[2] = pPart.Position(i_s).Z();
                 fstepPositions.push_back(this_xyz);
                 fstepTimes.push_back(pPart.Position(i_s).T());
               }
               //filling the tree per track
               fLightAnalysisTree->Fill();
             }
           } // if fMakeLightAnalysisTree
         }
       }
     }
     if (fUseLitePhotons) {
 
       //Get *ALL* SimPhotonsCollection from Event
       auto photon_handles = evt.getMany<std::vector<sim::SimPhotonsLite>>();
       if (photon_handles.size() == 0)
         throw art::Exception(art::errors::ProductNotFound)
           << "sim SimPhotons retrieved and you requested them.";
 
       //Get SimPhotonsLite from Event
       for (auto const& mod : fInputModule) {
         // Loop over direct/reflected photons
         for (auto const& ph_handle : photon_handles) {
           // Do some checking before we proceed
           if (!ph_handle.isValid()) continue;
           if (ph_handle.provenance()->moduleLabel() != mod)
             continue; //not the most efficient way of doing this, but preserves the logic of the module. Andrzej
 
           bool Reflected = (ph_handle.provenance()->productInstanceName() == "Reflected");
 
           //Reset counters
           fCountEventAll = 0;
           fCountEventDetected = 0;
 
           if (fVerbosity > 0)
             std::cout << "Found OpDet hit collection of size " << (*ph_handle).size() << std::endl;
 
           if ((*ph_handle).size() > 0) {
 
             for (auto const& photon : (*ph_handle)) {
               //Get data from HitCollection entry
               fOpChannel = photon.OpChannel;
               std::map<int, int> PhotonsMap = photon.DetectedPhotons;
 
               //Reset Counters
               fCountOpDetAll = 0;
               fCountOpDetDetected = 0;
               fCountOpDetReflDetected = 0;
 
               for (auto it = PhotonsMap.begin(); it != PhotonsMap.end(); it++) {
                 // Calculate wavelength in nm
                 if (Reflected) { fWavelength = kVisibleThreshold; }
                 else {
                   fWavelength = kVUVWavelength; // original
                 }
 
                 //Get arrival time from phot
                 fTime = it->first;
 
                 for (int i = 0; i < it->second; i++) {
                   // Increment per OpDet counters and fill per phot trees
                   fCountOpDetAll++;
                   if (fMakeAllPhotonsTree) fThePhotonTreeAll->Fill();
 
                   if (odresponse->detectedLite(fOpChannel)) {
                     if (fMakeDetectedPhotonsTree) fThePhotonTreeDetected->Fill();
                     // direct light
                     if (!Reflected) { fCountOpDetDetected++; }
                     else if (Reflected) {
                       fCountOpDetReflDetected++;
                     }
                     if (fVerbosity > 3)
                       std::cout << "OpDetResponseInterface PerPhoton : Event " << fEventID
                                 << " OpChannel " << fOpChannel << " Wavelength " << fWavelength
                                 << " Detected 1 " << std::endl;
                   }
                   else if (fVerbosity > 3) {
                     std::cout << "OpDetResponseInterface PerPhoton : Event " << fEventID
                               << " OpChannel " << fOpChannel << " Wavelength " << fWavelength
                               << " Detected 0 " << std::endl;
                   }
                 }
               }
 
               // Incremenent per event and fill Per OpDet trees
               if (fMakeOpDetsTree) fTheOpDetTree->Fill();
               fCountEventAll += fCountOpDetAll;
               fCountEventDetected += fCountOpDetDetected;
 
               if (fPVS->IsBuildJob())
                 storeVisibility(fOpChannel, fCountOpDetDetected, fCountOpDetReflDetected, fT0_vis);
 
               // Give per OpDet output
               if (fVerbosity > 2)
                 std::cout << "OpDetResponseInterface PerOpDet : Event " << fEventID << " OpDet "
                           << fOpChannel << " All " << fCountOpDetAll << " Det "
                           << fCountOpDetDetected << std::endl;
             }
             // Fill per event tree
             if (fMakeOpDetEventsTree) fTheEventTree->Fill();
 
             // Give per event output
             if (fVerbosity > 1)
               std::cout << "OpDetResponseInterface PerEvent : Event " << fEventID << " All "
                         << fCountOpDetAll << " Det " << fCountOpDetDetected << std::endl;
           }
           else {
             // if empty OpDet hit collection,
             // add an empty record to the per event tree
             if (fMakeOpDetEventsTree) fTheEventTree->Fill();
           }
         }
       }
     }
   } // SimPhotonCounter::analyze()

void opdet::SimPhotonCounter::beginJob ( )

virtual

Reimplemented from art::EDAnalyzer.

Definition at line 198 of file SimPhotonCounter_module.cc.

References geo::CryostatGeo::Boundaries(), geo::GeometryCore::Cryostat(), e, fCountEventAll, fCountEventDetected, fCountOpDetAll, fCountOpDetDetected, fCountOpDetReflDetected, fdEdx, fEnergy, fEventID, fLightAnalysisTree, fMakeAllPhotonsTree, fMakeDetectedPhotonsTree, fMakeLightAnalysisTree, fMakeOpDetEventsTree, fMakeOpDetsTree, fmotherTrackID, fOpChannel, foriginX, foriginY, foriginZ, fpdg, fProcess, fPVS, fRun, fSavePhotonOrigins, fSignalsvis, fSignalsvuv, fstepPositions, fstepTimes, fTheEventTree, fTheOpDetTree, fThePhotonTreeAll, fThePhotonTreeDetected, fTime, fTrackID, fWavelength, geo::OpDetGeo::GetCenter(), geo::BoxBoundedGeo::MinX(), geo::CryostatGeo::NOpDet(), geo::GeometryCore::OpDetGeoFromOpDet(), pi_serv, art::errors::ServiceNotFound, and phot::PhotonVisibilityService::StoreReflected().

   {
     // Get file service to store trees
     art::ServiceHandle<art::TFileService const> tfs;
     art::ServiceHandle<geo::Geometry const> geo;
 
     std::cout << "Optical Channels positions:  " << geo->Cryostat().NOpDet() << std::endl;
     for (int ch = 0; ch != int(geo->Cryostat().NOpDet()); ch++) {
       auto const OpDetCenter = geo->OpDetGeoFromOpDet(ch).GetCenter();
       std::cout << ch << "  " << OpDetCenter.X() << "  " << OpDetCenter.Y() << "  "
                 << OpDetCenter.Z() << std::endl;
     }
 
     auto const CryoBounds = geo->Cryostat().Boundaries();
     std::cout << "Cryo Boundaries" << std::endl;
     std::cout << "Xmin: " << CryoBounds.MinX() << " Xmax: " << CryoBounds.MaxX()
               << " Ymin: " << CryoBounds.MinY() << " Ymax: " << CryoBounds.MaxY()
               << " Zmin: " << CryoBounds.MinZ() << " Zmax: " << CryoBounds.MaxZ() << std::endl;
 
     try {
       pi_serv = &*(art::ServiceHandle<cheat::ParticleInventoryService const>());
     }
     catch (art::Exception const& e) {
       if (e.categoryCode() != art::errors::ServiceNotFound) throw;
       mf::LogError("SimPhotonCounter")
         << "ParticleInventoryService service is not configured!"
            " Please add it in the job configuration."
            " In the meanwhile, some checks to particles will be skipped.";
     }
 
     // Create and assign branch addresses to required tree
     if (fMakeAllPhotonsTree) {
       fThePhotonTreeAll = tfs->make<TTree>("AllPhotons", "AllPhotons");
       fThePhotonTreeAll->Branch("EventID", &fEventID, "EventID/I");
       fThePhotonTreeAll->Branch("Wavelength", &fWavelength, "Wavelength/F");
       fThePhotonTreeAll->Branch("OpChannel", &fOpChannel, "OpChannel/I");
       fThePhotonTreeAll->Branch("Time", &fTime, "Time/F");
       if (fSavePhotonOrigins) {
         fThePhotonTreeAll->Branch("originX", &foriginX, "originX/F");
         fThePhotonTreeAll->Branch("originY", &foriginY, "originY/F");
         fThePhotonTreeAll->Branch("originZ", &foriginZ, "originZ/F");
       }
     }
 
     if (fMakeDetectedPhotonsTree) {
       fThePhotonTreeDetected = tfs->make<TTree>("DetectedPhotons", "DetectedPhotons");
       fThePhotonTreeDetected->Branch("EventID", &fEventID, "EventID/I");
       fThePhotonTreeDetected->Branch("Wavelength", &fWavelength, "Wavelength/F");
       fThePhotonTreeDetected->Branch("OpChannel", &fOpChannel, "OpChannel/I");
       fThePhotonTreeDetected->Branch("Time", &fTime, "Time/F");
       if (fSavePhotonOrigins) {
         fThePhotonTreeDetected->Branch("originX", &foriginX, "originX/F");
         fThePhotonTreeDetected->Branch("originY", &foriginY, "originY/F");
         fThePhotonTreeDetected->Branch("originZ", &foriginZ, "originZ/F");
       }
     }
 
     if (fMakeOpDetsTree) {
       fTheOpDetTree = tfs->make<TTree>("OpDets", "OpDets");
       fTheOpDetTree->Branch("EventID", &fEventID, "EventID/I");
       fTheOpDetTree->Branch("OpChannel", &fOpChannel, "OpChannel/I");
       fTheOpDetTree->Branch("CountAll", &fCountOpDetAll, "CountAll/I");
       fTheOpDetTree->Branch("CountDetected", &fCountOpDetDetected, "CountDetected/I");
       if (fPVS->StoreReflected())
         fTheOpDetTree->Branch("CountReflDetected", &fCountOpDetReflDetected, "CountReflDetected/I");
       fTheOpDetTree->Branch("Time", &fTime, "Time/F");
     }
 
     if (fMakeOpDetEventsTree) {
       fTheEventTree = tfs->make<TTree>("OpDetEvents", "OpDetEvents");
       fTheEventTree->Branch("EventID", &fEventID, "EventID/I");
       fTheEventTree->Branch("CountAll", &fCountEventAll, "CountAll/I");
       fTheEventTree->Branch("CountDetected", &fCountEventDetected, "CountDetected/I");
       if (fPVS->StoreReflected())
         fTheOpDetTree->Branch("CountReflDetected", &fCountOpDetReflDetected, "CountReflDetected/I");
     }
 
     //generating the tree for the light analysis:
     if (fMakeLightAnalysisTree) {
       fLightAnalysisTree = tfs->make<TTree>("LightAnalysis", "LightAnalysis");
       fLightAnalysisTree->Branch("RunNumber", &fRun);
       fLightAnalysisTree->Branch("EventID", &fEventID);
       fLightAnalysisTree->Branch("TrackID", &fTrackID);
       fLightAnalysisTree->Branch("PdgCode", &fpdg);
       fLightAnalysisTree->Branch("MotherTrackID", &fmotherTrackID);
       fLightAnalysisTree->Branch("Energy", &fEnergy);
       fLightAnalysisTree->Branch("dEdx", &fdEdx);
       fLightAnalysisTree->Branch("StepPositions", &fstepPositions);
       fLightAnalysisTree->Branch("StepTimes", &fstepTimes);
       fLightAnalysisTree->Branch("SignalsVUV", &fSignalsvuv);
       fLightAnalysisTree->Branch("SignalsVisible", &fSignalsvis);
       fLightAnalysisTree->Branch("Process", &fProcess);
     }
   }

void opdet::SimPhotonCounter::ClearVectors ( )

Definition at line 90 of file SimPhotonCounter.cxx.

References _photonVector_late, _photonVector_prompt, and GetVectorSize().

 {
   for (size_t i = 0; i < GetVectorSize(); i++) {
     _photonVector_prompt[i] = 0.0;
     _photonVector_late[i] = 0.0;
   }
 }

template<typename T , BranchType BT>

ProductToken< T > art::ModuleBase::consumes ( InputTag const & tag )

protectedinherited

Definition at line 61 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::consumes().

   {
     return collector_.consumes<T, BT>(tag);
   }

ConsumesCollector & art::ModuleBase::consumesCollector ( )

protectedinherited

Definition at line 57 of file ModuleBase.cc.

References art::ModuleBase::collector_.

   {
     return collector_;
   }

template<typename T , BranchType BT>

void art::ModuleBase::consumesMany ( )

protectedinherited

Definition at line 75 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::consumesMany().

   {
     collector_.consumesMany<T, BT>();
   }

template<typename Element , BranchType = InEvent>

ViewToken<Element> art::ModuleBase::consumesView ( InputTag const & )

protectedinherited

template<typename T , BranchType BT>

ViewToken<T> art::ModuleBase::consumesView ( InputTag const & tag )

inherited

Definition at line 68 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::consumesView().

   {
     return collector_.consumesView<T, BT>(tag);
   }

void art::detail::Analyzer::doBeginJob ( SharedResources const & resources )

inherited

Definition at line 25 of file Analyzer.cc.

Referenced by art::detail::Analyzer::Analyzer().

   {
     setupQueues(resources);
     ProcessingFrame const frame{ScheduleID{}};
     beginJobWithFrame(frame);
   }

bool art::detail::Analyzer::doBeginRun	(	RunPrincipal &	rp,
		ModuleContext const &	mc
	)

inherited

Definition at line 68 of file Analyzer.cc.

References art::ModuleContext::scheduleID().

Referenced by art::detail::Analyzer::Analyzer().

   {
     ProcessingFrame const frame{mc.scheduleID()};
     beginRunWithFrame(std::as_const(rp).makeRun(mc), frame);
     return true;
   }

bool art::detail::Analyzer::doBeginSubRun	(	SubRunPrincipal &	srp,
		ModuleContext const &	mc
	)

inherited

Definition at line 84 of file Analyzer.cc.

References art::ModuleContext::scheduleID().

Referenced by art::detail::Analyzer::Analyzer().

   {
     ProcessingFrame const frame{mc.scheduleID()};
     beginSubRunWithFrame(std::as_const(srp).makeSubRun(mc), frame);
     return true;
   }

void art::detail::Analyzer::doEndJob ( )

inherited

Definition at line 33 of file Analyzer.cc.

Referenced by art::detail::Analyzer::Analyzer().

   {
     ProcessingFrame const frame{ScheduleID{}};
     endJobWithFrame(frame);
   }

bool art::detail::Analyzer::doEndRun	(	RunPrincipal &	rp,
		ModuleContext const &	mc
	)

inherited

Definition at line 76 of file Analyzer.cc.

References art::ModuleContext::scheduleID().

Referenced by art::detail::Analyzer::Analyzer().

   {
     ProcessingFrame const frame{mc.scheduleID()};
     endRunWithFrame(std::as_const(rp).makeRun(mc), frame);
     return true;
   }

bool art::detail::Analyzer::doEndSubRun	(	SubRunPrincipal &	srp,
		ModuleContext const &	mc
	)

inherited

Definition at line 92 of file Analyzer.cc.

References art::ModuleContext::scheduleID().

Referenced by art::detail::Analyzer::Analyzer().

   {
     ProcessingFrame const frame{mc.scheduleID()};
     endSubRunWithFrame(std::as_const(srp).makeSubRun(mc), frame);
     return true;
   }

bool art::detail::Analyzer::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
	)

inherited

Definition at line 100 of file Analyzer.cc.

References e, and art::ModuleContext::scheduleID().

Referenced by art::detail::Analyzer::Analyzer().

   {
     auto const e = std::as_const(ep).makeEvent(mc);
     if (wantEvent(mc.scheduleID(), e)) {
       ++counts_run;
       ProcessingFrame const frame{mc.scheduleID()};
       analyzeWithFrame(e, frame);
       ++counts_passed;
     }
     return true;
   }

void art::detail::Analyzer::doRespondToCloseInputFile ( FileBlock const & fb )

inherited

Definition at line 47 of file Analyzer.cc.

Referenced by art::detail::Analyzer::Analyzer().

   {
     ProcessingFrame const frame{ScheduleID{}};
     respondToCloseInputFileWithFrame(fb, frame);
   }

void art::detail::Analyzer::doRespondToCloseOutputFiles ( FileBlock const & fb )

inherited

Definition at line 61 of file Analyzer.cc.

Referenced by art::detail::Analyzer::Analyzer().

   {
     ProcessingFrame const frame{ScheduleID{}};
     respondToCloseOutputFilesWithFrame(fb, frame);
   }

void art::detail::Analyzer::doRespondToOpenInputFile ( FileBlock const & fb )

inherited

Definition at line 40 of file Analyzer.cc.

Referenced by art::detail::Analyzer::Analyzer().

   {
     ProcessingFrame const frame{ScheduleID{}};
     respondToOpenInputFileWithFrame(fb, frame);
   }

void art::detail::Analyzer::doRespondToOpenOutputFiles ( FileBlock const & fb )

inherited

Definition at line 54 of file Analyzer.cc.

Referenced by art::detail::Analyzer::Analyzer().

   {
     ProcessingFrame const frame{ScheduleID{}};
     respondToOpenOutputFilesWithFrame(fb, frame);
   }

void opdet::SimPhotonCounter::endJob ( )

virtual

Reimplemented from art::EDAnalyzer.

Definition at line 293 of file SimPhotonCounter_module.cc.

References fPVS, and phot::PhotonVisibilityService::IsBuildJob().

   {
     if (fPVS->IsBuildJob()) { art::ServiceHandle<phot::PhotonVisibilityService>()->StoreLibrary(); }
   }

std::array< std::vector< ProductInfo >, NumBranchTypes > const & art::ModuleBase::getConsumables ( ) const

inherited

Definition at line 43 of file ModuleBase.cc.

References art::ModuleBase::collector_, and art::ConsumesCollector::getConsumables().

   {
     return collector_.getConsumables();
   }

Handle< TriggerResults > art::Observer::getTriggerResults ( Event const & e ) const

protectedinherited

Definition at line 75 of file Observer.cc.

References art::ProductRetriever::get(), and art::Observer::selectors_.

Referenced by art::OutputModule::doWriteEvent(), and art::Observer::wantAllEvents().

   {
     if (selectors_) {
       return selectors_->getOneTriggerResults(e);
     }
 
     // The following applies for cases where no SelectEvents entries
     // exist.
     Handle<TriggerResults> h;
     if (e.get(empty_process_name, h)) {
       return h;
     }
     return Handle<TriggerResults>{};
   }

size_t opdet::SimPhotonCounter::GetVectorSize ( ) const

inline

Definition at line 50 of file SimPhotonCounter.h.

Referenced by AddOnePhoton(), ClearVectors(), opdet::FlashHypothesisComparison::FillSimPhotonCounterInfo(), Print(), and TotalPhotonVector().

50 { return _photonVector_prompt.size(); }

opdet::SimPhotonCounter::_photonVector_prompt

std::vector< float > _photonVector_prompt

Definition: SimPhotonCounter.h:102

bool opdet::SimPhotonCounter::isVisible ( double wavelength ) const

inlineprivate

Returns if we label as "visibile" a photon with specified wavelength [nm].

Definition at line 163 of file SimPhotonCounter_module.cc.

References art::errors::Configuration, art::EDAnalyzer::EDAnalyzer(), fInputModule, fMakeAllPhotonsTree, fMakeDetectedPhotonsTree, fMakeLightAnalysisTree, fMakeOpDetEventsTree, fMakeOpDetsTree, fPVS, fSavePhotonOrigins, fUseLitePhotons, fVerbosity, fhicl::ParameterSet::get(), phot::PhotonVisibilityService::IsBuildJob(), kVisibleThreshold, SimPhotonCounter(), phot::PhotonVisibilityService::StoreReflected(), and phot::PhotonVisibilityService::StoreReflT0().

Referenced by analyze().

163 { return fWavelength < kVisibleThreshold; }

opdet::SimPhotonCounter::fWavelength

Float_t fWavelength

Definition: SimPhotonCounter_module.cc:119

opdet::SimPhotonCounter::kVisibleThreshold

static constexpr double kVisibleThreshold

Threshold used to resolve between visible and ultraviolet light.

Definition: SimPhotonCounter_module.cc:82

float opdet::SimPhotonCounter::LatePhotonTotal ( ) const

inline

Definition at line 93 of file SimPhotonCounter.h.

Referenced by opdet::FlashHypothesisComparison::FillSimPhotonCounterInfo().

     {
       return std::accumulate(_photonVector_late.begin(), _photonVector_late.end(), 0.0);
     }

const std::vector<float>& opdet::SimPhotonCounter::LatePhotonVector ( ) const

inline

Definition at line 79 of file SimPhotonCounter.h.

Referenced by opdet::FlashHypothesisComparison::FillComparisonInfo(), opdet::FlashHypothesisComparison::FillSimPhotonCounterInfo(), and TotalPhotonVector().

79 { return _photonVector_late; }

opdet::SimPhotonCounter::_photonVector_late

std::vector< float > _photonVector_late

Definition: SimPhotonCounter.h:103

float opdet::SimPhotonCounter::LatePhotonVector ( size_t i ) const

inline

Definition at line 81 of file SimPhotonCounter.h.

81 { return _photonVector_late.at(i); }

opdet::SimPhotonCounter::_photonVector_late

std::vector< float > _photonVector_late

Definition: SimPhotonCounter.h:103

std::unique_ptr< Worker > art::ModuleBase::makeWorker ( WorkerParams const & wp )

inherited

Definition at line 37 of file ModuleBase.cc.

References art::ModuleBase::doMakeWorker(), and art::NumBranchTypes.

   {
     return doMakeWorker(wp);
   }

float opdet::SimPhotonCounter::MaxLateTime ( ) const

inline

Definition at line 62 of file SimPhotonCounter.h.

Referenced by Print().

62 { return _max_late_time; }

opdet::SimPhotonCounter::_max_late_time

float _max_late_time

Definition: SimPhotonCounter.h:108

float opdet::SimPhotonCounter::MaxPromptTime ( ) const

inline

Definition at line 60 of file SimPhotonCounter.h.

Referenced by Print().

60 { return _max_prompt_time; }

opdet::SimPhotonCounter::_max_prompt_time

float _max_prompt_time

Definition: SimPhotonCounter.h:106

float opdet::SimPhotonCounter::MaxWavelength ( ) const

inline

Definition at line 55 of file SimPhotonCounter.h.

55 { return _max_wavelength; }

opdet::SimPhotonCounter::_max_wavelength

float _max_wavelength

Definition: SimPhotonCounter.h:112

template<typename T , BranchType BT>

ProductToken< T > art::ModuleBase::mayConsume ( InputTag const & tag )

protectedinherited

Definition at line 82 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::mayConsume().

   {
     return collector_.mayConsume<T, BT>(tag);
   }

template<typename T , BranchType BT>

void art::ModuleBase::mayConsumeMany ( )

protectedinherited

Definition at line 96 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::mayConsumeMany().

   {
     collector_.mayConsumeMany<T, BT>();
   }

template<typename Element , BranchType = InEvent>

ViewToken<Element> art::ModuleBase::mayConsumeView ( InputTag const & )

protectedinherited

template<typename T , BranchType BT>

ViewToken<T> art::ModuleBase::mayConsumeView ( InputTag const & tag )

inherited

Definition at line 89 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::mayConsumeView().

   {
     return collector_.mayConsumeView<T, BT>(tag);
   }

float opdet::SimPhotonCounter::MinLateTime ( ) const

inline

Definition at line 61 of file SimPhotonCounter.h.

Referenced by Print().

61 { return _min_late_time; }

opdet::SimPhotonCounter::_min_late_time

float _min_late_time

Definition: SimPhotonCounter.h:107

float opdet::SimPhotonCounter::MinPromptTime ( ) const

inline

Definition at line 59 of file SimPhotonCounter.h.

Referenced by Print().

59 { return _min_prompt_time; }

opdet::SimPhotonCounter::_min_prompt_time

float _min_prompt_time

Definition: SimPhotonCounter.h:105

float opdet::SimPhotonCounter::MinWavelength ( ) const

inline

Definition at line 54 of file SimPhotonCounter.h.

54 { return _min_wavelength; }

opdet::SimPhotonCounter::_min_wavelength

float _min_wavelength

Definition: SimPhotonCounter.h:111

ModuleDescription const & art::ModuleBase::moduleDescription ( ) const

inherited

Definition at line 13 of file ModuleBase.cc.

References art::errors::LogicError.

Referenced by art::OutputModule::doRespondToOpenInputFile(), art::OutputModule::doWriteEvent(), art::Modifier::fillProductDescriptions(), art::OutputModule::makePlugins_(), art::OutputWorker::OutputWorker(), reco::shower::LArPandoraModularShowerCreation::produce(), art::Modifier::registerProducts(), and art::OutputModule::registerProducts().

   {
     if (md_.has_value()) {
       return *md_;
     }
 
     throw Exception{errors::LogicError,
                     "There was an error while calling moduleDescription().\n"}
       << "The moduleDescription() base-class member function cannot be called\n"
          "during module construction.  To determine which module is "
          "responsible\n"
          "for calling it, find the '<module type>:<module "
          "label>@Construction'\n"
          "tag in the message prefix above.  Please contact artists@fnal.gov\n"
          "for guidance.\n";
   }

float opdet::SimPhotonCounter::PhotonTotal ( ) const

inline

Definition at line 97 of file SimPhotonCounter.h.

References pyG4Element::Print().

97 { return (PromptPhotonTotal() + LatePhotonTotal()); }

opdet::SimPhotonCounter::PromptPhotonTotal

float PromptPhotonTotal() const

Definition: SimPhotonCounter.h:89

opdet::SimPhotonCounter::LatePhotonTotal

float LatePhotonTotal() const

Definition: SimPhotonCounter.h:93

void opdet::SimPhotonCounter::Print ( )

Definition at line 110 of file SimPhotonCounter.cxx.

References _photonVector_late, _photonVector_prompt, _qeVector, GetVectorSize(), MaxLateTime(), MaxPromptTime(), MinLateTime(), MinPromptTime(), and TotalPhotonVector().

 {
   std::cout << "Vector size: " << GetVectorSize() << std::endl;
   std::cout << "Time cut ranges: (" << MinPromptTime() << "," << MaxPromptTime() << ") , ("
             << MinLateTime() << "," << MaxLateTime() << ")" << std::endl;
   std::cout << "\t"
             << "i : QE / Prompt / Late / Total" << std::endl;
   for (size_t i = 0; i < GetVectorSize(); i++)
     std::cout << "\t" << i << ": " << _qeVector[i] << " / " << _photonVector_prompt[i] << " / "
               << _photonVector_late[i] << " / " << TotalPhotonVector(i) << std::endl;
 }

string const & art::Observer::processName ( ) const

protectedinherited

Definition at line 57 of file Observer.cc.

References art::Observer::process_name_.

Referenced by art::FileDumperOutput::printPrincipal().

   {
     return process_name_;
   }

float opdet::SimPhotonCounter::PromptPhotonTotal ( ) const

inline

Definition at line 89 of file SimPhotonCounter.h.

Referenced by opdet::FlashHypothesisComparison::FillSimPhotonCounterInfo().

     {
       return std::accumulate(_photonVector_prompt.begin(), _photonVector_prompt.end(), 0.0);
     }

const std::vector<float>& opdet::SimPhotonCounter::PromptPhotonVector ( ) const

inline

Definition at line 78 of file SimPhotonCounter.h.

Referenced by opdet::FlashHypothesisComparison::FillComparisonInfo(), opdet::FlashHypothesisComparison::FillSimPhotonCounterInfo(), opdet::FlashHypothesisComparison::RunComparison(), and TotalPhotonVector().

78 { return _photonVector_prompt; }

opdet::SimPhotonCounter::_photonVector_prompt

std::vector< float > _photonVector_prompt

Definition: SimPhotonCounter.h:102

float opdet::SimPhotonCounter::PromptPhotonVector ( size_t i ) const

inline

Definition at line 80 of file SimPhotonCounter.h.

80 { return _photonVector_prompt.at(i); }

opdet::SimPhotonCounter::_photonVector_prompt

std::vector< float > _photonVector_prompt

Definition: SimPhotonCounter.h:102

float opdet::SimPhotonCounter::QE ( size_t i ) const

inline

Definition at line 65 of file SimPhotonCounter.h.

65 { return _qeVector.at(i); }

opdet::SimPhotonCounter::_qeVector

std::vector< float > _qeVector

Definition: SimPhotonCounter.h:109

std::vector<float> const& opdet::SimPhotonCounter::QEVector ( ) const

inline

Definition at line 72 of file SimPhotonCounter.h.

72 { return _qeVector; }

opdet::SimPhotonCounter::_qeVector

std::vector< float > _qeVector

Definition: SimPhotonCounter.h:109

void art::ModuleBase::setModuleDescription ( ModuleDescription const & md )

inherited

Definition at line 31 of file ModuleBase.cc.

References art::ModuleBase::md_.

   {
     md_ = md;
   }

void opdet::SimPhotonCounter::SetQE	(	size_t	i,
		float	e
	)

inline

Definition at line 64 of file SimPhotonCounter.h.

References e.

64 { _qeVector.at(i) = e; }

opdet::SimPhotonCounter::_qeVector

std::vector< float > _qeVector

Definition: SimPhotonCounter.h:109

e

Float_t e

Definition: plot.C:35

void opdet::SimPhotonCounter::SetQEVector ( const std::vector< float > & eV )

inline

Definition at line 67 of file SimPhotonCounter.h.

     {
       SetVectorSize(eV.size());
       _qeVector = eV;
     }

void opdet::SimPhotonCounter::SetTimeRanges	(	float	t_p1,
		float	t_p2,
		float	t_l1,
		float	t_l2
	)

Definition at line 60 of file SimPhotonCounter.cxx.

References _max_late_time, _max_prompt_time, _min_late_time, and _min_prompt_time.

Referenced by SimPhotonCounter().

 {
   if (t_p2 < t_p1 || t_l2 < t_l1 || t_p2 > t_l1)
     throw std::runtime_error("ERROR in SimPhotonCounter: bad time ranges");
 
   _min_prompt_time = t_p1;
   _max_prompt_time = t_p2;
   _min_late_time = t_l1;
   _max_late_time = t_l2;
 }

void opdet::SimPhotonCounter::SetVectorSize ( size_t s )

inline

Definition at line 44 of file SimPhotonCounter.h.

     {
       _photonVector_prompt.resize(s);
       _photonVector_late.resize(s);
       _qeVector.resize(s);
     }

void opdet::SimPhotonCounter::SetWavelengthRanges	(	float	min_w,
		float	max_w
	)

Definition at line 44 of file SimPhotonCounter.cxx.

References _max_wavelength, and _min_wavelength.

Referenced by SimPhotonCounter().

 {
   if (min_w >= max_w) throw std::runtime_error("ERROR in SimPhotonCounter: bad wavelength range");
 
   _min_wavelength = min_w;
   _max_wavelength = max_w;
 }

void art::ModuleBase::sortConsumables ( std::string const & current_process_name )

inherited

Definition at line 49 of file ModuleBase.cc.

References art::ModuleBase::collector_, and art::ConsumesCollector::sortConsumables().

   {
     // Now that we know we have seen all the consumes declarations,
     // sort the results for fast lookup later.
     collector_.sortConsumables(current_process_name);
   }

void opdet::SimPhotonCounter::storeVisibility	(	int	channel,
		int	nDirectPhotons,
		int	nReflectedPhotons,
		double	reflectedT0 = `0.0`
	)		const

private

Definition at line 682 of file SimPhotonCounter_module.cc.

References DEFINE_ART_MODULE, fPVS, phot::PhotonVisibilityService::RetrieveLightProd(), phot::PhotonVisibilityService::SetLibraryEntry(), phot::PhotonVisibilityService::SetLibraryReflT0Entry(), phot::PhotonVisibilityService::StoreReflected(), and phot::PhotonVisibilityService::StoreReflT0().

Referenced by analyze().

   {
     phot::PhotonVisibilityService& pvs = *(art::ServiceHandle<phot::PhotonVisibilityService>());
 
     // ask PhotonVisibilityService which voxel was being served,
     // and how many photons where there generated (yikes!!);
     // this value was put there by LightSource (yikes!!)
     int VoxID;
     double NProd;
     fPVS->RetrieveLightProd(VoxID, NProd);
 
     pvs.SetLibraryEntry(VoxID, channel, double(nDirectPhotons) / NProd);
 
     //store reflected light
     if (fPVS->StoreReflected()) {
       pvs.SetLibraryEntry(VoxID, channel, double(nReflectedPhotons) / NProd, true);
 
       //store reflected first arrival time
       if (fPVS->StoreReflT0()) pvs.SetLibraryReflT0Entry(VoxID, channel, reflectedT0);
 
     } // if reflected
 
   } // SimPhotonCounter::storeVisibility()

std::vector< float > opdet::SimPhotonCounter::TotalPhotonVector ( ) const

Definition at line 98 of file SimPhotonCounter.cxx.

References util::begin(), util::end(), GetVectorSize(), LatePhotonVector(), and PromptPhotonVector().

Referenced by opdet::FlashHypothesisComparison::FillComparisonInfo(), opdet::FlashHypothesisComparison::FillSimPhotonCounterInfo(), and Print().

 {
 
   std::vector<float> totalPhotonVector(GetVectorSize());
   std::transform(PromptPhotonVector().begin(),
                  PromptPhotonVector().end(),
                  LatePhotonVector().begin(),
                  totalPhotonVector.begin(),
                  std::plus<float>());
   return totalPhotonVector;
 }

float opdet::SimPhotonCounter::TotalPhotonVector ( size_t i ) const

inline

Definition at line 84 of file SimPhotonCounter.h.

     {
       return (PromptPhotonVector(i) + LatePhotonVector(i));
     }

bool art::Observer::wantAllEvents ( ) const

inlineprotectednoexceptinherited

Definition at line 31 of file Observer.h.

References e, art::Observer::getTriggerResults(), art::Observer::wantAllEvents_, and art::Observer::wantEvent().

     {
       return wantAllEvents_;
     }

bool art::Observer::wantEvent	(	ScheduleID	id,
		Event const &	e
	)		const

protectedinherited

Definition at line 63 of file Observer.cc.

References art::Observer::rejectors_, art::Observer::selectors_, and art::Observer::wantAllEvents_.

Referenced by art::OutputModule::doEvent(), art::OutputModule::doWriteEvent(), and art::Observer::wantAllEvents().

   {
     if (wantAllEvents_) {
       return true;
     }
     bool const select_event = selectors_ ? selectors_->matchEvent(id, e) : true;
     bool const reject_event =
       rejectors_ ? rejectors_->matchEvent(id, e) : false;
     return select_event and not reject_event;
   }

float opdet::SimPhotonCounter::Wavelength ( const sim::OnePhoton & ph )

private

Definition at line 52 of file SimPhotonCounter.cxx.

References sim::OnePhoton::Energy.

Referenced by AddOnePhoton().

 {
   if (ph.Energy < std::numeric_limits<float>::epsilon())
     throw std::runtime_error("ERROR in SimPhotonCounter: photon energy is zero.");
 
   return 0.00124 / ph.Energy;
 }

Member Data Documentation

float opdet::SimPhotonCounter::_max_late_time

private

Definition at line 108 of file SimPhotonCounter.h.

Referenced by AddOnePhoton(), and SetTimeRanges().

float opdet::SimPhotonCounter::_max_prompt_time

private

Definition at line 106 of file SimPhotonCounter.h.

Referenced by AddOnePhoton(), and SetTimeRanges().

float opdet::SimPhotonCounter::_max_wavelength

private

Definition at line 112 of file SimPhotonCounter.h.

Referenced by AddOnePhoton(), and SetWavelengthRanges().

float opdet::SimPhotonCounter::_min_late_time

private

Definition at line 107 of file SimPhotonCounter.h.

Referenced by AddOnePhoton(), and SetTimeRanges().

float opdet::SimPhotonCounter::_min_prompt_time

private

Definition at line 105 of file SimPhotonCounter.h.

Referenced by AddOnePhoton(), and SetTimeRanges().

float opdet::SimPhotonCounter::_min_wavelength

private

Definition at line 111 of file SimPhotonCounter.h.

Referenced by AddOnePhoton(), and SetWavelengthRanges().

std::vector<float> opdet::SimPhotonCounter::_photonVector_late

private

Definition at line 103 of file SimPhotonCounter.h.

Referenced by AddOnePhoton(), ClearVectors(), Print(), and SimPhotonCounter().

std::vector<float> opdet::SimPhotonCounter::_photonVector_prompt

private

Definition at line 102 of file SimPhotonCounter.h.

Referenced by AddOnePhoton(), ClearVectors(), Print(), and SimPhotonCounter().

std::vector<float> opdet::SimPhotonCounter::_qeVector

private

Definition at line 109 of file SimPhotonCounter.h.

Referenced by AddOnePhoton(), Print(), and SimPhotonCounter().

Int_t opdet::SimPhotonCounter::fCountEventAll

private

Definition at line 130 of file SimPhotonCounter_module.cc.

Referenced by analyze(), and beginJob().

Int_t opdet::SimPhotonCounter::fCountEventDetected

private

Definition at line 131 of file SimPhotonCounter_module.cc.

Referenced by analyze(), and beginJob().

Int_t opdet::SimPhotonCounter::fCountOpDetAll

private

Definition at line 125 of file SimPhotonCounter_module.cc.

Referenced by analyze(), and beginJob().

Int_t opdet::SimPhotonCounter::fCountOpDetDetected

private

Definition at line 126 of file SimPhotonCounter_module.cc.

Referenced by analyze(), and beginJob().

Int_t opdet::SimPhotonCounter::fCountOpDetReflDetected

private

Definition at line 127 of file SimPhotonCounter_module.cc.

Referenced by analyze(), and beginJob().

double opdet::SimPhotonCounter::fdEdx

private

Definition at line 144 of file SimPhotonCounter_module.cc.

Referenced by analyze(), and beginJob().

double opdet::SimPhotonCounter::fEnergy

private

Definition at line 144 of file SimPhotonCounter_module.cc.

Referenced by analyze(), and beginJob().

Int_t opdet::SimPhotonCounter::fEventID

private

Definition at line 134 of file SimPhotonCounter_module.cc.

Referenced by analyze(), and beginJob().

TVector3 opdet::SimPhotonCounter::finalPhotonPosition

private

Definition at line 115 of file SimPhotonCounter_module.cc.

std::vector<std::string> opdet::SimPhotonCounter::fInputModule

private

Definition at line 99 of file SimPhotonCounter_module.cc.

Referenced by analyze(), and isVisible().

TTree* opdet::SimPhotonCounter::fLightAnalysisTree = nullptr

private

Definition at line 142 of file SimPhotonCounter_module.cc.

Referenced by analyze(), and beginJob().

bool opdet::SimPhotonCounter::fMakeAllPhotonsTree

private

Definition at line 104 of file SimPhotonCounter_module.cc.

Referenced by analyze(), beginJob(), and isVisible().

bool opdet::SimPhotonCounter::fMakeDetectedPhotonsTree

private

Definition at line 103 of file SimPhotonCounter_module.cc.

Referenced by analyze(), beginJob(), and isVisible().

bool opdet::SimPhotonCounter::fMakeLightAnalysisTree

private

Definition at line 138 of file SimPhotonCounter_module.cc.

Referenced by analyze(), beginJob(), and isVisible().

bool opdet::SimPhotonCounter::fMakeOpDetEventsTree

private

Definition at line 106 of file SimPhotonCounter_module.cc.

Referenced by analyze(), beginJob(), and isVisible().

bool opdet::SimPhotonCounter::fMakeOpDetsTree

private

Definition at line 105 of file SimPhotonCounter_module.cc.

Referenced by analyze(), beginJob(), and isVisible().

int opdet::SimPhotonCounter::fmotherTrackID

private

Definition at line 143 of file SimPhotonCounter_module.cc.

Referenced by analyze(), and beginJob().

Int_t opdet::SimPhotonCounter::fOpChannel

private

Definition at line 135 of file SimPhotonCounter_module.cc.

Referenced by analyze(), and beginJob().

Float_t opdet::SimPhotonCounter::foriginX

private

Definition at line 121 of file SimPhotonCounter_module.cc.

Referenced by analyze(), and beginJob().

Float_t opdet::SimPhotonCounter::foriginY

private

Definition at line 122 of file SimPhotonCounter_module.cc.

Referenced by analyze(), and beginJob().

Float_t opdet::SimPhotonCounter::foriginZ

private

Definition at line 123 of file SimPhotonCounter_module.cc.

Referenced by analyze(), and beginJob().

int opdet::SimPhotonCounter::fpdg

private

Definition at line 143 of file SimPhotonCounter_module.cc.

Referenced by analyze(), and beginJob().

std::vector<double> opdet::SimPhotonCounter::fPosition0

private

Definition at line 145 of file SimPhotonCounter_module.cc.

std::string opdet::SimPhotonCounter::fProcess

private

Definition at line 150 of file SimPhotonCounter_module.cc.

Referenced by analyze(), and beginJob().

phot::PhotonVisibilityService const* opdet::SimPhotonCounter::fPVS = nullptr

private

Definition at line 153 of file SimPhotonCounter_module.cc.

Referenced by analyze(), beginJob(), endJob(), isVisible(), and storeVisibility().

int opdet::SimPhotonCounter::fRun

private

Definition at line 143 of file SimPhotonCounter_module.cc.

Referenced by analyze(), and beginJob().

bool opdet::SimPhotonCounter::fSavePhotonOrigins

private

Definition at line 107 of file SimPhotonCounter_module.cc.

Referenced by beginJob(), and isVisible().

std::vector<std::vector<std::vector<double> > > opdet::SimPhotonCounter::fSignals_vis

private

Definition at line 140 of file SimPhotonCounter_module.cc.

Referenced by analyze().

std::vector<std::vector<std::vector<double> > > opdet::SimPhotonCounter::fSignals_vuv

private

Definition at line 139 of file SimPhotonCounter_module.cc.

Referenced by analyze().

std::vector<std::vector<double> > opdet::SimPhotonCounter::fSignalsvis

private

Definition at line 149 of file SimPhotonCounter_module.cc.

Referenced by analyze(), and beginJob().

std::vector<std::vector<double> > opdet::SimPhotonCounter::fSignalsvuv

private

Definition at line 148 of file SimPhotonCounter_module.cc.

Referenced by analyze(), and beginJob().

std::vector<std::vector<double> > opdet::SimPhotonCounter::fstepPositions

private

Definition at line 146 of file SimPhotonCounter_module.cc.

Referenced by analyze(), and beginJob().

std::vector<double> opdet::SimPhotonCounter::fstepTimes

private

Definition at line 147 of file SimPhotonCounter_module.cc.

Referenced by analyze(), and beginJob().

Float_t opdet::SimPhotonCounter::fT0_vis

private

Definition at line 128 of file SimPhotonCounter_module.cc.

Referenced by analyze().

TTree* opdet::SimPhotonCounter::fTheEventTree

private

Definition at line 95 of file SimPhotonCounter_module.cc.

Referenced by analyze(), and beginJob().

TTree* opdet::SimPhotonCounter::fTheOpDetTree

private

Definition at line 94 of file SimPhotonCounter_module.cc.

Referenced by analyze(), and beginJob().

TTree* opdet::SimPhotonCounter::fThePhotonTreeAll

private

Definition at line 92 of file SimPhotonCounter_module.cc.

Referenced by analyze(), and beginJob().

TTree* opdet::SimPhotonCounter::fThePhotonTreeDetected

private

Definition at line 93 of file SimPhotonCounter_module.cc.

Referenced by analyze(), and beginJob().

Float_t opdet::SimPhotonCounter::fTime

private

Definition at line 120 of file SimPhotonCounter_module.cc.

Referenced by analyze(), and beginJob().

int opdet::SimPhotonCounter::fTrackID

private

Definition at line 143 of file SimPhotonCounter_module.cc.

Referenced by analyze(), and beginJob().

bool const opdet::SimPhotonCounter::fUseLitePhotons

private

Definition at line 155 of file SimPhotonCounter_module.cc.

Referenced by analyze(), and isVisible().

int opdet::SimPhotonCounter::fVerbosity

private

Definition at line 101 of file SimPhotonCounter_module.cc.

Referenced by analyze(), and isVisible().

Float_t opdet::SimPhotonCounter::fWavelength

private

Definition at line 119 of file SimPhotonCounter_module.cc.

Referenced by analyze(), and beginJob().

TVector3 opdet::SimPhotonCounter::initialPhotonPosition

private

Definition at line 114 of file SimPhotonCounter_module.cc.

constexpr double opdet::SimPhotonCounter::kVisibleThreshold = 200.0

staticprivate

Threshold used to resolve between visible and ultraviolet light.

Definition at line 82 of file SimPhotonCounter_module.cc.

Referenced by analyze(), and isVisible().

constexpr double opdet::SimPhotonCounter::kVisibleWavelength = 450.0

staticprivate

Value used when a typical visible light wavelength is needed.

Definition at line 85 of file SimPhotonCounter_module.cc.

constexpr double opdet::SimPhotonCounter::kVUVWavelength = 128.0

staticprivate

Value used when a typical ultraviolet light wavelength is needed.

Definition at line 88 of file SimPhotonCounter_module.cc.

Referenced by analyze().

cheat::ParticleInventoryService const* opdet::SimPhotonCounter::pi_serv = nullptr

private

Definition at line 152 of file SimPhotonCounter_module.cc.

Referenced by analyze(), and beginJob().

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

larana/v10_00_11/source/larana/OpticalDetector/SimPhotonCounter.h
larana/v10_00_11/source/larana/OpticalDetector/SimPhotonCounter_module.cc
larana/v10_00_11/source/larana/OpticalDetector/SimPhotonCounter.cxx

Public Types

Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Static Private Attributes

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation