icarus::crt::CRTDetSim Class Reference

Inheritance diagram for icarus::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
string	fG4ModuleLabel

Private Attributes
CLHEP::HepRandomEngine &	fRandEngine
CRTDetSimAlg	detAlg

Detailed Description

Definition at line 109 of file icaruscode/icaruscode/CRT/CRTDetSim_module.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 136 of file icaruscode/icaruscode/CRT/CRTDetSim_module.cc.

                                                         : EDProducer{p},
    fRandEngine(art::ServiceHandle<rndm::NuRandomService>()->createEngine(*this, "HepJamesRandom", "crt", p, "Seed")),
    detAlg(p.get<fhicl::ParameterSet>("DetSimAlg"),fRandEngine)
{
    this->reconfigure(p);

    produces< vector<CRTData> >();
    produces< vector<sim::AuxDetIDE> >();
    produces< art::Assns<CRTData, sim::AuxDetIDE> >();
}

icarus::crt::CRTDetSim::CRTDetSim ( CRTDetSim const &  )

delete

icarus::crt::CRTDetSim::CRTDetSim ( CRTDetSim &&  )

delete

Member Function Documentation

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

delete

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

delete

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

override

Definition at line 148 of file icaruscode/icaruscode/CRT/CRTDetSim_module.cc.

                                                {

    //pointer to vector of CRT data products to be pushed to event
    std::unique_ptr< vector<CRTData> > triggeredCRTHits (
        new vector<CRTData>);
    //pointer to vector of AuxDetIDE products to be pushed to event
    std::unique_ptr< vector<sim::AuxDetIDE> > ides (
        new vector<sim::AuxDetIDE>);
    //pointer associations between CRT data products and AuxDetIDEs to be pushed to event
    std::unique_ptr< art::Assns<CRTData,sim::AuxDetIDE> > dataAssn (
        new art::Assns<CRTData,sim::AuxDetIDE> );
    art::PtrMaker<CRTData> makeDataPtr(event); 
    art::PtrMaker<sim::AuxDetIDE> makeIDEPtr(event); 

    //clear detAlg member data
    detAlg.ClearTaggers();

    // Handle for (truth) AuxDetSimChannels
    art::Handle< vector<sim::AuxDetSimChannel> > adChanHandle;;
    vector< art::Ptr<sim::AuxDetSimChannel> > adChanList;
    if (event.getByLabel(fG4ModuleLabel, adChanHandle) )
        art::fill_ptr_vector(adChanList, adChanHandle);

    mf::LogInfo("CRTDetSimProducer")
        <<"Number of AuxDetChannels = " << adChanList.size();

    int nide=0;
    auto const clockData = art::ServiceHandle<detinfo::DetectorClocksService>()->DataFor(event);
    for(auto const& adsc : adChanList) {

        auto const& auxDetIDEs = adsc->AuxDetIDEs();
        if(auxDetIDEs.size()>0) {
            nide++;
            detAlg.FillTaggers(clockData, adsc->AuxDetID(), adsc->AuxDetSensitiveID(), auxDetIDEs);
        }

    }//loop over AuxDetSimChannels

    //generate CRTData products, associates from filled detAlg member data
    int nData = 0;
    if(nide>0) {
        vector<std::pair<CRTData,vector<sim::AuxDetIDE>>> data = detAlg.CreateData();
        //time sort CRT data (ascending)
        std::sort(data.begin(),data.end(),
            [](const std::pair<CRTData,vector<sim::AuxDetIDE>>& d1, 
               const std::pair<CRTData,vector<sim::AuxDetIDE>>& d2) {
                return d1.first.fTs0 < d2.first.fTs0;
            });

        for(auto const& dataPair : data){
        
          triggeredCRTHits->push_back(dataPair.first);
          art::Ptr<CRTData> dataPtr = makeDataPtr(triggeredCRTHits->size()-1);
          nData++;

          for(auto const& ide : dataPair.second){
            ides->push_back(ide);
            art::Ptr<sim::AuxDetIDE> idePtr = makeIDEPtr(ides->size()-1);
            dataAssn->addSingle(dataPtr, idePtr);
          }
        }

    }

    if(nData==0)
        mf::LogWarning("CRTDetSimProducer")
          << "0 CRTData produced (expected for most neutrino events, never for cosmics)";

    //push products to event
    event.put(std::move(triggeredCRTHits));
    event.put(std::move(ides));
    event.put(std::move(dataAssn));

    mf::LogInfo("CRTDetSimProducer")
      << "Number of CRT data produced = "<< nData << '\n'
      << "  from " << nide << " AuxDetIDEs";
}

void icarus::crt::CRTDetSim::reconfigure

(

fhicl::ParameterSet const &

p

)

Definition at line 131 of file icaruscode/icaruscode/CRT/CRTDetSim_module.cc.

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

Member Data Documentation

CRTDetSimAlg icarus::crt::CRTDetSim::detAlg

private

Definition at line 126 of file icaruscode/icaruscode/CRT/CRTDetSim_module.cc.

string icarus::crt::CRTDetSim::fG4ModuleLabel

Definition at line 121 of file icaruscode/icaruscode/CRT/CRTDetSim_module.cc.

CLHEP::HepRandomEngine& icarus::crt::CRTDetSim::fRandEngine

private

Definition at line 125 of file icaruscode/icaruscode/CRT/CRTDetSim_module.cc.

---

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

• icaruscode/icaruscode/CRT/CRTDetSim_module.cc