Inheritance diagram for SiStripHitEffFromCalibTree:

Public Member Functions
	SiStripHitEffFromCalibTree (const edm::ParameterSet &)

	~SiStripHitEffFromCalibTree () override=default

Public Member Functions inherited from ConditionDBWriter< SiStripBadStrip >
	ConditionDBWriter (const edm::ParameterSet &iConfig)

	~ConditionDBWriter () override

Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::WatchRuns, edm::one::WatchLuminosityBlocks, edm::one::SharedResources >
	EDAnalyzer ()=default

	EDAnalyzer (const EDAnalyzer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

const EDAnalyzer &	operator= (const EDAnalyzer &)=delete

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const final

Public Member Functions inherited from edm::one::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDAnalyzerBase () override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

std::vector< ESProxyIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const *> const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	algoAnalyze (const edm::Event &e, const edm::EventSetup &c) override

void	algoBeginJob (const edm::EventSetup &) override

void	algoEndJob () override

void	computeEff (vector< TH1F > &vhfound, vector< TH1F > &vhtotal, string name)

TString	getLayerSideName (Long_t k)

std::unique_ptr< SiStripBadStrip >	getNewObject () override

void	makeHotColdMaps ()

void	makeSQLite ()

void	makeSummary ()

void	makeSummaryVsBx ()

void	makeSummaryVsCM ()

void	makeSummaryVsLumi ()

void	makeTKMap (bool autoTagging)

void	setBadComponents (int i, int component, SiStripQuality::BadComponent &BC, std::stringstream ssV[4][19], int NBadComponent[4][19][4])

void	totalStatistics ()

Private Attributes
bool	_autoIneffModTagging

string	_badModulesFile

unsigned int	_bunchx

unsigned int	_clusterMatchingMethod

float	_clusterTrajDist

SiStripDetInfo	_detInfo

float	_effPlotMin

float	_ResXSig

bool	_showEndcapSides

bool	_showOnlyGoodModules

bool	_showRings

bool	_showTOB6TEC9

unsigned int	_spaceBetweenTrains

float	_stripsApvEdge

TString	_title

edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	_tkGeomToken

float	_tkMapMin

edm::ESGetToken< TrackerTopology, TrackerTopologyRcd >	_tTopoToken

bool	_useCM

bool	_useOnlyHighPurityTracks

int	alllayerfound [35]

int	alllayertotal [35]

map< unsigned int, double >	BadModules

TH1F *	bxHisto

TTree *	CalibTree

vector< string >	CalibTreeFilenames

unsigned int	doSummary

map< pair< unsigned int, unsigned int >, array< double, 3 > >	eventInfos

edm::FileInPath	FileInPath_

edm::Service< TFileService >	fs

int	goodlayerfound [35]

int	goodlayertotal [35]

vector< hit >	hits [23]

vector< TH2F * >	HotColdMaps

TH1F *	instLumiHisto

long	layerfound [23]

map< unsigned int, vector< int > >	layerfound_perBx

vector< TH1F * >	layerfound_vsCM

vector< TH1F * >	layerfound_vsLumi

vector< TH1F * >	layerfound_vsPU

long	layertotal [23]

map< unsigned int, vector< int > >	layertotal_perBx

vector< TH1F * >	layertotal_vsCM

vector< TH1F * >	layertotal_vsLumi

vector< TH1F * >	layertotal_vsPU

map< unsigned int, pair< unsigned int, unsigned int > >	modCounter [23]

unsigned int	nModsMin

unsigned int	nTEClayers

TH1F *	PUHisto

SiStripQuality *	quality_

float	threshold

TrackerMap *	tkmap

TrackerMap *	tkmapbad

TrackerMap *	tkmapden

TrackerMap *	tkmapeff

TrackerMap *	tkmapnum

Static Private Attributes
static constexpr int	SiStripLayers = 22

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from ConditionDBWriter< SiStripBadStrip >
static void	fillPSetDescription (edm::ParameterSetDescription &desc)

Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from ConditionDBWriter< SiStripBadStrip >
void	setDoStore (const bool doStore)
	When set to false the payload will not be written to the db. More...

void	storeOnDbNow ()

cond::Time_t	timeOfLastIOV ()

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Definition at line 95 of file SiStripHitEffFromCalibTree.cc.

Constructor & Destructor Documentation

◆ SiStripHitEffFromCalibTree()

SiStripHitEffFromCalibTree::SiStripHitEffFromCalibTree ( const edm::ParameterSet & conf )

explicit

Definition at line 189 of file SiStripHitEffFromCalibTree.cc.

References _autoIneffModTagging, _badModulesFile, _bunchx, _clusterMatchingMethod, _clusterTrajDist, _detInfo, _effPlotMin, _ResXSig, _showEndcapSides, _showOnlyGoodModules, _showRings, _showTOB6TEC9, _spaceBetweenTrains, _stripsApvEdge, _title, _tkGeomToken, _tkMapMin, _tTopoToken, _useCM, _useOnlyHighPurityTracks, CalibTreeFilenames, doSummary, edm::EDConsumerBase::esConsumes(), FileInPath_, edm::FileInPath::fullPath(), edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), TFileService::kSharedResource, nModsMin, nTEClayers, quality_, SiStripDetInfoFileReader::read(), AlCaHarvesting_cff::SiStripQuality, AlCaHLTBitMon_QueryRunRegistry::string, and threshold.

     : ConditionDBWriter<SiStripBadStrip>(conf), FileInPath_(SiStripDetInfoFileReader::kDefaultFile) {
   usesResource(TFileService::kSharedResource);
   CalibTreeFilenames = conf.getUntrackedParameter<vector<std::string> >("CalibTreeFilenames");
   threshold = conf.getParameter<double>("Threshold");
   nModsMin = conf.getParameter<int>("nModsMin");
   doSummary = conf.getParameter<int>("doSummary");
   _badModulesFile = conf.getUntrackedParameter<std::string>("BadModulesFile", "");
   _autoIneffModTagging = conf.getUntrackedParameter<bool>("AutoIneffModTagging", false);
   _clusterMatchingMethod = conf.getUntrackedParameter<int>("ClusterMatchingMethod", 0);
   _ResXSig = conf.getUntrackedParameter<double>("ResXSig", -1);
   _clusterTrajDist = conf.getUntrackedParameter<double>("ClusterTrajDist", 64.0);
   _stripsApvEdge = conf.getUntrackedParameter<double>("StripsApvEdge", 10.0);
   _useOnlyHighPurityTracks = conf.getUntrackedParameter<bool>("UseOnlyHighPurityTracks", true);
   _bunchx = conf.getUntrackedParameter<int>("BunchCrossing", 0);
   _spaceBetweenTrains = conf.getUntrackedParameter<int>("SpaceBetweenTrains", 25);
   _useCM = conf.getUntrackedParameter<bool>("UseCommonMode", false);
   _showEndcapSides = conf.getUntrackedParameter<bool>("ShowEndcapSides", true);
   _showRings = conf.getUntrackedParameter<bool>("ShowRings", false);
   _showTOB6TEC9 = conf.getUntrackedParameter<bool>("ShowTOB6TEC9", false);
   _showOnlyGoodModules = conf.getUntrackedParameter<bool>("ShowOnlyGoodModules", false);
   _tkMapMin = conf.getUntrackedParameter<double>("TkMapMin", 0.9);
   _effPlotMin = conf.getUntrackedParameter<double>("EffPlotMin", 0.9);
   _title = conf.getParameter<std::string>("Title");
   _tkGeomToken = esConsumes();
   _tTopoToken = esConsumes();
   _detInfo = SiStripDetInfoFileReader::read(FileInPath_.fullPath());
 
   nTEClayers = 9;  // number of wheels
   if (_showRings)
     nTEClayers = 7;  // number of rings
 
   quality_ = new SiStripQuality(_detInfo);
 }

◆ ~SiStripHitEffFromCalibTree()

SiStripHitEffFromCalibTree::~SiStripHitEffFromCalibTree ( )

overridedefault

Member Function Documentation

◆ algoAnalyze()

void SiStripHitEffFromCalibTree::algoAnalyze	(	const edm::Event &	e,
		const edm::EventSetup &	c
	)

overrideprivatevirtual

Reimplemented from ConditionDBWriter< SiStripBadStrip >.

Definition at line 228 of file SiStripHitEffFromCalibTree.cc.

References _autoIneffModTagging, _badModulesFile, _bunchx, _clusterMatchingMethod, _clusterTrajDist, _detInfo, _ResXSig, _showRings, _showTOB6TEC9, _stripsApvEdge, _tkGeomToken, _tTopoToken, _useCM, _useOnlyHighPurityTracks, a, funct::abs(), accept(), alllayerfound, alllayertotal, electrons_cff::bool, Surface::bounds(), simKBmtfDigis_cfi::bx, bxHisto, c, CalibTree, CalibTreeFilenames, SiStripBadStrip::decode(), MillePedeFileConverter_cfg::e, mps_fire::end, eventInfos, cppFunctionSkipper::exception, fs, SiStripQuality::getBadComponentList(), SiStripBadStrip::getDataVectorBegin(), SiStripDetInfo::getNumberOfApvsAndStripLength(), SiStripBadStrip::getRegistryVectorBegin(), SiStripBadStrip::getRegistryVectorEnd(), goodlayerfound, goodlayertotal, muons_cff::highPurity, hits, mps_fire::i, l1ctLayer2EG_cff::id, hit::id, compare_using_db::ifile, muonGEMDigis_cfi::instLumi, instLumiHisto, createfilelist::int, dqmiolumiharvest::j, dqmdumpme::k, cmsLHEtoEOSManager::l, phase1PixelTopology::layer, layerfound, layerfound_perBx, layerfound_vsCM, layerfound_vsLumi, layerfound_vsPU, phase1PixelTopology::layerName, layertotal, layertotal_perBx, layertotal_vsCM, layertotal_vsLumi, layertotal_vsPU, mps_splice::line, LOGPRINT, TFileService::make(), makeHotColdMaps(), makeSQLite(), makeSummary(), makeSummaryVsBx(), makeSummaryVsCM(), makeSummaryVsLumi(), makeTKMap(), modCounter, mps_splice::mods, createIOVlist::nevents, nevt, nModsMin, StripTopology::nstrips(), nTEClayers, MLClient_cfi::PU, PUHisto, quality, quality_, SiStripBadStrip::data::range, FastTimerService_cff::range, writedatasetfile::run, setBadComponents(), SiStripLayers, StripGeomDetUnit::specificTopology(), contentValuesCheck::ss, str, AlCaHLTBitMon_QueryRunRegistry::string, TrackingMonitor_cfi::stripCluster, GeomDet::surface(), SiStripDetId::TEC, threshold, SiStripDetId::TIB, SiStripDetId::TID, SiStripDetId::TOB, totalStatistics(), Bounds::width(), x, hit::x, testProducerWithPsetDescEmpty_cfi::x1, y, hit::y, z, and hit::z.

                                                                                       {
   const auto& tkgeom = c.getData(_tkGeomToken);
   const auto& tTopo = c.getData(_tTopoToken);
 
   // read bad modules to mask
   ifstream badModules_file;
   set<uint32_t> badModules_list;
   if (!_badModulesFile.empty()) {
     badModules_file.open(_badModulesFile.c_str());
     uint32_t badmodule_detid;
     int mods, fiber1, fiber2, fiber3;
     if (badModules_file.is_open()) {
       string line;
       while (getline(badModules_file, line)) {
         if (badModules_file.eof())
           continue;
         stringstream ss(line);
         ss >> badmodule_detid >> mods >> fiber1 >> fiber2 >> fiber3;
         if (badmodule_detid != 0 && mods == 1 && (fiber1 == 1 || fiber2 == 1 || fiber3 == 1))
           badModules_list.insert(badmodule_detid);
       }
       badModules_file.close();
     }
   }
   if (!badModules_list.empty())
     LOGPRINT << "Remove additionnal bad modules from the analysis: ";
   set<uint32_t>::iterator itBadMod;
   for (itBadMod = badModules_list.begin(); itBadMod != badModules_list.end(); ++itBadMod)
     LOGPRINT << " " << *itBadMod;
 
   // initialze counters and histos
 
   bxHisto = fs->make<TH1F>("bx", "bx", 3600, 0, 3600);
   instLumiHisto = fs->make<TH1F>("instLumi", "inst. lumi.", 250, 0, 25000);
   PUHisto = fs->make<TH1F>("PU", "PU", 200, 0, 200);
 
   for (int l = 0; l < 35; l++) {
     goodlayertotal[l] = 0;
     goodlayerfound[l] = 0;
     alllayertotal[l] = 0;
     alllayerfound[l] = 0;
   }
 
   TH1F* resolutionPlots[23];
   for (Long_t ilayer = 0; ilayer < 23; ilayer++) {
     std::string lyrName = ::layerName(ilayer, _showRings, nTEClayers);
 
     resolutionPlots[ilayer] = fs->make<TH1F>(Form("resol_layer_%i", (int)(ilayer)), lyrName.c_str(), 125, -125, 125);
     resolutionPlots[ilayer]->GetXaxis()->SetTitle("trajX-clusX [strip unit]");
 
     layerfound_vsLumi.push_back(
         fs->make<TH1F>(Form("layerfound_vsLumi_layer_%i", (int)(ilayer)), lyrName.c_str(), 100, 0, 25000));
     layertotal_vsLumi.push_back(
         fs->make<TH1F>(Form("layertotal_vsLumi_layer_%i", (int)(ilayer)), lyrName.c_str(), 100, 0, 25000));
     layerfound_vsPU.push_back(
         fs->make<TH1F>(Form("layerfound_vsPU_layer_%i", (int)(ilayer)), lyrName.c_str(), 45, 0, 90));
     layertotal_vsPU.push_back(
         fs->make<TH1F>(Form("layertotal_vsPU_layer_%i", (int)(ilayer)), lyrName.c_str(), 45, 0, 90));
 
     if (_useCM) {
       layerfound_vsCM.push_back(
           fs->make<TH1F>(Form("layerfound_vsCM_layer_%i", (int)(ilayer)), lyrName.c_str(), 20, 0, 400));
       layertotal_vsCM.push_back(
           fs->make<TH1F>(Form("layertotal_vsCM_layer_%i", (int)(ilayer)), lyrName.c_str(), 20, 0, 400));
     }
     layertotal[ilayer] = 0;
     layerfound[ilayer] = 0;
   }
 
   if (!_autoIneffModTagging)
     LOGPRINT << "A module is bad if efficiency < " << threshold << " and has at least " << nModsMin << " nModsMin.";
   else
     LOGPRINT << "A module is bad if the upper limit on the efficiency is < to the avg in the layer - " << threshold
              << " and has at least " << nModsMin << " nModsMin.";
 
   unsigned int run, evt, bx{0};
   double instLumi, PU;
 
   //Open the ROOT Calib Tree
   for (unsigned int ifile = 0; ifile < CalibTreeFilenames.size(); ifile++) {
     LOGPRINT << "Loading file: " << CalibTreeFilenames[ifile];
     TFile* CalibTreeFile = TFile::Open(CalibTreeFilenames[ifile].c_str(), "READ");
 
     // Get event infos
     bool foundEventInfos = false;
     try {
       CalibTreeFile->cd("eventInfo");
     } catch (exception& e) {
       LOGPRINT << "No event infos tree";
     }
     TTree* EventTree = (TTree*)(gDirectory->Get("tree"));
 
     TLeaf* runLf;
     TLeaf* evtLf;
     TLeaf* BunchLf;
     TLeaf* InstLumiLf;
     TLeaf* PULf;
     if (EventTree) {
       LOGPRINT << "Found event infos tree";
 
       runLf = EventTree->GetLeaf("run");
       evtLf = EventTree->GetLeaf("event");
 
       BunchLf = EventTree->GetLeaf("bx");
       InstLumiLf = EventTree->GetLeaf("instLumi");
       PULf = EventTree->GetLeaf("PU");
 
       int nevt = EventTree->GetEntries();
       if (nevt)
         foundEventInfos = true;
 
       for (int j = 0; j < nevt; j++) {
         EventTree->GetEntry(j);
         run = runLf->GetValue();
         evt = evtLf->GetValue();
         bx = BunchLf->GetValue();
         instLumi = InstLumiLf->GetValue();
         PU = PULf->GetValue();
 
         bxHisto->Fill(bx);
         instLumiHisto->Fill(instLumi);
         PUHisto->Fill(PU);
 
         eventInfos[make_pair(run, evt)] = array<double, 3>{{(double)bx, instLumi, PU}};
       }
     }
 
     // Get hit infos
     CalibTreeFile->cd("anEff");
     CalibTree = (TTree*)(gDirectory->Get("traj"));
 
     runLf = CalibTree->GetLeaf("run");
     evtLf = CalibTree->GetLeaf("event");
     TLeaf* BadLf = CalibTree->GetLeaf("ModIsBad");
     TLeaf* sistripLf = CalibTree->GetLeaf("SiStripQualBad");
     TLeaf* idLf = CalibTree->GetLeaf("Id");
     TLeaf* acceptLf = CalibTree->GetLeaf("withinAcceptance");
     TLeaf* layerLf = CalibTree->GetLeaf("layer");
     //TLeaf* nHitsLf = CalibTree->GetLeaf("nHits");
     TLeaf* highPurityLf = CalibTree->GetLeaf("highPurity");
     TLeaf* xLf = CalibTree->GetLeaf("TrajGlbX");
     TLeaf* yLf = CalibTree->GetLeaf("TrajGlbY");
     TLeaf* zLf = CalibTree->GetLeaf("TrajGlbZ");
     TLeaf* ResXSigLf = CalibTree->GetLeaf("ResXSig");
     TLeaf* TrajLocXLf = CalibTree->GetLeaf("TrajLocX");
     TLeaf* TrajLocYLf = CalibTree->GetLeaf("TrajLocY");
     TLeaf* ClusterLocXLf = CalibTree->GetLeaf("ClusterLocX");
     BunchLf = CalibTree->GetLeaf("bunchx");
     InstLumiLf = CalibTree->GetLeaf("instLumi");
     PULf = CalibTree->GetLeaf("PU");
     TLeaf* CMLf = nullptr;
     if (_useCM)
       CMLf = CalibTree->GetLeaf("commonMode");
 
     int nevents = CalibTree->GetEntries();
     LOGPRINT << "Successfully loaded analyze function with " << nevents << " events!\n";
 
     map<pair<unsigned int, unsigned int>, array<double, 3> >::iterator itEventInfos;
 
     //Loop through all of the events
     for (int j = 0; j < nevents; j++) {
       CalibTree->GetEntry(j);
       run = (unsigned int)runLf->GetValue();
       evt = (unsigned int)evtLf->GetValue();
       unsigned int isBad = (unsigned int)BadLf->GetValue();
       unsigned int quality = (unsigned int)sistripLf->GetValue();
       unsigned int id = (unsigned int)idLf->GetValue();
       unsigned int accept = (unsigned int)acceptLf->GetValue();
       unsigned int layer_wheel = (unsigned int)layerLf->GetValue();
       unsigned int layer = layer_wheel;
       if (_showRings && layer > 10) {  // use rings instead of wheels
         if (layer < 14)
           layer = 10 + ((id >> 9) & 0x3);  //TID   3 disks and also 3 rings -> use the same container
         else
           layer = 13 + ((id >> 5) & 0x7);  //TEC
       }
       //unsigned int nHits = (unsigned int)nHitsLf->GetValue();
       bool highPurity = (bool)highPurityLf->GetValue();
       double x = xLf->GetValue();
       double y = yLf->GetValue();
       double z = zLf->GetValue();
       double resxsig = ResXSigLf->GetValue();
       double TrajLocX = TrajLocXLf->GetValue();
       double TrajLocY = TrajLocYLf->GetValue();
       double ClusterLocX = ClusterLocXLf->GetValue();
       double TrajLocXMid;
       double stripTrajMid;
       double stripCluster;
       bool badquality = false;
 
       instLumi = 0;
       PU = 0;
 
       // if no special tree with event infos, they may be stored in the hit eff tree
       if (!foundEventInfos) {
         bx = (unsigned int)BunchLf->GetValue();
         if (InstLumiLf != nullptr)
           instLumi = InstLumiLf->GetValue();  // branch not filled by default
         if (PULf != nullptr)
           PU = PULf->GetValue();  // branch not filled by default
       }
       int CM = -100;
       if (_useCM)
         CM = CMLf->GetValue();
 
       // Get infos from eventInfos if they exist
       if (foundEventInfos) {
         itEventInfos = eventInfos.find(make_pair(run, evt));
         if (itEventInfos != eventInfos.end()) {
           bx = itEventInfos->second[0];
           instLumi = itEventInfos->second[1];
           PU = itEventInfos->second[2];
         }
       }
 
       //We have two things we want to do, both an XY color plot, and the efficiency measurement
       //First, ignore anything that isn't in acceptance and isn't good quality
 
       if (_bunchx > 0 && _bunchx != bx)
         continue;
 
       //if(quality == 1 || accept != 1 || nHits < 8) continue;
       if (accept != 1)
         continue;
       if (_useOnlyHighPurityTracks && !highPurity)
         continue;
       if (quality == 1)
         badquality = true;
 
       // don't compute efficiencies in modules from TOB6 and TEC9
       if (!_showTOB6TEC9 && (layer_wheel == 10 || layer_wheel == SiStripLayers))
         continue;
 
       // don't use bad modules given in the bad module list
       itBadMod = badModules_list.find(id);
       if (itBadMod != badModules_list.end())
         continue;
 
       //Now that we have a good event, we need to look at if we expected it or not, and the location
       //if we didn't
       //Fill the missing hit information first
       bool badflag = false;
 
       // By default uses the old matching method
       if (_ResXSig < 0) {
         if (isBad == 1)
           badflag = true;  // isBad set to false in the tree when resxsig<999.0
       } else {
         if (isBad == 1 || resxsig > _ResXSig)
           badflag = true;
       }
 
       // Conversion of positions in strip unit
       int nstrips = -9;
       float Pitch = -9.0;
 
       if (resxsig == 1000.0) {  // special treatment, no GeomDetUnit associated in some cases when no cluster found
         Pitch = 0.0205;         // maximum
         nstrips = 768;          // maximum
         stripTrajMid = TrajLocX / Pitch + nstrips / 2.0;
         stripCluster = ClusterLocX / Pitch + nstrips / 2.0;
       } else {
         DetId ClusterDetId(id);
         const StripGeomDetUnit* stripdet = (const StripGeomDetUnit*)tkgeom.idToDetUnit(ClusterDetId);
         const StripTopology& Topo = stripdet->specificTopology();
         nstrips = Topo.nstrips();
         Pitch = stripdet->surface().bounds().width() / Topo.nstrips();
         stripTrajMid = TrajLocX / Pitch + nstrips / 2.0;  //layer01->10
         stripCluster = ClusterLocX / Pitch + nstrips / 2.0;
 
         // For trapezoidal modules: extrapolation of x trajectory position to the y middle of the module
         //  for correct comparison with cluster position
         float hbedge = 0;
         float htedge = 0;
         float hapoth = 0;
         if (layer >= 11) {
           const BoundPlane& plane = stripdet->surface();
           const TrapezoidalPlaneBounds* trapezoidalBounds(
               dynamic_cast<const TrapezoidalPlaneBounds*>(&(plane.bounds())));
           std::array<const float, 4> const& parameters = (*trapezoidalBounds).parameters();
           hbedge = parameters[0];
           htedge = parameters[1];
           hapoth = parameters[3];
           TrajLocXMid = TrajLocX / (1 + (htedge - hbedge) * TrajLocY / (htedge + hbedge) /
                                             hapoth);  // radialy extrapolated x loc position at middle
           stripTrajMid = TrajLocXMid / Pitch + nstrips / 2.0;
         }
       }
 
       if (!badquality && layer < 23) {
         if (resxsig != 1000.0)
           resolutionPlots[layer]->Fill(stripTrajMid - stripCluster);
         else
           resolutionPlots[layer]->Fill(1000);
       }
 
       // New matching methods
       int tapv = -9;
       int capv = -9;
       float stripInAPV = 64.;
 
       if (_clusterMatchingMethod >= 1) {
         badflag = false;          // reset
         if (resxsig == 1000.0) {  // default value when no cluster found in the module
           badflag = true;         // consider the module inefficient in this case
         } else {
           if (_clusterMatchingMethod == 2 ||
               _clusterMatchingMethod == 4) {  // check the distance between cluster and trajectory position
             if (abs(stripCluster - stripTrajMid) > _clusterTrajDist)
               badflag = true;
           }
           if (_clusterMatchingMethod == 3 ||
               _clusterMatchingMethod ==
                   4) {  // cluster and traj have to be in the same APV (don't take edges into accounts)
             tapv = (int)stripTrajMid / 128;
             capv = (int)stripCluster / 128;
             stripInAPV = stripTrajMid - tapv * 128;
 
             if (stripInAPV < _stripsApvEdge || stripInAPV > 128 - _stripsApvEdge)
               continue;
             if (tapv != capv)
               badflag = true;
           }
         }
       }
 
       if (badflag && !badquality) {
         hit temphit;
         temphit.x = x;
         temphit.y = y;
         temphit.z = z;
         temphit.id = id;
         hits[layer].push_back(temphit);
       }
       pair<unsigned int, unsigned int> newgoodpair(1, 1);
       pair<unsigned int, unsigned int> newbadpair(1, 0);
       //First, figure out if the module already exists in the map of maps
       map<unsigned int, pair<unsigned int, unsigned int> >::iterator it = modCounter[layer].find(id);
       if (!badquality) {
         if (it == modCounter[layer].end()) {
           if (badflag)
             modCounter[layer][id] = newbadpair;
           else
             modCounter[layer][id] = newgoodpair;
         } else {
           ((*it).second.first)++;
           if (!badflag)
             ((*it).second.second)++;
         }
 
         if (layerfound_perBx.find(bx) == layerfound_perBx.end()) {
           layerfound_perBx[bx] = vector<int>(23, 0);
           layertotal_perBx[bx] = vector<int>(23, 0);
         }
         if (!badflag)
           layerfound_perBx[bx][layer]++;
         layertotal_perBx[bx][layer]++;
 
         if (!badflag)
           layerfound_vsLumi[layer]->Fill(instLumi);
         layertotal_vsLumi[layer]->Fill(instLumi);
         if (!badflag)
           layerfound_vsPU[layer]->Fill(PU);
         layertotal_vsPU[layer]->Fill(PU);
 
         if (_useCM) {
           if (!badflag)
             layerfound_vsCM[layer]->Fill(CM);
           layertotal_vsCM[layer]->Fill(CM);
         }
 
         //Have to do the decoding for which side to go on (ugh)
         if (layer <= 10) {
           if (!badflag)
             goodlayerfound[layer]++;
           goodlayertotal[layer]++;
         } else if (layer > 10 && layer < 14) {
           if (((id >> 13) & 0x3) == 1) {
             if (!badflag)
               goodlayerfound[layer]++;
             goodlayertotal[layer]++;
           } else if (((id >> 13) & 0x3) == 2) {
             if (!badflag)
               goodlayerfound[layer + 3]++;
             goodlayertotal[layer + 3]++;
           }
         } else if (layer > 13 && layer <= SiStripLayers) {
           if (((id >> 18) & 0x3) == 1) {
             if (!badflag)
               goodlayerfound[layer + 3]++;
             goodlayertotal[layer + 3]++;
           } else if (((id >> 18) & 0x3) == 2) {
             if (!badflag)
               goodlayerfound[layer + 3 + nTEClayers]++;
             goodlayertotal[layer + 3 + nTEClayers]++;
           }
         }
       }
       //Do the one where we don't exclude bad modules!
       if (layer <= 10) {
         if (!badflag)
           alllayerfound[layer]++;
         alllayertotal[layer]++;
       } else if (layer > 10 && layer < 14) {
         if (((id >> 13) & 0x3) == 1) {
           if (!badflag)
             alllayerfound[layer]++;
           alllayertotal[layer]++;
         } else if (((id >> 13) & 0x3) == 2) {
           if (!badflag)
             alllayerfound[layer + 3]++;
           alllayertotal[layer + 3]++;
         }
       } else if (layer > 13 && layer <= SiStripLayers) {
         if (((id >> 18) & 0x3) == 1) {
           if (!badflag)
             alllayerfound[layer + 3]++;
           alllayertotal[layer + 3]++;
         } else if (((id >> 18) & 0x3) == 2) {
           if (!badflag)
             alllayerfound[layer + 3 + nTEClayers]++;
           alllayertotal[layer + 3 + nTEClayers]++;
         }
       }
       //At this point, both of our maps are loaded with the correct information
     }
   }  // go to next CalibTreeFile
 
   makeHotColdMaps();
   makeTKMap(_autoIneffModTagging);
   makeSQLite();
   totalStatistics();
   makeSummary();
   makeSummaryVsBx();
   makeSummaryVsLumi();
   if (_useCM)
     makeSummaryVsCM();
 
   //try to write out what's in the quality record
   int NTkBadComponent[4];  //k: 0=BadModule, 1=BadFiber, 2=BadApv, 3=BadStrips
   int NBadComponent[4][19][4];
   //legend: NBadComponent[i][j][k]= SubSystem i, layer/disk/wheel j, BadModule/Fiber/Apv k
   //     i: 0=TIB, 1=TID, 2=TOB, 3=TEC
   //     k: 0=BadModule, 1=BadFiber, 2=BadApv, 3=BadStrips
   std::stringstream ssV[4][19];
 
   for (int i = 0; i < 4; ++i) {
     NTkBadComponent[i] = 0;
     for (int j = 0; j < 19; ++j) {
       ssV[i][j].str("");
       for (int k = 0; k < 4; ++k)
         NBadComponent[i][j][k] = 0;
     }
   }
 
   std::vector<SiStripQuality::BadComponent> BC = quality_->getBadComponentList();
 
   for (size_t i = 0; i < BC.size(); ++i) {
     //&&&&&&&&&&&&&
     //Full Tk
     //&&&&&&&&&&&&&
 
     if (BC[i].BadModule)
       NTkBadComponent[0]++;
     if (BC[i].BadFibers)
       NTkBadComponent[1] += ((BC[i].BadFibers >> 2) & 0x1) + ((BC[i].BadFibers >> 1) & 0x1) + ((BC[i].BadFibers) & 0x1);
     if (BC[i].BadApvs)
       NTkBadComponent[2] += ((BC[i].BadApvs >> 5) & 0x1) + ((BC[i].BadApvs >> 4) & 0x1) + ((BC[i].BadApvs >> 3) & 0x1) +
                             ((BC[i].BadApvs >> 2) & 0x1) + ((BC[i].BadApvs >> 1) & 0x1) + ((BC[i].BadApvs) & 0x1);
 
     //&&&&&&&&&&&&&&&&&
     //Single SubSystem
     //&&&&&&&&&&&&&&&&&
 
     int component;
     SiStripDetId a(BC[i].detid);
     if (a.subdetId() == SiStripDetId::TIB) {
       //&&&&&&&&&&&&&&&&&
       //TIB
       //&&&&&&&&&&&&&&&&&
 
       component = tTopo.tibLayer(BC[i].detid);
       setBadComponents(0, component, BC[i], ssV, NBadComponent);
 
     } else if (a.subdetId() == SiStripDetId::TID) {
       //&&&&&&&&&&&&&&&&&
       //TID
       //&&&&&&&&&&&&&&&&&
 
       component = tTopo.tidSide(BC[i].detid) == 2 ? tTopo.tidWheel(BC[i].detid) : tTopo.tidWheel(BC[i].detid) + 3;
       setBadComponents(1, component, BC[i], ssV, NBadComponent);
 
     } else if (a.subdetId() == SiStripDetId::TOB) {
       //&&&&&&&&&&&&&&&&&
       //TOB
       //&&&&&&&&&&&&&&&&&
 
       component = tTopo.tobLayer(BC[i].detid);
       setBadComponents(2, component, BC[i], ssV, NBadComponent);
 
     } else if (a.subdetId() == SiStripDetId::TEC) {
       //&&&&&&&&&&&&&&&&&
       //TEC
       //&&&&&&&&&&&&&&&&&
 
       component = tTopo.tecSide(BC[i].detid) == 2 ? tTopo.tecWheel(BC[i].detid) : tTopo.tecWheel(BC[i].detid) + 9;
       setBadComponents(3, component, BC[i], ssV, NBadComponent);
     }
   }
 
   //&&&&&&&&&&&&&&&&&&
   // Single Strip Info
   //&&&&&&&&&&&&&&&&&&
   float percentage = 0;
 
   SiStripQuality::RegistryIterator rbegin = quality_->getRegistryVectorBegin();
   SiStripQuality::RegistryIterator rend = quality_->getRegistryVectorEnd();
 
   for (SiStripBadStrip::RegistryIterator rp = rbegin; rp != rend; ++rp) {
     unsigned int detid = rp->detid;
 
     int subdet = -999;
     int component = -999;
     SiStripDetId a(detid);
     if (a.subdetId() == 3) {
       subdet = 0;
       component = tTopo.tibLayer(detid);
     } else if (a.subdetId() == 4) {
       subdet = 1;
       component = tTopo.tidSide(detid) == 2 ? tTopo.tidWheel(detid) : tTopo.tidWheel(detid) + 3;
     } else if (a.subdetId() == 5) {
       subdet = 2;
       component = tTopo.tobLayer(detid);
     } else if (a.subdetId() == 6) {
       subdet = 3;
       component = tTopo.tecSide(detid) == 2 ? tTopo.tecWheel(detid) : tTopo.tecWheel(detid) + 9;
     }
 
     SiStripQuality::Range sqrange =
         SiStripQuality::Range(quality_->getDataVectorBegin() + rp->ibegin, quality_->getDataVectorBegin() + rp->iend);
 
     percentage = 0;
     for (int it = 0; it < sqrange.second - sqrange.first; it++) {
       unsigned int range = quality_->decode(*(sqrange.first + it)).range;
       NTkBadComponent[3] += range;
       NBadComponent[subdet][0][3] += range;
       NBadComponent[subdet][component][3] += range;
       percentage += range;
     }
     if (percentage != 0)
       percentage /= 128. * _detInfo.getNumberOfApvsAndStripLength(detid).first;
     if (percentage > 1)
       edm::LogError("SiStripQualityStatistics") << "PROBLEM detid " << detid << " value " << percentage << std::endl;
   }
   //&&&&&&&&&&&&&&&&&&
   // printout
   //&&&&&&&&&&&&&&&&&&
   std::ostringstream ss;
 
   ss << "\n-----------------\nNew IOV starting from run " << e.id().run() << " event " << e.id().event()
      << " lumiBlock " << e.luminosityBlock() << " time " << e.time().value() << "\n-----------------\n";
   ss << "\n-----------------\nGlobal Info\n-----------------";
   ss << "\nBadComponent \t  Modules \tFibers "
         "\tApvs\tStrips\n----------------------------------------------------------------";
   ss << "\nTracker:\t\t" << NTkBadComponent[0] << "\t" << NTkBadComponent[1] << "\t" << NTkBadComponent[2] << "\t"
      << NTkBadComponent[3];
   ss << "\nTIB:\t\t\t" << NBadComponent[0][0][0] << "\t" << NBadComponent[0][0][1] << "\t" << NBadComponent[0][0][2]
      << "\t" << NBadComponent[0][0][3];
   ss << "\nTID:\t\t\t" << NBadComponent[1][0][0] << "\t" << NBadComponent[1][0][1] << "\t" << NBadComponent[1][0][2]
      << "\t" << NBadComponent[1][0][3];
   ss << "\nTOB:\t\t\t" << NBadComponent[2][0][0] << "\t" << NBadComponent[2][0][1] << "\t" << NBadComponent[2][0][2]
      << "\t" << NBadComponent[2][0][3];
   ss << "\nTEC:\t\t\t" << NBadComponent[3][0][0] << "\t" << NBadComponent[3][0][1] << "\t" << NBadComponent[3][0][2]
      << "\t" << NBadComponent[3][0][3];
   ss << "\n";
 
   for (int i = 1; i < 5; ++i)
     ss << "\nTIB Layer " << i << " :\t\t" << NBadComponent[0][i][0] << "\t" << NBadComponent[0][i][1] << "\t"
        << NBadComponent[0][i][2] << "\t" << NBadComponent[0][i][3];
   ss << "\n";
   for (int i = 1; i < 4; ++i)
     ss << "\nTID+ Disk " << i << " :\t\t" << NBadComponent[1][i][0] << "\t" << NBadComponent[1][i][1] << "\t"
        << NBadComponent[1][i][2] << "\t" << NBadComponent[1][i][3];
   for (int i = 4; i < 7; ++i)
     ss << "\nTID- Disk " << i - 3 << " :\t\t" << NBadComponent[1][i][0] << "\t" << NBadComponent[1][i][1] << "\t"
        << NBadComponent[1][i][2] << "\t" << NBadComponent[1][i][3];
   ss << "\n";
   for (int i = 1; i < 7; ++i)
     ss << "\nTOB Layer " << i << " :\t\t" << NBadComponent[2][i][0] << "\t" << NBadComponent[2][i][1] << "\t"
        << NBadComponent[2][i][2] << "\t" << NBadComponent[2][i][3];
   ss << "\n";
   for (int i = 1; i < 10; ++i)
     ss << "\nTEC+ Disk " << i << " :\t\t" << NBadComponent[3][i][0] << "\t" << NBadComponent[3][i][1] << "\t"
        << NBadComponent[3][i][2] << "\t" << NBadComponent[3][i][3];
   for (int i = 10; i < 19; ++i)
     ss << "\nTEC- Disk " << i - 9 << " :\t\t" << NBadComponent[3][i][0] << "\t" << NBadComponent[3][i][1] << "\t"
        << NBadComponent[3][i][2] << "\t" << NBadComponent[3][i][3];
   ss << "\n";
 
   ss << "\n----------------------------------------------------------------\n\t\t   Detid  \tModules Fibers "
         "Apvs\n----------------------------------------------------------------";
   for (int i = 1; i < 5; ++i)
     ss << "\nTIB Layer " << i << " :" << ssV[0][i].str();
   ss << "\n";
   for (int i = 1; i < 4; ++i)
     ss << "\nTID+ Disk " << i << " :" << ssV[1][i].str();
   for (int i = 4; i < 7; ++i)
     ss << "\nTID- Disk " << i - 3 << " :" << ssV[1][i].str();
   ss << "\n";
   for (int i = 1; i < 7; ++i)
     ss << "\nTOB Layer " << i << " :" << ssV[2][i].str();
   ss << "\n";
   for (int i = 1; i < 10; ++i)
     ss << "\nTEC+ Disk " << i << " :" << ssV[3][i].str();
   for (int i = 10; i < 19; ++i)
     ss << "\nTEC- Disk " << i - 9 << " :" << ssV[3][i].str();
 
   LOGPRINT << ss.str();
 
   // store also bad modules in log file
   ofstream badModules;
   badModules.open("BadModules.log");
   badModules << "\n----------------------------------------------------------------\n\t\t   Detid  \tModules Fibers "
                 "Apvs\n----------------------------------------------------------------";
   for (int i = 1; i < 5; ++i)
     badModules << "\nTIB Layer " << i << " :" << ssV[0][i].str();
   badModules << "\n";
   for (int i = 1; i < 4; ++i)
     badModules << "\nTID+ Disk " << i << " :" << ssV[1][i].str();
   for (int i = 4; i < 7; ++i)
     badModules << "\nTID- Disk " << i - 3 << " :" << ssV[1][i].str();
   badModules << "\n";
   for (int i = 1; i < 7; ++i)
     badModules << "\nTOB Layer " << i << " :" << ssV[2][i].str();
   badModules << "\n";
   for (int i = 1; i < 10; ++i)
     badModules << "\nTEC+ Disk " << i << " :" << ssV[3][i].str();
   for (int i = 10; i < 19; ++i)
     badModules << "\nTEC- Disk " << i - 9 << " :" << ssV[3][i].str();
   badModules.close();
 }

◆ algoBeginJob()

void SiStripHitEffFromCalibTree::algoBeginJob ( const edm::EventSetup & )

overrideprivatevirtual

Reimplemented from ConditionDBWriter< SiStripBadStrip >.

Definition at line 224 of file SiStripHitEffFromCalibTree.cc.

224 {}

◆ algoEndJob()

void SiStripHitEffFromCalibTree::algoEndJob ( )

overrideprivatevirtual

Reimplemented from ConditionDBWriter< SiStripBadStrip >.

Definition at line 226 of file SiStripHitEffFromCalibTree.cc.

226 {}

◆ computeEff()

void SiStripHitEffFromCalibTree::computeEff	(	vector< TH1F *> &	vhfound,
		vector< TH1F *> &	vhtotal,
		string	name
	)

private

Definition at line 1451 of file SiStripHitEffFromCalibTree.cc.

References _showRings, MillePedeFileConverter_cfg::e, dqm-mbProfile::format, fs, mps_fire::i, phase1PixelTopology::layerName, TFileService::make(), Skims_PA_cff::name, MuonTCMETValueMapProducer_cff::nLayers, nTEClayers, and SiStripLayers.

Referenced by makeSummaryVsCM(), and makeSummaryVsLumi().

                                                                                                        {
   unsigned int nLayers = SiStripLayers;
   if (_showRings)
     nLayers = 20;
 
   TH1F* hfound;
   TH1F* htotal;
 
   for (unsigned int ilayer = 1; ilayer < nLayers; ilayer++) {
     hfound = vhfound[ilayer];
     htotal = vhtotal[ilayer];
 
     hfound->Sumw2();
     htotal->Sumw2();
 
     // new ROOT version: TGraph::Divide don't handle null or negative values
     for (Long_t i = 0; i < hfound->GetNbinsX() + 1; ++i) {
       if (hfound->GetBinContent(i) == 0)
         hfound->SetBinContent(i, 1e-6);
       if (htotal->GetBinContent(i) == 0)
         htotal->SetBinContent(i, 1);
     }
 
     TGraphAsymmErrors* geff = fs->make<TGraphAsymmErrors>(hfound->GetNbinsX());
     geff->SetName(Form("%s_layer%i", name.c_str(), ilayer));
     geff->BayesDivide(hfound, htotal);
     if (name == "effVsLumi")
       geff->SetTitle(
           fmt::format("Hit Efficiency vs inst. lumi. - {}", ::layerName(ilayer, _showRings, nTEClayers)).c_str());
     if (name == "effVsPU")
       geff->SetTitle(fmt::format("Hit Efficiency vs pileup - {}", ::layerName(ilayer, _showRings, nTEClayers)).c_str());
     if (name == "effVsCM")
       geff->SetTitle(
           fmt::format("Hit Efficiency vs common Mode - {}", ::layerName(ilayer, _showRings, nTEClayers)).c_str());
     geff->SetMarkerStyle(20);
   }
 }

◆ getLayerSideName()

TString SiStripHitEffFromCalibTree::getLayerSideName ( Long_t k )

private

Definition at line 1532 of file SiStripHitEffFromCalibTree.cc.

References _showRings, dqmdumpme::k, and nTEClayers.

Referenced by makeSummary().

                                                              {
   TString layername = "";
   TString ringlabel = "D";
   if (_showRings)
     ringlabel = "R";
   if (k > 0 && k < 5) {
     layername = TString("TIB L") + k;
   } else if (k > 4 && k < 11) {
     layername = TString("TOB L") + (k - 4);
   } else if (k > 10 && k < 14) {
     layername = TString("TID- ") + ringlabel + (k - 10);
   } else if (k > 13 && k < 17) {
     layername = TString("TID+ ") + ringlabel + (k - 13);
   } else if (k > 16 && k < 17 + nTEClayers) {
     layername = TString("TEC- ") + ringlabel + (k - 16);
   } else if (k > 16 + nTEClayers) {
     layername = TString("TEC+ ") + ringlabel + (k - 16 - nTEClayers);
   }
 
   return layername;
 }

◆ getNewObject()

std::unique_ptr< SiStripBadStrip > SiStripHitEffFromCalibTree::getNewObject ( )

overrideprivatevirtual

Implements ConditionDBWriter< SiStripBadStrip >.

Definition at line 1554 of file SiStripHitEffFromCalibTree.cc.

References SiStripBadStrip::getDataVectorBegin(), SiStripBadStrip::getRegistryVectorBegin(), SiStripBadStrip::getRegistryVectorEnd(), getGTfromDQMFile::obj, quality_, and FastTimerService_cff::range.

                                                                         {
   //Need this for a Condition DB Writer
   //Initialize a return variable
   auto obj = std::make_unique<SiStripBadStrip>();
 
   SiStripBadStrip::RegistryIterator rIter = quality_->getRegistryVectorBegin();
   SiStripBadStrip::RegistryIterator rIterEnd = quality_->getRegistryVectorEnd();
 
   for (; rIter != rIterEnd; ++rIter) {
     SiStripBadStrip::Range range(quality_->getDataVectorBegin() + rIter->ibegin,
                                  quality_->getDataVectorBegin() + rIter->iend);
     if (!obj->put(rIter->detid, range))
       edm::LogError("SiStripHitEffFromCalibTree")
           << "[SiStripHitEffFromCalibTree::getNewObject] detid already exists" << std::endl;
   }
 
   return obj;
 }

◆ makeHotColdMaps()

void SiStripHitEffFromCalibTree::makeHotColdMaps ( )

private

Definition at line 872 of file SiStripHitEffFromCalibTree.cc.

References fs, hits, HotColdMaps, LOGPRINT, TFileService::make(), phi, and SiStripLayers.

Referenced by algoAnalyze().

                                                  {
   LOGPRINT << "Entering hot cold map generation!\n";
   TStyle* gStyle = new TStyle("gStyle", "myStyle");
   gStyle->cd();
   gStyle->SetPalette(1);
   gStyle->SetCanvasColor(kWhite);
   gStyle->SetOptStat(0);
   //Here we make the hot/cold color maps that we love so very much
   //Already have access to the data as a private variable
   //Create all of the histograms in the TFileService
   TH2F* temph2;
   for (Long_t maplayer = 1; maplayer <= SiStripLayers; maplayer++) {
     //Initialize all of the histograms
     if (maplayer > 0 && maplayer <= 4) {
       //We are in the TIB
       temph2 = fs->make<TH2F>(Form("%s%i", "TIB", (int)(maplayer)), "TIB", 100, -1, 361, 100, -100, 100);
       temph2->GetXaxis()->SetTitle("Phi");
       temph2->GetXaxis()->SetBinLabel(1, TString("360"));
       temph2->GetXaxis()->SetBinLabel(50, TString("180"));
       temph2->GetXaxis()->SetBinLabel(100, TString("0"));
       temph2->GetYaxis()->SetTitle("Global Z");
       temph2->SetOption("colz");
       HotColdMaps.push_back(temph2);
     } else if (maplayer > 4 && maplayer <= 10) {
       //We are in the TOB
       temph2 = fs->make<TH2F>(Form("%s%i", "TOB", (int)(maplayer - 4)), "TOB", 100, -1, 361, 100, -120, 120);
       temph2->GetXaxis()->SetTitle("Phi");
       temph2->GetXaxis()->SetBinLabel(1, TString("360"));
       temph2->GetXaxis()->SetBinLabel(50, TString("180"));
       temph2->GetXaxis()->SetBinLabel(100, TString("0"));
       temph2->GetYaxis()->SetTitle("Global Z");
       temph2->SetOption("colz");
       HotColdMaps.push_back(temph2);
     } else if (maplayer > 10 && maplayer <= 13) {
       //We are in the TID
       //Split by +/-
       temph2 = fs->make<TH2F>(Form("%s%i", "TID-", (int)(maplayer - 10)), "TID-", 100, -100, 100, 100, -100, 100);
       temph2->GetXaxis()->SetTitle("Global Y");
       temph2->GetXaxis()->SetBinLabel(1, TString("+Y"));
       temph2->GetXaxis()->SetBinLabel(50, TString("0"));
       temph2->GetXaxis()->SetBinLabel(100, TString("-Y"));
       temph2->GetYaxis()->SetTitle("Global X");
       temph2->GetYaxis()->SetBinLabel(1, TString("-X"));
       temph2->GetYaxis()->SetBinLabel(50, TString("0"));
       temph2->GetYaxis()->SetBinLabel(100, TString("+X"));
       temph2->SetOption("colz");
       HotColdMaps.push_back(temph2);
       temph2 = fs->make<TH2F>(Form("%s%i", "TID+", (int)(maplayer - 10)), "TID+", 100, -100, 100, 100, -100, 100);
       temph2->GetXaxis()->SetTitle("Global Y");
       temph2->GetXaxis()->SetBinLabel(1, TString("+Y"));
       temph2->GetXaxis()->SetBinLabel(50, TString("0"));
       temph2->GetXaxis()->SetBinLabel(100, TString("-Y"));
       temph2->GetYaxis()->SetTitle("Global X");
       temph2->GetYaxis()->SetBinLabel(1, TString("-X"));
       temph2->GetYaxis()->SetBinLabel(50, TString("0"));
       temph2->GetYaxis()->SetBinLabel(100, TString("+X"));
       temph2->SetOption("colz");
       HotColdMaps.push_back(temph2);
     } else if (maplayer > 13) {
       //We are in the TEC
       //Split by +/-
       temph2 = fs->make<TH2F>(Form("%s%i", "TEC-", (int)(maplayer - 13)), "TEC-", 100, -120, 120, 100, -120, 120);
       temph2->GetXaxis()->SetTitle("Global Y");
       temph2->GetXaxis()->SetBinLabel(1, TString("+Y"));
       temph2->GetXaxis()->SetBinLabel(50, TString("0"));
       temph2->GetXaxis()->SetBinLabel(100, TString("-Y"));
       temph2->GetYaxis()->SetTitle("Global X");
       temph2->GetYaxis()->SetBinLabel(1, TString("-X"));
       temph2->GetYaxis()->SetBinLabel(50, TString("0"));
       temph2->GetYaxis()->SetBinLabel(100, TString("+X"));
       temph2->SetOption("colz");
       HotColdMaps.push_back(temph2);
       temph2 = fs->make<TH2F>(Form("%s%i", "TEC+", (int)(maplayer - 13)), "TEC+", 100, -120, 120, 100, -120, 120);
       temph2->GetXaxis()->SetTitle("Global Y");
       temph2->GetXaxis()->SetBinLabel(1, TString("+Y"));
       temph2->GetXaxis()->SetBinLabel(50, TString("0"));
       temph2->GetXaxis()->SetBinLabel(100, TString("-Y"));
       temph2->GetYaxis()->SetTitle("Global X");
       temph2->GetYaxis()->SetBinLabel(1, TString("-X"));
       temph2->GetYaxis()->SetBinLabel(50, TString("0"));
       temph2->GetYaxis()->SetBinLabel(100, TString("+X"));
       temph2->SetOption("colz");
       HotColdMaps.push_back(temph2);
     }
   }
   for (Long_t mylayer = 1; mylayer <= SiStripLayers; mylayer++) {
     //Determine what kind of plot we want to write out
     //Loop through the entirety of each layer
     //Create an array of the histograms
     vector<hit>::const_iterator iter;
     for (iter = hits[mylayer].begin(); iter != hits[mylayer].end(); iter++) {
       //Looping over the particular layer
       //Fill by 360-x to get the proper location to compare with TKMaps of phi
       //Also global xy is messed up
       if (mylayer > 0 && mylayer <= 4) {
         //We are in the TIB
         float phi = ::calcPhi(iter->x, iter->y);
         HotColdMaps[mylayer - 1]->Fill(360. - phi, iter->z, 1.);
       } else if (mylayer > 4 && mylayer <= 10) {
         //We are in the TOB
         float phi = ::calcPhi(iter->x, iter->y);
         HotColdMaps[mylayer - 1]->Fill(360. - phi, iter->z, 1.);
       } else if (mylayer > 10 && mylayer <= 13) {
         //We are in the TID
         //There are 2 different maps here
         int side = (((iter->id) >> 13) & 0x3);
         if (side == 1)
           HotColdMaps[(mylayer - 1) + (mylayer - 11)]->Fill(-iter->y, iter->x, 1.);
         else if (side == 2)
           HotColdMaps[(mylayer - 1) + (mylayer - 10)]->Fill(-iter->y, iter->x, 1.);
         //if(side == 1) HotColdMaps[(mylayer - 1) + (mylayer - 11)]->Fill(iter->x,iter->y,1.);
         //else if(side == 2) HotColdMaps[(mylayer - 1) + (mylayer - 10)]->Fill(iter->x,iter->y,1.);
       } else if (mylayer > 13) {
         //We are in the TEC
         //There are 2 different maps here
         int side = (((iter->id) >> 18) & 0x3);
         if (side == 1)
           HotColdMaps[(mylayer + 2) + (mylayer - 14)]->Fill(-iter->y, iter->x, 1.);
         else if (side == 2)
           HotColdMaps[(mylayer + 2) + (mylayer - 13)]->Fill(-iter->y, iter->x, 1.);
         //if(side == 1) HotColdMaps[(mylayer + 2) + (mylayer - 14)]->Fill(iter->x,iter->y,1.);
         //else if(side == 2) HotColdMaps[(mylayer + 2) + (mylayer - 13)]->Fill(iter->x,iter->y,1.);
       }
     }
   }
   LOGPRINT << "Finished HotCold Map Generation\n";
 }

◆ makeSQLite()

void SiStripHitEffFromCalibTree::makeSQLite ( )

private

Definition at line 1110 of file SiStripHitEffFromCalibTree.cc.

References _detInfo, BadModules, SiStripQuality::compact(), SiStripQuality::fillBadComponents(), SiStripDetInfo::getNumberOfApvsAndStripLength(), globals_cff::id1, createfilelist::int, LOGPRINT, SiStripBadStrip::put(), quality_, and FastTimerService_cff::range.

Referenced by algoAnalyze().

                                             {
   //Generate the SQLite file for use in the Database of the bad modules!
   LOGPRINT << "Entering SQLite file generation!\n";
   std::vector<unsigned int> BadStripList;
   unsigned short NStrips;
   unsigned int id1;
   std::unique_ptr<SiStripQuality> pQuality = std::make_unique<SiStripQuality>(_detInfo);
   //This is the list of the bad strips, use to mask out entire APVs
   //Now simply go through the bad hit list and mask out things that
   //are bad!
   map<unsigned int, double>::const_iterator it;
   for (it = BadModules.begin(); it != BadModules.end(); it++) {
     //We need to figure out how many strips are in this particular module
     //To Mask correctly!
     NStrips = _detInfo.getNumberOfApvsAndStripLength((*it).first).first * 128;
     LOGPRINT << "Number of strips module " << (*it).first << " is " << NStrips;
     BadStripList.push_back(pQuality->encode(0, NStrips, 0));
     //Now compact into a single bad module
     id1 = (unsigned int)(*it).first;
     LOGPRINT << "ID1 shoudl match list of modules above " << id1;
     quality_->compact(id1, BadStripList);
     SiStripQuality::Range range(BadStripList.begin(), BadStripList.end());
     quality_->put(id1, range);
     BadStripList.clear();
   }
   //Fill all the bad components now
   quality_->fillBadComponents();
 }

◆ makeSummary()

void SiStripHitEffFromCalibTree::makeSummary ( )

private

Definition at line 1187 of file SiStripHitEffFromCalibTree.cc.

References _effPlotMin, _showEndcapSides, _showOnlyGoodModules, _showRings, _showTOB6TEC9, _title, python.cmstools::all(), alllayerfound, alllayertotal, MillePedeFileConverter_cfg::e, newFWLiteAna::found, fs, getLayerSideName(), goodlayerfound, goodlayertotal, mps_fire::i, dqmiolumiharvest::j, dqmdumpme::k, label, phase1PixelTopology::layerName, LOGPRINT, TFileService::make(), MuonTCMETValueMapProducer_cff::nLayers, nTEClayers, compare::overlay(), and SiStripLayers.

Referenced by algoAnalyze().

                                              {
   //setTDRStyle();
 
   int nLayers = 34;
   if (_showRings)
     nLayers = 30;
   if (!_showEndcapSides) {
     if (!_showRings)
       nLayers = SiStripLayers;
     else
       nLayers = 20;
   }
 
   TH1F* found = fs->make<TH1F>("found", "found", nLayers + 1, 0, nLayers + 1);
   TH1F* all = fs->make<TH1F>("all", "all", nLayers + 1, 0, nLayers + 1);
   TH1F* found2 = fs->make<TH1F>("found2", "found2", nLayers + 1, 0, nLayers + 1);
   TH1F* all2 = fs->make<TH1F>("all2", "all2", nLayers + 1, 0, nLayers + 1);
   // first bin only to keep real data off the y axis so set to -1
   found->SetBinContent(0, -1);
   all->SetBinContent(0, 1);
 
   // new ROOT version: TGraph::Divide don't handle null or negative values
   for (Long_t i = 1; i < nLayers + 2; ++i) {
     found->SetBinContent(i, 1e-6);
     all->SetBinContent(i, 1);
     found2->SetBinContent(i, 1e-6);
     all2->SetBinContent(i, 1);
   }
 
   TCanvas* c7 = new TCanvas("c7", " test ", 10, 10, 800, 600);
   c7->SetFillColor(0);
   c7->SetGrid();
 
   int nLayers_max = nLayers + 1;  // barrel+endcap
   if (!_showEndcapSides)
     nLayers_max = 11;  // barrel
   for (Long_t i = 1; i < nLayers_max; ++i) {
     LOGPRINT << "Fill only good modules layer " << i << ":  S = " << goodlayerfound[i]
              << "    B = " << goodlayertotal[i];
     if (goodlayertotal[i] > 5) {
       found->SetBinContent(i, goodlayerfound[i]);
       all->SetBinContent(i, goodlayertotal[i]);
     }
 
     LOGPRINT << "Filling all modules layer " << i << ":  S = " << alllayerfound[i] << "    B = " << alllayertotal[i];
     if (alllayertotal[i] > 5) {
       found2->SetBinContent(i, alllayerfound[i]);
       all2->SetBinContent(i, alllayertotal[i]);
     }
   }
 
   // endcap - merging sides
   if (!_showEndcapSides) {
     for (Long_t i = 11; i < 14; ++i) {  // TID disks
       LOGPRINT << "Fill only good modules layer " << i << ":  S = " << goodlayerfound[i] + goodlayerfound[i + 3]
                << "    B = " << goodlayertotal[i] + goodlayertotal[i + 3];
       if (goodlayertotal[i] + goodlayertotal[i + 3] > 5) {
         found->SetBinContent(i, goodlayerfound[i] + goodlayerfound[i + 3]);
         all->SetBinContent(i, goodlayertotal[i] + goodlayertotal[i + 3]);
       }
       LOGPRINT << "Filling all modules layer " << i << ":  S = " << alllayerfound[i] + alllayerfound[i + 3]
                << "    B = " << alllayertotal[i] + alllayertotal[i + 3];
       if (alllayertotal[i] + alllayertotal[i + 3] > 5) {
         found2->SetBinContent(i, alllayerfound[i] + alllayerfound[i + 3]);
         all2->SetBinContent(i, alllayertotal[i] + alllayertotal[i + 3]);
       }
     }
     for (Long_t i = 17; i < 17 + nTEClayers; ++i) {  // TEC disks
       LOGPRINT << "Fill only good modules layer " << i - 3
                << ":  S = " << goodlayerfound[i] + goodlayerfound[i + nTEClayers]
                << "    B = " << goodlayertotal[i] + goodlayertotal[i + nTEClayers];
       if (goodlayertotal[i] + goodlayertotal[i + nTEClayers] > 5) {
         found->SetBinContent(i - 3, goodlayerfound[i] + goodlayerfound[i + nTEClayers]);
         all->SetBinContent(i - 3, goodlayertotal[i] + goodlayertotal[i + nTEClayers]);
       }
       LOGPRINT << "Filling all modules layer " << i - 3 << ":  S = " << alllayerfound[i] + alllayerfound[i + nTEClayers]
                << "    B = " << alllayertotal[i] + alllayertotal[i + nTEClayers];
       if (alllayertotal[i] + alllayertotal[i + nTEClayers] > 5) {
         found2->SetBinContent(i - 3, alllayerfound[i] + alllayerfound[i + nTEClayers]);
         all2->SetBinContent(i - 3, alllayertotal[i] + alllayertotal[i + nTEClayers]);
       }
     }
   }
 
   found->Sumw2();
   all->Sumw2();
 
   found2->Sumw2();
   all2->Sumw2();
 
   TGraphAsymmErrors* gr = fs->make<TGraphAsymmErrors>(nLayers + 1);
   gr->SetName("eff_good");
   gr->BayesDivide(found, all);
 
   TGraphAsymmErrors* gr2 = fs->make<TGraphAsymmErrors>(nLayers + 1);
   gr2->SetName("eff_all");
   gr2->BayesDivide(found2, all2);
 
   for (int j = 0; j < nLayers + 1; j++) {
     gr->SetPointError(j, 0., 0., gr->GetErrorYlow(j), gr->GetErrorYhigh(j));
     gr2->SetPointError(j, 0., 0., gr2->GetErrorYlow(j), gr2->GetErrorYhigh(j));
   }
 
   gr->GetXaxis()->SetLimits(0, nLayers);
   gr->SetMarkerColor(2);
   gr->SetMarkerSize(1.2);
   gr->SetLineColor(2);
   gr->SetLineWidth(4);
   gr->SetMarkerStyle(20);
   gr->SetMinimum(_effPlotMin);
   gr->SetMaximum(1.001);
   gr->GetYaxis()->SetTitle("Efficiency");
   gStyle->SetTitleFillColor(0);
   gStyle->SetTitleBorderSize(0);
   gr->SetTitle(_title);
 
   gr2->GetXaxis()->SetLimits(0, nLayers);
   gr2->SetMarkerColor(1);
   gr2->SetMarkerSize(1.2);
   gr2->SetLineColor(1);
   gr2->SetLineWidth(4);
   gr2->SetMarkerStyle(21);
   gr2->SetMinimum(_effPlotMin);
   gr2->SetMaximum(1.001);
   gr2->GetYaxis()->SetTitle("Efficiency");
   gr2->SetTitle(_title);
 
   for (Long_t k = 1; k < nLayers + 1; k++) {
     TString label;
     if (_showEndcapSides)
       label = getLayerSideName(k);
     else
       label = ::layerName(k, _showRings, nTEClayers);
     if (!_showTOB6TEC9) {
       if (k == 10)
         label = "";
       if (!_showRings && k == nLayers)
         label = "";
       if (!_showRings && _showEndcapSides && k == 25)
         label = "";
     }
     if (!_showRings) {
       if (_showEndcapSides) {
         gr->GetXaxis()->SetBinLabel(((k + 1) * 100 + 2) / (nLayers)-4, label);
         gr2->GetXaxis()->SetBinLabel(((k + 1) * 100 + 2) / (nLayers)-4, label);
       } else {
         gr->GetXaxis()->SetBinLabel((k + 1) * 100 / (nLayers)-6, label);
         gr2->GetXaxis()->SetBinLabel((k + 1) * 100 / (nLayers)-6, label);
       }
     } else {
       if (_showEndcapSides) {
         gr->GetXaxis()->SetBinLabel((k + 1) * 100 / (nLayers)-4, label);
         gr2->GetXaxis()->SetBinLabel((k + 1) * 100 / (nLayers)-4, label);
       } else {
         gr->GetXaxis()->SetBinLabel((k + 1) * 100 / (nLayers)-7, label);
         gr2->GetXaxis()->SetBinLabel((k + 1) * 100 / (nLayers)-7, label);
       }
     }
   }
 
   gr->Draw("AP");
   gr->GetXaxis()->SetNdivisions(36);
 
   c7->cd();
   TPad* overlay = new TPad("overlay", "", 0, 0, 1, 1);
   overlay->SetFillStyle(4000);
   overlay->SetFillColor(0);
   overlay->SetFrameFillStyle(4000);
   overlay->Draw("same");
   overlay->cd();
   if (!_showOnlyGoodModules)
     gr2->Draw("AP");
 
   TLegend* leg = new TLegend(0.70, 0.27, 0.88, 0.40);
   leg->AddEntry(gr, "Good Modules", "p");
   if (!_showOnlyGoodModules)
     leg->AddEntry(gr2, "All Modules", "p");
   leg->SetTextSize(0.020);
   leg->SetFillColor(0);
   leg->Draw("same");
 
   c7->SaveAs("Summary.png");
 }

◆ makeSummaryVsBx()

void SiStripHitEffFromCalibTree::makeSummaryVsBx ( )

private

Definition at line 1371 of file SiStripHitEffFromCalibTree.cc.

References _showRings, _spaceBetweenTrains, MillePedeFileConverter_cfg::e, dqmMemoryStats::float, dqm-mbProfile::format, newFWLiteAna::found, fs, layerfound_perBx, phase1PixelTopology::layerName, layertotal_perBx, LOGPRINT, LaserClient_cfi::low, TFileService::make(), MuonTCMETValueMapProducer_cff::nLayers, nTEClayers, SiStripLayers, dqmMemoryStats::total, and up.

Referenced by algoAnalyze().

                                                  {
   LOGPRINT << "Computing efficiency vs bx";
 
   unsigned int nLayers = SiStripLayers;
   if (_showRings)
     nLayers = 20;
 
   for (unsigned int ilayer = 1; ilayer < nLayers; ilayer++) {
     TH1F* hfound = fs->make<TH1F>(Form("foundVsBx_layer%i", ilayer), Form("layer %i", ilayer), 3565, 0, 3565);
     TH1F* htotal = fs->make<TH1F>(Form("totalVsBx_layer%i", ilayer), Form("layer %i", ilayer), 3565, 0, 3565);
 
     for (unsigned int ibx = 0; ibx < 3566; ibx++) {
       hfound->SetBinContent(ibx, 1e-6);
       htotal->SetBinContent(ibx, 1);
     }
     map<unsigned int, vector<int> >::iterator iterMapvsBx;
     for (iterMapvsBx = layerfound_perBx.begin(); iterMapvsBx != layerfound_perBx.end(); ++iterMapvsBx)
       hfound->SetBinContent(iterMapvsBx->first, iterMapvsBx->second[ilayer]);
     for (iterMapvsBx = layertotal_perBx.begin(); iterMapvsBx != layertotal_perBx.end(); ++iterMapvsBx)
       if (iterMapvsBx->second[ilayer] > 0)
         htotal->SetBinContent(iterMapvsBx->first, iterMapvsBx->second[ilayer]);
 
     hfound->Sumw2();
     htotal->Sumw2();
 
     TGraphAsymmErrors* geff = fs->make<TGraphAsymmErrors>(3564);
     geff->SetName(Form("effVsBx_layer%i", ilayer));
     geff->SetTitle(fmt::format("Hit Efficiency vs bx - {}", ::layerName(ilayer, _showRings, nTEClayers)).c_str());
     geff->BayesDivide(hfound, htotal);
 
     //Average over trains
     TGraphAsymmErrors* geff_avg = fs->make<TGraphAsymmErrors>();
     geff_avg->SetName(Form("effVsBxAvg_layer%i", ilayer));
     geff_avg->SetTitle(fmt::format("Hit Efficiency vs bx - {}", ::layerName(ilayer, _showRings, nTEClayers)).c_str());
     geff_avg->SetMarkerStyle(20);
     int ibx = 0;
     int previous_bx = -80;
     int delta_bx = 0;
     int nbx = 0;
     int found = 0;
     int total = 0;
     double sum_bx = 0;
     int ipt = 0;
     float low, up, eff;
     int firstbx = 0;
     for (iterMapvsBx = layertotal_perBx.begin(); iterMapvsBx != layertotal_perBx.end(); ++iterMapvsBx) {
       ibx = iterMapvsBx->first;
       delta_bx = ibx - previous_bx;
       // consider a new train
       if (delta_bx > (int)_spaceBetweenTrains && nbx > 0 && total > 0) {
         eff = found / (float)total;
         //LOGPRINT<<"new train "<<ipt<<" "<<sum_bx/nbx<<" "<<eff<<endl;
         geff_avg->SetPoint(ipt, sum_bx / nbx, eff);
         low = TEfficiency::Bayesian(total, found, .683, 1, 1, false);
         up = TEfficiency::Bayesian(total, found, .683, 1, 1, true);
         geff_avg->SetPointError(ipt, sum_bx / nbx - firstbx, previous_bx - sum_bx / nbx, eff - low, up - eff);
         ipt++;
         sum_bx = 0;
         found = 0;
         total = 0;
         nbx = 0;
         firstbx = ibx;
       }
       sum_bx += ibx;
       found += hfound->GetBinContent(ibx);
       total += htotal->GetBinContent(ibx);
       nbx++;
 
       previous_bx = ibx;
     }
     //last train
     eff = found / (float)total;
     //LOGPRINT<<"new train "<<ipt<<" "<<sum_bx/nbx<<" "<<eff<<endl;
     geff_avg->SetPoint(ipt, sum_bx / nbx, eff);
     low = TEfficiency::Bayesian(total, found, .683, 1, 1, false);
     up = TEfficiency::Bayesian(total, found, .683, 1, 1, true);
     geff_avg->SetPointError(ipt, sum_bx / nbx - firstbx, previous_bx - sum_bx / nbx, eff - low, up - eff);
   }
 }

◆ makeSummaryVsCM()

void SiStripHitEffFromCalibTree::makeSummaryVsCM ( )

private

Definition at line 1527 of file SiStripHitEffFromCalibTree.cc.

References computeEff(), layerfound_vsCM, layertotal_vsCM, and LOGPRINT.

Referenced by algoAnalyze().

                                                  {
   LOGPRINT << "Computing efficiency vs CM";
   computeEff(layerfound_vsCM, layertotal_vsCM, "effVsCM");
 }

◆ makeSummaryVsLumi()

void SiStripHitEffFromCalibTree::makeSummaryVsLumi ( )

private

Definition at line 1489 of file SiStripHitEffFromCalibTree.cc.

References _showRings, computeEff(), instLumiHisto, layerfound_vsLumi, layerfound_vsPU, layertotal_vsLumi, layertotal_vsPU, LOGPRINT, MuonTCMETValueMapProducer_cff::nLayers, PUHisto, and SiStripLayers.

Referenced by algoAnalyze().

                                                    {
   LOGPRINT << "Computing efficiency vs lumi";
 
   if (instLumiHisto->GetEntries())  // from infos per event
     LOGPRINT << "Avg conditions (avg+/-rms):   lumi :" << instLumiHisto->GetMean() << "+/-" << instLumiHisto->GetRMS()
              << "   pu: " << PUHisto->GetMean() << "+/-" << PUHisto->GetRMS();
 
   else {  // from infos per hit
 
     unsigned int nLayers = SiStripLayers;
     if (_showRings)
       nLayers = 20;
     unsigned int nLayersForAvg = 0;
     float layerLumi = 0;
     float layerPU = 0;
     float avgLumi = 0;
     float avgPU = 0;
 
     LOGPRINT << "Lumi summary:  (avg over trajectory measurements)";
     for (unsigned int ilayer = 1; ilayer < nLayers; ilayer++) {
       layerLumi = layertotal_vsLumi[ilayer]->GetMean();
       layerPU = layertotal_vsPU[ilayer]->GetMean();
       //LOGPRINT<<" layer "<<ilayer<<"  lumi: "<<layerLumi<<"  pu: "<<layerPU<<endl;
       if (layerLumi != 0 && layerPU != 0) {
         avgLumi += layerLumi;
         avgPU += layerPU;
         nLayersForAvg++;
       }
     }
     avgLumi /= nLayersForAvg;
     avgPU /= nLayersForAvg;
     LOGPRINT << "Avg conditions:   lumi :" << avgLumi << "  pu: " << avgPU;
   }
 
   computeEff(layerfound_vsLumi, layertotal_vsLumi, "effVsLumi");
   computeEff(layerfound_vsPU, layertotal_vsPU, "effVsPU");
 }

◆ makeTKMap()

void SiStripHitEffFromCalibTree::makeTKMap ( bool autoTagging = false )

private

Definition at line 1000 of file SiStripHitEffFromCalibTree.cc.

References _showRings, _title, _tkMapMin, BadModules, TrackerMap::fill(), TrackerMap::fillc(), fs, mps_fire::i, layerfound, phase1PixelTopology::layerName, layertotal, LOGPRINT, TFileService::make(), modCounter, nModsMin, nTEClayers, TrackerMap::save(), SiStripLayers, threshold, tkmap, tkmapbad, tkmapden, tkmapeff, and tkmapnum.

Referenced by algoAnalyze().

                                                                    {
   LOGPRINT << "Entering TKMap generation!\n";
   tkmap = new TrackerMap("  Detector Inefficiency  ");
   tkmapbad = new TrackerMap("  Inefficient Modules  ");
   tkmapeff = new TrackerMap(_title.Data());
   tkmapnum = new TrackerMap(" Detector numerator   ");
   tkmapden = new TrackerMap(" Detector denominator ");
 
   double myeff, mynum, myden, myeff_up;
   double layer_min_eff = 0;
 
   for (Long_t i = 1; i <= SiStripLayers; i++) {
     //Loop over every layer, extracting the information from
     //the map of the efficiencies
     layertotal[i] = 0;
     layerfound[i] = 0;
     TH1F* hEffInLayer =
         fs->make<TH1F>(Form("eff_layer%i", int(i)), Form("Module efficiency in layer %i", int(i)), 201, 0, 1.005);
 
     map<unsigned int, pair<unsigned int, unsigned int> >::const_iterator ih;
     for (ih = modCounter[i].begin(); ih != modCounter[i].end(); ih++) {
       //We should be in the layer in question, and looping over all of the modules in said layer
       //Generate the list for the TKmap, and the bad module list
       mynum = (double)(((*ih).second).second);
       myden = (double)(((*ih).second).first);
       if (myden > 0)
         myeff = mynum / myden;
       else
         myeff = 0;
       hEffInLayer->Fill(myeff);
 
       if (!autoTagging) {
         if ((myden >= nModsMin) && (myeff < threshold)) {
           //We have a bad module, put it in the list!
           BadModules[(*ih).first] = myeff;
           tkmapbad->fillc((*ih).first, 255, 0, 0);
           LOGPRINT << "Layer " << i << " (" << ::layerName(i, _showRings, nTEClayers) << ")  module " << (*ih).first
                    << " efficiency: " << myeff << " , " << mynum << "/" << myden;
         } else {
           //Fill the bad list with empty results for every module
           tkmapbad->fillc((*ih).first, 255, 255, 255);
         }
         if (myeff < threshold)
           LOGPRINT << "Layer " << i << " (" << ::layerName(i, _showRings, nTEClayers) << ")  module " << (*ih).first
                    << " efficiency: " << myeff << " , " << mynum << "/" << myden;
         if (myden < nModsMin) {
           LOGPRINT << "Layer " << i << " (" << ::layerName(i, _showRings, nTEClayers) << ")  module " << (*ih).first
                    << " is under occupancy at " << myden;
         }
       }
 
       //Put any module into the TKMap
       tkmap->fill((*ih).first, 1. - myeff);
       tkmapeff->fill((*ih).first, myeff);
       tkmapnum->fill((*ih).first, mynum);
       tkmapden->fill((*ih).first, myden);
 
       //Add the number of hits in the layer
       layertotal[i] += long(myden);
       layerfound[i] += long(mynum);
     }
 
     if (autoTagging) {
       //Compute threshold to use for each layer
       hEffInLayer->GetXaxis()->SetRange(3, hEffInLayer->GetNbinsX() + 1);  // Remove from the avg modules below 1%
       layer_min_eff =
           hEffInLayer->GetMean() - 2.5 * hEffInLayer->GetRMS();  // uses RMS in case the distribution is wide
       if (threshold > 2.5 * hEffInLayer->GetRMS())
         layer_min_eff = hEffInLayer->GetMean() - threshold;  // otherwise uses the parameter 'threshold'
       LOGPRINT << "Layer " << i << " threshold for bad modules: <" << layer_min_eff
                << "  (layer mean: " << hEffInLayer->GetMean() << " rms: " << hEffInLayer->GetRMS() << ")";
 
       hEffInLayer->GetXaxis()->SetRange(1, hEffInLayer->GetNbinsX() + 1);
 
       for (ih = modCounter[i].begin(); ih != modCounter[i].end(); ih++) {
         // Second loop over modules to tag inefficient ones
         mynum = (double)(((*ih).second).second);
         myden = (double)(((*ih).second).first);
         if (myden > 0)
           myeff = mynum / myden;
         else
           myeff = 0;
         // upper limit on the efficiency
         myeff_up = TEfficiency::Bayesian(myden, mynum, .99, 1, 1, true);
         if ((myden >= nModsMin) && (myeff_up < layer_min_eff)) {
           //We have a bad module, put it in the list!
           BadModules[(*ih).first] = myeff;
           tkmapbad->fillc((*ih).first, 255, 0, 0);
         } else {
           //Fill the bad list with empty results for every module
           tkmapbad->fillc((*ih).first, 255, 255, 255);
         }
         if (myeff_up < layer_min_eff + 0.08)  // printing message also for modules slighly above (8%) the limit
           LOGPRINT << "Layer " << i << " (" << ::layerName(i, _showRings, nTEClayers) << ")  module " << (*ih).first
                    << " efficiency: " << myeff << " , " << mynum << "/" << myden << " , upper limit: " << myeff_up;
         if (myden < nModsMin) {
           LOGPRINT << "Layer " << i << " (" << ::layerName(i, _showRings, nTEClayers) << ")  module " << (*ih).first
                    << " layer " << i << " is under occupancy at " << myden;
         }
       }
     }
   }
   tkmap->save(true, 0, 0, "SiStripHitEffTKMap.png");
   tkmapbad->save(true, 0, 0, "SiStripHitEffTKMapBad.png");
   tkmapeff->save(true, _tkMapMin, 1., "SiStripHitEffTKMapEff.png");
   tkmapnum->save(true, 0, 0, "SiStripHitEffTKMapNum.png");
   tkmapden->save(true, 0, 0, "SiStripHitEffTKMapDen.png");
   LOGPRINT << "Finished TKMap Generation\n";
 }

◆ setBadComponents()

void SiStripHitEffFromCalibTree::setBadComponents	(	int	i,
		int	component,
		SiStripQuality::BadComponent &	BC,
		std::stringstream	ssV[4][19],
		int	NBadComponent[4][19][4]
	)

private

Definition at line 1573 of file SiStripHitEffFromCalibTree.cc.

References _detInfo, SiStripQuality::BadComponent::BadApvs, SiStripQuality::BadComponent::BadFibers, SiStripQuality::BadComponent::BadModule, SiStripQuality::BadComponent::detid, SiStripDetInfo::getNumberOfApvsAndStripLength(), mps_fire::i, and testProducerWithPsetDescEmpty_cfi::x1.

Referenced by algoAnalyze().

                                                                                                                    {
   int napv = _detInfo.getNumberOfApvsAndStripLength(BC.detid).first;
 
   ssV[i][component] << "\n\t\t " << BC.detid << " \t " << BC.BadModule << " \t " << ((BC.BadFibers) & 0x1) << " ";
   if (napv == 4)
     ssV[i][component] << "x " << ((BC.BadFibers >> 1) & 0x1);
 
   if (napv == 6)
     ssV[i][component] << ((BC.BadFibers >> 1) & 0x1) << " " << ((BC.BadFibers >> 2) & 0x1);
   ssV[i][component] << " \t " << ((BC.BadApvs) & 0x1) << " " << ((BC.BadApvs >> 1) & 0x1) << " ";
   if (napv == 4)
     ssV[i][component] << "x x " << ((BC.BadApvs >> 2) & 0x1) << " " << ((BC.BadApvs >> 3) & 0x1);
   if (napv == 6)
     ssV[i][component] << ((BC.BadApvs >> 2) & 0x1) << " " << ((BC.BadApvs >> 3) & 0x1) << " "
                       << ((BC.BadApvs >> 4) & 0x1) << " " << ((BC.BadApvs >> 5) & 0x1) << " ";
 
   if (BC.BadApvs) {
     NBadComponent[i][0][2] += ((BC.BadApvs >> 5) & 0x1) + ((BC.BadApvs >> 4) & 0x1) + ((BC.BadApvs >> 3) & 0x1) +
                               ((BC.BadApvs >> 2) & 0x1) + ((BC.BadApvs >> 1) & 0x1) + ((BC.BadApvs) & 0x1);
     NBadComponent[i][component][2] += ((BC.BadApvs >> 5) & 0x1) + ((BC.BadApvs >> 4) & 0x1) +
                                       ((BC.BadApvs >> 3) & 0x1) + ((BC.BadApvs >> 2) & 0x1) +
                                       ((BC.BadApvs >> 1) & 0x1) + ((BC.BadApvs) & 0x1);
   }
   if (BC.BadFibers) {
     NBadComponent[i][0][1] += ((BC.BadFibers >> 2) & 0x1) + ((BC.BadFibers >> 1) & 0x1) + ((BC.BadFibers) & 0x1);
     NBadComponent[i][component][1] +=
         ((BC.BadFibers >> 2) & 0x1) + ((BC.BadFibers >> 1) & 0x1) + ((BC.BadFibers) & 0x1);
   }
   if (BC.BadModule) {
     NBadComponent[i][0][0]++;
     NBadComponent[i][component][0]++;
   }
 }

◆ totalStatistics()

void SiStripHitEffFromCalibTree::totalStatistics ( )

private

Definition at line 1139 of file SiStripHitEffFromCalibTree.cc.

References _showRings, mps_fire::i, layerfound, phase1PixelTopology::layerName, layertotal, LOGPRINT, nTEClayers, and SiStripLayers.

Referenced by algoAnalyze().

                                                  {
   //Calculate the statistics by layer
   int totalfound = 0;
   int totaltotal = 0;
   double layereff;
   int subdetfound[5];
   int subdettotal[5];
 
   for (Long_t i = 1; i < 5; i++) {
     subdetfound[i] = 0;
     subdettotal[i] = 0;
   }
 
   for (Long_t i = 1; i <= SiStripLayers; i++) {
     layereff = double(layerfound[i]) / double(layertotal[i]);
     LOGPRINT << "Layer " << i << " (" << ::layerName(i, _showRings, nTEClayers) << ") has total efficiency " << layereff
              << " " << layerfound[i] << "/" << layertotal[i];
     totalfound += layerfound[i];
     totaltotal += layertotal[i];
     if (i < 5) {
       subdetfound[1] += layerfound[i];
       subdettotal[1] += layertotal[i];
     }
     if (i >= 5 && i < 11) {
       subdetfound[2] += layerfound[i];
       subdettotal[2] += layertotal[i];
     }
     if (i >= 11 && i < 14) {
       subdetfound[3] += layerfound[i];
       subdettotal[3] += layertotal[i];
     }
     if (i >= 14) {
       subdetfound[4] += layerfound[i];
       subdettotal[4] += layertotal[i];
     }
   }
 
   LOGPRINT << "The total efficiency is " << double(totalfound) / double(totaltotal);
   LOGPRINT << "      TIB: " << double(subdetfound[1]) / subdettotal[1] << " " << subdetfound[1] << "/"
            << subdettotal[1];
   LOGPRINT << "      TOB: " << double(subdetfound[2]) / subdettotal[2] << " " << subdetfound[2] << "/"
            << subdettotal[2];
   LOGPRINT << "      TID: " << double(subdetfound[3]) / subdettotal[3] << " " << subdetfound[3] << "/"
            << subdettotal[3];
   LOGPRINT << "      TEC: " << double(subdetfound[4]) / subdettotal[4] << " " << subdetfound[4] << "/"
            << subdettotal[4];
 }

Member Data Documentation

◆ _autoIneffModTagging

bool SiStripHitEffFromCalibTree::_autoIneffModTagging

private

Definition at line 135 of file SiStripHitEffFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripHitEffFromCalibTree().

◆ _badModulesFile

string SiStripHitEffFromCalibTree::_badModulesFile

private

Definition at line 134 of file SiStripHitEffFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripHitEffFromCalibTree().

◆ _bunchx

unsigned int SiStripHitEffFromCalibTree::_bunchx

private

Definition at line 141 of file SiStripHitEffFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripHitEffFromCalibTree().

◆ _clusterMatchingMethod

unsigned int SiStripHitEffFromCalibTree::_clusterMatchingMethod

private

Definition at line 136 of file SiStripHitEffFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripHitEffFromCalibTree().

◆ _clusterTrajDist

float SiStripHitEffFromCalibTree::_clusterTrajDist

private

Definition at line 138 of file SiStripHitEffFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripHitEffFromCalibTree().

◆ _detInfo

SiStripDetInfo SiStripHitEffFromCalibTree::_detInfo

private

Definition at line 124 of file SiStripHitEffFromCalibTree.cc.

Referenced by algoAnalyze(), makeSQLite(), setBadComponents(), and SiStripHitEffFromCalibTree().

◆ _effPlotMin

float SiStripHitEffFromCalibTree::_effPlotMin

private

Definition at line 149 of file SiStripHitEffFromCalibTree.cc.

Referenced by makeSummary(), and SiStripHitEffFromCalibTree().

◆ _ResXSig

float SiStripHitEffFromCalibTree::_ResXSig

private

Definition at line 137 of file SiStripHitEffFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripHitEffFromCalibTree().

◆ _showEndcapSides

bool SiStripHitEffFromCalibTree::_showEndcapSides

private

Definition at line 144 of file SiStripHitEffFromCalibTree.cc.

Referenced by makeSummary(), and SiStripHitEffFromCalibTree().

◆ _showOnlyGoodModules

bool SiStripHitEffFromCalibTree::_showOnlyGoodModules

private

Definition at line 147 of file SiStripHitEffFromCalibTree.cc.

Referenced by makeSummary(), and SiStripHitEffFromCalibTree().

◆ _showRings

bool SiStripHitEffFromCalibTree::_showRings

private

Definition at line 145 of file SiStripHitEffFromCalibTree.cc.

Referenced by algoAnalyze(), computeEff(), getLayerSideName(), makeSummary(), makeSummaryVsBx(), makeSummaryVsLumi(), makeTKMap(), SiStripHitEffFromCalibTree(), and totalStatistics().

◆ _showTOB6TEC9

bool SiStripHitEffFromCalibTree::_showTOB6TEC9

private

Definition at line 146 of file SiStripHitEffFromCalibTree.cc.

Referenced by algoAnalyze(), makeSummary(), and SiStripHitEffFromCalibTree().

◆ _spaceBetweenTrains

unsigned int SiStripHitEffFromCalibTree::_spaceBetweenTrains

private

Definition at line 142 of file SiStripHitEffFromCalibTree.cc.

Referenced by makeSummaryVsBx(), and SiStripHitEffFromCalibTree().

◆ _stripsApvEdge

float SiStripHitEffFromCalibTree::_stripsApvEdge

private

Definition at line 139 of file SiStripHitEffFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripHitEffFromCalibTree().

◆ _title

TString SiStripHitEffFromCalibTree::_title

private

Definition at line 150 of file SiStripHitEffFromCalibTree.cc.

Referenced by html.PageSet::_getPage(), plotting.Subtract::create(), plotting.Transform::create(), plotting.FakeDuplicate::create(), plotting.CutEfficiency::create(), plotting.Plot::create(), makeSummary(), makeTKMap(), and SiStripHitEffFromCalibTree().

◆ _tkGeomToken

edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> SiStripHitEffFromCalibTree::_tkGeomToken

private

Definition at line 152 of file SiStripHitEffFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripHitEffFromCalibTree().

◆ _tkMapMin

float SiStripHitEffFromCalibTree::_tkMapMin

private

Definition at line 148 of file SiStripHitEffFromCalibTree.cc.

Referenced by makeTKMap(), and SiStripHitEffFromCalibTree().

◆ _tTopoToken

edm::ESGetToken<TrackerTopology, TrackerTopologyRcd> SiStripHitEffFromCalibTree::_tTopoToken

private

Definition at line 153 of file SiStripHitEffFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripHitEffFromCalibTree().

◆ _useCM

bool SiStripHitEffFromCalibTree::_useCM

private

Definition at line 143 of file SiStripHitEffFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripHitEffFromCalibTree().

◆ _useOnlyHighPurityTracks

bool SiStripHitEffFromCalibTree::_useOnlyHighPurityTracks

private

Definition at line 140 of file SiStripHitEffFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripHitEffFromCalibTree().

◆ alllayerfound

int SiStripHitEffFromCalibTree::alllayerfound[35]

private

Definition at line 185 of file SiStripHitEffFromCalibTree.cc.

Referenced by algoAnalyze(), and makeSummary().

◆ alllayertotal

int SiStripHitEffFromCalibTree::alllayertotal[35]

private

Definition at line 184 of file SiStripHitEffFromCalibTree.cc.

Referenced by algoAnalyze(), and makeSummary().

◆ BadModules

map<unsigned int, double> SiStripHitEffFromCalibTree::BadModules

private

Definition at line 186 of file SiStripHitEffFromCalibTree.cc.

Referenced by makeSQLite(), and makeTKMap().

◆ bxHisto

TH1F* SiStripHitEffFromCalibTree::bxHisto

private

Definition at line 157 of file SiStripHitEffFromCalibTree.cc.

Referenced by algoAnalyze().

◆ CalibTree

TTree* SiStripHitEffFromCalibTree::CalibTree

private

Definition at line 129 of file SiStripHitEffFromCalibTree.cc.

Referenced by algoAnalyze().

◆ CalibTreeFilenames

vector<string> SiStripHitEffFromCalibTree::CalibTreeFilenames

private

Definition at line 130 of file SiStripHitEffFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripHitEffFromCalibTree().

◆ doSummary

unsigned int SiStripHitEffFromCalibTree::doSummary

private

Definition at line 133 of file SiStripHitEffFromCalibTree.cc.

Referenced by SiStripHitEffFromCalibTree().

◆ eventInfos

map<pair<unsigned int, unsigned int>, array<double, 3> > SiStripHitEffFromCalibTree::eventInfos

private

Definition at line 162 of file SiStripHitEffFromCalibTree.cc.

Referenced by algoAnalyze().

◆ FileInPath_

edm::FileInPath SiStripHitEffFromCalibTree::FileInPath_

private

Definition at line 125 of file SiStripHitEffFromCalibTree.cc.

Referenced by SiStripHitEffFromCalibTree().

◆ fs

edm::Service<TFileService> SiStripHitEffFromCalibTree::fs

private

Definition at line 123 of file SiStripHitEffFromCalibTree.cc.

Referenced by algoAnalyze(), computeEff(), makeHotColdMaps(), makeSummary(), makeSummaryVsBx(), and makeTKMap().

◆ goodlayerfound

int SiStripHitEffFromCalibTree::goodlayerfound[35]

private

Definition at line 183 of file SiStripHitEffFromCalibTree.cc.

Referenced by algoAnalyze(), and makeSummary().

◆ goodlayertotal

int SiStripHitEffFromCalibTree::goodlayertotal[35]

private

Definition at line 182 of file SiStripHitEffFromCalibTree.cc.

Referenced by algoAnalyze(), and makeSummary().

◆ hits

vector<hit> SiStripHitEffFromCalibTree::hits[23]

private

Definition at line 164 of file SiStripHitEffFromCalibTree.cc.

Referenced by algoAnalyze(), and makeHotColdMaps().

◆ HotColdMaps

vector<TH2F*> SiStripHitEffFromCalibTree::HotColdMaps

private

Definition at line 165 of file SiStripHitEffFromCalibTree.cc.

Referenced by makeHotColdMaps().

◆ instLumiHisto

TH1F* SiStripHitEffFromCalibTree::instLumiHisto

private

Definition at line 158 of file SiStripHitEffFromCalibTree.cc.

Referenced by algoAnalyze(), and makeSummaryVsLumi().

◆ layerfound

long SiStripHitEffFromCalibTree::layerfound[23]

private

Definition at line 172 of file SiStripHitEffFromCalibTree.cc.

Referenced by algoAnalyze(), makeTKMap(), and totalStatistics().

◆ layerfound_perBx

map<unsigned int, vector<int> > SiStripHitEffFromCalibTree::layerfound_perBx

private

Definition at line 174 of file SiStripHitEffFromCalibTree.cc.

Referenced by algoAnalyze(), and makeSummaryVsBx().

◆ layerfound_vsCM

vector<TH1F*> SiStripHitEffFromCalibTree::layerfound_vsCM

private

Definition at line 180 of file SiStripHitEffFromCalibTree.cc.

Referenced by algoAnalyze(), and makeSummaryVsCM().

◆ layerfound_vsLumi

vector<TH1F*> SiStripHitEffFromCalibTree::layerfound_vsLumi

private

Definition at line 176 of file SiStripHitEffFromCalibTree.cc.

Referenced by algoAnalyze(), and makeSummaryVsLumi().

◆ layerfound_vsPU

vector<TH1F*> SiStripHitEffFromCalibTree::layerfound_vsPU

private

Definition at line 178 of file SiStripHitEffFromCalibTree.cc.

Referenced by algoAnalyze(), and makeSummaryVsLumi().

◆ layertotal

long SiStripHitEffFromCalibTree::layertotal[23]

private

Definition at line 173 of file SiStripHitEffFromCalibTree.cc.

Referenced by algoAnalyze(), makeTKMap(), and totalStatistics().

◆ layertotal_perBx

map<unsigned int, vector<int> > SiStripHitEffFromCalibTree::layertotal_perBx

private

Definition at line 175 of file SiStripHitEffFromCalibTree.cc.

Referenced by algoAnalyze(), and makeSummaryVsBx().

◆ layertotal_vsCM

vector<TH1F*> SiStripHitEffFromCalibTree::layertotal_vsCM

private

Definition at line 181 of file SiStripHitEffFromCalibTree.cc.

Referenced by algoAnalyze(), and makeSummaryVsCM().

◆ layertotal_vsLumi

vector<TH1F*> SiStripHitEffFromCalibTree::layertotal_vsLumi

private

Definition at line 177 of file SiStripHitEffFromCalibTree.cc.

Referenced by algoAnalyze(), and makeSummaryVsLumi().

◆ layertotal_vsPU

vector<TH1F*> SiStripHitEffFromCalibTree::layertotal_vsPU

private

Definition at line 179 of file SiStripHitEffFromCalibTree.cc.

Referenced by algoAnalyze(), and makeSummaryVsLumi().

◆ modCounter

map<unsigned int, pair<unsigned int, unsigned int> > SiStripHitEffFromCalibTree::modCounter[23]

private

Definition at line 166 of file SiStripHitEffFromCalibTree.cc.

Referenced by algoAnalyze(), and makeTKMap().

◆ nModsMin

unsigned int SiStripHitEffFromCalibTree::nModsMin

private

Definition at line 132 of file SiStripHitEffFromCalibTree.cc.

Referenced by algoAnalyze(), makeTKMap(), and SiStripHitEffFromCalibTree().

◆ nTEClayers

unsigned int SiStripHitEffFromCalibTree::nTEClayers

private

Definition at line 155 of file SiStripHitEffFromCalibTree.cc.

Referenced by algoAnalyze(), computeEff(), getLayerSideName(), makeSummary(), makeSummaryVsBx(), makeTKMap(), SiStripHitEffFromCalibTree(), and totalStatistics().

◆ PUHisto

TH1F* SiStripHitEffFromCalibTree::PUHisto

private

Definition at line 159 of file SiStripHitEffFromCalibTree.cc.

Referenced by algoAnalyze(), and makeSummaryVsLumi().

◆ quality_

SiStripQuality* SiStripHitEffFromCalibTree::quality_

private

Definition at line 126 of file SiStripHitEffFromCalibTree.cc.

Referenced by algoAnalyze(), getNewObject(), makeSQLite(), and SiStripHitEffFromCalibTree().

◆ SiStripLayers

constexpr int SiStripHitEffFromCalibTree::SiStripLayers = 22

staticprivate

Definition at line 121 of file SiStripHitEffFromCalibTree.cc.

Referenced by algoAnalyze(), computeEff(), makeHotColdMaps(), makeSummary(), makeSummaryVsBx(), makeSummaryVsLumi(), makeTKMap(), and totalStatistics().

◆ threshold

float SiStripHitEffFromCalibTree::threshold

private

Definition at line 131 of file SiStripHitEffFromCalibTree.cc.

Referenced by algoAnalyze(), utils.StatisticalTest::get_status(), makeTKMap(), and SiStripHitEffFromCalibTree().

◆ tkmap

TrackerMap* SiStripHitEffFromCalibTree::tkmap

private

Definition at line 167 of file SiStripHitEffFromCalibTree.cc.

Referenced by makeTKMap().

◆ tkmapbad

TrackerMap* SiStripHitEffFromCalibTree::tkmapbad

private

Definition at line 168 of file SiStripHitEffFromCalibTree.cc.

Referenced by makeTKMap().

◆ tkmapden

TrackerMap* SiStripHitEffFromCalibTree::tkmapden

private

Definition at line 171 of file SiStripHitEffFromCalibTree.cc.

Referenced by makeTKMap().

◆ tkmapeff

TrackerMap* SiStripHitEffFromCalibTree::tkmapeff

private

Definition at line 169 of file SiStripHitEffFromCalibTree.cc.

Referenced by makeTKMap().

◆ tkmapnum

TrackerMap* SiStripHitEffFromCalibTree::tkmapnum

private

Definition at line 170 of file SiStripHitEffFromCalibTree.cc.

Referenced by makeTKMap().

Public Member Functions

Private Member Functions

Private Attributes

Static Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ SiStripHitEffFromCalibTree()

◆ ~SiStripHitEffFromCalibTree()

Member Function Documentation

◆ algoAnalyze()

◆ algoBeginJob()

◆ algoEndJob()

◆ computeEff()

◆ getLayerSideName()

◆ getNewObject()

◆ makeHotColdMaps()

◆ makeSQLite()

◆ makeSummary()

◆ makeSummaryVsBx()

◆ makeSummaryVsCM()

◆ makeSummaryVsLumi()

◆ makeTKMap()

◆ setBadComponents()

◆ totalStatistics()

Member Data Documentation

◆ _autoIneffModTagging

◆ _badModulesFile

◆ _bunchx

◆ _clusterMatchingMethod

◆ _clusterTrajDist

◆ _detInfo

◆ _effPlotMin

◆ _ResXSig

◆ _showEndcapSides

◆ _showOnlyGoodModules

◆ _showRings

◆ _showTOB6TEC9

◆ _spaceBetweenTrains

◆ _stripsApvEdge

◆ _title

◆ _tkGeomToken

◆ _tkMapMin

◆ _tTopoToken

◆ _useCM

◆ _useOnlyHighPurityTracks

◆ alllayerfound

◆ alllayertotal

◆ BadModules

◆ bxHisto

◆ CalibTree

◆ CalibTreeFilenames

◆ doSummary

◆ eventInfos

◆ FileInPath_

◆ fs

◆ goodlayerfound

◆ goodlayertotal

◆ hits

◆ HotColdMaps

◆ instLumiHisto

◆ layerfound

◆ layerfound_perBx

◆ layerfound_vsCM

◆ layerfound_vsLumi

◆ layerfound_vsPU

◆ layertotal

◆ layertotal_perBx

◆ layertotal_vsCM

◆ layertotal_vsLumi

◆ layertotal_vsPU

◆ modCounter

◆ nModsMin

◆ nTEClayers

◆ PUHisto

◆ quality_

◆ SiStripLayers

◆ threshold

◆ tkmap

◆ tkmapbad