#include <SoftLeptonTagPlotter.h>

Inheritance diagram for SoftLeptonTagPlotter:

Public Member Functions
void	analyzeTag (const reco::BaseTagInfo *baseTagInfo, const double &jec, const int &jetFlavour)

void	analyzeTag (const reco::BaseTagInfo *baseTagInfo, const double &jec, const int &jetFlavour, const float &w)

void	epsPlot (const std::string &name)

virtual void	finalize (DQMStore::IBooker &ibook_, DQMStore::IGetter &igetter_)

void	psPlot (const std::string &name)

	SoftLeptonTagPlotter (const std::string &tagName, const EtaPtBin &etaPtBin, const edm::ParameterSet &pSet, const unsigned int &mc, const bool &willFinalize, DQMStore::IBooker &ibook)

	~SoftLeptonTagPlotter (void)

Public Member Functions inherited from BaseTagInfoPlotter
virtual void	analyzeTag (const std::vector< const reco::BaseTagInfo * > &tagInfos, const double &jec, const int &jetFlavour)

virtual void	analyzeTag (const std::vector< const reco::BaseTagInfo * > &tagInfos, const double &jec, const int &jetFlavour, const float &w)

	BaseTagInfoPlotter (const std::string &tagName, const EtaPtBin &etaPtBin)

virtual void	setEventSetup (const edm::EventSetup &setup)

virtual std::vector< std::string >	tagInfoRequirements () const

virtual	~BaseTagInfoPlotter ()

Public Member Functions inherited from BaseBTagPlotter
	BaseBTagPlotter (const std::string &tagName, const EtaPtBin &etaPtBin)

const EtaPtBin &	etaPtBin ()

virtual	~BaseBTagPlotter ()

Private Attributes
FlavourHistograms< double > *	m_deltaR [s_leptons]

FlavourHistograms< double > *	m_etaRel [s_leptons]

FlavourHistograms< double > *	m_leptonId [s_leptons]

FlavourHistograms< double > *	m_leptonPt [s_leptons]

FlavourHistograms< double > *	m_p0Par [s_leptons]

FlavourHistograms< double > *	m_ptRel [s_leptons]

FlavourHistograms< double > *	m_ratio [s_leptons]

FlavourHistograms< double > *	m_ratioRel [s_leptons]

FlavourHistograms< double > *	m_sip2d [s_leptons]

FlavourHistograms< double > *	m_sip2dsig [s_leptons]

FlavourHistograms< double > *	m_sip3d [s_leptons]

FlavourHistograms< double > *	m_sip3dsig [s_leptons]

unsigned int	mcPlots_

bool	willFinalize_

Static Private Attributes
static const int	s_leptons = 2

Additional Inherited Members
Protected Attributes inherited from BaseBTagPlotter
const EtaPtBin	etaPtBin_

const std::string	tagName_

const std::string	theExtensionString

Detailed Description

Definition at line 10 of file SoftLeptonTagPlotter.h.

Constructor & Destructor Documentation

SoftLeptonTagPlotter::SoftLeptonTagPlotter	(	const std::string &	tagName,
		const EtaPtBin &	etaPtBin,
		const edm::ParameterSet &	pSet,
		const unsigned int &	mc,
		const bool &	willFinalize,
		DQMStore::IBooker &	ibook
	)

Definition at line 13 of file SoftLeptonTagPlotter.cc.

References i, m_deltaR, m_etaRel, m_leptonId, m_leptonPt, m_p0Par, m_ptRel, m_ratio, m_ratioRel, m_sip2d, m_sip2dsig, m_sip3d, m_sip3dsig, mcPlots_, ordinal, alignCSCRings::s, s_leptons, AlCaHLTBitMon_QueryRunRegistry::string, BaseBTagPlotter::theExtensionString, and willFinalize_.

                                                                                         :
   BaseTagInfoPlotter(tagName, etaPtBin), mcPlots_(mc), willFinalize_(wf)
 {
   const std::string softLepDir(theExtensionString.substr(1));
   
   if(willFinalize_) return;
   
   for (int i = 0; i < s_leptons; i++) {
     std::ostringstream s("");
     s << ordinal[i] << " lepton ";
     m_leptonId[i] = new FlavourHistograms<double> (
                            s.str() + "id",
                            "Lepton identification discriminaint",
                            60, -0.1, 1.1, false, false, true, "b", softLepDir,mcPlots_, ibook);
     m_leptonPt[i] = new FlavourHistograms<double> (
                            s.str() + "pT",
                            "Lepton transverse moementum",
                            100, 0.0, 20.0, false, false, true, "b", softLepDir,mcPlots_, ibook);
     m_sip2dsig[i] = new FlavourHistograms<double> (
         s.str() + "sip2dsig",
         "Lepton signed 2D impact parameter significance",
         100, -20.0, 30.0, false, false, true, "b", softLepDir,mcPlots_, ibook);
     m_sip3dsig[i] = new FlavourHistograms<double> (
         s.str() + "sip3dsig",
         "Lepton signed 3D impact parameter significance",
         100, -20.0, 30.0, false, false, true, "b", softLepDir,mcPlots_, ibook);
     m_sip2d[i] = new FlavourHistograms<double> (
         s.str() + "sip2d",
         "Lepton signed 2D impact parameter",
         100, -20.0, 30.0, false, false, true, "b", softLepDir,mcPlots_, ibook);
     m_sip3d[i] = new FlavourHistograms<double> (
         s.str() + "sip3d",
         "Lepton signed 3D impact parameter",
         100, -20.0, 30.0, false, false, true, "b", softLepDir,mcPlots_, ibook);
     m_ptRel[i] = new FlavourHistograms<double> (
         s.str() +  "pT rel",
         "Lepton transverse moementum relative to jet axis",
         100, 0.0, 10.0, false, false, true, "b", softLepDir,mcPlots_, ibook);
     m_p0Par[i] = new FlavourHistograms<double> (
         s.str() + "p0 par",
         "Lepton moementum along jet axis in the B rest frame",
         100, 0.0, 10.0, false, false, true, "b", softLepDir,mcPlots_, ibook);
     m_etaRel[i] = new FlavourHistograms<double> (
         s.str() + "eta rel",
         "Lepton pseudorapidity relative to jet axis",
         100, -5.0, 25.0, false, false, true, "b", softLepDir,mcPlots_, ibook);
     m_deltaR[i] = new FlavourHistograms<double> (
         s.str() + "delta R",
         "Lepton pseudoangular distance from jet axis",
         100, 0.0, 0.6, false, false, true, "b", softLepDir,mcPlots_, ibook);
     m_ratio[i] = new FlavourHistograms<double> (
         s.str() + "energy ratio",
         "Ratio of lepton momentum to jet energy",
         100, 0.0, 2.0, false, false, true, "b", softLepDir,mcPlots_, ibook);
     m_ratioRel[i] = new FlavourHistograms<double> (
         s.str() + "parallel energy ratio",
         "Ratio of lepton momentum along the jet axis to jet energy",
         100, 0.0, 2.0, false, false, true, "b", softLepDir,mcPlots_, ibook);
   }
 }

SoftLeptonTagPlotter::~SoftLeptonTagPlotter ( void )

Definition at line 76 of file SoftLeptonTagPlotter.cc.

References i, m_deltaR, m_etaRel, m_leptonId, m_leptonPt, m_p0Par, m_ptRel, m_ratio, m_ratioRel, m_sip2d, m_sip2dsig, m_sip3d, m_sip3dsig, s_leptons, and willFinalize_.

 {
   if(willFinalize_) return;
 
   for (int i = 0; i != s_leptons; ++i) {
     delete m_leptonId[i];
     delete m_leptonPt[i];
     delete m_sip2dsig[i];
     delete m_sip3dsig[i];
     delete m_sip2d[i];
     delete m_sip3d[i];
     delete m_ptRel[i];
     delete m_p0Par[i];
     delete m_etaRel[i];
     delete m_deltaR[i];
     delete m_ratio[i];
     delete m_ratioRel[i];
   }
 }

Member Function Documentation

void SoftLeptonTagPlotter::analyzeTag	(	const reco::BaseTagInfo *	baseTagInfo,
		const double &	jec,
		const int &	jetFlavour
	)

virtual

Reimplemented from BaseTagInfoPlotter.

Definition at line 96 of file SoftLeptonTagPlotter.cc.

 {
   analyzeTag(baseTagInfo,jetFlavour,1.);
 }

void SoftLeptonTagPlotter::analyzeTag	(	const reco::BaseTagInfo *	baseTagInfo,
		const double &	jec,
		const int &	jetFlavour,
		const float &	w
	)

virtual

Reimplemented from BaseTagInfoPlotter.

Definition at line 102 of file SoftLeptonTagPlotter.cc.

 {
 
   const reco::SoftLeptonTagInfo * tagInfo = 
     dynamic_cast<const reco::SoftLeptonTagInfo *>(baseTagInfo);
 
   if (!tagInfo) {
     throw cms::Exception("Configuration")
       << "BTagPerformanceAnalyzer: Extended TagInfo not of type SoftLeptonTagInfo. " << endl;
   }
 
   int n_leptons = tagInfo->leptons();
 
   for (int i = 0; i != n_leptons && i != s_leptons; ++i) {
     const reco::SoftLeptonProperties& properties = tagInfo->properties(i);
     m_leptonPt[i]->fill( jetFlavour, tagInfo->lepton(i)->pt() ,w);
     m_leptonId[i]->fill( jetFlavour, properties.quality() ,w);
     m_sip2dsig[i]->fill(    jetFlavour, properties.sip2dsig ,w);
     m_sip3dsig[i]->fill(    jetFlavour, properties.sip3dsig ,w);
     m_sip2d[i]->fill(    jetFlavour, properties.sip2d ,w);
     m_sip3d[i]->fill(    jetFlavour, properties.sip3d ,w);
     m_ptRel[i]->fill(    jetFlavour, properties.ptRel ,w);
     m_p0Par[i]->fill(    jetFlavour, properties.p0Par ,w);
     m_etaRel[i]->fill(   jetFlavour, properties.etaRel ,w);
     m_deltaR[i]->fill(   jetFlavour, properties.deltaR ,w);
     m_ratio[i]->fill(    jetFlavour, properties.ratio ,w);
     m_ratioRel[i]->fill( jetFlavour, properties.ratioRel ,w);
   }
 }

void SoftLeptonTagPlotter::epsPlot ( const std::string & name )

virtual

Implements BaseBTagPlotter.

Definition at line 178 of file SoftLeptonTagPlotter.cc.

References FlavourHistograms< T >::epsPlot(), i, m_deltaR, m_etaRel, m_leptonId, m_leptonPt, m_p0Par, m_ptRel, m_ratio, m_ratioRel, m_sip2d, m_sip2dsig, m_sip3d, m_sip3dsig, s_leptons, and willFinalize_.

 {
   if(willFinalize_) return;
 
   for (int i=0; i != s_leptons; ++i) {
     m_leptonId[i]->epsPlot( name );
     m_leptonPt[i]->epsPlot( name );
     m_sip2d[i]->epsPlot( name );
     m_sip3d[i]->epsPlot( name );
     m_sip2dsig[i]->epsPlot( name );
     m_sip3dsig[i]->epsPlot( name );
     m_ptRel[i]->epsPlot( name );
     m_p0Par[i]->epsPlot( name );
     m_etaRel[i]->epsPlot( name );
     m_deltaR[i]->epsPlot( name );
     m_ratio[i]->epsPlot( name );
     m_ratioRel[i]->epsPlot( name );
   }
 }

virtual void SoftLeptonTagPlotter::finalize	(	DQMStore::IBooker &	ibook_,
		DQMStore::IGetter &	igetter_
	)

inlinevirtual

Implements BaseBTagPlotter.

Definition at line 22 of file SoftLeptonTagPlotter.h.

22 {}

void SoftLeptonTagPlotter::psPlot ( const std::string & name )

virtual

Implements BaseBTagPlotter.

Definition at line 135 of file SoftLeptonTagPlotter.cc.

References svgfig::canvas(), i, m_deltaR, m_etaRel, m_leptonId, m_leptonPt, m_p0Par, m_ptRel, m_ratio, m_ratioRel, m_sip2d, m_sip2dsig, m_sip3d, m_sip3dsig, FlavourHistograms< T >::plot(), s_leptons, setTDRStyle(), AlCaHLTBitMon_QueryRunRegistry::string, BaseBTagPlotter::theExtensionString, and willFinalize_.

 {
   if(willFinalize_) return;
 
   const std::string cName("SoftLeptonPlots" + theExtensionString);
   setTDRStyle()->cd();
   TCanvas canvas(cName.c_str(), cName.c_str(), 600, 900);
   canvas.UseCurrentStyle();
   canvas.Divide(2,3);
   canvas.Print((name + cName + ".ps[").c_str());
   for (int i = 0; i < s_leptons; i++) {
     canvas.cd(1)->Clear();
     m_leptonId[i]->plot();
     canvas.cd(2)->Clear();
     m_leptonPt[i]->plot();
     canvas.cd(3)->Clear();
     m_sip2d[i]->plot();
     canvas.cd(4)->Clear();
     m_sip3d[i]->plot();
     canvas.cd(5)->Clear();
     m_sip2dsig[i]->plot();
     canvas.cd(6)->Clear();
     m_sip3dsig[i]->plot();
     canvas.Print((name + cName + ".ps").c_str());
 
     canvas.cd(1)->Clear();
     m_etaRel[i]->plot();
     canvas.cd(2)->Clear();
     m_deltaR[i]->plot();
     canvas.cd(3)->Clear();
     m_ratio[i]->plot();
     canvas.cd(4)->Clear();
     m_ratioRel[i]->plot();
     canvas.cd(5)->Clear();
     m_ptRel[i]->plot();
     canvas.cd(6)->Clear();
     m_p0Par[i]->plot();
     canvas.Print((name + cName + ".ps").c_str());
   }
   canvas.Print((name + cName + ".ps]").c_str());
 }