PFClient_JetRes Class Reference

#include <PFClient_JetRes.h>

Inheritance diagram for PFClient_JetRes:

Public Member Functions
	PFClient_JetRes (const edm::ParameterSet &parameterSet)
Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzer ()
ModuleDescription const &	moduleDescription () const
std::string	workerType () const
virtual	~EDAnalyzer ()
Public Member Functions inherited from edm::EDConsumerBase
	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 (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

Private Member Functions
void	analyze (edm::Event const &, edm::EventSetup const &)
void	beginJob ()
void	createEfficiencyPlots (std::string &folder, std::string &name)
void	createResolutionPlots (std::string &folder, std::string &name)
void	doEfficiency ()
void	doSummaries ()
void	endJob ()
void	endRun (edm::Run const &run, edm::EventSetup const &eSetup)
void	getHistogramParameters (MonitorElement *me_slice, double &avarage, double &rms, double &mean, double &sigma)

Private Attributes
DQMStore *	dqmStore_
std::vector< std::string >	effHistogramNames_
bool	efficiencyFlag_
std::vector< std::string >	folderNames_
std::vector< std::string >	histogramNames_
std::vector< int >	PtBins_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)
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 12 of file PFClient_JetRes.h.

Constructor & Destructor Documentation

PFClient_JetRes::PFClient_JetRes

(

const edm::ParameterSet &

parameterSet

)

Definition at line 18 of file PFClient_JetRes.cc.

References effHistogramNames_, efficiencyFlag_, folderNames_, edm::ParameterSet::getParameter(), histogramNames_, and PtBins_.

{
  folderNames_ = parameterSet.getParameter< std::vector<std::string> >( "FolderNames" );
  histogramNames_  = parameterSet.getParameter< std::vector<std::string> >( "HistogramNames" );
  efficiencyFlag_ =  parameterSet.getParameter< bool> ("CreateEfficiencyPlots" );
  effHistogramNames_  = parameterSet.getParameter< std::vector<std::string> >( "HistogramNamesForEfficiencyPlots" );
  PtBins_ = parameterSet.getParameter< std::vector<int> >( "VariablePtBins" ) ;

}

Member Function Documentation

void PFClient_JetRes::analyze ( edm::Event const &  , edm::EventSetup const &    )

inlineprivatevirtual

Implements edm::EDAnalyzer.

Definition at line 19 of file PFClient_JetRes.h.

{;}

void PFClient_JetRes::beginJob ( void  )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 30 of file PFClient_JetRes.cc.

References dqmStore_, and cppFunctionSkipper::operator.

                               {

  dqmStore_ = edm::Service<DQMStore>().operator->();
}

void PFClient_JetRes::createEfficiencyPlots ( std::string &  folder, std::string &  name  )

private

Definition at line 200 of file PFClient_JetRes.cc.

References DQMStore::book1D(), MonitorElement::DQM_KIND_TH1F, dqmStore_, postValidation_cfi::efficiency, DQMStore::get(), MonitorElement::getAxisTitle(), MonitorElement::getBinContent(), MonitorElement::getNbinsX(), MonitorElement::getTH1F(), MonitorElement::kind(), MonitorElement::Reset(), MonitorElement::setBinContent(), DQMStore::setCurrentFolder(), MonitorElement::setEfficiencyFlag(), AlCaHLTBitMon_QueryRunRegistry::string, SiStripMonitorClusterAlca_cfi::xmax, and SiStripMonitorClusterAlca_cfi::xmin.

Referenced by doEfficiency().

                                                                              {     
  MonitorElement* me1 = dqmStore_->get(folder+"/"+name);
  MonitorElement* me2 = dqmStore_->get(folder+"/"+name+"ref_");
  if (!me1 || !me2) return;
  MonitorElement* me_eff;
  if ( (me1->kind() == MonitorElement::DQM_KIND_TH1F) &&
       (me1->kind() == MonitorElement::DQM_KIND_TH1F) ) {
    TH1* th1 = me1->getTH1F();
    size_t nbinx = me1->getNbinsX();
    
    float xmin = th1->GetXaxis()->GetXmin();
    float xmax = th1->GetXaxis()->GetXmax();
    std::string xtit = me1->getAxisTitle(1);
    std::string tit_new;
    tit_new = ";"+xtit+";Efficiency"; 

    dqmStore_->setCurrentFolder(folder);
    me_eff = dqmStore_->book1D("efficiency_"+name,tit_new, nbinx, xmin, xmax); 
                 
    double efficiency;
    me_eff->Reset(); me_eff->setEfficiencyFlag();
    for (size_t ix = 1; ix < nbinx+1; ++ix) {
      float val1 = me1->getBinContent(ix);
      float val2 = me2->getBinContent(ix);
      if (val2 > 0.0) efficiency = val1/val2;
      else efficiency = 0;   
      me_eff->setBinContent(ix,efficiency);
    }
  }
}

void PFClient_JetRes::createResolutionPlots ( std::string &  folder, std::string &  name  )

private

Definition at line 81 of file PFClient_JetRes.cc.

References PDRates::average, DQMStore::book1D(), MonitorElement::DQM_KIND_TH2D, MonitorElement::DQM_KIND_TH2F, MonitorElement::DQM_KIND_TH2S, dqmStore_, DQMStore::get(), getHistogramParameters(), MonitorElement::getNbinsY(), MonitorElement::getTH2F(), MonitorElement::kind(), timingPdfMaker::mean, PtBins_, plotscripts::rms(), MonitorElement::setBinContent(), MonitorElement::setBinError(), DQMStore::setCurrentFolder(), MonitorElement::setEfficiencyFlag(), MonitorElement::setEntries(), AlCaHLTBitMon_QueryRunRegistry::string, fw3dlego::xbins, SiStripMonitorClusterAlca_cfi::ymax, and SiStripMonitorClusterAlca_cfi::ymin.

Referenced by doSummaries().

                                                                              {     
  MonitorElement* me = dqmStore_->get(folder+"/"+name);
  if (!me) return;

  MonitorElement* pT[PtBins_.size() -1];
  std::vector<double> pTEntries(PtBins_.size()-1, 0) ;

  //std::vector<std::string> pTRange (PtBins_.size() -1) ;
  //char* pTRange[PtBins_.size() -1] ;
  TString pTRange[PtBins_.size() -1] ;
  //float pTCenter[PtBins_.size() -1] ;

  MonitorElement* me_average; 
  MonitorElement* me_rms; 
  MonitorElement* me_mean; 
  MonitorElement* me_sigma;

  if ( (me->kind() == MonitorElement::DQM_KIND_TH2F) ||
       (me->kind() == MonitorElement::DQM_KIND_TH2S) ||
       (me->kind() == MonitorElement::DQM_KIND_TH2D) ) {
    TH2* th = me->getTH2F();
    //size_t nbinx = me->getNbinsX();
    size_t nbinx = PtBins_.size() -1;
    size_t nbiny = me->getNbinsY();
    
    float ymin = th->GetYaxis()->GetXmin();
    float ymax = th->GetYaxis()->GetXmax();
    std::string xtit = th->GetXaxis()->GetTitle();
    //std::string ytit = th->GetYaxis()->GetTitle();
    std::string ytit = "#Deltap_{T}/p_{T}";

    float* xbins = new float[nbinx+1];
    for (size_t ix = 1; ix < nbinx+1; ++ix) {
      //xbins[ix-1] = th->GetBinLowEdge(ix);
       xbins[ix-1] = PtBins_[ix-1];
       //snprintf(pTRange[ix-1].data(), 15, "Pt%d_%d", PtBins_[ix-1], PtBins_[ix]);
       pTRange[ix-1] = TString::Format("Pt%d_%d", PtBins_[ix-1], PtBins_[ix]) ;
       if (name == "BRdelta_et_Over_et_VS_et_") pTRange[ix-1] = TString::Format("BRPt%d_%d", PtBins_[ix-1], PtBins_[ix]) ;
       else if (name == "ERdelta_et_Over_et_VS_et_") pTRange[ix-1] = TString::Format("ERPt%d_%d", PtBins_[ix-1], PtBins_[ix]) ;

       //pTCenter[ix-1] = (PtBins_[ix] - PtBins_[ix-1]) / 2. ;
       if (ix == nbinx) {
     //xbins[ix] = th->GetXaxis()->GetBinUpEdge(ix);
     xbins[ix] = PtBins_[ix];
       }
    }    

    std::string tit_new;
    dqmStore_->setCurrentFolder(folder);
    //MonitorElement* me_slice = dqmStore_->book1D("PFlowSlice", "PFlowSlice", nbiny, ymin, ymax); 
    
    tit_new = "Average "+ytit+";"+xtit+";Average_"+ytit ; 
    me_average = dqmStore_->book1D("average_"+name,tit_new, nbinx, xbins); 
    me_average->setEfficiencyFlag();
    tit_new = "RMS "+ytit+";"+xtit+";RMS_"+ytit ; 
    me_rms     = dqmStore_->book1D("rms_"+name,tit_new, nbinx, xbins); 
    me_rms->setEfficiencyFlag();
    tit_new = ";"+xtit+";Mean_"+ytit; 
    me_mean    = dqmStore_->book1D("mean_"+name,tit_new, nbinx, xbins); 
    me_mean->setEfficiencyFlag();
    tit_new = ";"+xtit+";Sigma_"+ytit;               
    me_sigma   = dqmStore_->book1D("sigma_"+name,tit_new, nbinx, xbins); 
    me_sigma->setEfficiencyFlag();
    
    double  average, rms, mean, sigma;
    for (size_t ix = 1; ix < nbinx+1; ++ix) {
      //me_slice->Reset();
      if (name == "delta_et_Over_et_VS_et_") pT[ix-1] = dqmStore_->book1D(pTRange[ix-1], TString::Format("Total %s;%s;Events", ytit.data(), ytit.data() ), nbiny, ymin, ymax) ;
      if (name == "BRdelta_et_Over_et_VS_et_") pT[ix-1] = dqmStore_->book1D(pTRange[ix-1], TString::Format("Barrel %s;%s;Events", ytit.data(), ytit.data()), nbiny, ymin, ymax) ;
      else if (name == "ERdelta_et_Over_et_VS_et_") pT[ix-1] = dqmStore_->book1D(pTRange[ix-1], TString::Format("Endcap %s;%s;Events", ytit.data(), ytit.data() ), nbiny, ymin, ymax) ;

      for (size_t iy = 0; iy <= nbiny+1; ++iy)  // add under and overflow
    if (th->GetBinContent(ix,iy)) {
      //me_slice->setBinContent(iy,th->GetBinContent(ix,iy)); 
      pT[ix-1]->setBinContent(iy, th->GetBinContent(ix,iy));
      pT[ix-1]->setBinError(iy, th->GetBinError(ix,iy));
      pTEntries[ix-1] += th->GetBinContent(ix,iy); }
      
      pT[ix-1]->setEntries( pTEntries[ix-1] ) ;

      //getHistogramParameters(me_slice, average, rms, mean, sigma);
      getHistogramParameters(pT[ix-1], average, rms, mean, sigma);
      me_average->setBinContent(ix,average);
      me_rms->setBinContent(ix,rms);
      me_mean->setBinContent(ix,mean);
      me_sigma->setBinContent(ix,sigma);
    }
    //if (me_slice) dqmStore_->removeElement(me_slice->getName());
    delete [] xbins;
  }
}

void PFClient_JetRes::doEfficiency ( )

private

Definition at line 66 of file PFClient_JetRes.cc.

References createEfficiencyPlots(), effHistogramNames_, folderNames_, cmsHarvester::path, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by endRun().

                                   {
  for (std::vector<std::string>::const_iterator ifolder = folderNames_.begin();
                                    ifolder != folderNames_.end(); ifolder++) {
    std::string path  = "ParticleFlow/"+(*ifolder);

    for (std::vector<std::string>::const_iterator ihist = effHistogramNames_.begin();
     ihist != effHistogramNames_.end(); ihist++) {
      std::string hname = (*ihist);
      createEfficiencyPlots(path, hname);
    }
  }
}

void PFClient_JetRes::doSummaries ( )

private

Definition at line 50 of file PFClient_JetRes.cc.

References createResolutionPlots(), folderNames_, histogramNames_, cmsHarvester::path, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by endRun().

                                  {

  for (std::vector<std::string>::const_iterator ifolder = folderNames_.begin();
                                    ifolder != folderNames_.end(); ifolder++) {
    std::string path  = "ParticleFlow/"+(*ifolder);

    for (std::vector<std::string>::const_iterator ihist = histogramNames_.begin();
       ihist != histogramNames_.end(); ihist++) {
      std::string hname = (*ihist); 
      createResolutionPlots(path, hname);
    }
  }
}

void PFClient_JetRes::endJob ( void  )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 44 of file PFClient_JetRes.cc.

                             {

}

void PFClient_JetRes::endRun ( edm::Run const &  run, edm::EventSetup const &  eSetup  )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 37 of file PFClient_JetRes.cc.

References doEfficiency(), doSummaries(), and efficiencyFlag_.

                                                                           {
  doSummaries();
  if (efficiencyFlag_) doEfficiency();
}

void PFClient_JetRes::getHistogramParameters ( MonitorElement *  me_slice, double &  avarage, double &  rms, double &  mean, double &  sigma  )

private

Definition at line 175 of file PFClient_JetRes.cc.

References MonitorElement::DQM_KIND_TH1F, MonitorElement::getMean(), MonitorElement::getRMS(), MonitorElement::getTH1F(), and MonitorElement::kind().

Referenced by createResolutionPlots().

                                                                  {
  average = 0.0;
  rms     = 0.0;
  mean    = 0.0;
  sigma   = 0.0;

  if (!me_slice) return;
  if  (me_slice->kind() == MonitorElement::DQM_KIND_TH1F) {
    average = me_slice->getMean();
    rms     = me_slice->getRMS();  
    TH1F* th_slice = me_slice->getTH1F();
    if (th_slice && th_slice->GetEntries() > 0) {
      th_slice->Fit( "gaus","Q0");
      TF1* gaus = th_slice->GetFunction( "gaus" );
      if (gaus) {
    sigma = gaus->GetParameter(2);
        mean  = gaus->GetParameter(1);
      }
    }
  }
}

Member Data Documentation

DQMStore* PFClient_JetRes::dqmStore_

private

Definition at line 37 of file PFClient_JetRes.h.

Referenced by beginJob(), createEfficiencyPlots(), and createResolutionPlots().

std::vector<std::string> PFClient_JetRes::effHistogramNames_

private

Definition at line 32 of file PFClient_JetRes.h.

Referenced by doEfficiency(), and PFClient_JetRes().

bool PFClient_JetRes::efficiencyFlag_

private

Definition at line 35 of file PFClient_JetRes.h.

Referenced by endRun(), and PFClient_JetRes().

std::vector<std::string> PFClient_JetRes::folderNames_

private

Definition at line 30 of file PFClient_JetRes.h.

Referenced by doEfficiency(), doSummaries(), and PFClient_JetRes().

std::vector<std::string> PFClient_JetRes::histogramNames_

private

Definition at line 31 of file PFClient_JetRes.h.

Referenced by doSummaries(), and PFClient_JetRes().

std::vector<int> PFClient_JetRes::PtBins_

private

Definition at line 33 of file PFClient_JetRes.h.

Referenced by createResolutionPlots(), and PFClient_JetRes().