ICARUSMCOpHit Class Reference

Inheritance diagram for ICARUSMCOpHit:

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

Private Attributes
double	_merge_period
std::string	_simph_producer
double	_spe_area
double	_spe_amp

Detailed Description

Definition at line 30 of file ICARUSMCOpHit_module.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 55 of file ICARUSMCOpHit_module.cc.

  : EDProducer{p}
{
  _merge_period = p.get<double>("MergePeriod");
  _simph_producer = p.get<std::string>("SimPhotonsProducer");
  _spe_area = p.get<double>("SPEArea");
  _spe_amp  = p.get<double>("SPEAmplitude");
  produces<std::vector<recob::OpHit> >();
}

ICARUSMCOpHit::ICARUSMCOpHit ( ICARUSMCOpHit const &  )

delete

ICARUSMCOpHit::ICARUSMCOpHit ( ICARUSMCOpHit &&  )

delete

Member Function Documentation

ICARUSMCOpHit& ICARUSMCOpHit::operator= ( ICARUSMCOpHit const &  )

delete

ICARUSMCOpHit& ICARUSMCOpHit::operator= ( ICARUSMCOpHit &&  )

delete

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

override

Definition at line 65 of file ICARUSMCOpHit_module.cc.

{
  auto oph_v = std::unique_ptr<std::vector<recob::OpHit> >(new std::vector<recob::OpHit>());

  // Get SimPhotons.

  art::Handle< std::vector< sim::SimPhotons > > simph_h;
  e.getByLabel(_simph_producer,simph_h);

  // Get SimPhotonsLite.

  art::Handle< std::vector< sim::SimPhotonsLite > > simphlite_h;
  e.getByLabel(_simph_producer,simphlite_h);

  // At least one, but not both, of the handles should be valid.

  if(!simph_h.isValid() && !simphlite_h.isValid()) {
    std::cerr << "Could not retrieve sim::SimPhotons or sim::SimPhotonsLite from producer label: " << _simph_producer << std::endl;
    throw std::exception();
  }
  if(simph_h.isValid() && simphlite_h.isValid()) {
    std::cerr << "Found both sim::SimPhotons and sim::SimPhotonsLite from producer label: " << _simph_producer << std::endl;
    throw std::exception();
  }

  // Combine the two handles into a single vector so that they can be processed in one loop.

  typedef std::variant<const sim::SimPhotons*, const sim::SimPhotonsLite*> EitherSimPhoton;
  std::vector<EitherSimPhoton> sim_photons;
  if(simph_h.isValid()) {
    sim_photons.reserve(simph_h->size());
    for(auto const& simph : *simph_h) {
      sim_photons.push_back(&simph);
    }
  }
  else if(simphlite_h.isValid()) {
    sim_photons.reserve(simphlite_h->size());
    for(auto const& simphlite : *simphlite_h) {
      sim_photons.push_back(&simphlite);
    }
  }

  // Loop over SimPhotons[Lite]

  std::vector<bool> processed_v;
  auto const clockData = art::ServiceHandle<detinfo::DetectorClocksService const>()->DataFor(e);
  for(auto const& simph : sim_photons) {
    // Make sure channel number is unique (e.g. one sim::SimPhotons per op channel)
    size_t opch = -1;
    if(simph_h.isValid())
      opch = std::get<const sim::SimPhotons*>(simph)->OpChannel();
    else if(simphlite_h.isValid())
      opch = std::get<const sim::SimPhotonsLite*>(simph)->OpChannel;

    if(opch >= processed_v.size()) processed_v.resize(opch+1,false);
    if(processed_v[opch]) {
      std::cerr << "Found duplicate channels in std::vector<sim::SimPhotons>! not expected (logic will fail).."<<std::endl;
      throw std::exception();
    }

    processed_v[opch] = true;
    bool in_window  = false;
    double oph_time = -1.e9;
    double pe = 0.;
    // Insert photon times into a sorted set
    std::map<double,size_t> time_m;
    if(simph_h.isValid()) {
      for(auto const& oneph : *std::get<const sim::SimPhotons*>(simph)) {
        double this_time = clockData.G4ToElecTime(oneph.Time) - clockData.TriggerTime();
        time_m[this_time] += 1;
      }
    }
    else if(simphlite_h.isValid()) {
      for(auto const& [ time_ns, nphotons ] : std::get<const sim::SimPhotonsLite*>(simph)->DetectedPhotons) {
        double this_time = clockData.G4ToElecTime(time_ns + 0.5) - clockData.TriggerTime();
        time_m[this_time] += nphotons;
      }
    }

    // Loop over the time vector, emplace photons
    for(auto const& time_photon_pair : time_m) {

      auto const& this_time = time_photon_pair.first;

      if(this_time > (oph_time + _merge_period) && in_window) {
        recob::OpHit oph(opch, 
                         oph_time,
                         oph_time + clockData.TriggerTime(),
                         0, // frame
                         1., // width
                         pe * _spe_area, // area,
                         pe * _spe_amp, // peakheight,
                         pe, // pe
                         0.);
        oph_v->emplace_back(std::move(oph));
        in_window = false;
        pe = 0;
      }

      if(!in_window) oph_time = this_time;
      in_window = true;
      pe += time_photon_pair.second;
    }
    if(in_window) {
      recob::OpHit oph(opch, 
                       oph_time,
                       oph_time + clockData.TriggerTime(),
                       0, // frame
                       1., // width
                       pe * _spe_area, // area,
                       pe * _spe_amp, // peakheight,
                       pe, // pe
                       0.);
      oph_v->emplace_back(std::move(oph));
    }
  }

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

Member Data Documentation

double ICARUSMCOpHit::_merge_period

private

Definition at line 49 of file ICARUSMCOpHit_module.cc.

std::string ICARUSMCOpHit::_simph_producer

private

Definition at line 50 of file ICARUSMCOpHit_module.cc.

double ICARUSMCOpHit::_spe_amp

private

Definition at line 51 of file ICARUSMCOpHit_module.cc.

double ICARUSMCOpHit::_spe_area

private

Definition at line 51 of file ICARUSMCOpHit_module.cc.

---

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

• ICARUSMCOpHit_module.cc