sim::MergeSimSourcesUtility Class Reference

#include "MergeSimSources.h"

Public Member Functions
	MergeSimSourcesUtility (std::vector< int > const &)
void	Reset ()
void	MergeMCParticles (std::vector< simb::MCParticle > &, const std::vector< simb::MCParticle > &, size_t)
void	MergeSimChannels (std::vector< sim::SimChannel > &, const std::vector< sim::SimChannel > &, size_t)
void	MergeAuxDetSimChannels (std::vector< sim::AuxDetSimChannel > &, const std::vector< sim::AuxDetSimChannel > &, size_t)
void	MergeSimPhotons (std::vector< sim::SimPhotons > &, const std::vector< sim::SimPhotons > &)
void	MergeSimPhotonsLite (std::vector< sim::SimPhotonsLite > &, const std::vector< sim::SimPhotonsLite > &)
const std::vector< std::vector< size_t > > &	GetMCParticleListMap ()

Private Member Functions
void	UpdateG4TrackIDRange (std::pair< int, int >, size_t)

Private Attributes
std::vector< int >	fG4TrackIDOffsets
std::vector< std::pair< int, int > >	fG4TrackIDRanges
std::vector< std::vector< size_t > >	fMCParticleListMap

Detailed Description

Definition at line 23 of file MergeSimSources.h.

Constructor & Destructor Documentation

sim::MergeSimSourcesUtility::MergeSimSourcesUtility

(

std::vector< int > const &

)

Title: MergeSimSources Utility Class Author: Wes Ketchum (wketc.nosp@m.hum@.nosp@m.lanl..nosp@m.gov)

Description: Class that merges different simulation sources together to created a combined sim list. Typically just merges vectors/maps/etc together. But, if anything as a G4 trackID, applies a user-defined offset to those IDs.

Definition at line 18 of file MergeSimSources.cxx.

References fG4TrackIDOffsets, and Reset().

{
  fG4TrackIDOffsets = offsets;
  Reset();
}

Member Function Documentation

const std::vector< std::vector<size_t> >& sim::MergeSimSourcesUtility::GetMCParticleListMap ( )

inline

Definition at line 49 of file MergeSimSources.h.

References fMCParticleListMap.

Referenced by sim::MergeSimSources::produce().

{ return fMCParticleListMap; }

void sim::MergeSimSourcesUtility::MergeAuxDetSimChannels	(	std::vector< sim::AuxDetSimChannel > &	merged_vector,
		const std::vector< sim::AuxDetSimChannel > &	input_vector,
		size_t	source_index
	)

Definition at line 90 of file MergeSimSources.cxx.

References fG4TrackIDOffsets, max, min, and UpdateG4TrackIDRange().

Referenced by sim::MergeSimSources::produce().

{
  if(source_index >= fG4TrackIDOffsets.size())
    std::runtime_error("ERROR in MergeSimSourcesUtility: Source index out of range!");

  merged_vector.reserve( merged_vector.size() + input_vector.size() );

  std::pair<int,int> range_trackID(std::numeric_limits<int>::max(),
                                   std::numeric_limits<int>::min());
  
  for(auto const& simchannel : input_vector){
    std::vector<sim::AuxDetSimChannel>::iterator it = std::find(merged_vector.begin(),merged_vector.end(),simchannel);

    if(it==merged_vector.end()){
      merged_vector.emplace_back(simchannel.AuxDetID(), simchannel.AuxDetSensitiveID());
      it = merged_vector.end() - 1;
    }
    
    std::pair<int,int> thisrange = it->MergeAuxDetSimChannel(simchannel,fG4TrackIDOffsets[source_index]);
    if(thisrange.first < range_trackID.first) range_trackID.first = thisrange.first;
    if(thisrange.second > range_trackID.second) range_trackID.second = thisrange.second;
  }
  
  UpdateG4TrackIDRange(range_trackID,source_index);
}

void sim::MergeSimSourcesUtility::MergeMCParticles	(	std::vector< simb::MCParticle > &	merged_vector,
		const std::vector< simb::MCParticle > &	input_vector,
		size_t	source_index
	)

Definition at line 33 of file MergeSimSources.cxx.

References fG4TrackIDOffsets, fMCParticleListMap, max, min, and UpdateG4TrackIDRange().

Referenced by sim::MergeSimSources::produce().

{

  if(source_index >= fG4TrackIDOffsets.size())
    std::runtime_error("ERROR in MergeSimSourcesUtility: Source index out of range!");

  fMCParticleListMap[source_index].resize(input_vector.size());
  merged_vector.reserve(merged_vector.size() + input_vector.size());

  std::pair<int,int> range_trackID(std::numeric_limits<int>::max(),
                                   std::numeric_limits<int>::min());

  for(size_t i_p=0; i_p<input_vector.size(); i_p++){
    merged_vector.emplace_back(input_vector[i_p],fG4TrackIDOffsets[source_index]);

    fMCParticleListMap[source_index][i_p] = merged_vector.size() - 1;

    if(merged_vector.back().TrackId() < range_trackID.first) 
      range_trackID.first = merged_vector.back().TrackId();
    if(merged_vector.back().TrackId() > range_trackID.second) 
      range_trackID.second = merged_vector.back().TrackId();

  }

  UpdateG4TrackIDRange(range_trackID,source_index);
}

void sim::MergeSimSourcesUtility::MergeSimChannels	(	std::vector< sim::SimChannel > &	merged_vector,
		const std::vector< sim::SimChannel > &	input_vector,
		size_t	source_index
	)

Definition at line 62 of file MergeSimSources.cxx.

References fG4TrackIDOffsets, max, min, and UpdateG4TrackIDRange().

Referenced by sim::MergeSimSources::produce().

{
  if(source_index >= fG4TrackIDOffsets.size())
    std::runtime_error("ERROR in MergeSimSourcesUtility: Source index out of range!");

  merged_vector.reserve( merged_vector.size() + input_vector.size() );

  std::pair<int,int> range_trackID(std::numeric_limits<int>::max(),
                                   std::numeric_limits<int>::min());
  
  for(auto const& simchannel : input_vector){
    std::vector<sim::SimChannel>::iterator it = std::find(merged_vector.begin(),merged_vector.end(),simchannel);

    if(it==merged_vector.end()){
      merged_vector.emplace_back(simchannel.Channel());
      it = merged_vector.end() - 1;
    }

    std::pair<int,int> thisrange = it->MergeSimChannel(simchannel,fG4TrackIDOffsets[source_index]);
    if(thisrange.first < range_trackID.first) range_trackID.first = thisrange.first;
    if(thisrange.second > range_trackID.second) range_trackID.second = thisrange.second;
  }

  UpdateG4TrackIDRange(range_trackID,source_index);
}

void sim::MergeSimSourcesUtility::MergeSimPhotons	(	std::vector< sim::SimPhotons > &	merged_vector,
		const std::vector< sim::SimPhotons > &	input_vector
	)

Definition at line 118 of file MergeSimSources.cxx.

Referenced by sim::MergeSimSources::produce().

{

  merged_vector.reserve( merged_vector.size() + input_vector.size() );

  for(auto const& simphotons : input_vector){
    std::vector<sim::SimPhotons>::iterator it = std::find(merged_vector.begin(),merged_vector.end(),simphotons);

    if(it==merged_vector.end()){
      merged_vector.emplace_back(simphotons.OpChannel());
      it = merged_vector.end() - 1;
    }

    *it += simphotons;
  }
}

void sim::MergeSimSourcesUtility::MergeSimPhotonsLite	(	std::vector< sim::SimPhotonsLite > &	merged_vector,
		const std::vector< sim::SimPhotonsLite > &	input_vector
	)

Definition at line 136 of file MergeSimSources.cxx.

Referenced by sim::MergeSimSources::produce().

{

  merged_vector.reserve( merged_vector.size() + input_vector.size() );

  for(auto const& simphotons : input_vector){
    std::vector<sim::SimPhotonsLite>::iterator it = std::find(merged_vector.begin(),merged_vector.end(),simphotons);
    
    if(it==merged_vector.end()){
      merged_vector.emplace_back(simphotons.OpChannel);
      it = merged_vector.end() - 1;
    }
    
    *it += simphotons;
  }
}

void sim::MergeSimSourcesUtility::Reset

(

)

Definition at line 24 of file MergeSimSources.cxx.

References fG4TrackIDOffsets, fG4TrackIDRanges, fMCParticleListMap, max, and min.

Referenced by MergeSimSourcesUtility(), and sim::MergeSimSources::produce().

{
  fG4TrackIDRanges.resize(fG4TrackIDOffsets.size(),
                          std::make_pair(std::numeric_limits<int>::max(),
                                         std::numeric_limits<int>::min()));
  fMCParticleListMap.resize(fG4TrackIDOffsets.size(),
                            std::vector<size_t>());
}

void sim::MergeSimSourcesUtility::UpdateG4TrackIDRange ( std::pair< int, int >  newrange, size_t  source_index  )

private

Definition at line 154 of file MergeSimSources.cxx.

References fG4TrackIDOffsets, fG4TrackIDRanges, and ss.

Referenced by MergeAuxDetSimChannels(), MergeMCParticles(), and MergeSimChannels().

{
  if(source_index >= fG4TrackIDOffsets.size())
    std::runtime_error("ERROR in MergeSimSourcesUtility: Source index out of range!");

  if( newrange.first >= fG4TrackIDRanges[source_index].first &&
      newrange.second <= fG4TrackIDRanges[source_index].second)
    return;
  
  for(size_t i=0; i<fG4TrackIDRanges.size(); i++){
    if(i==source_index) continue;

    if( (newrange.first >= fG4TrackIDRanges[i].first && newrange.first <= fG4TrackIDRanges[i].second) ||
        (newrange.second >= fG4TrackIDRanges[i].first && newrange.second <= fG4TrackIDRanges[i].second) )
      {
        std::stringstream ss;
        ss << "ERROR in MergeSimSourcesUtility: Source trackIDs overlap!"
           << "\n\t" << i << "\t" << fG4TrackIDRanges[i].first << " " << fG4TrackIDRanges[i].second
           << "\n\t" << "n\t" << newrange.first << " " << newrange.second;
        throw std::runtime_error(ss.str());
      }
  }

  if(newrange.first < fG4TrackIDRanges[source_index].first)
    fG4TrackIDRanges[source_index].first = newrange.first;
  if(newrange.second > fG4TrackIDRanges[source_index].second)
    fG4TrackIDRanges[source_index].second = newrange.second;
  

}

Member Data Documentation

std::vector<int> sim::MergeSimSourcesUtility::fG4TrackIDOffsets

private

Definition at line 53 of file MergeSimSources.h.

Referenced by MergeAuxDetSimChannels(), MergeMCParticles(), MergeSimChannels(), MergeSimSourcesUtility(), Reset(), and UpdateG4TrackIDRange().

std::vector< std::pair<int,int> > sim::MergeSimSourcesUtility::fG4TrackIDRanges

private

Definition at line 54 of file MergeSimSources.h.

Referenced by Reset(), and UpdateG4TrackIDRange().

std::vector< std::vector<size_t> > sim::MergeSimSourcesUtility::fMCParticleListMap

private

Definition at line 56 of file MergeSimSources.h.

Referenced by GetMCParticleListMap(), MergeMCParticles(), and Reset().

---

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

• larsim/v06_53_00/source/larsim/LArG4/MergeSimSources.h
• larsim/v06_53_00/source/larsim/LArG4/MergeSimSources.cxx