PhotonOfflineClient Class Reference

#include <PhotonOfflineClient.h>

Inheritance diagram for PhotonOfflineClient:

Public Member Functions
void	book2DHistoVector (DQMStore::IBooker &iBooker, std::vector< std::vector< MonitorElement * > > &vecOfHist, std::string histoType, std::string histoName, std::string title, int xbin, double xmin, double xmax, int ybin=1, double ymin=1, double ymax=2)
void	book3DHistoVector (DQMStore::IBooker &iBooker, std::vector< std::vector< std::vector< MonitorElement * > > > &vecOfHist, std::string histoType, std::string histoName, std::string title, int xbin, double xmin, double xmax, int ybin=1, double ymin=1, double ymax=2)
MonitorElement *	bookHisto (DQMStore::IBooker &iBooker, std::string histoName, std::string title, int bin, double min, double max)
void	dqmEndJob (DQMStore::IBooker &, DQMStore::IGetter &) override
	PhotonOfflineClient (const edm::ParameterSet &pset)
MonitorElement *	retrieveHisto (DQMStore::IGetter &iGetter, std::string dir, std::string name)
virtual void	runClient (DQMStore::IBooker &iBooker, DQMStore::IGetter &iGetter)
	~PhotonOfflineClient () override
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
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)
virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	dividePlots (MonitorElement *dividend, MonitorElement *numerator, MonitorElement *denominator)
void	dividePlots (MonitorElement *dividend, MonitorElement *numerator, double denominator)

Private Attributes
std::string	analyzerName_
bool	batch_
std::stringstream	currentFolder_
double	cutStep_
DQMStore *	dbe_
int	etaBin
double	etaMax
double	etaMin
int	etBin
double	etMax
double	etMin
bool	excludeBkgHistos_
int	histo_index_conversions_
int	histo_index_efficiency_
int	histo_index_invMass_
int	histo_index_photons_
std::string	inputFileName_
bool	minimalSetOfHistos_
int	numberOfSteps_
std::string	outputFileName_
std::vector< std::vector< MonitorElement * > >	p_badChannelsFractionVsEt_
std::vector< std::vector< MonitorElement * > >	p_badChannelsFractionVsEta_
std::vector< std::vector< MonitorElement * > >	p_badChannelsFractionVsPhi_
std::vector< std::vector< MonitorElement * > >	p_convFractionVsEt_
std::vector< std::vector< MonitorElement * > >	p_convFractionVsEta_
MonitorElement *	p_convFractionVsEtaLoose_
MonitorElement *	p_convFractionVsEtaTight_
MonitorElement *	p_convFractionVsEtLoose_
MonitorElement *	p_convFractionVsEtTight_
std::vector< std::vector< std::vector< MonitorElement * > > >	p_convFractionVsPhi_
MonitorElement *	p_efficiencyVsEtaHLT_
MonitorElement *	p_efficiencyVsEtaLoose_
MonitorElement *	p_efficiencyVsEtaTight_
MonitorElement *	p_efficiencyVsEtHLT_
MonitorElement *	p_efficiencyVsEtLoose_
MonitorElement *	p_efficiencyVsEtTight_
MonitorElement *	p_vertexReconstructionEfficiencyVsEta_
edm::ParameterSet	parameters_
std::vector< std::string >	parts_
int	phiBin
double	phiMax
double	phiMin
bool	standAlone_
std::vector< std::string >	types_
int	verbosity_

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

$Id: PhotonOfflineClient authors: Nancy Marinelli, U. of Notre Dame, US Jamie Antonelli, U. of Notre Dame, US

Definition at line 64 of file PhotonOfflineClient.h.

Constructor & Destructor Documentation

PhotonOfflineClient::PhotonOfflineClient ( const edm::ParameterSet &  pset )

explicit

Definition at line 25 of file PhotonOfflineClient.cc.

References conversionPostprocessing_cfi::etaBin, ALCARECOTkAlBeamHalo_cff::etaMax, ALCARECOTkAlBeamHalo_cff::etaMin, conversionPostprocessing_cfi::etBin, conversionPostprocessing_cfi::etMax, btagDijet_cfi::etMin, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), conversionPostprocessing_cfi::phiBin, AlignmentTrackSelector_cfi::phiMax, AlignmentTrackSelector_cfi::phiMin, and muonDTDigis_cfi::pset.

{
  
  //dbe_ = 0;
  //dbe_ = edm::Service<DQMStore>().operator->();
  // dbe_->setVerbose(0);
  parameters_ = pset;

  analyzerName_    = pset.getParameter<string>("analyzerName");
  cutStep_       = pset.getParameter<double>("cutStep");
  numberOfSteps_ = pset.getParameter<int>("numberOfSteps");

  etMin  = pset.getParameter<double>("etMin");
  etMax  = pset.getParameter<double>("etMax");
  etBin  = pset.getParameter<int>("etBin");
  etaMin = pset.getParameter<double>("etaMin");
  etaMax = pset.getParameter<double>("etaMax");
  etaBin = pset.getParameter<int>("etaBin");
  phiMin = pset.getParameter<double>("phiMin");
  phiMax = pset.getParameter<double>("phiMax");
  phiBin = pset.getParameter<int>("phiBin");

  standAlone_ = pset.getParameter<bool>("standAlone");
  batch_ = pset.getParameter<bool>("batch");
  excludeBkgHistos_ =  pset.getParameter<bool>("excludeBkgHistos");

  outputFileName_ = pset.getParameter<string>("OutputFileName");
  inputFileName_  = pset.getUntrackedParameter<string>("InputFileName");

  histo_index_photons_ = 0;
  histo_index_conversions_ = 0;
  histo_index_efficiency_ = 0;
  histo_index_invMass_ = 0;

  types_.push_back("All");
  types_.push_back("GoodCandidate");
  if (!excludeBkgHistos_)  types_.push_back("Background");

  parts_.push_back("AllEcal");
  parts_.push_back("Barrel");
  parts_.push_back("Endcaps");

}

PhotonOfflineClient::~PhotonOfflineClient ( )

override

Definition at line 69 of file PhotonOfflineClient.cc.

{}

Member Function Documentation

void PhotonOfflineClient::book2DHistoVector	(	DQMStore::IBooker &	iBooker,
		std::vector< std::vector< MonitorElement * > > &	vecOfHist,
		std::string	histoType,
		std::string	histoName,
		std::string	title,
		int	xbin,
		double	xmin,
		double	xmax,
		int	ybin = 1,
		double	ymin = 1,
		double	ymax = 2
	)

Definition at line 359 of file PhotonOfflineClient.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), DQMStore::IBooker::bookProfile(), TkAlMuonSelectors_cfi::cut, callgraph::kind, DQMStore::IBooker::pwd(), DQMStore::IBooker::setCurrentFolder(), fftjetcommon_cfi::title, parallelization::uint(), TrackerOfflineValidation_Dqm_cff::xmax, TrackerOfflineValidation_Dqm_cff::xmin, Phase2TrackerMonitorDigi_cff::ymax, and Phase2TrackerMonitorDigi_cff::ymin.

{
  int histo_index = 0;

  vector<MonitorElement*> temp1DVector;

  //determining which folder we're in
  bool conversionPlot = false;
  if(iBooker.pwd().find( "Conversions" ) != string::npos) conversionPlot = true;


  if(conversionPlot){
    histo_index_conversions_++;
    histo_index = histo_index_conversions_;
  }
  else{
    histo_index_photons_++;
    histo_index = histo_index_photons_;
  }


  stringstream histo_number_stream;
  histo_number_stream << "h_";
  if(histo_index<10)   histo_number_stream << "0";
  histo_number_stream << histo_index << "_";


  for(int cut = 0; cut != numberOfSteps_; ++cut){ //looping over Et cut values

    for(uint type=0;type!=types_.size();++type){  //looping over isolation type

      currentFolder_.str("");
      currentFolder_ << "Egamma/"+analyzerName_+"/" << types_[type] << "Photons/Et above " << (cut+1)*cutStep_ << " GeV";
      if(conversionPlot) currentFolder_ << "/Conversions";

      iBooker.setCurrentFolder(currentFolder_.str());

      string kind;
      if(conversionPlot) kind = " Conversions: ";
      else kind = " Photons: ";

      if(histoType=="1D")           temp1DVector.push_back(iBooker.book1D(histo_number_stream.str()+histoName,types_[type]+kind+title,xbin,xmin,xmax));
      else if(histoType=="2D")      temp1DVector.push_back(iBooker.book2D(histo_number_stream.str()+histoName,types_[type]+kind+title,xbin,xmin,xmax,ybin,ymin,ymax));
      else if(histoType=="Profile") temp1DVector.push_back(iBooker.bookProfile(histo_number_stream.str()+histoName,types_[type]+kind+title,xbin,xmin,xmax,ybin,ymin,ymax,""));
      //else cout << "bad histoType\n";
    }

    temp2DVector.push_back(temp1DVector);
    temp1DVector.clear();
  }



}

void PhotonOfflineClient::book3DHistoVector	(	DQMStore::IBooker &	iBooker,
		std::vector< std::vector< std::vector< MonitorElement * > > > &	vecOfHist,
		std::string	histoType,
		std::string	histoName,
		std::string	title,
		int	xbin,
		double	xmin,
		double	xmax,
		int	ybin = 1,
		double	ymin = 1,
		double	ymax = 2
	)

Definition at line 420 of file PhotonOfflineClient.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), DQMStore::IBooker::bookProfile(), TkAlMuonSelectors_cfi::cut, callgraph::kind, DQMStore::IBooker::pwd(), DQMStore::IBooker::setCurrentFolder(), fftjetcommon_cfi::title, parallelization::uint(), TrackerOfflineValidation_Dqm_cff::xmax, TrackerOfflineValidation_Dqm_cff::xmin, Phase2TrackerMonitorDigi_cff::ymax, and Phase2TrackerMonitorDigi_cff::ymin.

{
  int histo_index = 0;

  vector<MonitorElement*> temp1DVector;
  vector<vector<MonitorElement*> > temp2DVector;

  //determining which folder we're in
  bool conversionPlot = false;
  if(iBooker.pwd().find( "Conversions" ) != string::npos) conversionPlot = true;


  if(conversionPlot){
    histo_index_conversions_++;
    histo_index = histo_index_conversions_;
  }
  else{
    histo_index_photons_++;
    histo_index = histo_index_photons_;
  }

  stringstream histo_number_stream;
  histo_number_stream << "h_";
  if(histo_index<10)   histo_number_stream << "0";
  histo_number_stream << histo_index << "_";


  for(int cut = 0; cut != numberOfSteps_; ++cut){     //looping over Et cut values

    for(uint type=0;type!=types_.size();++type){      //looping over isolation type

      for(uint part=0;part!=parts_.size();++part){    //looping over different parts of the ecal

    currentFolder_.str("");
    currentFolder_ << "Egamma/"+analyzerName_+"/" << types_[type] << "Photons/Et above " << (cut+1)*cutStep_ << " GeV";
    if(conversionPlot) currentFolder_ << "/Conversions";

    iBooker.setCurrentFolder(currentFolder_.str());

    string kind;
    if(conversionPlot) kind = " Conversions: ";
    else kind = " Photons: ";

    if(histoType=="1D")           temp1DVector.push_back(iBooker.book1D(      histo_number_stream.str()+histoName+parts_[part],types_[type]+kind+parts_[part]+": "+title,xbin,xmin,xmax));
    else if(histoType=="2D")      temp1DVector.push_back(iBooker.book2D(      histo_number_stream.str()+histoName+parts_[part],types_[type]+kind+parts_[part]+": "+title,xbin,xmin,xmax,ybin,ymin,ymax));
    else if(histoType=="Profile") temp1DVector.push_back(iBooker.bookProfile( histo_number_stream.str()+histoName+parts_[part],types_[type]+kind+parts_[part]+": "+title,xbin,xmin,xmax,ybin,ymin,ymax,""));
    //else cout << "bad histoType\n";


      }

      temp2DVector.push_back(temp1DVector);
      temp1DVector.clear();
    }

    temp3DVector.push_back(temp2DVector);
    temp2DVector.clear();
  }


}

MonitorElement * PhotonOfflineClient::bookHisto	(	DQMStore::IBooker &	iBooker,
		std::string	histoName,
		std::string	title,
		int	bin,
		double	min,
		double	max
	)

Definition at line 333 of file PhotonOfflineClient.cc.

References stringResolutionProvider_cfi::bin, DQMStore::IBooker::book1D(), SiStripPI::max, min(), DQMStore::IBooker::pwd(), and fftjetcommon_cfi::title.

{

  int histo_index = 0;
  stringstream histo_number_stream;


  //determining which folder we're in
  if(iBooker.pwd().find( "InvMass" ) != string::npos){
    histo_index_invMass_++;
    histo_index = histo_index_invMass_;
  }
  if(iBooker.pwd().find( "Efficiencies" ) != string::npos){
    histo_index_efficiency_++;
    histo_index = histo_index_efficiency_;
  }
  histo_number_stream << "h_";
  if(histo_index<10)   histo_number_stream << "0";
  histo_number_stream << histo_index;

  return iBooker.book1D(histo_number_stream.str()+"_"+histoName,title,bin,min,max);

}

void PhotonOfflineClient::dividePlots ( MonitorElement *  dividend, MonitorElement *  numerator, MonitorElement *  denominator  )

private

Definition at line 293 of file PhotonOfflineClient.cc.

References MonitorElement::getBinContent(), MonitorElement::getEntries(), MonitorElement::getNbinsX(), MonitorElement::setBinContent(), MonitorElement::setBinError(), MonitorElement::setEfficiencyFlag(), MonitorElement::setEntries(), mathSSE::sqrt(), and relativeConstraints::value.

                                                                                                                      {
  double value,err;

  dividend->setEfficiencyFlag(); 
  if(denominator->getEntries()==0) return;

  for (int j=1; j<=numerator->getNbinsX(); j++){
    if (denominator->getBinContent(j)!=0){
      value = ((double) numerator->getBinContent(j))/((double) denominator->getBinContent(j));
      err = sqrt( value*(1-value) / ((double) denominator->getBinContent(j)) );
      dividend->setBinContent(j, value);
      dividend->setBinError(j,err);
    }
    else {
      dividend->setBinContent(j, 0);
      dividend->setBinError(j,0);
    }
    dividend->setEntries(numerator->getEntries());
  }
}

void PhotonOfflineClient::dividePlots ( MonitorElement *  dividend, MonitorElement *  numerator, double  denominator  )

private

Definition at line 315 of file PhotonOfflineClient.cc.

References MonitorElement::getBinContent(), MonitorElement::getNbinsX(), MonitorElement::setBinContent(), MonitorElement::setBinError(), MonitorElement::setEfficiencyFlag(), mathSSE::sqrt(), and relativeConstraints::value.

                                                                                                             {
  double value,err;

  dividend->setEfficiencyFlag(); 
  for (int j=1; j<=numerator->getNbinsX(); j++){
    if (denominator!=0){
      value = ((double) numerator->getBinContent(j))/denominator;
      err = sqrt( value*(1-value) / denominator);
      dividend->setBinContent(j, value);
      dividend->setBinError(j,err);
    }
    else {
      dividend->setBinContent(j, 0);
    }
  }

}

void PhotonOfflineClient::dqmEndJob ( DQMStore::IBooker &  iBooker, DQMStore::IGetter &  iGetter  )

overridevirtual

Implements DQMEDHarvester.

Definition at line 74 of file PhotonOfflineClient.cc.

{
  if(!standAlone_) runClient(iBooker, iGetter);
}

MonitorElement * PhotonOfflineClient::retrieveHisto	(	DQMStore::IGetter &	iGetter,
		std::string	dir,
		std::string	name
	)

Definition at line 486 of file PhotonOfflineClient.cc.

References DQMStore::IGetter::getContents(), and parallelization::uint().

                                                                                                  {
  //cout << "dir = " << dir << endl;
  //cout << "name = " << name << endl;
  vector<MonitorElement*> histoVector;
  uint indexOfRelevantHistogram=0;
  string fullMEName = "";
  histoVector = iGetter.getContents(dir);
  for(uint index=0;index!=histoVector.size();index++){
    string MEName = histoVector[index]->getName();
    if(MEName.find( name ) != string::npos){
      indexOfRelevantHistogram = index;
      break;
    }
  }
  return histoVector[indexOfRelevantHistogram];
}

void PhotonOfflineClient::runClient ( DQMStore::IBooker &  iBooker, DQMStore::IGetter &  iGetter  )

virtual

Definition at line 80 of file PhotonOfflineClient.cc.

References gather_cfg::cout, TkAlMuonSelectors_cfi::cut, pfDeepCMVADiscriminatorsJetTags_cfi::denominator, DQMStore::IGetter::dirExists(), conversionPostprocessing_cfi::etaBin, ALCARECOTkAlBeamHalo_cff::etaMax, ALCARECOTkAlBeamHalo_cff::etaMin, conversionPostprocessing_cfi::etBin, conversionPostprocessing_cfi::etMax, btagDijet_cfi::etMin, DQMStore::IGetter::get(), MonitorElement::getIntValue(), pfDeepCMVADiscriminatorsJetTags_cfi::numerator, conversionPostprocessing_cfi::phiBin, AlignmentTrackSelector_cfi::phiMax, AlignmentTrackSelector_cfi::phiMin, DQMStore::IGetter::removeElement(), DQMStore::IGetter::setCurrentFolder(), and parallelization::uint().

{
  //  if(!dbe_) return;

  //if(batch_)  dbe_->open(inputFileName_);

  //std::cout << " PostProcessing analyzer name " << analyzerName_ << std::endl;
  if(!iGetter.dirExists("Egamma/"+analyzerName_)){
    std::cout << "Folder Egamma/"+analyzerName_+" does not exist - Abort the efficiency calculation "  << std::endl;
    return;
  }

  //find out how many histograms are in the various folders
  histo_index_photons_     = iGetter.get("Egamma/"+analyzerName_+"/numberOfHistogramsInPhotonsFolder")->getIntValue();
  histo_index_conversions_ = iGetter.get("Egamma/"+analyzerName_+"/numberOfHistogramsInConversionsFolder")->getIntValue();
  histo_index_efficiency_  = iGetter.get("Egamma/"+analyzerName_+"/numberOfHistogramsInEfficiencyFolder")->getIntValue();
  histo_index_invMass_     = iGetter.get("Egamma/"+analyzerName_+"/numberOfHistogramsInInvMassFolder")->getIntValue();

  iGetter.setCurrentFolder("Egamma/"+analyzerName_+"/");
  iGetter.removeElement("numberOfHistogramsInPhotonsFolder");
  iGetter.removeElement("numberOfHistogramsInConversionsFolder");
  iGetter.removeElement("numberOfHistogramsInEfficiencyFolder");
  iGetter.removeElement("numberOfHistogramsInInvMassFolder");

  string AllPath    = "Egamma/"+analyzerName_+"/AllPhotons/";
  string IsoPath    = "Egamma/"+analyzerName_+"/GoodCandidatePhotons/";
  string NonisoPath = "Egamma/"+analyzerName_+"/BackgroundPhotons/";
  string EffPath    = "Egamma/"+analyzerName_+"/Efficiencies/";

  //booking efficiency histograms
  iGetter.setCurrentFolder(EffPath);

  p_efficiencyVsEtaLoose_ = bookHisto(iBooker,"EfficiencyVsEtaLoose","Fraction of Photons passing Loose Isolation vs #eta;#eta",etaBin,etaMin, etaMax);
  p_efficiencyVsEtLoose_  = bookHisto(iBooker,"EfficiencyVsEtLoose", "Fraction of Photons passing Loose Isolation vs E_{T};E_{T} (GeV)",etBin,etMin,etMax);
  p_efficiencyVsEtaTight_ = bookHisto(iBooker,"EfficiencyVsEtaTight","Fraction of Photons passing Tight Isolation vs #eta;#eta",etaBin,etaMin, etaMax);
  p_efficiencyVsEtTight_  = bookHisto(iBooker,"EfficiencyVsEtTight", "Fraction of Photons passing Tight Isolation vs E_{T};E_{T} (GeV)",etBin,etMin,etMax);

  p_efficiencyVsEtaHLT_ = bookHisto(iBooker,"EfficiencyVsEtaHLT","Fraction of Photons firing HLT vs #eta;#eta",etaBin,etaMin, etaMax);
  p_efficiencyVsEtHLT_  = bookHisto(iBooker,"EfficiencyVsEtHLT", "Fraction of Photons firing HLT vs E_{T};E_{T} (GeV)",etBin,etMin,etMax);

  p_convFractionVsEtaLoose_ = bookHisto(iBooker,"ConvFractionVsEtaLoose","Fraction of Loosely Isolated Photons which are matched to two tracks vs #eta;#eta",etaBin,etaMin, etaMax);
  p_convFractionVsEtLoose_  = bookHisto(iBooker,"ConvFractionVsEtLoose", "Fraction of Loosely Isolated Photons which are matched to two tracks vs E_{T};E_{T} (GeV)",etBin,etMin,etMax);
  p_convFractionVsEtaTight_ = bookHisto(iBooker,"ConvFractionVsEtaTight","Fraction of Tightly Isolated Photons which are matched to two tracks vs #eta;#eta",etaBin,etaMin, etaMax);
  p_convFractionVsEtTight_  = bookHisto(iBooker,"ConvFractionVsEtTight", "Fraction of Tightly Isolated Photons which are matched to two tracks vs E_{T};E_{T} (GeV)",etBin,etMin,etMax);

  p_vertexReconstructionEfficiencyVsEta_ = bookHisto(iBooker,"VertexReconstructionEfficiencyVsEta","Fraction of Converted Photons which have a valid vertex vs #eta;#eta",etaBin,etaMin, etaMax);

  //booking conversion fraction histograms
  iGetter.setCurrentFolder(AllPath+"Et above 20 GeV/Conversions");
  book2DHistoVector(iBooker, p_convFractionVsEt_, "1D","convFractionVsEt", "Fraction of Converted Photons vs E_{T};E_{T} (GeV)",etBin,etMin,etMax);
  book3DHistoVector(iBooker, p_convFractionVsPhi_,"1D","convFractionVsPhi","Fraction of Converted Photons vs #phi;#phi",phiBin,phiMin,phiMax);
  book2DHistoVector(iBooker, p_convFractionVsEta_, "1D","convFractionVsEta","Fraction of Converted Photons vs #eta;#eta",etaBin,etaMin,etaMax);

  //booking bad channel fraction histograms
  iGetter.setCurrentFolder(AllPath+"Et above 20 GeV/");
  book2DHistoVector(iBooker, p_badChannelsFractionVsPhi_, "1D","badChannelsFractionVsPhi","Fraction of Photons which have at least one bad channel vs #phi;#phi",phiBin,phiMin,phiMax);
  book2DHistoVector(iBooker, p_badChannelsFractionVsEta_, "1D","badChannelsFractionVsEta","Fraction of Photons which have at least one bad channel vs #eta;#eta",etaBin,etaMin, etaMax);
  book2DHistoVector(iBooker, p_badChannelsFractionVsEt_, "1D","badChannelsFractionVsEt", "Fraction of Photons which have at least one bad channel vs E_{T};E_{T} (GeV)",etBin,etMin,etMax);

  //making efficiency plots
  MonitorElement * dividend;
  MonitorElement * numerator;
  MonitorElement * denominator;

  currentFolder_.str("");
  currentFolder_ << AllPath << "Et above 20 GeV/";

  //HLT efficiency plots
  dividend    = retrieveHisto(iGetter,EffPath,"EfficiencyVsEtaHLT");
  numerator   = retrieveHisto(iGetter,EffPath,"phoEtaPostHLT");
  denominator = retrieveHisto(iGetter,EffPath,"phoEtaPreHLT");
  dividePlots(dividend,numerator,denominator);

  dividend    = retrieveHisto(iGetter,EffPath,"EfficiencyVsEtHLT");
  numerator   = retrieveHisto(iGetter,EffPath,"phoEtPostHLT");
  denominator = retrieveHisto(iGetter,EffPath,"phoEtPreHLT");
  dividePlots(dividend,numerator,denominator);

  //efficiencies vs Eta
  denominator = retrieveHisto(iGetter,currentFolder_.str(),"phoEta");

  dividend    = retrieveHisto(iGetter,EffPath,"EfficiencyVsEtaLoose");
  numerator   = retrieveHisto(iGetter,EffPath,"phoEtaLoose");
  dividePlots(dividend,numerator,denominator);

  dividend    = retrieveHisto(iGetter,EffPath,"EfficiencyVsEtaTight");
  numerator   = retrieveHisto(iGetter,EffPath,"phoEtaTight");
  dividePlots(dividend,numerator,denominator);

  //efficiencies vs Et
  denominator = retrieveHisto(iGetter,currentFolder_.str(),"phoEtAllEcal");

  dividend    = retrieveHisto(iGetter,EffPath,"EfficiencyVsEtLoose");
  numerator   = retrieveHisto(iGetter,EffPath,"phoEtLoose");
  dividePlots(dividend,numerator,denominator);

  dividend    = retrieveHisto(iGetter,EffPath,"EfficiencyVsEtTight");
  numerator   = retrieveHisto(iGetter,EffPath,"phoEtTight");
  dividePlots(dividend,numerator,denominator);

  //conversion fractions vs Eta
  dividend    = retrieveHisto(iGetter,EffPath,"ConvFractionVsEtaLoose");
  numerator   = retrieveHisto(iGetter,EffPath,"convEtaLoose");
  denominator = retrieveHisto(iGetter,EffPath,"phoEtaLoose");
  dividePlots(dividend,numerator,denominator);

  dividend    = retrieveHisto(iGetter,EffPath,"ConvFractionVsEtaTight");
  numerator   = retrieveHisto(iGetter,EffPath,"convEtaTight");
  denominator = retrieveHisto(iGetter,EffPath,"phoEtaTight");
  dividePlots(dividend,numerator,denominator);

  //conversion fractions vs Et
  dividend    = retrieveHisto(iGetter,EffPath,"ConvFractionVsEtLoose");
  numerator   = retrieveHisto(iGetter,EffPath,"convEtLoose");
  denominator = retrieveHisto(iGetter,EffPath,"phoEtLoose");
  dividePlots(dividend,numerator,denominator);

  dividend    = retrieveHisto(iGetter,EffPath,"ConvFractionVsEtTight");
  numerator   = retrieveHisto(iGetter,EffPath,"convEtTight");
  denominator = retrieveHisto(iGetter,EffPath,"phoEtTight");
  dividePlots(dividend,numerator,denominator);

  //conversion vertex recontruction efficiency
  dividend    = retrieveHisto(iGetter,EffPath,"VertexReconstructionEfficiencyVsEta");
  numerator   = retrieveHisto(iGetter,currentFolder_.str()+"Conversions/","phoConvEta");
  denominator = retrieveHisto(iGetter,EffPath,"phoEtaVertex");
  dividePlots(dividend,numerator,denominator);

  iGetter.setCurrentFolder(EffPath);

  iGetter.removeElement("phoEtaPreHLT");
  iGetter.removeElement("phoEtPreHLT");
  iGetter.removeElement("phoEtaPostHLT");
  iGetter.removeElement("phoEtPostHLT");
  iGetter.removeElement("phoEtaLoose");
  iGetter.removeElement("phoEtaTight");
  iGetter.removeElement("phoEtLoose");
  iGetter.removeElement("phoEtTight");
  iGetter.removeElement("phoEtaVertex");
  iGetter.removeElement("convEtaLoose");
  iGetter.removeElement("convEtaTight");
  iGetter.removeElement("convEtLoose");
  iGetter.removeElement("convEtTight");

  for(uint type=0;type!=types_.size();++type){

    for (int cut = 0; cut !=numberOfSteps_; ++cut) {

      currentFolder_.str("");
      currentFolder_ << "Egamma/"+analyzerName_+"/" << types_[type] << "Photons/Et above " << (cut+1)*cutStep_ << " GeV/";

      //making bad channel histograms

      //vs Et
      dividend    = retrieveHisto(iGetter,currentFolder_.str(),"badChannelsFractionVsEt");
      numerator   = retrieveHisto(iGetter,currentFolder_.str(),"phoEtBadChannels");
      denominator = retrieveHisto(iGetter,currentFolder_.str(),"phoEtAllEcal");
      dividePlots(dividend,numerator,denominator);

      //vs eta
      dividend    = retrieveHisto(iGetter,currentFolder_.str(),"badChannelsFractionVsEta");
      numerator   = retrieveHisto(iGetter,currentFolder_.str(),"phoEtaBadChannels");
      denominator = retrieveHisto(iGetter,currentFolder_.str(),"phoEta");
      dividePlots(dividend,numerator,denominator);

      //vs phi
      dividend    = retrieveHisto(iGetter,currentFolder_.str(),"badChannelsFractionVsPhi");
      numerator   = retrieveHisto(iGetter,currentFolder_.str(),"phoPhiBadChannels");
      denominator = retrieveHisto(iGetter,currentFolder_.str(),"phoPhiAllEcal");
      dividePlots(dividend,numerator,denominator);

      //making conversion fraction histograms

      //vs Et
      dividend    = retrieveHisto(iGetter,currentFolder_.str()+"Conversions/","convFractionVsEt");
      numerator   = retrieveHisto(iGetter,currentFolder_.str()+"Conversions/","phoConvEtAllEcal");
      denominator = retrieveHisto(iGetter,currentFolder_.str(),"phoEtAllEcal");
      dividePlots(dividend,numerator,denominator);

      //vs eta
      dividend    = retrieveHisto(iGetter,currentFolder_.str()+"Conversions/","convFractionVsEta");
      numerator   = retrieveHisto(iGetter,currentFolder_.str()+"Conversions/","phoConvEtaForEfficiency");
      denominator = retrieveHisto(iGetter,currentFolder_.str(),"phoEta");
      dividePlots(dividend,numerator,denominator);

      //vs phi
      dividend    = retrieveHisto(iGetter,currentFolder_.str()+"Conversions/","convFractionVsPhiAllEcal");
      numerator   = retrieveHisto(iGetter,currentFolder_.str()+"Conversions/","phoConvPhiForEfficiencyAllEcal");
      denominator = retrieveHisto(iGetter,currentFolder_.str(),"phoPhiAllEcal");
      dividePlots(dividend,numerator,denominator);
      dividend    = retrieveHisto(iGetter,currentFolder_.str()+"Conversions/","convFractionVsPhiBarrel");
      numerator   = retrieveHisto(iGetter,currentFolder_.str()+"Conversions/","phoConvPhiForEfficiencyBarrel");
      denominator = retrieveHisto(iGetter,currentFolder_.str(),"phoPhiBarrel");
      dividePlots(dividend,numerator,denominator);
      dividend    = retrieveHisto(iGetter,currentFolder_.str()+"Conversions/","convFractionVsPhiEndcaps");
      numerator   = retrieveHisto(iGetter,currentFolder_.str()+"Conversions/","phoConvPhiForEfficiencyEndcaps");
      denominator = retrieveHisto(iGetter,currentFolder_.str(),"phoPhiEndcaps");
      dividePlots(dividend,numerator,denominator);


      iGetter.setCurrentFolder(currentFolder_.str()+"Conversions/");
      iGetter.removeElement("phoConvEtaForEfficiency");
      iGetter.removeElement("phoConvPhiForEfficiencyAllEcal");
      iGetter.removeElement("phoConvPhiForEfficiencyBarrel");
      iGetter.removeElement("phoConvPhiForEfficiencyEndcaps");
    }
  }  
  
  // if(standAlone_) dbe_->save(outputFileName_);
  //else if(batch_) dbe_->save(inputFileName_);

}

Member Data Documentation

std::string PhotonOfflineClient::analyzerName_

private

Definition at line 102 of file PhotonOfflineClient.h.

bool PhotonOfflineClient::batch_

private

Definition at line 149 of file PhotonOfflineClient.h.

std::stringstream PhotonOfflineClient::currentFolder_

private

Definition at line 154 of file PhotonOfflineClient.h.

double PhotonOfflineClient::cutStep_

private

Definition at line 133 of file PhotonOfflineClient.h.

DQMStore* PhotonOfflineClient::dbe_

private

Definition at line 128 of file PhotonOfflineClient.h.

int PhotonOfflineClient::etaBin

private

Definition at line 143 of file PhotonOfflineClient.h.

double PhotonOfflineClient::etaMax

private

Definition at line 142 of file PhotonOfflineClient.h.

double PhotonOfflineClient::etaMin

private

Definition at line 141 of file PhotonOfflineClient.h.

int PhotonOfflineClient::etBin

private

Definition at line 140 of file PhotonOfflineClient.h.

double PhotonOfflineClient::etMax

private

Definition at line 139 of file PhotonOfflineClient.h.

double PhotonOfflineClient::etMin

private

Definition at line 138 of file PhotonOfflineClient.h.

bool PhotonOfflineClient::excludeBkgHistos_

private

Definition at line 136 of file PhotonOfflineClient.h.

int PhotonOfflineClient::histo_index_conversions_

private

Definition at line 157 of file PhotonOfflineClient.h.

int PhotonOfflineClient::histo_index_efficiency_

private

Definition at line 158 of file PhotonOfflineClient.h.

int PhotonOfflineClient::histo_index_invMass_

private

Definition at line 159 of file PhotonOfflineClient.h.

int PhotonOfflineClient::histo_index_photons_

private

Definition at line 156 of file PhotonOfflineClient.h.

std::string PhotonOfflineClient::inputFileName_

private

Definition at line 152 of file PhotonOfflineClient.h.

bool PhotonOfflineClient::minimalSetOfHistos_

private

Definition at line 135 of file PhotonOfflineClient.h.

int PhotonOfflineClient::numberOfSteps_

private

Definition at line 134 of file PhotonOfflineClient.h.

std::string PhotonOfflineClient::outputFileName_

private

Definition at line 151 of file PhotonOfflineClient.h.

std::vector<std::vector<MonitorElement*> > PhotonOfflineClient::p_badChannelsFractionVsEt_

private

Definition at line 121 of file PhotonOfflineClient.h.

std::vector<std::vector<MonitorElement*> > PhotonOfflineClient::p_badChannelsFractionVsEta_

private

Definition at line 119 of file PhotonOfflineClient.h.

std::vector<std::vector<MonitorElement*> > PhotonOfflineClient::p_badChannelsFractionVsPhi_

private

Definition at line 120 of file PhotonOfflineClient.h.

std::vector<std::vector<MonitorElement*> > PhotonOfflineClient::p_convFractionVsEt_

private

Definition at line 117 of file PhotonOfflineClient.h.

std::vector<std::vector<MonitorElement*> > PhotonOfflineClient::p_convFractionVsEta_

private

Definition at line 115 of file PhotonOfflineClient.h.

MonitorElement* PhotonOfflineClient::p_convFractionVsEtaLoose_

private

Definition at line 110 of file PhotonOfflineClient.h.

MonitorElement* PhotonOfflineClient::p_convFractionVsEtaTight_

private

Definition at line 112 of file PhotonOfflineClient.h.

MonitorElement* PhotonOfflineClient::p_convFractionVsEtLoose_

private

Definition at line 111 of file PhotonOfflineClient.h.

MonitorElement* PhotonOfflineClient::p_convFractionVsEtTight_

private

Definition at line 113 of file PhotonOfflineClient.h.

std::vector<std::vector<std::vector<MonitorElement*> > > PhotonOfflineClient::p_convFractionVsPhi_

private

Definition at line 116 of file PhotonOfflineClient.h.

MonitorElement* PhotonOfflineClient::p_efficiencyVsEtaHLT_

private

Definition at line 107 of file PhotonOfflineClient.h.

MonitorElement* PhotonOfflineClient::p_efficiencyVsEtaLoose_

private

Definition at line 103 of file PhotonOfflineClient.h.

MonitorElement* PhotonOfflineClient::p_efficiencyVsEtaTight_

private

Definition at line 105 of file PhotonOfflineClient.h.

MonitorElement* PhotonOfflineClient::p_efficiencyVsEtHLT_

private

Definition at line 108 of file PhotonOfflineClient.h.

MonitorElement* PhotonOfflineClient::p_efficiencyVsEtLoose_

private

Definition at line 104 of file PhotonOfflineClient.h.

MonitorElement* PhotonOfflineClient::p_efficiencyVsEtTight_

private

Definition at line 106 of file PhotonOfflineClient.h.

MonitorElement* PhotonOfflineClient::p_vertexReconstructionEfficiencyVsEta_

private

Definition at line 123 of file PhotonOfflineClient.h.

edm::ParameterSet PhotonOfflineClient::parameters_

private

Definition at line 131 of file PhotonOfflineClient.h.

Referenced by Mixins._TypedParameterizable::clone(), Types.PSet::clone(), and Mixins._TypedParameterizable::copy().

std::vector<std::string> PhotonOfflineClient::parts_

private

Definition at line 162 of file PhotonOfflineClient.h.

int PhotonOfflineClient::phiBin

private

Definition at line 146 of file PhotonOfflineClient.h.

double PhotonOfflineClient::phiMax

private

Definition at line 145 of file PhotonOfflineClient.h.

double PhotonOfflineClient::phiMin

private

Definition at line 144 of file PhotonOfflineClient.h.

bool PhotonOfflineClient::standAlone_

private

Definition at line 148 of file PhotonOfflineClient.h.

std::vector<std::string> PhotonOfflineClient::types_

private

Definition at line 161 of file PhotonOfflineClient.h.

int PhotonOfflineClient::verbosity_

private

Definition at line 129 of file PhotonOfflineClient.h.