Inheritance diagram for FakeFlash:

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

	FakeFlash (FakeFlash const &)=delete

	FakeFlash (FakeFlash &&)=delete

FakeFlash &	operator= (FakeFlash const &)=delete

FakeFlash &	operator= (FakeFlash &&)=delete

void	produce (art::Event &e) override

void	beginRun (art::Run &run) override

Private Member Functions
void	GenPosition (double &x, double &y, double &z)

void	FillSimPhotons (std::vector< sim::SimPhotons > &simph_v, int nphotons, size_t mother_trackid, const TLorentzVector &pos)

std::vector< double >	GenerateTime (size_t numphotons)

Private Attributes
bool	_verbose
	verbosity for debugging More...

double	_frequency
	[MHz] More...

double	_duration
	[us] More...

double	_tstart
	[ns] More...

size_t	_min_photons
	[photons] More...

size_t	_max_photons
	[photons] More...

std::vector< size_t >	_tpc_v
	List of TPC ID to be used. More...

double	_fast_frac
	fraction of prompt light More...

double	_fast_tau
	scintillation emission time constant for fast component More...

double	_slow_tau
	scintillation emission time constant for slow component More...

size_t	_ch_min
	channel range min to produce SimPhotons More...

size_t	_ch_max
	channel range max to produce SimPhotons More...

double	_xmax
	0.0-1.0 the x-position range in fraction of a TPC volume More...

double	_ymax
	0.0-1.0 the y-position range in fraction of a TPC volume More...

double	_zmax
	0.0-1.0 the z-position range in fraction of a TPC volume More...

double	_xmin
	0.0-1.0 the x-position range in fraction of a TPC volume More...

double	_ymin
	0.0-1.0 the y-position range in fraction of a TPC volume More...

double	_zmin
	0.0-1.0 the z-position range in fraction of a TPC volume More...

CLHEP::HepRandomEngine &	fFlatEngine

CLHEP::RandFlat *	fFlatRandom

CLHEP::RandExponential *	fExpoRandom

CLHEP::RandPoisson *	fPoisRandom

Detailed Description

Definition at line 40 of file FakeFlash_module.cc.

Constructor & Destructor Documentation

FakeFlash::FakeFlash ( fhicl::ParameterSet const & p )

explicit

Definition at line 89 of file FakeFlash_module.cc.

   : EDProducer{p}
   , fFlatEngine(art::ServiceHandle<rndm::NuRandomService>()->createEngine(*this, "HepJamesRandom", "Gen", p, "Seed"))
   // More initializers here.
 {
   _verbose = p.get<bool>("Verbose",false); // If you want someone to talk to you
   auto min_photons = p.get<int>("MinPhotons",24000);     // Min of the range of photons to be injected in one shot
   auto max_photons = p.get<int>("MaxPhotons",2400000);   // Max of the range of photons to be injected in one shot
   assert(min_photons < max_photons && min_photons>0 && max_photons>0);
   _min_photons = min_photons;
   _max_photons = max_photons;
 
   // scintillation params (yeah fine you can replace these w/ service)
   _fast_frac = p.get<double>("PromptLightFraction",0.23);
   _fast_tau  = p.get<double>("FastTimeConstant",0.006);
   _slow_tau  = p.get<double>("SlowTimeConstant",1.5);
 
   // channel range to simulate
   auto geop = lar::providerFrom<geo::Geometry>();
   _ch_min = 0;
   _ch_max = geop->NOpChannels() - 1;
   _ch_min = p.get<size_t>("ChannelMin",_ch_min);
   _ch_max = p.get<size_t>("ChannelMax",_ch_max);
   assert(_ch_min<_ch_max);
   // Given in micro-seconds as larsoft default time unit, but then we convert to ns to record as photon time
   _fast_tau *= 1.e3;
   _slow_tau *= 1.e3;
   // TPC list + range
   _tpc_v = p.get<std::vector<size_t> >("TPCList"); // The list of TPC in which flash can happen
   _xmax  = p.get<double>("XMax",1.0);
   _xmin  = p.get<double>("XMin",0.0);
   _ymax  = p.get<double>("YMax",1.0);
   _ymin  = p.get<double>("YMin",0.0);
   _zmax  = p.get<double>("ZMax",1.0);
   _zmin  = p.get<double>("ZMin",0.0);
   assert(_xmax>=_xmin && _ymax>=_ymin && _zmax>=_zmin &&
          _xmax<=1.0 && _ymax<=1.0 && _zmax<=1.0 &&
          _xmin>=0.0 && _ymin>=0.0 && _zmin>=0.0 );
   // generation timing info
   _frequency = p.get<double>("Frequency"); // The frequency of photon(s) injection
   _duration  = p.get<double>("Duration");  // The duration of photon(s) injection period
   _tstart    = p.get<double>("G4TStart");  // The start time of photon injection period
   produces<std::vector<sim::SimPhotons> >();
   produces<std::vector<simb::MCTruth> >();
   produces< sumdata::RunData, art::InRun >();
 
   fFlatRandom = new CLHEP::RandFlat(fFlatEngine,0,1);
   fExpoRandom = new CLHEP::RandExponential(fFlatEngine);
   fPoisRandom = new CLHEP::RandPoisson(fFlatEngine);
 }

FakeFlash::FakeFlash ( FakeFlash const & )

delete

FakeFlash::FakeFlash ( FakeFlash && )

delete

Member Function Documentation

void FakeFlash::beginRun ( art::Run & run )

override

Definition at line 140 of file FakeFlash_module.cc.

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

void FakeFlash::FillSimPhotons	(	std::vector< sim::SimPhotons > &	simph_v,
		int	nphotons,
		size_t	mother_trackid,
		const TLorentzVector &	pos
	)

private

Definition at line 178 of file FakeFlash_module.cc.

 {
   art::ServiceHandle<phot::PhotonVisibilityService const> pvs;
   double xyz[3];
   xyz[0] = pos.X();
   xyz[1] = pos.Y();
   xyz[2] = pos.Z();
   auto const& Visibilities = pvs->GetAllVisibilities(xyz);
   // Now loop over optical channels and fill SimPhotons array
   if(simph_v.empty()) {
     for(size_t opch=_ch_min; opch<=_ch_max; ++opch)
       simph_v.push_back(sim::SimPhotons(opch));
   }
   for(size_t opch=_ch_min; opch <= _ch_max; ++opch) {
     // get visibility
     size_t detected = fPoisRandom->fire((double)(nphotons) * Visibilities[opch]);
     //std::cout<<opch<<","<<detected<<std::endl;
     auto time_array = this->GenerateTime(detected);
     assert(time_array.size() == detected);
     // record
     auto& simph = simph_v[opch];
     for(size_t idx=0; idx<detected; ++idx) {
         sim::OnePhoton phot;
         phot.Time = pos.T() + time_array[idx];
         phot.InitialPosition.SetX(pos.X());
         phot.InitialPosition.SetY(pos.Y());
         phot.InitialPosition.SetZ(pos.Z());
         phot.MotherTrackID = mother_trackid;
         simph.emplace_back(phot);
     }
   }
 }

std::vector< double > FakeFlash::GenerateTime ( size_t numphotons )

private

Definition at line 214 of file FakeFlash_module.cc.

                                                            {
 
   std::vector<double> res(numphotons);
   for(int i=0; i<((int)numphotons); ++i) {
     if(fFlatRandom->fire(0.,1.) < _fast_frac)
       res[i] = fExpoRandom->fire(_fast_tau);
     else res[i] = fExpoRandom->fire(_slow_tau);
   }
   return res;
 }

void FakeFlash::GenPosition	(	double &	x,
		double &	y,
		double &	z
	)

private

Definition at line 152 of file FakeFlash_module.cc.

                                                            {
 
     size_t tpc_id = (size_t)(fFlatRandom->fire(0,_tpc_v.size()));
     bool found = false;
     // Implementation of required member function here.
     auto geop = lar::providerFrom<geo::Geometry>();
     for(size_t c=0; c<geop->Ncryostats(); ++c) {
         auto const& cryostat = geop->Cryostat(c);
         if(!cryostat.HasTPC(tpc_id)) continue;
         auto const& tpc = cryostat.TPC(tpc_id);
         auto const& tpcabox = tpc.ActiveBoundingBox();
         double xmin = tpcabox.MinX() + (tpcabox.MaxX() - tpcabox.MinX()) * _xmin;
         double xmax = tpcabox.MinX() + (tpcabox.MaxX() - tpcabox.MinX()) * _xmax;
         double ymin = tpcabox.MinY() + (tpcabox.MaxY() - tpcabox.MinY()) * _ymin;
         double ymax = tpcabox.MinY() + (tpcabox.MaxY() - tpcabox.MinY()) * _ymax;
         double zmin = tpcabox.MinZ() + (tpcabox.MaxZ() - tpcabox.MinZ()) * _zmin;
         double zmax = tpcabox.MinZ() + (tpcabox.MaxZ() - tpcabox.MinZ()) * _zmax;
         x = fFlatRandom->fire(xmin,xmax);
         y = fFlatRandom->fire(ymin,ymax);
         z = fFlatRandom->fire(zmin,zmax);
         found = true;
         break;
     }
     if(!found) std::cerr<< "\033[93mTPC " << tpc_id << " not found...\033[00m" << std::endl;
 }

FakeFlash& FakeFlash::operator= ( FakeFlash const & )

delete

FakeFlash& FakeFlash::operator= ( FakeFlash && )

delete

void FakeFlash::produce ( art::Event & e )

override

Definition at line 225 of file FakeFlash_module.cc.

 {
 
   auto simph_v = std::unique_ptr<std::vector<sim::SimPhotons> >(new std::vector<sim::SimPhotons>());
   auto mct_v   = std::unique_ptr<std::vector<simb::MCTruth> >(new std::vector<simb::MCTruth>());
   std::vector<simb::MCParticle> mcp_v;
   double clock = 0.;
   while(clock <= _duration) {
     // Determine photon count
     int nphotons = _min_photons + (_max_photons - _min_photons) * (fFlatRandom->fire(0,1));
 
     // Generate position 4 vector
     double time = _tstart + clock * 1.e3;
     double x,y,z;
     this->GenPosition(x,y,z);
     TLorentzVector pos(x,y,z,time);
 
     // Fill SimPhotons
     this->FillSimPhotons(*simph_v, (int)(mcp_v.size()), nphotons, pos);
 
     // Record a fake photon particle
     simb::MCParticle part(mcp_v.size(), 22, "primary", 0, 0., 1);
     TLorentzVector mom(0.,0.,0.,(double)nphotons);
     part.AddTrajectoryPoint(pos,mom);
     mcp_v.emplace_back(std::move(part));
 
     // count time
     clock += (1./_frequency);
   }
 
   simb::MCTruth mct;
   for(auto& part : mcp_v)
     mct.Add(part);
   mct_v->emplace_back(std::move(mct));
 
   e.put(std::move(simph_v));
   e.put(std::move(mct_v));
 }

Member Data Documentation

size_t FakeFlash::_ch_max

private

channel range max to produce SimPhotons

Definition at line 75 of file FakeFlash_module.cc.

size_t FakeFlash::_ch_min

private

channel range min to produce SimPhotons

Definition at line 74 of file FakeFlash_module.cc.

double FakeFlash::_duration

private

[us]

Definition at line 66 of file FakeFlash_module.cc.

double FakeFlash::_fast_frac

private

fraction of prompt light

Definition at line 71 of file FakeFlash_module.cc.

double FakeFlash::_fast_tau

private

scintillation emission time constant for fast component

Definition at line 72 of file FakeFlash_module.cc.

double FakeFlash::_frequency

private

[MHz]

Definition at line 65 of file FakeFlash_module.cc.

size_t FakeFlash::_max_photons

private

[photons]

Definition at line 69 of file FakeFlash_module.cc.

size_t FakeFlash::_min_photons

private

[photons]

Definition at line 68 of file FakeFlash_module.cc.

double FakeFlash::_slow_tau

private

scintillation emission time constant for slow component

Definition at line 73 of file FakeFlash_module.cc.

std::vector<size_t> FakeFlash::_tpc_v

private

List of TPC ID to be used.

Definition at line 70 of file FakeFlash_module.cc.

double FakeFlash::_tstart

private

[ns]

Definition at line 67 of file FakeFlash_module.cc.

bool FakeFlash::_verbose

private

verbosity for debugging

Definition at line 64 of file FakeFlash_module.cc.

double FakeFlash::_xmax

private

0.0-1.0 the x-position range in fraction of a TPC volume

Definition at line 76 of file FakeFlash_module.cc.

double FakeFlash::_xmin

private

0.0-1.0 the x-position range in fraction of a TPC volume

Definition at line 79 of file FakeFlash_module.cc.

double FakeFlash::_ymax

private

0.0-1.0 the y-position range in fraction of a TPC volume

Definition at line 77 of file FakeFlash_module.cc.

double FakeFlash::_ymin

private

0.0-1.0 the y-position range in fraction of a TPC volume

Definition at line 80 of file FakeFlash_module.cc.

double FakeFlash::_zmax

private

0.0-1.0 the z-position range in fraction of a TPC volume

Definition at line 78 of file FakeFlash_module.cc.

double FakeFlash::_zmin

private

0.0-1.0 the z-position range in fraction of a TPC volume

Definition at line 81 of file FakeFlash_module.cc.

CLHEP::RandExponential* FakeFlash::fExpoRandom

private

Definition at line 84 of file FakeFlash_module.cc.

CLHEP::HepRandomEngine& FakeFlash::fFlatEngine

private

Definition at line 82 of file FakeFlash_module.cc.

CLHEP::RandFlat* FakeFlash::fFlatRandom

private

Definition at line 83 of file FakeFlash_module.cc.

CLHEP::RandPoisson* FakeFlash::fPoisRandom

private

Definition at line 85 of file FakeFlash_module.cc.

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

FakeFlash_module.cc

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation