MergeSimSources.h

Go to the documentation of this file.

#ifndef MERGESIMSOURCES_H
#define MERGESIMSOURCES_H

#include "nusimdata/SimulationBase/MCParticle.h"

#include "lardataobj/Simulation/SimChannel.h"
#include "lardataobj/Simulation/SimPhotons.h"
#include "lardataobj/Simulation/AuxDetSimChannel.h"

namespace sim{

  class MergeSimSourcesUtility{

  public:

    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() { return fMCParticleListMap; }
    
  private:

    std::vector<int>                   fG4TrackIDOffsets;
    std::vector< std::pair<int,int> >  fG4TrackIDRanges;

    std::vector< std::vector<size_t> > fMCParticleListMap;

    void UpdateG4TrackIDRange(std::pair<int,int>,size_t);
    
  }; //end MergeSimSourcesUtility class

} //end namespace sim

#endif