CMS 3D CMS Logo

List of all members | Public Types | Public Member Functions | Public Attributes | Static Public Attributes
SiPixelDigitizerAlgorithm::PixelEfficiencies Struct Reference

Public Types

enum  shiftEnumerator { FPixRocIdShift = 3, BPixRocIdShift = 6 }
 

Public Member Functions

void init_from_db (const edm::ESHandle< TrackerGeometry > &, const edm::ESHandle< SiPixelDynamicInefficiency > &)
 
bool matches (const DetId &, const DetId &, const std::vector< uint32_t > &)
 
 PixelEfficiencies (const edm::ParameterSet &conf, bool AddPixelInefficiency, int NumberOfBarrelLayers, int NumberOfEndcapDisks)
 

Public Attributes

std::map< uint32_t, double > ChipGeomFactors
 
std::map< uint32_t, double > ColGeomFactors
 
unsigned int FPixIndex
 
bool FromConfig
 
std::map< uint32_t, size_t > iPU
 
std::map< uint32_t, double > PixelGeomFactors
 
std::map< uint32_t, std::vector< double > > PixelGeomFactorsROCBigPixels
 
std::map< uint32_t, std::vector< double > > PixelGeomFactorsROCStdPixels
 
std::vector< double > pu_scale
 
double theInnerEfficiency_FPix [20]
 
double theInstLumiScaleFactor
 
std::vector< double > theLadderEfficiency_BPix [20]
 
std::vector< double > theModuleEfficiency_BPix [20]
 
double theOuterEfficiency_FPix [20]
 
double thePixelChipEfficiency [20]
 
double thePixelColEfficiency [20]
 
double thePixelEfficiency [20]
 
std::vector< std::vector< double > > thePUEfficiency
 

Static Public Attributes

static const int rocIdMaskBits = 0x1F
 

Detailed Description

Internal use only.

Definition at line 241 of file SiPixelDigitizerAlgorithm.h.

Member Enumeration Documentation

Constructor & Destructor Documentation

SiPixelDigitizerAlgorithm::PixelEfficiencies::PixelEfficiencies ( const edm::ParameterSet conf,
bool  AddPixelInefficiency,
int  NumberOfBarrelLayers,
int  NumberOfEndcapDisks 
)

Definition at line 398 of file SiPixelDigitizerAlgorithm.cc.

References relativeConstraints::empty, Exception, edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), mps_fire::i, SiPixelDigitizerAlgorithm::NumberOfBarrelLayers, SiPixelDigitizerAlgorithm::NumberOfEndcapDisks, and findQualityFiles::size.

398  {
399  // pixel inefficiency
400  // Don't use Hard coded values, read inefficiencies in from DB/python config or don't use any
401  int NumberOfTotLayers = NumberOfBarrelLayers + NumberOfEndcapDisks;
404  FromConfig =
405  conf.exists("thePixelColEfficiency_BPix1") && conf.exists("thePixelColEfficiency_BPix2") && conf.exists("thePixelColEfficiency_BPix3") &&
406  conf.exists("thePixelColEfficiency_FPix1") && conf.exists("thePixelColEfficiency_FPix2") &&
407  conf.exists("thePixelEfficiency_BPix1") && conf.exists("thePixelEfficiency_BPix2") && conf.exists("thePixelEfficiency_BPix3") &&
408  conf.exists("thePixelEfficiency_FPix1") && conf.exists("thePixelEfficiency_FPix2") &&
409  conf.exists("thePixelChipEfficiency_BPix1") && conf.exists("thePixelChipEfficiency_BPix2") && conf.exists("thePixelChipEfficiency_BPix3") &&
410  conf.exists("thePixelChipEfficiency_FPix1") && conf.exists("thePixelChipEfficiency_FPix2");
411  if (NumberOfBarrelLayers==3) FromConfig = FromConfig && conf.exists("theLadderEfficiency_BPix1") && conf.exists("theLadderEfficiency_BPix2") && conf.exists("theLadderEfficiency_BPix3") &&
412  conf.exists("theModuleEfficiency_BPix1") && conf.exists("theModuleEfficiency_BPix2") && conf.exists("theModuleEfficiency_BPix3") &&
413  conf.exists("thePUEfficiency_BPix1") && conf.exists("thePUEfficiency_BPix2") && conf.exists("thePUEfficiency_BPix3") &&
414  conf.exists("theInnerEfficiency_FPix1") && conf.exists("theInnerEfficiency_FPix2") &&
415  conf.exists("theOuterEfficiency_FPix1") && conf.exists("theOuterEfficiency_FPix2") &&
416  conf.exists("thePUEfficiency_FPix_Inner") && conf.exists("thePUEfficiency_FPix_Outer") &&
417  conf.exists("theInstLumiScaleFactor");
418  if (NumberOfBarrelLayers>=4) FromConfig = FromConfig && conf.exists("thePixelColEfficiency_BPix4") &&
419  conf.exists("thePixelEfficiency_BPix4") && conf.exists("thePixelChipEfficiency_BPix4");
420  if (NumberOfEndcapDisks>=3) FromConfig = FromConfig && conf.exists("thePixelColEfficiency_FPix4") &&
421  conf.exists("thePixelEfficiency_FPix3") && conf.exists("thePixelChipEfficiency_FPix3");
422  if (FromConfig) {
423  LogInfo ("PixelDigitizer ") <<"The PixelDigitizer inefficiency configuration is read from the config file.\n";
424  theInstLumiScaleFactor = conf.getParameter<double>("theInstLumiScaleFactor");
425  int i=0;
426  thePixelColEfficiency[i++] = conf.getParameter<double>("thePixelColEfficiency_BPix1");
427  thePixelColEfficiency[i++] = conf.getParameter<double>("thePixelColEfficiency_BPix2");
428  thePixelColEfficiency[i++] = conf.getParameter<double>("thePixelColEfficiency_BPix3");
429  if (NumberOfBarrelLayers>=4){thePixelColEfficiency[i++] = conf.getParameter<double>("thePixelColEfficiency_BPix4");}
430  //
431  i=0;
432  thePixelEfficiency[i++] = conf.getParameter<double>("thePixelEfficiency_BPix1");
433  thePixelEfficiency[i++] = conf.getParameter<double>("thePixelEfficiency_BPix2");
434  thePixelEfficiency[i++] = conf.getParameter<double>("thePixelEfficiency_BPix3");
435  if (NumberOfBarrelLayers>=4){thePixelEfficiency[i++] = conf.getParameter<double>("thePixelEfficiency_BPix4");}
436  //
437  i=0;
438  thePixelChipEfficiency[i++] = conf.getParameter<double>("thePixelChipEfficiency_BPix1");
439  thePixelChipEfficiency[i++] = conf.getParameter<double>("thePixelChipEfficiency_BPix2");
440  thePixelChipEfficiency[i++] = conf.getParameter<double>("thePixelChipEfficiency_BPix3");
441  if (NumberOfBarrelLayers>=4){thePixelChipEfficiency[i++] = conf.getParameter<double>("thePixelChipEfficiency_BPix4");}
442  //
443  if (NumberOfBarrelLayers==3){
444  i=0;
445  theLadderEfficiency_BPix[i++] = conf.getParameter<std::vector<double> >("theLadderEfficiency_BPix1");
446  theLadderEfficiency_BPix[i++] = conf.getParameter<std::vector<double> >("theLadderEfficiency_BPix2");
447  theLadderEfficiency_BPix[i++] = conf.getParameter<std::vector<double> >("theLadderEfficiency_BPix3");
448  if ( ((theLadderEfficiency_BPix[0].size()!=20) || (theLadderEfficiency_BPix[1].size()!=32) ||
449  (theLadderEfficiency_BPix[2].size()!=44)) && (NumberOfBarrelLayers==3) )
450  throw cms::Exception("Configuration") << "Wrong ladder number in efficiency config!";
451  //
452  i=0;
453  theModuleEfficiency_BPix[i++] = conf.getParameter<std::vector<double> >("theModuleEfficiency_BPix1");
454  theModuleEfficiency_BPix[i++] = conf.getParameter<std::vector<double> >("theModuleEfficiency_BPix2");
455  theModuleEfficiency_BPix[i++] = conf.getParameter<std::vector<double> >("theModuleEfficiency_BPix3");
456  if ( ((theModuleEfficiency_BPix[0].size()!=4) || (theModuleEfficiency_BPix[1].size()!=4) ||
457  (theModuleEfficiency_BPix[2].size()!=4)) && (NumberOfBarrelLayers==3) )
458  throw cms::Exception("Configuration") << "Wrong module number in efficiency config!";
459  //
460  thePUEfficiency.push_back(conf.getParameter<std::vector<double> >("thePUEfficiency_BPix1"));
461  thePUEfficiency.push_back(conf.getParameter<std::vector<double> >("thePUEfficiency_BPix2"));
462  thePUEfficiency.push_back(conf.getParameter<std::vector<double> >("thePUEfficiency_BPix3"));
463  if ( ((thePUEfficiency[0].empty()) || (thePUEfficiency[1].empty()) ||
464  (thePUEfficiency[2].empty())) && (NumberOfBarrelLayers==3) )
465  throw cms::Exception("Configuration") << "At least one PU efficiency (BPix) number is needed in efficiency config!";
466  }
467  // The next is needed for Phase2 Tracker studies
468  if (NumberOfBarrelLayers>=5){
469  if (NumberOfTotLayers>20){throw cms::Exception("Configuration") <<"SiPixelDigitizer was given more layers than it can handle";}
470  // For Phase2 tracker layers just set the outermost BPix inefficiency to 99.9% THESE VALUES ARE HARDCODED ALSO ELSEWHERE IN THIS FILE
471  for (int j=5 ; j<=NumberOfBarrelLayers ; j++){
472  thePixelColEfficiency[j-1]=0.999;
473  thePixelEfficiency[j-1]=0.999;
474  thePixelChipEfficiency[j-1]=0.999;
475  }
476  }
477  //
478  i=FPixIndex;
479  thePixelColEfficiency[i++] = conf.getParameter<double>("thePixelColEfficiency_FPix1");
480  thePixelColEfficiency[i++] = conf.getParameter<double>("thePixelColEfficiency_FPix2");
481  if (NumberOfEndcapDisks>=3){thePixelColEfficiency[i++] = conf.getParameter<double>("thePixelColEfficiency_FPix3");}
482  i=FPixIndex;
483  thePixelEfficiency[i++] = conf.getParameter<double>("thePixelEfficiency_FPix1");
484  thePixelEfficiency[i++] = conf.getParameter<double>("thePixelEfficiency_FPix2");
485  if (NumberOfEndcapDisks>=3){thePixelEfficiency[i++] = conf.getParameter<double>("thePixelEfficiency_FPix3");}
486  i=FPixIndex;
487  thePixelChipEfficiency[i++] = conf.getParameter<double>("thePixelChipEfficiency_FPix1");
488  thePixelChipEfficiency[i++] = conf.getParameter<double>("thePixelChipEfficiency_FPix2");
489  if (NumberOfEndcapDisks>=3){thePixelChipEfficiency[i++] = conf.getParameter<double>("thePixelChipEfficiency_FPix3");}
490  // The next is needed for Phase2 Tracker studies
491  if (NumberOfEndcapDisks>=4){
492  if (NumberOfTotLayers>20){throw cms::Exception("Configuration") <<"SiPixelDigitizer was given more layers than it can handle";}
493  // For Phase2 tracker layers just set the extra FPix disk inefficiency to 99.9% THESE VALUES ARE HARDCODED ALSO ELSEWHERE IN THIS FILE
494  for (int j=4+FPixIndex ; j<=NumberOfEndcapDisks+NumberOfBarrelLayers ; j++){
495  thePixelColEfficiency[j-1]=0.999;
496  thePixelEfficiency[j-1]=0.999;
497  thePixelChipEfficiency[j-1]=0.999;
498  }
499  }
500  //FPix Dynamic Inefficiency
501  if (NumberOfBarrelLayers==3){
502  i=FPixIndex;
503  theInnerEfficiency_FPix[i++] = conf.getParameter<double>("theInnerEfficiency_FPix1");
504  theInnerEfficiency_FPix[i++] = conf.getParameter<double>("theInnerEfficiency_FPix2");
505  i=FPixIndex;
506  theOuterEfficiency_FPix[i++] = conf.getParameter<double>("theOuterEfficiency_FPix1");
507  theOuterEfficiency_FPix[i++] = conf.getParameter<double>("theOuterEfficiency_FPix2");
508  thePUEfficiency.push_back(conf.getParameter<std::vector<double> >("thePUEfficiency_FPix_Inner"));
509  thePUEfficiency.push_back(conf.getParameter<std::vector<double> >("thePUEfficiency_FPix_Outer"));
510  if ( ((thePUEfficiency[3].empty()) || (thePUEfficiency[4].empty())) && (NumberOfEndcapDisks==2) )
511  throw cms::Exception("Configuration") << "At least one (FPix) PU efficiency number is needed in efficiency config!";
512  pu_scale.resize(thePUEfficiency.size());
513  }
514  }
515  else LogInfo ("PixelDigitizer ") <<"The PixelDigitizer inefficiency configuration is read from the database.\n";
516  }
517  // the first "NumberOfBarrelLayers" settings [0],[1], ... , [NumberOfBarrelLayers-1] are for the barrel pixels
518  // the next "NumberOfEndcapDisks" settings [NumberOfBarrelLayers],[NumberOfBarrelLayers+1], ... [NumberOfEndcapDisks+NumberOfBarrelLayers-1]
519 }
size
Write out results.
T getParameter(std::string const &) const
bool exists(std::string const &parameterName) const
checks if a parameter exists
std::vector< std::vector< double > > thePUEfficiency

Member Function Documentation

void SiPixelDigitizerAlgorithm::PixelEfficiencies::init_from_db ( const edm::ESHandle< TrackerGeometry > &  geom,
const edm::ESHandle< SiPixelDynamicInefficiency > &  SiPixelDynamicInefficiency 
)

Definition at line 530 of file SiPixelDigitizerAlgorithm.cc.

References PixelTopology::colsperroc(), TrackerGeometry::detUnits(), Exception, SiPixelDynamicInefficiency::getChipGeomFactors(), SiPixelDynamicInefficiency::getColGeomFactors(), SiPixelDynamicInefficiency::getDetIdmasks(), SiPixelDynamicInefficiency::getPixelGeomFactors(), SiPixelDynamicInefficiency::getPUFactors(), SiPixelDynamicInefficiency::gettheInstLumiScaleFactor(), mps_fire::i, TrackerGeometry::idToDet(), TrackerGeometry::isThere(), patRefSel_triggerMatching_cfi::matches, SiStripPI::max, min(), GeomDetEnumerators::P1PXB, GeomDetEnumerators::P1PXEC, PixelSubdetector::PixelBarrel, DetId::rawId(), PixelTopology::rowsperroc(), edm::shift, PixelGeomDetUnit::specificTopology(), DetId::subdetId(), and ecaldqm::topology().

Referenced by SiPixelDigitizerAlgorithm::init_DynIneffDB().

530  {
531 
532  theInstLumiScaleFactor = SiPixelDynamicInefficiency->gettheInstLumiScaleFactor();
533  const std::map<uint32_t, double>& PixelGeomFactorsDBIn = SiPixelDynamicInefficiency->getPixelGeomFactors();
534  const std::map<uint32_t, double>& ColGeomFactorsDB = SiPixelDynamicInefficiency->getColGeomFactors();
535  const std::map<uint32_t, double>& ChipGeomFactorsDB = SiPixelDynamicInefficiency->getChipGeomFactors();
536  const std::map<uint32_t, std::vector<double> >& PUFactors = SiPixelDynamicInefficiency->getPUFactors();
537  std::vector<uint32_t > DetIdmasks = SiPixelDynamicInefficiency->getDetIdmasks();
538 
539  // Loop on all modules, initialize map for easy access
540  for( const auto& it_module : geom->detUnits()) {
541  if( dynamic_cast<PixelGeomDetUnit const*>(it_module)==nullptr) continue;
542  const DetId detid = it_module->geographicalId();
543  uint32_t rawid = detid.rawId();
544  PixelGeomFactors[rawid] = 1;
545  ColGeomFactors[rawid] = 1;
546  ChipGeomFactors[rawid] = 1;
547  PixelGeomFactorsROCStdPixels[rawid] = std::vector<double>(16,1);
548  PixelGeomFactorsROCBigPixels[rawid] = std::vector<double>(16,1);
549  }
550 
551  // ROC level inefficiency for phase 1 (disentangle scale factors for big and std size pixels)
552  std::map<uint32_t, double> PixelGeomFactorsDB;
553 
556  for (auto db_factor : PixelGeomFactorsDBIn){
557  int shift = DetId(db_factor.first).subdetId() ==
559  unsigned int rocMask = rocIdMaskBits << shift;
560  unsigned int rocId = (((db_factor.first) & rocMask) >> shift);
561  if (rocId != 0) {
562  rocId--;
563  unsigned int rawid = db_factor.first & (~rocMask);
564  const PixelGeomDetUnit * theGeomDet = dynamic_cast<const PixelGeomDetUnit*> (geom->idToDet(rawid));
565  PixelTopology const * topology = &(theGeomDet->specificTopology());
566  const int nPixelsInROC = topology->rowsperroc()*topology->colsperroc();
567  const int nBigPixelsInROC = 2*topology->rowsperroc()+topology->colsperroc()-2;
568  double factor = db_factor.second;
569  double badFraction = 1 - factor;
570  double bigPixelFraction = static_cast<double> (nBigPixelsInROC)/nPixelsInROC;
571  double stdPixelFraction = 1. - bigPixelFraction;
572 
573  double badFractionBig = std::min(bigPixelFraction, badFraction);
574  double badFractionStd = std::max(0., badFraction - badFractionBig);
575  double badFractionBigReNormalized = badFractionBig/bigPixelFraction;
576  double badFractionStdReNormalized = badFractionStd/stdPixelFraction;
577  PixelGeomFactorsROCStdPixels[rawid][rocId] *= (1. - badFractionStdReNormalized);
578  PixelGeomFactorsROCBigPixels[rawid][rocId] *= (1. - badFractionBigReNormalized);
579  }
580  else{
581  PixelGeomFactorsDB[db_factor.first] = db_factor.second;
582  }
583  }
584  } // is Phase 1 geometry
585  else{
586  PixelGeomFactorsDB = PixelGeomFactorsDBIn;
587  }
588 
589  // Loop on all modules, store module level geometrical scale factors
590  for( const auto& it_module : geom->detUnits()) {
591  if( dynamic_cast<PixelGeomDetUnit const*>(it_module)==nullptr) continue;
592  const DetId detid = it_module->geographicalId();
593  uint32_t rawid = detid.rawId();
594  for (auto db_factor : PixelGeomFactorsDB) if (matches(detid, DetId(db_factor.first), DetIdmasks)) PixelGeomFactors[rawid] *= db_factor.second;
595  for (auto db_factor : ColGeomFactorsDB) if (matches(detid, DetId(db_factor.first), DetIdmasks)) ColGeomFactors[rawid] *= db_factor.second;
596  for (auto db_factor : ChipGeomFactorsDB) if (matches(detid, DetId(db_factor.first), DetIdmasks)) ChipGeomFactors[rawid] *= db_factor.second;
597  }
598 
599  // piluep scale factors are calculated once per event
600  // therefore vector index is stored in a map for each module that matches to a db_id
601  size_t i=0;
602  for (auto factor : PUFactors) {
603  const DetId db_id = DetId(factor.first);
604  for( const auto& it_module : geom->detUnits()) {
605  if( dynamic_cast<PixelGeomDetUnit const*>(it_module)==nullptr) continue;
606  const DetId detid = it_module->geographicalId();
607  if (!matches(detid, db_id, DetIdmasks)) continue;
608  if (iPU.count(detid.rawId())) {
609  throw cms::Exception("Database")<<"Multiple db_ids match to same module in SiPixelDynamicInefficiency DB Object";
610  } else {
611  iPU[detid.rawId()] = i;
612  }
613  }
614  thePUEfficiency.push_back(factor.second);
615  ++i;
616  }
617  pu_scale.resize(thePUEfficiency.size());
618 }
virtual int rowsperroc() const =0
CaloTopology const * topology(0)
const DetContainer & detUnits() const override
Returm a vector of all GeomDet.
constexpr uint32_t rawId() const
get the raw id
Definition: DetId.h:47
std::vector< std::vector< double > > thePUEfficiency
bool isThere(GeomDetEnumerators::SubDetector subdet) const
const std::vector< uint32_t > getDetIdmasks() const
virtual int colsperroc() const =0
std::map< uint32_t, std::vector< double > > PixelGeomFactorsROCBigPixels
constexpr int subdetId() const
get the contents of the subdetector field (not cast into any detector&#39;s numbering enum) ...
Definition: DetId.h:41
T min(T a, T b)
Definition: MathUtil.h:58
Definition: DetId.h:18
std::map< uint32_t, std::vector< double > > PixelGeomFactorsROCStdPixels
virtual const PixelTopology & specificTopology() const
Returns a reference to the pixel proxy topology.
bool matches(const DetId &, const DetId &, const std::vector< uint32_t > &)
const TrackerGeomDet * idToDet(DetId) const override
static unsigned int const shift
const std::map< unsigned int, std::vector< double > > & getPUFactors() const
const std::map< unsigned int, double > & getColGeomFactors() const
const std::map< unsigned int, double > & getPixelGeomFactors() const
const std::map< unsigned int, double > & getChipGeomFactors() const
bool SiPixelDigitizerAlgorithm::PixelEfficiencies::matches ( const DetId detid,
const DetId db_id,
const std::vector< uint32_t > &  DetIdmasks 
)

Definition at line 620 of file SiPixelDigitizerAlgorithm.cc.

References DetId::det(), mps_fire::i, DetId::rawId(), and DetId::subdetId().

620  {
621  if (detid.subdetId() != db_id.subdetId()) return false;
622  for (size_t i=0; i<DetIdmasks.size(); ++i) {
623  DetId maskid = DetId(DetIdmasks.at(i));
624  if (maskid.subdetId() != db_id.subdetId()) continue;
625  if ((detid.rawId()&maskid.rawId()) != (db_id.rawId()&maskid.rawId()) &&
626  (db_id.rawId()&maskid.rawId()) != DetId(db_id.det(), db_id.subdetId()).rawId()) return false;
627  }
628  return true;
629 }
constexpr uint32_t rawId() const
get the raw id
Definition: DetId.h:47
constexpr int subdetId() const
get the contents of the subdetector field (not cast into any detector&#39;s numbering enum) ...
Definition: DetId.h:41
Definition: DetId.h:18
constexpr Detector det() const
get the detector field from this detid
Definition: DetId.h:39

Member Data Documentation

std::map<uint32_t, double> SiPixelDigitizerAlgorithm::PixelEfficiencies::ChipGeomFactors
std::map<uint32_t, double> SiPixelDigitizerAlgorithm::PixelEfficiencies::ColGeomFactors
unsigned int SiPixelDigitizerAlgorithm::PixelEfficiencies::FPixIndex
bool SiPixelDigitizerAlgorithm::PixelEfficiencies::FromConfig
std::map<uint32_t, size_t > SiPixelDigitizerAlgorithm::PixelEfficiencies::iPU
std::map<uint32_t, double> SiPixelDigitizerAlgorithm::PixelEfficiencies::PixelGeomFactors
std::map<uint32_t, std::vector<double> > SiPixelDigitizerAlgorithm::PixelEfficiencies::PixelGeomFactorsROCBigPixels
std::map<uint32_t, std::vector<double> > SiPixelDigitizerAlgorithm::PixelEfficiencies::PixelGeomFactorsROCStdPixels
std::vector<double> SiPixelDigitizerAlgorithm::PixelEfficiencies::pu_scale
const int SiPixelDigitizerAlgorithm::PixelEfficiencies::rocIdMaskBits = 0x1F
static

Definition at line 269 of file SiPixelDigitizerAlgorithm.h.

double SiPixelDigitizerAlgorithm::PixelEfficiencies::theInnerEfficiency_FPix[20]
double SiPixelDigitizerAlgorithm::PixelEfficiencies::theInstLumiScaleFactor
std::vector<double> SiPixelDigitizerAlgorithm::PixelEfficiencies::theLadderEfficiency_BPix[20]
std::vector<double> SiPixelDigitizerAlgorithm::PixelEfficiencies::theModuleEfficiency_BPix[20]
double SiPixelDigitizerAlgorithm::PixelEfficiencies::theOuterEfficiency_FPix[20]
double SiPixelDigitizerAlgorithm::PixelEfficiencies::thePixelChipEfficiency[20]
double SiPixelDigitizerAlgorithm::PixelEfficiencies::thePixelColEfficiency[20]
double SiPixelDigitizerAlgorithm::PixelEfficiencies::thePixelEfficiency[20]
std::vector<std::vector<double> > SiPixelDigitizerAlgorithm::PixelEfficiencies::thePUEfficiency