sim::SimListUtils Class Reference

#include "SimListUtils.h"

Static Public Member Functions
static sim::LArVoxelList	GetLArVoxelList (const art::Event &evt, std::string moduleLabel)
static sim::SimPhotonsCollection	GetSimPhotonsCollection (const art::Event &evt, std::string moduleLabel)

Detailed Description

Definition at line 22 of file SimListUtils.h.

Member Function Documentation

sim::LArVoxelList sim::SimListUtils::GetLArVoxelList ( const art::Event &  evt, std::string  moduleLabel  )

static

Definition at line 24 of file SimListUtils.cxx.

References sim::LArVoxelList::Add(), sim::LArVoxelList::find(), art::ProductRetriever::getView(), sim::LArG4Parameters::GeVToElectrons(), detinfo::sampling_rate(), detinfo::trigger_offset(), x, y, and z.

Referenced by evdb_tool::DrawLArVoxel3D::Draw(), evd::SimulationDrawer::MCTruth3D(), and evd::SimulationDrawer::MCTruthOrtho().

  {
    art::ServiceHandle<sim::LArG4Parameters const> lgp;
    auto const clocks = art::ServiceHandle<detinfo::DetectorClocksService const>()->DataFor(evt);

    // get the sim::SimChannels
    std::vector<const sim::SimChannel*> sccol;
    evt.getView(moduleLabel, sccol);

    sim::LArVoxelList voxList;

    // loop over the voxels and put them into the list
    for (auto itr = sccol.begin(); itr != sccol.end(); ++itr) {

      // get all sim::IDE associated with this channel
      const auto& idemap = (*itr)->TDCIDEMap();

      // loop over all the sim::IDE values
      for (auto mitr = idemap.begin(); mitr != idemap.end(); mitr++) {

        double time = (*mitr).first - trigger_offset(clocks);
        time *= sampling_rate(clocks);

        // loop over the sim::IDE objects
        const std::vector<sim::IDE>& ide = (*mitr).second;
        for (size_t i = 0; i < ide.size(); ++i) {

          sim::LArVoxelID larVoxelID(ide[i].x, ide[i].y, ide[i].z, time);

          // if energy is unassigned the TrackId is sim::kNoParticleId
          voxList.Add(larVoxelID, ide[i].numElectrons / lgp->GeVToElectrons(), ide[i].trackID);

          // set the voxel id for the just added LArVoxelData
          (*voxList.find(larVoxelID)).second.SetVoxelID(larVoxelID);

        } // end loop over ide for this tdc
      }   // end loop over map
    }     // end loop over sim::SimChannels

    return voxList;
  }

sim::SimPhotonsCollection sim::SimListUtils::GetSimPhotonsCollection ( const art::Event &  evt, std::string  moduleLabel  )

static

get the voxels from the event handle

loop over the pmthits and put them into the list

make an entry in the list for this pmt id

add the photons to the entry

Definition at line 69 of file SimListUtils.cxx.

References art::ProductRetriever::getByLabel(), and sim::SimPhotons::SetChannel().

Referenced by opdet::OptDetDigitizer::produce(), and opdet::OpMCDigi::produce().

  {
    art::Handle<std::vector<sim::SimPhotons>> pmtHandle;
    evt.getByLabel(moduleLabel, pmtHandle);
    const std::vector<sim::SimPhotons>& pmt(*pmtHandle);

    sim::SimPhotonsCollection pmtList;
    pmtList.clear();
    // std::cout << "Building SimPhotonsCollection" << std::endl;

    for (auto itr = pmt.begin(); itr != pmt.end(); ++itr) {

      int ch = (*itr).OpChannel();
      if (pmtList.find(ch) == pmtList.end()) {
        // Create new photon object
        sim::SimPhotons new_photons;
        new_photons.clear();
        new_photons.SetChannel(ch);
        new_photons.reserve((*itr).size());
        pmtList.insert(std::pair<int, sim::SimPhotons>(ch, new_photons));
      }

      for (auto pitr = (*itr).begin(); pitr != (*itr).end(); ++pitr)
        pmtList[ch].push_back(sim::OnePhoton((*pitr)));
    }

    return pmtList;
    // return std::move(pmtList);
  }

---

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

• larsim/v09_43_00/source/larsim/Simulation/SimListUtils.h
• larsim/v09_43_00/source/larsim/Simulation/SimListUtils.cxx