#include <L1TOccupancyClient.h>

Inheritance diagram for L1TOccupancyClient:

Public Member Functions
	L1TOccupancyClient (const edm::ParameterSet &ps)
	Constructor. More...

	~L1TOccupancyClient () override
	Destructor. More...

Public Member Functions inherited from DQMEDHarvester
void	accumulate (edm::Event const &ev, edm::EventSetup const &es) final

void	beginJob () override

void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &) final

void	beginRun (edm::Run const &run, edm::EventSetup const &) override

virtual void	dqmAnalyze (DQMStore::IBooker &, DQMStore::IGetter &, edm::Event const &, edm::EventSetup const &)

	DQMEDHarvester (edm::ParameterSet const &iConfig)

	DQMEDHarvester ()

virtual void	dqmEndRun (DQMStore::IBooker &, DQMStore::IGetter &, edm::Run const &, edm::EventSetup const &)

void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) final

void	endLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &es) final

void	endProcessBlockProduce (edm::ProcessBlock &) final

void	endRun (edm::Run const &, edm::EventSetup const &) override

void	endRunProduce (edm::Run &run, edm::EventSetup const &es) final

	~DQMEDHarvester () override=default

Public Member Functions inherited from edm::one::EDProducer< T >
	EDProducer ()=default

	EDProducer (const EDProducer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const final

Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDProducerBase () override

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

std::vector < edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector < edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

	~ProducerBase () noexcept(false) 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)

Protected Member Functions
void	book (DQMStore::IBooker &ibooker, DQMStore::IGetter &igetter)

void	dqmEndJob (DQMStore::IBooker &ibooker, DQMStore::IGetter &igetter) override

void	dqmEndLuminosityBlock (DQMStore::IBooker &ibooker, DQMStore::IGetter &igetter, const edm::LuminosityBlock &lumiSeg, const edm::EventSetup &c) override

double	xySymmetry (const edm::ParameterSet &ps, std::string test_name, std::vector< std::pair< int, double > > &deadChannels, std::vector< std::pair< int, double > > &statDev, bool &enoughStats)

Protected Member Functions inherited from edm::ProducerBase
ProducesCollector	producesCollector ()

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 () noexcept

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

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)

Private Member Functions
int	compareWithStrip (TH2F *histo, std::string test, int binStrip, int nBins, int axis, double avg, const edm::ParameterSet &ps, std::vector< std::pair< int, double > > &deadChannels)

double	getAvrg (TH2F *h2f, std::string test, int axis, int nBins, int binStrip, int avrgMode)

void	getBinCoordinateOnAxisWithValue (TH2F *h2f, double content, int &coord, int axis)

void	printDeadChannels (const std::vector< std::pair< int, double > > &deadChannels, TH2F *h2f, const std::vector< std::pair< int, double > > &statDev, std::string test_name)

Private Attributes
TFile *	file_

L1TOccupancyClientHistogramService *	hservice_

std::map< std::string, MonitorElement * >	meCertification

std::map< std::string, MonitorElement * >	meDifferential

std::map< std::string, MonitorElement * >	meResults

std::vector< edm::ParameterSet * >	mValidTests

edm::ParameterSet	parameters_

std::vector< edm::ParameterSet >	tests_

bool	verbose_

Additional Inherited Members
Public Types inherited from DQMEDHarvester
typedef dqm::harvesting::DQMStore	DQMStore

typedef dqm::harvesting::MonitorElement	MonitorElement

Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Attributes inherited from DQMEDHarvester
DQMStore *	dqmstore_

edm::GetterOfProducts< DQMToken >	jobmegetter_

edm::EDPutTokenT< DQMToken >	jobToken_

edm::GetterOfProducts< DQMToken >	lumimegetter_

edm::EDPutTokenT< DQMToken >	lumiToken_

edm::GetterOfProducts< DQMToken >	runmegetter_

edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

Definition at line 30 of file L1TOccupancyClient.h.

Constructor & Destructor Documentation

L1TOccupancyClient::L1TOccupancyClient ( const edm::ParameterSet & ps )

Constructor.

Definition at line 28 of file L1TOccupancyClient.cc.

References gather_cfg::cout, and edm::ParameterSet::getParameter().

                                                                 {
   // Get parameters
   parameters_ = ps;
   verbose_ = ps.getParameter<bool>("verbose");
   tests_ = ps.getParameter<std::vector<ParameterSet> >("testParams");
 
   if (verbose_) {
     cout << "[L1TOccupancyClient:] Called constructor" << endl;
   }
 }

L1TOccupancyClient::~L1TOccupancyClient ( )

override

Destructor.

Definition at line 43 of file L1TOccupancyClient.cc.

References gather_cfg::cout.

                                         {
   if (verbose_) {
     cout << "[L1TOccupancyClient:] Called destructor" << endl;
   }
 }

Member Function Documentation

void L1TOccupancyClient::book	(	DQMStore::IBooker &	ibooker,
		DQMStore::IGetter &	igetter
	)

protected

Definition at line 56 of file L1TOccupancyClient.cc.

References dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), gather_cfg::cout, edm::ParameterSet::getParameter(), visualization-live-secondInstance_cfg::m, dqm::impl::MonitorElement::Reset(), dqm::implementation::NavigatorBase::setCurrentFolder(), dqm::impl::MonitorElement::setTitle(), and runGCPTkAlMap::title.

                                                                                 {
   hservice_ = new L1TOccupancyClientHistogramService(parameters_, ibooker, verbose_);
 
   if (verbose_) {
     cout << "[L1TOccupancyClient:] Called beginRun" << endl;
 
     // In verbose mode we will produce an extra output file with several tests
     file_ = TFile::Open("DQM_L1TOccupancyClient_Snapshots_LS.root", "RECREATE");
   }
 
   ibooker.setCurrentFolder("L1T/L1TOccupancy/");
   //dbe_->setCurrentFolder("L1T/L1TOccupancy/Results");
   //dbe_->setCurrentFolder("L1T/L1TOccupancy/BadCellValues");
   //dbe_->setCurrentFolder("L1T/L1TOccupancy/Certification");
 
   // Loop over all tests in defined
   for (vector<ParameterSet>::iterator it = tests_.begin(); it != tests_.end(); it++) {
     // If the test algorithm is XYSymmetry we create the necessary histograms
     if ((*it).getUntrackedParameter<string>("algoName", "XYSymmetry") == "XYSymmetry") {
       // Getting Parameters for the test
       string testName = (*it).getParameter<string>("testName");
       ParameterSet algoParameters = (*it).getParameter<ParameterSet>("algoParams");
       string histPath = algoParameters.getParameter<string>("histPath");
 
       if (verbose_) {
         cout << "[L1TOccupancyClient:] Monitored histogram path: " << histPath << endl;
 
         // Creating verbose file directory structure
         // test_name/test_name_Results,
         // test_name/test_name_Histos
         // TDirectory *td  = file_->mkdir(testName.c_str()             ,testName.c_str());
         //FIXME: sub never used gcc361 warning
         //TDirectory *sub = td   ->mkdir((testName+"_Results").c_str(),string("_Results").c_str());
 
         //sub = td->mkdir((testName+"_Histos").c_str()      ,(testName+"_Histos").c_str());
         //sub = td->mkdir((testName+"_Histos_AllLS").c_str(),(testName+"_Histos_AllLS").c_str());
       }
 
       // Load histograms in service instance
       if (hservice_->loadHisto(igetter, testName, histPath)) {
         // Mask channels specified in python file
         hservice_->setMaskedBins(testName, algoParameters.getParameter<vector<ParameterSet> >("maskedAreas"));
 
         // Book MonitorElements
         // * Test results
         ibooker.setCurrentFolder("L1T/L1TOccupancy/Results");
         string title = testName;
         MonitorElement* m = ibooker.book2D(title.c_str(), hservice_->getDifferentialHistogram(testName));
         m->setTitle(title);
         m->Reset();
         meResults[title] = m;
 
         // * Which cells are masked as bad
         ibooker.setCurrentFolder("L1T/L1TOccupancy/HistogramDiff");
         title = testName;
         m = ibooker.book2D(title.c_str(), hservice_->getDifferentialHistogram(testName));
         m->Reset();
         m->setTitle(title);
         meDifferential[title] = m;
 
         // * Fraction of bad cells
         ibooker.setCurrentFolder("L1T/L1TOccupancy/Certification");
         title = testName;
         m = ibooker.book1D(title.c_str(), title.c_str(), 2500, -.5, 2500. - .5);
         m->setTitle(title);
         meCertification[title] = m;
 
         mValidTests.push_back(&(*it));
       }
     }
   }
 }

int L1TOccupancyClient::compareWithStrip	(	TH2F *	histo,
		std::string	test,
		int	binStrip,
		int	nBins,
		int	axis,
		double	avg,
		const edm::ParameterSet &	ps,
		std::vector< std::pair< int, double > > &	deadChannels
	)

private

Definition at line 731 of file L1TOccupancyClient.cc.

References funct::abs(), gather_cfg::cout, alignCSCRings::e, edm::ParameterSet::getUntrackedParameter(), mps_fire::i, infinity, and submitPVValidationJobs::params.

                                                                                 {
   int dead = 0;
 
   //
   if (iAxis == 1) {
     // Get and set parameters for working curves
     TF1* fmuup = new TF1("fmuup", "TMath::Log(TMath::PoissonI(x,[0])/TMath::PoissonI(x,[1]))", -10000., 10000.);
     TF1* fmulow = new TF1("fmulow", "TMath::Log(TMath::PoissonI(x,[0])/TMath::PoissonI(x,[1]))", -10000., 10000.);
     fmuup->SetParameter(0, iAvg * iPS.getUntrackedParameter<double>("factorup", 2.0));
     fmuup->SetParameter(1, iAvg);
     fmulow->SetParameter(0, iAvg * iPS.getUntrackedParameter<double>("factorlow", 0.1));
     fmulow->SetParameter(1, iAvg);
 
     TF1* fchi = new TF1("fchi", "[0]*x**2+[1]*x+[2]", 0., 1500.);
 
     // Evaluate sigma up
     vector<double> defaultChi2up;
     defaultChi2up.push_back(5.45058e-05);
     defaultChi2up.push_back(0.268756);
     defaultChi2up.push_back(-11.7515);
 
     vector<double> params = iPS.getUntrackedParameter<vector<double> >("params_chi2_up", defaultChi2up);
     for (unsigned int i = 0; i < params.size(); i++) {
       fchi->SetParameter(i, params[i]);
     }
     double sigma_up = fchi->Eval(iAvg);
 
     // Evaluate sigma low
     vector<double> defaultChi2low;
     defaultChi2low.push_back(4.11095e-05);
     defaultChi2low.push_back(0.577451);
     defaultChi2low.push_back(-10.378);
 
     params = iPS.getUntrackedParameter<vector<double> >("params_chi2_low", defaultChi2low);
     for (unsigned int i = 0; i < params.size(); i++) {
       fchi->SetParameter(i, params[i]);
     }
     double sigma_low = fchi->Eval(iAvg);
 
     if (verbose_) {
       cout << "binstrip= " << iBinStrip << ", sigmaup= " << sigma_up << ", sigmalow= " << sigma_low << endl;
     }
 
     for (int i = 1; i <= iNBins; i++) {
       if (verbose_) {
         cout << "    " << i << " binContent: up:" << fmuup->Eval(iHist->GetBinContent(iBinStrip, i))
              << " low: " << fmulow->Eval(iHist->GetBinContent(iBinStrip, i)) << endl;
       }
 
       // Evaluate chi2 for cells
       double muup = fmuup->Eval(iHist->GetBinContent(iBinStrip, i));
       double mulow = fmulow->Eval(iHist->GetBinContent(iBinStrip, i));
 
       // If channel is masked -> set it to value -1
       if (hservice_->isMasked(iTestName, iBinStrip, i)) {
         oChannels.push_back(pair<int, double>(iHist->GetBin(iBinStrip, i), -1.0));
       }
       //else perform test
       else if (muup > sigma_up || mulow > sigma_low ||
                ((fabs(muup) == std::numeric_limits<double>::infinity()) &&
                 (fabs(mulow) == std::numeric_limits<double>::infinity()))) {
         dead++;
         oChannels.push_back(
             pair<int, double>(iHist->GetBin(iBinStrip, i), abs(iHist->GetBinContent(iBinStrip, i) - iAvg) / iAvg));
       }
     }
   }
   //
   else if (iAxis == 2) {
     //get and set parameters for working curves
     TF1* fmuup = new TF1("fmuup", "TMath::Log(TMath::PoissonI(x,[0])/TMath::PoissonI(x,[1]))", -10000., 10000.);
     TF1* fmulow = new TF1("fmulow", "TMath::Log(TMath::PoissonI(x,[0])/TMath::PoissonI(x,[1]))", -10000., 10000.);
     fmuup->SetParameter(0, iAvg * iPS.getUntrackedParameter<double>("factorup", 2.0));
     fmuup->SetParameter(1, iAvg);
     fmulow->SetParameter(0, iAvg * iPS.getUntrackedParameter<double>("factorlow", 0.1));
     fmulow->SetParameter(1, iAvg);
 
     TF1* fchi = new TF1("fchi", "[0]*x**2+[1]*x+[2]", 0., 1500.);
 
     // Evaluate sigma up
     vector<double> defaultChi2up;
     defaultChi2up.push_back(5.45058e-05);
     defaultChi2up.push_back(0.268756);
     defaultChi2up.push_back(-11.7515);
 
     vector<double> params = iPS.getUntrackedParameter<vector<double> >("params_chi2_up", defaultChi2up);
     for (unsigned int i = 0; i < params.size(); i++) {
       fchi->SetParameter(i, params[i]);
     }
     double sigma_up = fchi->Eval(iAvg);
 
     // Evaluate sigma low
     vector<double> defaultChi2low;
     defaultChi2low.push_back(4.11095e-05);
     defaultChi2low.push_back(0.577451);
     defaultChi2low.push_back(-10.378);
 
     params = iPS.getUntrackedParameter<vector<double> >("params_chi2_low", defaultChi2low);
     for (unsigned int i = 0; i < params.size(); i++) {
       fchi->SetParameter(i, params[i]);
     }
     double sigma_low = fchi->Eval(iAvg);
 
     if (verbose_) {
       cout << "binstrip= " << iBinStrip << ", sigmaup= " << sigma_up << ", sigmalow= " << sigma_low << endl;
     }
 
     for (int i = 1; i <= iNBins; i++) {
       if (verbose_) {
         cout << "    " << i << " binContent: up:" << fmuup->Eval(iHist->GetBinContent(i, iBinStrip))
              << " low: " << fmulow->Eval(iHist->GetBinContent(i, iBinStrip)) << endl;
       }
 
       //evaluate chi2 for cells
       double muup = fmuup->Eval(iHist->GetBinContent(i, iBinStrip));
       double mulow = fmulow->Eval(iHist->GetBinContent(i, iBinStrip));
 
       //if channel is masked -> set it to value -1
       if (hservice_->isMasked(iTestName, i, iBinStrip)) {
         oChannels.push_back(pair<int, double>(iHist->GetBin(iBinStrip, i), -1.0));
       }
       //else perform test
       else if (muup > sigma_up || mulow > sigma_low ||
                ((fabs(muup) == std::numeric_limits<double>::infinity()) &&
                 (fabs(mulow) == std::numeric_limits<double>::infinity()))) {
         dead++;
         oChannels.push_back(
             pair<int, double>(iHist->GetBin(i, iBinStrip), abs(iHist->GetBinContent(i, iBinStrip) - iAvg) / iAvg));
       }
     }
   } else {
     if (verbose_) {
       cout << "invalid axis" << endl;
     }
   }
 
   return dead;
 }

void L1TOccupancyClient::dqmEndJob	(	DQMStore::IBooker &	ibooker,
		DQMStore::IGetter &	igetter
	)

overrideprotectedvirtual

Implements DQMEDHarvester.

Definition at line 136 of file L1TOccupancyClient.cc.

References newFWLiteAna::bin, gather_cfg::cout, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), mps_fire::i, str, and ctpps_dqm_sourceclient-live_cfg::test.

                                                                                      {
   book(ibooker, igetter);
 
   if (verbose_) {
     cout << "[L1TOccupancyClient:] Called endRun()" << endl;
   }
 
   // Loop over every test in python
   for (std::vector<ParameterSet*>::iterator it = mValidTests.begin(); it != mValidTests.end(); it++) {
     ParameterSet& test = (**it);
     string algo_name = test.getUntrackedParameter<string>("algoName", "XYSymmetry");
     string test_name = test.getParameter<string>("testName");
 
     if (verbose_) {
       cout << "[L1TOccupancyClient:] Starting calculations for: " << algo_name << " on: " << test_name << endl;
     }
 
     if (algo_name == "XYSymmetry") {
       ParameterSet ps = (**it).getParameter<ParameterSet>("algoParams");
       string histPath = ps.getParameter<string>("histPath");
 
       vector<pair<int, double> > deadChannels;
       vector<pair<int, double> > statDev;
       bool enoughStats = false;
 
       // Make final block
       hservice_->updateHistogramEndRun(test_name);
 
       // Perform the test
       double dead = xySymmetry(ps, test_name, deadChannels, statDev, enoughStats);
       stringstream str;
       str << test_name << "_cumu_LS_EndRun";
 
       if (verbose_) {
         TH2F* cumulative_save = (TH2F*)hservice_->getDifferentialHistogram(test_name)->Clone(str.str().c_str());
 
         cumulative_save->SetTitle(str.str().c_str());
 
         TDirectory* td = file_->GetDirectory(test_name.c_str());
 
         td->cd(string(test_name + "_Histos_AllLS").c_str());
 
         cumulative_save->Write();
       }
       // If we have enough statistics, we can write test result
       if (enoughStats) {
         // Make the result histogram
         printDeadChannels(deadChannels, meResults[test_name]->getTH2F(), statDev, test_name);
 
         if (verbose_) {
           TH2F* cumulative_save = (TH2F*)hservice_->getDifferentialHistogram(test_name)->Clone(str.str().c_str());
           cumulative_save->SetTitle(str.str().c_str());
           TDirectory* td = file_->GetDirectory(("DQM_L1TOccupancyClient_Snapshots_LS.root:/" + test_name).c_str());
           td->cd(string(test_name + "_Histos").c_str());
           cumulative_save->Write();
 
           // save the result histo
           TH2F* h2f = meResults[test_name]->getTH2F();
           stringstream str2;
           str2 << test_name << "_result_LS_EndRun";
           TH2F* dead_save = (TH2F*)h2f->Clone(str2.str().c_str());
 
           td->cd(string(test_name + "_Results").c_str());
           dead_save->SetTitle(str2.str().c_str());
           dead_save->Write();
         }
 
         // Updating test results
         meDifferential[test_name]->Reset();
         meDifferential[test_name]->getTH2F()->Add(hservice_->getDifferentialHistogram(test_name));
 
         vector<int> lsCertification = hservice_->getLSCertification(test_name);
 
         // Fill fraction of dead channels
         for (unsigned int i = 0; i < lsCertification.size(); i++) {
           int bin = meCertification[test_name]->getTH1()->FindBin(lsCertification[i]);
           meCertification[test_name]->getTH1()->SetBinContent(bin, 1 - dead);
         }
 
         // Reset differential histo
         hservice_->resetHisto(test_name);
 
         if (verbose_) {
           cout << "Now we have enough statstics for " << test_name << endl;
         }
 
       } else {
         if (verbose_) {
           cout << "we don't have enough statstics for " << test_name << endl;
         }
 
         // Getting LS which this test monitored
         vector<int> lsCertification = hservice_->getLSCertification(test_name);
 
         // Fill fraction of dead channels
         for (unsigned int i = 0; i < lsCertification.size(); i++) {
           int bin = meCertification[test_name]->getTH1()->FindBin(lsCertification[i]);
           meCertification[test_name]->getTH1()->SetBinContent(bin, -1);
         }
       }
     } else {
       if (verbose_) {
         cout << "No valid algorithm" << std::endl;
       }
     }
   }
 
   if (verbose_) {
     file_->Close();
   }
 
   delete hservice_;
 }

void L1TOccupancyClient::dqmEndLuminosityBlock	(	DQMStore::IBooker &	ibooker,
		DQMStore::IGetter &	igetter,
		const edm::LuminosityBlock &	lumiSeg,
		const edm::EventSetup &	c
	)

overrideprotectedvirtual

Reimplemented from DQMEDHarvester.

Definition at line 265 of file L1TOccupancyClient.cc.

References newFWLiteAna::bin, gather_cfg::cout, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), mps_fire::i, edm::LuminosityBlockBase::id(), edm::LuminosityBlockID::luminosityBlock(), str, and ctpps_dqm_sourceclient-live_cfg::test.

                                                                        {
   book(ibooker, igetter);
 
   int eventLS = lumiSeg.id().luminosityBlock();
 
   if (verbose_) {
     cout << "[L1TOccupancyClient:] Called endLuminosityBlock()" << endl;
     cout << "[L1TOccupancyClient:] Lumisection: " << eventLS << endl;
   }
 
   // Loop over every test in python
   for (std::vector<ParameterSet*>::const_iterator it = mValidTests.begin(); it != mValidTests.end(); it++) {
     ParameterSet& test = (**it);
     string algo_name = test.getUntrackedParameter<string>("algoName", "XYSymmetry");
     string test_name = test.getParameter<string>("testName");
 
     if (verbose_) {
       cout << "[L1TOccupancyClient:] Starting calculations for " << algo_name << " on:" << test_name << endl;
     }
 
     if (algo_name == "XYSymmetry") {
       ParameterSet ps = (**it).getParameter<ParameterSet>("algoParams");
       string histPath = ps.getParameter<string>("histPath");
 
       vector<pair<int, double> > deadChannels;
       vector<pair<int, double> > statDev;
       bool enoughStats = false;
 
       // Update histo's data with data of this LS
       hservice_->updateHistogramEndLS(igetter, test_name, histPath, eventLS);
 
       // Perform the test
       double dead = xySymmetry(ps, test_name, deadChannels, statDev, enoughStats);
       stringstream str;
       str << test_name << "_cumu_LS_" << eventLS;
 
       if (verbose_) {
         TH2F* cumulative_save = (TH2F*)hservice_->getDifferentialHistogram(test_name)->Clone(str.str().c_str());
         cumulative_save->SetTitle(str.str().c_str());
         TDirectory* td = file_->GetDirectory(test_name.c_str());
         td->cd(string(test_name + "_Histos_AllLS").c_str());
         cumulative_save->Write();
       }
 
       // If we have enough statistics, we can write test result
       if (enoughStats) {
         // Make the result histogram
         printDeadChannels(deadChannels, meResults[test_name]->getTH2F(), statDev, test_name);
 
         if (verbose_) {
           TH2F* cumulative_save = (TH2F*)hservice_->getDifferentialHistogram(test_name)->Clone(str.str().c_str());
           cumulative_save->SetTitle(str.str().c_str());
           TDirectory* td = file_->GetDirectory(("DQM_L1TOccupancyClient_Snapshots_LS.root:/" + test_name).c_str());
           td->cd(string(test_name + "_Histos").c_str());
           cumulative_save->Write();
 
           // save the result histo
           TH2F* h2f = meResults[test_name]->getTH2F();
           stringstream str2;
           str2 << test_name << "_result_LS_" << eventLS;
           TH2F* dead_save = (TH2F*)h2f->Clone(str2.str().c_str());
 
           td->cd(string(test_name + "_Results").c_str());
           dead_save->SetTitle(str2.str().c_str());
           dead_save->Write();
         }
 
         // Updating test results
         meDifferential[test_name]->Reset();
         meDifferential[test_name]->getTH2F()->Add(hservice_->getDifferentialHistogram(test_name));
 
         vector<int> lsCertification = hservice_->getLSCertification(test_name);
 
         // Fill fraction of dead channels
         for (unsigned int i = 0; i < lsCertification.size(); i++) {
           int bin = meCertification[test_name]->getTH1()->FindBin(lsCertification[i]);
           meCertification[test_name]->getTH1()->SetBinContent(bin, 1 - dead);
         }
 
         // Reset differential histo
         hservice_->resetHisto(test_name);
 
         if (verbose_) {
           cout << "Now we have enough statstics for " << test_name << endl;
         }
 
       } else {
         if (verbose_) {
           cout << "we don't have enough statstics for " << test_name << endl;
         }
       }
     } else {
       if (verbose_) {
         cout << "No valid algorithm" << std::endl;
       }
     }
   }
 }

double L1TOccupancyClient::getAvrg	(	TH2F *	h2f,
		std::string	test,
		int	axis,
		int	nBins,
		int	binStrip,
		int	avrgMode
	)

private

Definition at line 597 of file L1TOccupancyClient.cc.

References gather_cfg::cout, timingPdfMaker::histo, mps_fire::i, and makeHLTPrescaleTable::values.

                                                                                                                     {
   double avg = 0.0;
   TH1D* proj = nullptr;
   TH2F* histo = new TH2F(*iHist);
 
   std::vector<double> values;
   int marked;
 
   if (iAxis == 1) {
     switch (iAvgMode) {
       // arithmetic average
       case 1:
         marked = hservice_->maskBins(iTestName, histo, iBinStrip, iAxis);
         proj = histo->ProjectionX();
         avg = proj->GetBinContent(iBinStrip) / (iNBins - marked);
         break;
 
       // median
       case 2:
         marked = hservice_->maskBins(iTestName, histo, iBinStrip, iAxis);
         proj = histo->ProjectionY("_py", iBinStrip, iBinStrip);
         for (int i = 0; i < iNBins; i++) {
           values.push_back(proj->GetBinContent(i + 1));
         }
         avg = TMath::Median(iNBins, &values[0]);
         break;
       default:
         if (verbose_) {
           cout << "Invalid averaging mode!" << endl;
         }
         break;
     }
   } else if (iAxis == 2) {
     switch (iAvgMode) {
       // arithmetic average
       case 1:
         marked = hservice_->maskBins(iTestName, histo, iBinStrip, iAxis);
         proj = histo->ProjectionY();
         avg = proj->GetBinContent(iBinStrip) / (iNBins - marked);
         break;
       // median
       case 2:
         marked = hservice_->maskBins(iTestName, histo, iBinStrip, iAxis);
         proj = histo->ProjectionX("_px", iBinStrip, iBinStrip);
         for (int i = 0; i < iNBins; i++) {
           values.push_back(proj->GetBinContent(i + 1));
         }
 
         avg = TMath::Median(iNBins, &values[0]);
         break;
       default:
         if (verbose_) {
           cout << "invalid averaging mode!" << endl;
         }
         break;
     }
   } else {
     if (verbose_) {
       cout << "invalid axis" << endl;
     }
   }
   delete proj;
   delete histo;
   return avg;
 }

void L1TOccupancyClient::getBinCoordinateOnAxisWithValue	(	TH2F *	h2f,
		double	content,
		int &	coord,
		int	axis
	)

private

Definition at line 886 of file L1TOccupancyClient.cc.

                                                                                                                    {
   int nBinsX = iHist->GetNbinsX();  //actual number of bins x
   int nBinsY = iHist->GetNbinsY();  //actual number of bins y
 
   if (iAxis == 1) {
     int global = iHist->GetXaxis()->FindFixBin(iValue);
 
     // If parameter exceeds axis' value: set to maximum number of bins in x-axis
     if (global > nBinsX * nBinsY) {
       global = iHist->GetXaxis()->GetLast();
     }
 
     // Get coordinates of bin
     int y, z;
     iHist->GetBinXYZ(global, oBinCoordinate, y, z);
   } else if (iAxis == 2) {
     int global = iHist->GetYaxis()->FindFixBin(iValue);
 
     // If parameter exceeds axis' value: set to maximum number of bins in x-axis
     if (global > nBinsX * nBinsY) {
       global = iHist->GetYaxis()->GetLast();
     }
 
     // Get coordinates of bin
     int x, z;
     iHist->GetBinXYZ(global, x, oBinCoordinate, z);
   }
 }

void L1TOccupancyClient::printDeadChannels	(	const std::vector< std::pair< int, double > > &	deadChannels,
		TH2F *	h2f,
		const std::vector< std::pair< int, double > > &	statDev,
		std::string	test_name
	)

private

Definition at line 672 of file L1TOccupancyClient.cc.

References newFWLiteAna::bin, HLT_FULL_cff::chi2, gather_cfg::cout, gpuVertexFinder::printf(), and trackerHitRTTI::vector.

                                                              {
   // Reset the dead channels histogram
   oHistDeadChannels->Reset();
   if (verbose_) {
     cout << "suspect or masked channels of " << iTestName << ": ";
   }
 
   int x, y, z;
   float chi2 = 0.0;
 
   // put all bad (value=1) and masked (value=-1) cells in histo
   for (std::vector<pair<int, double> >::const_iterator it = iDeadChannels.begin(); it != iDeadChannels.end(); it++) {
     int bin = (*it).first;
     oHistDeadChannels->GetBinXYZ(bin, x, y, z);
 
     if (hservice_->isMasked(iTestName, x, y)) {
       oHistDeadChannels->SetBinContent(bin, -1);
       if (verbose_) {
         printf("(%4i,%4i) Masked\n", x, y);
       }
     } else {
       oHistDeadChannels->SetBinContent(bin, 1);
       if (verbose_) {
         printf("(%4i,%4i) Failed test\n", x, y);
       }
     }
   }
 
   // FIXME: Is this needed?
   for (std::vector<pair<int, double> >::const_iterator it = statDev.begin(); it != statDev.end(); it++) {
     double dev = (*it).second;
     chi2 += dev;
   }
   //put total chi2 in float
 
   if (verbose_) {
     cout << "total number of suspect channels: " << (iDeadChannels.size() - (hservice_->getNBinsMasked(iTestName)))
          << endl;
   }
 }

double L1TOccupancyClient::xySymmetry	(	const edm::ParameterSet &	ps,
		std::string	test_name,
		std::vector< std::pair< int, double > > &	deadChannels,
		std::vector< std::pair< int, double > > &	statDev,
		bool &	enoughStats
	)

protected

Definition at line 388 of file L1TOccupancyClient.cc.

References gather_cfg::cout, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), mps_fire::i, dqmiolumiharvest::j, isotrackApplyRegressor::k, Max(), and submitPVValidationJobs::params.

                                                          {
   // Getting differential histogram for this this thes
   TH2F* diffHist = hservice_->getDifferentialHistogram(iTestName);
 
   int pAxis = ps.getUntrackedParameter<int>("axis", 1);
   int pAverageMode = ps.getUntrackedParameter<int>("averageMode", 2);  // 1=arith. mean, 2=median
   int nBinsX = diffHist->GetNbinsX();                                  // actual number of bins x
   int nBinsY = diffHist->GetNbinsY();                                  // actual number of bins y
 
   // Axis==1 : Means symmetry axis is vertical
   if (pAxis == 1) {
     int maxBinStrip, centralBinStrip;  // x-coordinate of strips
 
     maxBinStrip = nBinsX;
 
     // If takeCenter=true  determine central bin of the pAxis
     // If takeCenter=false determine the bin to use based user input
     if (ps.getUntrackedParameter<bool>("takeCenter", true)) {
       centralBinStrip = nBinsX / 2 + 1;
     } else {
       double pAxisSymmetryValue = ps.getParameter<double>("axisSymmetryValue");
       getBinCoordinateOnAxisWithValue(diffHist, pAxisSymmetryValue, centralBinStrip, 1);
     }
 
     // Assuming odd number of strips --> first comparison is middle strip to itself
     int upBinStrip = centralBinStrip;
     int lowBinStrip = centralBinStrip;
 
     // If even number decrease lowBinstrip by one
     if (nBinsX % 2 == 0) {
       lowBinStrip--;
     }
 
     // Do we have enough statistics? Min(Max(strip_i,strip_j))>threshold
     std::unique_ptr<double[]> maxAvgs(new double[maxBinStrip - upBinStrip + 1]);
 
     int nActualStrips = 0;  //number of strips that are not fully masked
     for (int i = 0, j = upBinStrip, k = lowBinStrip; j <= maxBinStrip; i++, j++, k--) {
       double avg1 = getAvrg(diffHist, iTestName, pAxis, nBinsY, j, pAverageMode);
       double avg2 = getAvrg(diffHist, iTestName, pAxis, nBinsY, k, pAverageMode);
 
       // Protection for when both strips are masked
       if (!hservice_->isStripMasked(iTestName, j, pAxis) && !hservice_->isStripMasked(iTestName, k, pAxis)) {
         maxAvgs[i] = TMath::Max(avg1, avg2);
         nActualStrips++;
       }
     }
 
     vector<double> defaultMu0up;
     defaultMu0up.push_back(13.7655);
     defaultMu0up.push_back(184.742);
     defaultMu0up.push_back(50735.3);
     defaultMu0up.push_back(-97.6793);
 
     TF1 tf("myFunc", "[0]*(TMath::Log(x*[1]+[2]))+[3]", 10., 11000.);
     vector<double> params = ps.getUntrackedParameter<vector<double> >("params_mu0_up", defaultMu0up);
     for (unsigned int i = 0; i < params.size(); i++) {
       tf.SetParameter(i, params[i]);
     }
     int statsup = (int)tf.Eval(hservice_->getNBinsHistogram(iTestName));
 
     vector<double> defaultMu0low;
     defaultMu0low.push_back(2.19664);
     defaultMu0low.push_back(1.94546);
     defaultMu0low.push_back(-99.3263);
     defaultMu0low.push_back(19.388);
 
     params = ps.getUntrackedParameter<vector<double> >("params_mu0_low", defaultMu0low);
     for (unsigned int i = 0; i < params.size(); i++) {
       tf.SetParameter(i, params[i]);
     }
     int statslow = (int)tf.Eval(hservice_->getNBinsHistogram(iTestName));
 
     if (verbose_) {
       cout << "nbins: " << hservice_->getNBinsHistogram(iTestName) << endl;
       cout << "statsup= " << statsup << ", statslow= " << statslow << endl;
     }
 
     enoughStats = TMath::MinElement(nActualStrips, maxAvgs.get()) > TMath::Max(statsup, statslow);
     if (verbose_) {
       cout << "stats: " << TMath::MinElement(nActualStrips, maxAvgs.get())
            << ", statsAvg: " << diffHist->GetEntries() / hservice_->getNBinsHistogram(iTestName)
            << ", threshold: " << TMath::Max(statsup, statslow) << endl;
     }
 
     //if enough statistics
     //make the test
     if (enoughStats) {
       for (; upBinStrip <= maxBinStrip; upBinStrip++, lowBinStrip--) {
         double avg = getAvrg(diffHist, iTestName, pAxis, nBinsY, upBinStrip, pAverageMode);
         compareWithStrip(
             diffHist, iTestName, lowBinStrip, nBinsY, pAxis, avg, ps, deadChannels);  //compare with lower side
 
         avg = getAvrg(diffHist, iTestName, pAxis, nBinsY, lowBinStrip, pAverageMode);
         compareWithStrip(
             diffHist, iTestName, upBinStrip, nBinsY, pAxis, avg, ps, deadChannels);  //compare with upper side
       }
     }
   }
 
   // pAxis==2 : Means symetry pAxis is horizontal
   else if (pAxis == 2) {
     int maxBinStrip, centralBinStrip;  //x-coordinate of strips
 
     maxBinStrip = nBinsY;
 
     // Determine center of diagram: either with set pAxis or middle of diagram
     if (ps.getUntrackedParameter<bool>("takeCenter", true)) {
       centralBinStrip = nBinsY / 2 + 1;
     } else {
       double pAxisSymmetryValue = ps.getParameter<double>("axisSymmetryValue");
       getBinCoordinateOnAxisWithValue(diffHist, pAxisSymmetryValue, centralBinStrip, 2);
     }
 
     //assuming odd number of strips --> first comparison is middle strip to itself
     int lowBinStrip = centralBinStrip, upBinStrip = centralBinStrip;
 
     //if even number
     if (nBinsX % 2 == 0) {
       //decrease lowBinstrip by one
       lowBinStrip--;
     }
 
     //do we have enough statistics? Min(Max(strip_i,strip_j))>threshold
     std::unique_ptr<double[]> maxAvgs(new double[maxBinStrip - upBinStrip + 1]);
     int nActualStrips = 0;
     for (int i = 0, j = upBinStrip, k = lowBinStrip; j <= maxBinStrip; i++, j++, k--) {
       double avg1 = getAvrg(diffHist, iTestName, pAxis, nBinsX, j, pAverageMode);
       double avg2 = getAvrg(diffHist, iTestName, pAxis, nBinsX, k, pAverageMode);
       if (!hservice_->isStripMasked(iTestName, j, pAxis) && !hservice_->isStripMasked(iTestName, k, pAxis)) {
         maxAvgs[i] = TMath::Max(avg1, avg2);
         nActualStrips++;
       }
     }
 
     vector<double> defaultMu0up;
     defaultMu0up.push_back(13.7655);
     defaultMu0up.push_back(184.742);
     defaultMu0up.push_back(50735.3);
     defaultMu0up.push_back(-97.6793);
 
     vector<double> params = ps.getUntrackedParameter<std::vector<double> >("params_mu0_up", defaultMu0up);
     TF1 tf("myFunc", "[0]*(TMath::Log(x*[1]+[2]))+[3]", 10., 11000.);
     for (unsigned int i = 0; i < params.size(); i++) {
       tf.SetParameter(i, params[i]);
     }
     int statsup = (int)tf.Eval(hservice_->getNBinsHistogram(iTestName));
 
     vector<double> defaultMu0low;
     defaultMu0low.push_back(2.19664);
     defaultMu0low.push_back(1.94546);
     defaultMu0low.push_back(-99.3263);
     defaultMu0low.push_back(19.388);
 
     params = ps.getUntrackedParameter<std::vector<double> >("params_mu0_low", defaultMu0low);
     for (unsigned int i = 0; i < params.size(); i++) {
       tf.SetParameter(i, params[i]);
     }
     int statslow = (int)tf.Eval(hservice_->getNBinsHistogram(iTestName));
     if (verbose_) {
       cout << "statsup= " << statsup << ", statslow= " << statslow << endl;
     }
     enoughStats = TMath::MinElement(nActualStrips, maxAvgs.get()) > TMath::Max(statsup, statslow);
     if (verbose_) {
       cout << "stats: " << TMath::MinElement(nActualStrips, maxAvgs.get())
            << ", statsAvg: " << diffHist->GetEntries() / hservice_->getNBinsHistogram(iTestName)
            << ", threshold: " << TMath::Max(statsup, statslow) << endl;
     }
 
     //if we have enough statistics
     //make the test
     if (enoughStats) {
       for (; upBinStrip <= maxBinStrip; upBinStrip++, lowBinStrip--) {
         double avg = getAvrg(diffHist, iTestName, pAxis, nBinsX, upBinStrip, pAverageMode);
         compareWithStrip(
             diffHist, iTestName, lowBinStrip, nBinsX, pAxis, avg, ps, deadChannels);  //compare with lower side
 
         avg = getAvrg(diffHist, iTestName, pAxis, nBinsX, lowBinStrip, pAverageMode);
         compareWithStrip(
             diffHist, iTestName, upBinStrip, nBinsX, pAxis, avg, ps, deadChannels);  //compare with upper side
       }
     }
   } else {
     if (verbose_) {
       cout << "Invalid axis" << endl;
     }
   }
 
   return (deadChannels.size() - hservice_->getNBinsMasked(iTestName)) * 1.0 / hservice_->getNBinsHistogram(iTestName);
 }

Member Data Documentation

TFile* L1TOccupancyClient::file_

private

Definition at line 56 of file L1TOccupancyClient.h.

L1TOccupancyClientHistogramService* L1TOccupancyClient::hservice_

private

Definition at line 55 of file L1TOccupancyClient.h.

std::map<std::string, MonitorElement*> L1TOccupancyClient::meCertification

private

Definition at line 67 of file L1TOccupancyClient.h.

std::map<std::string, MonitorElement*> L1TOccupancyClient::meDifferential

private

Definition at line 66 of file L1TOccupancyClient.h.

std::map<std::string, MonitorElement*> L1TOccupancyClient::meResults

private

Definition at line 65 of file L1TOccupancyClient.h.

std::vector<edm::ParameterSet*> L1TOccupancyClient::mValidTests

private

Definition at line 63 of file L1TOccupancyClient.h.

edm::ParameterSet L1TOccupancyClient::parameters_

private

Definition at line 54 of file L1TOccupancyClient.h.

Referenced by Mixins._TypedParameterizable::clone(), Types.PSet::clone(), Types.EDAlias::clone(), and Mixins._TypedParameterizable::copy().

std::vector<edm::ParameterSet> L1TOccupancyClient::tests_

private

Definition at line 62 of file L1TOccupancyClient.h.

bool L1TOccupancyClient::verbose_

private

Definition at line 59 of file L1TOccupancyClient.h.

Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation