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 399 of file SiPixelDigitizerAlgorithm.cc.

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

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

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

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

621  {
622  if (detid.subdetId() != db_id.subdetId()) return false;
623  for (size_t i=0; i<DetIdmasks.size(); ++i) {
624  DetId maskid = DetId(DetIdmasks.at(i));
625  if (maskid.subdetId() != db_id.subdetId()) continue;
626  if ((detid.rawId()&maskid.rawId()) != (db_id.rawId()&maskid.rawId()) &&
627  (db_id.rawId()&maskid.rawId()) != DetId(db_id.det(), db_id.subdetId()).rawId()) return false;
628  }
629  return true;
630 }
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