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 25 of file SimListUtils.cxx.

  {
    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 71 of file SimListUtils.cxx.

  {
    /// get the voxels from the event handle
    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;

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

      int ch = (*itr).OpChannel();
      /// make an entry in the list for this pmt id
      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));
      }

      /// add the photons to the entry
      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:

• SimListUtils.h
• SimListUtils.cxx