#include <MuonGEMDigisHarvestor.h>

Inheritance diagram for MuonGEMDigisHarvestor:

Public Member Functions
TProfile *	ComputeEff (TH1F num, TH1F denum)

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

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

void	ProcessBooking (DQMStore::IBooker &, DQMStore::IGetter &, const char label, TString suffix, TH1F track_hist, TH1F *sh_hist)

	~MuonGEMDigisHarvestor () 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	compareable_dbe_hist_prefix_

std::string	compareable_dbe_path_

std::string	dbe_hist_prefix_

std::string	dbe_path_

std::string	outputFile_

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 20 of file MuonGEMDigisHarvestor.h.

Constructor & Destructor Documentation

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

explicit

constructor

Data Format.

Geometry Log messages

Definition at line 48 of file MuonGEMDigisHarvestor.cc.

References compareable_dbe_hist_prefix_, compareable_dbe_path_, dbe_hist_prefix_, dbe_path_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), outputFile_, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                       {
   dbe_path_ = ps.getParameter<std::string>("dbePath");
   dbe_hist_prefix_ = ps.getParameter<std::string>("dbeHistPrefix");
   compareable_dbe_path_ = ps.getParameter<std::string>("compareDBEPath");
   compareable_dbe_hist_prefix_ = ps.getParameter<std::string>("compareDBEHistPrefix");
 
   outputFile_ = ps.getUntrackedParameter<std::string>("outputFile", "myfile.root");
 }

MuonGEMDigisHarvestor::~MuonGEMDigisHarvestor ( )

override

destructor

Definition at line 57 of file MuonGEMDigisHarvestor.cc.

57 {}

Member Function Documentation

TProfile * MuonGEMDigisHarvestor::ComputeEff	(	TH1F *	num,
		TH1F *	denum
	)

Definition at line 58 of file MuonGEMDigisHarvestor.cc.

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

Referenced by ProcessBooking().

                                                                   {
   std::string name = "eff_" + std::string(num->GetName());
   std::string title = "Eff. " + 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;
       LogDebug("MuonGEMDigisHarvestor") << "Alert! specific bin's num is bigger than denum";
     }
     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 MuonGEMDigisHarvestor::dqmEndJob	(	DQMStore::IBooker &	ibooker,
		DQMStore::IGetter &	ig
	)

overridevirtual

Implements DQMEDHarvester.

Definition at line 125 of file MuonGEMDigisHarvestor.cc.

References GEMDetLabel::c_suffix, dbe_path_, DEFINE_FWK_MODULE, DQMStore::IGetter::get(), mps_fire::i, gen::k, GEMDetLabel::l_suffix, LogDebug, ProcessBooking(), GEMDetLabel::s_suffix, DQMStore::IGetter::setCurrentFolder(), and createPayload::suffix.

                                                                                    {
   ig.setCurrentFolder(dbe_path_);
   TH1F *gem_trk_eta[3];
   TH1F *gem_trk_phi[3][2];
 
   TH1F *sh_eta[3][4];
   TH1F *sh_phi[3][4][3];
 
   // simplePlots
   /*
   for( int region = -1 ; region <= 1 ; region = region+2) {
     for ( int station = 0 ; station <2 ; station++) {
       if ( station ==1 ) station=2 ;
       TString dcEta_label =
   TString::Format("%s%s_r%d%s",dbe_path_.c_str(),dbe_hist_prefix_.c_str(),
   region, s_suffix[station].c_str()); TString denum_dcEta_label =
   TString::Format("%s%s_r%d%s",compareable_dbe_path_.c_str(),compareable_dbe_hist_prefix_.c_str(),
   region, s_suffix[station].c_str());
 
 
         if ( ig.get( dcEta_label.Data()) != nullptr && ig.get(
   denum_dcEta_label.Data()) != nullptr) { TH2F* dcEta = (TH2F*)ig.get(
   dcEta_label.Data())->getTH2F()->Clone(); TH2F* denum_dcEta = (TH2F*)ig.get(
   denum_dcEta_label.Data())->getTH2F()->Clone(); dcEta->Divide(denum_dcEta);
         TH2F* eff_dcEta = (TH2F*)dcEta->Clone();
 
         TString eff_dcEta_title = TString::Format("Hits Efficiency on detector
   component at r%d%s",region,s_suffix[station].c_str()); TString eff_dcEta_label
   = TString::Format("eff_DigiHit_r%d%s",region,s_suffix[station].c_str());
 
 
         eff_dcEta->SetName( eff_dcEta_label.Data());
         eff_dcEta->SetTitle( eff_dcEta_title.Data());
 
         ibooker.book2D(eff_dcEta->GetName(), eff_dcEta);
       }
       else {
         std::cout<<"Failed to get histograms"<<std::endl;
         std::cout<<dcEta_label<<std::endl;
         std::cout<<denum_dcEta_label<<std::endl;
       }
     }
   }
   */
 
   using namespace GEMDetLabel;
 
   // detailPlots
   for (unsigned int i = 0; i < s_suffix.size(); i++) {
     TString eta_label = TString(dbe_path_) + "track_eta" + s_suffix[i];
     TString phi_label;
     if (ig.get(eta_label.Data()) != nullptr) {
       gem_trk_eta[i] = (TH1F *)ig.get(eta_label.Data())->getTH1F()->Clone();
       gem_trk_eta[i]->Sumw2();
     } else
       LogDebug("MuonGEMDigisHarvestor") << "Can not found track_eta";
     for (unsigned int k = 0; k < c_suffix.size(); k++) {
       phi_label = TString(dbe_path_.c_str()) + "track_phi" + s_suffix[i] + c_suffix[k];
       if (ig.get(phi_label.Data()) != nullptr) {
         gem_trk_phi[i][k] = (TH1F *)ig.get(phi_label.Data())->getTH1F()->Clone();
         gem_trk_phi[i][k]->Sumw2();
       } else
         LogDebug("MuonGEMDigisHarvestor") << "Can not found track_phi";
     }
 
     if (ig.get(eta_label.Data()) != nullptr && ig.get(phi_label.Data()) != nullptr) {
       for (unsigned int j = 0; j < l_suffix.size(); j++) {
         TString suffix = TString(s_suffix[i]) + TString(l_suffix[j]);
         TString eta_label = TString(dbe_path_) + "dg_sh_eta" + suffix;
         if (ig.get(eta_label.Data()) != nullptr) {
           sh_eta[i][j] = (TH1F *)ig.get(eta_label.Data())->getTH1F()->Clone();
           sh_eta[i][j]->Sumw2();
         } else
           LogDebug("MuonGEMDigisHarvestor") << "Can not found eta histogram : " << eta_label;
         ProcessBooking(ibooker, ig, "dg_eta", suffix, gem_trk_eta[i], sh_eta[i][j]);
         ProcessBooking(ibooker, ig, "pad_eta", suffix, gem_trk_eta[i], sh_eta[i][j]);
         ProcessBooking(ibooker, ig, "copad_eta", suffix, gem_trk_eta[i], sh_eta[i][j]);
         for (unsigned int k = 0; k < c_suffix.size(); k++) {
           suffix = TString(s_suffix[i]) + TString(l_suffix[j]) + TString(c_suffix[k]);
           TString phi_label = TString(dbe_path_) + "dg_sh_phi" + suffix;
           if (ig.get(phi_label.Data()) != nullptr) {
             sh_phi[i][j][k] = (TH1F *)ig.get(phi_label.Data())->getTH1F()->Clone();
             sh_phi[i][j][k]->Sumw2();
           } else {
             LogDebug("MuonGEMDigisHarvestor") << "Can not found phi plots : " << phi_label;
             continue;
           }
           ProcessBooking(ibooker, ig, "dg_phi", suffix, gem_trk_phi[i][k], sh_phi[i][j][k]);
           ProcessBooking(ibooker, ig, "pad_phi", suffix, gem_trk_phi[i][k], sh_phi[i][j][k]);
           ProcessBooking(ibooker, ig, "copad_phi", suffix, gem_trk_phi[i][k], sh_phi[i][j][k]);
         }
       }
     } else
       LogDebug("MuonGEMDigisHarvestor") << "Can not find eta or phi of all track";
   }
 }

void MuonGEMDigisHarvestor::ProcessBooking	(	DQMStore::IBooker &	ibooker,
		DQMStore::IGetter &	ig,
		const char *	label,
		TString	suffix,
		TH1F *	track_hist,
		TH1F *	sh_hist
	)

Definition at line 91 of file MuonGEMDigisHarvestor.cc.

References DQMStore::IBooker::bookProfile(), ComputeEff(), dbe_path_, DQMStore::IGetter::get(), create_public_lumi_plots::hist, LogDebug, PostProcessor_cff::profile, AlCaHLTBitMon_QueryRunRegistry::string, createPayload::suffix, and runGCPTkAlMap::title.

Referenced by dqmEndJob().

                                                           {
   TString dbe_label = TString(dbe_path_) + label + suffix;
   if (ig.get(dbe_label.Data()) != nullptr && sh_hist != nullptr && track_hist != nullptr) {
     TH1F *hist = (TH1F *)ig.get(dbe_label.Data())->getTH1F()->Clone();
     TProfile *profile = ComputeEff(hist, track_hist);
     TProfile *profile_sh = ComputeEff(hist, sh_hist);
     profile_sh->SetName((profile->GetName() + std::string("_sh")).c_str());
     TString x_axis_title = TString(hist->GetXaxis()->GetTitle());
     TString title = TString::Format(
         "Eff. for a SimTrack to have an associated GEM digi in %s;%s;Eff.", suffix.Data(), x_axis_title.Data());
     TString title2 = TString::Format(
         "Eff. for a SimTrack to have an associated GEM digi in "
         "%s with a matched SimHit;%s;Eff.",
         suffix.Data(),
         x_axis_title.Data());
     profile->SetTitle(title.Data());
     profile_sh->SetTitle(title2.Data());
     ibooker.bookProfile(profile->GetName(), profile);
     ibooker.bookProfile(profile_sh->GetName(), profile_sh);
   } else {
     LogDebug("MuonGEMDigisHarvestor") << "Can not found histogram of " << dbe_label;
     if (track_hist == nullptr)
       LogDebug("MuonGEMDigisHarvestor") << "track not found";
     if (sh_hist == nullptr)
       LogDebug("MuonGEMDigisHarvestor") << "sh_hist not found";
   }
   return;
 }