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)
	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 > &thePixelTemp_)
static bool	pushfile (const SiPixelGenErrorDBObject &dbobject, std::vector< SiPixelGenErrorStore > &thePixelTemp_)

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 103 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 105 of file SiPixelGenError.h.

References SiPixelGenErrorEntry::cotalpha, and SiPixelGenErrorEntry::cotbeta.

Member Function Documentation

float SiPixelGenError::fbin ( int  i )

inline

Return lower bound of Qbin definition.

Definition at line 130 of file SiPixelGenError.h.

References Exception, and mps_fire::i.

float SiPixelGenError::lorxbias ( )

inline

signed lorentz x-bias (microns)

Definition at line 128 of file SiPixelGenError.h.

float SiPixelGenError::lorxwidth ( )

inline

signed lorentz x-width (microns)

Definition at line 126 of file SiPixelGenError.h.

Referenced by PixelCPEGeneric::localPosition().

float SiPixelGenError::lorybias ( )

inline

signed lorentz y-bias (microns)

Definition at line 127 of file SiPixelGenError.h.

float SiPixelGenError::lorywidth ( )

inline

signed lorentz y-width (microns)

Definition at line 125 of file SiPixelGenError.h.

Referenced by PixelCPEGeneric::localPosition().

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

static

Definition at line 396 of file SiPixelGenError.cc.

                                                                                {
   
   for (auto & templ : thePixelTemp_) {
      for ( auto iy=0; iy<templ.head.NTy; ++iy ) templ.cotbetaY[iy]=templ.enty[iy].cotbeta;
      for ( auto iy=0; iy<templ.head.NTyx; ++iy )  templ.cotbetaX[iy]= templ.entx[iy][0].cotbeta;
      for ( auto ix=0; ix<templ.head.NTxx; ++ix ) templ.cotalphaX[ix]=templ.entx[0][ix].cotalpha;
   }
   
}

bool SiPixelGenError::pushfile ( int  filenum, std::vector< SiPixelGenErrorStore > &  thePixelTemp_  )

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 53 of file SiPixelGenError.cc.

References SiPixelGenErrorHeader::Bfield, EnergyCorrector::c, SiPixelGenErrorEntry::cotalpha, SiPixelGenErrorEntry::cotbeta, 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().

{
   // Add info stored in external file numbered filenum to theGenErrorStore
   
   // Local variables
   int i, j, k;
   float costrk[3]={0,0,0};
   const char *tempfile;
   //   char title[80]; remove this
   char c;
   const int code_version={1};
   
   
   
   //  Create a filename for this run
   
   std::ostringstream tout;
   
   //  Create different path in CMSSW than standalone
   
#ifndef SI_PIXEL_TEMPLATE_STANDALONE
   tout << "CalibTracker/SiPixelESProducers/data/generror_summary_zp"
   << std::setw(4) << std::setfill('0') << std::right << filenum << ".out" << std::ends;
   std::string tempf = tout.str();
   edm::FileInPath file( tempf.c_str() );
   tempfile = (file.fullPath()).c_str();
#else
   tout << "generror_summary_zp" << std::setw(4) << std::setfill('0') << std::right << filenum << ".out" << std::ends;
   std::string tempf = tout.str();
   tempfile = tempf.c_str();
#endif
   
   //  open the Generic file
   
   std::ifstream in_file(tempfile, std::ios::in);
   
   if(in_file.is_open()) {
      
      // Create a local GenError storage entry
      
      SiPixelGenErrorStore theCurrentTemp;
      
      // Read-in a header string first and print it
      
      for (i=0; (c=in_file.get()) != '\n'; ++i) {
         if(i < 79) {theCurrentTemp.head.title[i] = c;}
      }
      if(i > 78) {i=78;}
      theCurrentTemp.head.title[i+1] ='\0';
      LOGINFO("SiPixelGenError") << "Loading Pixel GenError File - " << theCurrentTemp.head.title << ENDL;
      
      // next, the header information
      
      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 >>
      theCurrentTemp.head.temperature >> theCurrentTemp.head.fluence >> theCurrentTemp.head.qscale >> theCurrentTemp.head.s50 >>
      theCurrentTemp.head.lorywidth >> theCurrentTemp.head.lorxwidth >> theCurrentTemp.head.ysize >> theCurrentTemp.head.xsize >>
      theCurrentTemp.head.zsize >> theCurrentTemp.head.ss50 >> theCurrentTemp.head.lorybias >> theCurrentTemp.head.lorxbias >>
      theCurrentTemp.head.fbin[0] >> theCurrentTemp.head.fbin[1] >> theCurrentTemp.head.fbin[2];
      
      if(in_file.fail()) {LOGERROR("SiPixelGenError") << "Error reading file, no GenError load" << ENDL; return false;}
      
      LOGINFO("SiPixelGenError") << "GenError ID = " << theCurrentTemp.head.ID << ", GenError Version " << theCurrentTemp.head.templ_version << ", Bfield = " << theCurrentTemp.head.Bfield
      << ", NTy = " << theCurrentTemp.head.NTy << ", NTyx = " << theCurrentTemp.head.NTyx<< ", NTxx = " << theCurrentTemp.head.NTxx << ", Dtype = " << theCurrentTemp.head.Dtype
      << ", Bias voltage " << theCurrentTemp.head.Vbias << ", temperature "
      << theCurrentTemp.head.temperature << ", fluence " << theCurrentTemp.head.fluence << ", Q-scaling factor " << theCurrentTemp.head.qscale
      << ", 1/2 multi dcol threshold " << theCurrentTemp.head.s50 << ", 1/2 single dcol threshold " << theCurrentTemp.head.ss50
      << ", y Lorentz Width " << theCurrentTemp.head.lorywidth << ", y Lorentz Bias " << theCurrentTemp.head.lorybias
      << ", x Lorentz width " << theCurrentTemp.head.lorxwidth << ", x Lorentz Bias " << theCurrentTemp.head.lorxbias
      << ", Q/Q_avg fractions for Qbin defs " << theCurrentTemp.head.fbin[0] << ", " << theCurrentTemp.head.fbin[1]
      << ", " << theCurrentTemp.head.fbin[2]
      << ", pixel x-size " << theCurrentTemp.head.xsize << ", y-size " << theCurrentTemp.head.ysize << ", zsize " << theCurrentTemp.head.zsize << ENDL;
      
      if(theCurrentTemp.head.templ_version < code_version) {LOGERROR("SiPixelGenError") << "code expects version " << code_version << ", no GenError load" << ENDL; return false;}
      
#ifdef SI_PIXEL_TEMPLATE_USE_BOOST
      
      // next, layout the 1-d/2-d structures needed to store GenError info
      
      theCurrentTemp.enty.resize(boost::extents[theCurrentTemp.head.NTy]);
      
      theCurrentTemp.entx.resize(boost::extents[theCurrentTemp.head.NTyx][theCurrentTemp.head.NTxx]);
      
#endif
      
      // next, loop over all y-angle entries
      
      for (i=0; i < theCurrentTemp.head.NTy; ++i) {
         
         in_file >> theCurrentTemp.enty[i].runnum >> costrk[0] >> costrk[1] >> costrk[2];
         
         if(in_file.fail()) {LOGERROR("SiPixelGenError") << "Error reading file 1, no GenError load, run # " << theCurrentTemp.enty[i].runnum << ENDL; return false;}
         
         // Calculate the alpha, beta, and cot(beta) for this entry
         
         theCurrentTemp.enty[i].cotalpha = costrk[0]/costrk[2];
         
         theCurrentTemp.enty[i].cotbeta = costrk[1]/costrk[2];
         
         in_file >> theCurrentTemp.enty[i].qavg >> theCurrentTemp.enty[i].pixmax >> theCurrentTemp.enty[i].dyone
         >> theCurrentTemp.enty[i].syone >> theCurrentTemp.enty[i].dxone >> theCurrentTemp.enty[i].sxone;
         
         if(in_file.fail()) {LOGERROR("SiPixelGenError") << "Error reading file 2, no GenError load, run # " << theCurrentTemp.enty[i].runnum << ENDL; return false;}
         
         in_file >> theCurrentTemp.enty[i].dytwo >> theCurrentTemp.enty[i].sytwo >> theCurrentTemp.enty[i].dxtwo
         >> theCurrentTemp.enty[i].sxtwo >> theCurrentTemp.enty[i].qmin >>  theCurrentTemp.enty[i].qmin2;
         
         if(in_file.fail()) {LOGERROR("SiPixelGenError") << "Error reading file 3, no GenError load, run # " << theCurrentTemp.enty[i].runnum << ENDL; return false;}
         
         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;}
         
         
         for (j=0; j<4; ++j) {
            
            in_file >> theCurrentTemp.enty[i].yavggen[j] >> theCurrentTemp.enty[i].yrmsgen[j] >> theCurrentTemp.enty[i].xavggen[j] >> theCurrentTemp.enty[i].xrmsgen[j];
            
            if(in_file.fail()) {LOGERROR("SiPixelGenError") << "Error reading file 14a, no GenError load, run # " << theCurrentTemp.enty[i].runnum << ENDL; return false;}
         }
         
         
      }
      
      // next, loop over all barrel x-angle entries
      
      for (k=0; k < theCurrentTemp.head.NTyx; ++k) {
         
         for (i=0; i < theCurrentTemp.head.NTxx; ++i) {
            
            in_file >> theCurrentTemp.entx[k][i].runnum >> costrk[0] >> costrk[1] >> costrk[2];
            
            if(in_file.fail()) {LOGERROR("SiPixelGenError") << "Error reading file 17, no GenError load, run # " << theCurrentTemp.entx[k][i].runnum << ENDL; return false;}
            
            // Calculate the alpha, beta, and cot(beta) for this entry
            
            theCurrentTemp.entx[k][i].cotalpha = costrk[0]/costrk[2];
            
            theCurrentTemp.entx[k][i].cotbeta = costrk[1]/costrk[2];
            
            in_file >> theCurrentTemp.entx[k][i].qavg >> theCurrentTemp.entx[k][i].pixmax >> theCurrentTemp.entx[k][i].dyone
            >> theCurrentTemp.entx[k][i].syone >> theCurrentTemp.entx[k][i].dxone >> theCurrentTemp.entx[k][i].sxone;
            
            if(in_file.fail()) {LOGERROR("SiPixelGenError") << "Error reading file 18, no GenError load, run # " << theCurrentTemp.entx[k][i].runnum << ENDL; return false;}
            
            in_file >> theCurrentTemp.entx[k][i].dytwo >> theCurrentTemp.entx[k][i].sytwo >> theCurrentTemp.entx[k][i].dxtwo
            >> theCurrentTemp.entx[k][i].sxtwo >> theCurrentTemp.entx[k][i].qmin >> theCurrentTemp.entx[k][i].qmin2;
            //             >> theCurrentTemp.entx[k][i].mpvvav >> theCurrentTemp.entx[k][i].sigmavav >> theCurrentTemp.entx[k][i].kappavav;
            
            if(in_file.fail()) {LOGERROR("SiPixelGenError") << "Error reading file 19, no GenError load, run # " << theCurrentTemp.entx[k][i].runnum << ENDL; return false;}
            
            
            for (j=0; j<4; ++j) {
               
               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];
               
               if(in_file.fail()) {LOGERROR("SiPixelGenError") << "Error reading file 30a, no GenError load, run # " << theCurrentTemp.entx[k][i].runnum << ENDL; return false;}
            }
                                                
         }
      }
      
      
      in_file.close();
      
      // Add this info to the store
      
      thePixelTemp_.push_back(theCurrentTemp);
      
      postInit(thePixelTemp_);
      
      return true;
      
   } else {
      
      // If file didn't open, report this
      
      LOGERROR("SiPixelGenError") << "Error opening File" << tempfile << ENDL;
      return false;
      
   }
   
} // TempInit

bool SiPixelGenError::pushfile ( const SiPixelGenErrorDBObject &  dbobject, std::vector< SiPixelGenErrorStore > &  thePixelTemp_  )

static

This routine initializes the global GenError structures from an SiPixelGenErrorDBObject

Parameters

dbobject

- db storing multiple generic calibrations

Definition at line 242 of file SiPixelGenError.cc.

References SiPixelGenErrorHeader::Bfield, SiPixelGenErrorDBObject::char2float::c, SiPixelGenErrorEntry::cotalpha, SiPixelGenErrorEntry::cotbeta, SiPixelGenErrorHeader::Dtype, SiPixelGenErrorEntry::dxone, SiPixelGenErrorEntry::dxtwo, SiPixelGenErrorEntry::dyone, SiPixelGenErrorEntry::dytwo, ENDL, SiPixelGenErrorStore::entx, SiPixelGenErrorStore::enty, SiPixelGenErrorDBObject::char2float::f, SiPixelGenErrorDBObject::fail(), SiPixelGenErrorHeader::fbin, SiPixelGenErrorHeader::fluence, SiPixelGenErrorStore::head, mps_fire::i, SiPixelGenErrorHeader::ID, SiPixelGenErrorDBObject::incrementIndex(), SiPixelGenErrorDBObject::index(), gen::k, LOGERROR, LOGINFO, SiPixelGenErrorHeader::lorxbias, SiPixelGenErrorHeader::lorxwidth, SiPixelGenErrorHeader::lorybias, SiPixelGenErrorHeader::lorywidth, funct::m, SiPixelGenErrorHeader::NTxx, SiPixelGenErrorHeader::NTy, SiPixelGenErrorHeader::NTyx, SiPixelGenErrorDBObject::numOfTempl(), SiPixelGenErrorEntry::pixmax, SiPixelGenErrorEntry::qavg, SiPixelGenErrorEntry::qmin, SiPixelGenErrorEntry::qmin2, SiPixelGenErrorHeader::qscale, SiPixelGenErrorEntry::runnum, SiPixelGenErrorHeader::s50, SiPixelGenErrorHeader::ss50, SiPixelGenErrorDBObject::sVector(), SiPixelGenErrorEntry::sxone, SiPixelGenErrorEntry::sxtwo, SiPixelGenErrorEntry::syone, SiPixelGenErrorEntry::sytwo, groupFilesInBlocks::temp, SiPixelGenErrorHeader::temperature, SiPixelGenErrorHeader::templ_version, SiPixelGenErrorHeader::title, SiPixelGenErrorHeader::Vbias, SiPixelGenErrorEntry::xavggen, SiPixelGenErrorEntry::xrmsgen, SiPixelGenErrorHeader::xsize, SiPixelGenErrorEntry::yavggen, SiPixelGenErrorEntry::yrmsgen, SiPixelGenErrorHeader::ysize, and SiPixelGenErrorHeader::zsize.

                                                                                  {
   // Add GenError stored in external dbobject to theGenErrorStore
   
   // Local variables
   int i, j, k;
   float costrk[3]={0,0,0};
   //   const char *tempfile;
   const int code_version={1};
   
   // We must create a new object because dbobject must be a const and our stream must not be
   SiPixelGenErrorDBObject db = dbobject;
   
   // Create a local GenError storage entry
   SiPixelGenErrorStore theCurrentTemp;
   
   // Fill the GenError storage for each GenError calibration stored in the db
   for(int m=0; m<db.numOfTempl(); ++m) {
      
      // Read-in a header string first and print it
      
      SiPixelGenErrorDBObject::char2float temp;
      for (i=0; i<20; ++i) {
         temp.f = db.sVector()[db.index()];
         theCurrentTemp.head.title[4*i] = temp.c[0];
         theCurrentTemp.head.title[4*i+1] = temp.c[1];
         theCurrentTemp.head.title[4*i+2] = temp.c[2];
         theCurrentTemp.head.title[4*i+3] = temp.c[3];
         db.incrementIndex(1);
      }
      
      theCurrentTemp.head.title[79] = '\0';
      LOGINFO("SiPixelGenError") << "Loading Pixel GenError File - "
      << theCurrentTemp.head.title << ENDL;
      
      // next, the header information
      
      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 >>
      theCurrentTemp.head.temperature >> theCurrentTemp.head.fluence >> theCurrentTemp.head.qscale >> theCurrentTemp.head.s50 >>
      theCurrentTemp.head.lorywidth >> theCurrentTemp.head.lorxwidth >> theCurrentTemp.head.ysize >> theCurrentTemp.head.xsize >>
      theCurrentTemp.head.zsize >> theCurrentTemp.head.ss50 >> theCurrentTemp.head.lorybias >> theCurrentTemp.head.lorxbias >>
      theCurrentTemp.head.fbin[0] >> theCurrentTemp.head.fbin[1] >> theCurrentTemp.head.fbin[2];
      
      if(db.fail()) {LOGERROR("SiPixelGenError")
         << "Error reading file, no GenError load" << ENDL; return false;}
      
      LOGINFO("SiPixelGenError") << "GenError ID = " << theCurrentTemp.head.ID << ", GenError Version " << theCurrentTemp.head.templ_version << ", Bfield = " << theCurrentTemp.head.Bfield
      << ", NTy = " << theCurrentTemp.head.NTy << ", NTyx = " << theCurrentTemp.head.NTyx<< ", NTxx = " << theCurrentTemp.head.NTxx << ", Dtype = " << theCurrentTemp.head.Dtype
      << ", Bias voltage " << theCurrentTemp.head.Vbias << ", temperature "
      << theCurrentTemp.head.temperature << ", fluence " << theCurrentTemp.head.fluence << ", Q-scaling factor " << theCurrentTemp.head.qscale
      << ", 1/2 multi dcol threshold " << theCurrentTemp.head.s50 << ", 1/2 single dcol threshold " << theCurrentTemp.head.ss50
      << ", y Lorentz Width " << theCurrentTemp.head.lorywidth << ", y Lorentz Bias " << theCurrentTemp.head.lorybias
      << ", x Lorentz width " << theCurrentTemp.head.lorxwidth << ", x Lorentz Bias " << theCurrentTemp.head.lorxbias
      << ", Q/Q_avg fractions for Qbin defs " << theCurrentTemp.head.fbin[0] << ", " << theCurrentTemp.head.fbin[1]
      << ", " << theCurrentTemp.head.fbin[2]
      << ", pixel x-size " << theCurrentTemp.head.xsize << ", y-size " << theCurrentTemp.head.ysize << ", zsize " << theCurrentTemp.head.zsize << ENDL;
      
      LOGINFO("SiPixelGenError") << "Loading Pixel GenError - "
                  << theCurrentTemp.head.title << " version "
                  <<theCurrentTemp.head.templ_version <<" code v."
                  <<code_version<<ENDL;
      if(theCurrentTemp.head.templ_version < code_version) {
         LOGERROR("SiPixelGenError") << "code expects version " << code_version
         << ", no GenError load" << ENDL; return false;}
            
#ifdef SI_PIXEL_TEMPLATE_USE_BOOST
      
      // next, layout the 1-d/2-d structures needed to store GenError
      
      theCurrentTemp.enty.resize(boost::extents[theCurrentTemp.head.NTy]);
      
      theCurrentTemp.entx.resize(boost::extents[theCurrentTemp.head.NTyx][theCurrentTemp.head.NTxx]);
      
#endif
      
      // next, loop over all barrel y-angle entries
      
      for (i=0; i < theCurrentTemp.head.NTy; ++i) {
         
         db >> theCurrentTemp.enty[i].runnum >> costrk[0] >> costrk[1] >> costrk[2];
         
         if(db.fail()) {LOGERROR("SiPixelGenError") << "Error reading file 1, no GenError load, run # " << theCurrentTemp.enty[i].runnum << ENDL; return false;}
         
         // Calculate the alpha, beta, and cot(beta) for this entry
         
         theCurrentTemp.enty[i].cotalpha = costrk[0]/costrk[2];
         
         theCurrentTemp.enty[i].cotbeta = costrk[1]/costrk[2];
         
         db >> theCurrentTemp.enty[i].qavg >> theCurrentTemp.enty[i].pixmax >> theCurrentTemp.enty[i].dyone
         >> theCurrentTemp.enty[i].syone >> theCurrentTemp.enty[i].dxone >> theCurrentTemp.enty[i].sxone;
         
         if(db.fail()) {LOGERROR("SiPixelGenError") << "Error reading file 2, no GenError load, run # " << theCurrentTemp.enty[i].runnum << ENDL; return false;}
         
         db >> theCurrentTemp.enty[i].dytwo >> theCurrentTemp.enty[i].sytwo >> theCurrentTemp.enty[i].dxtwo
         >> theCurrentTemp.enty[i].sxtwo >> theCurrentTemp.enty[i].qmin >> theCurrentTemp.enty[i].qmin2;
         
         for (j=0; j<4; ++j) {
            
            db >> theCurrentTemp.enty[i].yavggen[j] >> theCurrentTemp.enty[i].yrmsgen[j] >> theCurrentTemp.enty[i].xavggen[j] >> theCurrentTemp.enty[i].xrmsgen[j];
            
            if(db.fail()) {LOGERROR("SiPixelGenError") << "Error reading file 14a, no GenError load, run # " << theCurrentTemp.enty[i].runnum << ENDL; return false;}
         }
         
      }
      
      // next, loop over all barrel x-angle entries
      
      for (k=0; k < theCurrentTemp.head.NTyx; ++k) {
         
         for (i=0; i < theCurrentTemp.head.NTxx; ++i) {
            
            db >> theCurrentTemp.entx[k][i].runnum >> costrk[0] >> costrk[1] >> costrk[2];
            
            if(db.fail()) {LOGERROR("SiPixelGenError") << "Error reading file 17, no GenError load, run # " << theCurrentTemp.entx[k][i].runnum << ENDL; return false;}
            
            // Calculate the alpha, beta, and cot(beta) for this entry
            
            theCurrentTemp.entx[k][i].cotalpha = costrk[0]/costrk[2];
            
            theCurrentTemp.entx[k][i].cotbeta = costrk[1]/costrk[2];
            
            db >> theCurrentTemp.entx[k][i].qavg >> theCurrentTemp.entx[k][i].pixmax >> theCurrentTemp.entx[k][i].dyone
            >> theCurrentTemp.entx[k][i].syone >> theCurrentTemp.entx[k][i].dxone >> theCurrentTemp.entx[k][i].sxone;
            
            if(db.fail()) {LOGERROR("SiPixelGenError") << "Error reading file 18, no GenError load, run # " << theCurrentTemp.entx[k][i].runnum << ENDL; return false;}
            
            db >> theCurrentTemp.entx[k][i].dytwo >> theCurrentTemp.entx[k][i].sytwo >> theCurrentTemp.entx[k][i].dxtwo
            >> theCurrentTemp.entx[k][i].sxtwo >> theCurrentTemp.entx[k][i].qmin >> theCurrentTemp.entx[k][i].qmin2;
            
            for (j=0; j<4; ++j) {
               
               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];
               
               if(db.fail()) {LOGERROR("SiPixelGenError") << "Error reading file 30a, no GenError load, run # " << theCurrentTemp.entx[k][i].runnum << ENDL; return false;}
            }
            
         }
      }
      
      
      // Add this GenError to the store
      
      thePixelTemp_.push_back(theCurrentTemp);
      
      postInit(thePixelTemp_);
      
   }
   return true;
   
} // TempInit

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 464 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().

{
   // Interpolate for a new set of track angles
   
   
   // Find the index corresponding to id
   
   
   if(id != id_current_) {
      
      index_id_ = -1;
      for( int i=0; i<(int)thePixelTemp_.size(); ++i) {
         if(id == thePixelTemp_[i].head.ID) {
            index_id_ = i;
            id_current_ = id;
            lorywidth_ = thePixelTemp_[i].head.lorywidth;
            lorxwidth_ = thePixelTemp_[i].head.lorxwidth;
            lorybias_ = thePixelTemp_[i].head.lorybias;
            lorxbias_ = thePixelTemp_[i].head.lorxbias;
            for(int j=0; j<3; ++j) {fbin_[j] = thePixelTemp_[i].head.fbin[j];}
            
            // Pixel sizes to the private variables
            
            xsize_ = thePixelTemp_[i].head.xsize;
            ysize_ = thePixelTemp_[i].head.ysize;
            zsize_ = thePixelTemp_[i].head.zsize;
            
            break;
         }
      }
   }
   
   int index = index_id_;
   
#ifndef SI_PIXEL_TEMPLATE_STANDALONE
   if(index < 0 || index >= (int)thePixelTemp_.size()) {
      throw cms::Exception("DataCorrupt") << "SiPixelGenError::qbin can't find needed GenError ID = " << id << std::endl;
   }
#else
   assert(index >= 0 && index < (int)thePixelTemp_.size());
#endif
   
   //
   
   
   auto const & templ = thePixelTemp_[index];
   
   // Interpolate the absolute value of cot(beta)
   
   auto acotb = std::abs(cotbeta);
   
   //   qcorrect corrects the cot(alpha)=0 cluster charge for non-zero cot(alpha)
   
   auto cotalpha0 =  thePixelTemp_[index].enty[0].cotalpha;
   auto qcorrect=std::sqrt((1.f+cotbeta*cotbeta+cotalpha*cotalpha)/(1.f+cotbeta*cotbeta+cotalpha0*cotalpha0));
   
   // for some cosmics, the ususal gymnastics are incorrect
   
   float cota = cotalpha;
   float cotb = acotb;
   bool flip_y; bool flip_x;
   // for some cosmics, the ususal gymnastics are incorrect
   flip_x = false;
   flip_y = false;
   switch(thePixelTemp_[index_id_].head.Dtype) {
      case 0:
         if(cotbeta < 0.f) {flip_y = true;}
         break;
      case 1:
         if(locBz < 0.f) {
            cotb = cotbeta;
         } else {
            cotb = -cotbeta;
            flip_y = true;
         }
         break;
      case 2:
      case 3:
      case 4:
         if(locBx*locBz < 0.f) {
            cota = -cotalpha;
            flip_x = true;
         }
         if(locBx > 0.f) {
            cotb = cotbeta;
         } else {
            cotb = -cotbeta;
            flip_y = true;
         }
         break;
      default:
#ifndef SI_PIXEL_TEMPLATE_STANDALONE
         throw cms::Exception("DataCorrupt") << "SiPixelTemplate::illegal subdetector ID = " << thePixelTemp_[index_id_].head.Dtype << std::endl;
#else
         std::cout << "SiPixelTemplate::illegal subdetector ID = " << thePixelTemp_[index_id_].head.Dtype << std::endl;
#endif
   }
   
   // Copy the charge scaling factor to the private variable
   
   if(flip_y) {
      lorybias_ = -lorybias_;
      lorywidth_ = -lorywidth_;
   }
   if(flip_x) {
      lorxbias_ = -lorxbias_;
      lorxwidth_ = -lorxwidth_;
   }
   
   
   auto qscale = thePixelTemp_[index].head.qscale;
   
   
   /*
    lorywidth = thePixelTemp_[index].head.lorywidth;
    if(locBz > 0.f) {lorywidth = -lorywidth;}
    lorxwidth = thePixelTemp_[index].head.lorxwidth;
    */
   
   
   auto Ny = thePixelTemp_[index].head.NTy;
   auto Nyx = thePixelTemp_[index].head.NTyx;
   auto Nxx = thePixelTemp_[index].head.NTxx;
   
#ifndef SI_PIXEL_TEMPLATE_STANDALONE
   if(Ny < 2 || Nyx < 1 || Nxx < 2) {
      throw cms::Exception("DataCorrupt") << "GenError ID = " << id_current_ << "has too few entries: Ny/Nyx/Nxx = " << Ny << "/" << Nyx << "/" << Nxx << std::endl;
   }
#else
   assert(Ny > 1 && Nyx > 0 && Nxx > 1);
#endif
   
   // next, loop over all y-angle entries
   
   auto ilow = 0;
   auto ihigh = 0;
   auto yratio = 0.f;
   
   {
      auto j = std::lower_bound(templ.cotbetaY,templ.cotbetaY+Ny,cotb);
      if (j==templ.cotbetaY+Ny) { --j;  yratio = 1.f; }
      else if (j==templ.cotbetaY) { ++j; yratio = 0.f;}
      else { yratio = (cotb - (*(j-1)) )/( (*j) - (*(j-1)) ) ; }
      
      ihigh = j-templ.cotbetaY;
      ilow = ihigh-1;
   }
   
   
   
   // Interpolate/store all y-related quantities (flip displacements when flip_y)
   
   auto qavg = (1.f - yratio)*thePixelTemp_[index].enty[ilow].qavg + yratio*thePixelTemp_[index].enty[ihigh].qavg;
   qavg *= qcorrect;
   auto qmin = (1.f - yratio)*thePixelTemp_[index].enty[ilow].qmin + yratio*thePixelTemp_[index].enty[ihigh].qmin;
   qmin *= qcorrect;
   auto qmin2 = (1.f - yratio)*thePixelTemp_[index].enty[ilow].qmin2 + yratio*thePixelTemp_[index].enty[ihigh].qmin2;
   qmin2 *= qcorrect;
   
   
#ifndef SI_PIXEL_TEMPLATE_STANDALONE
   if(qavg <= 0.f || qmin <= 0.f) {
      throw cms::Exception("DataCorrupt") << "SiPixelGenError::qbin, qavg or qmin <= 0,"
      << " Probably someone called the generic pixel reconstruction with an illegal trajectory state" << std::endl;
   }
#else
   assert(qavg > 0.f && qmin > 0.f);
#endif
   
   //  Scale the input charge to account for differences between pixelav and CMSSW simulation or data
   
   auto qtotal = qscale*qclus;
   
   // uncertainty and final corrections depend upon total charge bin
   
   auto fq = qtotal/qavg;
   int  binq;
   if(fq > fbin_[0]) {
      binq=0;
   } else {
      if(fq > fbin_[1]) {
         binq=1;
      } else {
         if(fq > fbin_[2]) {
            binq=2;
         } else {
            binq=3;
         }
      }
   }
   
   auto yrmsgen =(1.f - yratio)*thePixelTemp_[index].enty[ilow].yrmsgen[binq] + yratio*thePixelTemp_[index].enty[ihigh].yrmsgen[binq];
   sy1 = (1.f - yratio)*thePixelTemp_[index].enty[ilow].syone + yratio*thePixelTemp_[index].enty[ihigh].syone;
   sy2 = (1.f - yratio)*thePixelTemp_[index].enty[ilow].sytwo + yratio*thePixelTemp_[index].enty[ihigh].sytwo;
   
   
   // next, loop over all x-angle entries, first, find relevant y-slices
   
   auto iylow = 0;
   auto iyhigh = 0;
   auto yxratio = 0.f;
   
   
   {
      auto j = std::lower_bound(templ.cotbetaX,templ.cotbetaX+Nyx,acotb);
      if (j==templ.cotbetaX+Nyx) { --j;  yxratio = 1.f; }
      else if (j==templ.cotbetaX) { ++j; yxratio = 0.f;}
      else { yxratio = (acotb - (*(j-1)) )/( (*j) - (*(j-1)) ) ; }
      
      iyhigh = j-templ.cotbetaX;
      iylow = iyhigh -1;
   }
   
   
   
   auto xxratio = 0.f;
   
   {
      auto j = std::lower_bound(templ.cotalphaX,templ.cotalphaX+Nxx,cotalpha);
      if (j==templ.cotalphaX+Nxx) { --j;  xxratio = 1.f; }
      else if (j==templ.cotalphaX) { ++j; xxratio = 0.f;}
      else { xxratio = (cota - (*(j-1)) )/( (*j) - (*(j-1)) ) ; }
      
      ihigh = j-templ.cotalphaX;
      ilow = ihigh-1;
   }
   
   
   
   sx1 = (1.f - xxratio)*thePixelTemp_[index].entx[0][ilow].sxone + xxratio*thePixelTemp_[index].entx[0][ihigh].sxone;
   sx2 = (1.f - xxratio)*thePixelTemp_[index].entx[0][ilow].sxtwo + xxratio*thePixelTemp_[index].entx[0][ihigh].sxtwo;
   
   // pixmax is the maximum allowed pixel charge (used for truncation)
   
   pixmx=(1.f - yxratio)*((1.f - xxratio)*thePixelTemp_[index].entx[iylow][ilow].pixmax + xxratio*thePixelTemp_[index].entx[iylow][ihigh].pixmax)
   +yxratio*((1.f - xxratio)*thePixelTemp_[index].entx[iyhigh][ilow].pixmax + xxratio*thePixelTemp_[index].entx[iyhigh][ihigh].pixmax);
   
   auto xrmsgen = (1.f - yxratio)*((1.f - xxratio)*thePixelTemp_[index].entx[iylow][ilow].xrmsgen[binq] + xxratio*thePixelTemp_[index].entx[iylow][ihigh].xrmsgen[binq])
   +yxratio*((1.f - xxratio)*thePixelTemp_[index].entx[iyhigh][ilow].xrmsgen[binq] + xxratio*thePixelTemp_[index].entx[iyhigh][ihigh].xrmsgen[binq]);
   
   
   if(irradiationCorrections) {
      auto yavggen =(1.f - yratio)*thePixelTemp_[index].enty[ilow].yavggen[binq] + yratio*thePixelTemp_[index].enty[ihigh].yavggen[binq];
      if(flip_y) {yavggen = -yavggen;}
      deltay = yavggen;
      dy1 = (1.f - yratio)*thePixelTemp_[index].enty[ilow].dyone + yratio*thePixelTemp_[index].enty[ihigh].dyone;
      if(flip_y) {dy1 = -dy1;}
      dy2 = (1.f - yratio)*thePixelTemp_[index].enty[ilow].dytwo + yratio*thePixelTemp_[index].enty[ihigh].dytwo;
      if(flip_y) {dy2 = -dy2;}

      auto xavggen = (1.f - yxratio)*((1.f - xxratio)*thePixelTemp_[index].entx[iylow][ilow].xavggen[binq] + xxratio*thePixelTemp_[index].entx[iylow][ihigh].xavggen[binq])
      +yxratio*((1.f - xxratio)*thePixelTemp_[index].entx[iyhigh][ilow].xavggen[binq] + xxratio*thePixelTemp_[index].entx[iyhigh][ihigh].xavggen[binq]);
      if(flip_x) {xavggen = -xavggen;}
      deltax = xavggen;
      dx1 = (1.f - xxratio)*thePixelTemp_[index].entx[0][ilow].dxone + xxratio*thePixelTemp_[index].entx[0][ihigh].dxone;
      if(flip_x) {dx1 = -dx1;}
      dx2 = (1.f - xxratio)*thePixelTemp_[index].entx[0][ilow].dxtwo + xxratio*thePixelTemp_[index].entx[0][ihigh].dxtwo;
      if(flip_x) {dx2 = -dx2;}
   }

   
   //  Take the errors and bias from the correct charge bin
   
   sigmay = yrmsgen;
   
   sigmax = xrmsgen;
   
   // If the charge is too small (then flag it)
   
   if(qtotal < 0.95f*qmin) {binq = 5;} else {if(qtotal < 0.95f*qmin2) {binq = 4;}}
   
   return binq;
   
} // qbin

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 433 of file SiPixelGenError.cc.

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

                                {
   // Find the index corresponding to id
   
   if(id != id_current_) {
      
      index_id_ = -1;
      for( int i=0; i<(int)thePixelTemp_.size(); ++i) {
         if(id == thePixelTemp_[i].head.ID) {
            index_id_ = i;
            id_current_ = id;
            //
            lorywidth_ = thePixelTemp_[i].head.lorywidth;
            lorxwidth_ = thePixelTemp_[i].head.lorxwidth;
            lorybias_ = thePixelTemp_[i].head.lorybias;
            lorxbias_ = thePixelTemp_[i].head.lorxbias;
            
            //for(int j=0; j<3; ++j) {fbin_[j] = thePixelTemp_[i].head.fbin[j];}
            
            // Pixel sizes to the private variables
            xsize_ = thePixelTemp_[i].head.xsize;
            ysize_ = thePixelTemp_[i].head.ysize;
            zsize_ = thePixelTemp_[i].head.zsize;
            
            break;
         }
      }
   }
   return index_id_;
}

float SiPixelGenError::xsize ( )

inline

pixel x-size (microns)

Definition at line 137 of file SiPixelGenError.h.

float SiPixelGenError::ysize ( )

inline

pixel y-size (microns)

Definition at line 138 of file SiPixelGenError.h.

float SiPixelGenError::zsize ( )

inline

pixel z-size or thickness (microns)

Definition at line 139 of file SiPixelGenError.h.

Member Data Documentation

float SiPixelGenError::fbin_[3]

private

The QBin definitions in Q_clus/Q_avg.

Definition at line 155 of file SiPixelGenError.h.

int SiPixelGenError::id_current_

private

current id

Definition at line 145 of file SiPixelGenError.h.

int SiPixelGenError::index_id_

private

current index

Definition at line 146 of file SiPixelGenError.h.

float SiPixelGenError::lorxbias_

private

Lorentz x-width.

Definition at line 153 of file SiPixelGenError.h.

float SiPixelGenError::lorxwidth_

private

Lorentz x-width.

Definition at line 151 of file SiPixelGenError.h.

float SiPixelGenError::lorybias_

private

Lorentz y-width (sign corrected for fpix frame)

Definition at line 154 of file SiPixelGenError.h.

float SiPixelGenError::lorywidth_

private

Lorentz y-width (sign corrected for fpix frame)

Definition at line 152 of file SiPixelGenError.h.

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

private

Definition at line 163 of file SiPixelGenError.h.

float SiPixelGenError::xsize_

private

Pixel x-size.

Definition at line 156 of file SiPixelGenError.h.

float SiPixelGenError::ysize_

private

Pixel y-size.

Definition at line 157 of file SiPixelGenError.h.

float SiPixelGenError::zsize_

private

Pixel z-size (thickness)

Definition at line 158 of file SiPixelGenError.h.