#include <TrackProbabilityTagPlotter.h>

Inheritance diagram for TrackProbabilityTagPlotter:

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

void	analyzeTag (const reco::BaseTagInfo *tagInfo, 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)

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

	~TrackProbabilityTagPlotter ()

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
EffPurFromHistos *	effPurFromHistos [4]

double	endEffPur_

unsigned int	mcPlots_

int	nBinEffPur_

double	startEffPur_

FlavourHistograms< double > *	tkcntHistosSig2D [5]

FlavourHistograms< double > *	tkcntHistosSig3D [5]

bool	willFinalize_

Additional Inherited Members
Protected Attributes inherited from BaseBTagPlotter
const EtaPtBin	etaPtBin_

const std::string	tagName_

const std::string	theExtensionString

Detailed Description

Definition at line 12 of file TrackProbabilityTagPlotter.h.

Constructor & Destructor Documentation

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

Definition at line 7 of file TrackProbabilityTagPlotter.cc.

References dir, AlCaHLTBitMon_QueryRunRegistry::string, BaseBTagPlotter::theExtensionString, tkcntHistosSig2D, tkcntHistosSig3D, and willFinalize_.

                                                                                                 :
   BaseTagInfoPlotter(tagName, etaPtBin),
   nBinEffPur_(pSet.getParameter<int>("nBinEffPur")),
   startEffPur_(pSet.getParameter<double>("startEffPur")),
   endEffPur_(pSet.getParameter<double>("endEffPur")),
   mcPlots_(mc), willFinalize_(wf)
 {
   const std::string dir(theExtensionString.substr(1));
   if (willFinalize_) return;
 
   tkcntHistosSig3D[4] = new FlavourHistograms<double>
        ("ips_3D" + theExtensionString, "3D Probability of impact parameter",
     50, -1.0, 1.0, false, true, true, "b", dir, mc, ibook) ;
 
   tkcntHistosSig3D[0] = new FlavourHistograms<double>
        ("ips1_3D" + theExtensionString, "3D Probability of impact parameter 1st trk",
     50, -1.0, 1.0, false, true, true, "b", dir, mc, ibook) ;
 
   tkcntHistosSig3D[1] = new FlavourHistograms<double>
        ("ips2_3D" + theExtensionString, "3D Probability of impact parameter 2nd trk",
     50, -1.0, 1.0, false, true, true, "b", dir, mc, ibook) ;
 
   tkcntHistosSig3D[2] = new FlavourHistograms<double>
        ("ips3_3D" + theExtensionString, "3D Probability of impact parameter 3rd trk",
     50, -1.0, 1.0, false, true, true, "b", dir, mc, ibook) ;
 
   tkcntHistosSig3D[3] = new FlavourHistograms<double>
        ("ips4_3D" + theExtensionString, "3D Probability of impact parameter 4th trk",
     50, -1.0, 1.0, false, true, true, "b", dir, mc, ibook) ;
 
   tkcntHistosSig2D[4] = new FlavourHistograms<double>
        ("ips_2D" + theExtensionString, "2D Probability of impact parameter",
     50, -1.0, 1.0, false, true, true, "b", dir, mc, ibook) ;
 
   tkcntHistosSig2D[0] = new FlavourHistograms<double>
        ("ips1_2D" + theExtensionString, "2D Probability of impact parameter 1st trk",
     50, -1.0, 1.0, false, true, true, "b", dir, mc, ibook) ;
 
   tkcntHistosSig2D[1] = new FlavourHistograms<double>
        ("ips2_2D" + theExtensionString, "2D Probability of impact parameter 2nd trk",
     50, -1.0, 1.0, false, true, true, "b", dir, mc, ibook) ;
 
   tkcntHistosSig2D[2] = new FlavourHistograms<double>
        ("ips3_2D" + theExtensionString, "2D Probability of impact parameter 3rd trk",
     50, -1.0, 1.0, false, true, true, "b", dir, mc, ibook) ;
 
   tkcntHistosSig2D[3] = new FlavourHistograms<double>
        ("ips4" + theExtensionString, "2D Probability of impact parameter 4th trk",
     50, -1.0, 1.0, false, true, true, "b", dir, mc, ibook) ;
 }

TrackProbabilityTagPlotter::~TrackProbabilityTagPlotter ( )

Definition at line 61 of file TrackProbabilityTagPlotter.cc.

References effPurFromHistos, gen::n, tkcntHistosSig2D, tkcntHistosSig3D, and willFinalize_.

 {
 
   if (willFinalize_) {
     for(int n=1; n != 3; ++n) {
       delete tkcntHistosSig3D[n];
       delete effPurFromHistos[n-1];
     }
     for(int n=1; n != 3; ++n) {
       delete tkcntHistosSig2D[n];
       delete effPurFromHistos[n+1];
     }
     return;
   }
   for(int n=0; n != 5; ++n) {
     delete tkcntHistosSig2D[n];
     delete tkcntHistosSig3D[n];
   }
 }

Member Function Documentation

void TrackProbabilityTagPlotter::analyzeTag	(	const reco::BaseTagInfo *	tagInfo,
		const double &	jec,
		const int &	jetFlavour
	)

virtual

Reimplemented from BaseTagInfoPlotter.

Definition at line 106 of file TrackProbabilityTagPlotter.cc.

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

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

virtual

Reimplemented from BaseTagInfoPlotter.

Definition at line 82 of file TrackProbabilityTagPlotter.cc.

References Exception, lumiContext::fill, gen::n, reco::TrackProbabilityTagInfo::probability(), reco::TrackProbabilityTagInfo::selectedTracks(), dumpRecoGeometry_cfg::tagInfo, tkcntHistosSig2D, tkcntHistosSig3D, and w.

 {
   const reco::TrackProbabilityTagInfo * tagInfo = 
     dynamic_cast<const reco::TrackProbabilityTagInfo *>(baseTagInfo);
 
   if (!tagInfo) {
     throw cms::Exception("Configuration")
       << "BTagPerformanceAnalyzer: Extended TagInfo not of type TrackProbabilityTagInfo. " << endl;
   }
 
   for(int n=0; n != tagInfo->selectedTracks(1) && n != 4; ++n)
     tkcntHistosSig2D[n]->fill(jetFlavour, tagInfo->probability(n,1),w);
   for(int n=0; n != tagInfo->selectedTracks(0) && n != 4; ++n)
     tkcntHistosSig3D[n]->fill(jetFlavour, tagInfo->probability(n,0),w);
 
   for(int n=0; n != tagInfo->selectedTracks(1); ++n)
     tkcntHistosSig2D[4]->fill(jetFlavour, tagInfo->probability(n,1),w);
   for(int n=0; n != tagInfo->selectedTracks(0); ++n)
     tkcntHistosSig3D[4]->fill(jetFlavour, tagInfo->probability(n,0),w);
 }

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

virtual

Implements BaseBTagPlotter.

Definition at line 201 of file TrackProbabilityTagPlotter.cc.

References effPurFromHistos, FlavourHistograms< T >::epsPlot(), gen::n, tkcntHistosSig2D, tkcntHistosSig3D, and willFinalize_.

 {
   if (willFinalize_) {
     for(int n=0; n != 4; ++n) effPurFromHistos[n]->epsPlot(name);
     return;
   }
   for(int n=0; n != 5; ++n) {
     tkcntHistosSig2D[n]->epsPlot(name);
     tkcntHistosSig3D[n]->epsPlot(name);
   }
 }

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

virtual

Implements BaseBTagPlotter.

Definition at line 112 of file TrackProbabilityTagPlotter.cc.

References bookConverter::compute(), dir, effPurFromHistos, endEffPur_, mcPlots_, gen::n, nBinEffPur_, startEffPur_, AlCaHLTBitMon_QueryRunRegistry::string, BaseBTagPlotter::theExtensionString, tkcntHistosSig2D, and tkcntHistosSig3D.

 {
   //
   // final processing:
   // produce the misid. vs. eff histograms
   //
   const std::string dir("TrackProbability"+theExtensionString);
 
   tkcntHistosSig3D[1] = new FlavourHistograms<double>
        ("ips2_3D" + theExtensionString, "3D Probability of impact parameter 2nd trk",
     50, -1.0, 1.0, "b", dir, mcPlots_, igetter_) ;
   effPurFromHistos[0] = new EffPurFromHistos (tkcntHistosSig3D[1],dir,mcPlots_, ibook, 
                           nBinEffPur_, startEffPur_, endEffPur_);
 
   tkcntHistosSig3D[2] = new FlavourHistograms<double>
        ("ips3_3D" + theExtensionString, "3D Probability of impact parameter 3rd trk",
     50, -1.0, 1.0, "b", dir, mcPlots_, igetter_) ;
   effPurFromHistos[1] = new EffPurFromHistos (tkcntHistosSig3D[2],dir,mcPlots_, ibook, 
                           nBinEffPur_, startEffPur_, endEffPur_);
 
   tkcntHistosSig2D[1] = new FlavourHistograms<double>
        ("ips2_2D" + theExtensionString, "2D Probability of impact parameter 2nd trk",
     50, -1.0, 1.0, "b", dir, mcPlots_, igetter_) ;
   effPurFromHistos[2] = new EffPurFromHistos (tkcntHistosSig2D[1],dir,mcPlots_, ibook, 
                           nBinEffPur_, startEffPur_, endEffPur_);
 
   tkcntHistosSig2D[2] = new FlavourHistograms<double>
        ("ips3_2D" + theExtensionString, "2D Probability of impact parameter 3rd trk",
     50, -1.0, 1.0, "b", dir, mcPlots_, igetter_) ;
   effPurFromHistos[3] = new EffPurFromHistos (tkcntHistosSig2D[2],dir,mcPlots_, ibook, 
                           nBinEffPur_, startEffPur_, endEffPur_);
 
   for(int n=0; n != 4; ++n) effPurFromHistos[n]->compute(ibook);
 }

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

virtual

Implements BaseBTagPlotter.

Definition at line 147 of file TrackProbabilityTagPlotter.cc.

References svgfig::canvas(), EffPurFromHistos::discriminatorCutEfficScan(), EffPurFromHistos::discriminatorNoCutEffic(), effPurFromHistos, gen::n, EffPurFromHistos::plot(), FlavourHistograms< T >::plot(), setTDRStyle(), AlCaHLTBitMon_QueryRunRegistry::string, BaseBTagPlotter::theExtensionString, tkcntHistosSig2D, tkcntHistosSig3D, and willFinalize_.

 {
   const std::string cName("TrackProbabilityPlots"+ theExtensionString);
   setTDRStyle()->cd();
   TCanvas canvas(cName.c_str(), cName.c_str(), 600, 900);
   canvas.UseCurrentStyle();
   if (willFinalize_) {
     for(int n=0; n != 2; ++n) {
       canvas.Print((name + cName + ".ps").c_str());
       canvas.Clear();
       canvas.Divide(2,3);
       canvas.cd(1);
       effPurFromHistos[0+n]->discriminatorNoCutEffic()->plot();
       canvas.cd(2);
       effPurFromHistos[0+n]->discriminatorCutEfficScan()->plot();
       canvas.cd(3);
       effPurFromHistos[0+n]->plot();
       canvas.cd(4);
       effPurFromHistos[1+n]->discriminatorNoCutEffic()->plot();
       canvas.cd(5);
       effPurFromHistos[1+n]->discriminatorCutEfficScan()->plot();
       canvas.cd(6);
       effPurFromHistos[1+n]->plot();
     }
     return;
   }
 
   canvas.Clear();
   canvas.Divide(2,3);
   canvas.Print((name + cName + ".ps[").c_str());
   canvas.cd(1);
 
   tkcntHistosSig3D[4]->plot();
   for(int n=0; n != 4; ++n) {
     canvas.cd(2+n);
     tkcntHistosSig3D[n]->plot();
   }
 
   canvas.Print((name + cName + ".ps").c_str());
   canvas.Clear();
   canvas.Divide(2,3);
 
   canvas.cd(1);
   tkcntHistosSig2D[4]->plot();
   for(int n=0; n != 4; ++n) {
     canvas.cd(2+n);
     tkcntHistosSig2D[n]->plot();
   }
 
   canvas.Print((name + cName + ".ps").c_str());
   canvas.Print((name + cName + ".ps]").c_str());
 }