---

SamplerRouter.cxx

Go to the documentation of this file.

// ================================================
// Implementation of SamplerRouter
// ================================================

#include "FastShowerUtils/SamplerRouter.h"
#include "FastShowerUtils/ParticleParameters.h"
#include <assert.h>
#include <iostream>

namespace FastShower{
  SamplerRouter::SamplerRouter(const 
                               std::vector<std::pair<double, ISampler*> >& v,
                               double x, std::string label):
    DebugBase(label), m_floor(x){
    m_samplers.insert(v.begin(), v.end());
    cout<<"Sampler router finishing construction"<<endl;
  }
  //
  SamplerRouter::SamplerRouter(const SamplerRouter& c):ISampler(), DebugBase(){
    SamplerTable::const_iterator iter = c.m_samplers.begin();
    SamplerTable::const_iterator end  = c.m_samplers.end();
    
    for(; iter!=end; ++iter){
      m_samplers[(*iter).first]=(*iter).second->clone();
    }
    m_floor = c.m_floor;
    cout<<"CopyConstructed SamplersRouter, size "<<m_samplers.size()<<endl;
  }
  //
  void SamplerRouter::sample(const PolyArgs& pa, CoreSamples& cs) const{

    ParticleParameters pp = *pa.pp();
    double indQ = indexQ(pp, cs);
    assert ( indQ >= m_floor); 
    assert ( m_samplers.upper_bound(indQ) != m_samplers.end());

    ISampler* sampler = m_samplers.upper_bound(indQ)->second;
    sampler->sample(pa, cs);
  }
  //
  double SamplerRouter::lastValue(const CoreSamples& cs) const {return 0;}
  //   
  void SamplerRouter::components(std::vector<const IDebug*>& v) const{
    v.clear();
    std::map<double, ISampler*>::const_iterator iter = m_samplers.begin();
    std::map<double, ISampler*>::const_iterator end = m_samplers.end();
    for(;iter!=end;++iter){ v.push_back( (*iter).second ); }
   }
 }//namespace

---