---

DistRandomiser1D.cxx

Go to the documentation of this file.

#include "FastShowerUtils/Normalisers/DistRandomiser1D.h"

// CLHEP
#include "CLHEP/Matrix/Vector.h"

#include <cmath>
#include <string>
#include <algorithm>
#include <iostream>
#include <fstream>
namespace FastShower{

  DistRandomiser1D::DistRandomiser1D(const std::string& f){

    // get data form file (histogram)
    std::vector<double> v;
    readData(f,v);
    
    // needed to initialise RandGeneral
    std::vector<double>::size_type n =  v.size(); // # of bins
    double dist[50];
    assert(n == 50);
    // 
    int i;
    std::vector<double>::const_iterator itr = v.begin();
    std::vector<double>::const_iterator end = v.end();
    for (i=0, itr; itr!=end; ++i, ++itr){
      dist[i] = (*itr);
    }
    m_Dist = new RandGeneral(dist, n, 0);
  }
  //
  void
  DistRandomiser1D::readData(const std::string& f, std::vector<double>& v) const {
    ifstream data(f.c_str());
    istream_iterator<int> dataItr(data);
    istream_iterator<int> endOfData;
    copy(dataItr,endOfData,back_inserter(v));
  }
  //
  double DistRandomiser1D::sample() const {
    return 2.0*m_Dist->fire();
  }
  //
  DistRandomiser1D::~DistRandomiser1D(){
    delete m_Dist;
  }
}//namespace

---