MuonGEMDigisHarvestor Class Reference

#include <MuonGEMDigisHarvestor.h>

Inheritance diagram for MuonGEMDigisHarvestor:

Public Member Functions
TProfile *	ComputeEff (TH1F *num, TH1F *denum)
virtual void	dqmEndJob (DQMStore::IBooker &, DQMStore::IGetter &)
	MuonGEMDigisHarvestor (const edm::ParameterSet &)
	constructor More...
void	ProcessBooking (DQMStore::IBooker &, DQMStore::IGetter &, const char *label, TString suffix, TH1F *track_hist, TH1F *sh_hist)
virtual	~MuonGEMDigisHarvestor ()
	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 ()
ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesDependentUpon (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &oModuleLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

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

Constructor & Destructor Documentation

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

explicit

constructor

Definition at line 50 of file MuonGEMDigisHarvestor.cc.

References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), 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 ( )

virtual

destructor

Definition at line 61 of file MuonGEMDigisHarvestor.cc.

{
}

Member Function Documentation

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

Definition at line 64 of file MuonGEMDigisHarvestor.cc.

References i, LogDebug, mergeVDriftHistosByStation::name, mathSSE::sqrt(), AlCaHLTBitMon_QueryRunRegistry::string, groupFilesInBlocks::temp, and indexGen::title.

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

virtual

Implements DQMEDHarvester.

Definition at line 121 of file MuonGEMDigisHarvestor.cc.

References GEMDetLabel::c_suffix, DQMStore::IGetter::get(), getTH1F(), i, j, relval_2017::k, GEMDetLabel::l_suffix, LogDebug, GEMDetLabel::s_suffix, DQMStore::IGetter::setCurrentFolder(), and createPayload::suffix.

{
  ig.setCurrentFolder(dbe_path_.c_str());
  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;
      }
    }
  }
  */

  // detailPlots
  for( int i = 0 ; i < 3 ; 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 ( int k=0 ; k <3 ; 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( int j = 0; j < 4 ; 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 ( int k= 0 ; k< 3 ; 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 95 of file MuonGEMDigisHarvestor.cc.

References DQMStore::IBooker::bookProfile(), DQMStore::IGetter::get(), getTH1F(), estimatePileup::hist, LogDebug, AlCaHLTBitMon_QueryRunRegistry::string, createPayload::suffix, and indexGen::title.

{
  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;
}

Member Data Documentation

std::string MuonGEMDigisHarvestor::compareable_dbe_hist_prefix_

private

Definition at line 34 of file MuonGEMDigisHarvestor.h.

std::string MuonGEMDigisHarvestor::compareable_dbe_path_

private

Definition at line 34 of file MuonGEMDigisHarvestor.h.

std::string MuonGEMDigisHarvestor::dbe_hist_prefix_

private

Definition at line 34 of file MuonGEMDigisHarvestor.h.

std::string MuonGEMDigisHarvestor::dbe_path_

private

Definition at line 34 of file MuonGEMDigisHarvestor.h.

std::string MuonGEMDigisHarvestor::outputFile_

private

Definition at line 34 of file MuonGEMDigisHarvestor.h.