SurveyPxbDicer.cc

Go to the documentation of this file.

#include "Alignment/SurveyAnalysis/interface/SurveyPxbImage.h"
#include "Alignment/SurveyAnalysis/interface/SurveyPxbDicer.h"

//#include <stdexcept>
#include <vector>
#include <cmath>
#include <string>
#include <sstream>
#include <functional>
#include <algorithm>
#include <fstream>
#include "DataFormats/GeometryVector/interface/Point3DBase.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>

SurveyPxbDicer::SurveyPxbDicer(std::vector<edm::ParameterSet> pars, unsigned int seed)
{
    CLHEP::HepRandom::setTheSeed( seed );
    mean_a0 = getParByName("a0", "mean", pars);
    sigma_a0 = getParByName("a0", "sigma", pars);
    mean_a1 = getParByName("a1", "mean", pars);
    sigma_a1 = getParByName("a1", "sigma", pars);
    mean_scale = getParByName("scale", "mean", pars);
    sigma_scale = getParByName("scale", "sigma", pars);
    mean_phi = getParByName("phi", "mean", pars);
    sigma_phi = getParByName("phi", "sigma", pars);
    mean_x = getParByName("x", "mean", pars);
    sigma_x = getParByName("x", "sigma", pars);
    mean_y = getParByName("y", "mean", pars);
    sigma_y = getParByName("y", "sigma", pars);
}

std::string SurveyPxbDicer::doDice(const fidpoint_t &fidpointvec, const idPair_t &id, const bool rotate)
{
    // Dice the local parameters
    const value_t a0 = ranGauss(mean_a0, sigma_a0);
    const value_t a1 = ranGauss(mean_a1, sigma_a1);
    const value_t scale = ranGauss(mean_scale, sigma_scale);
    const value_t phi = ranGauss(mean_phi, sigma_phi);
    const value_t a2 = scale * cos(phi);
    const value_t a3 = scale * sin(phi);
    const coord_t p0 = transform(fidpointvec[0],a0,a1,a2,a3);
    const coord_t p1 = transform(fidpointvec[1],a0,a1,a2,a3);
    const coord_t p2 = transform(fidpointvec[2],a0,a1,a2,a3);
    const coord_t p3 = transform(fidpointvec[3],a0,a1,a2,a3);
    const value_t sign = rotate ? -1 : 1;
    std::ostringstream oss;
    oss << id.first << " " 
         << sign*ranGauss(p0.x(),sigma_x) << " " << -sign*ranGauss(p0.y(),sigma_y) << " "
         << sign*ranGauss(p1.x(),sigma_x) << " " << -sign*ranGauss(p1.y(),sigma_y) << " "
         << id.second << " "
         << sign*ranGauss(p2.x(),sigma_x) << " " << -sign*ranGauss(p2.y(),sigma_y) << " "
         << sign*ranGauss(p3.x(),sigma_x) << " " << -sign*ranGauss(p3.y(),sigma_y) << " "
         << sigma_x << " " << sigma_y << " "
         << rotate
         << " # MC-truth:"
         << " a0-a3: " << a0 << " " << a1 << " " << a2 << " " << a3
         << " S: " << scale
         << " phi: " << phi
         << " x0: " << fidpointvec[0].x() << " " << fidpointvec[0].y()
         << " x1: " << fidpointvec[1].x() << " " << fidpointvec[1].y()
         << " x2: " << fidpointvec[2].x() << " " << fidpointvec[2].y()
         << " x3: " << fidpointvec[3].x() << " " << fidpointvec[3].y()
         << std::endl;
    return oss.str();
}

void SurveyPxbDicer::doDice(const fidpoint_t &fidpointvec, const idPair_t &id, std::ofstream &outfile, const bool rotate)
{
    outfile << doDice(fidpointvec, id, rotate);
}

SurveyPxbDicer::coord_t SurveyPxbDicer::transform(const coord_t &x, const value_t &a0, const value_t &a1, const value_t &a2, const value_t &a3)
{
    return coord_t(a0+a2*x.x()+a3*x.y(), a1-a3*x.x()+a2*x.y());
}

SurveyPxbDicer::value_t SurveyPxbDicer::getParByName(const std::string &name, const std::string &par, const std::vector<edm::ParameterSet>& pars)
{
    std::vector<edm::ParameterSet>::const_iterator it;
    it = std::find_if(pars.begin(), pars.end(), std::bind1st(findParByName(),name));
    if (it == pars.end()) { throw std::runtime_error("Parameter not found in SurveyPxbDicer::getParByName"); }
    return (*it).getParameter<value_t>(par);
}