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 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 > &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 102 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 104 of file SiPixelGenError.h.

References id_current_, and index_id_.

Member Function Documentation

float SiPixelGenError::fbin ( int  i )

inline

Return lower bound of Qbin definition.

Definition at line 128 of file SiPixelGenError.h.

References assert(), edm::hlt::Exception, fbin_, and i.

float SiPixelGenError::lorxbias ( )

inline

signed lorentz x-bias (microns)

Definition at line 126 of file SiPixelGenError.h.

References lorxbias_.

float SiPixelGenError::lorxwidth ( )

inline

signed lorentz x-width (microns)

Definition at line 124 of file SiPixelGenError.h.

References lorxwidth_.

Referenced by PixelCPEGeneric::localPosition().

float SiPixelGenError::lorybias ( )

inline

signed lorentz y-bias (microns)

Definition at line 125 of file SiPixelGenError.h.

References lorybias_.

float SiPixelGenError::lorywidth ( )

inline

signed lorentz y-width (microns)

Definition at line 123 of file SiPixelGenError.h.

References lorywidth_.

Referenced by PixelCPEGeneric::localPosition().

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

static

Definition at line 395 of file SiPixelGenError.cc.

References citk::for().

                                                                                {
   
   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 52 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, mergeVDriftHistosByStation::file, SiPixelGenErrorHeader::fluence, SiPixelGenErrorStore::head, i, SiPixelGenErrorHeader::ID, recoMuon::in, create_public_lumi_plots::in_file, j, relval_steps::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 241 of file SiPixelGenError.cc.

References SiPixelGenErrorHeader::Bfield, SiPixelGenErrorDBObject::char2float::c, SiPixelGenErrorEntry::cotalpha, SiPixelGenErrorEntry::cotbeta, EcalCondDB::db, 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, i, SiPixelGenErrorHeader::ID, SiPixelGenErrorDBObject::incrementIndex(), SiPixelGenErrorDBObject::index(), j, relval_steps::k, LOGERROR, LOGINFO, SiPixelGenErrorHeader::lorxbias, SiPixelGenErrorHeader::lorxwidth, SiPixelGenErrorHeader::lorybias, SiPixelGenErrorHeader::lorywidth, visualization-live-secondInstance_cfg::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	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 463 of file SiPixelGenError.cc.

References funct::abs(), assert(), edm::hlt::Exception, f, i, cmsHarvester::index, j, 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 cotb; bool flip_y;
   if(thePixelTemp_[index].head.Dtype == 0) {
      cotb = acotb;
      flip_y = false;
      if(cotbeta < 0.f) {flip_y = true;}
   } else {
      if(locBz < 0.f) {
         cotb = cotbeta;
         flip_y = false;
      } else {
         cotb = -cotbeta;
         flip_y = true;
      }
   }
   
   // Copy the charge scaling factor to the private variable
   
   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)
   
   dy1 = (1.f - yratio)*thePixelTemp_[index].enty[ilow].dyone + yratio*thePixelTemp_[index].enty[ihigh].dyone;
   if(flip_y) {dy1 = -dy1;}
   sy1 = (1.f - yratio)*thePixelTemp_[index].enty[ilow].syone + yratio*thePixelTemp_[index].enty[ihigh].syone;
   dy2 = (1.f - yratio)*thePixelTemp_[index].enty[ilow].dytwo + yratio*thePixelTemp_[index].enty[ihigh].dytwo;
   if(flip_y) {dy2 = -dy2;}
   sy2 = (1.f - yratio)*thePixelTemp_[index].enty[ilow].sytwo + yratio*thePixelTemp_[index].enty[ihigh].sytwo;
   
   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 yavggen =(1.f - yratio)*thePixelTemp_[index].enty[ilow].yavggen[binq] + yratio*thePixelTemp_[index].enty[ihigh].yavggen[binq];
   if(flip_y) {yavggen = -yavggen;}
   auto yrmsgen =(1.f - yratio)*thePixelTemp_[index].enty[ilow].yrmsgen[binq] + yratio*thePixelTemp_[index].enty[ihigh].yrmsgen[binq];
   
   
   // 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;
   }
   
   
   
   ilow = ihigh = 0;
   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 = (cotalpha - (*(j-1)) )/( (*j) - (*(j-1)) ) ; }
      
      ihigh = j-templ.cotalphaX;
      ilow = ihigh-1;
   }
   
   
   
   dx1 = (1.f - xxratio)*thePixelTemp_[index].entx[0][ilow].dxone + xxratio*thePixelTemp_[index].entx[0][ihigh].dxone;
   sx1 = (1.f - xxratio)*thePixelTemp_[index].entx[0][ilow].sxone + xxratio*thePixelTemp_[index].entx[0][ihigh].sxone;
   dx2 = (1.f - xxratio)*thePixelTemp_[index].entx[0][ilow].dxtwo + xxratio*thePixelTemp_[index].entx[0][ihigh].dxtwo;
   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 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]);
   
   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]);
   
   
   
   //  Take the errors and bias from the correct charge bin
    
   sigmay = yrmsgen; deltay = yavggen;
   
   sigmax = xrmsgen; deltax = xavggen;
   
   // 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 432 of file SiPixelGenError.cc.

References i.

                                {
   // 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 135 of file SiPixelGenError.h.

References xsize_.

float SiPixelGenError::ysize ( )

inline

pixel y-size (microns)

Definition at line 136 of file SiPixelGenError.h.

References ysize_.

float SiPixelGenError::zsize ( )

inline

pixel z-size or thickness (microns)

Definition at line 137 of file SiPixelGenError.h.

References zsize_.

Member Data Documentation

float SiPixelGenError::fbin_[3]

private

The QBin definitions in Q_clus/Q_avg.

Definition at line 153 of file SiPixelGenError.h.

Referenced by fbin().

int SiPixelGenError::id_current_

private

current id

Definition at line 143 of file SiPixelGenError.h.

Referenced by SiPixelGenError().

int SiPixelGenError::index_id_

private

current index

Definition at line 144 of file SiPixelGenError.h.

Referenced by SiPixelGenError().

float SiPixelGenError::lorxbias_

private

Lorentz x-width.

Definition at line 151 of file SiPixelGenError.h.

Referenced by lorxbias().

float SiPixelGenError::lorxwidth_

private

Lorentz x-width.

Definition at line 149 of file SiPixelGenError.h.

Referenced by lorxwidth().

float SiPixelGenError::lorybias_

private

Lorentz y-width (sign corrected for fpix frame)

Definition at line 152 of file SiPixelGenError.h.

Referenced by lorybias().

float SiPixelGenError::lorywidth_

private

Lorentz y-width (sign corrected for fpix frame)

Definition at line 150 of file SiPixelGenError.h.

Referenced by lorywidth().

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

private

Definition at line 161 of file SiPixelGenError.h.

float SiPixelGenError::xsize_

private

Pixel x-size.

Definition at line 154 of file SiPixelGenError.h.

Referenced by xsize().

float SiPixelGenError::ysize_

private

Pixel y-size.

Definition at line 155 of file SiPixelGenError.h.

Referenced by ysize().

float SiPixelGenError::zsize_

private

Pixel z-size (thickness)

Definition at line 156 of file SiPixelGenError.h.

Referenced by zsize().