#include <MuonME0SegHarvestor.h>

Inheritance diagram for MuonME0SegHarvestor:

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

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

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

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

	~MuonME0SegHarvestor ()
	destructor More...

Public Member Functions inherited from DQMEDHarvester
virtual void	analyze (edm::Event const &, edm::EventSetup const &) final

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

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

	DQMEDHarvester (void)

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

virtual void	endJob () final

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

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

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

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

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDAnalyzerBase ()

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

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

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

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

virtual	~EDConsumerBase () noexcept(false)

Private Attributes
std::string	dbe_path_

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 edm::one::EDAnalyzerBase
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 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)

Detailed Description

Definition at line 18 of file MuonME0SegHarvestor.h.

Constructor & Destructor Documentation

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

explicit

constructor

Geometry.

Log messages

Definition at line 34 of file MuonME0SegHarvestor.cc.

References dbe_path_, and AlCaHLTBitMon_QueryRunRegistry::string.

 {
   dbe_path_ = std::string("MuonME0RecHitsV/ME0SegmentsTask/");
 }

MuonME0SegHarvestor::~MuonME0SegHarvestor ( )

destructor

Definition at line 40 of file MuonME0SegHarvestor.cc.

41 {

42 }

Member Function Documentation

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

Definition at line 45 of file MuonME0SegHarvestor.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 = "Segment 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("MuonME0SegHarvestor")<<"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 MuonME0SegHarvestor::dqmEndJob	(	DQMStore::IBooker &	ibooker,
		DQMStore::IGetter &	ig
	)

virtual

Implements DQMEDHarvester.

Definition at line 106 of file MuonME0SegHarvestor.cc.

References dbe_path_, DEFINE_FWK_MODULE, DQMStore::IGetter::get(), ProcessBooking(), and DQMStore::IGetter::setCurrentFolder().

 {
   ig.setCurrentFolder(dbe_path_.c_str());
     
   TString eta_label_den = TString(dbe_path_)+"me0_simsegment_eta";
   TString eta_label_num = TString(dbe_path_)+"me0_matchedsimsegment_eta";
   TString pt_label_den = TString(dbe_path_)+"me0_simsegment_pt";
   TString pt_label_num = TString(dbe_path_)+"me0_matchedsimsegment_pt";
   TString phi_label_den = TString(dbe_path_)+"me0_simsegment_phi";
   TString phi_label_num = TString(dbe_path_)+"me0_matchedsimsegment_phi";
 
   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();
       
       ProcessBooking( ibooker, ig, "me0segment_eff_vs_eta", num_vs_eta, den_vs_eta );
       
       delete num_vs_eta;
       delete den_vs_eta;
       
   }
   else edm::LogWarning("MuonME0SegHarvestor")<<"Can not find histograms: "<<eta_label_num<<" or "<<eta_label_den;
  
   if( ig.get(pt_label_num.Data()) !=nullptr && ig.get(pt_label_den.Data()) !=nullptr ) {
         
       TH1F* num_vs_pt = (TH1F*)ig.get(pt_label_num.Data())->getTH1F()->Clone();
       num_vs_pt->Sumw2();
       TH1F *den_vs_pt = (TH1F*)ig.get(pt_label_den.Data())->getTH1F()->Clone();
       den_vs_pt->Sumw2();
       
       ProcessBooking( ibooker, ig, "me0segment_eff_vs_pt", num_vs_pt, den_vs_pt );
       
       delete num_vs_pt;
       delete den_vs_pt;
       
   }
   else edm::LogWarning("MuonME0SegHarvestor")<<"Can not find histograms: "<<pt_label_num<<" or "<<pt_label_den;
 
   if( ig.get(phi_label_num.Data()) !=nullptr && ig.get(phi_label_den.Data()) !=nullptr ) {
         
       TH1F* num_vs_phi = (TH1F*)ig.get(phi_label_num.Data())->getTH1F()->Clone();
       num_vs_phi->Sumw2();
       TH1F *den_vs_phi = (TH1F*)ig.get(phi_label_den.Data())->getTH1F()->Clone();
       den_vs_phi->Sumw2();
       
       ProcessBooking( ibooker, ig, "me0segment_eff_vs_phi", num_vs_phi, den_vs_phi );
       
       delete num_vs_phi;
       delete den_vs_phi;
       
   }
   else edm::LogWarning("MuonME0SegHarvestor")<<"Can not find histograms: "<<phi_label_num<<" or "<<phi_label_den;
 
     
 }

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

Definition at line 77 of file MuonME0SegHarvestor.cc.

References DQMStore::IBooker::bookProfile(), ComputeEff(), NewPostProcessor_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("Segment Efficiency;%s;Eff.",x_axis_title.Data());
 
     profile->SetTitle( title.Data());
     ibooker.bookProfile( profile->GetName(),profile);
       
     delete profile;
       
   }
   else {
       
     edm::LogWarning("MuonME0SegHarvestor")<<"Can not find histograms";
     if ( num == nullptr) edm::LogWarning("MuonME0SegHarvestor")<<"num not found";
     if ( den == nullptr) edm::LogWarning("MuonME0SegHarvestor")<<"den not found";
       
   }
   return;
     
 }

Member Data Documentation

std::string MuonME0SegHarvestor::dbe_path_

private

Definition at line 31 of file MuonME0SegHarvestor.h.

Referenced by dqmEndJob(), and MuonME0SegHarvestor().

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation