A module for optical MC testing and library building. More...

Inheritance diagram for evgen::LightSource:

Public Types
using	ModuleType = EDProducer

using	WorkerType = WorkerT< EDProducer >

template<typename UserConfig , typename KeysToIgnore = void>
using	Table = ProducerBase::Table< UserConfig, KeysToIgnore >

Public Member Functions
	LightSource (fhicl::ParameterSet const &pset)

void	produce (art::Event &evt)

void	beginRun (art::Run &run)

template<typename PROD , BranchType B = InEvent>
ProductID	getProductID (std::string const &instanceName={}) const

template<typename PROD , BranchType B>
ProductID	getProductID (ModuleDescription const &moduleDescription, std::string const &instanceName) const

bool	modifiesEvent () const

template<typename T , BranchType = InEvent>
ProductToken< T >	consumes (InputTag const &)

template<typename T , art::BranchType BT>
art::ProductToken< T >	consumes (InputTag const &it)

template<typename T , BranchType = InEvent>
void	consumesMany ()

template<typename Element , BranchType = InEvent>
ViewToken< Element >	consumesView (InputTag const &)

template<typename T , art::BranchType BT>
art::ViewToken< T >	consumesView (InputTag const &it)

template<typename T , BranchType = InEvent>
ProductToken< T >	mayConsume (InputTag const &)

template<typename T , art::BranchType BT>
art::ProductToken< T >	mayConsume (InputTag const &it)

template<typename T , BranchType = InEvent>
void	mayConsumeMany ()

template<typename Element , BranchType = InEvent>
ViewToken< Element >	mayConsumeView (InputTag const &)

template<typename T , art::BranchType BT>
art::ViewToken< T >	mayConsumeView (InputTag const &it)

base_engine_t &	createEngine (seed_t seed)

base_engine_t &	createEngine (seed_t seed, std::string const &kind_of_engine_to_make)

base_engine_t &	createEngine (seed_t seed, std::string const &kind_of_engine_to_make, label_t const &engine_label)

seed_t	get_seed_value (fhicl::ParameterSet const &pset, char const key[]="seed", seed_t const implicit_seed=-1)

Static Public Member Functions
static cet::exempt_ptr< Consumer >	non_module_context ()

Protected Member Functions
CurrentProcessingContext const *	currentContext () const

void	validateConsumedProduct (BranchType const bt, ProductInfo const &pi)

void	prepareForJob (fhicl::ParameterSet const &pset)

void	showMissingConsumes () const

Private Member Functions
void	Sample (simb::MCTruth &truth)

Private Attributes
std::string	fVersion

std::ifstream	fInputFile

std::string	fFileName

char	fDummyString [256]

TTree *	fPhotonsGenerated

TLorentzVector	fShotPos

TLorentzVector	fShotMom

Int_t	fEvID

int	fSourceMode

bool	fFillTree

int	fPosDist

int	fTDist

int	fPDist

int	fXSteps

int	fYSteps

int	fZSteps

sim::PhotonVoxelDef	fThePhotonVoxelDef

int	fVoxelCount

int	fCurrentVoxel

TVector3	fTPCCenter

TVector3	fTPCDimensions

std::vector< double >	fRegionMin

std::vector< double >	fRegionMax

bool	fUseCustomRegion

double	fX

double	fY

double	fZ

bool	fPointSource

double	fT

double	fSigmaX

double	fSigmaY

double	fSigmaZ

double	fSigmaT

double	fP

double	fSigmaP

int	fN

int	fFirstVoxel

int	fLastVoxel

Static Private Attributes
static const int	kUNIF = 0

static const int	kGAUS = 1

static const int	kFILE = 0

static const int	kSCAN = 1

Detailed Description

A module for optical MC testing and library building.

Light source event generator which simulate an extended isotropic photon source.

The light source can be run in two modes, file mode or scan mode. Each requires the specification of a different set of parameters.

File mode

Light source position, intensity and shape are supplied on an event by event basis in a text file. See the example provided for the format. Pararmeters required:

int32   SourceMode = 0      - sets light source to file mode
string  FileName            - file of per event light source specifications
int32   PosDist             - how to distribute production points sampled in momentum, position
int32   PDist                   and time ranges specified.  For all of these :
int32   TDist                   0 = uniform and 1 = gauss
bool    FillTree            - whether to write a tree of photon production points to fileservice

Upon reaching the end of the file, the light source will loop back to the first point. hence a one line text file will give a constant light source size, position and intensity.

Scan mode

Divide volume into cuboidal regions and produce an isotropic light source in each, using one region per event. User can specify either to use the full detector volume or some custom specified volume.

This mode is used when building a fast photon sim library, and performing volume scan sensitivity studies.

int32   SourceMode = 1      - sets light source to scan mode
int32   N                   - number of photons to shoot from each point
double  P                   - peak photon momentum (or energy) in eV
double  SigmaP              - momentum distribution width
double  XSteps              - Number of regions to divide volume into in each direction 
double  YSteps
double  ZSteps
double  T0                  - Peak time of photon production
double  SigmaT              - time distribution width
int32   PosDist             - how to distribute production points sampled in momentum, position
int32   PDist                 and time ranges specified.  For all of these :
int32   TDist                   0 = uniform and 1 = gaussian
bool    FillTree            - whether to write a tree of photon production points to fileservice
bool    UseCustomRegion     - supply our own volme specification or use the full detector volume?
vdouble[3]  RegionMin       - bounding corners of the custom volume specification 
vdouble[3]  RegionMax           (only used if UseCustomRegion=true)

Definition at line 106 of file LightSource_module.cc.

Member Typedef Documentation

using art::EDProducer::ModuleType = EDProducer

inherited

Definition at line 34 of file EDProducer.h.

template<typename UserConfig , typename KeysToIgnore = void>

using art::EDProducer::Table = ProducerBase::Table<UserConfig, KeysToIgnore>

inherited

Definition at line 43 of file EDProducer.h.

using art::EDProducer::WorkerType = WorkerT<EDProducer>

inherited

Definition at line 35 of file EDProducer.h.

Constructor & Destructor Documentation

evgen::LightSource::LightSource ( fhicl::ParameterSet const & pset )

explicit

Definition at line 188 of file LightSource_module.cc.

References art::EngineCreator::createEngine(), fCurrentVoxel, fDummyString, fEvID, fFileName, fFillTree, fFirstVoxel, fInputFile, fLastVoxel, fN, fP, fPDist, fPhotonsGenerated, fPointSource, fPosDist, fRegionMax, fRegionMin, fShotMom, fShotPos, fSigmaP, fSigmaT, fSigmaX, fSigmaY, fSigmaZ, fSourceMode, fT, fTDist, fThePhotonVoxelDef, fUseCustomRegion, fVoxelCount, fXSteps, fYSteps, fZSteps, fhicl::ParameterSet::get(), sim::PhotonVoxelDef::GetNVoxels(), phot::PhotonVisibilityService::GetVoxelDef(), sim::PhotonVoxelDef::GetVoxelSize(), kFILE, kSCAN, and art::TFileDirectory::make().

   {
 
     // get the random number seed, use a random default if not specified    
     // in the configuration file.  
  
     fSourceMode   =     (pset.get<int >("SourceMode")  );
     fFillTree     =     (pset.get<bool>("FillTree")    );
     fPosDist      =     (pset.get<int >("PosDist")     );
     fPDist        =     (pset.get<int >("PDist")       );
     fTDist        =     (pset.get<int >("TDist")       );
     
     // create a default random engine; obtain the random seed from NuRandomService,
     // unless overridden in configuration with key "Seed"
     art::ServiceHandle<rndm::NuRandomService>()
       ->createEngine(*this, pset, "Seed");
 
     // load optional parameters in function
     produces< sumdata::RunData, art::InRun >();
     produces< std::vector<simb::MCTruth> >();
  
     if(fSourceMode==kFILE)
       {
         fFileName  = pset.get<std::string>("SteeringFile");
         fInputFile.open(fFileName.c_str());
         fInputFile.getline(fDummyString,256);
 
       }
     else if (fSourceMode==kSCAN)
       {
         fT      = pset.get<double>("T0");
         fSigmaT = pset.get<double>("SigmaT");
         fN      = pset.get<int   >("N");
 
         fFirstVoxel     = pset.get<int   >("FirstVoxel");
         fLastVoxel      = pset.get<int   >("LastVoxel");
         
         fP      = pset.get<double>("P");
         fSigmaP = pset.get<double>("SigmaP");
 
         fUseCustomRegion = pset.get<bool>("UseCustomRegion");
         fPointSource = pset.get<bool>("PointSource",false);
 
         if(fUseCustomRegion)
           {     
             fRegionMin = pset.get< std::vector<double> >("RegionMin");
             fRegionMax = pset.get< std::vector<double> >("RegionMax");
             fXSteps = pset.get<int >("XSteps");
             fYSteps = pset.get<int >("YSteps");
             fZSteps = pset.get<int >("ZSteps");
           }
 
         art::ServiceHandle<geo::Geometry> geo;
         // get TPC dimensions removed. -TA
 
 
         fCurrentVoxel=0;
 
         // define voxelization based on parameters read from config.
         // There are two modes - either read the dimensions of the TPC from
         // the geometry, or use values specified by the user.
         if(!fUseCustomRegion)
           {
             art::ServiceHandle<phot::PhotonVisibilityService> vis;   
             fThePhotonVoxelDef = vis->GetVoxelDef();
           }
         else
           {
             fThePhotonVoxelDef = sim::PhotonVoxelDef(fRegionMin[0], 
                                                      fRegionMax[0],
                                                      fXSteps,
                                                      fRegionMin[1],
                                                      fRegionMax[1],
                                                      fYSteps,
                                                      fRegionMin[2],
                                                      fRegionMax[2],
                                                      fZSteps);   
           }
         
 
         // Set distribution widths to voxel size
 
         fSigmaX = fThePhotonVoxelDef.GetVoxelSize().X()/2.0;
         fSigmaY = fThePhotonVoxelDef.GetVoxelSize().Y()/2.0;
         fSigmaZ = fThePhotonVoxelDef.GetVoxelSize().Z()/2.0;
 
         // Get number of voxels we will step through
         
         fVoxelCount = fThePhotonVoxelDef.GetNVoxels();
         
         if(fLastVoxel<0) fLastVoxel = fVoxelCount;
         
         mf::LogVerbatim("LightSource") << "Light Source : Determining voxel params : " 
                                        << fVoxelCount << " " 
                                        << fSigmaX     << " "
                                        << fSigmaY     << " " 
                                        <<fSigmaZ;
 
       }
     else{
       throw cet::exception("LightSource") << "EVGEN Light Source : Unrecognised light source mode\n";
     }
     
 
 
 
     if(fFillTree)
       {
         art::ServiceHandle<art::TFileService> tfs;
         fPhotonsGenerated = tfs->make<TTree>("PhotonsGenerated","PhotonsGenerated");
         fPhotonsGenerated->Branch("X",&(fShotPos[0]),"X/D");
         fPhotonsGenerated->Branch("Y",&(fShotPos[1]),"Y/D");
         fPhotonsGenerated->Branch("Z",&(fShotPos[2]),"Z/D");
         fPhotonsGenerated->Branch("T",&(fShotPos[3]),"T/D");
         fPhotonsGenerated->Branch("PX",&(fShotMom[0]),"PX/D");
         fPhotonsGenerated->Branch("PY",&(fShotMom[1]),"PY/D");
         fPhotonsGenerated->Branch("PZ",&(fShotMom[2]),"PZ/D");
         fPhotonsGenerated->Branch("PT",&(fShotMom[3]),"PT/D");
         fPhotonsGenerated->Branch("EventID",&fEvID,"EventID/I");
       }
   }

Member Function Documentation

void evgen::LightSource::beginRun ( art::Run & run )

virtual

Reimplemented from art::EDProducer.

Definition at line 312 of file LightSource_module.cc.

References geo::GeometryCore::DetectorName(), fCurrentVoxel, fFirstVoxel, and art::Run::put().

   {
 
     // grab the geometry object to see what geometry we are using
     art::ServiceHandle<geo::Geometry> geo;
     std::unique_ptr<sumdata::RunData> runcol(new sumdata::RunData(geo->DetectorName()));
 
     run.put(std::move(runcol));
     
     fCurrentVoxel=fFirstVoxel;
 
     return;
   }

template<typename T , BranchType = InEvent>

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

inherited

template<typename T , art::BranchType BT>

art::ProductToken<T> art::Consumer::consumes ( InputTag const & it )

inherited

Definition at line 147 of file Consumer.h.

References art::InputTag::instance(), art::InputTag::label(), and art::InputTag::process().

 {
   if (!moduleContext_)
     return ProductToken<T>::invalid();
 
   consumables_[BT].emplace_back(ConsumableType::Product,
                                 TypeID{typeid(T)},
                                 it.label(),
                                 it.instance(),
                                 it.process());
   return ProductToken<T>{it};
 }

template<typename T , art::BranchType BT>

void art::Consumer::consumesMany ( )

inherited

Definition at line 162 of file Consumer.h.

 {
   if (!moduleContext_)
     return;
 
   consumables_[BT].emplace_back(ConsumableType::Many, TypeID{typeid(T)});
 }

template<typename Element , BranchType = InEvent>

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

inherited

template<typename T , art::BranchType BT>

art::ViewToken<T> art::Consumer::consumesView ( InputTag const & it )

inherited

Definition at line 172 of file Consumer.h.

References art::InputTag::instance(), art::InputTag::label(), and art::InputTag::process().

 {
   if (!moduleContext_)
     return ViewToken<T>::invalid();
 
   consumables_[BT].emplace_back(ConsumableType::ViewElement,
                                 TypeID{typeid(T)},
                                 it.label(),
                                 it.instance(),
                                 it.process());
   return ViewToken<T>{it};
 }

EngineCreator::base_engine_t & EngineCreator::createEngine ( seed_t seed )

inherited

Definition at line 26 of file EngineCreator.cc.

References art::EngineCreator::rng().

Referenced by evgen::CosmicsGen::CosmicsGen(), rndm::NuRandomService::createEngine(), cluster::fuzzyCluster::fuzzyCluster(), cluster::HoughLineFinder::HoughLineFinder(), art::MixFilter< T >::initEngine_(), larg4::LArG4::LArG4(), LightSource(), evgen::NeutronOsc::NeutronOsc(), evgen::NucleonDecay::NucleonDecay(), opdet::OpMCDigi::OpMCDigi(), opdet::OptDetDigitizer::OptDetDigitizer(), phot::PhotonLibraryPropagation::PhotonLibraryPropagation(), detsim::SimDriftElectrons::SimDriftElectrons(), evgen::SingleGen::SingleGen(), evgen::SNNueAr40CCGen::SNNueAr40CCGen(), ToyOneShowerGen::ToyOneShowerGen(), and trkf::Track3DKalman::Track3DKalman().

 {
   return rng()->createEngine(placeholder_schedule_id(), seed);
 }

EngineCreator::base_engine_t & EngineCreator::createEngine	(	seed_t	seed,
		std::string const &	kind_of_engine_to_make
	)

inherited

Definition at line 32 of file EngineCreator.cc.

References art::EngineCreator::rng().

 {
   return rng()->createEngine(
     placeholder_schedule_id(), seed, kind_of_engine_to_make);
 }

EngineCreator::base_engine_t & EngineCreator::createEngine	(	seed_t	seed,
		std::string const &	kind_of_engine_to_make,
		label_t const &	engine_label
	)

inherited

Definition at line 40 of file EngineCreator.cc.

References art::EngineCreator::rng().

 {
   return rng()->createEngine(
     placeholder_schedule_id(), seed, kind_of_engine_to_make, engine_label);
 }

CurrentProcessingContext const * art::EDProducer::currentContext ( ) const

protectedinherited

Definition at line 120 of file EDProducer.cc.

References art::EDProducer::current_context_.

   {
     return current_context_.get();
   }

EngineCreator::seed_t EngineCreator::get_seed_value	(	fhicl::ParameterSet const &	pset,
		char const	key[] = `"seed"`,
		seed_t const	implicit_seed = `-1`
	)

inherited

Definition at line 49 of file EngineCreator.cc.

References fhicl::ParameterSet::get().

Referenced by art::MixFilter< T >::initEngine_().

 {
   auto const& explicit_seeds = pset.get<std::vector<int>>(key, {});
   return explicit_seeds.empty() ? implicit_seed : explicit_seeds.front();
 }

template<typename PROD , BranchType B>

ProductID art::EDProducer::getProductID ( std::string const & instanceName = {} ) const

inlineinherited

Definition at line 123 of file EDProducer.h.

References art::EDProducer::moduleDescription_.

   {
     return ProducerBase::getProductID<PROD, B>(moduleDescription_,
                                                instanceName);
   }

template<typename PROD , BranchType B>

ProductID art::ProducerBase::getProductID	(	ModuleDescription const &	moduleDescription,
		std::string const &	instanceName
	)		const

inherited

Definition at line 56 of file ProducerBase.h.

References B, and art::ModuleDescription::moduleLabel().

Referenced by art::ProducerBase::modifiesEvent().

   {
     auto const& pd =
       get_ProductDescription<PROD>(B, md.moduleLabel(), instanceName);
     return pd.productID();
   }

template<typename T , BranchType = InEvent>

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

inherited

template<typename T , art::BranchType BT>

art::ProductToken<T> art::Consumer::mayConsume ( InputTag const & it )

inherited

Definition at line 190 of file Consumer.h.

References art::InputTag::instance(), art::InputTag::label(), and art::InputTag::process().

 {
   if (!moduleContext_)
     return ProductToken<T>::invalid();
 
   consumables_[BT].emplace_back(ConsumableType::Product,
                                 TypeID{typeid(T)},
                                 it.label(),
                                 it.instance(),
                                 it.process());
   return ProductToken<T>{it};
 }

template<typename T , art::BranchType BT>

void art::Consumer::mayConsumeMany ( )

inherited

Definition at line 205 of file Consumer.h.

 {
   if (!moduleContext_)
     return;
 
   consumables_[BT].emplace_back(ConsumableType::Many, TypeID{typeid(T)});
 }

template<typename Element , BranchType = InEvent>

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

inherited

template<typename T , art::BranchType BT>

art::ViewToken<T> art::Consumer::mayConsumeView ( InputTag const & it )

inherited

Definition at line 215 of file Consumer.h.

References art::InputTag::instance(), art::InputTag::label(), and art::InputTag::process().

 {
   if (!moduleContext_)
     return ViewToken<T>::invalid();
 
   consumables_[BT].emplace_back(ConsumableType::ViewElement,
                                 TypeID{typeid(T)},
                                 it.label(),
                                 it.instance(),
                                 it.process());
   return ViewToken<T>{it};
 }

bool art::ProducerBase::modifiesEvent ( ) const

inlineinherited

Definition at line 40 of file ProducerBase.h.

References art::ProducerBase::getProductID().

     {
       return true;
     }

cet::exempt_ptr< art::Consumer > art::Consumer::non_module_context ( )

staticinherited

Definition at line 76 of file Consumer.cc.

 {
   static Consumer invalid{InvalidTag{}};
   return &invalid;
 }

void art::Consumer::prepareForJob ( fhicl::ParameterSet const & pset )

protectedinherited

Definition at line 89 of file Consumer.cc.

References fhicl::ParameterSet::get_if_present().

Referenced by art::EDProducer::doBeginJob(), art::EDFilter::doBeginJob(), and art::EDAnalyzer::doBeginJob().

 {
   if (!moduleContext_)
     return;
 
   pset.get_if_present("errorOnMissingConsumes", requireConsumes_);
   for (auto& consumablesPerBranch : consumables_) {
     cet::sort_all(consumablesPerBranch);
   }
 }

void evgen::LightSource::produce ( art::Event & evt )

virtual

Implements art::EDProducer.

Definition at line 327 of file LightSource_module.cc.

References e, fCurrentVoxel, fDummyString, fFileName, fFirstVoxel, fInputFile, fLastVoxel, fN, fP, fSigmaP, fSigmaT, fSigmaX, fSigmaY, fSigmaZ, fSourceMode, fT, fThePhotonVoxelDef, fX, fY, fZ, sim::PhotonVoxel::GetCenter(), sim::PhotonVoxelDef::GetPhotonVoxel(), phot::PhotonVisibilityService::IsBuildJob(), kFILE, kSCAN, simb::kSingleParticle, art::Event::put(), Sample(), art::errors::ServiceNotFound, simb::MCTruth::SetOrigin(), and phot::PhotonVisibilityService::StoreLightProd().

   {
     
     // FILE MODE -
     //  Each event, read coordinates of gun and number of photons to shoot from file
     
     if(fSourceMode==kFILE)
       {
         // Loop file if required
         if(fInputFile.eof()){
           mf::LogWarning("LightSource") << "EVGEN Light Source : Warning, reached end of file,"
                                         << " looping back to beginning";
           fInputFile.seekg(0,std::ios::beg);
           fInputFile.clear();
         }
         
         if(!fInputFile.is_open() || fInputFile.fail() ){
           throw cet::exception("LightSource") << "EVGEN Light Source : File error in " 
                                               << fFileName << "\n";
         }
         else{ 
           // read in one line
           fInputFile >> fX >> fY >> fZ >> fT 
                      >> fSigmaX >> fSigmaY >> fSigmaZ >> fSigmaT 
                      >> fP >> fSigmaP >> fN;
           fInputFile.getline(fDummyString,256);
           fThePhotonVoxelDef = sim::PhotonVoxelDef(fX - fSigmaX, 
                                                    fX + fSigmaX,
                                                    1,
                                                    fY - fSigmaY,
                                                    fY + fSigmaY,
                                                    1,
                                                    fZ - fSigmaZ,
                                                    fZ + fSigmaZ,
                                                    1);   
             
           fCurrentVoxel=0;
         }
       }
 
 
     // SCAN MODE -
     //  Step through detector using a number of steps provided in the config file
     //  firing a constant number of photons from each point
     else if(fSourceMode==kSCAN)
       {
 
 
         TVector3 VoxelCenter = fThePhotonVoxelDef.GetPhotonVoxel(fCurrentVoxel).GetCenter();
         
         fX = VoxelCenter.X();
         fY = VoxelCenter.Y();
         fZ = VoxelCenter.Z();
         
       }
         
     
     // UNRECOGNISED MODE 
     //  - neither file or scan mode, probably a config file error
 
     else{
       throw cet::exception("LightSource") <<"EVGEN : Light Source, unrecognised source mode\n";
     }
     
 
 
     //    std::cout<<"EVGEN Light source to be placed at  (x, y, z, t, dx, dy, dz, dt, n) " <<
     //      fX << " " << fY << " " <<fZ << " " << fT << " " << 
     //      fSigmaX << " " << fSigmaY << " " << fSigmaZ << " " << fSigmaT << " " <<
     //      fN << std::endl<<std::endl;
     
       
     std::unique_ptr< std::vector<simb::MCTruth> > truthcol(new std::vector<simb::MCTruth>);
 
     simb::MCTruth truth;
     truth.SetOrigin(simb::kSingleParticle);
     Sample(truth);
     
     //     std::cout << "put mctruth into the vector" << std::endl;
     truthcol->push_back(truth);
 
     //     std::cout << "add vector to the event " << truthcol->size() << std::endl;
     evt.put(std::move(truthcol));
     
     phot::PhotonVisibilityService* vis = nullptr;
     try {
       vis = art::ServiceHandle<phot::PhotonVisibilityService>().get();
     }
     catch (art::Exception const& e) {
       // if the service is not configured, then this is not a build job
       // (it is a simple generation job instead)
       if (e.categoryCode() != art::errors::ServiceNotFound) throw;
     }
 
     if(vis && vis->IsBuildJob())
       {
         mf::LogVerbatim("LightSource") << "Light source : Stowing voxel params ";
         vis->StoreLightProd(fCurrentVoxel,fN);
       }
 
     if(fCurrentVoxel!=fLastVoxel) 
       {
         ++fCurrentVoxel;
       }
     else
       {
         mf::LogVerbatim("LightSource") << "EVGEN Light Source fully scanned detector.  Starting over.";
         fCurrentVoxel=fFirstVoxel;
       }
     
     return;
 
 
   }

void evgen::LightSource::Sample ( simb::MCTruth & truth )

private

Definition at line 443 of file LightSource_module.cc.

References simb::MCTruth::Add(), simb::MCParticle::AddTrajectoryPoint(), DEFINE_ART_MODULE, fFillTree, fN, fP, fPDist, fPhotonsGenerated, fPointSource, fPosDist, fShotMom, fShotPos, fSigmaP, fSigmaT, fSigmaX, fSigmaY, fSigmaZ, fT, fTDist, fX, fY, fZ, art::RandomNumberGenerator::getEngine(), kGAUS, part, art::EngineCreator::rng(), and x.

Referenced by produce().

   {
     mf::LogVerbatim("LightSource") <<"Light source debug : Shooting at " << fX <<" " << fY<<" "<< fZ;
     
     // get the random number generator service and make some CLHEP generators
     art::ServiceHandle<art::RandomNumberGenerator> rng;
     CLHEP::HepRandomEngine &engine = rng->getEngine();
     CLHEP::RandFlat   flat(engine);
     CLHEP::RandGaussQ gauss(engine);
 
     for(int j=0; j!=fN; ++j){
       // Choose momentum (supplied in eV, convert to GeV)
       double p = fP;
       if (fPDist == kGAUS) {
         p = gauss.fire(fP, fSigmaP);
       }
       else {
         p = fP + fSigmaP*(2.0*flat.fire()-1.0);
       }
       p /= 1000000000.;
       
       // Choose position
       TVector3 x;
       if(fPointSource) {
         x[0] = fX;
         x[1] = fY;
         x[2] = fZ;
       }
       else {
         if (fPosDist == kGAUS) {
           x[0] = gauss.fire(fX, fSigmaX);
           x[1] = gauss.fire(fY, fSigmaY);
           x[2] = gauss.fire(fZ, fSigmaZ);
         }
         else {
           x[0] = fX + fSigmaX*(2.0*flat.fire()-1.0);
           x[1] = fY + fSigmaY*(2.0*flat.fire()-1.0);
           x[2] = fZ + fSigmaZ*(2.0*flat.fire()-1.0);
         }
 
       }
       
       // Choose time
       double t;
       if (fTDist == kGAUS) {
         t = gauss.fire(fT, fSigmaT);
       }
       else {
         t = fT + fSigmaT * (2.0 * flat.fire()-1.0);
       }
         
         
       //assume the position is relative to the center of the TPC
       //x += fTPCCenter;
         
       fShotPos = TLorentzVector(x[0], x[1], x[2], t);
         
 
       // Choose angles
       double costh = 2 * flat.fire() - 1;
       double sinth = pow(1-pow(costh,2),0.5);
       double phi   = 2 * M_PI * flat.fire();
       
       // Generate momentum 4-vector
       
       fShotMom = TLorentzVector( p*sinth*cos(phi), 
                                  p*sinth*sin(phi), 
                                  p*costh, 
                                  p                  );
         
       int trackid = -1*(j+1); // set track id to -i as these are all primary particles and have id <= 0
       std::string primary("primary");
       int PDG=0; //optical photons have PDG 0
 
       simb::MCParticle part(trackid, PDG, primary);
       part.AddTrajectoryPoint(fShotPos, fShotMom);
 
       if(fFillTree)
         fPhotonsGenerated->Fill();
         
       mct.Add(part);
     }
 
   }

void art::Consumer::showMissingConsumes ( ) const

protectedinherited

Definition at line 125 of file Consumer.cc.

Referenced by art::EDProducer::doEndJob(), art::EDFilter::doEndJob(), art::EDAnalyzer::doEndJob(), and art::RootOutput::endJob().

 {
   if (!moduleContext_)
     return;
 
   // If none of the branches have missing consumes statements, exit early.
   if (std::all_of(cbegin(missingConsumes_),
                   cend(missingConsumes_),
                   [](auto const& perBranch) { return perBranch.empty(); }))
     return;
 
   constexpr cet::HorizontalRule rule{60};
   mf::LogPrint log{"MTdiagnostics"};
   log << '\n'
       << rule('=') << '\n'
       << "The following consumes (or mayConsume) statements are missing from\n"
       << module_context(moduleDescription_) << '\n'
       << rule('-') << '\n';
 
   cet::for_all_with_index(
     missingConsumes_, [&log](std::size_t const i, auto const& perBranch) {
       for (auto const& pi : perBranch) {
         log << "  "
             << assemble_consumes_statement(static_cast<BranchType>(i), pi)
             << '\n';
       }
     });
   log << rule('=');
 }

void art::Consumer::validateConsumedProduct	(	BranchType const	bt,
		ProductInfo const &	pi
	)

protectedinherited

Definition at line 101 of file Consumer.cc.

References art::errors::ProductRegistrationFailure.

 {
   // Early exits if consumes tracking has been disabled or if the
   // consumed product is an allowed consumable.
   if (!moduleContext_)
     return;
 
   if (cet::binary_search_all(consumables_[bt], pi))
     return;
 
   if (requireConsumes_) {
     throw Exception(errors::ProductRegistrationFailure,
                     "Consumer: an error occurred during validation of a "
                     "retrieved product\n\n")
       << "The following consumes (or mayConsume) statement is missing from\n"
       << module_context(moduleDescription_) << ":\n\n"
       << "  " << assemble_consumes_statement(bt, pi) << "\n\n";
   }
 
   missingConsumes_[bt].insert(pi);
 }

Member Data Documentation

int evgen::LightSource::fCurrentVoxel

private

Definition at line 153 of file LightSource_module.cc.

Referenced by beginRun(), LightSource(), and produce().

char evgen::LightSource::fDummyString[256]

private

Definition at line 130 of file LightSource_module.cc.

Referenced by LightSource(), and produce().

Int_t evgen::LightSource::fEvID

private

Definition at line 136 of file LightSource_module.cc.

Referenced by LightSource().

std::string evgen::LightSource::fFileName

private

Definition at line 129 of file LightSource_module.cc.

Referenced by LightSource(), and produce().

bool evgen::LightSource::fFillTree

private

Definition at line 140 of file LightSource_module.cc.

Referenced by LightSource(), and Sample().

int evgen::LightSource::fFirstVoxel

private

Definition at line 180 of file LightSource_module.cc.

Referenced by beginRun(), LightSource(), and produce().

std::ifstream evgen::LightSource::fInputFile

private

Definition at line 128 of file LightSource_module.cc.

Referenced by LightSource(), and produce().

int evgen::LightSource::fLastVoxel

private

Definition at line 181 of file LightSource_module.cc.

Referenced by LightSource(), and produce().

int evgen::LightSource::fN

private

Definition at line 178 of file LightSource_module.cc.

Referenced by LightSource(), produce(), and Sample().

double evgen::LightSource::fP

private

Definition at line 174 of file LightSource_module.cc.

Referenced by LightSource(), produce(), and Sample().

int evgen::LightSource::fPDist

private

Definition at line 143 of file LightSource_module.cc.

Referenced by LightSource(), and Sample().

TTree* evgen::LightSource::fPhotonsGenerated

private

Definition at line 133 of file LightSource_module.cc.

Referenced by LightSource(), and Sample().

bool evgen::LightSource::fPointSource

private

Definition at line 168 of file LightSource_module.cc.

Referenced by LightSource(), and Sample().

int evgen::LightSource::fPosDist

private

Definition at line 141 of file LightSource_module.cc.

Referenced by LightSource(), and Sample().

std::vector<double> evgen::LightSource::fRegionMax

private

Definition at line 160 of file LightSource_module.cc.

Referenced by LightSource().

std::vector<double> evgen::LightSource::fRegionMin

private

Definition at line 159 of file LightSource_module.cc.

Referenced by LightSource().

TLorentzVector evgen::LightSource::fShotMom

private

Definition at line 135 of file LightSource_module.cc.

Referenced by LightSource(), and Sample().

TLorentzVector evgen::LightSource::fShotPos

private

Definition at line 134 of file LightSource_module.cc.

Referenced by LightSource(), and Sample().

double evgen::LightSource::fSigmaP

private

Definition at line 175 of file LightSource_module.cc.

Referenced by LightSource(), produce(), and Sample().

double evgen::LightSource::fSigmaT

private

Definition at line 173 of file LightSource_module.cc.

Referenced by LightSource(), produce(), and Sample().

double evgen::LightSource::fSigmaX

private

Definition at line 170 of file LightSource_module.cc.

Referenced by LightSource(), produce(), and Sample().

double evgen::LightSource::fSigmaY

private

Definition at line 171 of file LightSource_module.cc.

Referenced by LightSource(), produce(), and Sample().

double evgen::LightSource::fSigmaZ

private

Definition at line 172 of file LightSource_module.cc.

Referenced by LightSource(), produce(), and Sample().

int evgen::LightSource::fSourceMode

private

Definition at line 139 of file LightSource_module.cc.

Referenced by LightSource(), and produce().

double evgen::LightSource::fT

private

Definition at line 169 of file LightSource_module.cc.

Referenced by LightSource(), produce(), and Sample().

int evgen::LightSource::fTDist

private

Definition at line 142 of file LightSource_module.cc.

Referenced by LightSource(), and Sample().

sim::PhotonVoxelDef evgen::LightSource::fThePhotonVoxelDef

private

Definition at line 150 of file LightSource_module.cc.

Referenced by LightSource(), and produce().

TVector3 evgen::LightSource::fTPCCenter

private

Definition at line 157 of file LightSource_module.cc.

TVector3 evgen::LightSource::fTPCDimensions

private

Definition at line 158 of file LightSource_module.cc.

bool evgen::LightSource::fUseCustomRegion

private

Definition at line 161 of file LightSource_module.cc.

Referenced by LightSource().

std::string evgen::LightSource::fVersion

private

Definition at line 119 of file LightSource_module.cc.

int evgen::LightSource::fVoxelCount

private

Definition at line 152 of file LightSource_module.cc.

Referenced by LightSource().

double evgen::LightSource::fX

private

Definition at line 165 of file LightSource_module.cc.

Referenced by produce(), and Sample().

int evgen::LightSource::fXSteps

private

Definition at line 146 of file LightSource_module.cc.

Referenced by LightSource().

double evgen::LightSource::fY

private

Definition at line 166 of file LightSource_module.cc.

Referenced by produce(), and Sample().

int evgen::LightSource::fYSteps

private

Definition at line 147 of file LightSource_module.cc.

Referenced by LightSource().

double evgen::LightSource::fZ

private

Definition at line 167 of file LightSource_module.cc.

Referenced by produce(), and Sample().

int evgen::LightSource::fZSteps

private

Definition at line 148 of file LightSource_module.cc.

Referenced by LightSource().

const int evgen::LightSource::kFILE = 0

staticprivate

Definition at line 124 of file LightSource_module.cc.

Referenced by LightSource(), and produce().

const int evgen::LightSource::kGAUS = 1

staticprivate

Definition at line 123 of file LightSource_module.cc.

Referenced by Sample().

const int evgen::LightSource::kSCAN = 1

staticprivate

Definition at line 125 of file LightSource_module.cc.

Referenced by LightSource(), and produce().

const int evgen::LightSource::kUNIF = 0

staticprivate

Definition at line 122 of file LightSource_module.cc.

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

/cvmfs/larsoft.opensciencegrid.org/products/larsim/v07_09_00/source/larsim/EventGenerator/LightSource_module.cc

Public Types

Public Member Functions

Static Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Static Private Attributes

Detailed Description

File mode

Scan mode

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation