SurveyPxbDicer.h

Go to the documentation of this file.

#ifndef GUARD_surveypxbdicer_h
#define GUARD_surveypxbdicer_h
 
#include <vector>
#include <functional>
#include <utility>
#include <fstream>
#include "DataFormats/GeometryVector/interface/Point3DBase.h"
#include "Alignment/SurveyAnalysis/interface/SurveyPxbImage.h"
#include "Math/SMatrix.h"
#include "Math/SVector.h"
#include "CLHEP/Random/RandGauss.h"
#include "CLHEP/Random/Random.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
 
#include <iostream>
#include <stdexcept>
 
class SurveyPxbDicer {
public:
  typedef SurveyPxbImage::coord_t coord_t;
  typedef double value_t;
  typedef std::vector<coord_t> fidpoint_t;
  typedef unsigned int count_t;
  typedef SurveyPxbImage::id_t id_t;
  typedef std::pair<id_t, id_t> idPair_t;
 
  // Constructors
  SurveyPxbDicer(){};
  // Constructor from VPSet
  SurveyPxbDicer(const std::vector<edm::ParameterSet> &pars, unsigned int seed);
 
  std::string doDice(const fidpoint_t &fidpointvec, const idPair_t &id, const bool rotate = false);
  void doDice(const fidpoint_t &fidpointvec, const idPair_t &id, std::ofstream &outfile, const bool rotate = false);
 
private:
  value_t ranGauss(value_t mean, value_t sigma) { return CLHEP::RandGauss::shoot(mean, sigma); };
 
  value_t mean_a0, sigma_a0;
  value_t mean_a1, sigma_a1;
  value_t mean_scale, sigma_scale;
  value_t mean_phi, sigma_phi;
  value_t mean_x, sigma_x;
  value_t mean_y, sigma_y;
 
  value_t getParByName(const std::string &name, const std::string &par, const std::vector<edm::ParameterSet> &pars);
 
  coord_t transform(const coord_t &x, const value_t &a0, const value_t &a1, const value_t &a2, const value_t &a3);
 
  struct findParByName {
    bool operator()(const std::string &name, const edm::ParameterSet &parset) const {
      return (parset.getParameter<std::string>("name") == name);
    }
  };
};
 
#endif