sim::SimListUtils Class Reference

#include "SimListUtils.h"

Public Member Functions
	SimListUtils ()
virtual	~SimListUtils ()

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 21 of file SimListUtils.h.

Constructor & Destructor Documentation

sim::SimListUtils::SimListUtils

(

)

Definition at line 23 of file SimListUtils.cxx.

  {
  }

sim::SimListUtils::~SimListUtils ( )

virtual

Definition at line 28 of file SimListUtils.cxx.

  {
  }

Member Function Documentation

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

static

Definition at line 35 of file SimListUtils.cxx.

References sim::LArVoxelList::Add(), sim::LArVoxelList::find(), art::DataViewImpl::getView(), sim::LArG4Parameters::GeVToElectrons(), detinfo::DetectorProperties::SamplingRate(), detinfo::DetectorProperties::TriggerOffset(), x, y, and z.

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

  {
    art::ServiceHandle<sim::LArG4Parameters> lgp;
    const detinfo::DetectorProperties* detprop = lar::providerFrom<detinfo::DetectorPropertiesService>();

    // 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();
      //std::map<unsigned short, std::vector<sim::IDE> >::const_iterator mitr;

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

        double time = (*mitr).first - detprop->TriggerOffset();
        time *= detprop->SamplingRate();
        
        // 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 84 of file SimListUtils.cxx.

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

Referenced by opdet::SimPhotonCounter::analyze(), opdet::OpMCDigi::produce(), and opdet::OptDetDigitizer::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/v06_53_00/source/larsim/Simulation/SimListUtils.h
• larsim/v06_53_00/source/larsim/Simulation/SimListUtils.cxx