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