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::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

Detailed Description

Internal use only.

Definition at line 230 of file SiPixelDigitizerAlgorithm.h.

Constructor & Destructor Documentation

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

Definition at line 366 of file SiPixelDigitizerAlgorithm.cc.

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

                                                                                                                                                                        {
   // pixel inefficiency
   // Don't use Hard coded values, read inefficiencies in from DB/python config or don't use any
   int NumberOfTotLayers = NumberOfBarrelLayers + NumberOfEndcapDisks;
   FPixIndex=NumberOfBarrelLayers;
   if (AddPixelInefficiency){
     FromConfig = 
       conf.exists("thePixelColEfficiency_BPix1") && conf.exists("thePixelColEfficiency_BPix2") && conf.exists("thePixelColEfficiency_BPix3") &&
       conf.exists("thePixelColEfficiency_FPix1") && conf.exists("thePixelColEfficiency_FPix2") &&
       conf.exists("thePixelEfficiency_BPix1") && conf.exists("thePixelEfficiency_BPix2") && conf.exists("thePixelEfficiency_BPix3") &&
       conf.exists("thePixelEfficiency_FPix1") && conf.exists("thePixelEfficiency_FPix2") &&
       conf.exists("thePixelChipEfficiency_BPix1") && conf.exists("thePixelChipEfficiency_BPix2") && conf.exists("thePixelChipEfficiency_BPix3") &&
       conf.exists("thePixelChipEfficiency_FPix1") && conf.exists("thePixelChipEfficiency_FPix2");
     if (NumberOfBarrelLayers==3) FromConfig = FromConfig && conf.exists("theLadderEfficiency_BPix1") && conf.exists("theLadderEfficiency_BPix2") && conf.exists("theLadderEfficiency_BPix3") &&
       conf.exists("theModuleEfficiency_BPix1") && conf.exists("theModuleEfficiency_BPix2") && conf.exists("theModuleEfficiency_BPix3") &&
       conf.exists("thePUEfficiency_BPix1") && conf.exists("thePUEfficiency_BPix2") && conf.exists("thePUEfficiency_BPix3") &&
       conf.exists("theInnerEfficiency_FPix1") && conf.exists("theInnerEfficiency_FPix2") &&
       conf.exists("theOuterEfficiency_FPix1") && conf.exists("theOuterEfficiency_FPix2") &&
       conf.exists("thePUEfficiency_FPix_Inner") && conf.exists("thePUEfficiency_FPix_Outer") &&
       conf.exists("theInstLumiScaleFactor");
     if (NumberOfBarrelLayers>=4) FromConfig = FromConfig && conf.exists("thePixelColEfficiency_BPix4") &&
       conf.exists("thePixelEfficiency_BPix4") && conf.exists("thePixelChipEfficiency_BPix4");
     if (NumberOfEndcapDisks>=3) FromConfig = FromConfig && conf.exists("thePixelColEfficiency_FPix4") &&
       conf.exists("thePixelEfficiency_FPix3") && conf.exists("thePixelChipEfficiency_FPix3");
     if (FromConfig) {
       LogInfo ("PixelDigitizer ") <<"The PixelDigitizer inefficiency configuration is read from the config file.\n";
       theInstLumiScaleFactor = conf.getParameter<double>("theInstLumiScaleFactor");
       int i=0;
       thePixelColEfficiency[i++] = conf.getParameter<double>("thePixelColEfficiency_BPix1");
       thePixelColEfficiency[i++] = conf.getParameter<double>("thePixelColEfficiency_BPix2");
       thePixelColEfficiency[i++] = conf.getParameter<double>("thePixelColEfficiency_BPix3");
       if (NumberOfBarrelLayers>=4){thePixelColEfficiency[i++] = conf.getParameter<double>("thePixelColEfficiency_BPix4");}
       //
       i=0;
       thePixelEfficiency[i++] = conf.getParameter<double>("thePixelEfficiency_BPix1");
       thePixelEfficiency[i++] = conf.getParameter<double>("thePixelEfficiency_BPix2");
       thePixelEfficiency[i++] = conf.getParameter<double>("thePixelEfficiency_BPix3");
       if (NumberOfBarrelLayers>=4){thePixelEfficiency[i++] = conf.getParameter<double>("thePixelEfficiency_BPix4");}
       //
       i=0;
       thePixelChipEfficiency[i++] = conf.getParameter<double>("thePixelChipEfficiency_BPix1");
       thePixelChipEfficiency[i++] = conf.getParameter<double>("thePixelChipEfficiency_BPix2");
       thePixelChipEfficiency[i++] = conf.getParameter<double>("thePixelChipEfficiency_BPix3");
       if (NumberOfBarrelLayers>=4){thePixelChipEfficiency[i++] = conf.getParameter<double>("thePixelChipEfficiency_BPix4");}
       //
       if (NumberOfBarrelLayers==3){
         i=0;
         theLadderEfficiency_BPix[i++] = conf.getParameter<std::vector<double> >("theLadderEfficiency_BPix1");
         theLadderEfficiency_BPix[i++] = conf.getParameter<std::vector<double> >("theLadderEfficiency_BPix2");
         theLadderEfficiency_BPix[i++] = conf.getParameter<std::vector<double> >("theLadderEfficiency_BPix3");
         if ( ((theLadderEfficiency_BPix[0].size()!=20) || (theLadderEfficiency_BPix[1].size()!=32) ||
               (theLadderEfficiency_BPix[2].size()!=44)) && (NumberOfBarrelLayers==3) )  
           throw cms::Exception("Configuration") << "Wrong ladder number in efficiency config!";
         //           
         i=0;
         theModuleEfficiency_BPix[i++] = conf.getParameter<std::vector<double> >("theModuleEfficiency_BPix1");
         theModuleEfficiency_BPix[i++] = conf.getParameter<std::vector<double> >("theModuleEfficiency_BPix2");
         theModuleEfficiency_BPix[i++] = conf.getParameter<std::vector<double> >("theModuleEfficiency_BPix3");
         if ( ((theModuleEfficiency_BPix[0].size()!=4) || (theModuleEfficiency_BPix[1].size()!=4) ||
               (theModuleEfficiency_BPix[2].size()!=4)) && (NumberOfBarrelLayers==3) )  
           throw cms::Exception("Configuration") << "Wrong module number in efficiency config!";
         //
         thePUEfficiency.push_back(conf.getParameter<std::vector<double> >("thePUEfficiency_BPix1"));
         thePUEfficiency.push_back(conf.getParameter<std::vector<double> >("thePUEfficiency_BPix2"));
         thePUEfficiency.push_back(conf.getParameter<std::vector<double> >("thePUEfficiency_BPix3"));                        
         if ( ((thePUEfficiency[0].empty()) || (thePUEfficiency[1].empty()) || 
               (thePUEfficiency[2].empty())) && (NumberOfBarrelLayers==3) )
           throw cms::Exception("Configuration") << "At least one PU efficiency (BPix) number is needed in efficiency config!";
       }
       // The next is needed for Phase2 Tracker studies
       if (NumberOfBarrelLayers>=5){
         if (NumberOfTotLayers>20){throw cms::Exception("Configuration") <<"SiPixelDigitizer was given more layers than it can handle";}
         // For Phase2 tracker layers just set the outermost BPix inefficiency to 99.9% THESE VALUES ARE HARDCODED ALSO ELSEWHERE IN THIS FILE
         for (int j=5 ; j<=NumberOfBarrelLayers ; j++){
           thePixelColEfficiency[j-1]=0.999;
           thePixelEfficiency[j-1]=0.999;
           thePixelChipEfficiency[j-1]=0.999;
         }
       }
       //
       i=FPixIndex;
       thePixelColEfficiency[i++]   = conf.getParameter<double>("thePixelColEfficiency_FPix1");
       thePixelColEfficiency[i++]   = conf.getParameter<double>("thePixelColEfficiency_FPix2");
       if (NumberOfEndcapDisks>=3){thePixelColEfficiency[i++]   = conf.getParameter<double>("thePixelColEfficiency_FPix3");}
       i=FPixIndex;
       thePixelEfficiency[i++]      = conf.getParameter<double>("thePixelEfficiency_FPix1");
       thePixelEfficiency[i++]      = conf.getParameter<double>("thePixelEfficiency_FPix2");
       if (NumberOfEndcapDisks>=3){thePixelEfficiency[i++]      = conf.getParameter<double>("thePixelEfficiency_FPix3");}
       i=FPixIndex;
       thePixelChipEfficiency[i++]  = conf.getParameter<double>("thePixelChipEfficiency_FPix1");
       thePixelChipEfficiency[i++]  = conf.getParameter<double>("thePixelChipEfficiency_FPix2");
       if (NumberOfEndcapDisks>=3){thePixelChipEfficiency[i++]  = conf.getParameter<double>("thePixelChipEfficiency_FPix3");}
       // The next is needed for Phase2 Tracker studies
       if (NumberOfEndcapDisks>=4){
         if (NumberOfTotLayers>20){throw cms::Exception("Configuration") <<"SiPixelDigitizer was given more layers than it can handle";}
         // For Phase2 tracker layers just set the extra FPix disk inefficiency to 99.9% THESE VALUES ARE HARDCODED ALSO ELSEWHERE IN THIS FILE
         for (int j=4+FPixIndex ; j<=NumberOfEndcapDisks+NumberOfBarrelLayers ; j++){
           thePixelColEfficiency[j-1]=0.999;
           thePixelEfficiency[j-1]=0.999;
           thePixelChipEfficiency[j-1]=0.999;
         }
       }
       //FPix Dynamic Inefficiency
       if (NumberOfBarrelLayers==3){
         i=FPixIndex;
         theInnerEfficiency_FPix[i++] = conf.getParameter<double>("theInnerEfficiency_FPix1");
         theInnerEfficiency_FPix[i++] = conf.getParameter<double>("theInnerEfficiency_FPix2");
         i=FPixIndex;
         theOuterEfficiency_FPix[i++] = conf.getParameter<double>("theOuterEfficiency_FPix1");
         theOuterEfficiency_FPix[i++] = conf.getParameter<double>("theOuterEfficiency_FPix2");
         thePUEfficiency.push_back(conf.getParameter<std::vector<double> >("thePUEfficiency_FPix_Inner"));
         thePUEfficiency.push_back(conf.getParameter<std::vector<double> >("thePUEfficiency_FPix_Outer"));
         if ( ((thePUEfficiency[3].empty()) || (thePUEfficiency[4].empty())) && (NumberOfEndcapDisks==2) )
           throw cms::Exception("Configuration") << "At least one (FPix) PU efficiency number is needed in efficiency config!";
         pu_scale.resize(thePUEfficiency.size());
       }
     }
     else LogInfo ("PixelDigitizer ") <<"The PixelDigitizer inefficiency configuration is read from the database.\n";
   }
   // the first "NumberOfBarrelLayers" settings [0],[1], ... , [NumberOfBarrelLayers-1] are for the barrel pixels
   // the next  "NumberOfEndcapDisks"  settings [NumberOfBarrelLayers],[NumberOfBarrelLayers+1], ... [NumberOfEndcapDisks+NumberOfBarrelLayers-1]
 }

Member Function Documentation

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

Definition at line 498 of file SiPixelDigitizerAlgorithm.cc.

References TrackerGeometry::detUnits(), Exception, SiPixelDynamicInefficiency::getChipGeomFactors(), SiPixelDynamicInefficiency::getColGeomFactors(), SiPixelDynamicInefficiency::getDetIdmasks(), SiPixelDynamicInefficiency::getPixelGeomFactors(), SiPixelDynamicInefficiency::getPUFactors(), SiPixelDynamicInefficiency::gettheInstLumiScaleFactor(), mps_fire::i, patRefSel_triggerMatching_cfi::matches, and DetId::rawId().

Referenced by SiPixelDigitizerAlgorithm::init_DynIneffDB().

                                                                                                                                                                                  {
 
   theInstLumiScaleFactor = SiPixelDynamicInefficiency->gettheInstLumiScaleFactor();
   const std::map<uint32_t, double>& PixelGeomFactorsDB = SiPixelDynamicInefficiency->getPixelGeomFactors();
   const std::map<uint32_t, double>& ColGeomFactorsDB = SiPixelDynamicInefficiency->getColGeomFactors();
   const std::map<uint32_t, double>& ChipGeomFactorsDB = SiPixelDynamicInefficiency->getChipGeomFactors();
   const std::map<uint32_t, std::vector<double> >& PUFactors = SiPixelDynamicInefficiency->getPUFactors();
   std::vector<uint32_t > DetIdmasks = SiPixelDynamicInefficiency->getDetIdmasks();
   
   // Loop on all modules, calculate geometrical scale factors and store in map for easy access
   for( const auto& it_module : geom->detUnits()) {
     if( dynamic_cast<PixelGeomDetUnit const*>(it_module)==nullptr) continue;
     const DetId detid = it_module->geographicalId();
     uint32_t rawid = detid.rawId();
     PixelGeomFactors[rawid] = 1;
     ColGeomFactors[rawid] = 1;
     ChipGeomFactors[rawid] = 1;
     for (auto db_factor : PixelGeomFactorsDB) if (matches(detid, DetId(db_factor.first), DetIdmasks)) PixelGeomFactors[rawid] *= db_factor.second;
     for (auto db_factor : ColGeomFactorsDB) if (matches(detid, DetId(db_factor.first), DetIdmasks)) ColGeomFactors[rawid] *= db_factor.second;
     for (auto db_factor : ChipGeomFactorsDB) if (matches(detid, DetId(db_factor.first), DetIdmasks)) ChipGeomFactors[rawid] *= db_factor.second;
   }
 
   // piluep scale factors are calculated once per event
   // therefore vector index is stored in a map for each module that matches to a db_id
   size_t i=0;
   for (auto factor : PUFactors) {
     const DetId db_id = DetId(factor.first);
     for( const auto& it_module : geom->detUnits()) {
       if( dynamic_cast<PixelGeomDetUnit const*>(it_module)==nullptr) continue;
       const DetId detid = it_module->geographicalId();
       if (!matches(detid, db_id, DetIdmasks)) continue;
       if (iPU.count(detid.rawId())) {
     throw cms::Exception("Database")<<"Multiple db_ids match to same module in SiPixelDynamicInefficiency DB Object";
       } else {
     iPU[detid.rawId()] = i;
       }
     }
     thePUEfficiency.push_back(factor.second);
     ++i;
   }
   pu_scale.resize(thePUEfficiency.size());
 }

bool SiPixelDigitizerAlgorithm::PixelEfficiencies::matches	(	const DetId &	detid,
		const DetId &	db_id,
		const std::vector< uint32_t > &	DetIdmasks
	)

Definition at line 541 of file SiPixelDigitizerAlgorithm.cc.

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

                                                                                                                                        {
   if (detid.subdetId() != db_id.subdetId()) return false;
   for (size_t i=0; i<DetIdmasks.size(); ++i) {
     DetId maskid = DetId(DetIdmasks.at(i));
     if (maskid.subdetId() != db_id.subdetId()) continue;
     if ((detid.rawId()&maskid.rawId()) != (db_id.rawId()&maskid.rawId()) && 
     (db_id.rawId()&maskid.rawId()) != DetId(db_id.det(), db_id.subdetId()).rawId()) return false;
   }
   return true;
 }