/data/refman/pasoursint/CMSSW_6_1_1/src/Alignment/SurveyAnalysis/interface/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(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: public std::binary_function< std::string, edm::ParameterSet, bool >
        {
                bool operator () ( const std::string &name, const edm::ParameterSet &parset )
                        const { return (parset.getParameter<std::string>("name") == name);}
        };

};

#endif