SiPixelTemplate2D Class Reference

#include <SiPixelTemplate2D.h>

Public Member Functions
float	chi2avgone ()
	average y chi^2 for 1 pixel clusters More...
float	chi2minone ()
	minimum of y chi^2 for 1 pixel clusters More...
float	chi2ppix ()
	average chi^2 per struck pixel More...
float	chi2scale ()
	scale factor for chi^2 distribution More...
float	clslenx ()
	cluster x-size More...
float	clsleny ()
	cluster y-size More...
float	delyavg ()
	average difference between clsleny_ and cluster length [with threshold effects] More...
float	delysig ()
	rms difference between clsleny_ and cluster length [with threshold effects] More...
float	fbin (int i)
	Return lower bound of Qbin definition. More...
bool	getid (int id)
bool	interpolate (int id, float cotalpha, float cotbeta, float locBz, float locBx)
float	kappavav ()
	kappa parameter in Vavilov distribution More...
void	landau_par (float lanpar[2][5])
	Return the Landau probability parameters for this set of cot(alpha, cot(beta) More...
float	lorxdrift ()
	signed lorentz x-width (microns) More...
float	lorydrift ()
	signed lorentz y-width (microns) More...
float	mpvvav ()
	most probable Q in Vavilov distribution More...
float	offsetx (int i)
	x-offset in 4 charge bins More...
float	offsety (int i)
	y-offset in 4 charge bins More...
float	pixmax ()
	maximum pixel charge More...
float	qavg ()
	average cluster charge for this set of track angles More...
float	qscale ()
	charge scaling factor More...
float	s50 ()
	1/2 of the pixel threshold signal in adc units More...
float	scalex (int i)
	x-error scale factor in 4 charge bins More...
float	scalexavg ()
	x-reco error scaling factor More...
float	scaley (int i)
	y-error scale factor in 4 charge bins More...
float	scaleyavg ()
	y-reco error scaling factor More...
void	sideload (SiPixelTemplateEntry2D *entry, int iDtype, float locBx, float locBz, float lorwdy, float lorwdx, float q50, float fbin[3], float xsize, float ysize, float zsize)
float	sigmavav ()
	scale factor in Vavilov distribution More...
	SiPixelTemplate2D (const std::vector< SiPixelTemplateStore2D > &thePixelTemp)
	Default constructor. More...
float	sizex ()
	return x size of template cluster More...
float	sizey ()
	return y size of template cluster More...
int	storesize ()
	return the size of the template store (the number of stored IDs More...
float	sxymax ()
	max pixel signal for pixel error calculation More...
float	xsize ()
	pixel x-size (microns) More...
void	xysigma2 (float qpixel, int index, float &xysig2)
bool	xytemp (float xhit, float yhit, bool ydouble[21+2], bool xdouble[13+2], float template2d[13+2][21+2], bool dervatives, float dpdx2d[2][13+2][21+2], float &QTemplate)
bool	xytemp (float xhit, float yhit, bool ydouble[21+2], bool xdouble[13+2], float template2d[13+2][21+2])
bool	xytemp (int id, float cotalpha, float cotbeta, float xhit, float yhit, std::vector< bool > &ydouble, std::vector< bool > &xdouble, float template2d[13+2][21+2])
float	ysize ()
	pixel y-size (microns) More...
float	zsize ()
	pixel z-size or thickness (microns) More...

Static Public Member Functions
static bool	pushfile (int filenum, std::vector< SiPixelTemplateStore2D > &pixelTemp, std::string dir="")
static bool	pushfile (const SiPixel2DTemplateDBObject &dbobject, std::vector< SiPixelTemplateStore2D > &pixelTemp)

Private Member Functions
bool	checkIllegalIndex (const std::string whichMethod, int indMax, int i)

Private Attributes
float	adcota_
	fractional pixel distance of cot(alpha) from jx0_ More...
float	adcotb_
	fractional pixel distance of cot(beta) from iy0_ More...
float	chi2avgone_
	average chi^2 for 1 pixel clusters More...
float	chi2minone_
	minimum of chi^2 for 1 pixel clusters More...
float	chi2ppix_
	average chi^2 per struck pixel More...
float	chi2scale_
	scale factor for chi2 distribution More...
float	clslenx_
	projected x-length of cluster More...
float	clsleny_
	projected y-length of cluster More...
float	cota_current_
	current cot alpha More...
float	cotalpha0_
	minimum cot(alpha) covered More...
float	cotalpha1_
	maximum cot(alpha) covered More...
float	cotb_current_
	current cot beta More...
float	cotbeta0_
	minimum cot(beta) covered More...
float	cotbeta1_
	maximum cot(beta) covered More...
float	deltacota_
	cot(alpha) bin size More...
float	deltacotb_
	cot(beta) bin size More...
float	delyavg_
	average difference between clsleny_ and cluster length [with threshold effects] More...
float	delysig_
	rms of difference between clsleny_ and cluster length [with threshold effects] More...
int	Dtype_
	flags BPix (=0) or FPix (=1) More...
const SiPixelTemplateEntry2D *	entry00_
const SiPixelTemplateEntry2D *	entry01_
const SiPixelTemplateEntry2D *	entry10_
float	fbin_ [3]
	The QBin definitions in Q_clus/Q_avg. More...
bool	flip_x_
	flip x sign-sensitive quantities More...
bool	flip_y_
	flip y sign-sensitive quantities More...
int	id_current_
	current id More...
int	imax_
	max y index of templated cluster More...
int	imin_
	min y index of templated cluster More...
int	index_id_
	current index More...
int	iy0_
	index of nearest cot(beta) bin More...
int	iy1_
	index of next-nearest cot(beta) bin More...
int	jmax_
	max x index of templated cluster More...
int	jmin_
	min x index of templated cluster More...
int	jx0_
	index of nearest cot(alpha) bin More...
int	jx1_
	index of next-nearest cot(alpha) bin More...
float	kappavav_
	kappa parameter in Vavilov distribution More...
float	lanpar_ [2][5]
	Interpolated Landau parameters. More...
float	lorxdrift_
	Lorentz x-drift. More...
float	lorxwidth_
	Lorentz x-width. More...
float	lorydrift_
	Lorentz y-drift. More...
float	lorywidth_
	Lorentz y-width (sign corrected for fpix frame) More...
float	mpvvav_
	most probable Q in Vavilov distribution More...
int	Nxx_
	number of cot(alpha)-entries (rows) in template More...
int	Nyx_
	number of cot(beta)-entries (columns) in template More...
float	offsetx_ [4]
	x-offset in charge bins More...
float	offsety_ [4]
	y-offset in charge bins More...
float	pixmax_
	maximum pixel charge More...
float	qavg_
	average cluster charge for this set of track angles More...
float	qscale_
	charge scaling factor More...
float	s50_
	1/2 of the pixel threshold signal in adc units More...
float	scalex_ [4]
	x-error scale factor in charge bins More...
float	scalexavg_
	average x-error scale factor More...
float	scaley_ [4]
	y-error scale factor in charge bins More...
float	scaleyavg_
	average y-error scale factor More...
float	sigmavav_
	scale factor in Vavilov distribution More...
bool	success_
	true if cotalpha, cotbeta are inside of the acceptance (dynamically loaded) More...
float	sxymax_
	average pixel signal for y-projection of cluster More...
const std::vector< SiPixelTemplateStore2D > &	thePixelTemp_
float	xsize_
	Pixel x-size. More...
float	xypary0x0_ [2][5]
	Polynomial error parameterization at ix0,iy0. More...
float	xypary0x1_ [2][5]
	Polynomial error parameterization at ix1,iy0. More...
float	xypary1x0_ [2][5]
	Polynomial error parameterization at ix0,iy1. More...
float	xytemp_ [13+2][21+2]
	template for xy-reconstruction More...
float	ysize_
	Pixel y-size. More...
float	zsize_
	Pixel z-size (thickness) More...

Detailed Description

A template management class. SiPixelTemplate contains thePixelTemp (a std::vector of SiPixelTemplateStore, each of which is a collection of many SiPixelTemplateEntries). Each SiPixelTemplateStore 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 SiPixelTemplate is xytemp(), which produces a template on the fly, given a specific track's alpha and beta. The results are kept in data members and accessed via inline getters.

The resulting template is then used by PixelTempReco2D() (a global function) which get the reference for SiPixelTemplate & templ and uses the current template to reconstruct the SiPixelRecHit.

Definition at line 141 of file SiPixelTemplate2D.h.

Constructor & Destructor Documentation

SiPixelTemplate2D::SiPixelTemplate2D ( const std::vector< SiPixelTemplateStore2D > &  thePixelTemp )

inline

Default constructor.

Definition at line 143 of file SiPixelTemplate2D.h.

References BXM2, BYM2, SiPixelTemplateEntry2D::cotalpha, SiPixelTemplateEntry2D::cotbeta, dir, mps_splice::entry, HiCaloJetParameters_cff::interpolate, SiPixelTemplateEntry2D::lanpar, AlCaHLTBitMon_QueryRunRegistry::string, xsize, SiPixelTemplateEntry2D::xytemp, and ysize.

Member Function Documentation

bool SiPixelTemplate2D::checkIllegalIndex ( const std::string  whichMethod, int  indMax, int  i  )

inlineprivate

Definition at line 215 of file SiPixelTemplate2D.h.

References Exception.

      {
#ifndef SI_PIXEL_TEMPLATE_STANDALONE
         if(i < 0 || i > indMax) {throw cms::Exception("DataCorrupt") << "SiPixelTemplate2D::"<< whichMethod <<" called with illegal index = " << i << std::endl;}
#else
         assert(i>=0 && i<indMax+1);

#endif
         return true;
   }

float SiPixelTemplate2D::chi2avgone ( )

inline

average y chi^2 for 1 pixel clusters

Definition at line 195 of file SiPixelTemplate2D.h.

float SiPixelTemplate2D::chi2minone ( )

inline

minimum of y chi^2 for 1 pixel clusters

Definition at line 196 of file SiPixelTemplate2D.h.

float SiPixelTemplate2D::chi2ppix ( )

inline

average chi^2 per struck pixel

Definition at line 193 of file SiPixelTemplate2D.h.

Referenced by SiPixelTemplateReco2D::PixelTempReco2D().

float SiPixelTemplate2D::chi2scale ( )

inline

scale factor for chi^2 distribution

Definition at line 194 of file SiPixelTemplate2D.h.

Referenced by SiPixelTemplateReco2D::PixelTempReco2D().

float SiPixelTemplate2D::clslenx ( )

inline

cluster x-size

Definition at line 203 of file SiPixelTemplate2D.h.

float SiPixelTemplate2D::clsleny ( )

inline

cluster y-size

Definition at line 202 of file SiPixelTemplate2D.h.

Referenced by SiPixelTemplateReco2D::PixelTempReco2D().

float SiPixelTemplate2D::delyavg ( )

inline

average difference between clsleny_ and cluster length [with threshold effects]

Definition at line 206 of file SiPixelTemplate2D.h.

float SiPixelTemplate2D::delysig ( )

inline

rms difference between clsleny_ and cluster length [with threshold effects]

Definition at line 207 of file SiPixelTemplate2D.h.

float SiPixelTemplate2D::fbin ( int  i )

inline

Return lower bound of Qbin definition.

Definition at line 190 of file SiPixelTemplate2D.h.

References mps_fire::i.

Referenced by SiPixelTemplateReco2D::PixelTempReco2D().

bool SiPixelTemplate2D::getid

(

int

id

)

Definition at line 522 of file SiPixelTemplate2D.cc.

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

Referenced by SiPixelDigitizerAlgorithm::PixelTempRewgt2D().

{
   if(id != id_current_) {
      
      // Find the index corresponding to id
      
      index_id_ = -1;
      for(int i=0; i<(int)thePixelTemp_.size(); ++i) {
         
         if(id == thePixelTemp_[i].head.ID) {
            
            index_id_ = i;
            id_current_ = id;
            
            // Copy the charge scaling factor to the private variable
            
            Dtype_ = thePixelTemp_[index_id_].head.Dtype;
            
            // Copy the charge scaling factor to the private variable
            
            qscale_ = thePixelTemp_[index_id_].head.qscale;
            
            // Copy the pseudopixel signal size to the private variable
            
            s50_ = thePixelTemp_[index_id_].head.s50;
            
            // Copy Qbinning info to private variables
            
            for(int j=0; j<3; ++j) {fbin_[j] = thePixelTemp_[index_id_].head.fbin[j];}
            
            // Copy the Lorentz widths to private variables
            
            lorywidth_ = thePixelTemp_[index_id_].head.lorywidth;
            lorxwidth_ = thePixelTemp_[index_id_].head.lorxwidth;
            
            // Copy the pixel sizes private variables
            
            xsize_ = thePixelTemp_[index_id_].head.xsize;
            ysize_ = thePixelTemp_[index_id_].head.ysize;
            zsize_ = thePixelTemp_[index_id_].head.zsize;
            
            // Determine the size of this template
            
            Nyx_ = thePixelTemp_[index_id_].head.NTyx;
            Nxx_ = thePixelTemp_[index_id_].head.NTxx;
#ifndef SI_PIXEL_TEMPLATE_STANDALONE
            if(Nyx_ < 2 || Nxx_ < 2) {
               throw cms::Exception("DataCorrupt") << "template ID = " << id_current_ << "has too few entries: Nyx/Nxx = " << Nyx_ << "/" << Nxx_ << std::endl;
            }
#else
            assert(Nyx_ > 1 && Nxx_ > 1);
#endif
            int imidx = Nxx_/2;
            
            cotalpha0_ =  thePixelTemp_[index_id_].entry[0][0].cotalpha;
            cotalpha1_ =  thePixelTemp_[index_id_].entry[0][Nxx_-1].cotalpha;
            deltacota_ = (cotalpha1_-cotalpha0_)/(float)(Nxx_-1);
            
            cotbeta0_ =  thePixelTemp_[index_id_].entry[0][imidx].cotbeta;
            cotbeta1_ =  thePixelTemp_[index_id_].entry[Nyx_-1][imidx].cotbeta;
            deltacotb_ = (cotbeta1_-cotbeta0_)/(float)(Nyx_-1);
            
            break;
         }
      }
   }

#ifndef SI_PIXEL_TEMPLATE_STANDALONE
   if(index_id_ < 0 || index_id_ >= (int)thePixelTemp_.size()) {
      throw cms::Exception("DataCorrupt") << "SiPixelTemplate2D::interpolate can't find needed template ID = " << id
      << ", Are you using the correct global tag?" << std::endl;
   }
#else
   assert(index_id_ >= 0 && index_id_ < (int)thePixelTemp_.size());
#endif
   return true;
}

bool SiPixelTemplate2D::interpolate	(	int	id,
		float	cotalpha,
		float	cotbeta,
		float	locBz,
		float	locBx
	)

Interpolate stored 2-D information for input angles

Parameters

id	- (input) the id of the template
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 Phase 0 FPix IP-related tracks, locBz < 0 for cot(beta) > 0 and locBz > 0 for cot(beta) < 0 for Phase 1 FPix IP-related tracks, see next comment
locBx	- (input) the sign of this quantity is used to determine whether to flip cot(alpha/beta)<0 quantities from cot(alpha/beta)>0 (FPix only) for Phase 1 FPix IP-related tracks, locBx/locBz > 0 for cot(alpha) > 0 and locBx/locBz < 0 for cot(alpha) < 0 for Phase 1 FPix IP-related tracks, locBx > 0 for cot(beta) > 0 and locBx < 0 for cot(beta) < 0

Definition at line 613 of file SiPixelTemplate2D.cc.

References gather_cfg::cout, Exception, f, objects.autophobj::float, mps_fire::i, and createfilelist::int.

Referenced by SiPixelTemplateReco2D::PixelTempReco2D(), and xytemp().

{
   
   
   // Interpolate for a new set of track angles
   
   // Local variables

   float acotb, dcota, dcotb;
   
   // Check to see if interpolation is valid
   
   if(id != id_current_ || cotalpha != cota_current_ || cotbeta != cotb_current_) {
      cota_current_ = cotalpha; cotb_current_ = cotbeta;
      success_ = getid(id);
   }
   
#ifndef SI_PIXEL_TEMPLATE_STANDALONE
   if(index_id_ < 0 || index_id_ >= (int)thePixelTemp_.size()) {
      throw cms::Exception("DataCorrupt") << "SiPixelTemplate2D::interpolate can't find needed template ID = " << id
      << ", Are you using the correct global tag?" << std::endl;
   }
#else
   assert(index_id_ >= 0 && index_id_ < (int)thePixelTemp_.size());
#endif
   
   // Check angle limits and et up interpolation parameters
   
   float cota = cotalpha;
   flip_x_ = false;
   flip_y_ = false;
   switch(Dtype_) {
      case 0:
         if(cotbeta < 0.f) {flip_y_ = true;}
         break;
      case 1:
         if(locBz > 0.f) {flip_y_ = true;}
         break;
      case 2:
      case 3:
      case 4:
      case 5:
         if(locBx*locBz < 0.f) {
            cota = fabs(cotalpha);
            flip_x_ = true;
         }
         if(locBx < 0.f) {flip_y_ = true;}
         break;
      default:
#ifndef SI_PIXEL_TEMPLATE_STANDALONE
         throw cms::Exception("DataCorrupt") << "SiPixelTemplate2D::illegal subdetector ID = " << thePixelTemp_[index_id_].head.Dtype << std::endl;
#else
         std::cout << "SiPixelTemplate:2D:illegal subdetector ID = " << thePixelTemp_[index_id_].head.Dtype << std::endl;
#endif
      }
   
   if(cota < cotalpha0_) {
      success_ = false;
      jx0_ = 0;
      jx1_ = 1;
      adcota_ = 0.f;
   } else if(cota > cotalpha1_) {
      success_ = false;
      jx0_ = Nxx_ - 1;
      jx1_ = jx0_ - 1;
      adcota_ = 0.f;
   } else {
      jx0_ = (int)((cota-cotalpha0_)/deltacota_+0.5f);
      dcota = (cota - (cotalpha0_ + jx0_*deltacota_))/deltacota_;
      adcota_ = fabs(dcota);
      if(dcota > 0.f) {jx1_ = jx0_ + 1;if(jx1_ > Nxx_-1) jx1_ = jx0_-1;} else {jx1_ = jx0_ - 1; if(jx1_ < 0) jx1_ = jx0_+1;}
   }
   
   // Interpolate the absolute value of cot(beta)
   
   acotb = fabs(cotbeta);
   
   if(acotb < cotbeta0_) {
      success_ = false;
      iy0_ = 0;
      iy1_ = 1;
      adcotb_ = 0.f;
   } else if(acotb > cotbeta1_) {
      success_ = false;
      iy0_ = Nyx_ - 1;
      iy1_ = iy0_ - 1;
      adcotb_ = 0.f;
   } else {
      iy0_ = (int)((acotb-cotbeta0_)/deltacotb_+0.5f);
      dcotb = (acotb - (cotbeta0_ + iy0_*deltacotb_))/deltacotb_;
      adcotb_ = fabs(dcotb);
      if(dcotb > 0.f) {iy1_ = iy0_ + 1; if(iy1_ > Nyx_-1) iy1_ = iy0_-1;} else {iy1_ = iy0_ - 1; if(iy1_ < 0) iy1_ = iy0_+1;}
   }
   
   //  Calculate signed quantities
   
   lorydrift_ = lorywidth_/2.;
   if(flip_y_) lorydrift_ = -lorydrift_;
   lorxdrift_ = lorxwidth_/2.;
   if(flip_x_) lorxdrift_ = -lorxdrift_;
   
   // Use pointers to the three angle pairs used in the interpolation
   
   
   entry00_ = &thePixelTemp_[index_id_].entry[iy0_][jx0_];
   entry10_ = &thePixelTemp_[index_id_].entry[iy1_][jx0_];
   entry01_ = &thePixelTemp_[index_id_].entry[iy0_][jx1_];
   
   // Interpolate things in cot(alpha)-cot(beta)
   
   qavg_ = entry00_->qavg
   +adcota_*(entry01_->qavg - entry00_->qavg)
   +adcotb_*(entry10_->qavg - entry00_->qavg);
   
   pixmax_ = entry00_->pixmax
   +adcota_*(entry01_->pixmax - entry00_->pixmax)
   +adcotb_*(entry10_->pixmax - entry00_->pixmax);
   
   sxymax_ = entry00_->sxymax
   +adcota_*(entry01_->sxymax - entry00_->sxymax)
   +adcotb_*(entry10_->sxymax - entry00_->sxymax);
   
   chi2avgone_ = entry00_->chi2avgone
   +adcota_*(entry01_->chi2avgone - entry00_->chi2avgone)
   +adcotb_*(entry10_->chi2avgone - entry00_->chi2avgone);
   
   chi2minone_ = entry00_->chi2minone
   +adcota_*(entry01_->chi2minone - entry00_->chi2minone)
   +adcotb_*(entry10_->chi2minone - entry00_->chi2minone);
   
   clsleny_ = entry00_->clsleny
   +adcota_*(entry01_->clsleny - entry00_->clsleny)
   +adcotb_*(entry10_->clsleny - entry00_->clsleny);
   
   clslenx_ = entry00_->clslenx
   +adcota_*(entry01_->clslenx - entry00_->clslenx)
   +adcotb_*(entry10_->clslenx - entry00_->clslenx);
   
   
   chi2ppix_ = entry00_->chi2ppix
   +adcota_*(entry01_->chi2ppix - entry00_->chi2ppix)
   +adcotb_*(entry10_->chi2ppix - entry00_->chi2ppix);
   
   chi2scale_ = entry00_->chi2scale
   +adcota_*(entry01_->chi2scale - entry00_->chi2scale)
   +adcotb_*(entry10_->chi2scale - entry00_->chi2scale);
   
   scaleyavg_ = entry00_->scaleyavg
   +adcota_*(entry01_->scaleyavg - entry00_->scaleyavg)
   +adcotb_*(entry10_->scaleyavg - entry00_->scaleyavg);
   
   scalexavg_ = entry00_->scalexavg
   +adcota_*(entry01_->scalexavg - entry00_->scalexavg)
   +adcotb_*(entry10_->scalexavg - entry00_->scalexavg);
   
   delyavg_ = entry00_->delyavg
   +adcota_*(entry01_->delyavg - entry00_->delyavg)
   +adcotb_*(entry10_->delyavg - entry00_->delyavg);
   
   delysig_ = entry00_->delysig
   +adcota_*(entry01_->delysig - entry00_->delysig)
   +adcotb_*(entry10_->delysig - entry00_->delysig);
   
   mpvvav_ = entry00_->mpvvav
   +adcota_*(entry01_->mpvvav - entry00_->mpvvav)
   +adcotb_*(entry10_->mpvvav - entry00_->mpvvav);
   
   sigmavav_ = entry00_->sigmavav
   +adcota_*(entry01_->sigmavav - entry00_->sigmavav)
   +adcotb_*(entry10_->sigmavav - entry00_->sigmavav);
   
   kappavav_ = entry00_->kappavav
   +adcota_*(entry01_->kappavav - entry00_->kappavav)
   +adcotb_*(entry10_->kappavav - entry00_->kappavav);
   
   for(int i=0; i<4 ; ++i) {
      scalex_[i] = entry00_->scalex[i]
      +adcota_*(entry01_->scalex[i] - entry00_->scalex[i])
      +adcotb_*(entry10_->scalex[i] - entry00_->scalex[i]);
      
      scaley_[i] = entry00_->scaley[i]
      +adcota_*(entry01_->scaley[i] - entry00_->scaley[i])
      +adcotb_*(entry10_->scaley[i] - entry00_->scaley[i]);
      
      offsetx_[i] = entry00_->offsetx[i]
      +adcota_*(entry01_->offsetx[i] - entry00_->offsetx[i])
      +adcotb_*(entry10_->offsetx[i] - entry00_->offsetx[i]);
      if(flip_x_) offsetx_[i] = -offsetx_[i];
      
      offsety_[i] = entry00_->offsety[i]
      +adcota_*(entry01_->offsety[i] - entry00_->offsety[i])
      +adcotb_*(entry10_->offsety[i] - entry00_->offsety[i]);
      if(flip_y_) offsety_[i] = -offsety_[i];
   }
   
   for(int i=0; i<2 ; ++i) {
      for(int j=0; j<5 ; ++j) {
         // Charge loss switches sides when cot(beta) changes sign
         if(flip_y_) {
            xypary0x0_[1-i][j] = (float)entry00_->xypar[i][j];
            xypary1x0_[1-i][j] = (float)entry10_->xypar[i][j];
            xypary0x1_[1-i][j] = (float)entry01_->xypar[i][j];
            lanpar_[1-i][j] = entry00_->lanpar[i][j]
            +adcota_*(entry01_->lanpar[i][j] - entry00_->lanpar[i][j])
            +adcotb_*(entry10_->lanpar[i][j] - entry00_->lanpar[i][j]);
         } else {
            xypary0x0_[i][j] = (float)entry00_->xypar[i][j];
            xypary1x0_[i][j] = (float)entry10_->xypar[i][j];
            xypary0x1_[i][j] = (float)entry01_->xypar[i][j];
            lanpar_[i][j] = entry00_->lanpar[i][j]
            +adcota_*(entry01_->lanpar[i][j] - entry00_->lanpar[i][j])
            +adcotb_*(entry10_->lanpar[i][j] - entry00_->lanpar[i][j]);
         }
      }
   }
   
   return success_;
} // interpolate

float SiPixelTemplate2D::kappavav ( )

inline

kappa parameter in Vavilov distribution

Definition at line 199 of file SiPixelTemplate2D.h.

Referenced by SiPixelTemplateReco2D::PixelTempReco2D().

void SiPixelTemplate2D::landau_par

(

float

lanpar[2][5]

)

Return the Landau probability parameters for this set of cot(alpha, cot(beta)

Definition at line 1563 of file SiPixelTemplate2D.cc.

References mps_fire::i.

Referenced by SiPixelTemplateSplit::PixelTempSplit().

{
   // Interpolate using quantities already stored in the private variables

   for(int i=0; i<2; ++i) {
      for(int j=0; j<5; ++j) {
         lanpar[i][j] = lanpar_[i][j];
      }
   }
   return;
   
} // End lan_par

float SiPixelTemplate2D::lorxdrift ( )

inline

signed lorentz x-width (microns)

Definition at line 201 of file SiPixelTemplate2D.h.

Referenced by SiPixelTemplateReco2D::PixelTempReco2D().

float SiPixelTemplate2D::lorydrift ( )

inline

signed lorentz y-width (microns)

Definition at line 200 of file SiPixelTemplate2D.h.

Referenced by SiPixelTemplateReco2D::PixelTempReco2D().

float SiPixelTemplate2D::mpvvav ( )

inline

most probable Q in Vavilov distribution

Definition at line 197 of file SiPixelTemplate2D.h.

Referenced by SiPixelTemplateReco2D::PixelTempReco2D().

float SiPixelTemplate2D::offsetx ( int  i )

inline

x-offset in 4 charge bins

Definition at line 188 of file SiPixelTemplate2D.h.

References mps_fire::i.

Referenced by SiPixelTemplateReco2D::PixelTempReco2D().

float SiPixelTemplate2D::offsety ( int  i )

inline

y-offset in 4 charge bins

Definition at line 189 of file SiPixelTemplate2D.h.

References mps_fire::i.

Referenced by SiPixelTemplateReco2D::PixelTempReco2D().

float SiPixelTemplate2D::pixmax ( )

inline

maximum pixel charge

Definition at line 182 of file SiPixelTemplate2D.h.

Referenced by SiPixelTemplateReco2D::PixelTempReco2D().

bool SiPixelTemplate2D::pushfile ( int  filenum, std::vector< SiPixelTemplateStore2D > &  pixelTemp, std::string  dir = ""  )

static

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

Parameters

filenum

- an integer NNNN used in the filename template_summary_zpNNNN

Alt implementation: for (unsigned cnt=4-tempfile.length(); cnt > 0; cnt– ){ tempfile = "0" + tempfile; }

Definition at line 62 of file SiPixelTemplate2D.cc.

References SiPixelTemplateHeader2D::Bfield, EnergyCorrector::c, SiPixelTemplateEntry2D::chi2ppix, SiPixelTemplateEntry2D::chi2scale, SiPixelTemplateEntry2D::clslenx, SiPixelTemplateEntry2D::clsleny, SiPixelTemplateEntry2D::costrk, SiPixelTemplateEntry2D::cotalpha, SiPixelTemplateEntry2D::cotbeta, SiPixelTemplateEntry2D::delyavg, SiPixelTemplateEntry2D::delysig, SiPixelTemplateHeader2D::Dtype, ENDL, SiPixelTemplateStore2D::entry, SiPixelTemplateHeader2D::fbin, FrontierConditions_GlobalTag_cff::file, SiPixelTemplateHeader2D::fluence, edm::FileInPath::fullPath(), SiPixelTemplateStore2D::head, mps_fire::i, SiPixelTemplateHeader2D::ID, runGCPTkAlMap::in_file, createfilelist::int, SiPixelTemplateEntry2D::iymax, SiPixelTemplateEntry2D::iymin, SiPixelTemplateEntry2D::jxmax, SiPixelTemplateEntry2D::jxmin, gen::k, SiPixelTemplateEntry2D::kappavav, checklumidiff::l, SiPixelTemplateEntry2D::lanpar, LOGERROR, LOGINFO, SiPixelTemplateHeader2D::lorxbias, SiPixelTemplateHeader2D::lorxwidth, SiPixelTemplateHeader2D::lorybias, SiPixelTemplateHeader2D::lorywidth, SiPixelTemplateEntry2D::mpvvav, SiPixelTemplateHeader2D::NTxx, SiPixelTemplateHeader2D::NTy, SiPixelTemplateHeader2D::NTyx, SiPixelTemplateEntry2D::offsetx, SiPixelTemplateEntry2D::offsety, SiPixelTemplateEntry2D::pixmax, SiPixelTemplateEntry2D::qavg, SiPixelTemplateHeader2D::qscale, SiPixelTemplateEntry2D::runnum, SiPixelTemplateHeader2D::s50, SiPixelTemplateEntry2D::scalex, SiPixelTemplateEntry2D::scalexavg, SiPixelTemplateEntry2D::scaley, SiPixelTemplateEntry2D::scaleyavg, SiPixelTemplateEntry2D::sigmavav, SiPixelTemplateEntry2D::spare, SiPixelTemplateHeader2D::ss50, AlCaHLTBitMon_QueryRunRegistry::string, SiPixelTemplateEntry2D::sxymax, T2XSIZE, T2YSIZE, SiPixelTemplateHeader2D::temperature, SiPixelTemplateHeader2D::templ_version, SiPixelTemplateHeader2D::title, SiPixelTemplateHeader2D::Vbias, SiPixelTemplateHeader2D::xsize, SiPixelTemplateEntry2D::xypar, SiPixelTemplateEntry2D::xytemp, SiPixelTemplateHeader2D::ysize, and SiPixelTemplateHeader2D::zsize.

Referenced by SiPixelDigitizerAlgorithm::init(), and PixelCPEClusterRepair::PixelCPEClusterRepair().

{
   // Add template stored in external file numbered filenum to theTemplateStore
   
   // Local variables
   const int code_version={21};
   
   //  Create a filename for this run
   std::string tempfile = std::to_string(filenum);
   
   
   //  Create different path in CMSSW than standalone
   
#ifndef SI_PIXEL_TEMPLATE_STANDALONE
   // If integer filenum has less than 4 digits, prepend 0's until we have four numerical characters, e.g. "0292"
   int nzeros = 4-tempfile.length();
   if (nzeros > 0)
     tempfile = std::string(nzeros, '0') + tempfile;

   tempfile = dir + "template_summary2D_zp" + tempfile + ".out";
   edm::FileInPath file( tempfile );         // Find the file in CMSSW
   tempfile = file.fullPath();               // Put it back with the whole path.

#else
   // This is the same as above, but more elegant.
   std::ostringstream tout;
   tout << "template_summary2D_zp" << std::setw(4) << std::setfill('0') << std::right << filenum << ".out" << std::ends;
   tempfile = tout.str();

#endif
   
   //  Open the template file
   //
   std::ifstream in_file(tempfile);
   if(in_file.is_open() && in_file.good()) {

      // Create a local template storage entry
      SiPixelTemplateStore2D theCurrentTemp;
      
      // Read-in a header string first and print it
      char c;
      for (int i=0; (c=in_file.get()) != '\n'; ++i) {
         if(i < 79) {theCurrentTemp.head.title[i] = c;}
         else       {theCurrentTemp.head.title[79] ='\0';}
      }
      LOGINFO("SiPixelTemplate2D") << "Loading Pixel Template 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;
      
      if(in_file.fail()) {LOGERROR("SiPixelTemplate2D") << "Error reading file 0A, no template load" << ENDL; return false;}
      
      if(theCurrentTemp.head.templ_version > 17) {
         in_file >> 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("SiPixelTemplate2D") << "Error reading file 0B, no template load"
            << ENDL; return false;}
      } else {
// This is for older [legacy] payloads 
         theCurrentTemp.head.ss50 = theCurrentTemp.head.s50;
         theCurrentTemp.head.lorybias = theCurrentTemp.head.lorywidth/2.f;
         theCurrentTemp.head.lorxbias = theCurrentTemp.head.lorxwidth/2.f;
         theCurrentTemp.head.fbin[0] = 1.5f;
         theCurrentTemp.head.fbin[1] = 1.00f;
         theCurrentTemp.head.fbin[2] = 0.85f;
      }
      
      LOGINFO("SiPixelTemplate2D") << "Template ID = " << theCurrentTemp.head.ID << ", Template 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("SiPixelTemplate2D") << "code expects version " << code_version << ", no template load" << ENDL; return false;}
      
      if(theCurrentTemp.head.NTy != 0) {LOGERROR("SiPixelTemplate2D") << "Trying to load 1-d template info into the 2-d template object, check your DB/global tag!" << ENDL; return false;}
      
      // next, layout the 2-d structure needed to store template
      
      theCurrentTemp.entry = new SiPixelTemplateEntry2D*[theCurrentTemp.head.NTyx];
      theCurrentTemp.entry[0] = new SiPixelTemplateEntry2D[theCurrentTemp.head.NTyx*theCurrentTemp.head.NTxx];
      for(int i = 1; i < theCurrentTemp.head.NTyx; ++i) theCurrentTemp.entry[i] = theCurrentTemp.entry[i-1] + theCurrentTemp.head.NTxx;      

      // Read in the file info
      
      for (int iy=0; iy < theCurrentTemp.head.NTyx; ++iy) {
         for(int jx=0; jx < theCurrentTemp.head.NTxx; ++jx) {
            
            in_file >> theCurrentTemp.entry[iy][jx].runnum >> theCurrentTemp.entry[iy][jx].costrk[0]
            >> theCurrentTemp.entry[iy][jx].costrk[1] >> theCurrentTemp.entry[iy][jx].costrk[2];
            
            if(in_file.fail()) {
          LOGERROR("SiPixelTemplate2D") << "Error reading file 1, no template load, run # " << theCurrentTemp.entry[iy][jx].runnum << ENDL;
          delete [] theCurrentTemp.entry[0];
          delete [] theCurrentTemp.entry;
          return false;
        }
            
            // Calculate cot(alpha) and cot(beta) for this entry
            
            theCurrentTemp.entry[iy][jx].cotalpha = theCurrentTemp.entry[iy][jx].costrk[0]/theCurrentTemp.entry[iy][jx].costrk[2];
            
            theCurrentTemp.entry[iy][jx].cotbeta = theCurrentTemp.entry[iy][jx].costrk[1]/theCurrentTemp.entry[iy][jx].costrk[2];
            
            in_file >> theCurrentTemp.entry[iy][jx].qavg >> theCurrentTemp.entry[iy][jx].pixmax >> theCurrentTemp.entry[iy][jx].sxymax >> theCurrentTemp.entry[iy][jx].iymin
            >> theCurrentTemp.entry[iy][jx].iymax >> theCurrentTemp.entry[iy][jx].jxmin >> theCurrentTemp.entry[iy][jx].jxmax;
            
            if(in_file.fail()) {
          LOGERROR("SiPixelTemplate2D") << "Error reading file 2, no template load, run # " << theCurrentTemp.entry[iy][jx].runnum << ENDL;
          delete [] theCurrentTemp.entry[0];
          delete [] theCurrentTemp.entry; 
          return false;
        }
            
            for (int k=0; k<2; ++k) {
               
               in_file >> theCurrentTemp.entry[iy][jx].xypar[k][0] >> theCurrentTemp.entry[iy][jx].xypar[k][1]
               >> theCurrentTemp.entry[iy][jx].xypar[k][2] >> theCurrentTemp.entry[iy][jx].xypar[k][3] >> theCurrentTemp.entry[iy][jx].xypar[k][4];
               
               if(in_file.fail()) {
         LOGERROR("SiPixelTemplate2D") << "Error reading file 3, no template load, run # " << theCurrentTemp.entry[iy][jx].runnum << ENDL; 
         delete [] theCurrentTemp.entry[0];
         delete [] theCurrentTemp.entry;
         return false;
           }               
            }
            
            for (int k=0; k<2; ++k) {
               
               in_file >> theCurrentTemp.entry[iy][jx].lanpar[k][0] >> theCurrentTemp.entry[iy][jx].lanpar[k][1]
               >> theCurrentTemp.entry[iy][jx].lanpar[k][2] >> theCurrentTemp.entry[iy][jx].lanpar[k][3] >> theCurrentTemp.entry[iy][jx].lanpar[k][4];
               
               if(in_file.fail()) {
         LOGERROR("SiPixelTemplate2D") << "Error reading file 4, no template load, run # " << theCurrentTemp.entry[iy][jx].runnum << ENDL; 
         delete [] theCurrentTemp.entry[0];
         delete [] theCurrentTemp.entry;
         return false;
           }
               
            }
            
      
//  Read the 2D template entries as floats [they are formatted that way] and cast to short ints

            float dummy[T2YSIZE];
            for (int l=0; l<7; ++l) {
               for (int k=0; k<7; ++k) {
                  for (int j=0; j<T2XSIZE; ++j) {
                     for (int i=0; i<T2YSIZE; ++i) {
                        in_file >> dummy[i];                       
                     }                     
                     if(in_file.fail()) {
               LOGERROR("SiPixelTemplate2D") << "Error reading file 5, no template load, run # " << theCurrentTemp.entry[iy][jx].runnum << ENDL; 
               delete [] theCurrentTemp.entry[0];
               delete [] theCurrentTemp.entry;
               return false;
             }
                     for (int i=0; i<T2YSIZE; ++i) {
                        theCurrentTemp.entry[iy][jx].xytemp[k][l][i][j] = (short int)dummy[i];                        
                     }                    
                  }
               }
            }
            
            
            in_file >> theCurrentTemp.entry[iy][jx].chi2ppix >> theCurrentTemp.entry[iy][jx].chi2scale >> theCurrentTemp.entry[iy][jx].offsetx[0] >> theCurrentTemp.entry[iy][jx].offsetx[1]
            >> theCurrentTemp.entry[iy][jx].offsetx[2] >> theCurrentTemp.entry[iy][jx].offsetx[3]>> theCurrentTemp.entry[iy][jx].offsety[0] >> theCurrentTemp.entry[iy][jx].offsety[1]
            >> theCurrentTemp.entry[iy][jx].offsety[2] >> theCurrentTemp.entry[iy][jx].offsety[3];
            
            if(in_file.fail()) {
          LOGERROR("SiPixelTemplate2D") << "Error reading file 6, no template load, run # " << theCurrentTemp.entry[iy][jx].runnum << ENDL; 
          delete [] theCurrentTemp.entry[0];
          delete [] theCurrentTemp.entry;
          return false;
        }
            
            in_file >> theCurrentTemp.entry[iy][jx].clsleny >> theCurrentTemp.entry[iy][jx].clslenx >> theCurrentTemp.entry[iy][jx].mpvvav >> theCurrentTemp.entry[iy][jx].sigmavav
            >> theCurrentTemp.entry[iy][jx].kappavav >> theCurrentTemp.entry[iy][jx].scalexavg >> theCurrentTemp.entry[iy][jx].scaleyavg >> theCurrentTemp.entry[iy][jx].delyavg >> theCurrentTemp.entry[iy][jx].delysig >> theCurrentTemp.entry[iy][jx].spare[0];
            
            if(in_file.fail()) {
          LOGERROR("SiPixelTemplate2D") << "Error reading file 7, no template load, run # " << theCurrentTemp.entry[iy][jx].runnum << ENDL; 
          delete [] theCurrentTemp.entry[0];
          delete [] theCurrentTemp.entry;         
          return false;
        }
            
            in_file >> theCurrentTemp.entry[iy][jx].scalex[0] >> theCurrentTemp.entry[iy][jx].scalex[1] >> theCurrentTemp.entry[iy][jx].scalex[2] >> theCurrentTemp.entry[iy][jx].scalex[3]
            >> theCurrentTemp.entry[iy][jx].scaley[0] >> theCurrentTemp.entry[iy][jx].scaley[1] >> theCurrentTemp.entry[iy][jx].scaley[2] >> theCurrentTemp.entry[iy][jx].scaley[3]
            >> theCurrentTemp.entry[iy][jx].spare[1]  >> theCurrentTemp.entry[iy][jx].spare[2];
            
            if(in_file.fail()) {
          LOGERROR("SiPixelTemplate2D") << "Error reading file 8, no template load, run # " << theCurrentTemp.entry[iy][jx].runnum << ENDL; 
          delete [] theCurrentTemp.entry[0];
          delete [] theCurrentTemp.entry;         
          return false;
        }
            
         }
         
      }
      
      
      in_file.close();
      
      // Add this template to the store
      
      pixelTemp.push_back(theCurrentTemp);
      
      return true;
      
   } else {
      
      // If file didn't open, report this
      
      LOGERROR("SiPixelTemplate2D") << "Error opening File" << tempfile << ENDL;
      return false;
      
   }
   
} // TempInit

bool SiPixelTemplate2D::pushfile ( const SiPixel2DTemplateDBObject &  dbobject, std::vector< SiPixelTemplateStore2D > &  pixelTemp  )

static

This routine initializes the global template structures from an external file template_summary_zpNNNN where NNNN are four digits

Parameters

dbobject

- db storing multiple template calibrations

Definition at line 301 of file SiPixelTemplate2D.cc.

References SiPixelTemplateHeader2D::Bfield, SiPixel2DTemplateDBObject::char2float::c, SiPixelTemplateEntry2D::chi2ppix, SiPixelTemplateEntry2D::chi2scale, SiPixelTemplateEntry2D::clslenx, SiPixelTemplateEntry2D::clsleny, SiPixelTemplateEntry2D::costrk, SiPixelTemplateEntry2D::cotalpha, SiPixelTemplateEntry2D::cotbeta, SiPixelTemplateEntry2D::delyavg, SiPixelTemplateEntry2D::delysig, SiPixelTemplateHeader2D::Dtype, ENDL, SiPixelTemplateStore2D::entry, SiPixel2DTemplateDBObject::char2float::f, SiPixel2DTemplateDBObject::fail(), SiPixelTemplateHeader2D::fbin, SiPixelTemplateHeader2D::fluence, SiPixelTemplateStore2D::head, mps_fire::i, SiPixelTemplateHeader2D::ID, SiPixel2DTemplateDBObject::incrementIndex(), SiPixel2DTemplateDBObject::index(), createfilelist::int, SiPixelTemplateEntry2D::iymax, SiPixelTemplateEntry2D::iymin, SiPixelTemplateEntry2D::jxmax, SiPixelTemplateEntry2D::jxmin, gen::k, SiPixelTemplateEntry2D::kappavav, checklumidiff::l, SiPixelTemplateEntry2D::lanpar, LOGERROR, LOGINFO, SiPixelTemplateHeader2D::lorxbias, SiPixelTemplateHeader2D::lorxwidth, SiPixelTemplateHeader2D::lorybias, SiPixelTemplateHeader2D::lorywidth, funct::m, SiPixelTemplateEntry2D::mpvvav, SiPixelTemplateHeader2D::NTxx, SiPixelTemplateHeader2D::NTy, SiPixelTemplateHeader2D::NTyx, SiPixel2DTemplateDBObject::numOfTempl(), SiPixelTemplateEntry2D::offsetx, SiPixelTemplateEntry2D::offsety, SiPixelTemplateEntry2D::pixmax, SiPixelTemplateEntry2D::qavg, SiPixelTemplateHeader2D::qscale, SiPixelTemplateEntry2D::runnum, SiPixelTemplateHeader2D::s50, SiPixelTemplateEntry2D::scalex, SiPixelTemplateEntry2D::scalexavg, SiPixelTemplateEntry2D::scaley, SiPixelTemplateEntry2D::scaleyavg, SiPixelTemplateEntry2D::sigmavav, SiPixelTemplateEntry2D::spare, SiPixelTemplateHeader2D::ss50, SiPixel2DTemplateDBObject::sVector(), SiPixelTemplateEntry2D::sxymax, T2XSIZE, T2YSIZE, groupFilesInBlocks::temp, SiPixelTemplateHeader2D::temperature, SiPixelTemplateHeader2D::templ_version, SiPixelTemplateHeader2D::title, SiPixelTemplateHeader2D::Vbias, SiPixelTemplateHeader2D::xsize, SiPixelTemplateEntry2D::xypar, SiPixelTemplateEntry2D::xytemp, SiPixelTemplateHeader2D::ysize, and SiPixelTemplateHeader2D::zsize.

{
   // Add template stored in external dbobject to theTemplateStore
   
   const int code_version={21};
   
   // We must create a new object because dbobject must be a const and our stream must not be
   SiPixel2DTemplateDBObject db = dbobject;
   
   // Create a local template storage entry
   SiPixelTemplateStore2D theCurrentTemp;
   
   // Fill the template storage for each template calibration stored in the db
   for(int m=0; m<db.numOfTempl(); ++m)
   {
      
      // Read-in a header string first and print it
      
      SiPixel2DTemplateDBObject::char2float temp;
      for (int 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("SiPixelTemplate2D") << "Loading Pixel Template 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;
      
      if(db.fail()) {LOGERROR("SiPixelTemplate2D") << "Error reading file 0A, no template load" << ENDL; return false;}
      
      LOGINFO("SiPixelTemplate2D") << "Loading Pixel Template File - " << theCurrentTemp.head.title
      <<" code version = "<<code_version
      <<" object version "<<theCurrentTemp.head.templ_version
      << ENDL;
      
      if(theCurrentTemp.head.templ_version > 17) {
         db >> 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("SiPixelTemplate2D") << "Error reading file 0B, no template load"
            << ENDL; return false;}
      } else {
// This is for older [legacy] payloads and the numbers are indeed magic [they are part of the payload for v>17]
         theCurrentTemp.head.ss50 = theCurrentTemp.head.s50;
         theCurrentTemp.head.lorybias = theCurrentTemp.head.lorywidth/2.f;
         theCurrentTemp.head.lorxbias = theCurrentTemp.head.lorxwidth/2.f;
         theCurrentTemp.head.fbin[0] = 1.50f;
         theCurrentTemp.head.fbin[1] = 1.00f;
         theCurrentTemp.head.fbin[2] = 0.85f;
      }
      
      LOGINFO("SiPixelTemplate2D")
      << "Template ID = " << theCurrentTemp.head.ID << ", Template 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) {
         LOGINFO("SiPixelTemplate2D") << "code expects version "<< code_version << " finds "
         << theCurrentTemp.head.templ_version <<", load anyway " << ENDL;
      }
      
      if(theCurrentTemp.head.NTy != 0) {LOGERROR("SiPixelTemplate2D") << "Trying to load 1-d template info into the 2-d template object, check your DB/global tag!" << ENDL; return false;}
      
      
      // next, layout the 2-d structure needed to store template
      
      theCurrentTemp.entry = new SiPixelTemplateEntry2D*[theCurrentTemp.head.NTyx];
      theCurrentTemp.entry[0] = new SiPixelTemplateEntry2D[theCurrentTemp.head.NTyx*theCurrentTemp.head.NTxx];
      for(int i = 1; i < theCurrentTemp.head.NTyx; ++i) theCurrentTemp.entry[i] = theCurrentTemp.entry[i-1] + theCurrentTemp.head.NTxx;
      
      // Read in the file info
      
      for (int iy=0; iy < theCurrentTemp.head.NTyx; ++iy) {
         for(int jx=0; jx < theCurrentTemp.head.NTxx; ++jx) {
            
            db >> theCurrentTemp.entry[iy][jx].runnum >> theCurrentTemp.entry[iy][jx].costrk[0]
            >> theCurrentTemp.entry[iy][jx].costrk[1] >> theCurrentTemp.entry[iy][jx].costrk[2];
            
            if(db.fail()) {
          LOGERROR("SiPixelTemplate2D") << "Error reading file 1, no template load, run # " << theCurrentTemp.entry[iy][jx].runnum << ENDL; 
          delete [] theCurrentTemp.entry[0];
          delete [] theCurrentTemp.entry;         
          return false;
        }
            
            // Calculate cot(alpha) and cot(beta) for this entry
            
            theCurrentTemp.entry[iy][jx].cotalpha = theCurrentTemp.entry[iy][jx].costrk[0]/theCurrentTemp.entry[iy][jx].costrk[2];
            
            theCurrentTemp.entry[iy][jx].cotbeta = theCurrentTemp.entry[iy][jx].costrk[1]/theCurrentTemp.entry[iy][jx].costrk[2];
            
            db >> theCurrentTemp.entry[iy][jx].qavg >> theCurrentTemp.entry[iy][jx].pixmax >> theCurrentTemp.entry[iy][jx].sxymax >> theCurrentTemp.entry[iy][jx].iymin
            >> theCurrentTemp.entry[iy][jx].iymax >> theCurrentTemp.entry[iy][jx].jxmin >> theCurrentTemp.entry[iy][jx].jxmax;
            
            if(db.fail()) {
          LOGERROR("SiPixelTemplate2D") << "Error reading file 2, no template load, run # " << theCurrentTemp.entry[iy][jx].runnum << ENDL; 
          delete [] theCurrentTemp.entry[0];
          delete [] theCurrentTemp.entry;         
          return false;
        }
            
            for (int k=0; k<2; ++k) {
               
               db >> theCurrentTemp.entry[iy][jx].xypar[k][0] >> theCurrentTemp.entry[iy][jx].xypar[k][1]
               >> theCurrentTemp.entry[iy][jx].xypar[k][2] >> theCurrentTemp.entry[iy][jx].xypar[k][3] >> theCurrentTemp.entry[iy][jx].xypar[k][4];
               
               if(db.fail()) {
         LOGERROR("SiPixelTemplate2D") << "Error reading file 3, no template load, run # " << theCurrentTemp.entry[iy][jx].runnum << ENDL; 
         delete [] theCurrentTemp.entry[0];
         delete [] theCurrentTemp.entry;          
         return false;
           }
               
            }
            
            for (int k=0; k<2; ++k) {
               
               db >> theCurrentTemp.entry[iy][jx].lanpar[k][0] >> theCurrentTemp.entry[iy][jx].lanpar[k][1]
               >> theCurrentTemp.entry[iy][jx].lanpar[k][2] >> theCurrentTemp.entry[iy][jx].lanpar[k][3] >> theCurrentTemp.entry[iy][jx].lanpar[k][4];
               
               if(db.fail()) {
         LOGERROR("SiPixelTemplate2D") << "Error reading file 4, no template load, run # " << theCurrentTemp.entry[iy][jx].runnum << ENDL; 
         delete [] theCurrentTemp.entry[0];
         delete [] theCurrentTemp.entry;          
         return false;
           }
               
            }
            
//  Read the 2D template entries as floats [they are formatted that way] and cast to short ints

            float dummy[T2YSIZE];
            for (int l=0; l<7; ++l) {
               for (int k=0; k<7; ++k) {
                  for (int j=0; j<T2XSIZE; ++j) {
                     for (int i=0; i<T2YSIZE; ++i) {
                        db >> dummy[i];                       
                     }                     
                     if(db.fail()) {
               LOGERROR("SiPixelTemplate2D") << "Error reading file 5, no template load, run # " << theCurrentTemp.entry[iy][jx].runnum << ENDL; 
               delete [] theCurrentTemp.entry[0];
               delete [] theCurrentTemp.entry;        
               return false;
             }
                     for (int i=0; i<T2YSIZE; ++i) {
                        theCurrentTemp.entry[iy][jx].xytemp[k][l][i][j] = (short int)dummy[i];                        
                     }                    
                  }
               }
            }
            
            db >> theCurrentTemp.entry[iy][jx].chi2ppix >> theCurrentTemp.entry[iy][jx].chi2scale >> theCurrentTemp.entry[iy][jx].offsetx[0] >> theCurrentTemp.entry[iy][jx].offsetx[1]
            >> theCurrentTemp.entry[iy][jx].offsetx[2] >> theCurrentTemp.entry[iy][jx].offsetx[3]>> theCurrentTemp.entry[iy][jx].offsety[0] >> theCurrentTemp.entry[iy][jx].offsety[1]
            >> theCurrentTemp.entry[iy][jx].offsety[2] >> theCurrentTemp.entry[iy][jx].offsety[3];
            
            if(db.fail()) {
          LOGERROR("SiPixelTemplate2D") << "Error reading file 6, no template load, run # " << theCurrentTemp.entry[iy][jx].runnum << ENDL; 
          delete [] theCurrentTemp.entry[0];
          delete [] theCurrentTemp.entry;         
          return false;
        }
            
            db >> theCurrentTemp.entry[iy][jx].clsleny >> theCurrentTemp.entry[iy][jx].clslenx >> theCurrentTemp.entry[iy][jx].mpvvav >> theCurrentTemp.entry[iy][jx].sigmavav
            >> theCurrentTemp.entry[iy][jx].kappavav >> theCurrentTemp.entry[iy][jx].scalexavg >> theCurrentTemp.entry[iy][jx].scaleyavg >> theCurrentTemp.entry[iy][jx].delyavg >> theCurrentTemp.entry[iy][jx].delysig >> theCurrentTemp.entry[iy][jx].spare[0];
            
            if(db.fail()) {
          LOGERROR("SiPixelTemplate2D") << "Error reading file 7, no template load, run # " << theCurrentTemp.entry[iy][jx].runnum << ENDL; 
          delete [] theCurrentTemp.entry[0];
          delete [] theCurrentTemp.entry;         
          return false;
        }
            
            db >> theCurrentTemp.entry[iy][jx].scalex[0] >> theCurrentTemp.entry[iy][jx].scalex[1] >> theCurrentTemp.entry[iy][jx].scalex[2] >> theCurrentTemp.entry[iy][jx].scalex[3]
            >> theCurrentTemp.entry[iy][jx].scaley[0] >> theCurrentTemp.entry[iy][jx].scaley[1] >> theCurrentTemp.entry[iy][jx].scaley[2] >> theCurrentTemp.entry[iy][jx].scaley[3]
            >> theCurrentTemp.entry[iy][jx].spare[1]  >> theCurrentTemp.entry[iy][jx].spare[2];
            
            if(db.fail()) {
          LOGERROR("SiPixelTemplate2D") << "Error reading file 8, no template load, run # " << theCurrentTemp.entry[iy][jx].runnum << ENDL; 
          delete [] theCurrentTemp.entry[0];
          delete [] theCurrentTemp.entry;         
          return false;
        }
            
         }
      }
         
// Add this template to the store
   
   pixelTemp.push_back(theCurrentTemp);    
     
   }
   
   
   return true;
   
} // TempInit

float SiPixelTemplate2D::qavg ( )

inline

average cluster charge for this set of track angles

Definition at line 181 of file SiPixelTemplate2D.h.

Referenced by SiPixelTemplateReco2D::PixelTempReco2D().

float SiPixelTemplate2D::qscale ( )

inline

charge scaling factor

Definition at line 183 of file SiPixelTemplate2D.h.

Referenced by SiPixelTemplateReco2D::PixelTempReco2D().

float SiPixelTemplate2D::s50 ( )

inline

1/2 of the pixel threshold signal in adc units

Definition at line 184 of file SiPixelTemplate2D.h.

Referenced by SiPixelTemplateReco2D::PixelTempReco2D(), and SiPixelDigitizerAlgorithm::PixelTempRewgt2D().

float SiPixelTemplate2D::scalex ( int  i )

inline

x-error scale factor in 4 charge bins

Definition at line 186 of file SiPixelTemplate2D.h.

References mps_fire::i.

Referenced by SiPixelTemplateReco2D::PixelTempReco2D().

float SiPixelTemplate2D::scalexavg ( )

inline

x-reco error scaling factor

Definition at line 205 of file SiPixelTemplate2D.h.

float SiPixelTemplate2D::scaley ( int  i )

inline

y-error scale factor in 4 charge bins

Definition at line 187 of file SiPixelTemplate2D.h.

References mps_fire::i.

Referenced by SiPixelTemplateReco2D::PixelTempReco2D().

float SiPixelTemplate2D::scaleyavg ( )

inline

y-reco error scaling factor

Definition at line 204 of file SiPixelTemplate2D.h.

void SiPixelTemplate2D::sideload	(	SiPixelTemplateEntry2D *	entry,
		int	iDtype,
		float	locBx,
		float	locBz,
		float	lorwdy,
		float	lorwdx,
		float	q50,
		float	fbin[3],
		float	xsize,
		float	ysize,
		float	zsize
	)

Load template info for single angle point to invoke template reco for template generation

Parameters

entry	- (input) pointer to template entry
sizex	- (input) pixel x-size
sizey	- (input) pixel y-size
sizez	- (input) pixel z-size

Definition at line 841 of file SiPixelTemplate2D.cc.

References gather_cfg::cout, mps_splice::entry, Exception, f, objects.autophobj::float, mps_fire::i, xsize, and ysize.

{
   // Set class variables to the input parameters
   
   entry00_ = entry;
   entry01_ = entry;
   entry10_ = entry;
   Dtype_ = iDtype;
   lorywidth_ = lorwdy;
   lorxwidth_ = lorwdx;
   xsize_ = xsize;
   ysize_ = ysize;
   zsize_ = zsize;
   s50_ = q50;
   qscale_ = 1.f;
   for(int i=0; i<3; ++i) {fbin_[i] = fbin[i];}
   
   // Set other class variables
   
   adcota_ = 0.f;
   adcotb_ = 0.f;
   
   // Interpolate things in cot(alpha)-cot(beta)
   
   qavg_ = entry00_->qavg;
   
   pixmax_ = entry00_->pixmax;
   
   sxymax_ = entry00_->sxymax;
   
   clsleny_ = entry00_->clsleny;
   
   clslenx_ = entry00_->clslenx;
   
   scaleyavg_ = 1.f;
   
   scalexavg_ = 1.f;
   
   delyavg_ = 0.f;
   
   delysig_ = 0.f;
   
   for(int i=0; i<4 ; ++i) {
      scalex_[i] = 1.f;
      scaley_[i] = 1.f;
      offsetx_[i] = 0.f;
      offsety_[i] = 0.f;
   }
   
   // This works only for IP-related tracks
   
   flip_x_ = false;
   flip_y_ = false;
   float cotbeta = entry00_->cotbeta;
   switch(Dtype_) {
      case 0:
         if(cotbeta < 0.f) {flip_y_ = true;}
         break;
      case 1:
         if(locBz > 0.f) {
            flip_y_ = true;
         }
         break;
      case 2:
      case 3:
      case 4:
      case 5:
         if(locBx*locBz < 0.f) {
            flip_x_ = true;
         }
         if(locBx < 0.f) {
            flip_y_ = true;
         }
         break;
      default:
#ifndef SI_PIXEL_TEMPLATE_STANDALONE
         throw cms::Exception("DataCorrupt") << "SiPixelTemplate2D::illegal subdetector ID = " << iDtype << std::endl;
#else
         std::cout << "SiPixelTemplate:2D:illegal subdetector ID = " << iDtype << std::endl;
#endif
   }
   
   //  Calculate signed quantities
   
   lorydrift_ = lorywidth_/2.;
   if(flip_y_) lorydrift_ = -lorydrift_;
   lorxdrift_ = lorxwidth_/2.;
   if(flip_x_) lorxdrift_ = -lorxdrift_;
   
   for(int i=0; i<2 ; ++i) {
      for(int j=0; j<5 ; ++j) {
         // Charge loss switches sides when cot(beta) changes sign
         if(flip_y_) {
            xypary0x0_[1-i][j] = (float)entry00_->xypar[i][j];
            xypary1x0_[1-i][j] = (float)entry00_->xypar[i][j];
            xypary0x1_[1-i][j] = (float)entry00_->xypar[i][j];
            lanpar_[1-i][j] = entry00_->lanpar[i][j];
         } else {
            xypary0x0_[i][j] = (float)entry00_->xypar[i][j];
            xypary1x0_[i][j] = (float)entry00_->xypar[i][j];
            xypary0x1_[i][j] = (float)entry00_->xypar[i][j];
            lanpar_[i][j] = entry00_->lanpar[i][j];
         }
      }
   }
   return;
}

float SiPixelTemplate2D::sigmavav ( )

inline

scale factor in Vavilov distribution

Definition at line 198 of file SiPixelTemplate2D.h.

Referenced by SiPixelTemplateReco2D::PixelTempReco2D().

float SiPixelTemplate2D::sizex ( )

inline

return x size of template cluster

Definition at line 191 of file SiPixelTemplate2D.h.

float SiPixelTemplate2D::sizey ( )

inline

return y size of template cluster

Definition at line 192 of file SiPixelTemplate2D.h.

int SiPixelTemplate2D::storesize ( )

inline

return the size of the template store (the number of stored IDs

Definition at line 211 of file SiPixelTemplate2D.h.

float SiPixelTemplate2D::sxymax ( )

inline

max pixel signal for pixel error calculation

Definition at line 185 of file SiPixelTemplate2D.h.

Referenced by SiPixelTemplateReco2D::PixelTempReco2D().

float SiPixelTemplate2D::xsize ( )

inline

pixel x-size (microns)

Definition at line 208 of file SiPixelTemplate2D.h.

Referenced by SiPixelTemplateReco2D::PixelTempReco2D(), and SiPixelDigitizerAlgorithm::PixelTempRewgt2D().

void SiPixelTemplate2D::xysigma2	(	float	qpixel,
		int	index,
		float &	xysig2
	)

Return y error (squared) for an input signal and yindex Add large Q scaling for use in cluster splitting.

Parameters

qpixel	- (input) pixel charge
index	- (input) y-index index of pixel
xysig2	- (output) square error

Definition at line 1511 of file SiPixelTemplate2D.cc.

References BYM2, Exception, f, and T2HYP1.

Referenced by SiPixelTemplateReco2D::PixelTempReco2D().

{
   // Interpolate using quantities already stored in the private variables
   
   // Local variables
   float sigi, sigi2, sigi3, sigi4, qscale, err2, err00;
   
   // Make sure that input is OK
   
#ifndef SI_PIXEL_TEMPLATE_STANDALONE
   if(index < 1 || index >= BYM2) {
      throw cms::Exception("DataCorrupt") << "SiPixelTemplate2D::ysigma2 called with index = " << index << std::endl;
   }
#else
   assert(index > 0 && index < BYM2);
#endif
   
   // Define the maximum signal to use in the parameterization
   
   // Evaluate pixel-by-pixel uncertainties (weights) for the templ analysis
   
   if(qpixel < sxymax_) {
      sigi = qpixel;
      qscale = 1.f;
   } else {
      sigi = sxymax_;
      qscale = qpixel/sxymax_;
   }
   sigi2 = sigi*sigi; sigi3 = sigi2*sigi; sigi4 = sigi3*sigi;
   if(index <= T2HYP1) {
      err00 = xypary0x0_[0][0]+xypary0x0_[0][1]*sigi+xypary0x0_[0][2]*sigi2+xypary0x0_[0][3]*sigi3+xypary0x0_[0][4]*sigi4;
      err2 = err00
    +adcota_*(xypary0x1_[0][0]+xypary0x1_[0][1]*sigi+xypary0x1_[0][2]*sigi2+xypary0x1_[0][3]*sigi3+xypary0x1_[0][4]*sigi4 - err00)
    +adcotb_*(xypary1x0_[0][0]+xypary1x0_[0][1]*sigi+xypary1x0_[0][2]*sigi2+xypary1x0_[0][3]*sigi3+xypary1x0_[0][4]*sigi4 - err00);
   } else {
      err00 = xypary0x0_[1][0]+xypary0x0_[1][1]*sigi+xypary0x0_[1][2]*sigi2+xypary0x0_[1][3]*sigi3+xypary0x0_[1][4]*sigi4;
      err2 = err00
    +adcota_*(xypary0x1_[1][0]+xypary0x1_[1][1]*sigi+xypary0x1_[1][2]*sigi2+xypary0x1_[1][3]*sigi3+xypary0x1_[1][4]*sigi4 - err00)
    +adcotb_*(xypary1x0_[1][0]+xypary1x0_[1][1]*sigi+xypary1x0_[1][2]*sigi2+xypary1x0_[1][3]*sigi3+xypary1x0_[1][4]*sigi4 - err00);
   }
   xysig2 =qscale*err2;
   if(xysig2 <= 0.f) {xysig2 = s50_*s50_;}
   
   return;
   
} // End xysigma2

bool SiPixelTemplate2D::xytemp	(	float	xhit,
		float	yhit,
		bool	ydouble[21+2],
		bool	xdouble[13+2],
		float	template2d[13+2][21+2],
		bool	derivatives,
		float	dpdx2d[2][13+2][21+2],
		float &	QTemplate
	)

Parameters

xhit	- (input) x-position of hit relative to the lower left corner of pixel[1][1] (to allow for the "padding" of the two-d clusters in the splitter)
yhit	- (input) y-position of hit relative to the lower left corner of pixel[1][1]
ydouble	- (input) STL vector of 21 element array to flag a double-pixel starting at cluster[1][1]
xdouble	- (input) STL vector of 11 element array to flag a double-pixel starting at cluster[1][1]
template2d	- (output) 2d template of size matched to the cluster. Input must be zeroed since charge is added only.

Definition at line 958 of file SiPixelTemplate2D.cc.

References BXM2, BXM3, BYM2, BYM3, PVValHelper::dx, PVValHelper::dy, f, objects.autophobj::float, mps_fire::i, createfilelist::int, gen::k, reco::ParticleMasses::k0, funct::m, SiStripPI::max, min(), gen::n, T2HX, T2HY, T2XSIZE, and T2YSIZE.

Referenced by SiPixelTemplateReco2D::PixelTempReco2D(), SiPixelDigitizerAlgorithm::PixelTempRewgt2D(), SiPixelTemplateSplit::PixelTempSplit(), and xytemp().

{
   // Interpolate for a new set of track angles
   
   // Local variables
   int pixx, pixy, k0, k1, l0, l1, deltax, deltay, iflipy, jflipx, imin, imax, jmin, jmax;
   int m, n;
   float dx, dy, ddx, ddy, adx, ady;
   //   const float deltaxy[2] = {8.33f, 12.5f};
   const float deltaxy[2] = {16.67f, 25.0f};
   
   // Check to see if interpolation is valid
   
   
   // next, determine the indices of the closest point in k (y-displacement), l (x-displacement)
   // pixy and pixx are the indices of the struck pixel in the (Ty,Tx) system
   // k0,k1 are the k-indices of the closest and next closest point
   // l0,l1 are the l-indices of the closest and next closest point
   
   pixy = (int)floorf(yhit/ysize_);
   dy = yhit-(pixy+0.5f)*ysize_;
   if(flip_y_) {dy = -dy;}
   k0 = (int)(dy/ysize_*6.f+3.5f);
   if(k0 < 0) k0 = 0;
   if(k0 > 6) k0 = 6;
   ddy = 6.f*dy/ysize_ - (k0-3);
   ady = fabs(ddy);
   if(ddy > 0.f) {k1 = k0 + 1; if(k1 > 6) k1 = k0-1;} else {k1 = k0 - 1; if(k1 < 0) k1 = k0+1;}
   pixx = (int)floorf(xhit/xsize_);
   dx = xhit-(pixx+0.5f)*xsize_;
   if(flip_x_) {dx = -dx;}
   l0 = (int)(dx/xsize_*6.f+3.5f);
   if(l0 < 0) l0 = 0;
   if(l0 > 6) l0 = 6;
   ddx = 6.f*dx/xsize_ - (l0-3);
   adx = fabs(ddx);
   if(ddx > 0.f) {l1 = l0 + 1; if(l1 > 6) l1 = l0-1;} else {l1 = l0 - 1; if(l1 < 0) l1 = l0+1;}
   
   // OK, lets do the template interpolation.
   
   // First find the limits of the indices for non-zero pixels
   
   imin = std::min(entry00_->iymin,entry10_->iymin);
   imin_ = std::min(imin,entry01_->iymin);
   
   jmin = std::min(entry00_->jxmin,entry10_->jxmin);
   jmin_ = std::min(jmin,entry01_->jxmin);
   
   imax = std::max(entry00_->iymax,entry10_->iymax);
   imax_ = std::max(imax,entry01_->iymax);
   
   jmax = std::max(entry00_->jxmax,entry10_->jxmax);
   jmax_ = std::max(jmax,entry01_->jxmax);
   
   // Calculate the x and y offsets to make the new template
   
   // First, shift the struck pixel coordinates to the (Ty+2, Tx+2) system
   
   ++pixy; ++pixx;
   
   // In the template store, the struck pixel is always (THy,THx)
   
   deltax = pixx - T2HX;
   deltay = pixy - T2HY;
   
   //  First zero the local 2-d template
   
   for(int j=0; j<BXM2; ++j) {for(int i=0; i<BYM2; ++i) {xytemp_[j][i] = 0.f;}}
   
   // Loop over the non-zero part of the template index space and interpolate
   
   for(int j=jmin_; j<=jmax_; ++j) {
      // Flip indices as needed
      if(flip_x_) {jflipx=T2XSIZE-1-j; m = deltax+jflipx;} else {m = deltax+j;}
      for(int i=imin_; i<=imax_; ++i) {
         if(flip_y_) {iflipy=T2YSIZE-1-i; n = deltay+iflipy;} else {n = deltay+i;}
         if(m>=0 && m<=BXM3 && n>=0 && n<=BYM3) {
            xytemp_[m][n] = (float)entry00_->xytemp[k0][l0][i][j]
            + adx*(float)(entry00_->xytemp[k0][l1][i][j] - entry00_->xytemp[k0][l0][i][j])
            + ady*(float)(entry00_->xytemp[k1][l0][i][j] - entry00_->xytemp[k0][l0][i][j])
            + adcota_*(float)(entry01_->xytemp[k0][l0][i][j] - entry00_->xytemp[k0][l0][i][j])
            + adcotb_*(float)(entry10_->xytemp[k0][l0][i][j] - entry00_->xytemp[k0][l0][i][j]);
         }
      }
   }
   
   //combine rows and columns to simulate double pixels
   
   for(int n=1; n<BYM3; ++n) {
      if(ydouble[n]) {
         //  Combine the y-columns
         for(int m=1; m<BXM3; ++m) {
            xytemp_[m][n] += xytemp_[m][n+1];
         }
         //  Now shift the remaining pixels over by one column
         for(int i=n+1; i<BYM3; ++i) {
            for(int m=1; m<BXM3; ++m) {
               xytemp_[m][i] = xytemp_[m][i+1];
            }
         }
      }
   }
   
   //combine rows and columns to simulate double pixels
   
   for(int m=1; m<BXM3; ++m) {
      if(xdouble[m]) {
         //  Combine the x-rows
         for(int n=1; n<BYM3; ++n) {
            xytemp_[m][n] += xytemp_[m+1][n];
         }
         //  Now shift the remaining pixels over by one row
         for(int j=m+1; j<BXM3; ++j) {
            for(n=1; n<BYM3; ++n) {
               xytemp_[j][n] = xytemp_[j+1][n];
            }
         }
      }
   }
   
   //  Finally, loop through and increment the external template
   
   float qtemptot = 0.f;
   
   for(int n=1; n<BYM3; ++n) {
      for(int m=1; m<BXM3; ++m) {
         if(xytemp_[m][n] != 0.f) {template2d[m][n] += xytemp_[m][n]; qtemptot += xytemp_[m][n];}
      }
   }
   
   QTemplate = qtemptot;
   
   if(derivatives) {
      
      float dxytempdx[2][BXM2][BYM2], dxytempdy[2][BXM2][BYM2];
      
      for(int k = 0; k<2; ++k) {
         for(int i = 0; i<BXM2; ++i) {
            for(int j = 0; j<BYM2; ++j) {
               dxytempdx[k][i][j] = 0.f;
               dxytempdy[k][i][j] = 0.f;
               dpdx2d[k][i][j] = 0.f;
            }
         }
      }
      
      // First do shifted +x template
      
      pixx = (int)floorf((xhit+deltaxy[0])/xsize_);
      dx = (xhit+deltaxy[0])-(pixx+0.5f)*xsize_;
      if(flip_x_) {dx = -dx;}
      l0 = (int)(dx/xsize_*6.f+3.5f);
      if(l0 < 0) l0 = 0;
      if(l0 > 6) l0 = 6;
      ddx = 6.f*dx/xsize_ - (l0-3);
      adx = fabs(ddx);
      if(ddx > 0.f) {l1 = l0 + 1; if(l1 > 6) l1 = l0-1;} else {l1 = l0 - 1; if(l1 < 0) l1 = l0+1;}
      
      // OK, lets do the template interpolation.
      
      // Calculate the x and y offsets to make the new template
      
      // First, shift the struck pixel coordinates to the (Ty+2, Tx+2) system
      
      ++pixx;
      
      // In the template store, the struck pixel is always (THy,THx)
      
      deltax = pixx - T2HX;
      
      // Loop over the non-zero part of the template index space and interpolate
      
      for(int j=jmin_; j<=jmax_; ++j) {
         // Flip indices as needed
         if(flip_x_) {jflipx=T2XSIZE-1-j; m = deltax+jflipx;} else {m = deltax+j;}
         for(int i=imin_; i<=imax_; ++i) {
            if(flip_y_) {iflipy=T2YSIZE-1-i; n = deltay+iflipy;} else {n = deltay+i;}
            if(m>=0 && m<=BXM3 && n>=0 && n<=BYM3) {
               dxytempdx[1][m][n] = (float)entry00_->xytemp[k0][l0][i][j]
               + adx*(float)(entry00_->xytemp[k0][l1][i][j] -   entry00_->xytemp[k0][l0][i][j])
               + ady*(float)(entry00_->xytemp[k1][l0][i][j] - entry00_->xytemp[k0][l0][i][j])
               + adcota_*(float)(entry01_->xytemp[k0][l0][i][j] - entry00_->xytemp[k0][l0][i][j])
               + adcotb_*(float)(entry10_->xytemp[k0][l0][i][j] - entry00_->xytemp[k0][l0][i][j]);
            }
         }
      }
      
      //combine rows and columns to simulate double pixels
      
      for(int n=1; n<BYM3; ++n) {
         if(ydouble[n]) {
            //  Combine the y-columns
            for(int m=1; m<BXM3; ++m) {
               dxytempdx[1][m][n] += dxytempdx[1][m][n+1];
            }
            //  Now shift the remaining pixels over by one column
            for(int i=n+1; i<BYM3; ++i) {
               for(int m=1; m<BXM3; ++m) {
                  dxytempdx[1][m][i] = dxytempdx[1][m][i+1];
               }
            }
         }
      }
      
      //combine rows and columns to simulate double pixels
      
      for(int m=1; m<BXM3; ++m) {
         if(xdouble[m]) {
            //  Combine the x-rows
            for(int n=1; n<BYM3; ++n) {
               dxytempdx[1][m][n] += dxytempdx[1][m+1][n];
            }
            //  Now shift the remaining pixels over by one row
            for(int j=m+1; j<BXM3; ++j) {
               for(int n=1; n<BYM3; ++n) {
                  dxytempdx[1][j][n] = dxytempdx[1][j+1][n];
               }
            }
         }
      }
      
      // Next do shifted -x template
      
      pixx = (int)floorf((xhit-deltaxy[0])/xsize_);
      dx = (xhit-deltaxy[0])-(pixx+0.5f)*xsize_;
      if(flip_x_) {dx = -dx;}
      l0 = (int)(dx/xsize_*6.f+3.5f);
      if(l0 < 0) l0 = 0;
      if(l0 > 6) l0 = 6;
      ddx = 6.f*dx/xsize_ - (l0-3);
      adx = fabs(ddx);
      if(ddx > 0.f) {l1 = l0 + 1; if(l1 > 6) l1 = l0-1;} else {l1 = l0 - 1; if(l1 < 0) l1 = l0+1;}
      
      // OK, lets do the template interpolation.
      
      // Calculate the x and y offsets to make the new template
      
      // First, shift the struck pixel coordinates to the (Ty+2, Tx+2) system
      
      ++pixx;
      
      // In the template store, the struck pixel is always (THy,THx)
      
      deltax = pixx - T2HX;
      
      // Loop over the non-zero part of the template index space and interpolate
      
      for(int j=jmin_; j<=jmax_; ++j) {
         // Flip indices as needed
         if(flip_x_) {jflipx=T2XSIZE-1-j; m = deltax+jflipx;} else {m = deltax+j;}
         for(int i=imin_; i<=imax_; ++i) {
            if(flip_y_) {iflipy=T2YSIZE-1-i; n = deltay+iflipy;} else {n = deltay+i;}
            if(m>=0 && m<=BXM3 && n>=0 && n<=BYM3) {
               dxytempdx[0][m][n]  = (float)entry00_->xytemp[k0][l0][i][j]
               + adx*(float)(entry00_->xytemp[k0][l1][i][j] -   entry00_->xytemp[k0][l0][i][j])
               + ady*(float)(entry00_->xytemp[k1][l0][i][j] - entry00_->xytemp[k0][l0][i][j])
               + adcota_*(float)(entry01_->xytemp[k0][l0][i][j] - entry00_->xytemp[k0][l0][i][j])
               + adcotb_*(float)(entry10_->xytemp[k0][l0][i][j] - entry00_->xytemp[k0][l0][i][j]);
            }
         }
      }
      
      //combine rows and columns to simulate double pixels
      
      for(int n=1; n<BYM3; ++n) {
         if(ydouble[n]) {
            //  Combine the y-columns
            for(int m=1; m<BXM3; ++m) {
               dxytempdx[0][m][n] += dxytempdx[0][m][n+1];
            }
            //  Now shift the remaining pixels over by one column
            for(int i=n+1; i<BYM3; ++i) {
               for(int m=1; m<BXM3; ++m) {
                  dxytempdx[0][m][i] = dxytempdx[0][m][i+1];
               }
            }
         }
      }
      
      //combine rows and columns to simulate double pixels
      
      for(int m=1; m<BXM3; ++m) {
         if(xdouble[m]) {
            //  Combine the x-rows
            for(int n=1; n<BYM3; ++n) {
               dxytempdx[0][m][n] += dxytempdx[0][m+1][n];
            }
            //  Now shift the remaining pixels over by one row
            for(int j=m+1; j<BXM3; ++j) {
               for(int n=1; n<BYM3; ++n) {
                  dxytempdx[0][j][n] = dxytempdx[0][j+1][n];
               }
            }
         }
      }
      
      
      //  Finally, normalize the derivatives and copy the results to the output array
      
      for(int n=1; n<BYM3; ++n) {
         for(int m=1; m<BXM3; ++m) {
            dpdx2d[0][m][n] = (dxytempdx[1][m][n] - dxytempdx[0][m][n])/(2.*deltaxy[0]);
         }
      }
      
      // Next, do shifted y template
      
      pixy = (int)floorf((yhit+deltaxy[1])/ysize_);
      dy = (yhit+deltaxy[1])-(pixy+0.5f)*ysize_;
      if(flip_y_) {dy = -dy;}
      k0 = (int)(dy/ysize_*6.f+3.5f);
      if(k0 < 0) k0 = 0;
      if(k0 > 6) k0 = 6;
      ddy = 6.f*dy/ysize_ - (k0-3);
      ady = fabs(ddy);
      if(ddy > 0.f) {k1 = k0 + 1; if(k1 > 6) k1 = k0-1;} else {k1 = k0 - 1; if(k1 < 0) k1 = k0+1;}
      pixx = (int)floorf(xhit/xsize_);
      dx = xhit-(pixx+0.5f)*xsize_;
      if(flip_x_) {dx = -dx;}
      l0 = (int)(dx/xsize_*6.f+3.5f);
      if(l0 < 0) l0 = 0;
      if(l0 > 6) l0 = 6;
      ddx = 6.f*dx/xsize_ - (l0-3);
      adx = fabs(ddx);
      if(ddx > 0.f) {l1 = l0 + 1; if(l1 > 6) l1 = l0-1;} else {l1 = l0 - 1; if(l1 < 0) l1 = l0+1;}
      
      // OK, lets do the template interpolation.
      
      // Calculate the x and y offsets to make the new template
      
      // First, shift the struck pixel coordinates to the (Ty+2, Tx+2) system
      
      ++pixy; ++pixx;
      
      // In the template store, the struck pixel is always (THy,THx)
      
      deltax = pixx - T2HX;
      deltay = pixy - T2HY;
      
      // Loop over the non-zero part of the template index space and interpolate
      
      for(int j=jmin_; j<=jmax_; ++j) {
         // Flip indices as needed
         if(flip_x_) {jflipx=T2XSIZE-1-j; m = deltax+jflipx;} else {m = deltax+j;}
         for(int i=imin_; i<=imax_; ++i) {
            if(flip_y_) {iflipy=T2YSIZE-1-i; n = deltay+iflipy;} else {n = deltay+i;}
            if(m>=0 && m<=BXM3 && n>=0 && n<=BYM3) {
               dxytempdy[1][m][n] = (float)entry00_->xytemp[k0][l0][i][j]
               + adx*(float)(entry00_->xytemp[k0][l1][i][j] - entry00_->xytemp[k0][l0][i][j])
               + ady*(float)(entry00_->xytemp[k1][l0][i][j] - entry00_->xytemp[k0][l0][i][j])
               + adcota_*(float)(entry01_->xytemp[k0][l0][i][j] - entry00_->xytemp[k0][l0][i][j])
               + adcotb_*(float)(entry10_->xytemp[k0][l0][i][j] - entry00_->xytemp[k0][l0][i][j]);
            }
         }
      }
      
      //combine rows and columns to simulate double pixels
      
      for(int n=1; n<BYM3; ++n) {
         if(ydouble[n]) {
            //  Combine the y-columns
            for(int m=1; m<BXM3; ++m) {
               dxytempdy[1][m][n] += dxytempdy[1][m][n+1];
            }
            //  Now shift the remaining pixels over by one column
            for(int i=n+1; i<BYM3; ++i) {
               for(int m=1; m<BXM3; ++m) {
                  dxytempdy[1][m][i] = dxytempdy[1][m][i+1];
               }
            }
         }
      }
      
      //combine rows and columns to simulate double pixels
      
      for(int m=1; m<BXM3; ++m) {
         if(xdouble[m]) {
            //  Combine the x-rows
            for(int n=1; n<BYM3; ++n) {
               dxytempdy[1][m][n] += dxytempdy[1][m+1][n];
            }
            //  Now shift the remaining pixels over by one row
            for(int j=m+1; j<BXM3; ++j) {
               for(int n=1; n<BYM3; ++n) {
                  dxytempdy[1][j][n] = dxytempdy[1][j+1][n];
               }
            }
         }
      }
      
      // Next, do shifted -y template
      
      pixy = (int)floorf((yhit-deltaxy[1])/ysize_);
      dy = (yhit-deltaxy[1])-(pixy+0.5f)*ysize_;
      if(flip_y_) {dy = -dy;}
      k0 = (int)(dy/ysize_*6.f+3.5f);
      if(k0 < 0) k0 = 0;
      if(k0 > 6) k0 = 6;
      ddy = 6.f*dy/ysize_ - (k0-3);
      ady = fabs(ddy);
      if(ddy > 0.f) {k1 = k0 + 1; if(k1 > 6) k1 = k0-1;} else {k1 = k0 - 1; if(k1 < 0) k1 = k0+1;}
      
      // OK, lets do the template interpolation.
      
      // Calculate the x and y offsets to make the new template
      
      // First, shift the struck pixel coordinates to the (Ty+2, Tx+2) system
      
      ++pixy;
      
      // In the template store, the struck pixel is always (THy,THx)
      
      deltay = pixy - T2HY;
      
      // Loop over the non-zero part of the template index space and interpolate
      
      for(int j=jmin_; j<=jmax_; ++j) {
         // Flip indices as needed
         if(flip_x_) {jflipx=T2XSIZE-1-j; m = deltax+jflipx;} else {m = deltax+j;}
         for(int i=imin_; i<=imax_; ++i) {
            if(flip_y_) {iflipy=T2YSIZE-1-i; n = deltay+iflipy;} else {n = deltay+i;}
            if(m>=0 && m<=BXM3 && n>=0 && n<=BYM3) {
               dxytempdy[0][m][n] = (float)entry00_->xytemp[k0][l0][i][j]
               + adx*(float)(entry00_->xytemp[k0][l1][i][j] - entry00_->xytemp[k0][l0][i][j])
               + ady*(float)(entry00_->xytemp[k1][l0][i][j] - entry00_->xytemp[k0][l0][i][j])
               + adcota_*(float)(entry01_->xytemp[k0][l0][i][j] - entry00_->xytemp[k0][l0][i][j])
               + adcotb_*(float)(entry10_->xytemp[k0][l0][i][j] - entry00_->xytemp[k0][l0][i][j]);
            }
         }
      }
      
      //combine rows and columns to simulate double pixels
      
      for(int n=1; n<BYM3; ++n) {
         if(ydouble[n]) {
            //  Combine the y-columns
            for(int m=1; m<BXM3; ++m) {
               dxytempdy[0][m][n] += dxytempdy[0][m][n+1];
            }
            //  Now shift the remaining pixels over by one column
            for(int i=n+1; i<BYM3; ++i) {
               for(int m=1; m<BXM3; ++m) {
                  dxytempdy[0][m][i] = dxytempdy[0][m][i+1];
               }
            }
         }
      }
      
      //combine rows and columns to simulate double pixels
      
      for(int m=1; m<BXM3; ++m) {
         if(xdouble[m]) {
            //  Combine the x-rows
            for(int n=1; n<BYM3; ++n) {
               dxytempdy[0][m][n] += dxytempdy[0][m+1][n];
            }
            //  Now shift the remaining pixels over by one row
            for(int j=m+1; j<BXM3; ++j) {
               for(int n=1; n<BYM3; ++n) {
                  dxytempdy[0][j][n] = dxytempdy[0][j+1][n];
               }
            }
         }
      }
      
      //  Finally, normalize the derivatives and copy the results to the output array
      
      for(int n=1; n<BYM3; ++n) {
         for(int m=1; m<BXM3; ++m) {
            dpdx2d[1][m][n] = (dxytempdy[1][m][n] - dxytempdy[0][m][n])/(2.*deltaxy[1]);
         }
      }
   }
   
   return success_;
} // xytemp

bool SiPixelTemplate2D::xytemp	(	float	xhit,
		float	yhit,
		bool	ydouble[21+2],
		bool	xdouble[13+2],
		float	template2d[13+2][21+2]
	)

Interpolate stored 2-D information for input angles and hit position to make a 2-D template

Parameters

xhit	- (input) x-position of hit relative to the lower left corner of pixel[1][1] (to allow for the "padding" of the two-d clusters in the splitter)
yhit	- (input) y-position of hit relative to the lower left corner of pixel[1][1]
ydouble	- (input) STL vector of 21 element array to flag a double-pixel starting at cluster[1][1]
xdouble	- (input) STL vector of 11 element array to flag a double-pixel starting at cluster[1][1]
template2d	- (output) 2d template of size matched to the cluster. Input must be zeroed since charge is added only.

Definition at line 1446 of file SiPixelTemplate2D.cc.

References BXM2, BYM2, and xytemp().

{
   // Interpolate for a new set of track angles
   
   bool derivatives = false;
   float dpdx2d[2][BXM2][BYM2];
   float QTemplate;
   
   return SiPixelTemplate2D::xytemp(xhit, yhit, ydouble, xdouble, template2d, derivatives, dpdx2d, QTemplate);
   
} // xytemp

bool SiPixelTemplate2D::xytemp	(	int	id,
		float	cotalpha,
		float	cotbeta,
		float	xhit,
		float	yhit,
		std::vector< bool > &	ydouble,
		std::vector< bool > &	xdouble,
		float	template2d[13+2][21+2]
	)

Interpolate stored 2-D information for input angles and hit position to make a 2-D template

Parameters

id	- (input) the id of the template
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)
xhit	- (input) x-position of hit relative to the lower left corner of pixel[1][1] (to allow for the "padding" of the two-d clusters in the splitter)
yhit	- (input) y-position of hit relative to the lower left corner of pixel[1][1]
ydouble	- (input) STL vector of 21 element array to flag a double-pixel starting at cluster[1][1]
xdouble	- (input) STL vector of 11 element array to flag a double-pixel starting at cluster[1][1]
template2d	- (output) 2d template of size matched to the cluster. Input must be zeroed since charge is added only.

Definition at line 1472 of file SiPixelTemplate2D.cc.

References BXM2, BYM2, f, mps_fire::i, interpolate(), TXSIZE, TYSIZE, and xytemp().

{
   
   // Local variables
   
   bool derivatives = false;
   float dpdx2d[2][BXM2][BYM2];
   float QTemplate;
   float locBx = 1.f;
   if(cotbeta < 0.f) {locBx = -1.f;}
   float locBz = locBx;
   if(cotalpha < 0.f) {locBz = -locBx;}
   
   bool yd[BYM2], xd[BXM2];
   
   yd[0] = false; yd[BYM2-1] = false;
   for(int i=0; i<TYSIZE; ++i) { yd[i+1] = ydouble[i];}
   xd[0] = false; xd[BXM2-1] = false;
   for(int j=0; j<TXSIZE; ++j) { xd[j+1] = xdouble[j];}
   
   
   // Interpolate for a new set of track angles
   
   if(SiPixelTemplate2D::interpolate(id, cotalpha, cotbeta, locBz, locBx)) {
      return SiPixelTemplate2D::xytemp(xhit, yhit, yd, xd, template2d, derivatives, dpdx2d, QTemplate);
   } else {
      return false;
   }
   
} // xytemp

float SiPixelTemplate2D::ysize ( )

inline

pixel y-size (microns)

Definition at line 209 of file SiPixelTemplate2D.h.

Referenced by SiPixelTemplateReco2D::PixelTempReco2D(), and SiPixelDigitizerAlgorithm::PixelTempRewgt2D().

float SiPixelTemplate2D::zsize ( )

inline

pixel z-size or thickness (microns)

Definition at line 210 of file SiPixelTemplate2D.h.

Member Data Documentation

float SiPixelTemplate2D::adcota_

private

fractional pixel distance of cot(alpha) from jx0_

Definition at line 247 of file SiPixelTemplate2D.h.

float SiPixelTemplate2D::adcotb_

private

fractional pixel distance of cot(beta) from iy0_

Definition at line 244 of file SiPixelTemplate2D.h.

float SiPixelTemplate2D::chi2avgone_

private

average chi^2 for 1 pixel clusters

Definition at line 271 of file SiPixelTemplate2D.h.

float SiPixelTemplate2D::chi2minone_

private

minimum of chi^2 for 1 pixel clusters

Definition at line 272 of file SiPixelTemplate2D.h.

float SiPixelTemplate2D::chi2ppix_

private

average chi^2 per struck pixel

Definition at line 269 of file SiPixelTemplate2D.h.

float SiPixelTemplate2D::chi2scale_

private

scale factor for chi2 distribution

Definition at line 270 of file SiPixelTemplate2D.h.

float SiPixelTemplate2D::clslenx_

private

projected x-length of cluster

Definition at line 274 of file SiPixelTemplate2D.h.

float SiPixelTemplate2D::clsleny_

private

projected y-length of cluster

Definition at line 273 of file SiPixelTemplate2D.h.

float SiPixelTemplate2D::cota_current_

private

current cot alpha

Definition at line 231 of file SiPixelTemplate2D.h.

float SiPixelTemplate2D::cotalpha0_

private

minimum cot(alpha) covered

Definition at line 239 of file SiPixelTemplate2D.h.

float SiPixelTemplate2D::cotalpha1_

private

maximum cot(alpha) covered

Definition at line 240 of file SiPixelTemplate2D.h.

float SiPixelTemplate2D::cotb_current_

private

current cot beta

Definition at line 232 of file SiPixelTemplate2D.h.

float SiPixelTemplate2D::cotbeta0_

private

minimum cot(beta) covered

Definition at line 236 of file SiPixelTemplate2D.h.

float SiPixelTemplate2D::cotbeta1_

private

maximum cot(beta) covered

Definition at line 237 of file SiPixelTemplate2D.h.

float SiPixelTemplate2D::deltacota_

private

cot(alpha) bin size

Definition at line 241 of file SiPixelTemplate2D.h.

float SiPixelTemplate2D::deltacotb_

private

cot(beta) bin size

Definition at line 238 of file SiPixelTemplate2D.h.

float SiPixelTemplate2D::delyavg_

private

average difference between clsleny_ and cluster length [with threshold effects]

Definition at line 277 of file SiPixelTemplate2D.h.

float SiPixelTemplate2D::delysig_

private

rms of difference between clsleny_ and cluster length [with threshold effects]

Definition at line 278 of file SiPixelTemplate2D.h.

int SiPixelTemplate2D::Dtype_

private

flags BPix (=0) or FPix (=1)

Definition at line 235 of file SiPixelTemplate2D.h.

const SiPixelTemplateEntry2D* SiPixelTemplate2D::entry00_

private

Definition at line 294 of file SiPixelTemplate2D.h.

const SiPixelTemplateEntry2D* SiPixelTemplate2D::entry01_

private

Definition at line 296 of file SiPixelTemplate2D.h.

const SiPixelTemplateEntry2D* SiPixelTemplate2D::entry10_

private

Definition at line 295 of file SiPixelTemplate2D.h.

float SiPixelTemplate2D::fbin_[3]

private

The QBin definitions in Q_clus/Q_avg.

Definition at line 293 of file SiPixelTemplate2D.h.

bool SiPixelTemplate2D::flip_x_

private

flip x sign-sensitive quantities

Definition at line 253 of file SiPixelTemplate2D.h.

bool SiPixelTemplate2D::flip_y_

private

flip y sign-sensitive quantities

Definition at line 252 of file SiPixelTemplate2D.h.

int SiPixelTemplate2D::id_current_

private

current id

Definition at line 229 of file SiPixelTemplate2D.h.

int SiPixelTemplate2D::imax_

private

max y index of templated cluster

Definition at line 249 of file SiPixelTemplate2D.h.

int SiPixelTemplate2D::imin_

private

min y index of templated cluster

Definition at line 248 of file SiPixelTemplate2D.h.

int SiPixelTemplate2D::index_id_

private

current index

Definition at line 230 of file SiPixelTemplate2D.h.

int SiPixelTemplate2D::iy0_

private

index of nearest cot(beta) bin

Definition at line 242 of file SiPixelTemplate2D.h.

int SiPixelTemplate2D::iy1_

private

index of next-nearest cot(beta) bin

Definition at line 243 of file SiPixelTemplate2D.h.

int SiPixelTemplate2D::jmax_

private

max x index of templated cluster

Definition at line 251 of file SiPixelTemplate2D.h.

int SiPixelTemplate2D::jmin_

private

min x index of templated cluster

Definition at line 250 of file SiPixelTemplate2D.h.

int SiPixelTemplate2D::jx0_

private

index of nearest cot(alpha) bin

Definition at line 245 of file SiPixelTemplate2D.h.

int SiPixelTemplate2D::jx1_

private

index of next-nearest cot(alpha) bin

Definition at line 246 of file SiPixelTemplate2D.h.

float SiPixelTemplate2D::kappavav_

private

kappa parameter in Vavilov distribution

Definition at line 285 of file SiPixelTemplate2D.h.

float SiPixelTemplate2D::lanpar_[2][5]

private

Interpolated Landau parameters.

Definition at line 268 of file SiPixelTemplate2D.h.

float SiPixelTemplate2D::lorxdrift_

private

Lorentz x-drift.

Definition at line 289 of file SiPixelTemplate2D.h.

float SiPixelTemplate2D::lorxwidth_

private

Lorentz x-width.

Definition at line 287 of file SiPixelTemplate2D.h.

float SiPixelTemplate2D::lorydrift_

private

Lorentz y-drift.

Definition at line 288 of file SiPixelTemplate2D.h.

float SiPixelTemplate2D::lorywidth_

private

Lorentz y-width (sign corrected for fpix frame)

Definition at line 286 of file SiPixelTemplate2D.h.

float SiPixelTemplate2D::mpvvav_

private

most probable Q in Vavilov distribution

Definition at line 283 of file SiPixelTemplate2D.h.

int SiPixelTemplate2D::Nxx_

private

number of cot(alpha)-entries (rows) in template

Definition at line 234 of file SiPixelTemplate2D.h.

int SiPixelTemplate2D::Nyx_

private

number of cot(beta)-entries (columns) in template

Definition at line 233 of file SiPixelTemplate2D.h.

float SiPixelTemplate2D::offsetx_[4]

private

x-offset in charge bins

Definition at line 281 of file SiPixelTemplate2D.h.

float SiPixelTemplate2D::offsety_[4]

private

y-offset in charge bins

Definition at line 282 of file SiPixelTemplate2D.h.

float SiPixelTemplate2D::pixmax_

private

maximum pixel charge

Definition at line 260 of file SiPixelTemplate2D.h.

float SiPixelTemplate2D::qavg_

private

average cluster charge for this set of track angles

Definition at line 259 of file SiPixelTemplate2D.h.

float SiPixelTemplate2D::qscale_

private

charge scaling factor

Definition at line 261 of file SiPixelTemplate2D.h.

float SiPixelTemplate2D::s50_

private

1/2 of the pixel threshold signal in adc units

Definition at line 262 of file SiPixelTemplate2D.h.

float SiPixelTemplate2D::scalex_[4]

private

x-error scale factor in charge bins

Definition at line 279 of file SiPixelTemplate2D.h.

float SiPixelTemplate2D::scalexavg_

private

average x-error scale factor

Definition at line 275 of file SiPixelTemplate2D.h.

float SiPixelTemplate2D::scaley_[4]

private

y-error scale factor in charge bins

Definition at line 280 of file SiPixelTemplate2D.h.

float SiPixelTemplate2D::scaleyavg_

private

average y-error scale factor

Definition at line 276 of file SiPixelTemplate2D.h.

float SiPixelTemplate2D::sigmavav_

private

scale factor in Vavilov distribution

Definition at line 284 of file SiPixelTemplate2D.h.

bool SiPixelTemplate2D::success_

private

true if cotalpha, cotbeta are inside of the acceptance (dynamically loaded)

Definition at line 254 of file SiPixelTemplate2D.h.

float SiPixelTemplate2D::sxymax_

private

average pixel signal for y-projection of cluster

Definition at line 263 of file SiPixelTemplate2D.h.

const std::vector< SiPixelTemplateStore2D >& SiPixelTemplate2D::thePixelTemp_

private

Definition at line 300 of file SiPixelTemplate2D.h.

float SiPixelTemplate2D::xsize_

private

Pixel x-size.

Definition at line 290 of file SiPixelTemplate2D.h.

float SiPixelTemplate2D::xypary0x0_[2][5]

private

Polynomial error parameterization at ix0,iy0.

Definition at line 265 of file SiPixelTemplate2D.h.

float SiPixelTemplate2D::xypary0x1_[2][5]

private

Polynomial error parameterization at ix1,iy0.

Definition at line 267 of file SiPixelTemplate2D.h.

float SiPixelTemplate2D::xypary1x0_[2][5]

private

Polynomial error parameterization at ix0,iy1.

Definition at line 266 of file SiPixelTemplate2D.h.

float SiPixelTemplate2D::xytemp_[ 13+2][ 21+2]

private

template for xy-reconstruction

Definition at line 264 of file SiPixelTemplate2D.h.

float SiPixelTemplate2D::ysize_

private

Pixel y-size.

Definition at line 291 of file SiPixelTemplate2D.h.

float SiPixelTemplate2D::zsize_

private

Pixel z-size (thickness)

Definition at line 292 of file SiPixelTemplate2D.h.