edm::DataMixingSiPixelMCDigiWorker::PixelEfficiencies Struct Reference

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 83 of file DataMixingSiPixelMCDigiWorker.h.

Constructor & Destructor Documentation

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

Definition at line 92 of file DataMixingSiPixelMCDigiWorker.cc.

References Exception, edm::ParameterSet::exists(), FPixIndex, FromConfig, edm::ParameterSet::getParameter(), i, j, edm::DataMixingSiPixelMCDigiWorker::NumberOfBarrelLayers, edm::DataMixingSiPixelMCDigiWorker::NumberOfEndcapDisks, pu_scale, findQualityFiles::size, theInnerEfficiency_FPix, theInstLumiScaleFactor, theLadderEfficiency_BPix, theModuleEfficiency_BPix, theOuterEfficiency_FPix, thePixelChipEfficiency, thePixelColEfficiency, thePixelEfficiency, and thePUEfficiency.

                                                                                                                                                                             {
  // 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].size()==0) || (thePUEfficiency[1].size()==0) || 
              (thePUEfficiency[2].size()==0)) && (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].size()==0) || (thePUEfficiency[4].size()==0)) && (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 edm::DataMixingSiPixelMCDigiWorker::PixelEfficiencies::init_from_db	(	const edm::ESHandle< TrackerGeometry > &	geom,
		const edm::ESHandle< SiPixelDynamicInefficiency > &	SiPixelDynamicInefficiency
	)

Definition at line 225 of file DataMixingSiPixelMCDigiWorker.cc.

References cond::rpcobgas::detid, Exception, i, matches, and DetId::rawId().

Referenced by edm::DataMixingSiPixelMCDigiWorker::init_DynIneffDB().

                                                                                                                                                                                     {

  theInstLumiScaleFactor = SiPixelDynamicInefficiency->gettheInstLumiScaleFactor();
  std::map<uint32_t, double> PixelGeomFactorsDB = SiPixelDynamicInefficiency->getPixelGeomFactors();
  std::map<uint32_t, double> ColGeomFactorsDB = SiPixelDynamicInefficiency->getColGeomFactors();
  std::map<uint32_t, double> ChipGeomFactorsDB = SiPixelDynamicInefficiency->getChipGeomFactors();
  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(TrackerGeometry::DetUnitContainer::const_iterator it_module = geom->detUnits().begin(); it_module != geom->detUnits().end(); it_module++) {
    if( dynamic_cast<PixelGeomDetUnit const*>((*it_module))==0) 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(TrackerGeometry::DetUnitContainer::const_iterator it_module = geom->detUnits().begin(); it_module != geom->detUnits().end(); it_module++) {
      if( dynamic_cast<PixelGeomDetUnit const*>((*it_module))==0) 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 edm::DataMixingSiPixelMCDigiWorker::PixelEfficiencies::matches	(	const DetId &	detid,
		const DetId &	db_id,
		const std::vector< uint32_t > &	DetIdmasks
	)

Definition at line 268 of file DataMixingSiPixelMCDigiWorker.cc.

References DetId::det(), 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;
}

Member Data Documentation

std::map<uint32_t, double> edm::DataMixingSiPixelMCDigiWorker::PixelEfficiencies::ChipGeomFactors

Definition at line 102 of file DataMixingSiPixelMCDigiWorker.h.

Referenced by edm::DataMixingSiPixelMCDigiWorker::putSiPixel().

std::map<uint32_t, double> edm::DataMixingSiPixelMCDigiWorker::PixelEfficiencies::ColGeomFactors

Definition at line 101 of file DataMixingSiPixelMCDigiWorker.h.

Referenced by edm::DataMixingSiPixelMCDigiWorker::putSiPixel().

unsigned int edm::DataMixingSiPixelMCDigiWorker::PixelEfficiencies::FPixIndex

Definition at line 98 of file DataMixingSiPixelMCDigiWorker.h.

Referenced by PixelEfficiencies(), and edm::DataMixingSiPixelMCDigiWorker::putSiPixel().

bool edm::DataMixingSiPixelMCDigiWorker::PixelEfficiencies::FromConfig

Definition at line 85 of file DataMixingSiPixelMCDigiWorker.h.

Referenced by edm::DataMixingSiPixelMCDigiWorker::init_DynIneffDB(), PixelEfficiencies(), and edm::DataMixingSiPixelMCDigiWorker::putSiPixel().

std::map<uint32_t, size_t > edm::DataMixingSiPixelMCDigiWorker::PixelEfficiencies::iPU

Definition at line 103 of file DataMixingSiPixelMCDigiWorker.h.

Referenced by edm::DataMixingSiPixelMCDigiWorker::putSiPixel().

std::map<uint32_t, double> edm::DataMixingSiPixelMCDigiWorker::PixelEfficiencies::PixelGeomFactors

Definition at line 100 of file DataMixingSiPixelMCDigiWorker.h.

Referenced by edm::DataMixingSiPixelMCDigiWorker::putSiPixel().

std::vector<double> edm::DataMixingSiPixelMCDigiWorker::PixelEfficiencies::pu_scale

Definition at line 88 of file DataMixingSiPixelMCDigiWorker.h.

Referenced by PixelEfficiencies(), edm::DataMixingSiPixelMCDigiWorker::putSiPixel(), and edm::DataMixingSiPixelMCDigiWorker::setPileupInfo().

double edm::DataMixingSiPixelMCDigiWorker::PixelEfficiencies::theInnerEfficiency_FPix[20]

Definition at line 96 of file DataMixingSiPixelMCDigiWorker.h.

Referenced by PixelEfficiencies(), and edm::DataMixingSiPixelMCDigiWorker::putSiPixel().

double edm::DataMixingSiPixelMCDigiWorker::PixelEfficiencies::theInstLumiScaleFactor

Definition at line 87 of file DataMixingSiPixelMCDigiWorker.h.

Referenced by PixelEfficiencies(), and edm::DataMixingSiPixelMCDigiWorker::setPileupInfo().

std::vector<double> edm::DataMixingSiPixelMCDigiWorker::PixelEfficiencies::theLadderEfficiency_BPix[20]

Definition at line 93 of file DataMixingSiPixelMCDigiWorker.h.

Referenced by PixelEfficiencies(), and edm::DataMixingSiPixelMCDigiWorker::putSiPixel().

std::vector<double> edm::DataMixingSiPixelMCDigiWorker::PixelEfficiencies::theModuleEfficiency_BPix[20]

Definition at line 94 of file DataMixingSiPixelMCDigiWorker.h.

Referenced by PixelEfficiencies(), and edm::DataMixingSiPixelMCDigiWorker::putSiPixel().

double edm::DataMixingSiPixelMCDigiWorker::PixelEfficiencies::theOuterEfficiency_FPix[20]

Definition at line 97 of file DataMixingSiPixelMCDigiWorker.h.

Referenced by PixelEfficiencies(), and edm::DataMixingSiPixelMCDigiWorker::putSiPixel().

double edm::DataMixingSiPixelMCDigiWorker::PixelEfficiencies::thePixelChipEfficiency[20]

Definition at line 92 of file DataMixingSiPixelMCDigiWorker.h.

Referenced by PixelEfficiencies(), and edm::DataMixingSiPixelMCDigiWorker::putSiPixel().

double edm::DataMixingSiPixelMCDigiWorker::PixelEfficiencies::thePixelColEfficiency[20]

Definition at line 91 of file DataMixingSiPixelMCDigiWorker.h.

Referenced by PixelEfficiencies(), and edm::DataMixingSiPixelMCDigiWorker::putSiPixel().

double edm::DataMixingSiPixelMCDigiWorker::PixelEfficiencies::thePixelEfficiency[20]

Definition at line 90 of file DataMixingSiPixelMCDigiWorker.h.

Referenced by PixelEfficiencies(), and edm::DataMixingSiPixelMCDigiWorker::putSiPixel().

std::vector<std::vector<double> > edm::DataMixingSiPixelMCDigiWorker::PixelEfficiencies::thePUEfficiency

Definition at line 89 of file DataMixingSiPixelMCDigiWorker.h.

Referenced by PixelEfficiencies(), and edm::DataMixingSiPixelMCDigiWorker::setPileupInfo().