#include <PhotonOfflineClient.h>

Inheritance diagram for PhotonOfflineClient:

Public Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)

virtual void	beginJob ()

std::vector< std::vector < MonitorElement * > >	book2DHistoVector (std::string histoType, std::string histoName, std::string title, int xbin, double xmin, double xmax, int ybin=1, double ymin=1, double ymax=2)

std::vector< std::vector < std::vector< MonitorElement * > > >	book3DHistoVector (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 (std::string histoName, std::string title, int bin, double min, double max)

virtual void	endJob ()

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

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

	PhotonOfflineClient (const edm::ParameterSet &pset)

MonitorElement *	retrieveHisto (std::string dir, std::string name)

virtual void	runClient ()

virtual	~PhotonOfflineClient ()

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

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

Definition at line 65 of file PhotonOfflineClient.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 25 of file PhotonOfflineClient.cc.

References dbe_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), cppFunctionSkipper::operator, and DQMStore::setVerbose().

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

virtual

Definition at line 69 of file PhotonOfflineClient.cc.

70 {}

Member Function Documentation

void PhotonOfflineClient::analyze	(	const edm::Event &	e,
		const edm::EventSetup &	esup
	)

virtual

Implements edm::EDAnalyzer.

Definition at line 73 of file PhotonOfflineClient.cc.

74 {}

void PhotonOfflineClient::beginJob ( void )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 71 of file PhotonOfflineClient.cc.

72 {}

vector< vector< MonitorElement * > > PhotonOfflineClient::book2DHistoVector	(	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 364 of file PhotonOfflineClient.cc.

References DQMStore::book1D(), DQMStore::book2D(), DQMStore::bookProfile(), GOODCOLL_filter_cfg::cut, dbe_, DQMStore::pwd(), DQMStore::setCurrentFolder(), indexGen::title, SiStripMonitorClusterAlca_cfi::xmax, SiStripMonitorClusterAlca_cfi::xmin, SiStripMonitorClusterAlca_cfi::ymax, and SiStripMonitorClusterAlca_cfi::ymin.

 {
   int histo_index = 0;
 
   vector<MonitorElement*> temp1DVector;
   vector<vector<MonitorElement*> > temp2DVector;
 
   //determining which folder we're in
   bool conversionPlot = false;
   if(dbe_->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";
 
       dbe_->setCurrentFolder(currentFolder_.str());
 
       string kind;
       if(conversionPlot) kind = " Conversions: ";
       else kind = " Photons: ";
 
       if(histoType=="1D")           temp1DVector.push_back(dbe_->book1D(histo_number_stream.str()+histoName,types_[type]+kind+title,xbin,xmin,xmax));
       else if(histoType=="2D")      temp1DVector.push_back(dbe_->book2D(histo_number_stream.str()+histoName,types_[type]+kind+title,xbin,xmin,xmax,ybin,ymin,ymax));
       else if(histoType=="Profile") temp1DVector.push_back(dbe_->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();
   }
 
   return temp2DVector;
 
 }

vector< vector< vector< MonitorElement * > > > PhotonOfflineClient::book3DHistoVector	(	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 424 of file PhotonOfflineClient.cc.

References DQMStore::book1D(), DQMStore::book2D(), DQMStore::bookProfile(), GOODCOLL_filter_cfg::cut, dbe_, DQMStore::pwd(), DQMStore::setCurrentFolder(), indexGen::title, SiStripMonitorClusterAlca_cfi::xmax, SiStripMonitorClusterAlca_cfi::xmin, SiStripMonitorClusterAlca_cfi::ymax, and SiStripMonitorClusterAlca_cfi::ymin.

 {
   int histo_index = 0;
 
   vector<MonitorElement*> temp1DVector;
   vector<vector<MonitorElement*> > temp2DVector;
   vector<vector<vector<MonitorElement*> > > temp3DVector;
 
 
   //determining which folder we're in
   bool conversionPlot = false;
   if(dbe_->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";
 
     dbe_->setCurrentFolder(currentFolder_.str());
 
     string kind;
     if(conversionPlot) kind = " Conversions: ";
     else kind = " Photons: ";
 
     if(histoType=="1D")           temp1DVector.push_back(dbe_->book1D(      histo_number_stream.str()+histoName+parts_[part],types_[type]+kind+parts_[part]+": "+title,xbin,xmin,xmax));
     else if(histoType=="2D")      temp1DVector.push_back(dbe_->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(dbe_->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();
   }
 
   return temp3DVector;
 }

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

Definition at line 340 of file PhotonOfflineClient.cc.

References newFWLiteAna::bin, DQMStore::book1D(), dbe_, bookConverter::max, min(), DQMStore::pwd(), and indexGen::title.

 {
 
   int histo_index = 0;
   stringstream histo_number_stream;
 
 
   //determining which folder we're in
   if(dbe_->pwd().find( "InvMass" ) != string::npos){
     histo_index_invMass_++;
     histo_index = histo_index_invMass_;
   }
   if(dbe_->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 dbe_->book1D(histo_number_stream.str()+"_"+histoName,title,bin,min,max);
 
 }

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

private

Definition at line 300 of file PhotonOfflineClient.cc.

References MonitorElement::getBinContent(), MonitorElement::getEntries(), MonitorElement::getNbinsX(), j, 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 322 of file PhotonOfflineClient.cc.

References MonitorElement::getBinContent(), MonitorElement::getNbinsX(), j, 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::endJob ( void )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 78 of file PhotonOfflineClient.cc.

 {
   if(standAlone_) runClient();
 }

void PhotonOfflineClient::endLuminosityBlock	(	const edm::LuminosityBlock &	,
		const edm::EventSetup &	setup
	)

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 75 of file PhotonOfflineClient.cc.

76 {}

void PhotonOfflineClient::endRun	(	const edm::Run &	run,
		const edm::EventSetup &	setup
	)

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 82 of file PhotonOfflineClient.cc.

 {
   if(!standAlone_) runClient();
 }

MonitorElement * PhotonOfflineClient::retrieveHisto	(	std::string	dir,
		std::string	name
	)

Definition at line 491 of file PhotonOfflineClient.cc.

References dbe_, DQMStore::getContents(), and cmsHarvester::index.

                                                                          {
   //cout << "dir = " << dir << endl;
   //cout << "name = " << name << endl;
   vector<MonitorElement*> histoVector;
   uint indexOfRelevantHistogram=0;
   string fullMEName = "";
   histoVector = dbe_->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 ( )

virtual

Definition at line 87 of file PhotonOfflineClient.cc.

References gather_cfg::cout, GOODCOLL_filter_cfg::cut, dbe_, cuy::denominator, DQMStore::dirExists(), DQMStore::get(), MonitorElement::getIntValue(), cuy::numerator, DQMStore::open(), DQMStore::removeElement(), DQMStore::save(), and DQMStore::setCurrentFolder().

 {
   if(!dbe_) return;
 
   if(batch_)  dbe_->open(inputFileName_);
 
   //std::cout << " PostProcessing analyzer name " << analyzerName_ << std::endl;
   if(!dbe_->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_     = dbe_->get("Egamma/"+analyzerName_+"/numberOfHistogramsInPhotonsFolder")->getIntValue();
   histo_index_conversions_ = dbe_->get("Egamma/"+analyzerName_+"/numberOfHistogramsInConversionsFolder")->getIntValue();
   histo_index_efficiency_  = dbe_->get("Egamma/"+analyzerName_+"/numberOfHistogramsInEfficiencyFolder")->getIntValue();
   histo_index_invMass_     = dbe_->get("Egamma/"+analyzerName_+"/numberOfHistogramsInInvMassFolder")->getIntValue();
 
   dbe_->setCurrentFolder("Egamma/"+analyzerName_+"/");
   dbe_->removeElement("numberOfHistogramsInPhotonsFolder");
   dbe_->removeElement("numberOfHistogramsInConversionsFolder");
   dbe_->removeElement("numberOfHistogramsInEfficiencyFolder");
   dbe_->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
   dbe_->setCurrentFolder(EffPath);
 
   p_efficiencyVsEtaLoose_ = bookHisto("EfficiencyVsEtaLoose","Fraction of Photons passing Loose Isolation vs #eta;#eta",etaBin,etaMin, etaMax);
   p_efficiencyVsEtLoose_  = bookHisto("EfficiencyVsEtLoose", "Fraction of Photons passing Loose Isolation vs E_{T};E_{T} (GeV)",etBin,etMin,etMax);
   p_efficiencyVsEtaTight_ = bookHisto("EfficiencyVsEtaTight","Fraction of Photons passing Tight Isolation vs #eta;#eta",etaBin,etaMin, etaMax);
   p_efficiencyVsEtTight_  = bookHisto("EfficiencyVsEtTight", "Fraction of Photons passing Tight Isolation vs E_{T};E_{T} (GeV)",etBin,etMin,etMax);
 
   p_efficiencyVsEtaHLT_ = bookHisto("EfficiencyVsEtaHLT","Fraction of Photons firing HLT vs #eta;#eta",etaBin,etaMin, etaMax);
   p_efficiencyVsEtHLT_  = bookHisto("EfficiencyVsEtHLT", "Fraction of Photons firing HLT vs E_{T};E_{T} (GeV)",etBin,etMin,etMax);
 
   p_convFractionVsEtaLoose_ = bookHisto("ConvFractionVsEtaLoose","Fraction of Loosely Isolated Photons which are matched to two tracks vs #eta;#eta",etaBin,etaMin, etaMax);
   p_convFractionVsEtLoose_  = bookHisto("ConvFractionVsEtLoose", "Fraction of Loosely Isolated Photons which are matched to two tracks vs E_{T};E_{T} (GeV)",etBin,etMin,etMax);
   p_convFractionVsEtaTight_ = bookHisto("ConvFractionVsEtaTight","Fraction of Tightly Isolated Photons which are matched to two tracks vs #eta;#eta",etaBin,etaMin, etaMax);
   p_convFractionVsEtTight_  = bookHisto("ConvFractionVsEtTight", "Fraction of Tightly Isolated Photons which are matched to two tracks vs E_{T};E_{T} (GeV)",etBin,etMin,etMax);
 
   p_vertexReconstructionEfficiencyVsEta_ = bookHisto("VertexReconstructionEfficiencyVsEta","Fraction of Converted Photons which have a valid vertex vs #eta;#eta",etaBin,etaMin, etaMax);
 
   //booking conversion fraction histograms
   dbe_->setCurrentFolder(AllPath+"Et above 20 GeV/Conversions");
   p_convFractionVsEt_  = book2DHistoVector("1D","convFractionVsEt", "Fraction of Converted Photons vs E_{T};E_{T} (GeV)",etBin,etMin,etMax);
   p_convFractionVsPhi_ = book3DHistoVector("1D","convFractionVsPhi","Fraction of Converted Photons vs #phi;#phi",phiBin,phiMin,phiMax);
   p_convFractionVsEta_ = book2DHistoVector("1D","convFractionVsEta","Fraction of Converted Photons vs #eta;#eta",etaBin,etaMin,etaMax);
 
   //booking bad channel fraction histograms
   dbe_->setCurrentFolder(AllPath+"Et above 20 GeV/");
   p_badChannelsFractionVsPhi_ = book2DHistoVector("1D","badChannelsFractionVsPhi","Fraction of Photons which have at least one bad channel vs #phi;#phi",phiBin,phiMin,phiMax);
   p_badChannelsFractionVsEta_ = book2DHistoVector("1D","badChannelsFractionVsEta","Fraction of Photons which have at least one bad channel vs #eta;#eta",etaBin,etaMin, etaMax);
   p_badChannelsFractionVsEt_  = book2DHistoVector("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(EffPath,"EfficiencyVsEtaHLT");
   numerator   = retrieveHisto(EffPath,"phoEtaPostHLT");
   denominator = retrieveHisto(EffPath,"phoEtaPreHLT");
   dividePlots(dividend,numerator,denominator);
 
   dividend    = retrieveHisto(EffPath,"EfficiencyVsEtHLT");
   numerator   = retrieveHisto(EffPath,"phoEtPostHLT");
   denominator = retrieveHisto(EffPath,"phoEtPreHLT");
   dividePlots(dividend,numerator,denominator);
 
   //efficiencies vs Eta
   denominator = retrieveHisto(currentFolder_.str(),"phoEta");
 
   dividend    = retrieveHisto(EffPath,"EfficiencyVsEtaLoose");
   numerator   = retrieveHisto(EffPath,"phoEtaLoose");
   dividePlots(dividend,numerator,denominator);
 
   dividend    = retrieveHisto(EffPath,"EfficiencyVsEtaTight");
   numerator   = retrieveHisto(EffPath,"phoEtaTight");
   dividePlots(dividend,numerator,denominator);
 
   //efficiencies vs Et
   denominator = retrieveHisto(currentFolder_.str(),"phoEtAllEcal");
 
   dividend    = retrieveHisto(EffPath,"EfficiencyVsEtLoose");
   numerator   = retrieveHisto(EffPath,"phoEtLoose");
   dividePlots(dividend,numerator,denominator);
 
   dividend    = retrieveHisto(EffPath,"EfficiencyVsEtTight");
   numerator   = retrieveHisto(EffPath,"phoEtTight");
   dividePlots(dividend,numerator,denominator);
 
   //conversion fractions vs Eta
   dividend    = retrieveHisto(EffPath,"ConvFractionVsEtaLoose");
   numerator   = retrieveHisto(EffPath,"convEtaLoose");
   denominator = retrieveHisto(EffPath,"phoEtaLoose");
   dividePlots(dividend,numerator,denominator);
 
   dividend    = retrieveHisto(EffPath,"ConvFractionVsEtaTight");
   numerator   = retrieveHisto(EffPath,"convEtaTight");
   denominator = retrieveHisto(EffPath,"phoEtaTight");
   dividePlots(dividend,numerator,denominator);
 
   //conversion fractions vs Et
   dividend    = retrieveHisto(EffPath,"ConvFractionVsEtLoose");
   numerator   = retrieveHisto(EffPath,"convEtLoose");
   denominator = retrieveHisto(EffPath,"phoEtLoose");
   dividePlots(dividend,numerator,denominator);
 
   dividend    = retrieveHisto(EffPath,"ConvFractionVsEtTight");
   numerator   = retrieveHisto(EffPath,"convEtTight");
   denominator = retrieveHisto(EffPath,"phoEtTight");
   dividePlots(dividend,numerator,denominator);
 
   //conversion vertex recontruction efficiency
   dividend    = retrieveHisto(EffPath,"VertexReconstructionEfficiencyVsEta");
   numerator   = retrieveHisto(currentFolder_.str()+"Conversions/","phoConvEta");
   denominator = retrieveHisto(EffPath,"phoEtaVertex");
   dividePlots(dividend,numerator,denominator);
 
   dbe_->setCurrentFolder(EffPath);
 
   dbe_->removeElement("phoEtaPreHLT");
   dbe_->removeElement("phoEtPreHLT");
   dbe_->removeElement("phoEtaPostHLT");
   dbe_->removeElement("phoEtPostHLT");
   dbe_->removeElement("phoEtaLoose");
   dbe_->removeElement("phoEtaTight");
   dbe_->removeElement("phoEtLoose");
   dbe_->removeElement("phoEtTight");
   dbe_->removeElement("phoEtaVertex");
   dbe_->removeElement("convEtaLoose");
   dbe_->removeElement("convEtaTight");
   dbe_->removeElement("convEtLoose");
   dbe_->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(currentFolder_.str(),"badChannelsFractionVsEt");
       numerator   = retrieveHisto(currentFolder_.str(),"phoEtBadChannels");
       denominator = retrieveHisto(currentFolder_.str(),"phoEtAllEcal");
       dividePlots(dividend,numerator,denominator);
 
       //vs eta
       dividend    = retrieveHisto(currentFolder_.str(),"badChannelsFractionVsEta");
       numerator   = retrieveHisto(currentFolder_.str(),"phoEtaBadChannels");
       denominator = retrieveHisto(currentFolder_.str(),"phoEta");
       dividePlots(dividend,numerator,denominator);
 
       //vs phi
       dividend    = retrieveHisto(currentFolder_.str(),"badChannelsFractionVsPhi");
       numerator   = retrieveHisto(currentFolder_.str(),"phoPhiBadChannels");
       denominator = retrieveHisto(currentFolder_.str(),"phoPhiAllEcal");
       dividePlots(dividend,numerator,denominator);
 
       //making conversion fraction histograms
 
       //vs Et
       dividend    = retrieveHisto(currentFolder_.str()+"Conversions/","convFractionVsEt");
       numerator   = retrieveHisto(currentFolder_.str()+"Conversions/","phoConvEtAllEcal");
       denominator = retrieveHisto(currentFolder_.str(),"phoEtAllEcal");
       dividePlots(dividend,numerator,denominator);
 
       //vs eta
       dividend    = retrieveHisto(currentFolder_.str()+"Conversions/","convFractionVsEta");
       numerator   = retrieveHisto(currentFolder_.str()+"Conversions/","phoConvEtaForEfficiency");
       denominator = retrieveHisto(currentFolder_.str(),"phoEta");
       dividePlots(dividend,numerator,denominator);
 
       //vs phi
       dividend    = retrieveHisto(currentFolder_.str()+"Conversions/","convFractionVsPhiAllEcal");
       numerator   = retrieveHisto(currentFolder_.str()+"Conversions/","phoConvPhiForEfficiencyAllEcal");
       denominator = retrieveHisto(currentFolder_.str(),"phoPhiAllEcal");
       dividePlots(dividend,numerator,denominator);
       dividend    = retrieveHisto(currentFolder_.str()+"Conversions/","convFractionVsPhiBarrel");
       numerator   = retrieveHisto(currentFolder_.str()+"Conversions/","phoConvPhiForEfficiencyBarrel");
       denominator = retrieveHisto(currentFolder_.str(),"phoPhiBarrel");
       dividePlots(dividend,numerator,denominator);
       dividend    = retrieveHisto(currentFolder_.str()+"Conversions/","convFractionVsPhiEndcaps");
       numerator   = retrieveHisto(currentFolder_.str()+"Conversions/","phoConvPhiForEfficiencyEndcaps");
       denominator = retrieveHisto(currentFolder_.str(),"phoPhiEndcaps");
       dividePlots(dividend,numerator,denominator);
 
 
       dbe_->setCurrentFolder(currentFolder_.str()+"Conversions/");
       dbe_->removeElement("phoConvEtaForEfficiency");
       dbe_->removeElement("phoConvPhiForEfficiencyAllEcal");
       dbe_->removeElement("phoConvPhiForEfficiencyBarrel");
       dbe_->removeElement("phoConvPhiForEfficiencyEndcaps");
     }
   }  
   
   if(standAlone_) dbe_->save(outputFileName_);
   else if(batch_) dbe_->save(inputFileName_);
 
 }

Member Data Documentation

std::string PhotonOfflineClient::analyzerName_

private

Definition at line 97 of file PhotonOfflineClient.h.

bool PhotonOfflineClient::batch_

private

Definition at line 144 of file PhotonOfflineClient.h.

std::stringstream PhotonOfflineClient::currentFolder_

private

Definition at line 149 of file PhotonOfflineClient.h.

double PhotonOfflineClient::cutStep_

private

Definition at line 128 of file PhotonOfflineClient.h.

DQMStore* PhotonOfflineClient::dbe_

private

Definition at line 123 of file PhotonOfflineClient.h.

int PhotonOfflineClient::etaBin

private

Definition at line 138 of file PhotonOfflineClient.h.

double PhotonOfflineClient::etaMax

private

Definition at line 137 of file PhotonOfflineClient.h.

double PhotonOfflineClient::etaMin

private

Definition at line 136 of file PhotonOfflineClient.h.

int PhotonOfflineClient::etBin

private

Definition at line 135 of file PhotonOfflineClient.h.

double PhotonOfflineClient::etMax

private

Definition at line 134 of file PhotonOfflineClient.h.

double PhotonOfflineClient::etMin

private

Definition at line 133 of file PhotonOfflineClient.h.

bool PhotonOfflineClient::excludeBkgHistos_

private

Definition at line 131 of file PhotonOfflineClient.h.

int PhotonOfflineClient::histo_index_conversions_

private

Definition at line 152 of file PhotonOfflineClient.h.

int PhotonOfflineClient::histo_index_efficiency_

private

Definition at line 153 of file PhotonOfflineClient.h.

int PhotonOfflineClient::histo_index_invMass_

private

Definition at line 154 of file PhotonOfflineClient.h.

int PhotonOfflineClient::histo_index_photons_

private

Definition at line 151 of file PhotonOfflineClient.h.

std::string PhotonOfflineClient::inputFileName_

private

Definition at line 147 of file PhotonOfflineClient.h.

bool PhotonOfflineClient::minimalSetOfHistos_

private

Definition at line 130 of file PhotonOfflineClient.h.

int PhotonOfflineClient::numberOfSteps_

private

Definition at line 129 of file PhotonOfflineClient.h.

std::string PhotonOfflineClient::outputFileName_

private

Definition at line 146 of file PhotonOfflineClient.h.

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

private

Definition at line 116 of file PhotonOfflineClient.h.

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

private

Definition at line 114 of file PhotonOfflineClient.h.

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

private

Definition at line 115 of file PhotonOfflineClient.h.

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

private

Definition at line 112 of file PhotonOfflineClient.h.

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

private

Definition at line 110 of file PhotonOfflineClient.h.

MonitorElement* PhotonOfflineClient::p_convFractionVsEtaLoose_

private

Definition at line 105 of file PhotonOfflineClient.h.

MonitorElement* PhotonOfflineClient::p_convFractionVsEtaTight_

private

Definition at line 107 of file PhotonOfflineClient.h.

MonitorElement* PhotonOfflineClient::p_convFractionVsEtLoose_

private

Definition at line 106 of file PhotonOfflineClient.h.

MonitorElement* PhotonOfflineClient::p_convFractionVsEtTight_

private

Definition at line 108 of file PhotonOfflineClient.h.

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

private

Definition at line 111 of file PhotonOfflineClient.h.

MonitorElement* PhotonOfflineClient::p_efficiencyVsEtaHLT_

private

Definition at line 102 of file PhotonOfflineClient.h.

MonitorElement* PhotonOfflineClient::p_efficiencyVsEtaLoose_

private

Definition at line 98 of file PhotonOfflineClient.h.

MonitorElement* PhotonOfflineClient::p_efficiencyVsEtaTight_

private

Definition at line 100 of file PhotonOfflineClient.h.

MonitorElement* PhotonOfflineClient::p_efficiencyVsEtHLT_

private

Definition at line 103 of file PhotonOfflineClient.h.

MonitorElement* PhotonOfflineClient::p_efficiencyVsEtLoose_

private

Definition at line 99 of file PhotonOfflineClient.h.

MonitorElement* PhotonOfflineClient::p_efficiencyVsEtTight_

private

Definition at line 101 of file PhotonOfflineClient.h.

MonitorElement* PhotonOfflineClient::p_vertexReconstructionEfficiencyVsEta_

private

Definition at line 118 of file PhotonOfflineClient.h.

edm::ParameterSet PhotonOfflineClient::parameters_

private

Definition at line 126 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 157 of file PhotonOfflineClient.h.

int PhotonOfflineClient::phiBin

private

Definition at line 141 of file PhotonOfflineClient.h.

double PhotonOfflineClient::phiMax

private

Definition at line 140 of file PhotonOfflineClient.h.

double PhotonOfflineClient::phiMin

private

Definition at line 139 of file PhotonOfflineClient.h.

bool PhotonOfflineClient::standAlone_

private

Definition at line 143 of file PhotonOfflineClient.h.

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

private

Definition at line 156 of file PhotonOfflineClient.h.

int PhotonOfflineClient::verbosity_

private

Definition at line 124 of file PhotonOfflineClient.h.

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation