CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
List of all members | Classes | Public Types | Public Member Functions | Private Member Functions | Private Attributes
SurveyPxbDicer Class Reference

#include <SurveyPxbDicer.h>

Classes

struct  findParByName
 Function object for searching for a parameter in the VPSet. More...
 

Public Types

typedef SurveyPxbImage::coord_t coord_t
 
typedef unsigned int count_t
 
typedef std::vector< coord_tfidpoint_t
 
typedef SurveyPxbImage::id_t id_t
 
typedef std::pair< id_t, id_tidPair_t
 
typedef double value_t
 

Public Member Functions

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)
 
 SurveyPxbDicer ()
 
 SurveyPxbDicer (const std::vector< edm::ParameterSet > &pars, unsigned int seed)
 

Private Member Functions

value_t getParByName (const std::string &name, const std::string &par, const std::vector< edm::ParameterSet > &pars)
 
value_t ranGauss (value_t mean, value_t sigma)
 
coord_t transform (const coord_t &x, const value_t &a0, const value_t &a1, const value_t &a2, const value_t &a3)
 

Private Attributes

value_t mean_a0
 
value_t mean_a1
 
value_t mean_phi
 
value_t mean_scale
 
value_t mean_x
 
value_t mean_y
 
value_t sigma_a0
 
value_t sigma_a1
 
value_t sigma_phi
 
value_t sigma_scale
 
value_t sigma_x
 
value_t sigma_y
 

Detailed Description

Class to dice a picture from a given set of fiducial points This class has its use for toy MC simulations to validate the PXB survey

Definition at line 21 of file SurveyPxbDicer.h.

Member Typedef Documentation

Definition at line 24 of file SurveyPxbDicer.h.

typedef unsigned int SurveyPxbDicer::count_t

Definition at line 27 of file SurveyPxbDicer.h.

typedef std::vector<coord_t> SurveyPxbDicer::fidpoint_t

Definition at line 26 of file SurveyPxbDicer.h.

Definition at line 28 of file SurveyPxbDicer.h.

typedef std::pair<id_t,id_t> SurveyPxbDicer::idPair_t

Definition at line 29 of file SurveyPxbDicer.h.

typedef double SurveyPxbDicer::value_t

Definition at line 25 of file SurveyPxbDicer.h.

Constructor & Destructor Documentation

SurveyPxbDicer::SurveyPxbDicer ( )
inline

Definition at line 32 of file SurveyPxbDicer.h.

32 {};
SurveyPxbDicer::SurveyPxbDicer ( const std::vector< edm::ParameterSet > &  pars,
unsigned int  seed 
)

Definition at line 21 of file SurveyPxbDicer.cc.

References getParByName(), mean_a0, mean_a1, mean_phi, mean_scale, mean_x, mean_y, sigma_a0, sigma_a1, sigma_phi, sigma_scale, sigma_x, and sigma_y.

22 {
23  CLHEP::HepRandom::setTheSeed( seed );
24  mean_a0 = getParByName("a0", "mean", pars);
25  sigma_a0 = getParByName("a0", "sigma", pars);
26  mean_a1 = getParByName("a1", "mean", pars);
27  sigma_a1 = getParByName("a1", "sigma", pars);
28  mean_scale = getParByName("scale", "mean", pars);
29  sigma_scale = getParByName("scale", "sigma", pars);
30  mean_phi = getParByName("phi", "mean", pars);
31  sigma_phi = getParByName("phi", "sigma", pars);
32  mean_x = getParByName("x", "mean", pars);
33  sigma_x = getParByName("x", "sigma", pars);
34  mean_y = getParByName("y", "mean", pars);
35  sigma_y = getParByName("y", "sigma", pars);
36 }
value_t mean_scale
value_t sigma_scale
value_t getParByName(const std::string &name, const std::string &par, const std::vector< edm::ParameterSet > &pars)

Member Function Documentation

std::string SurveyPxbDicer::doDice ( const fidpoint_t fidpointvec,
const idPair_t id,
const bool  rotate = false 
)

Invoke the dicer

Parameters
fidpointvecvector with fiducial points where values need to be diced for and transformed to the photo fram
idPairpair of the id values

Definition at line 38 of file SurveyPxbDicer.cc.

References funct::cos(), mean_a0, mean_a1, mean_phi, mean_scale, p1, p2, p3, phi, ranGauss(), pileupReCalc_HLTpaths::scale, edm::second(), sigma_a0, sigma_a1, sigma_phi, sigma_scale, sigma_x, sigma_y, jetcorrextractor::sign(), funct::sin(), transform(), x, PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().

Referenced by MillePedeAlignmentAlgorithm::addPxbSurvey(), and doDice().

39 {
40  // Dice the local parameters
41  const value_t a0 = ranGauss(mean_a0, sigma_a0);
42  const value_t a1 = ranGauss(mean_a1, sigma_a1);
45  const value_t a2 = scale * cos(phi);
46  const value_t a3 = scale * sin(phi);
47  const coord_t p0 = transform(fidpointvec[0],a0,a1,a2,a3);
48  const coord_t p1 = transform(fidpointvec[1],a0,a1,a2,a3);
49  const coord_t p2 = transform(fidpointvec[2],a0,a1,a2,a3);
50  const coord_t p3 = transform(fidpointvec[3],a0,a1,a2,a3);
51  const value_t sign = rotate ? -1 : 1;
52  std::ostringstream oss;
53  oss << id.first << " "
54  << sign*ranGauss(p0.x(),sigma_x) << " " << -sign*ranGauss(p0.y(),sigma_y) << " "
55  << sign*ranGauss(p1.x(),sigma_x) << " " << -sign*ranGauss(p1.y(),sigma_y) << " "
56  << id.second << " "
57  << sign*ranGauss(p2.x(),sigma_x) << " " << -sign*ranGauss(p2.y(),sigma_y) << " "
58  << sign*ranGauss(p3.x(),sigma_x) << " " << -sign*ranGauss(p3.y(),sigma_y) << " "
59  << sigma_x << " " << sigma_y << " "
60  << rotate
61  << " # MC-truth:"
62  << " a0-a3: " << a0 << " " << a1 << " " << a2 << " " << a3
63  << " S: " << scale
64  << " phi: " << phi
65  << " x0: " << fidpointvec[0].x() << " " << fidpointvec[0].y()
66  << " x1: " << fidpointvec[1].x() << " " << fidpointvec[1].y()
67  << " x2: " << fidpointvec[2].x() << " " << fidpointvec[2].y()
68  << " x3: " << fidpointvec[3].x() << " " << fidpointvec[3].y()
69  << std::endl;
70  return oss.str();
71 }
value_t mean_scale
double sign(double x)
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
value_t sigma_scale
coord_t transform(const coord_t &x, const value_t &a0, const value_t &a1, const value_t &a2, const value_t &a3)
U second(std::pair< T, U > const &p)
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
value_t ranGauss(value_t mean, value_t sigma)
double p2[4]
Definition: TauolaWrapper.h:90
SurveyPxbImage::coord_t coord_t
def rotate
Definition: svgfig.py:704
double p1[4]
Definition: TauolaWrapper.h:89
double p3[4]
Definition: TauolaWrapper.h:91
void SurveyPxbDicer::doDice ( const fidpoint_t fidpointvec,
const idPair_t id,
std::ofstream &  outfile,
const bool  rotate = false 
)

Definition at line 73 of file SurveyPxbDicer.cc.

References doDice().

74 {
75  outfile << doDice(fidpointvec, id, rotate);
76 }
std::string doDice(const fidpoint_t &fidpointvec, const idPair_t &id, const bool rotate=false)
list outfile
Definition: EdgesToViz.py:91
def rotate
Definition: svgfig.py:704
SurveyPxbDicer::value_t SurveyPxbDicer::getParByName ( const std::string &  name,
const std::string &  par,
const std::vector< edm::ParameterSet > &  pars 
)
private

Gets parameter by name from the VPSet

Parameters
namename of the parameter to be searched for in field 'name' of the VPSet
parselects the value, i.e. mean or sigma
parsreference to VPSet

Definition at line 83 of file SurveyPxbDicer.cc.

Referenced by SurveyPxbDicer().

84 {
85  std::vector<edm::ParameterSet>::const_iterator it;
86  it = std::find_if(pars.begin(), pars.end(), std::bind1st(findParByName(),name));
87  if (it == pars.end()) { throw std::runtime_error("Parameter not found in SurveyPxbDicer::getParByName"); }
88  return (*it).getParameter<value_t>(par);
89 }
value_t SurveyPxbDicer::ranGauss ( value_t  mean,
value_t  sigma 
)
inlineprivate

invoke the RNG to geat a gaussian smeared value

Parameters
meanmean value
sigma

Definition at line 46 of file SurveyPxbDicer.h.

Referenced by doDice().

46 {return CLHEP::RandGauss::shoot(mean,sigma);};
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 
)
private

Transform the diced values to the frame of a toy photograph

Parameters
xcoordinate to be transformed from local frame to photo frame
a0Transformation parameter, shift in x
a1Transformation parameter, shift in y
a2Transformation parameter, scale*cos(phi)
a3Transformation parameter, scale*sin(phi)

Definition at line 78 of file SurveyPxbDicer.cc.

References PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().

Referenced by Vispa.Views.LineDecayView.DecayLine::boundingRect(), Vispa.Views.LineDecayView.DecayLine::containsPoint(), doDice(), and Vispa.Views.LineDecayView.DecayLine::paint().

79 {
80  return coord_t(a0+a2*x.x()+a3*x.y(), a1-a3*x.x()+a2*x.y());
81 }
SurveyPxbImage::coord_t coord_t

Member Data Documentation

value_t SurveyPxbDicer::mean_a0
private

Definition at line 46 of file SurveyPxbDicer.h.

Referenced by doDice(), and SurveyPxbDicer().

value_t SurveyPxbDicer::mean_a1
private

Definition at line 49 of file SurveyPxbDicer.h.

Referenced by doDice(), and SurveyPxbDicer().

value_t SurveyPxbDicer::mean_phi
private

Definition at line 51 of file SurveyPxbDicer.h.

Referenced by doDice(), and SurveyPxbDicer().

value_t SurveyPxbDicer::mean_scale
private

Definition at line 50 of file SurveyPxbDicer.h.

Referenced by doDice(), and SurveyPxbDicer().

value_t SurveyPxbDicer::mean_x
private

Definition at line 52 of file SurveyPxbDicer.h.

Referenced by SurveyPxbDicer().

value_t SurveyPxbDicer::mean_y
private

Definition at line 53 of file SurveyPxbDicer.h.

Referenced by SurveyPxbDicer().

value_t SurveyPxbDicer::sigma_a0
private

Definition at line 46 of file SurveyPxbDicer.h.

Referenced by doDice(), and SurveyPxbDicer().

value_t SurveyPxbDicer::sigma_a1
private

Definition at line 49 of file SurveyPxbDicer.h.

Referenced by doDice(), and SurveyPxbDicer().

value_t SurveyPxbDicer::sigma_phi
private

Definition at line 51 of file SurveyPxbDicer.h.

Referenced by doDice(), and SurveyPxbDicer().

value_t SurveyPxbDicer::sigma_scale
private

Definition at line 50 of file SurveyPxbDicer.h.

Referenced by doDice(), and SurveyPxbDicer().

value_t SurveyPxbDicer::sigma_x
private

Definition at line 52 of file SurveyPxbDicer.h.

Referenced by doDice(), and SurveyPxbDicer().

value_t SurveyPxbDicer::sigma_y
private

Definition at line 53 of file SurveyPxbDicer.h.

Referenced by doDice(), and SurveyPxbDicer().