FakePhotoS Class Reference

Inheritance diagram for FakePhotoS:

Public Member Functions
	FakePhotoS (fhicl::ParameterSet const &p)
	FakePhotoS (FakePhotoS const &)=delete
	FakePhotoS (FakePhotoS &&)=delete
FakePhotoS &	operator= (FakePhotoS const &)=delete
FakePhotoS &	operator= (FakePhotoS &&)=delete
void	produce (art::Event &e) override
void	beginRun (art::Run &run) override

Private Member Functions
std::vector< double >	GenerateTime (size_t numphotons)

Private Attributes
bool	_verbose
double	_frequency
	[MHz] More...
double	_duration
	[us] More...
double	_tstart
	[ns] More...
size_t	_min_pe
	[p.e.] More...
size_t	_max_pe
	[p.e.] More...
bool	_simulate_time
	whether or not to simulate scintillation fall time 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...
std::vector< unsigned int >	_ch_v
CLHEP::HepRandomEngine &	fFlatEngine
CLHEP::RandFlat *	fFlatRandom
CLHEP::RandExponential *	fExpoRandom

Detailed Description

Definition at line 35 of file FakePhotoS_module.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 72 of file FakePhotoS_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_pe = p.get<int>("MinPE",1);     // Min of the range of PE count to be injected in one shot
  auto max_pe = p.get<int>("MaxPE",1);     // Max of the range of PE count to be injected in one shot
  assert(min_pe <= max_pe && min_pe>0 && max_pe>0); 
  _min_pe = min_pe;
  _max_pe = max_pe;
  _ch_v.clear();
  _ch_v = p.get<std::vector<unsigned int> >("Channels",_ch_v); // Specify if you want to use only select channels
  _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< sumdata::RunData, art::InRun >();
  // scintillation params (yeah fine you can replace these w/ service)
  _simulate_time = p.get<bool>("SimulateTime",false); // scintillation fall timing simulation
  _fast_frac = p.get<double>("PromptLightFraction",0.23);
  _fast_tau  = p.get<double>("FastTimeConstant",0.006);
  _slow_tau  = p.get<double>("SlowTimeConstant",1.5);

  fFlatRandom = new CLHEP::RandFlat(fFlatEngine,_min_pe,_max_pe);
  fExpoRandom = new CLHEP::RandExponential(fFlatEngine);
}

FakePhotoS::FakePhotoS ( FakePhotoS const &  )

delete

FakePhotoS::FakePhotoS ( FakePhotoS &&  )

delete

Member Function Documentation

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

override

Definition at line 100 of file FakePhotoS_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;
}

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

private

Definition at line 113 of file FakePhotoS_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;
}

FakePhotoS& FakePhotoS::operator= ( FakePhotoS const &  )

delete

FakePhotoS& FakePhotoS::operator= ( FakePhotoS &&  )

delete

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

override

Definition at line 124 of file FakePhotoS_module.cc.

{
  // Implementation of required member function here.
  auto simph_v = std::unique_ptr<std::vector<sim::SimPhotons> >(new std::vector<sim::SimPhotons>());

  if(_ch_v.empty()) {
    auto const geop = lar::providerFrom<geo::Geometry>();
    _ch_v.reserve(geop->NOpChannels());
    for(size_t opch=0; opch<geop->NOpChannels(); ++opch) 
      _ch_v.push_back(opch);
  }
  simph_v->reserve(_ch_v.size());
  for(auto const& opch : _ch_v) {
    if(_verbose) std::cout << "OpChannel " << opch << std::endl;
    sim::SimPhotons sphot(opch);
    sphot.reserve(int(_duration * _frequency));
    double clock = 0.;
    while(clock <= _duration) {
      
      size_t npe = fFlatRandom->fire(_min_pe,_max_pe+0.999999);
      if(_verbose) std::cout << npe << "@" << (int)(_tstart + clock * 1.e3) << "[ns] " << std::flush; 
      std::vector<double> time_array(npe,0.);
      if(_simulate_time) time_array = this->GenerateTime(npe);
      for(size_t ctr=0;ctr<npe;++ctr) {
        sim::OnePhoton phot;
        phot.Time = _tstart + clock * 1.e3 + time_array[ctr];
        sphot.emplace_back(std::move(phot));
      }
      clock += (1./_frequency);
    }
    if(_verbose) std::cout << std::endl;
    simph_v->emplace_back(std::move(sphot));
  }

  e.put(std::move(simph_v));
}

Member Data Documentation

std::vector<unsigned int> FakePhotoS::_ch_v

private

Definition at line 64 of file FakePhotoS_module.cc.

double FakePhotoS::_duration

private

[us]

Definition at line 56 of file FakePhotoS_module.cc.

double FakePhotoS::_fast_frac

private

fraction of prompt light

Definition at line 61 of file FakePhotoS_module.cc.

double FakePhotoS::_fast_tau

private

scintillation emission time constant for fast component

Definition at line 62 of file FakePhotoS_module.cc.

double FakePhotoS::_frequency

private

[MHz]

Definition at line 55 of file FakePhotoS_module.cc.

size_t FakePhotoS::_max_pe

private

[p.e.]

Definition at line 59 of file FakePhotoS_module.cc.

size_t FakePhotoS::_min_pe

private

[p.e.]

Definition at line 58 of file FakePhotoS_module.cc.

bool FakePhotoS::_simulate_time

private

whether or not to simulate scintillation fall time

Definition at line 60 of file FakePhotoS_module.cc.

double FakePhotoS::_slow_tau

private

scintillation emission time constant for slow component

Definition at line 63 of file FakePhotoS_module.cc.

double FakePhotoS::_tstart

private

[ns]

Definition at line 57 of file FakePhotoS_module.cc.

bool FakePhotoS::_verbose

private

Definition at line 54 of file FakePhotoS_module.cc.

CLHEP::RandExponential* FakePhotoS::fExpoRandom

private

Definition at line 68 of file FakePhotoS_module.cc.

CLHEP::HepRandomEngine& FakePhotoS::fFlatEngine

private

Definition at line 66 of file FakePhotoS_module.cc.

CLHEP::RandFlat* FakePhotoS::fFlatRandom

private

Definition at line 67 of file FakePhotoS_module.cc.

---

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

• FakePhotoS_module.cc