MuonME0DigisHarvestor Class Reference

#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)
virtual void	dqmEndJob (DQMStore::IBooker &, DQMStore::IGetter &)
	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 ()
	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 &	itemsToGetFromEvent () 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 MuonME0DigisHarvestor.h.

Constructor & Destructor Documentation

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

explicit

constructor

Geometry.

Log messages

Definition at line 34 of file MuonME0DigisHarvestor.cc.

References dbe_path_, and AlCaHLTBitMon_QueryRunRegistry::string.

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

MuonME0DigisHarvestor::~MuonME0DigisHarvestor

(

)

destructor

Definition at line 40 of file MuonME0DigisHarvestor.cc.

{
}

Member Function Documentation

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

Definition at line 105 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 45 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  )

virtual

Implements DQMEDHarvester.

Definition at line 157 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_.c_str());
 
  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 77 of file MuonME0DigisHarvestor.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("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 129 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().