#include <MuonME0DigisHarvestor.h>

Inheritance diagram for MuonME0DigisHarvestor:

Public Member Functions
TH1F *	ComputeBKG (TH1F hist1, TH1F hist2, std::string nameHist)

TProfile *	ComputeEff (TH1F num, TH1F denum, std::string nameHist)

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

	MuonME0DigisHarvestor (const edm::ParameterSet &)
	constructor More...

void	ProcessBooking (DQMStore::IBooker &, DQMStore::IGetter &, std::string nameHist, TH1F num, TH1F den)

void	ProcessBookingBKG (DQMStore::IBooker &ibooker, DQMStore::IGetter &ig, std::string nameHist, TH1F hist, TH1F hist2)

	~MuonME0DigisHarvestor () override
	destructor More...

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

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

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

	DQMEDHarvester ()

virtual void	dqmEndLuminosityBlock (DQMStore::IBooker &, DQMStore::IGetter &, edm::LuminosityBlock const &, edm::EventSetup const &)

void	endJob () final

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

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

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

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

	~DQMEDHarvester () override=default

Public Member Functions inherited from edm::one::EDProducer< edm::Accumulator, edm::EndLuminosityBlockProducer, edm::EndRunProducer, edm::one::WatchLuminosityBlocks, edm::one::WatchRuns, edm::one::SharedResources >
	EDProducer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () 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

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::vector< ModuleDescription const * > &modules, 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

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Attributes
std::string	dbe_path_

Additional Inherited Members
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 Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

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

Protected Attributes inherited from DQMEDHarvester
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

Definition at line 18 of file MuonME0DigisHarvestor.h.

Constructor & Destructor Documentation

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

explicit

constructor

Geometry.

Log messages

Definition at line 32 of file MuonME0DigisHarvestor.cc.

References dbe_path_, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                       {
   dbe_path_ = std::string("MuonME0DigisV/ME0DigisTask/");
 }

MuonME0DigisHarvestor::~MuonME0DigisHarvestor ( )

override

destructor

Definition at line 36 of file MuonME0DigisHarvestor.cc.

36 {}

Member Function Documentation

TH1F * MuonME0DigisHarvestor::ComputeBKG	(	TH1F *	hist1,
		TH1F *	hist2,
		std::string	nameHist
	)

Definition at line 95 of file MuonME0DigisHarvestor.cc.

References stringResolutionProvider_cfi::bin, MillePedeFileConverter_cfg::e, dataset::name, nEvts, Pi, Scenarios_cff::scale, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by ProcessBookingBKG().

                                                                                     {
   std::string name = "rate_" + nameHist;
   hist1->SetName(name.c_str());
   for (int bin = 1; bin <= hist1->GetNbinsX(); ++bin) {
     double R_min = hist1->GetBinCenter(bin) - 0.5 * hist1->GetBinWidth(bin);
     double R_max = hist1->GetBinCenter(bin) + 0.5 * hist1->GetBinWidth(bin);
 
     double Area = TMath::Pi() * (R_max * R_max - R_min * R_min);
     hist1->SetBinContent(bin, (hist1->GetBinContent(bin)) / Area);
     hist1->SetBinError(bin, (hist1->GetBinError(bin)) / Area);
   }
 
   int nEvts = hist2->GetEntries();
   float scale = 6 * 2 * nEvts * 3 * 25e-9;  // New redigitizer saves hits only in the BX range: [-1,+1], so the
                                             // number of background hits has to be divided by 3
   hist1->Scale(1.0 / scale);
   return hist1;
 }

TProfile * MuonME0DigisHarvestor::ComputeEff	(	TH1F *	num,
		TH1F *	denum,
		std::string	nameHist
	)

Definition at line 38 of file MuonME0DigisHarvestor.cc.

References mps_fire::i, dataset::name, mathSSE::sqrt(), AlCaHLTBitMon_QueryRunRegistry::string, groupFilesInBlocks::temp, and fftjetcommon_cfi::title.

Referenced by ProcessBooking().

                                                                                       {
   std::string name = "eff_" + nameHist;
   std::string title = "Digi Efficiency" + std::string(num->GetTitle());
   TProfile *efficHist = new TProfile(name.c_str(),
                                      title.c_str(),
                                      denum->GetXaxis()->GetNbins(),
                                      denum->GetXaxis()->GetXmin(),
                                      denum->GetXaxis()->GetXmax());
 
   for (int i = 1; i <= denum->GetNbinsX(); i++) {
     double nNum = num->GetBinContent(i);
     double nDenum = denum->GetBinContent(i);
     if (nDenum == 0 || nNum == 0) {
       continue;
     }
     if (nNum > nDenum) {
       double temp = nDenum;
       nDenum = nNum;
       nNum = temp;
       edm::LogWarning("MuonME0DigisHarvestor")
           << "Alert! specific bin's num is bigger than denum " << i << " " << nNum << " " << nDenum;
     }
     const double effVal = nNum / nDenum;
     efficHist->SetBinContent(i, effVal);
     efficHist->SetBinEntries(i, 1);
     efficHist->SetBinError(i, 0);
     const double errLo = TEfficiency::ClopperPearson((int)nDenum, (int)nNum, 0.683, false);
     const double errUp = TEfficiency::ClopperPearson((int)nDenum, (int)nNum, 0.683, true);
     const double errVal = (effVal - errLo > errUp - effVal) ? effVal - errLo : errUp - effVal;
     efficHist->SetBinError(i, sqrt(effVal * effVal + errVal * errVal));
   }
   return efficHist;
 }

void MuonME0DigisHarvestor::dqmEndJob	(	DQMStore::IBooker &	ibooker,
		DQMStore::IGetter &	ig
	)

overridevirtual

Implements DQMEDHarvester.

Definition at line 136 of file MuonME0DigisHarvestor.cc.

References dbe_path_, DEFINE_FWK_MODULE, DQMStore::IGetter::get(), mps_fire::i, GEMDetLabel::l_suffix, dataset::name, ProcessBooking(), ProcessBookingBKG(), DQMStore::IGetter::setCurrentFolder(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                    {
   ig.setCurrentFolder(dbe_path_);
 
   const char *l_suffix[6] = {"_l1", "_l2", "_l3", "_l4", "_l5", "_l6"};
   const char *r_suffix[2] = {"-1", "1"};
 
   TString eta_label_den_tot = TString(dbe_path_) + "me0_strip_dg_den_eta_tot";
   TString eta_label_num_tot = TString(dbe_path_) + "me0_strip_dg_num_eta_tot";
   if (ig.get(eta_label_num_tot.Data()) != nullptr && ig.get(eta_label_den_tot.Data()) != nullptr) {
     TH1F *num_vs_eta_tot = (TH1F *)ig.get(eta_label_num_tot.Data())->getTH1F()->Clone();
     num_vs_eta_tot->Sumw2();
     TH1F *den_vs_eta_tot = (TH1F *)ig.get(eta_label_den_tot.Data())->getTH1F()->Clone();
     den_vs_eta_tot->Sumw2();
 
     ProcessBooking(ibooker, ig, "me0_strip_dg_eta_tot", num_vs_eta_tot, den_vs_eta_tot);
 
     delete num_vs_eta_tot;
     delete den_vs_eta_tot;
 
   } else
     edm::LogWarning("MuonME0DigisHarvestor")
         << "Can not find histograms: " << eta_label_num_tot << " or " << eta_label_den_tot;
 
   for (int i = 0; i < 2; i++) {
     for (int j = 0; j < 6; j++) {
       TString eta_label_den = TString(dbe_path_) + "me0_strip_dg_den_eta" + r_suffix[i] + l_suffix[j];
       TString eta_label_num = TString(dbe_path_) + "me0_strip_dg_num_eta" + r_suffix[i] + l_suffix[j];
 
       if (ig.get(eta_label_num.Data()) != nullptr && ig.get(eta_label_den.Data()) != nullptr) {
         TH1F *num_vs_eta = (TH1F *)ig.get(eta_label_num.Data())->getTH1F()->Clone();
         num_vs_eta->Sumw2();
         TH1F *den_vs_eta = (TH1F *)ig.get(eta_label_den.Data())->getTH1F()->Clone();
         den_vs_eta->Sumw2();
 
         std::string r_s = r_suffix[i];
         std::string l_s = l_suffix[j];
         std::string name = "me0_strip_dg_eta" + r_s + l_s;
         ProcessBooking(ibooker, ig, name, num_vs_eta, den_vs_eta);
 
         delete num_vs_eta;
         delete den_vs_eta;
 
       } else
         edm::LogWarning("MuonME0DigisHarvestor")
             << "Can not find histograms: " << eta_label_num << " " << eta_label_den;
     }
   }
 
   TString label_eleBkg = TString(dbe_path_) + "me0_strip_dg_bkgElePos_radius";
   TString label_neuBkg = TString(dbe_path_) + "me0_strip_dg_bkgNeutral_radius";
   TString label_totBkg = TString(dbe_path_) + "me0_strip_dg_bkg_radius_tot";
   TString label_evts = TString(dbe_path_) + "num_evts";
 
   if (ig.get(label_evts.Data()) != nullptr) {
     TH1F *numEvts = (TH1F *)ig.get(label_evts.Data())->getTH1F()->Clone();
 
     if (ig.get(label_eleBkg.Data()) != nullptr) {
       TH1F *eleBkg = (TH1F *)ig.get(label_eleBkg.Data())->getTH1F()->Clone();
       eleBkg->Sumw2();
       ProcessBookingBKG(ibooker, ig, "me0_strip_dg_elePosBkg_rad", eleBkg, numEvts);
 
       delete eleBkg;
     }
     if (ig.get(label_neuBkg.Data()) != nullptr) {
       TH1F *neuBkg = (TH1F *)ig.get(label_neuBkg.Data())->getTH1F()->Clone();
       neuBkg->Sumw2();
       ProcessBookingBKG(ibooker, ig, "me0_strip_dg_neuBkg_rad", neuBkg, numEvts);
 
       delete neuBkg;
     }
     if (ig.get(label_totBkg.Data()) != nullptr) {
       TH1F *totBkg = (TH1F *)ig.get(label_totBkg.Data())->getTH1F()->Clone();
       totBkg->Sumw2();
       ProcessBookingBKG(ibooker, ig, "me0_strip_dg_totBkg_rad", totBkg, numEvts);
 
       delete totBkg;
     }
 
     delete numEvts;
   }
 }

void MuonME0DigisHarvestor::ProcessBooking	(	DQMStore::IBooker &	ibooker,
		DQMStore::IGetter &	ig,
		std::string	nameHist,
		TH1F *	num,
		TH1F *	den
	)

Definition at line 72 of file MuonME0DigisHarvestor.cc.

References DQMStore::IBooker::bookProfile(), ComputeEff(), PostProcessor_cff::profile, and fftjetcommon_cfi::title.

Referenced by dqmEndJob().

                                                                                              {
   if (num != nullptr && den != nullptr) {
     TProfile *profile = ComputeEff(num, den, nameHist);
 
     TString x_axis_title = TString(num->GetXaxis()->GetTitle());
     TString title = TString::Format("Digi Efficiency;%s;Eff.", x_axis_title.Data());
 
     profile->SetTitle(title.Data());
     ibooker.bookProfile(profile->GetName(), profile);
 
     delete profile;
 
   } else {
     edm::LogWarning("MuonME0DigisHarvestor") << "Can not find histograms";
     if (num == nullptr)
       edm::LogWarning("MuonME0DigisHarvestor") << "num not found";
     if (den == nullptr)
       edm::LogWarning("MuonME0DigisHarvestor") << "den not found";
   }
   return;
 }

void MuonME0DigisHarvestor::ProcessBookingBKG	(	DQMStore::IBooker &	ibooker,
		DQMStore::IGetter &	ig,
		std::string	nameHist,
		TH1F *	hist,
		TH1F *	hist2
	)

Definition at line 114 of file MuonME0DigisHarvestor.cc.

References DQMStore::IBooker::book1D(), ComputeBKG(), RPCpg::rate(), and fftjetcommon_cfi::title.

Referenced by dqmEndJob().

                                                                                                  {
   if (hist1 != nullptr && hist2 != nullptr) {
     TH1F *rate = ComputeBKG(hist1, hist2, nameHist);
 
     TString x_axis_title = TString(hist1->GetXaxis()->GetTitle());
     TString origTitle = TString(hist1->GetTitle());
     TString title = TString::Format((origTitle + ";%s;Rate [Hz/cm^{2}]").Data(), x_axis_title.Data());
 
     rate->SetTitle(title.Data());
     ibooker.book1D(rate->GetName(), rate);
 
   } else {
     edm::LogWarning("MuonME0DigisHarvestor") << "Can not find histograms";
     if (hist1 == nullptr)
       edm::LogWarning("MuonME0DigisHarvestor") << "num not found";
     if (hist2 == nullptr)
       edm::LogWarning("MuonME0DigisHarvestor") << "den not found";
   }
   return;
 }

Member Data Documentation

std::string MuonME0DigisHarvestor::dbe_path_

private

Definition at line 33 of file MuonME0DigisHarvestor.h.

Referenced by dqmEndJob(), and MuonME0DigisHarvestor().

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation