CMS 3D CMS Logo
List of all members | Public Member Functions | Static Public Member Functions | Private Attributes
SiPixelGenError Class Reference

#include <SiPixelGenError.h>

Public Member Functions

float fbin (int i)
 Return lower bound of Qbin definition. More...
 
float lorxbias ()
 signed lorentz x-bias (microns) More...
 
float lorxwidth ()
 signed lorentz x-width (microns) More...
 
float lorybias ()
 signed lorentz y-bias (microns) More...
 
float lorywidth ()
 signed lorentz y-width (microns) More...
 
int qbin (int id, float cotalpha, float cotbeta, float locBz, float locBx, float qclus, bool irradiationCorrections, int &pixmx, float &sigmay, float &deltay, float &sigmax, float &deltax, float &sy1, float &dy1, float &sy2, float &dy2, float &sx1, float &dx1, float &sx2, float &dx2)
 
int qbin (int id, float cotalpha, float cotbeta, float locBz, float locBx, float qclus, float &pixmx, float &sigmay, float &deltay, float &sigmax, float &deltax, float &sy1, float &dy1, float &sy2, float &dy2, float &sx1, float &dx1, float &sx2, float &dx2)
 
int qbin (int id)
 
 SiPixelGenError (const std::vector< SiPixelGenErrorStore > &thePixelTemp)
 Constructor for cases in which template store already exists. More...
 
float xsize ()
 pixel x-size (microns) More...
 
float ysize ()
 pixel y-size (microns) More...
 
float zsize ()
 pixel z-size or thickness (microns) More...
 

Static Public Member Functions

static void postInit (std::vector< SiPixelGenErrorStore > &thePixelTemp_)
 
static bool pushfile (int filenum, std::vector< SiPixelGenErrorStore > &pixelTemp, std::string dir="")
 
static bool pushfile (const SiPixelGenErrorDBObject &dbobject, std::vector< SiPixelGenErrorStore > &pixelTemp)
 

Private Attributes

float fbin_ [3]
 The QBin definitions in Q_clus/Q_avg. More...
 
int id_current_
 current id More...
 
int index_id_
 current index More...
 
float lorxbias_
 Lorentz x-width. More...
 
float lorxwidth_
 Lorentz x-width. More...
 
float lorybias_
 Lorentz y-width (sign corrected for fpix frame) More...
 
float lorywidth_
 Lorentz y-width (sign corrected for fpix frame) More...
 
const std::vector< SiPixelGenErrorStore > & thePixelTemp_
 
float xsize_
 Pixel x-size. More...
 
float ysize_
 Pixel y-size. More...
 
float zsize_
 Pixel z-size (thickness) More...
 

Detailed Description

A Generic Algorithm info management class. SiPixelGenError contains thePixelTemp (a std::vector of SiPixelGenErrorStore, each of which is a collection of many SiPixelGenErrorEntries). Each SiPixelGenErrorStore corresponds to a given detector condition, and is valid for a range of runs. We allow more than one Store since the may change over time.

This class reads templates from files via pushfile() method.

The main functionality of SiPixelGenError is qbin(), which produces algorithm info on the fly, given a specific track's alpha and beta.

Definition at line 117 of file SiPixelGenError.h.

Constructor & Destructor Documentation

SiPixelGenError::SiPixelGenError ( const std::vector< SiPixelGenErrorStore > &  thePixelTemp)
inline

Constructor for cases in which template store already exists.

Definition at line 119 of file SiPixelGenError.h.

References SiPixelGenErrorEntry::cotalpha, SiPixelGenErrorEntry::cotbeta, dir, and AlCaHLTBitMon_QueryRunRegistry::string.

Member Function Documentation

float SiPixelGenError::fbin ( int  i)
inline

Return lower bound of Qbin definition.

Definition at line 150 of file SiPixelGenError.h.

References Exception, and mps_fire::i.

float SiPixelGenError::lorxbias ( )
inline

signed lorentz x-bias (microns)

Definition at line 148 of file SiPixelGenError.h.

float SiPixelGenError::lorxwidth ( )
inline

signed lorentz x-width (microns)

Definition at line 146 of file SiPixelGenError.h.

Referenced by PixelCPEGeneric::localPosition().

float SiPixelGenError::lorybias ( )
inline

signed lorentz y-bias (microns)

Definition at line 147 of file SiPixelGenError.h.

float SiPixelGenError::lorywidth ( )
inline

signed lorentz y-width (microns)

Definition at line 145 of file SiPixelGenError.h.

Referenced by PixelCPEGeneric::localPosition().

void SiPixelGenError::postInit ( std::vector< SiPixelGenErrorStore > &  thePixelTemp_)
static

Definition at line 412 of file SiPixelGenError.cc.

412  {
413 
414  for (auto & templ : thePixelTemp_) {
415  for ( auto iy=0; iy<templ.head.NTy; ++iy ) templ.cotbetaY[iy]=templ.enty[iy].cotbeta;
416  for ( auto iy=0; iy<templ.head.NTyx; ++iy ) templ.cotbetaX[iy]= templ.entx[iy][0].cotbeta;
417  for ( auto ix=0; ix<templ.head.NTxx; ++ix ) templ.cotalphaX[ix]=templ.entx[0][ix].cotalpha;
418  }
419 
420 }
SiPixelGenError::thePixelTemp_
const std::vector< SiPixelGenErrorStore > & thePixelTemp_
Definition: SiPixelGenError.h:183
bool SiPixelGenError::pushfile ( int  filenum,
std::vector< SiPixelGenErrorStore > &  pixelTemp,
std::string  dir = "" 
)
static

This routine initializes the global GenError structures from an external file generror_summary_zpNNNN where NNNN are four digits of filenum.

Parameters
filenum- an integer NNNN used in the filename generror_summary_zpNNNN

Definition at line 57 of file SiPixelGenError.cc.

References SiPixelGenErrorHeader::Bfield, EnergyCorrector::c, SiPixelGenErrorEntry::cotalpha, SiPixelGenErrorStore::cotalphaX, SiPixelGenErrorEntry::cotbeta, SiPixelGenErrorStore::cotbetaX, SiPixelGenErrorStore::cotbetaY, SiPixelGenErrorHeader::Dtype, SiPixelGenErrorEntry::dxone, SiPixelGenErrorEntry::dxtwo, SiPixelGenErrorEntry::dyone, SiPixelGenErrorEntry::dytwo, ENDL, SiPixelGenErrorStore::entx, SiPixelGenErrorStore::enty, SiPixelGenErrorHeader::fbin, FrontierConditions_GlobalTag_cff::file, SiPixelGenErrorHeader::fluence, SiPixelGenErrorStore::head, mps_fire::i, SiPixelGenErrorHeader::ID, recoMuon::in, create_public_lumi_plots::in_file, gen::k, LOGERROR, LOGINFO, SiPixelGenErrorHeader::lorxbias, SiPixelGenErrorHeader::lorxwidth, SiPixelGenErrorHeader::lorybias, SiPixelGenErrorHeader::lorywidth, SiPixelGenErrorHeader::NTxx, SiPixelGenErrorHeader::NTy, SiPixelGenErrorHeader::NTyx, SiPixelGenErrorEntry::pixmax, SiPixelGenErrorEntry::qavg, SiPixelGenErrorEntry::qmin, SiPixelGenErrorEntry::qmin2, SiPixelGenErrorHeader::qscale, SiPixelGenErrorEntry::runnum, SiPixelGenErrorHeader::s50, SiPixelGenErrorHeader::ss50, AlCaHLTBitMon_QueryRunRegistry::string, SiPixelGenErrorEntry::sxone, SiPixelGenErrorEntry::sxtwo, SiPixelGenErrorEntry::syone, SiPixelGenErrorEntry::sytwo, SiPixelGenErrorHeader::temperature, SiPixelGenErrorHeader::templ_version, SiPixelGenErrorHeader::title, SiPixelGenErrorHeader::Vbias, SiPixelGenErrorEntry::xavggen, SiPixelGenErrorEntry::xrmsgen, SiPixelGenErrorHeader::xsize, SiPixelGenErrorEntry::yavggen, SiPixelGenErrorEntry::yrmsgen, SiPixelGenErrorHeader::ysize, and SiPixelGenErrorHeader::zsize.

Referenced by PixelCPEGeneric::PixelCPEGeneric().

58 {
59  // Add info stored in external file numbered filenum to theGenErrorStore
60 
61  // Local variables
62  int i, j, k;
63  float costrk[3]={0,0,0};
64  const char *tempfile;
65  // char title[80]; remove this
66  char c;
67  const int code_version={1};
68 
69 
70 
71  // Create a filename for this run
72 
73  std::ostringstream tout;
74 
75  // Create different path in CMSSW than standalone
76 
77 #ifndef SI_PIXEL_TEMPLATE_STANDALONE
78  tout << dir << "generror_summary_zp"
79  << std::setw(4) << std::setfill('0') << std::right << filenum << ".out" << std::ends;
80  std::string tempf = tout.str();
81  edm::FileInPath file( tempf.c_str() );
82  tempfile = (file.fullPath()).c_str();
83 #else
84  tout << "generror_summary_zp" << std::setw(4) << std::setfill('0') << std::right << filenum << ".out" << std::ends;
85  std::string tempf = tout.str();
86  tempfile = tempf.c_str();
87 #endif
88 
89  // open the Generic file
90 
91  std::ifstream in_file(tempfile, std::ios::in);
92 
93  if(in_file.is_open()) {
94 
95  // Create a local GenError storage entry
96 
97  SiPixelGenErrorStore theCurrentTemp;
98 
99  // Read-in a header string first and print it
100 
101  for (i=0; (c=in_file.get()) != '\n'; ++i) {
102  if(i < 79) {theCurrentTemp.head.title[i] = c;}
103  }
104  if(i > 78) {i=78;}
105  theCurrentTemp.head.title[i+1] ='\0';
106  LOGINFO("SiPixelGenError") << "Loading Pixel GenError File - " << theCurrentTemp.head.title << ENDL;
107 
108  // next, the header information
109 
110  in_file >> theCurrentTemp.head.ID >> theCurrentTemp.head.templ_version >> theCurrentTemp.head.Bfield >> theCurrentTemp.head.NTy >> theCurrentTemp.head.NTyx >> theCurrentTemp.head.NTxx >> theCurrentTemp.head.Dtype >> theCurrentTemp.head.Vbias >>
111  theCurrentTemp.head.temperature >> theCurrentTemp.head.fluence >> theCurrentTemp.head.qscale >> theCurrentTemp.head.s50 >>
112  theCurrentTemp.head.lorywidth >> theCurrentTemp.head.lorxwidth >> theCurrentTemp.head.ysize >> theCurrentTemp.head.xsize >>
113  theCurrentTemp.head.zsize >> theCurrentTemp.head.ss50 >> theCurrentTemp.head.lorybias >> theCurrentTemp.head.lorxbias >>
114  theCurrentTemp.head.fbin[0] >> theCurrentTemp.head.fbin[1] >> theCurrentTemp.head.fbin[2];
115 
116  if(in_file.fail()) {LOGERROR("SiPixelGenError") << "Error reading file, no GenError load" << ENDL; return false;}
117 
118  LOGINFO("SiPixelGenError") << "GenError ID = " << theCurrentTemp.head.ID << ", GenError Version " << theCurrentTemp.head.templ_version << ", Bfield = " << theCurrentTemp.head.Bfield
119  << ", NTy = " << theCurrentTemp.head.NTy << ", NTyx = " << theCurrentTemp.head.NTyx<< ", NTxx = " << theCurrentTemp.head.NTxx << ", Dtype = " << theCurrentTemp.head.Dtype
120  << ", Bias voltage " << theCurrentTemp.head.Vbias << ", temperature "
121  << theCurrentTemp.head.temperature << ", fluence " << theCurrentTemp.head.fluence << ", Q-scaling factor " << theCurrentTemp.head.qscale
122  << ", 1/2 multi dcol threshold " << theCurrentTemp.head.s50 << ", 1/2 single dcol threshold " << theCurrentTemp.head.ss50
123  << ", y Lorentz Width " << theCurrentTemp.head.lorywidth << ", y Lorentz Bias " << theCurrentTemp.head.lorybias
124  << ", x Lorentz width " << theCurrentTemp.head.lorxwidth << ", x Lorentz Bias " << theCurrentTemp.head.lorxbias
125  << ", Q/Q_avg fractions for Qbin defs " << theCurrentTemp.head.fbin[0] << ", " << theCurrentTemp.head.fbin[1]
126  << ", " << theCurrentTemp.head.fbin[2]
127  << ", pixel x-size " << theCurrentTemp.head.xsize << ", y-size " << theCurrentTemp.head.ysize << ", zsize " << theCurrentTemp.head.zsize << ENDL;
128 
129  if(theCurrentTemp.head.templ_version < code_version) {LOGERROR("SiPixelGenError") << "code expects version " << code_version << ", no GenError load" << ENDL; return false;}
130 
131 #ifdef SI_PIXEL_TEMPLATE_USE_BOOST
132 
133  // next, layout the 1-d/2-d structures needed to store GenError info
134 
135 
136  theCurrentTemp.cotbetaY = new float[theCurrentTemp.head.NTy];
137  theCurrentTemp.cotbetaX = new float[theCurrentTemp.head.NTyx];
138  theCurrentTemp.cotalphaX = new float[theCurrentTemp.head.NTxx];
139 
140  theCurrentTemp.enty.resize(boost::extents[theCurrentTemp.head.NTy]);
141 
142  theCurrentTemp.entx.resize(boost::extents[theCurrentTemp.head.NTyx][theCurrentTemp.head.NTxx]);
143 
144 #endif
145 
146  // next, loop over all y-angle entries
147 
148  for (i=0; i < theCurrentTemp.head.NTy; ++i) {
149 
150  in_file >> theCurrentTemp.enty[i].runnum >> costrk[0] >> costrk[1] >> costrk[2];
151 
152  if(in_file.fail()) {LOGERROR("SiPixelGenError") << "Error reading file 1, no GenError load, run # " << theCurrentTemp.enty[i].runnum << ENDL; return false;}
153 
154  // Calculate the alpha, beta, and cot(beta) for this entry
155 
156  theCurrentTemp.enty[i].cotalpha = costrk[0]/costrk[2];
157 
158  theCurrentTemp.enty[i].cotbeta = costrk[1]/costrk[2];
159 
160  in_file >> theCurrentTemp.enty[i].qavg >> theCurrentTemp.enty[i].pixmax >> theCurrentTemp.enty[i].dyone
161  >> theCurrentTemp.enty[i].syone >> theCurrentTemp.enty[i].dxone >> theCurrentTemp.enty[i].sxone;
162 
163  if(in_file.fail()) {LOGERROR("SiPixelGenError") << "Error reading file 2, no GenError load, run # " << theCurrentTemp.enty[i].runnum << ENDL; return false;}
164 
165  in_file >> theCurrentTemp.enty[i].dytwo >> theCurrentTemp.enty[i].sytwo >> theCurrentTemp.enty[i].dxtwo
166  >> theCurrentTemp.enty[i].sxtwo >> theCurrentTemp.enty[i].qmin >> theCurrentTemp.enty[i].qmin2;
167 
168  if(in_file.fail()) {LOGERROR("SiPixelGenError") << "Error reading file 3, no GenError load, run # " << theCurrentTemp.enty[i].runnum << ENDL; return false;}
169 
170  if(theCurrentTemp.enty[i].qmin <= 0.) {LOGERROR("SiPixelGenError") << "Error in GenError ID " << theCurrentTemp.head.ID << " qmin = " << theCurrentTemp.enty[i].qmin << ", run # " << theCurrentTemp.enty[i].runnum << ENDL; return false;}
171 
172 
173  for (j=0; j<4; ++j) {
174 
175  in_file >> theCurrentTemp.enty[i].yavggen[j] >> theCurrentTemp.enty[i].yrmsgen[j] >> theCurrentTemp.enty[i].xavggen[j] >> theCurrentTemp.enty[i].xrmsgen[j];
176 
177  if(in_file.fail()) {LOGERROR("SiPixelGenError") << "Error reading file 14a, no GenError load, run # " << theCurrentTemp.enty[i].runnum << ENDL; return false;}
178  }
179 
180 
181  }
182 
183  // next, loop over all barrel x-angle entries
184 
185  for (k=0; k < theCurrentTemp.head.NTyx; ++k) {
186 
187  for (i=0; i < theCurrentTemp.head.NTxx; ++i) {
188 
189  in_file >> theCurrentTemp.entx[k][i].runnum >> costrk[0] >> costrk[1] >> costrk[2];
190 
191  if(in_file.fail()) {LOGERROR("SiPixelGenError") << "Error reading file 17, no GenError load, run # " << theCurrentTemp.entx[k][i].runnum << ENDL; return false;}
192 
193  // Calculate the alpha, beta, and cot(beta) for this entry
194 
195  theCurrentTemp.entx[k][i].cotalpha = costrk[0]/costrk[2];
196 
197  theCurrentTemp.entx[k][i].cotbeta = costrk[1]/costrk[2];
198 
199  in_file >> theCurrentTemp.entx[k][i].qavg >> theCurrentTemp.entx[k][i].pixmax >> theCurrentTemp.entx[k][i].dyone
200  >> theCurrentTemp.entx[k][i].syone >> theCurrentTemp.entx[k][i].dxone >> theCurrentTemp.entx[k][i].sxone;
201 
202  if(in_file.fail()) {LOGERROR("SiPixelGenError") << "Error reading file 18, no GenError load, run # " << theCurrentTemp.entx[k][i].runnum << ENDL; return false;}
203 
204  in_file >> theCurrentTemp.entx[k][i].dytwo >> theCurrentTemp.entx[k][i].sytwo >> theCurrentTemp.entx[k][i].dxtwo
205  >> theCurrentTemp.entx[k][i].sxtwo >> theCurrentTemp.entx[k][i].qmin >> theCurrentTemp.entx[k][i].qmin2;
206  // >> theCurrentTemp.entx[k][i].mpvvav >> theCurrentTemp.entx[k][i].sigmavav >> theCurrentTemp.entx[k][i].kappavav;
207 
208  if(in_file.fail()) {LOGERROR("SiPixelGenError") << "Error reading file 19, no GenError load, run # " << theCurrentTemp.entx[k][i].runnum << ENDL; return false;}
209 
210 
211  for (j=0; j<4; ++j) {
212 
213  in_file >> theCurrentTemp.entx[k][i].yavggen[j] >> theCurrentTemp.entx[k][i].yrmsgen[j] >> theCurrentTemp.entx[k][i].xavggen[j] >> theCurrentTemp.entx[k][i].xrmsgen[j];
214 
215  if(in_file.fail()) {LOGERROR("SiPixelGenError") << "Error reading file 30a, no GenError load, run # " << theCurrentTemp.entx[k][i].runnum << ENDL; return false;}
216  }
217 
218  }
219  }
220 
221 
222  in_file.close();
223 
224  // Add this info to the store
225 
226  pixelTemp.push_back(theCurrentTemp);
227 
228  postInit(pixelTemp);
229 
230  return true;
231 
232  } else {
233 
234  // If file didn't open, report this
235 
236  LOGERROR("SiPixelGenError") << "Error opening File" << tempfile << ENDL;
237  return false;
238 
239  }
240 
241 } // TempInit
SiPixelGenErrorEntry::qavg
float qavg
average cluster charge for this set of track angles (now includes threshold effects) ...
Definition: SiPixelGenError.h:35
SiPixelGenErrorHeader::ss50
float ss50
1/2 of the single hit dcol threshold in electrons
Definition: SiPixelGenError.h:72
SiPixelGenErrorEntry::qmin
float qmin
minimum cluster charge for valid hit (keeps 99.9% of simulated hits)
Definition: SiPixelGenError.h:45
mps_fire.i
i
Definition: mps_fire.py:330
SiPixelGenErrorHeader::Bfield
float Bfield
Bfield in Tesla.
Definition: SiPixelGenError.h:75
SiPixelGenErrorEntry::runnum
int runnum
< Basic template entry corresponding to a single set of track angles
Definition: SiPixelGenError.h:32
SiPixelGenErrorEntry::yrmsgen
float yrmsgen[4]
generic algorithm: average y-rms of reconstruction binned in 4 charge bins
Definition: SiPixelGenError.h:48
SiPixelGenErrorEntry::dxone
float dxone
mean offset/correction for one pixel x-clusters
Definition: SiPixelGenError.h:39
LOGINFO
#define LOGINFO(x)
Definition: SiPixelGenError.cc:39
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
SiPixelGenErrorHeader::Vbias
float Vbias
detector bias potential in Volts
Definition: SiPixelGenError.h:68
SiPixelGenErrorEntry::xrmsgen
float xrmsgen[4]
generic algorithm: average x-rms of reconstruction binned in 4 charge bins
Definition: SiPixelGenError.h:50
create_public_lumi_plots.in_file
in_file
Definition: create_public_lumi_plots.py:864
SiPixelGenErrorHeader::NTxx
int NTxx
number of Template x-entries in each slice
Definition: SiPixelGenError.h:61
SiPixelGenErrorHeader::fbin
float fbin[3]
The QBin definitions in Q_clus/Q_avg.
Definition: SiPixelGenError.h:76
SiPixelGenErrorEntry::syone
float syone
rms for one pixel y-clusters
Definition: SiPixelGenError.h:38
SiPixelGenErrorHeader::zsize
float zsize
pixel size (for future use in upgraded geometry)
Definition: SiPixelGenError.h:79
SiPixelGenErrorStore::enty
SiPixelGenErrorEntry enty[60]
60 Barrel y templates spanning cluster lengths from 0px to +18px [28 entries for fpix] ...
Definition: SiPixelGenError.h:90
SiPixelGenErrorHeader::s50
float s50
1/2 of the multihit dcol threshold in electrons
Definition: SiPixelGenError.h:71
SiPixelGenErrorEntry::cotalpha
float cotalpha
cot(alpha) is proportional to cluster length in x and is basis of interpolation
Definition: SiPixelGenError.h:33
SiPixelGenErrorHeader::lorxbias
float lorxbias
estimate of x-lorentz bias
Definition: SiPixelGenError.h:67
SiPixelGenErrorHeader::temperature
float temperature
detector temperature in deg K
Definition: SiPixelGenError.h:69
SiPixelGenErrorHeader::xsize
float xsize
pixel size (for future use in upgraded geometry)
Definition: SiPixelGenError.h:77
SiPixelGenError::postInit
static void postInit(std::vector< SiPixelGenErrorStore > &thePixelTemp_)
Definition: SiPixelGenError.cc:412
SiPixelGenErrorHeader::templ_version
int templ_version
Version number of the template to ensure code compatibility.
Definition: SiPixelGenError.h:74
SiPixelGenErrorEntry::cotbeta
float cotbeta
cot(beta) is proportional to cluster length in y and is basis of interpolation
Definition: SiPixelGenError.h:34
FrontierConditions_GlobalTag_cff.file
file
Definition: FrontierConditions_GlobalTag_cff.py:13
SiPixelGenErrorHeader::ID
int ID
< template header structure
Definition: SiPixelGenError.h:58
SiPixelGenErrorEntry::sxone
float sxone
rms for one pixel x-clusters
Definition: SiPixelGenError.h:40
SiPixelGenErrorHeader::fluence
float fluence
radiation fluence in n_eq/cm^2
Definition: SiPixelGenError.h:70
LOGERROR
#define LOGERROR(x)
Definition: SiPixelGenError.cc:37
SiPixelGenErrorEntry::dytwo
float dytwo
mean offset/correction for one double-pixel y-clusters
Definition: SiPixelGenError.h:41
SiPixelGenErrorHeader::NTy
int NTy
number of Template y entries
Definition: SiPixelGenError.h:59
SiPixelGenErrorHeader::qscale
float qscale
Charge scaling to match cmssw and pixelav.
Definition: SiPixelGenError.h:63
SiPixelGenErrorEntry::dxtwo
float dxtwo
mean offset/correction for one double-pixel x-clusters
Definition: SiPixelGenError.h:43
SiPixelGenErrorStore::head
SiPixelGenErrorHeader head
< template storage structure
Definition: SiPixelGenError.h:85
gen::k
int k[5][pyjets_maxn]
Definition: Cascade2Hadronizer.cc:79
SiPixelGenErrorHeader::ysize
float ysize
pixel size (for future use in upgraded geometry)
Definition: SiPixelGenError.h:78
SiPixelGenErrorHeader::title
char title[80]
template title
Definition: SiPixelGenError.h:73
SiPixelGenErrorEntry::sytwo
float sytwo
rms for one double-pixel y-clusters
Definition: SiPixelGenError.h:42
SiPixelGenErrorHeader::lorywidth
float lorywidth
estimate of y-lorentz width for optimal resolution
Definition: SiPixelGenError.h:64
ENDL
#define ENDL
Definition: SiPixelGenError.cc:40
SiPixelGenErrorEntry::yavggen
float yavggen[4]
generic algorithm: average y-bias of reconstruction binned in 4 charge bins
Definition: SiPixelGenError.h:47
SiPixelGenErrorEntry::dyone
float dyone
mean offset/correction for one pixel y-clusters
Definition: SiPixelGenError.h:37
SiPixelGenErrorHeader::Dtype
int Dtype
detector type (0=BPix, 1=FPix)
Definition: SiPixelGenError.h:62
SiPixelGenErrorHeader::NTyx
int NTyx
number of Template y-slices of x entries
Definition: SiPixelGenError.h:60
dir
dbl *** dir
Definition: mlp_gen.cc:35
SiPixelGenErrorEntry::xavggen
float xavggen[4]
generic algorithm: average x-bias of reconstruction binned in 4 charge bins
Definition: SiPixelGenError.h:49
SiPixelGenErrorEntry::pixmax
float pixmax
maximum charge for individual pixels in cluster
Definition: SiPixelGenError.h:36
SiPixelGenErrorEntry::sxtwo
float sxtwo
rms for one double-pixel x-clusters
Definition: SiPixelGenError.h:44
SiPixelGenErrorStore::entx
SiPixelGenErrorEntry entx[5][29]
29 Barrel x templates spanning cluster lengths from -6px (-1.125Rad) to +6px (+1.125Rad) in each of 5...
Definition: SiPixelGenError.h:91
SiPixelGenErrorHeader::lorybias
float lorybias
estimate of y-lorentz bias
Definition: SiPixelGenError.h:66
SiPixelGenErrorHeader::lorxwidth
float lorxwidth
estimate of x-lorentz width for optimal resolution
Definition: SiPixelGenError.h:65
bool SiPixelGenError::pushfile ( const SiPixelGenErrorDBObject &  dbobject,
std::vector< SiPixelGenErrorStore > &  pixelTemp 
)
static

This routine initializes the global GenError structures from an SiPixelGenErrorDBObject

Parameters
dbobject- db storing multiple generic calibrations

SiPixelGenErrorStore theCurrentTemp; // not on the stack...

Definition at line 251 of file SiPixelGenError.cc.

References SiPixelGenErrorDBObject::char2float::c, ENDL, SiPixelGenErrorDBObject::char2float::f, SiPixelGenErrorDBObject::fail(), mps_fire::i, SiPixelGenErrorDBObject::incrementIndex(), SiPixelGenErrorDBObject::index(), gen::k, LOGERROR, LOGINFO, funct::m, SiPixelGenErrorDBObject::numOfTempl(), SiPixelGenErrorDBObject::sVector(), and groupFilesInBlocks::temp.

252  {
253  // Add GenError stored in external dbobject to theGenErrorStore
254 
255  // Local variables
256  int i, j, k;
257  float costrk[3]={0,0,0};
258  // const char *tempfile;
259  const int code_version={1};
260 
261  // We must create a new object because dbobject must be a const and our stream must not be
262  SiPixelGenErrorDBObject db = dbobject;
263 
264  // Create a local GenError storage entry
266  auto tmpPtr = std::make_unique<SiPixelGenErrorStore>(); // must be allocated on the heap instead
267  auto & theCurrentTemp = *tmpPtr;
268 
269  // Fill the GenError storage for each GenError calibration stored in the db
270  for(int m=0; m<db.numOfTempl(); ++m) {
271 
272  // Read-in a header string first and print it
273 
274  SiPixelGenErrorDBObject::char2float temp;
275  for (i=0; i<20; ++i) {
276  temp.f = db.sVector()[db.index()];
277  theCurrentTemp.head.title[4*i] = temp.c[0];
278  theCurrentTemp.head.title[4*i+1] = temp.c[1];
279  theCurrentTemp.head.title[4*i+2] = temp.c[2];
280  theCurrentTemp.head.title[4*i+3] = temp.c[3];
281  db.incrementIndex(1);
282  }
283 
284  theCurrentTemp.head.title[79] = '\0';
285  LOGINFO("SiPixelGenError") << "Loading Pixel GenError File - "
286  << theCurrentTemp.head.title << ENDL;
287 
288  // next, the header information
289 
290  db >> theCurrentTemp.head.ID >> theCurrentTemp.head.templ_version >> theCurrentTemp.head.Bfield >> theCurrentTemp.head.NTy >> theCurrentTemp.head.NTyx >> theCurrentTemp.head.NTxx >> theCurrentTemp.head.Dtype >> theCurrentTemp.head.Vbias >>
291  theCurrentTemp.head.temperature >> theCurrentTemp.head.fluence >> theCurrentTemp.head.qscale >> theCurrentTemp.head.s50 >>
292  theCurrentTemp.head.lorywidth >> theCurrentTemp.head.lorxwidth >> theCurrentTemp.head.ysize >> theCurrentTemp.head.xsize >>
293  theCurrentTemp.head.zsize >> theCurrentTemp.head.ss50 >> theCurrentTemp.head.lorybias >> theCurrentTemp.head.lorxbias >>
294  theCurrentTemp.head.fbin[0] >> theCurrentTemp.head.fbin[1] >> theCurrentTemp.head.fbin[2];
295 
296  if(db.fail()) {LOGERROR("SiPixelGenError")
297  << "Error reading file, no GenError load" << ENDL; return false;}
298 
299  LOGINFO("SiPixelGenError") << "GenError ID = " << theCurrentTemp.head.ID << ", GenError Version " << theCurrentTemp.head.templ_version << ", Bfield = " << theCurrentTemp.head.Bfield
300  << ", NTy = " << theCurrentTemp.head.NTy << ", NTyx = " << theCurrentTemp.head.NTyx<< ", NTxx = " << theCurrentTemp.head.NTxx << ", Dtype = " << theCurrentTemp.head.Dtype
301  << ", Bias voltage " << theCurrentTemp.head.Vbias << ", temperature "
302  << theCurrentTemp.head.temperature << ", fluence " << theCurrentTemp.head.fluence << ", Q-scaling factor " << theCurrentTemp.head.qscale
303  << ", 1/2 multi dcol threshold " << theCurrentTemp.head.s50 << ", 1/2 single dcol threshold " << theCurrentTemp.head.ss50
304  << ", y Lorentz Width " << theCurrentTemp.head.lorywidth << ", y Lorentz Bias " << theCurrentTemp.head.lorybias
305  << ", x Lorentz width " << theCurrentTemp.head.lorxwidth << ", x Lorentz Bias " << theCurrentTemp.head.lorxbias
306  << ", Q/Q_avg fractions for Qbin defs " << theCurrentTemp.head.fbin[0] << ", " << theCurrentTemp.head.fbin[1]
307  << ", " << theCurrentTemp.head.fbin[2]
308  << ", pixel x-size " << theCurrentTemp.head.xsize << ", y-size " << theCurrentTemp.head.ysize << ", zsize " << theCurrentTemp.head.zsize << ENDL;
309 
310  LOGINFO("SiPixelGenError") << "Loading Pixel GenError - "
311  << theCurrentTemp.head.title << " version "
312  <<theCurrentTemp.head.templ_version <<" code v."
313  <<code_version<<ENDL;
314  if(theCurrentTemp.head.templ_version < code_version) {
315  LOGERROR("SiPixelGenError") << "code expects version " << code_version
316  << ", no GenError load" << ENDL; return false;}
317 
318 #ifdef SI_PIXEL_TEMPLATE_USE_BOOST
319 
320  // next, layout the 1-d/2-d structures needed to store GenError
321 
322  // &&& Who is going to delete these? Are we leaking memory?
323  theCurrentTemp.cotbetaY = new float[theCurrentTemp.head.NTy];
324  theCurrentTemp.cotbetaX = new float[theCurrentTemp.head.NTyx];
325  theCurrentTemp.cotalphaX = new float[theCurrentTemp.head.NTxx];
326 
327  theCurrentTemp.enty.resize(boost::extents[theCurrentTemp.head.NTy]);
328 
329  theCurrentTemp.entx.resize(boost::extents[theCurrentTemp.head.NTyx][theCurrentTemp.head.NTxx]);
330 
331 #endif
332 
333  // next, loop over all barrel y-angle entries
334 
335  for (i=0; i < theCurrentTemp.head.NTy; ++i) {
336 
337  db >> theCurrentTemp.enty[i].runnum >> costrk[0] >> costrk[1] >> costrk[2];
338 
339  if(db.fail()) {LOGERROR("SiPixelGenError") << "Error reading file 1, no GenError load, run # " << theCurrentTemp.enty[i].runnum << ENDL; return false;}
340 
341  // Calculate the alpha, beta, and cot(beta) for this entry
342 
343  theCurrentTemp.enty[i].cotalpha = costrk[0]/costrk[2];
344 
345  theCurrentTemp.enty[i].cotbeta = costrk[1]/costrk[2];
346 
347  db >> theCurrentTemp.enty[i].qavg >> theCurrentTemp.enty[i].pixmax >> theCurrentTemp.enty[i].dyone
348  >> theCurrentTemp.enty[i].syone >> theCurrentTemp.enty[i].dxone >> theCurrentTemp.enty[i].sxone;
349 
350  if(db.fail()) {LOGERROR("SiPixelGenError") << "Error reading file 2, no GenError load, run # " << theCurrentTemp.enty[i].runnum << ENDL; return false;}
351 
352  db >> theCurrentTemp.enty[i].dytwo >> theCurrentTemp.enty[i].sytwo >> theCurrentTemp.enty[i].dxtwo
353  >> theCurrentTemp.enty[i].sxtwo >> theCurrentTemp.enty[i].qmin >> theCurrentTemp.enty[i].qmin2;
354 
355  for (j=0; j<4; ++j) {
356 
357  db >> theCurrentTemp.enty[i].yavggen[j] >> theCurrentTemp.enty[i].yrmsgen[j] >> theCurrentTemp.enty[i].xavggen[j] >> theCurrentTemp.enty[i].xrmsgen[j];
358 
359  if(db.fail()) {LOGERROR("SiPixelGenError") << "Error reading file 14a, no GenError load, run # " << theCurrentTemp.enty[i].runnum << ENDL; return false;}
360  }
361 
362  }
363 
364  // next, loop over all barrel x-angle entries
365 
366  for (k=0; k < theCurrentTemp.head.NTyx; ++k) {
367 
368  for (i=0; i < theCurrentTemp.head.NTxx; ++i) {
369 
370  db >> theCurrentTemp.entx[k][i].runnum >> costrk[0] >> costrk[1] >> costrk[2];
371 
372  if(db.fail()) {LOGERROR("SiPixelGenError") << "Error reading file 17, no GenError load, run # " << theCurrentTemp.entx[k][i].runnum << ENDL; return false;}
373 
374  // Calculate the alpha, beta, and cot(beta) for this entry
375 
376  theCurrentTemp.entx[k][i].cotalpha = costrk[0]/costrk[2];
377 
378  theCurrentTemp.entx[k][i].cotbeta = costrk[1]/costrk[2];
379 
380  db >> theCurrentTemp.entx[k][i].qavg >> theCurrentTemp.entx[k][i].pixmax >> theCurrentTemp.entx[k][i].dyone
381  >> theCurrentTemp.entx[k][i].syone >> theCurrentTemp.entx[k][i].dxone >> theCurrentTemp.entx[k][i].sxone;
382 
383  if(db.fail()) {LOGERROR("SiPixelGenError") << "Error reading file 18, no GenError load, run # " << theCurrentTemp.entx[k][i].runnum << ENDL; return false;}
384 
385  db >> theCurrentTemp.entx[k][i].dytwo >> theCurrentTemp.entx[k][i].sytwo >> theCurrentTemp.entx[k][i].dxtwo
386  >> theCurrentTemp.entx[k][i].sxtwo >> theCurrentTemp.entx[k][i].qmin >> theCurrentTemp.entx[k][i].qmin2;
387 
388  for (j=0; j<4; ++j) {
389 
390  db >> theCurrentTemp.entx[k][i].yavggen[j] >> theCurrentTemp.entx[k][i].yrmsgen[j] >> theCurrentTemp.entx[k][i].xavggen[j] >> theCurrentTemp.entx[k][i].xrmsgen[j];
391 
392  if(db.fail()) {LOGERROR("SiPixelGenError") << "Error reading file 30a, no GenError load, run # " << theCurrentTemp.entx[k][i].runnum << ENDL; return false;}
393  }
394 
395  }
396  }
397 
398 
399  // Add this GenError to the store
400 
401  pixelTemp.push_back(theCurrentTemp);
402 
403  postInit(pixelTemp);
404 
405  }
406  return true;
407 
408 } // TempInit
mps_fire.i
i
Definition: mps_fire.py:330
LOGINFO
#define LOGINFO(x)
Definition: SiPixelGenError.cc:39
SiPixelGenError::postInit
static void postInit(std::vector< SiPixelGenErrorStore > &thePixelTemp_)
Definition: SiPixelGenError.cc:412
SiPixelGenErrorDBObject::sVector
std::vector< float > sVector() const
Definition: SiPixelGenErrorDBObject.h:67
LOGERROR
#define LOGERROR(x)
Definition: SiPixelGenError.cc:37
gen::k
int k[5][pyjets_maxn]
Definition: Cascade2Hadronizer.cc:79
ENDL
#define ENDL
Definition: SiPixelGenError.cc:40
int SiPixelGenError::qbin ( int  id,
float  cotalpha,
float  cotbeta,
float  locBz,
float  locBx,
float  qclus,
bool  irradiationCorrections,
int &  pixmx,
float &  sigmay,
float &  deltay,
float &  sigmax,
float &  deltax,
float &  sy1,
float &  dy1,
float &  sy2,
float &  dy2,
float &  sx1,
float &  dx1,
float &  sx2,
float &  dx2 
)

Definition at line 480 of file SiPixelGenError.cc.

References funct::abs(), gather_cfg::cout, Exception, f, mps_fire::i, triggerObjects_cff::id, createfilelist::int, and mathSSE::sqrt().

Referenced by PixelCPEGeneric::localPosition(), and qbin().

484 {
485  // Interpolate for a new set of track angles
486 
487 
488  // Find the index corresponding to id
489 
490 
491  if(id != id_current_) {
492 
493  index_id_ = -1;
494  for( int i=0; i<(int)thePixelTemp_.size(); ++i) {
495  if(id == thePixelTemp_[i].head.ID) {
496  index_id_ = i;
497  id_current_ = id;
498  lorywidth_ = thePixelTemp_[i].head.lorywidth;
499  lorxwidth_ = thePixelTemp_[i].head.lorxwidth;
500  lorybias_ = thePixelTemp_[i].head.lorybias;
501  lorxbias_ = thePixelTemp_[i].head.lorxbias;
502  for(int j=0; j<3; ++j) {fbin_[j] = thePixelTemp_[i].head.fbin[j];}
503 
504  // Pixel sizes to the private variables
505 
506  xsize_ = thePixelTemp_[i].head.xsize;
507  ysize_ = thePixelTemp_[i].head.ysize;
508  zsize_ = thePixelTemp_[i].head.zsize;
509 
510  break;
511  }
512  }
513  }
514 
515  int index = index_id_;
516 
517 #ifndef SI_PIXEL_TEMPLATE_STANDALONE
518  if(index < 0 || index >= (int)thePixelTemp_.size()) {
519  throw cms::Exception("DataCorrupt") << "SiPixelGenError::qbin can't find needed GenError ID = " << id << std::endl;
520  }
521 #else
522  assert(index >= 0 && index < (int)thePixelTemp_.size());
523 #endif
524 
525  //
526 
527 
528  auto const & templ = thePixelTemp_[index];
529 
530  // Interpolate the absolute value of cot(beta)
531 
532  auto acotb = std::abs(cotbeta);
533 
534  // qcorrect corrects the cot(alpha)=0 cluster charge for non-zero cot(alpha)
535 
536  auto cotalpha0 = thePixelTemp_[index].enty[0].cotalpha;
537  auto qcorrect=std::sqrt((1.f+cotbeta*cotbeta+cotalpha*cotalpha)/(1.f+cotbeta*cotbeta+cotalpha0*cotalpha0));
538 
539  // for some cosmics, the ususal gymnastics are incorrect
540 
541  float cota = cotalpha;
542  float cotb = acotb;
543  bool flip_y; bool flip_x;
544  // for some cosmics, the ususal gymnastics are incorrect
545  flip_x = false;
546  flip_y = false;
547  switch(thePixelTemp_[index_id_].head.Dtype) {
548  case 0:
549  if(cotbeta < 0.f) {flip_y = true;}
550  break;
551  case 1:
552  if(locBz < 0.f) {
553  cotb = cotbeta;
554  } else {
555  cotb = -cotbeta;
556  flip_y = true;
557  }
558  break;
559  case 2:
560  case 3:
561  case 4:
562  case 5:
563  if(locBx*locBz < 0.f) {
564  cota = -cotalpha;
565  flip_x = true;
566  }
567  if(locBx > 0.f) {
568  cotb = cotbeta;
569  } else {
570  cotb = -cotbeta;
571  flip_y = true;
572  }
573  break;
574  default:
575 #ifndef SI_PIXEL_TEMPLATE_STANDALONE
576  throw cms::Exception("DataCorrupt") << "SiPixelGenError::illegal subdetector ID = " << thePixelTemp_[index_id_].head.Dtype << std::endl;
577 #else
578  std::cout << "SiPixelGenError::illegal subdetector ID = " << thePixelTemp_[index_id_].head.Dtype << std::endl;
579 #endif
580  }
581 
582  // Copy the charge scaling factor to the private variable
583 
584  if(flip_y) {
585  lorybias_ = -lorybias_;
586  lorywidth_ = -lorywidth_;
587  }
588  if(flip_x) {
589  lorxbias_ = -lorxbias_;
590  lorxwidth_ = -lorxwidth_;
591  }
592 
593 
594  auto qscale = thePixelTemp_[index].head.qscale;
595 
596 
597  /*
598  lorywidth = thePixelTemp_[index].head.lorywidth;
599  if(locBz > 0.f) {lorywidth = -lorywidth;}
600  lorxwidth = thePixelTemp_[index].head.lorxwidth;
601  */
602 
603 
604  auto Ny = thePixelTemp_[index].head.NTy;
605  auto Nyx = thePixelTemp_[index].head.NTyx;
606  auto Nxx = thePixelTemp_[index].head.NTxx;
607 
608 #ifndef SI_PIXEL_TEMPLATE_STANDALONE
609  if(Ny < 2 || Nyx < 1 || Nxx < 2) {
610  throw cms::Exception("DataCorrupt") << "GenError ID = " << id_current_ << "has too few entries: Ny/Nyx/Nxx = " << Ny << "/" << Nyx << "/" << Nxx << std::endl;
611  }
612 #else
613  assert(Ny > 1 && Nyx > 0 && Nxx > 1);
614 #endif
615 
616  // next, loop over all y-angle entries
617 
618  auto ilow = 0;
619  auto ihigh = 0;
620  auto yratio = 0.f;
621 
622  {
623  auto j = std::lower_bound(templ.cotbetaY,templ.cotbetaY+Ny,cotb);
624  if (j==templ.cotbetaY+Ny) { --j; yratio = 1.f; }
625  else if (j==templ.cotbetaY) { ++j; yratio = 0.f;}
626  else { yratio = (cotb - (*(j-1)) )/( (*j) - (*(j-1)) ) ; }
627 
628  ihigh = j-templ.cotbetaY;
629  ilow = ihigh-1;
630  }
631 
632 
633 
634  // Interpolate/store all y-related quantities (flip displacements when flip_y)
635 
636  auto qavg = (1.f - yratio)*thePixelTemp_[index].enty[ilow].qavg + yratio*thePixelTemp_[index].enty[ihigh].qavg;
637  qavg *= qcorrect;
638  auto qmin = (1.f - yratio)*thePixelTemp_[index].enty[ilow].qmin + yratio*thePixelTemp_[index].enty[ihigh].qmin;
639  qmin *= qcorrect;
640  auto qmin2 = (1.f - yratio)*thePixelTemp_[index].enty[ilow].qmin2 + yratio*thePixelTemp_[index].enty[ihigh].qmin2;
641  qmin2 *= qcorrect;
642 
643 
644 #ifndef SI_PIXEL_TEMPLATE_STANDALONE
645  if(qavg <= 0.f || qmin <= 0.f) {
646  throw cms::Exception("DataCorrupt") << "SiPixelGenError::qbin, qavg or qmin <= 0,"
647  << " Probably someone called the generic pixel reconstruction with an illegal trajectory state" << std::endl;
648  }
649 #else
650  assert(qavg > 0.f && qmin > 0.f);
651 #endif
652 
653  // Scale the input charge to account for differences between pixelav and CMSSW simulation or data
654 
655  auto qtotal = qscale*qclus;
656 
657  // uncertainty and final corrections depend upon total charge bin
658 
659  auto fq = qtotal/qavg;
660  int binq;
661  if(fq > fbin_[0]) {
662  binq=0;
663  } else {
664  if(fq > fbin_[1]) {
665  binq=1;
666  } else {
667  if(fq > fbin_[2]) {
668  binq=2;
669  } else {
670  binq=3;
671  }
672  }
673  }
674 
675  auto yrmsgen =(1.f - yratio)*thePixelTemp_[index].enty[ilow].yrmsgen[binq] + yratio*thePixelTemp_[index].enty[ihigh].yrmsgen[binq];
676  sy1 = (1.f - yratio)*thePixelTemp_[index].enty[ilow].syone + yratio*thePixelTemp_[index].enty[ihigh].syone;
677  sy2 = (1.f - yratio)*thePixelTemp_[index].enty[ilow].sytwo + yratio*thePixelTemp_[index].enty[ihigh].sytwo;
678 
679  if(irradiationCorrections){
680 
681  auto yavggen =(1.f - yratio)*thePixelTemp_[index].enty[ilow].yavggen[binq] + yratio*thePixelTemp_[index].enty[ihigh].yavggen[binq];
682  if(flip_y) {yavggen = -yavggen;}
683  deltay = yavggen;
684  dy1 = (1.f - yratio)*thePixelTemp_[index].enty[ilow].dyone + yratio*thePixelTemp_[index].enty[ihigh].dyone;
685  if(flip_y) {dy1 = -dy1;}
686  dy2 = (1.f - yratio)*thePixelTemp_[index].enty[ilow].dytwo + yratio*thePixelTemp_[index].enty[ihigh].dytwo;
687  if(flip_y) {dy2 = -dy2;}
688 
689  }
690 
691  // next, loop over all x-angle entries, first, find relevant y-slices
692 
693  auto iylow = 0;
694  auto iyhigh = 0;
695  auto yxratio = 0.f;
696 
697 
698  {
699  auto j = std::lower_bound(templ.cotbetaX,templ.cotbetaX+Nyx,acotb);
700  if (j==templ.cotbetaX+Nyx) { --j; yxratio = 1.f; }
701  else if (j==templ.cotbetaX) { ++j; yxratio = 0.f;}
702  else { yxratio = (acotb - (*(j-1)) )/( (*j) - (*(j-1)) ) ; }
703 
704  iyhigh = j-templ.cotbetaX;
705  iylow = iyhigh -1;
706  }
707 
708 
709 
710  auto xxratio = 0.f;
711 
712  {
713  auto j = std::lower_bound(templ.cotalphaX,templ.cotalphaX+Nxx,cota);
714  if (j==templ.cotalphaX+Nxx) { --j; xxratio = 1.f; }
715  else if (j==templ.cotalphaX) { ++j; xxratio = 0.f;}
716  else { xxratio = (cota - (*(j-1)) )/( (*j) - (*(j-1)) ) ; }
717 
718  ihigh = j-templ.cotalphaX;
719  ilow = ihigh-1;
720  }
721 
722 
723 
724  sx1 = (1.f - xxratio)*thePixelTemp_[index].entx[0][ilow].sxone + xxratio*thePixelTemp_[index].entx[0][ihigh].sxone;
725  sx2 = (1.f - xxratio)*thePixelTemp_[index].entx[0][ilow].sxtwo + xxratio*thePixelTemp_[index].entx[0][ihigh].sxtwo;
726 
727  // pixmax is the maximum allowed pixel charge (used for truncation)
728 
729  pixmx=(1.f - yxratio)*((1.f - xxratio)*thePixelTemp_[index].entx[iylow][ilow].pixmax + xxratio*thePixelTemp_[index].entx[iylow][ihigh].pixmax)
730  +yxratio*((1.f - xxratio)*thePixelTemp_[index].entx[iyhigh][ilow].pixmax + xxratio*thePixelTemp_[index].entx[iyhigh][ihigh].pixmax);
731 
732  auto xrmsgen = (1.f - yxratio)*((1.f - xxratio)*thePixelTemp_[index].entx[iylow][ilow].xrmsgen[binq] + xxratio*thePixelTemp_[index].entx[iylow][ihigh].xrmsgen[binq])
733  +yxratio*((1.f - xxratio)*thePixelTemp_[index].entx[iyhigh][ilow].xrmsgen[binq] + xxratio*thePixelTemp_[index].entx[iyhigh][ihigh].xrmsgen[binq]);
734 
735 
736  if(irradiationCorrections) {
737 
738  auto xavggen = (1.f - yxratio)*((1.f - xxratio)*thePixelTemp_[index].entx[iylow][ilow].xavggen[binq] + xxratio*thePixelTemp_[index].entx[iylow][ihigh].xavggen[binq])
739  +yxratio*((1.f - xxratio)*thePixelTemp_[index].entx[iyhigh][ilow].xavggen[binq] + xxratio*thePixelTemp_[index].entx[iyhigh][ihigh].xavggen[binq]);
740  if(flip_x) {xavggen = -xavggen;}
741  deltax = xavggen;
742  dx1 = (1.f - xxratio)*thePixelTemp_[index].entx[0][ilow].dxone + xxratio*thePixelTemp_[index].entx[0][ihigh].dxone;
743  if(flip_x) {dx1 = -dx1;}
744  dx2 = (1.f - xxratio)*thePixelTemp_[index].entx[0][ilow].dxtwo + xxratio*thePixelTemp_[index].entx[0][ihigh].dxtwo;
745  if(flip_x) {dx2 = -dx2;}
746  }
747 
748 
749  // Take the errors and bias from the correct charge bin
750 
751  sigmay = yrmsgen;
752 
753  sigmax = xrmsgen;
754 
755  // If the charge is too small (then flag it)
756 
757  if(qtotal < 0.95f*qmin) {binq = 5;} else {if(qtotal < 0.95f*qmin2) {binq = 4;}}
758 
759  return binq;
760 
761 } // qbin
SiPixelGenError::xsize_
float xsize_
Pixel x-size.
Definition: SiPixelGenError.h:176
SiPixelGenError::lorxbias_
float lorxbias_
Lorentz x-width.
Definition: SiPixelGenError.h:173
mps_fire.i
i
Definition: mps_fire.py:330
SiPixelGenError::lorywidth_
float lorywidth_
Lorentz y-width (sign corrected for fpix frame)
Definition: SiPixelGenError.h:172
SiPixelGenError::ysize_
float ysize_
Pixel y-size.
Definition: SiPixelGenError.h:177
triggerObjects_cff.id
id
Definition: triggerObjects_cff.py:31
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:18
SiPixelGenError::fbin_
float fbin_[3]
The QBin definitions in Q_clus/Q_avg.
Definition: SiPixelGenError.h:175
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
f
double f[11][100]
Definition: MuScleFitUtils.cc:78
SiPixelGenError::zsize_
float zsize_
Pixel z-size (thickness)
Definition: SiPixelGenError.h:178
SiPixelGenError::lorxwidth_
float lorxwidth_
Lorentz x-width.
Definition: SiPixelGenError.h:171
SiPixelGenError::id_current_
int id_current_
current id
Definition: SiPixelGenError.h:165
SiPixelGenError::thePixelTemp_
const std::vector< SiPixelGenErrorStore > & thePixelTemp_
Definition: SiPixelGenError.h:183
gather_cfg.cout
cout
Definition: gather_cfg.py:144
SiPixelGenError::lorybias_
float lorybias_
Lorentz y-width (sign corrected for fpix frame)
Definition: SiPixelGenError.h:174
SiPixelGenError::index_id_
int index_id_
current index
Definition: SiPixelGenError.h:166
int SiPixelGenError::qbin ( int  id,
float  cotalpha,
float  cotbeta,
float  locBz,
float  locBx,
float  qclus,
float &  pixmx,
float &  sigmay,
float &  deltay,
float &  sigmax,
float &  deltax,
float &  sy1,
float &  dy1,
float &  sy2,
float &  dy2,
float &  sx1,
float &  dx1,
float &  sx2,
float &  dx2 
)

Definition at line 763 of file SiPixelGenError.cc.

References objects.autophobj::float, and qbin().

767 {
768  // Interpolate for a new set of track angles
769 
770  bool irradiationCorrections = true;
771  int ipixmx, ibin;
772 
773  ibin = SiPixelGenError::qbin(id, cotalpha, cotbeta, locBz, locBx, qclus, irradiationCorrections,
774  ipixmx, sigmay, deltay, sigmax, deltax, sy1, dy1, sy2, dy2, sx1, dx1, sx2, dx2);
775 
776  pixmx = (float)ipixmx;
777 
778  return ibin;
779 
780 }
SiPixelGenError::qbin
int qbin(int id, float cotalpha, float cotbeta, float locBz, float locBx, float qclus, bool irradiationCorrections, int &pixmx, float &sigmay, float &deltay, float &sigmax, float &deltax, float &sy1, float &dy1, float &sy2, float &dy2, float &sx1, float &dx1, float &sx2, float &dx2)
Definition: SiPixelGenError.cc:480
int SiPixelGenError::qbin ( int  id)

Interpolate beta/alpha angles to produce an expected average charge. Return int (0-4) describing the charge of the cluster [0: 1.5<Q/Qavg, 1: 1<Q/Qavg<1.5, 2: 0.85<Q/Qavg<1, 3: 0.95Qmin<Q<0.85Qavg, 4: Q<0.95Qmin].

Parameters
id- (input) index of the GenError to use
cotalpha- (input) the cotangent of the alpha track angle (see CMS IN 2004/014)
cotbeta- (input) the cotangent of the beta track angle (see CMS IN 2004/014)
locBz- (input) the sign of this quantity is used to determine whether to flip cot(beta)<0 quantities from cot(beta)>0 (FPix only) for FPix IP-related tracks, locBz < 0 for cot(beta) > 0 and locBz > 0 for cot(beta) < 0
qclus- (input) the cluster charge in electrons
pixmax- (output) the maximum pixel charge in electrons (truncation value)
sigmay- (output) the estimated y-error for CPEGeneric in microns
deltay- (output) the estimated y-bias for CPEGeneric in microns
sigmax- (output) the estimated x-error for CPEGeneric in microns
deltax- (output) the estimated x-bias for CPEGeneric in microns
sy1- (output) the estimated y-error for 1 single-pixel clusters in microns
dy1- (output) the estimated y-bias for 1 single-pixel clusters in microns
sy2- (output) the estimated y-error for 1 double-pixel clusters in microns
dy2- (output) the estimated y-bias for 1 double-pixel clusters in microns
sx1- (output) the estimated x-error for 1 single-pixel clusters in microns
dx1- (output) the estimated x-bias for 1 single-pixel clusters in microns
sx2- (output) the estimated x-error for 1 double-pixel clusters in microns
dx2- (output) the estimated x-bias for 1 double-pixel clusters in microns

Definition at line 449 of file SiPixelGenError.cc.

References mps_fire::i, triggerObjects_cff::id, and createfilelist::int.

449  {
450  // Find the index corresponding to id
451 
452  if(id != id_current_) {
453 
454  index_id_ = -1;
455  for( int i=0; i<(int)thePixelTemp_.size(); ++i) {
456  if(id == thePixelTemp_[i].head.ID) {
457  index_id_ = i;
458  id_current_ = id;
459  //
460  lorywidth_ = thePixelTemp_[i].head.lorywidth;
461  lorxwidth_ = thePixelTemp_[i].head.lorxwidth;
462  lorybias_ = thePixelTemp_[i].head.lorybias;
463  lorxbias_ = thePixelTemp_[i].head.lorxbias;
464 
465  //for(int j=0; j<3; ++j) {fbin_[j] = thePixelTemp_[i].head.fbin[j];}
466 
467  // Pixel sizes to the private variables
468  xsize_ = thePixelTemp_[i].head.xsize;
469  ysize_ = thePixelTemp_[i].head.ysize;
470  zsize_ = thePixelTemp_[i].head.zsize;
471 
472  break;
473  }
474  }
475  }
476  return index_id_;
477 }
SiPixelGenError::xsize_
float xsize_
Pixel x-size.
Definition: SiPixelGenError.h:176
SiPixelGenError::lorxbias_
float lorxbias_
Lorentz x-width.
Definition: SiPixelGenError.h:173
mps_fire.i
i
Definition: mps_fire.py:330
SiPixelGenError::lorywidth_
float lorywidth_
Lorentz y-width (sign corrected for fpix frame)
Definition: SiPixelGenError.h:172
SiPixelGenError::ysize_
float ysize_
Pixel y-size.
Definition: SiPixelGenError.h:177
triggerObjects_cff.id
id
Definition: triggerObjects_cff.py:31
SiPixelGenError::zsize_
float zsize_
Pixel z-size (thickness)
Definition: SiPixelGenError.h:178
SiPixelGenError::lorxwidth_
float lorxwidth_
Lorentz x-width.
Definition: SiPixelGenError.h:171
SiPixelGenError::id_current_
int id_current_
current id
Definition: SiPixelGenError.h:165
SiPixelGenError::thePixelTemp_
const std::vector< SiPixelGenErrorStore > & thePixelTemp_
Definition: SiPixelGenError.h:183
SiPixelGenError::lorybias_
float lorybias_
Lorentz y-width (sign corrected for fpix frame)
Definition: SiPixelGenError.h:174
SiPixelGenError::index_id_
int index_id_
current index
Definition: SiPixelGenError.h:166
float SiPixelGenError::xsize ( )
inline

pixel x-size (microns)

Definition at line 157 of file SiPixelGenError.h.

float SiPixelGenError::ysize ( )
inline

pixel y-size (microns)

Definition at line 158 of file SiPixelGenError.h.

float SiPixelGenError::zsize ( )
inline

pixel z-size or thickness (microns)

Definition at line 159 of file SiPixelGenError.h.

Member Data Documentation

float SiPixelGenError::fbin_[3]
private

The QBin definitions in Q_clus/Q_avg.

Definition at line 175 of file SiPixelGenError.h.

int SiPixelGenError::id_current_
private

current id

Definition at line 165 of file SiPixelGenError.h.

int SiPixelGenError::index_id_
private

current index

Definition at line 166 of file SiPixelGenError.h.

float SiPixelGenError::lorxbias_
private

Lorentz x-width.

Definition at line 173 of file SiPixelGenError.h.

float SiPixelGenError::lorxwidth_
private

Lorentz x-width.

Definition at line 171 of file SiPixelGenError.h.

float SiPixelGenError::lorybias_
private

Lorentz y-width (sign corrected for fpix frame)

Definition at line 174 of file SiPixelGenError.h.

float SiPixelGenError::lorywidth_
private

Lorentz y-width (sign corrected for fpix frame)

Definition at line 172 of file SiPixelGenError.h.

const std::vector< SiPixelGenErrorStore >& SiPixelGenError::thePixelTemp_
private

Definition at line 183 of file SiPixelGenError.h.

float SiPixelGenError::xsize_
private

Pixel x-size.

Definition at line 176 of file SiPixelGenError.h.

float SiPixelGenError::ysize_
private

Pixel y-size.

Definition at line 177 of file SiPixelGenError.h.

float SiPixelGenError::zsize_
private

Pixel z-size (thickness)

Definition at line 178 of file SiPixelGenError.h.

Generated for CMSSW Reference Manual by doxygen 1.8.11