sbnd::crt::CRTDetSim Class Reference

#include <CRTDetSim.h>

Inheritance diagram for sbnd::crt::CRTDetSim:

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

Public Attributes
std::string	fG4ModuleLabel

Private Attributes
CLHEP::HepRandomEngine &	fEngine
	Reference to art-managed random-number engine. More...
double	fG4RefTime
	Stores the G4 reference time. More...
CRTDetSimAlg	fDetAlg
	Instance of the CRT detector simulation algorithm. More...

Detailed Description

Definition at line 31 of file CRTDetSim.h.

Constructor & Destructor Documentation

sbnd::crt::CRTDetSim::CRTDetSim ( fhicl::ParameterSet const &  p )

explicit

Definition at line 53 of file sbndcode/sbndcode/CRT/CRTDetSim_module.cc.

  : EDProducer{p}
  , fEngine(art::ServiceHandle<rndm::NuRandomService>{}->createEngine(*this, "HepJamesRandom", "crt", p, "Seed"))

sbnd::crt::CRTDetSim::CRTDetSim ( CRTDetSim const &  )

delete

sbnd::crt::CRTDetSim::CRTDetSim ( CRTDetSim &&  )

delete

Member Function Documentation

CRTDetSim& sbnd::crt::CRTDetSim::operator= ( CRTDetSim const &  )

delete

CRTDetSim& sbnd::crt::CRTDetSim::operator= ( CRTDetSim &&  )

delete

void sbnd::crt::CRTDetSim::produce ( art::Event &  e )

override

Definition at line 70 of file sbndcode/sbndcode/CRT/CRTDetSim_module.cc.

                                    {

  std::unique_ptr<std::vector<sbnd::crt::FEBData> > FEBDataOut(new std::vector<sbnd::crt::FEBData>);
  art::PtrMaker<sbnd::crt::FEBData> makeDataPtr(e);

  std::unique_ptr<std::vector<sim::AuxDetIDE> > auxDetIdes(new std::vector<sim::AuxDetIDE>);
  art::PtrMaker<sim::AuxDetIDE> makeIdePtr(e);

  std::unique_ptr<art::Assns<sbnd::crt::FEBData, sim::AuxDetIDE, sbnd::crt::FEBTruthInfo>> Dataassn
    (new art::Assns<sbnd::crt::FEBData, sim::AuxDetIDE, sbnd::crt::FEBTruthInfo>);

  // Handle for (truth) AuxDetSimChannels
  art::Handle<std::vector<sim::AuxDetSimChannel> > channels;
  e.getByLabel(fG4ModuleLabel, channels);

  if (!channels.isValid()) {
    mf::LogWarning("CRTDetSim") << "No AuxDetSimChannel..." << std::endl;
  }

  fDetAlg.ClearTaggers();


  //
  // Step 1: Construct Taggers
  //
  for(auto const& adsc : *channels) {

    if(adsc.AuxDetID() == UINT_MAX || adsc.AuxDetSensitiveID() == UINT_MAX) {
      mf::LogWarning("CRTDetSim") << "AuxDetSimChannel with ID: UINT_MAX\n"
                                  << "skipping channel...";
      continue;
    }

    if(adsc.AuxDetIDEs().size() > 0) {
      mf::LogDebug("CRTDetSim") << "Adding AuxDetSimChannel with AuxDetID " << adsc.AuxDetID()
                                << ", AuxDetSensitiveID " << adsc.AuxDetSensitiveID() << std::endl;
      fDetAlg.FillTaggers(adsc.AuxDetID(), adsc.AuxDetSensitiveID(), adsc.AuxDetIDEs());
    }

  } // loop over AuxDetSimChannels


  //
  // Step 1: Apply Coincidence, deadtime, etc.
  //

  fDetAlg.CreateData();
  std::vector<std::pair<sbnd::crt::FEBData, std::vector<sim::AuxDetIDE>>> data = fDetAlg.GetData();
  std::vector<std::vector<int>> auxdata = fDetAlg.GetAuxData();

  //
  // Step 3: Save output
  //

  // for(auto const& dataPair : data){
  for (size_t i = 0; i < data.size(); i++) {

    auto & dataPair = data[i];
    auto & feb = dataPair.first;
    auto & ides = dataPair.second;
    auto & idxs = auxdata.at(i);

    FEBDataOut->push_back(feb);
    art::Ptr<FEBData> dataPtr = makeDataPtr(FEBDataOut->size()-1);

    for(size_t j = 0; j < ides.size(); j++){
      auto const& ide = ides[j];
      auxDetIdes->push_back(ide);
      FEBTruthInfo const feb_truth_info {idxs[j]};
      mf::LogDebug("CRTDetSim") << "Adding " << idxs[j] << " to FEBTruthInfo." << std::endl;
      art::Ptr<sim::AuxDetIDE> idePtr = makeIdePtr(auxDetIdes->size()-1);
      Dataassn->addSingle(dataPtr, idePtr, feb_truth_info);
    }
  }

  mf::LogInfo("CRTDetSim") << "Number of FEBData objects produced: " << FEBDataOut->size() << "\n";

  e.put(std::move(FEBDataOut));
  e.put(std::move(auxDetIdes));
  e.put(std::move(Dataassn));
}

void sbnd::crt::CRTDetSim::reconfigure

(

fhicl::ParameterSet const &

p

)

Definition at line 48 of file sbndcode/sbndcode/CRT/CRTDetSim_module.cc.

                                                       {
  fG4ModuleLabel = p.get<std::string>("G4ModuleLabel");
}

Member Data Documentation

CRTDetSimAlg sbnd::crt::CRTDetSim::fDetAlg

private

Instance of the CRT detector simulation algorithm.

Definition at line 48 of file CRTDetSim.h.

CLHEP::HepRandomEngine& sbnd::crt::CRTDetSim::fEngine

private

Reference to art-managed random-number engine.

Definition at line 46 of file CRTDetSim.h.

std::string sbnd::crt::CRTDetSim::fG4ModuleLabel

Definition at line 42 of file CRTDetSim.h.

double sbnd::crt::CRTDetSim::fG4RefTime

private

Stores the G4 reference time.

Definition at line 47 of file CRTDetSim.h.

---

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

• CRTDetSim.h
• sbndcode/sbndcode/CRT/CRTDetSim_module.cc