tmtt::Histos Class Reference

#include <Histos.h>

Public Member Functions
virtual bool	available () const
virtual void	book ()
virtual bool	enabled () const
virtual void	endJobAnalysis (const HTrphi::ErrorMonitor *htRphiErrMon=nullptr)
virtual void	fill (const InputData &inputData, const Array2D< std::unique_ptr< Sector >> &mSectors, const Array2D< std::unique_ptr< HTrphi >> &mHtPhis, const Array2D< std::unique_ptr< Make3Dtracks >> &mGet3Dtrks, const std::map< std::string, std::list< const L1fittedTrack * >> &mapFinalTracks)
	Histos (const Settings *settings)
virtual void	trackerGeometryAnalysis (const std::list< TrackerModule > &listTrackerModule)
virtual	~Histos ()=default

Protected Member Functions
virtual TFileDirectory	bookEtaPhiSectors ()
virtual TFileDirectory	bookInputData ()
virtual TFileDirectory	bookRphiHT ()
virtual TFileDirectory	bookRZfilters ()
virtual TFileDirectory	bookTrackCands (const std::string &tName)
virtual std::map< std::string, TFileDirectory >	bookTrackFitting ()
virtual void	fillEtaPhiSectors (const InputData &inputData, const Array2D< std::unique_ptr< Sector >> &mSectors)
virtual void	fillInputData (const InputData &inputData)
virtual void	fillRphiHT (const Array2D< std::unique_ptr< HTrphi >> &mHtRphis)
virtual void	fillRZfilters (const Array2D< std::unique_ptr< Make3Dtracks >> &mMake3Dtrks)
virtual void	fillTrackCands (const InputData &inputData, const Array2D< std::unique_ptr< Make3Dtracks >> &mMake3Dtrks, const std::string &tName)
virtual void	fillTrackCands (const InputData &inputData, const std::vector< L1track3D > &tracks, const std::string &tName)
virtual void	fillTrackFitting (const InputData &inputData, const std::map< std::string, std::list< const L1fittedTrack * >> &mapFinalTracks)
virtual void	makeEfficiencyPlot (TFileDirectory &inputDir, TEfficiency *outputEfficiency, TH1F *pass, TH1F *all, TString name, TString title)
virtual void	plotHybridDupRemovalEfficiency ()
virtual TFileDirectory	plotTrackEffAfterFit (const std::string &fitName)
virtual TFileDirectory	plotTrackEfficiency (const std::string &tName)
virtual void	plotTrackletSeedEfficiency ()
virtual void	printFitTrackPerformance (const std::string &fitName)
virtual void	printHybridDupRemovalPerformance ()
virtual void	printTrackletSeedFindingPerformance ()
virtual void	printTrackPerformance (const std::string &tName)

Protected Attributes
bool	bApproxMistake_
float	chosenRofZ_
edm::Service< TFileService >	fs_
unsigned int	genMinStubLayers_
TGraph *	graphBVsZoverR_
TH1F *	hisBendResStub_
TH1F *	hisBendStub_
std::map< std::string, TProfile * >	hisD0ResVsTrueEta_
std::map< std::string, TProfile * >	hisD0ResVsTrueInvPt_
std::map< std::string, TH1F * >	hisEta_
std::map< std::string, TH1F * >	hisEtaRes_
std::map< std::string, TProfile * >	hisEtaResVsTrueEta_
std::map< std::string, TProfile * >	hisEtaResVsTrueInvPt_
std::map< std::string, TH1F * >	hisFitChi2DofRphiMatched_
std::map< std::string, TH1F * >	hisFitChi2DofRphiUnmatched_
std::map< std::string, TH1F * >	hisFitChi2DofRzMatched_
std::map< std::string, TH1F * >	hisFitChi2DofRzUnmatched_
std::map< std::string, TH1F * >	hisFitD0Matched_
std::map< std::string, TH1F * >	hisFitD0Unmatched_
std::map< std::string, TH1F * >	hisFitEtaMatched_
std::map< std::string, TH1F * >	hisFitEtaUnmatched_
std::map< std::string, TH1F * >	hisFitPhi0Matched_
std::map< std::string, TH1F * >	hisFitPhi0Unmatched_
std::map< std::string, TH1F * >	hisFitQinvPtMatched_
std::map< std::string, TH1F * >	hisFitQinvPtUnmatched_
std::map< std::string, TH1F * >	hisFitTPd0ForAlgEff_
std::map< std::string, TH1F * >	hisFitTPd0ForEff_
std::map< std::string, TH1F * >	hisFitTPetaForAlgEff_
std::map< std::string, TH1F * >	hisFitTPetaForEff_
std::map< std::string, TH1F * >	hisFitTPinvptForAlgEff_
std::map< std::string, TH1F * >	hisFitTPinvptForEff_
std::map< std::string, TH1F * >	hisFitTPphiForAlgEff_
std::map< std::string, TH1F * >	hisFitTPphiForEff_
std::map< std::string, TH1F * >	hisFitTPz0ForAlgEff_
std::map< std::string, TH1F * >	hisFitTPz0ForEff_
std::map< std::string, TH1F * >	hisFitZ0Matched_
std::map< std::string, TH1F * >	hisFitZ0Unmatched_
std::map< std::string, TH1F * >	hisFracMatchStubsOnTracks_
std::map< std::string, TH1F * >	hisKalmanChi2DofSkipLay0Matched_
std::map< std::string, TH1F * >	hisKalmanChi2DofSkipLay0Unmatched_
std::map< std::string, TH1F * >	hisKalmanChi2DofSkipLay1Matched_
std::map< std::string, TH1F * >	hisKalmanChi2DofSkipLay1Unmatched_
std::map< std::string, TH1F * >	hisKalmanChi2DofSkipLay2Matched_
std::map< std::string, TH1F * >	hisKalmanChi2DofSkipLay2Unmatched_
std::map< std::string, TH1F * >	hisKalmanNumUpdateCalls_
std::map< std::string, TH1F * >	hisLayersPerTrack_
TH1F *	hisNumEtaSecsPerStub_
std::map< std::string, TH1F * >	hisNumFitTrks_
std::map< std::string, TH1F * >	hisNumFitTrksPerNon_
std::map< std::string, TH1F * >	hisNumFitTrksPerSect_
TH1F *	hisNumLayersPerTP_
TH1F *	hisNumPhiSecsPerStub_
TH1F *	hisNumPSLayersPerTP_
std::map< std::string, TH1F * >	hisNumStubsPerLink_
TH1F *	hisNumStubsPerSec_
std::map< std::string, TH1F * >	hisNumTracksVsQoverPt_
std::map< std::string, TH1F * >	hisNumTrksPerNon_
std::map< std::string, TH1F * >	hisPerfFitTPetaForAlgEff_
std::map< std::string, TH1F * >	hisPerfFitTPetaForEff_
std::map< std::string, TH1F * >	hisPerfFitTPinvptForAlgEff_
std::map< std::string, TH1F * >	hisPerfFitTPinvptForEff_
std::map< std::string, TH1F * >	hisPerfRecoTPetaForAlgEff_
std::map< std::string, TH1F * >	hisPerfRecoTPetaForEff_
std::map< std::string, TH1F * >	hisPerfRecoTPinvptForAlgEff_
std::map< std::string, TH1F * >	hisPerfRecoTPinvptForEff_
std::map< std::string, TH1F * >	hisPhi0_
std::map< std::string, TH1F * >	hisPhi0Res_
std::map< std::string, TProfile * >	hisPhi0ResVsTrueEta_
std::map< std::string, TProfile * >	hisPhi0ResVsTrueInvPt_
std::map< std::string, TH1F * >	hisQoverPt_
std::map< std::string, TH1F * >	hisQoverPtRes_
std::map< std::string, TProfile * >	hisQoverPtResVsTrueEta_
std::map< std::string, TProfile * >	hisQoverPtResVsTrueInvPt_
std::map< std::string, TH1F * >	hisRecoTPd0ForAlgEff_
std::map< std::string, TH1F * >	hisRecoTPd0ForEff_
std::map< std::string, TH1F * >	hisRecoTPetaForAlgEff_
std::map< std::string, TH1F * >	hisRecoTPetaForEff_
std::map< std::string, TH1F * >	hisRecoTPinvptForAlgEff_
std::map< std::string, TH1F * >	hisRecoTPinvptForEff_
std::map< std::string, TH1F * >	hisRecoTPphiForAlgEff_
std::map< std::string, TH1F * >	hisRecoTPphiForEff_
std::map< std::string, TH1F * >	hisRecoTPz0ForAlgEff_
std::map< std::string, TH1F * >	hisRecoTPz0ForEff_
TProfile *	hisStubIneffiVsEta_
TProfile *	hisStubIneffiVsInvPt_
TProfile *	hisStubKillFE_
std::map< std::string, TH1F * >	hisStubsOnTracksPerNon_
std::map< std::string, TH1F * >	hisStubsPerFitTrack_
std::map< std::string, TH1F * >	hisStubsPerTrack_
TH1F *	hisStubsVsEta_
TH1F *	hisStubsVsR_
TH1F *	hisTPd0ForAlgEff_
TH1F *	hisTPd0ForEff_
TH1F *	hisTPetaForAlgEff_
TH1F *	hisTPetaForEff_
TH1F *	hisTPinvptForAlgEff_
TH1F *	hisTPinvptForEff_
TH1F *	hisTPphiForAlgEff_
TH1F *	hisTPphiForEff_
TH1F *	hisTPz0ForAlgEff_
TH1F *	hisTPz0ForEff_
std::map< std::string, TH1F * >	hisZ0_
std::map< std::string, TH1F * >	hisZ0Res_
std::map< std::string, TProfile * >	hisZ0ResVsTrueEta_
std::map< std::string, TProfile * >	hisZ0ResVsTrueInvPt_
float	houghMinPt_
unsigned int	houghNbinsPhi_
unsigned int	houghNbinsPt_
std::map< unsigned int, float >	mapBarrelLayerMaxR_
std::map< unsigned int, float >	mapBarrelLayerMinR_
std::map< unsigned int, float >	mapEndcapWheelMaxZ_
std::map< unsigned int, float >	mapEndcapWheelMinZ_
std::map< unsigned int, float >	mapExtraAModuleTypeMaxR_
std::map< unsigned int, float >	mapExtraAModuleTypeMaxZ_
std::map< unsigned int, float >	mapExtraAModuleTypeMinR_
std::map< unsigned int, float >	mapExtraAModuleTypeMinZ_
std::map< unsigned int, float >	mapExtraBModuleTypeMaxR_
std::map< unsigned int, float >	mapExtraBModuleTypeMaxZ_
std::map< unsigned int, float >	mapExtraBModuleTypeMinR_
std::map< unsigned int, float >	mapExtraBModuleTypeMinZ_
std::map< unsigned int, float >	mapExtraCModuleTypeMaxR_
std::map< unsigned int, float >	mapExtraCModuleTypeMaxZ_
std::map< unsigned int, float >	mapExtraCModuleTypeMinR_
std::map< unsigned int, float >	mapExtraCModuleTypeMinZ_
std::map< unsigned int, float >	mapExtraDModuleTypeMaxR_
std::map< unsigned int, float >	mapExtraDModuleTypeMaxZ_
std::map< unsigned int, float >	mapExtraDModuleTypeMinR_
std::map< unsigned int, float >	mapExtraDModuleTypeMinZ_
std::map< unsigned int, float >	mapModuleTypeMaxR_
std::map< unsigned int, float >	mapModuleTypeMaxZ_
std::map< unsigned int, float >	mapModuleTypeMinR_
std::map< unsigned int, float >	mapModuleTypeMinZ_
unsigned int	numEtaRegions_
std::map< std::string, unsigned int >	numFitAlgEff_
std::map< std::string, unsigned int >	numFitAlgEffPass_
std::map< std::string, unsigned int >	numFitPerfAlgEff_
std::map< std::string, unsigned int >	numFitPerfAlgEffPass_
unsigned int	numPhiSectors_
bool	oldSumW2opt_
bool	printedGeomAnalysis_
std::map< std::string, TProfile * >	profDupFitTrksVsEta_
std::map< std::string, TProfile * >	profDupFitTrksVsInvPt_
std::map< std::string, TProfile * >	profDupTracksVsEta_
std::map< std::string, TProfile * >	profDupTracksVsInvPt_
std::map< std::string, TProfile * >	profFitChi2DofRphiVsInvPtMatched_
std::map< std::string, TProfile * >	profFitChi2DofRphiVsInvPtUnmatched_
std::map< std::string, TProfile * >	profMeanStubsPerLink_
std::map< std::string, TProfile * >	profNumFitTracks_
std::map< std::string, TProfile * >	profNumTrackCands_
std::map< std::string, TProfile * >	profNumTracksVsEta_
std::map< std::string, TProfile * >	profStubsOnFitTracks_
std::map< std::string, TProfile * >	profStubsOnTracks_
bool	ranRZfilter_
bool	resPlotOpt_
const Settings *	settings_
std::map< std::string, TEfficiency * >	teffAlgEffFitVsD0_
std::map< std::string, TEfficiency * >	teffAlgEffFitVsEta_
std::map< std::string, TEfficiency * >	teffAlgEffFitVsInvPt_
std::map< std::string, TEfficiency * >	teffAlgEffFitVsPhi_
std::map< std::string, TEfficiency * >	teffAlgEffFitVsZ0_
std::map< std::string, TEfficiency * >	teffAlgEffVsD0_
std::map< std::string, TEfficiency * >	teffAlgEffVsEta_
std::map< std::string, TEfficiency * >	teffAlgEffVsInvPt_
std::map< std::string, TEfficiency * >	teffAlgEffVsPhi_
std::map< std::string, TEfficiency * >	teffAlgEffVsZ0_
std::map< std::string, TEfficiency * >	teffEffFitVsD0_
std::map< std::string, TEfficiency * >	teffEffFitVsEta_
std::map< std::string, TEfficiency * >	teffEffFitVsInvPt_
std::map< std::string, TEfficiency * >	teffEffFitVsPhi_
std::map< std::string, TEfficiency * >	teffEffFitVsZ0_
std::map< std::string, TEfficiency * >	teffEffVsD0_
std::map< std::string, TEfficiency * >	teffEffVsEta_
std::map< std::string, TEfficiency * >	teffEffVsInvPt_
std::map< std::string, TEfficiency * >	teffEffVsPhi_
std::map< std::string, TEfficiency * >	teffEffVsZ0_
std::map< std::string, TEfficiency * >	teffPerfAlgEffFitVsEta_
std::map< std::string, TEfficiency * >	teffPerfAlgEffFitVsInvPt_
std::map< std::string, TEfficiency * >	teffPerfAlgEffVsEta_
std::map< std::string, TEfficiency * >	teffPerfAlgEffVsInvPt_
std::map< std::string, TEfficiency * >	teffPerfAlgEffVsPt_
std::map< std::string, TEfficiency * >	teffPerfEffFitVsEta_
std::map< std::string, TEfficiency * >	teffPerfEffFitVsInvPt_
std::map< std::string, TEfficiency * >	teffPerfEffVsEta_
std::map< std::string, TEfficiency * >	teffPerfEffVsInvPt_
std::vector< std::string >	trackFitters_
std::vector< std::string >	useRZfilter_

Detailed Description

Definition at line 35 of file Histos.h.

Constructor & Destructor Documentation

Histos::Histos

(

const Settings *

settings

)

Definition at line 39 of file Histos.cc.

References tmtt::Settings::chosenRofZ(), chosenRofZ_, tmtt::Settings::genMinStubLayers(), genMinStubLayers_, tmtt::Settings::houghMinPt(), houghMinPt_, tmtt::Settings::houghNbinsPhi(), houghNbinsPhi_, tmtt::Settings::houghNbinsPt(), houghNbinsPt_, tmtt::Settings::numEtaRegions(), numEtaRegions_, tmtt::Settings::numPhiSectors(), numPhiSectors_, ranRZfilter_, tmtt::Settings::resPlotOpt(), resPlotOpt_, tmtt::Settings::trackFitters(), trackFitters_, tmtt::Settings::useRZfilter(), and useRZfilter_.

                                         : settings_(settings), oldSumW2opt_(false), bApproxMistake_(false) {
    genMinStubLayers_ = settings->genMinStubLayers();
    numPhiSectors_ = settings->numPhiSectors();
    numEtaRegions_ = settings->numEtaRegions();
    houghMinPt_ = settings->houghMinPt();
    houghNbinsPt_ = settings->houghNbinsPt();
    houghNbinsPhi_ = settings->houghNbinsPhi();
    chosenRofZ_ = settings->chosenRofZ();
    trackFitters_ = settings->trackFitters();
    useRZfilter_ = settings->useRZfilter();
    ranRZfilter_ = (not useRZfilter_.empty());  // Was any r-z track filter run?
    resPlotOpt_ = settings->resPlotOpt();       // Only use signal events for helix resolution plots?
  }

virtual tmtt::Histos::~Histos ( )

virtualdefault

Member Function Documentation

virtual bool tmtt::Histos::available ( ) const

inlinevirtual

Definition at line 57 of file Histos.h.

References fs_, and edm::Service< T >::isAvailable().

Referenced by enabled().

{ return fs_.isAvailable(); }

void Histos::book ( )

virtual

Definition at line 55 of file Histos.cc.

References bookEtaPhiSectors(), bookInputData(), bookRphiHT(), bookRZfilters(), bookTrackCands(), bookTrackFitting(), enabled(), oldSumW2opt_, and ranRZfilter_.

Referenced by tmtt::GlobalCacheTMTT::GlobalCacheTMTT().

                    {
    // Don't bother booking histograms if user didn't request them via TFileService in their cfg.
    if (not this->enabled())
      return;

    oldSumW2opt_ = TH1::GetDefaultSumw2();
    TH1::SetDefaultSumw2(true);

    // Book histograms about input data.
    this->bookInputData();
    // Book histograms checking if (eta,phi) sector definition choices are good.
    this->bookEtaPhiSectors();
    // Book histograms checking filling of r-phi HT array.
    this->bookRphiHT();
    // Book histograms about r-z track filters.
    if (ranRZfilter_)
      this->bookRZfilters();
    // Book histograms studying 3D track candidates found after HT.
    this->bookTrackCands("HT");
    // Book histograms studying 3D track candidates found after r-z track filter.
    if (ranRZfilter_)
      this->bookTrackCands("RZ");
    // Book histograms studying track fitting performance
    this->bookTrackFitting();
  }

TFileDirectory Histos::bookEtaPhiSectors ( )

protectedvirtual

Definition at line 369 of file Histos.cc.

References fs_, hisNumEtaSecsPerStub_, hisNumPhiSecsPerStub_, hisNumStubsPerSec_, TFileDirectory::make(), and TFileService::mkdir().

Referenced by book().

                                           {
    TFileDirectory inputDir = fs_->mkdir("CheckSectors");

    // Check if stubs excessively duplicated between overlapping sectors.
    hisNumEtaSecsPerStub_ =
        inputDir.make<TH1F>("NumEtaSecPerStub", "; No. of #eta sectors each stub in", 20, -0.5, 19.5);
    hisNumPhiSecsPerStub_ =
        inputDir.make<TH1F>("NumPhiSecPerStub", "; No. of #phi sectors each stub in", 20, -0.5, 19.5);

    // Count stubs per (eta,phi) sector.
    hisNumStubsPerSec_ = inputDir.make<TH1F>("NumStubsPerSec", "; No. of stubs per sector", 150, -0.5, 299.5);

    return inputDir;
  }

TFileDirectory Histos::bookInputData ( )

protectedvirtual

Definition at line 121 of file Histos.cc.

References fs_, hisBendResStub_, hisBendStub_, hisNumLayersPerTP_, hisNumPSLayersPerTP_, hisStubIneffiVsEta_, hisStubIneffiVsInvPt_, hisStubKillFE_, hisStubsVsEta_, hisStubsVsR_, hisTPd0ForAlgEff_, hisTPd0ForEff_, hisTPetaForAlgEff_, hisTPetaForEff_, hisTPinvptForAlgEff_, hisTPinvptForEff_, hisTPphiForAlgEff_, hisTPphiForEff_, hisTPz0ForAlgEff_, hisTPz0ForEff_, M_PI, TFileDirectory::make(), and TFileService::mkdir().

Referenced by book().

                                       {
    TFileDirectory inputDir = fs_->mkdir("InputData");

    hisStubsVsEta_ = inputDir.make<TH1F>("StubsVsEta", "; #eta; No. stubs in tracker", 30, -3.0, 3.0);
    hisStubsVsR_ = inputDir.make<TH1F>("StubsVsR", "; radius (cm); No. stubs in tracker", 1200, 0., 120.);

    hisNumLayersPerTP_ =
        inputDir.make<TH1F>("NumLayersPerTP", "; Number of layers per TP for alg. eff.", 20, -0.5, 19.5);
    hisNumPSLayersPerTP_ =
        inputDir.make<TH1F>("NumPSLayersPerTP", "; Number of PS layers per TP for alg. eff.", 20, -0.5, 19.5);

    // Study efficiency of tightened front end-electronics cuts.

    hisStubKillFE_ = inputDir.make<TProfile>(
        "StubKillFE", "; barrelLayer or 10+endcapRing; Stub fraction rejected by FE chip", 30, -0.5, 29.5);
    hisStubIneffiVsInvPt_ =
        inputDir.make<TProfile>("StubIneffiVsPt", "; 1/Pt; Inefficiency of FE chip for good stubs", 25, 0.0, 0.5);
    hisStubIneffiVsEta_ =
        inputDir.make<TProfile>("StubIneffiVsEta", "; |#eta|; Inefficiency of FE chip for good stubs", 15, 0.0, 3.0);

    // Study stub resolution.

    hisBendStub_ = inputDir.make<TH1F>("BendStub", "; Stub bend in units of strips", 59, -7.375, 7.375);
    hisBendResStub_ = inputDir.make<TH1F>("BendResStub", "; Stub bend minus TP bend in units of strips", 100, -5., 5.);

    // Histos for denominator of tracking efficiency
    hisTPinvptForEff_ = inputDir.make<TH1F>("TPinvptForEff", "; TP 1/Pt (for effi.);", 50, 0., 0.5);
    hisTPetaForEff_ = inputDir.make<TH1F>("TPetaForEff", "; TP #eta (for effi.);", 20, -3., 3.);
    hisTPphiForEff_ = inputDir.make<TH1F>("TPphiForEff", "; TP #phi (for effi.);", 20, -M_PI, M_PI);
    hisTPd0ForEff_ = inputDir.make<TH1F>("TPd0ForEff", "; TP d0 (for effi.);", 40, 0., 4.);
    hisTPz0ForEff_ = inputDir.make<TH1F>("TPz0ForEff", "; TP z0 (for effi.);", 50, 0., 25.);
    //
    hisTPinvptForAlgEff_ = inputDir.make<TH1F>("TPinvptForAlgEff", "; TP 1/Pt (for alg. effi.);", 50, 0., 0.5);
    hisTPetaForAlgEff_ = inputDir.make<TH1F>("TPetaForAlgEff", "; TP #eta (for alg. effi.);", 20, -3., 3.);
    hisTPphiForAlgEff_ = inputDir.make<TH1F>("TPphiForAlgEff", "; TP #phi (for alg. effi.);", 20, -M_PI, M_PI);
    hisTPd0ForAlgEff_ = inputDir.make<TH1F>("TPd0ForAlgEff", "; TP d0 (for alg. effi.);", 40, 0., 4.);
    hisTPz0ForAlgEff_ = inputDir.make<TH1F>("TPz0ForAlgEff", "; TP z0 (for alg. effi.);", 50, 0., 25.);

    return inputDir;
  }

TFileDirectory Histos::bookRphiHT ( )

protectedvirtual

Definition at line 437 of file Histos.cc.

References fs_, and TFileService::mkdir().

Referenced by book().

                                    {
    TFileDirectory inputDir = fs_->mkdir("HTrphi");

    return inputDir;
  }

TFileDirectory Histos::bookRZfilters ( )

protectedvirtual

Definition at line 455 of file Histos.cc.

References fs_, and TFileService::mkdir().

Referenced by book().

                                       {
    TFileDirectory inputDir = fs_->mkdir("RZfilters");

    return inputDir;
  }

TFileDirectory Histos::bookTrackCands ( const std::string &  tName )

protectedvirtual

Definition at line 471 of file Histos.cc.

References fs_, hisEta_, hisEtaRes_, hisFracMatchStubsOnTracks_, hisLayersPerTrack_, hisNumStubsPerLink_, hisNumTracksVsQoverPt_, hisNumTrksPerNon_, hisPerfRecoTPetaForAlgEff_, hisPerfRecoTPetaForEff_, hisPerfRecoTPinvptForAlgEff_, hisPerfRecoTPinvptForEff_, hisPhi0_, hisPhi0Res_, hisQoverPt_, hisQoverPtRes_, hisRecoTPd0ForAlgEff_, hisRecoTPd0ForEff_, hisRecoTPetaForAlgEff_, hisRecoTPetaForEff_, hisRecoTPinvptForAlgEff_, hisRecoTPinvptForEff_, hisRecoTPphiForAlgEff_, hisRecoTPphiForEff_, hisRecoTPz0ForAlgEff_, hisRecoTPz0ForEff_, hisStubsOnTracksPerNon_, hisStubsPerTrack_, hisZ0_, hisZ0Res_, M_PI, TFileDirectory::make(), TFileService::mkdir(), HLT_FULL_cff::nEta, numEtaRegions_, profDupTracksVsEta_, profDupTracksVsInvPt_, profMeanStubsPerLink_, profNumTrackCands_, profNumTracksVsEta_, profStubsOnTracks_, and alignCSCRings::s.

Referenced by book().

                                                           {
    // Now book histograms for studying tracking in general.

    auto addn = [tName](const string& s) { return TString::Format("%s_%s", s.c_str(), tName.c_str()); };

    TFileDirectory inputDir = fs_->mkdir(addn("TrackCands").Data());

    bool TMTT = (tName == "HT" || tName == "RZ");

    // Count tracks in various ways (including/excluding duplicates, excluding fakes ...)
    profNumTrackCands_[tName] =
        inputDir.make<TProfile>(addn("NumTrackCands"), "; class; N. of tracks in tracker", 7, 0.5, 7.5);
    profNumTrackCands_[tName]->GetXaxis()->SetBinLabel(7, "TP for eff recoed");
    profNumTrackCands_[tName]->GetXaxis()->SetBinLabel(6, "TP recoed");
    profNumTrackCands_[tName]->GetXaxis()->SetBinLabel(5, "TP recoed x #eta sector dups");
    profNumTrackCands_[tName]->GetXaxis()->SetBinLabel(4, "TP recoed x sector dups");
    profNumTrackCands_[tName]->GetXaxis()->SetBinLabel(2, "TP recoed x track dups");
    profNumTrackCands_[tName]->GetXaxis()->SetBinLabel(1, "reco tracks including fakes");
    profNumTrackCands_[tName]->LabelsOption("d");

    hisNumTrksPerNon_[tName] = inputDir.make<TH1F>(addn("NumTrksPerNon"), "; No. tracks per nonant;", 100, -0.5, 399.5);

    unsigned int nEta = numEtaRegions_;
    hisNumTracksVsQoverPt_[tName] =
        inputDir.make<TH1F>(addn("NumTracksVsQoverPt"), "; Q/Pt; No. of tracks in tracker", 100, -0.5, 0.5);
    if (TMTT) {
      profNumTracksVsEta_[tName] = inputDir.make<TProfile>(
          addn("NumTracksVsEta"), "; #eta region; No. of tracks in tracker", nEta, -0.5, nEta - 0.5);
    }

    // Count stubs per event assigned to tracks (determines HT data output rate)

    profStubsOnTracks_[tName] =
        inputDir.make<TProfile>(addn("StubsOnTracks"), "; ; No. of stubs on tracks per event", 1, 0.5, 1.5);
    hisStubsOnTracksPerNon_[tName] =
        inputDir.make<TH1F>(addn("StubsOnTracksPerNon"), "; No. of stubs on tracks per nonant", 100, -0.5, 4999.5);

    hisStubsPerTrack_[tName] = inputDir.make<TH1F>(addn("StubsPerTrack"), ";No. of stubs per track;", 50, -0.5, 49.5);
    hisLayersPerTrack_[tName] =
        inputDir.make<TH1F>(addn("LayersPerTrack"), ";No. of layers with stubs per track;", 20, -0.5, 19.5);

    if (TMTT) {
      hisNumStubsPerLink_[tName] =
          inputDir.make<TH1F>(addn("NumStubsPerLink"), "; Mean #stubs per MHT output opto-link;", 50, -0.5, 249.5);
      profMeanStubsPerLink_[tName] =
          inputDir.make<TProfile>(addn("MeanStubsPerLink"), "; Mean #stubs per MHT output opto-link;", 36, -0.5, 35.5);
    }

    hisFracMatchStubsOnTracks_[tName] = inputDir.make<TH1F>(
        addn("FracMatchStubsOnTracks"), "; Frac. of stubs per trk matching best TP;", 101, -0.005, 1.005);

    if (TMTT) {
      // Study duplication of tracks within an individual HT array.
      profDupTracksVsEta_[tName] =
          inputDir.make<TProfile>(addn("DupTracksVsTPeta"), "; #eta; No. of dup. trks per TP;", 15, 0.0, 3.0);
      profDupTracksVsInvPt_[tName] =
          inputDir.make<TProfile>(addn("DupTracksVsInvPt"), "; 1/Pt; No. of dup. trks per TP", 25, 0., 0.5);
    }

    // Histos of track params.
    hisQoverPt_[tName] = inputDir.make<TH1F>(addn("QoverPt"), "; track q/Pt", 100, -0.5, 0.5);
    hisPhi0_[tName] = inputDir.make<TH1F>(addn("Phi0"), "; track #phi0", 70, -3.5, 3.5);
    hisEta_[tName] = inputDir.make<TH1F>(addn("Eta"), "; track #eta", 60, -3.0, 3.0);
    hisZ0_[tName] = inputDir.make<TH1F>(addn("Z0"), "; track z0", 100, -25.0, 25.0);

    // Histos of track parameter resolution
    hisQoverPtRes_[tName] = inputDir.make<TH1F>(addn("QoverPtRes"), "; track resolution in q/Pt", 100, -0.06, 0.06);
    hisPhi0Res_[tName] = inputDir.make<TH1F>(addn("Phi0Res"), "; track resolution in #phi0", 100, -0.04, 0.04);
    hisEtaRes_[tName] = inputDir.make<TH1F>(addn("EtaRes"), "; track resolution in #eta", 100, -1.0, 1.0);
    hisZ0Res_[tName] = inputDir.make<TH1F>(addn("Z0Res"), "; track resolution in z0", 100, -10.0, 10.0);

    // Histos for tracking efficiency vs. TP kinematics
    hisRecoTPinvptForEff_[tName] =
        inputDir.make<TH1F>(addn("RecoTPinvptForEff"), "; TP 1/Pt of recoed tracks (for effi.);", 50, 0., 0.5);
    hisRecoTPetaForEff_[tName] =
        inputDir.make<TH1F>(addn("RecoTPetaForEff"), "; TP #eta of recoed tracks (for effi.);", 20, -3., 3.);
    hisRecoTPphiForEff_[tName] =
        inputDir.make<TH1F>(addn("RecoTPphiForEff"), "; TP #phi of recoed tracks (for effi.);", 20, -M_PI, M_PI);

    // Histo for efficiency to reconstruct track perfectly (no incorrect hits).
    hisPerfRecoTPinvptForEff_[tName] = inputDir.make<TH1F>(
        addn("PerfRecoTPinvptForEff"), "; TP 1/Pt of recoed tracks (for perf. effi.);", 50, 0., 0.5);
    hisPerfRecoTPetaForEff_[tName] =
        inputDir.make<TH1F>(addn("PerfRecoTPetaForEff"), "; TP #eta of recoed tracks (for perf. effi.);", 20, -3., 3.);

    // Histos for  tracking efficiency vs. TP production point
    hisRecoTPd0ForEff_[tName] =
        inputDir.make<TH1F>(addn("RecoTPd0ForEff"), "; TP d0 of recoed tracks (for effi.);", 40, 0., 4.);
    hisRecoTPz0ForEff_[tName] =
        inputDir.make<TH1F>(addn("RecoTPz0ForEff"), "; TP z0 of recoed tracks (for effi.);", 50, 0., 25.);

    // Histos for algorithmic tracking efficiency vs. TP kinematics
    hisRecoTPinvptForAlgEff_[tName] =
        inputDir.make<TH1F>(addn("RecoTPinvptForAlgEff"), "; TP 1/Pt of recoed tracks (for alg. effi.);", 50, 0., 0.5);
    hisRecoTPetaForAlgEff_[tName] =
        inputDir.make<TH1F>(addn("RecoTPetaForAlgEff"), "; TP #eta of recoed tracks (for alg. effi.);", 20, -3., 3.);
    hisRecoTPphiForAlgEff_[tName] = inputDir.make<TH1F>(
        addn("RecoTPphiForAlgEff"), "; TP #phi of recoed tracks (for alg. effi.);", 20, -M_PI, M_PI);

    // Histo for efficiency to reconstruct track perfectly (no incorrect hits).
    hisPerfRecoTPinvptForAlgEff_[tName] = inputDir.make<TH1F>(
        addn("PerfRecoTPinvptForAlgEff"), "; TP 1/Pt of recoed tracks (for perf. alg. effi.);", 50, 0., 0.5);
    hisPerfRecoTPetaForAlgEff_[tName] =
        inputDir.make<TH1F>(addn("PerfRecoTPetaForAlgEff"), "; TP #eta (for perf. alg. effi.);", 20, -3., 3.);

    // Histos for algorithmic tracking efficiency vs. TP production point
    hisRecoTPd0ForAlgEff_[tName] =
        inputDir.make<TH1F>(addn("RecoTPd0ForAlgEff"), "; TP d0 of recoed tracks (for alg. effi.);", 40, 0., 4.);
    hisRecoTPz0ForAlgEff_[tName] =
        inputDir.make<TH1F>(addn("RecoTPz0ForAlgEff"), "; TP z0 of recoed tracks (for alg. effi.);", 50, 0., 25.);

    return inputDir;
  }

map< string, TFileDirectory > Histos::bookTrackFitting ( )

protectedvirtual

Definition at line 865 of file Histos.cc.

References fs_, hisD0ResVsTrueEta_, hisD0ResVsTrueInvPt_, hisEtaResVsTrueEta_, hisEtaResVsTrueInvPt_, hisFitChi2DofRphiMatched_, hisFitChi2DofRphiUnmatched_, hisFitChi2DofRzMatched_, hisFitChi2DofRzUnmatched_, hisFitD0Matched_, hisFitD0Unmatched_, hisFitEtaMatched_, hisFitEtaUnmatched_, hisFitPhi0Matched_, hisFitPhi0Unmatched_, hisFitQinvPtMatched_, hisFitQinvPtUnmatched_, hisFitTPd0ForAlgEff_, hisFitTPd0ForEff_, hisFitTPetaForAlgEff_, hisFitTPetaForEff_, hisFitTPinvptForAlgEff_, hisFitTPinvptForEff_, hisFitTPphiForAlgEff_, hisFitTPphiForEff_, hisFitTPz0ForAlgEff_, hisFitTPz0ForEff_, hisFitZ0Matched_, hisFitZ0Unmatched_, hisKalmanChi2DofSkipLay0Matched_, hisKalmanChi2DofSkipLay0Unmatched_, hisKalmanChi2DofSkipLay1Matched_, hisKalmanChi2DofSkipLay1Unmatched_, hisKalmanChi2DofSkipLay2Matched_, hisKalmanChi2DofSkipLay2Unmatched_, hisKalmanNumUpdateCalls_, hisNumFitTrks_, hisNumFitTrksPerNon_, hisPerfFitTPetaForAlgEff_, hisPerfFitTPetaForEff_, hisPerfFitTPinvptForAlgEff_, hisPerfFitTPinvptForEff_, hisPhi0ResVsTrueEta_, hisPhi0ResVsTrueInvPt_, hisQoverPtResVsTrueEta_, hisQoverPtResVsTrueInvPt_, hisStubsPerFitTrack_, hisZ0ResVsTrueEta_, hisZ0ResVsTrueInvPt_, M_PI, TFileDirectory::make(), TFileService::mkdir(), profDupFitTrksVsEta_, profDupFitTrksVsInvPt_, profFitChi2DofRphiVsInvPtMatched_, profFitChi2DofRphiVsInvPtUnmatched_, profNumFitTracks_, profStubsOnFitTracks_, alignCSCRings::s, and trackFitters_.

Referenced by book().

                                                       {
    map<string, TFileDirectory> inputDirMap;

    for (const string& fitName : trackFitters_) {
      // Define lambda function to facilitate adding "fitName" histogram names.
      auto addn = [fitName](const string& s) { return TString::Format("%s_%s", s.c_str(), fitName.c_str()); };

      TFileDirectory inputDir = fs_->mkdir(fitName);
      inputDirMap[fitName] = inputDir;

      profNumFitTracks_[fitName] =
          inputDir.make<TProfile>(addn("NumFitTracks"), "; class; # of fitted tracks", 11, 0.5, 11.5, -0.5, 9.9e6);
      profNumFitTracks_[fitName]->GetXaxis()->SetBinLabel(7, "TP for eff fitted");
      profNumFitTracks_[fitName]->GetXaxis()->SetBinLabel(6, "TP fitted");
      profNumFitTracks_[fitName]->GetXaxis()->SetBinLabel(2, "Fit tracks that are genuine");
      profNumFitTracks_[fitName]->GetXaxis()->SetBinLabel(1, "Fit tracks including fakes");
      profNumFitTracks_[fitName]->LabelsOption("d");

      hisNumFitTrks_[fitName] =
          inputDir.make<TH1F>(addn("NumFitTrks"), "; No. fitted tracks in tracker;", 200, -0.5, 399.5);
      hisNumFitTrksPerNon_[fitName] =
          inputDir.make<TH1F>(addn("NumFitTrksPerNon"), "; No. fitted tracks per nonant;", 200, -0.5, 199.5);

      hisStubsPerFitTrack_[fitName] =
          inputDir.make<TH1F>(addn("StubsPerFitTrack"), "; No. of stubs per fitted track", 20, -0.5, 19.5);
      profStubsOnFitTracks_[fitName] = inputDir.make<TProfile>(
          addn("StubsOnFitTracks"), "; ; No. of stubs on all fitted tracks per event", 1, 0.5, 1.5);

      hisFitQinvPtMatched_[fitName] =
          inputDir.make<TH1F>(addn("FitQinvPtMatched"), "Fitted q/p_{T} for matched tracks", 120, -0.6, 0.6);
      hisFitPhi0Matched_[fitName] =
          inputDir.make<TH1F>(addn("FitPhi0Matched"), "Fitted #phi_{0} for matched tracks", 70, -3.5, 3.5);
      hisFitD0Matched_[fitName] =
          inputDir.make<TH1F>(addn("FitD0Matched"), "Fitted d_{0} for matched tracks", 100, -2., 2.);
      hisFitZ0Matched_[fitName] =
          inputDir.make<TH1F>(addn("FitZ0Matched"), "Fitted z_{0} for matched tracks", 100, -25., 25.);
      hisFitEtaMatched_[fitName] =
          inputDir.make<TH1F>(addn("FitEtaMatched"), "Fitted #eta for matched tracks", 70, -3.5, 3.5);

      hisFitQinvPtUnmatched_[fitName] =
          inputDir.make<TH1F>(addn("FitQinvPtUnmatched"), "Fitted q/p_{T} for unmatched tracks", 120, -0.6, 0.6);
      hisFitPhi0Unmatched_[fitName] =
          inputDir.make<TH1F>(addn("FitPhi0Unmatched"), "Fitted #phi_{0} for unmatched tracks", 70, -3.5, 3.5);
      hisFitD0Unmatched_[fitName] =
          inputDir.make<TH1F>(addn("FitD0Unmatched"), "Fitted d_{0} for unmatched tracks", 100, -2., 2.);
      hisFitZ0Unmatched_[fitName] =
          inputDir.make<TH1F>(addn("FitZ0Unmatched"), "Fitted z_{0} for unmatched tracks", 100, -25., 25.);
      hisFitEtaUnmatched_[fitName] =
          inputDir.make<TH1F>(addn("FitEtaUnmatched"), "Fitted #eta for unmatched tracks", 70, -3.5, 3.5);

      const unsigned int nBinsChi2 = 39;
      const float chi2dofBins[nBinsChi2 + 1] = {0.0,  0.2,  0.4,  0.6,   0.8,   1.0,   1.2,   1.4,   1.6,   1.8,
                                                2.0,  2.4,  2.8,  3.2,   3.6,   4.0,   4.5,   5.0,   6.0,   7.0,
                                                8.0,  9.0,  10.0, 12.0,  14.0,  16.0,  18.0,  20.0,  25.0,  30.0,
                                                40.0, 50.0, 75.0, 100.0, 150.0, 200.0, 250.0, 350.0, 500.0, 1000.0};

      hisFitChi2DofRphiMatched_[fitName] =
          inputDir.make<TH1F>(addn("FitChi2DofRphiMatched"), ";#chi^{2}rphi;", nBinsChi2, chi2dofBins);
      hisFitChi2DofRzMatched_[fitName] =
          inputDir.make<TH1F>(addn("FitChi2DofRzMatched"), ";#chi^{2}rz/DOF;", nBinsChi2, chi2dofBins);
      profFitChi2DofRphiVsInvPtMatched_[fitName] =
          inputDir.make<TProfile>(addn("FitChi2DofRphiVsInvPtMatched"), "; 1/p_{T}; Fit #chi^{2}rphi/dof", 25, 0., 0.5);

      hisFitChi2DofRphiUnmatched_[fitName] =
          inputDir.make<TH1F>(addn("FitChi2DofRphiUnmatched"), ";#chi^{2}rphi/DOF;", nBinsChi2, chi2dofBins);
      hisFitChi2DofRzUnmatched_[fitName] =
          inputDir.make<TH1F>(addn("FitChi2DofRzUnmatched"), ";#chi^{2}rz/DOF;", nBinsChi2, chi2dofBins);
      profFitChi2DofRphiVsInvPtUnmatched_[fitName] = inputDir.make<TProfile>(
          addn("FitChi2DofRphiVsInvPtUnmatched"), "; 1/p_{T}; Fit #chi^{2}rphi/dof", 25, 0., 0.5);

      // Monitoring specific track fit algorithms.
      if (fitName.find("KF") != string::npos) {
        hisKalmanNumUpdateCalls_[fitName] =
            inputDir.make<TH1F>(addn("KalmanNumUpdateCalls"), "; Calls to KF updator;", 100, -0.5, 99.5);

        hisKalmanChi2DofSkipLay0Matched_[fitName] = inputDir.make<TH1F>(
            addn("KalmanChi2DofSkipLay0Matched"), ";#chi^{2} for nSkippedLayers = 0;", nBinsChi2, chi2dofBins);
        hisKalmanChi2DofSkipLay1Matched_[fitName] = inputDir.make<TH1F>(
            addn("KalmanChi2DofSkipLay1Matched"), ";#chi^{2} for nSkippedLayers = 1;", nBinsChi2, chi2dofBins);
        hisKalmanChi2DofSkipLay2Matched_[fitName] = inputDir.make<TH1F>(
            addn("KalmanChi2DofSkipLay2Matched"), ";#chi^{2} for nSkippedLayers = 2;", nBinsChi2, chi2dofBins);
        hisKalmanChi2DofSkipLay0Unmatched_[fitName] = inputDir.make<TH1F>(
            addn("KalmanChi2DofSkipLay0Unmatched"), ";#chi^{2} for nSkippedLayers = 0;", nBinsChi2, chi2dofBins);
        hisKalmanChi2DofSkipLay1Unmatched_[fitName] = inputDir.make<TH1F>(
            addn("KalmanChi2DofSkipLay1Unmatched"), ";#chi^{2} for nSkippedLayers = 1;", nBinsChi2, chi2dofBins);
        hisKalmanChi2DofSkipLay2Unmatched_[fitName] = inputDir.make<TH1F>(
            addn("KalmanChi2DofSkipLay2Unmatched"), ";#chi^{2} for nSkippedLayers = 2;", nBinsChi2, chi2dofBins);
      }

      // Plots of helix param resolution.

      hisQoverPtResVsTrueEta_[fitName] = inputDir.make<TProfile>(
          addn("QoverPtResVsTrueEta"), "q/p_{T} resolution; |#eta|; q/p_{T} resolution", 30, 0.0, 3.0);
      hisPhi0ResVsTrueEta_[fitName] = inputDir.make<TProfile>(
          addn("PhiResVsTrueEta"), "#phi_{0} resolution; |#eta|; #phi_{0} resolution", 30, 0.0, 3.0);
      hisEtaResVsTrueEta_[fitName] =
          inputDir.make<TProfile>(addn("EtaResVsTrueEta"), "#eta resolution; |#eta|; #eta resolution", 30, 0.0, 3.0);
      hisZ0ResVsTrueEta_[fitName] =
          inputDir.make<TProfile>(addn("Z0ResVsTrueEta"), "z_{0} resolution; |#eta|; z_{0} resolution", 30, 0.0, 3.0);
      hisD0ResVsTrueEta_[fitName] =
          inputDir.make<TProfile>(addn("D0ResVsTrueEta"), "d_{0} resolution; |#eta|; d_{0} resolution", 30, 0.0, 3.0);

      hisQoverPtResVsTrueInvPt_[fitName] = inputDir.make<TProfile>(
          addn("QoverPtResVsTrueInvPt"), "q/p_{T} resolution; 1/p_{T}; q/p_{T} resolution", 25, 0.0, 0.5);
      hisPhi0ResVsTrueInvPt_[fitName] = inputDir.make<TProfile>(
          addn("PhiResVsTrueInvPt"), "#phi_{0} resolution; 1/p_{T}; #phi_{0} resolution", 25, 0.0, 0.5);
      hisEtaResVsTrueInvPt_[fitName] =
          inputDir.make<TProfile>(addn("EtaResVsTrueInvPt"), "#eta resolution; 1/p_{T}; #eta resolution", 25, 0.0, 0.5);
      hisZ0ResVsTrueInvPt_[fitName] = inputDir.make<TProfile>(
          addn("Z0ResVsTrueInvPt"), "z_{0} resolution; 1/p_{T}; z_{0} resolution", 25, 0.0, 0.5);
      hisD0ResVsTrueInvPt_[fitName] = inputDir.make<TProfile>(
          addn("D0ResVsTrueInvPt"), "d_{0} resolution; 1/p_{T}; d_{0} resolution", 25, 0.0, 0.5);

      // Duplicate track histos.
      profDupFitTrksVsEta_[fitName] =
          inputDir.make<TProfile>(addn("DupFitTrksVsEta"), "; #eta; No. of duplicate tracks per TP", 12, 0., 3.);
      profDupFitTrksVsInvPt_[fitName] =
          inputDir.make<TProfile>(addn("DupFitTrksVsInvPt"), "; 1/Pt; No. of duplicate tracks per TP", 25, 0., 0.5);

      // Histos for tracking efficiency vs. TP kinematics. (Binning must match similar histos in bookTrackCands()).
      hisFitTPinvptForEff_[fitName] =
          inputDir.make<TH1F>(addn("FitTPinvptForEff"), "; TP 1/Pt of fitted tracks (for effi.);", 50, 0., 0.5);
      hisFitTPetaForEff_[fitName] =
          inputDir.make<TH1F>(addn("FitTPetaForEff"), "; TP #eta of fitted tracks (for effi.);", 20, -3., 3.);
      hisFitTPphiForEff_[fitName] =
          inputDir.make<TH1F>(addn("FitTPphiForEff"), "; TP #phi of fitted tracks (for effi.);", 20, -M_PI, M_PI);

      // Histo for efficiency to reconstruct track perfectly (no incorrect hits). (Binning must match similar histos in bookTrackCands()).
      hisPerfFitTPinvptForEff_[fitName] = inputDir.make<TH1F>(
          addn("PerfFitTPinvptForEff"), "; TP 1/Pt of fitted tracks (for perf. effi.);", 50, 0., 0.5);
      hisPerfFitTPetaForEff_[fitName] = inputDir.make<TH1F>(
          addn("PerfFitTPetaForEff"), "; TP #eta of fitted tracks (for perfect effi.);", 20, -3., 3.);

      // Histos for tracking efficiency vs. TP production point. (Binning must match similar histos in bookTrackCands()).
      hisFitTPd0ForEff_[fitName] =
          inputDir.make<TH1F>(addn("FitTPd0ForEff"), "; TP d0 of fitted tracks (for effi.);", 40, 0., 4.);
      hisFitTPz0ForEff_[fitName] =
          inputDir.make<TH1F>(addn("FitTPz0ForEff"), "; TP z0 of fitted tracks (for effi.);", 50, 0., 25.);

      // Histos for algorithmic tracking efficiency vs. TP kinematics. (Binning must match similar histos in bookTrackCands()).
      hisFitTPinvptForAlgEff_[fitName] =
          inputDir.make<TH1F>(addn("FitTPinvptForAlgEff"), "; TP 1/Pt of fitted tracks (for alg. effi.);", 50, 0., 0.5);
      hisFitTPetaForAlgEff_[fitName] =
          inputDir.make<TH1F>(addn("FitTPetaForAlgEff"), "; TP #eta of fitted tracks (for alg. effi.);", 20, -3., 3.);
      hisFitTPphiForAlgEff_[fitName] = inputDir.make<TH1F>(
          addn("FitTPphiForAlgEff"), "; TP #phi of fitted tracks (for alg. effi.);", 20, -M_PI, M_PI);

      // Histo for efficiency to reconstruct track perfectly (no incorrect hits). (Binning must match similar histos in bookTrackCands()).
      hisPerfFitTPinvptForAlgEff_[fitName] = inputDir.make<TH1F>(
          addn("PerfFitTPinvptForAlgEff"), "; TP 1/Pt of fitted tracks (for perf. alg. effi.);", 50, 0., 0.5);
      hisPerfFitTPetaForAlgEff_[fitName] =
          inputDir.make<TH1F>(addn("PerfFitTPetaForAlgEff"), "; TP #eta (for perf. alg. effi.);", 20, -3., 3.);

      // Histos for algorithmic tracking efficiency vs. TP production point. (Binning must match similar histos in bookTrackCands()).
      hisFitTPd0ForAlgEff_[fitName] =
          inputDir.make<TH1F>(addn("FitTPd0ForAlgEff"), "; TP d0 of fitted tracks (for alg. effi.);", 40, 0., 4.);
      hisFitTPz0ForAlgEff_[fitName] =
          inputDir.make<TH1F>(addn("FitTPz0ForAlgEff"), "; TP z0 of fitted tracks (for alg. effi.);", 50, 0., 25.);
    }
    return inputDirMap;
  }

virtual bool tmtt::Histos::enabled ( ) const

inlinevirtual

Definition at line 60 of file Histos.h.

References available(), tmtt::Settings::enableHistos(), and settings_.

Referenced by book(), endJobAnalysis(), fill(), and trackerGeometryAnalysis().

{ return (settings_->enableHistos() && available()); }

void Histos::endJobAnalysis ( const HTrphi::ErrorMonitor *  htRphiErrMon = nullptr )

virtual

Definition at line 1563 of file Histos.cc.

References bApproxMistake_, enabled(), DivergingColor::frac, hisRecoTPinvptForEff_, tmtt::Settings::hybrid(), phase1PixelTopology::layer, mapBarrelLayerMaxR_, mapBarrelLayerMinR_, mapEndcapWheelMaxZ_, mapEndcapWheelMinZ_, mapExtraAModuleTypeMaxR_, mapExtraAModuleTypeMaxZ_, mapExtraAModuleTypeMinR_, mapExtraAModuleTypeMinZ_, mapExtraBModuleTypeMaxR_, mapExtraBModuleTypeMaxZ_, mapExtraBModuleTypeMinR_, mapExtraBModuleTypeMinZ_, mapExtraCModuleTypeMaxR_, mapExtraCModuleTypeMaxZ_, mapExtraCModuleTypeMinR_, mapExtraCModuleTypeMinZ_, mapExtraDModuleTypeMaxR_, mapExtraDModuleTypeMaxZ_, mapExtraDModuleTypeMinR_, mapExtraDModuleTypeMinZ_, mapModuleTypeMaxR_, mapModuleTypeMaxZ_, mapModuleTypeMinR_, mapModuleTypeMinZ_, tmtt::HTrphi::ErrorMonitor::maxLineGradient, tmtt::HTrphi::ErrorMonitor::numErrorsNorm, tmtt::HTrphi::ErrorMonitor::numErrorsTypeA, tmtt::HTrphi::ErrorMonitor::numErrorsTypeB, oldSumW2opt_, AlCaHLTBitMon_ParallelJobs::p, plotHybridDupRemovalEfficiency(), plotTrackEffAfterFit(), plotTrackEfficiency(), plotTrackletSeedEfficiency(), printFitTrackPerformance(), printHybridDupRemovalPerformance(), printTrackletSeedFindingPerformance(), printTrackPerformance(), ranRZfilter_, settings_, and trackFitters_.

Referenced by tmtt::TMTrackProducer::globalEndJob().

                                                                    {
    // Don't bother producing summary if user didn't request histograms via TFileService in their cfg.
    if (not this->enabled())
      return;

    // Protection when running in wierd mixed hybrid-TMTT modes.
    bool wierdMixedMode = (hisRecoTPinvptForEff_.find("TRACKLET") == hisRecoTPinvptForEff_.end());

    if (settings_->hybrid() && not wierdMixedMode) {
      // Produce plots of tracking efficieny after tracklet pattern reco.
      this->plotTrackletSeedEfficiency();
      this->plotTrackEfficiency("TRACKLET");
      this->plotHybridDupRemovalEfficiency();

    } else {
      // Produce plots of tracking efficiency using track candidates found after HT.
      this->plotTrackEfficiency("HT");

      // Optionally produce plots of tracking efficiency using track candidates found after r-z track filter.
      if (ranRZfilter_)
        this->plotTrackEfficiency("RZ");
    }

    // Produce more plots of tracking efficiency using track candidates after track fit.
    for (auto& fitName : trackFitters_) {
      this->plotTrackEffAfterFit(fitName);
    }

    PrintL1trk() << "=========================================================================";

    // Print r (z) range in which each barrel layer (endcap wheel) appears.
    // (Needed by firmware).
    PrintL1trk();
    PrintL1trk() << "--- r range in which stubs in each barrel layer appear ---";
    for (const auto& p : mapBarrelLayerMinR_) {
      unsigned int layer = p.first;
      PrintL1trk() << "   layer = " << layer << " : " << mapBarrelLayerMinR_[layer] << " < r < "
                   << mapBarrelLayerMaxR_[layer];
    }
    PrintL1trk() << "--- |z| range in which stubs in each endcap wheel appear ---";
    for (const auto& p : mapEndcapWheelMinZ_) {
      unsigned int layer = p.first;
      PrintL1trk() << "   wheel = " << layer << " : " << mapEndcapWheelMinZ_[layer] << " < |z| < "
                   << mapEndcapWheelMaxZ_[layer];
    }

    // Print (r,|z|) range in which each module type (defined in DigitalStub) appears.
    // (Needed by firmware).
    PrintL1trk();
    PrintL1trk() << "--- (r,|z|) range in which each module type (defined in DigitalStub) appears ---";
    for (const auto& p : mapModuleTypeMinR_) {
      unsigned int modType = p.first;
      PrintL1trk() << "   Module type = " << modType << setprecision(1) << " : r range = ("
                   << mapModuleTypeMinR_[modType] << "," << mapModuleTypeMaxR_[modType] << "); z range = ("
                   << mapModuleTypeMinZ_[modType] << "," << mapModuleTypeMaxZ_[modType] << ")";
    }
    // Ugly bodge to allow for modules in barrel layers 1-2 & endcap wheels 3-5 being different.
    PrintL1trk() << "and in addition";
    for (const auto& p : mapExtraAModuleTypeMinR_) {
      unsigned int modType = p.first;
      PrintL1trk() << "   Module type = " << modType << setprecision(1) << " : r range = ("
                   << mapExtraAModuleTypeMinR_[modType] << "," << mapExtraAModuleTypeMaxR_[modType] << "); z range = ("
                   << mapExtraAModuleTypeMinZ_[modType] << "," << mapExtraAModuleTypeMaxZ_[modType] << ")";
    }
    PrintL1trk() << "and in addition";
    for (const auto& p : mapExtraBModuleTypeMinR_) {
      unsigned int modType = p.first;
      PrintL1trk() << "   Module type = " << modType << setprecision(1) << " : r range = ("
                   << mapExtraBModuleTypeMinR_[modType] << "," << mapExtraBModuleTypeMaxR_[modType] << "); z range = ("
                   << mapExtraBModuleTypeMinZ_[modType] << "," << mapExtraBModuleTypeMaxZ_[modType] << ")";
    }
    PrintL1trk() << "and in addition";
    for (const auto& p : mapExtraCModuleTypeMinR_) {
      unsigned int modType = p.first;
      PrintL1trk() << "   Module type = " << modType << setprecision(1) << " : r range = ("
                   << mapExtraCModuleTypeMinR_[modType] << "," << mapExtraCModuleTypeMaxR_[modType] << "); z range = ("
                   << mapExtraCModuleTypeMinZ_[modType] << "," << mapExtraCModuleTypeMaxZ_[modType] << ")";
    }
    PrintL1trk() << "and in addition";
    for (const auto& p : mapExtraDModuleTypeMinR_) {
      unsigned int modType = p.first;
      PrintL1trk() << "   Module type = " << modType << setprecision(1) << " : r range = ("
                   << mapExtraDModuleTypeMinR_[modType] << "," << mapExtraDModuleTypeMaxR_[modType] << "); z range = ("
                   << mapExtraDModuleTypeMinZ_[modType] << "," << mapExtraDModuleTypeMaxZ_[modType] << ")";
    }
    PrintL1trk();

    if (settings_->hybrid() && not wierdMixedMode) {
      //--- Print summary of tracklet pattern reco
      this->printTrackletSeedFindingPerformance();
      this->printTrackPerformance("TRACKLET");
      this->printHybridDupRemovalPerformance();
    } else {
      //--- Print summary of track-finding performance after HT
      this->printTrackPerformance("HT");
      //--- Optionally print summary of track-finding performance after r-z track filter.
      if (ranRZfilter_)
        this->printTrackPerformance("RZ");
    }

    //--- Print summary of track-finding performance after helix fit, for each track fitting algorithm used.
    for (const string& fitName : trackFitters_) {
      this->printFitTrackPerformance(fitName);
    }
    PrintL1trk() << "=========================================================================";

    if (htRphiErrMon != nullptr && not settings_->hybrid()) {
      // Check that stub filling was consistent with known limitations of HT firmware design.

      PrintL1trk() << "Max. |gradients| of stub lines in HT array is: r-phi = " << htRphiErrMon->maxLineGradient;

      if (htRphiErrMon->maxLineGradient > 1.) {
        PrintL1trk()
            << "WARNING: Line |gradient| exceeds 1, which firmware will not be able to cope with! Please adjust HT "
               "array size to avoid this.";

      } else if (htRphiErrMon->numErrorsTypeA > 0.) {
        float frac = float(htRphiErrMon->numErrorsTypeA) / float(htRphiErrMon->numErrorsNorm);
        PrintL1trk()
            << "WARNING: Despite line gradients being less than one, some fraction of HT columns have filled cells "
               "with no filled neighbours in W, SW or NW direction. Firmware will object to this! ";
        PrintL1trk() << "This fraction = " << frac << " for r-phi HT";

      } else if (htRphiErrMon->numErrorsTypeB > 0.) {
        float frac = float(htRphiErrMon->numErrorsTypeB) / float(htRphiErrMon->numErrorsNorm);
        PrintL1trk()
            << "WARNING: Despite line gradients being less than one, some fraction of HT columns recorded individual "
               "stubs being added to more than two cells! Thomas firmware will object to this! ";
        PrintL1trk() << "This fraction = " << frac << " for r-phi HT";
      }
    }

    // Check if GP B approximation cfg params are inconsistent.
    if (bApproxMistake_)
      PrintL1trk() << "\n WARNING: BApprox cfg params are inconsistent - see printout above.";

    // Restore original ROOT default cfg.
    TH1::SetDefaultSumw2(oldSumW2opt_);
  }

void Histos::fill ( const InputData &  inputData, const Array2D< std::unique_ptr< Sector >> &  mSectors, const Array2D< std::unique_ptr< HTrphi >> &  mHtPhis, const Array2D< std::unique_ptr< Make3Dtracks >> &  mGet3Dtrks, const std::map< std::string, std::list< const L1fittedTrack * >> &  mapFinalTracks  )

virtual

Definition at line 83 of file Histos.cc.

References enabled(), fillEtaPhiSectors(), fillInputData(), fillRphiHT(), fillRZfilters(), fillTrackCands(), fillTrackFitting(), and ranRZfilter_.

Referenced by tmtt::TMTrackProducer::produce().

                                                                                            {
    // Each function here protected by a mytex lock, so only one thread can run it at a time.

    // Don't bother filling histograms if user didn't request them via TFileService in their cfg.
    if (not this->enabled())
      return;

    // Fill histograms about input data.
    this->fillInputData(inputData);

    // Fill histograms checking if (eta,phi) sector definition choices are good.
    this->fillEtaPhiSectors(inputData, mSectors);

    // Fill histograms checking filling of r-phi HT array.
    this->fillRphiHT(mHtRphis);

    // Fill histograms about r-z track filters.
    if (ranRZfilter_)
      this->fillRZfilters(mMake3Dtrks);

    // Fill histograms studying 3D track candidates found after HT.
    this->fillTrackCands(inputData, mMake3Dtrks, "HT");

    // Fill histograms studying 3D track candidates found after r-z track filter.
    if (ranRZfilter_) {
      this->fillTrackCands(inputData, mMake3Dtrks, "RZ");
    }

    // Fill histograms studying track fitting performance
    this->fillTrackFitting(inputData, mapFinalTracks);
  }

void Histos::fillEtaPhiSectors ( const InputData &  inputData, const Array2D< std::unique_ptr< Sector >> &  mSectors  )

protectedvirtual

Definition at line 386 of file Histos.cc.

References hisNumEtaSecsPerStub_, hisNumPhiSecsPerStub_, hisNumStubsPerSec_, tmtt::Sector::inside(), tmtt::Sector::insideEta(), tmtt::Sector::insidePhi(), mutex, numEtaRegions_, numPhiSectors_, and tmtt::InputData::stubsConst().

Referenced by fill().

                                                                                                        {
    // Allow only one thread to run this function at a time
    static std::mutex myMutex;
    std::lock_guard<std::mutex> myGuard(myMutex);

    const list<const Stub*>& vStubs = inputData.stubsConst();
    //const list<TP>& vTPs = inputData.getTPs();

    //=== Loop over all stubs, counting how many sectors each one appears in.

    for (const Stub* stub : vStubs) {
      // Number of (eta,phi), phi & eta sectors containing this stub.
      unsigned int nEtaSecs = 0;
      unsigned int nPhiSecs = 0;

      // Loop over (eta, phi) sectors.
      for (unsigned int iPhiSec = 0; iPhiSec < numPhiSectors_; iPhiSec++) {
        for (unsigned int iEtaReg = 0; iEtaReg < numEtaRegions_; iEtaReg++) {
          const Sector* sector = mSectors(iPhiSec, iEtaReg).get();

          // Check if sector contains stub stub, and if so count it.
          // Take care to just use one eta (phi) typical region when counting phi (eta) sectors.
          if (iPhiSec == 0 && sector->insideEta(stub))
            nEtaSecs++;
          if (iEtaReg == 0 && sector->insidePhi(stub))
            nPhiSecs++;
        }
      }

      // Plot number of sectors each stub appears in.
      hisNumEtaSecsPerStub_->Fill(nEtaSecs);
      hisNumPhiSecsPerStub_->Fill(nPhiSecs);
    }

    //=== Loop over all sectors, counting the stubs in each one.
    for (unsigned int iEtaReg = 0; iEtaReg < numEtaRegions_; iEtaReg++) {
      for (unsigned int iPhiSec = 0; iPhiSec < numPhiSectors_; iPhiSec++) {
        const Sector* sector = mSectors(iPhiSec, iEtaReg).get();

        unsigned int nStubs = 0;
        for (const Stub* stub : vStubs) {
          if (sector->inside(stub))
            nStubs++;
        }
        hisNumStubsPerSec_->Fill(nStubs);
      }
    }
  }

void Histos::fillInputData ( const InputData &  inputData )

protectedvirtual

Definition at line 164 of file Histos.cc.

References funct::abs(), tmtt::InputData::allStubs(), tmtt::Utility::countLayers(), tmtt::InputData::getTPs(), hisBendResStub_, hisBendStub_, hisNumLayersPerTP_, hisNumPSLayersPerTP_, hisStubIneffiVsEta_, hisStubIneffiVsInvPt_, hisStubKillFE_, hisStubsVsEta_, hisStubsVsR_, hisTPd0ForAlgEff_, hisTPd0ForEff_, hisTPetaForAlgEff_, hisTPetaForEff_, hisTPinvptForAlgEff_, hisTPinvptForEff_, hisTPphiForAlgEff_, hisTPphiForEff_, hisTPz0ForAlgEff_, hisTPz0ForEff_, phase1PixelTopology::layer, mapBarrelLayerMaxR_, mapBarrelLayerMinR_, mapEndcapWheelMaxZ_, mapEndcapWheelMinZ_, mapExtraAModuleTypeMaxR_, mapExtraAModuleTypeMaxZ_, mapExtraAModuleTypeMinR_, mapExtraAModuleTypeMinZ_, mapExtraBModuleTypeMaxR_, mapExtraBModuleTypeMaxZ_, mapExtraBModuleTypeMinR_, mapExtraBModuleTypeMinZ_, mapExtraCModuleTypeMaxR_, mapExtraCModuleTypeMaxZ_, mapExtraCModuleTypeMinR_, mapExtraCModuleTypeMinZ_, mapExtraDModuleTypeMaxR_, mapExtraDModuleTypeMaxZ_, mapExtraDModuleTypeMinR_, mapExtraDModuleTypeMinZ_, mapModuleTypeMaxR_, mapModuleTypeMaxZ_, mapModuleTypeMinR_, mapModuleTypeMinZ_, mutex, alignCSCRings::r, alignCSCRings::s, settings_, tmtt::InputData::stubsConst(), and cmsswSequenceInfo::tp.

Referenced by fill().

                                                       {
    // Allow only one thread to run this function at a time
    static std::mutex myMutex;
    std::lock_guard<std::mutex> myGuard(myMutex);

    const list<const Stub*>& vStubs = inputData.stubsConst();
    const list<TP>& vTPs = inputData.getTPs();

    for (const Stub* stub : vStubs) {
      hisStubsVsEta_->Fill(stub->eta());
      hisStubsVsR_->Fill(stub->r());
    }

    // Study efficiency of stubs to pass front-end electronics cuts.

    const list<Stub>& vAllStubs = inputData.allStubs();  // Get all stubs prior to FE cuts to do this.
    for (const Stub& s : vAllStubs) {
      unsigned int layerOrTenPlusRing = s.barrel() ? s.layerId() : 10 + s.trackerModule()->endcapRing();
      // Fraction of all stubs (good and bad) failing tightened front-end electronics cuts.
      hisStubKillFE_->Fill(layerOrTenPlusRing, (!s.frontendPass()));
    }

    // Study efficiency for good stubs of tightened front end-electronics cuts.
    for (const TP& tp : vTPs) {
      if (tp.useForAlgEff()) {  // Only bother for stubs that are on TP that we have a chance of reconstructing.
        const vector<const Stub*>& stubs = tp.assocStubs();
        for (const Stub* s : stubs) {
          hisStubIneffiVsInvPt_->Fill(1. / tp.pt(), (!s->frontendPass()));
          hisStubIneffiVsEta_->Fill(std::abs(tp.eta()), (!s->frontendPass()));
        }
      }
    }

    // Plot stub bend-derived information.
    for (const Stub* stub : vStubs) {
      hisBendStub_->Fill(stub->bend());
    }

    // Look at stub resolution.
    for (const TP& tp : vTPs) {
      if (tp.useForAlgEff()) {
        const vector<const Stub*>& assStubs = tp.assocStubs();

        for (const Stub* stub : assStubs) {
          hisBendResStub_->Fill(stub->bend() - tp.dphi(stub->r()) / stub->dphiOverBend());
        }

        if (std::abs(tp.eta()) < 0.5) {
          double nLayersOnTP = Utility::countLayers(settings_, assStubs, true, false);
          double nPSLayersOnTP = Utility::countLayers(settings_, assStubs, true, true);
          hisNumLayersPerTP_->Fill(nLayersOnTP);
          hisNumPSLayersPerTP_->Fill(nPSLayersOnTP);
        }
      }
    }

    // Determine r (z) range of each barrel layer (endcap wheel).

    for (const Stub* stub : vStubs) {
      unsigned int layer = stub->layerId();
      if (stub->barrel()) {
        // Get range in r of each barrel layer.
        float r = stub->r();
        if (mapBarrelLayerMinR_.find(layer) == mapBarrelLayerMinR_.end()) {
          mapBarrelLayerMinR_[layer] = r;
          mapBarrelLayerMaxR_[layer] = r;
        } else {
          if (mapBarrelLayerMinR_[layer] > r)
            mapBarrelLayerMinR_[layer] = r;
          if (mapBarrelLayerMaxR_[layer] < r)
            mapBarrelLayerMaxR_[layer] = r;
        }
      } else {
        layer = layer % 10;
        // Range in |z| of each endcap wheel.
        float z = std::abs(stub->z());
        if (mapEndcapWheelMinZ_.find(layer) == mapEndcapWheelMinZ_.end()) {
          mapEndcapWheelMinZ_[layer] = z;
          mapEndcapWheelMaxZ_[layer] = z;
        } else {
          if (mapEndcapWheelMinZ_[layer] > z)
            mapEndcapWheelMinZ_[layer] = z;
          if (mapEndcapWheelMaxZ_[layer] < z)
            mapEndcapWheelMaxZ_[layer] = z;
        }
      }
    }

    // Determine Range in (r,|z|) of each module type.

    for (const Stub* stub : vStubs) {
      float r = stub->r();
      float z = std::abs(stub->z());
      unsigned int modType = stub->trackerModule()->moduleTypeID();
      // Do something ugly, as modules in 1-2nd & 3-4th endcap wheels are different to those in wheel 5 ...
      // And boundary between flat & tilted modules in barrel layers 1-3 varies in z.
      if (stub->barrel() && stub->layerId() == 1) {  // barrel layer 1
        if (mapExtraAModuleTypeMinR_.find(modType) == mapExtraAModuleTypeMinR_.end()) {
          mapExtraAModuleTypeMinR_[modType] = r;
          mapExtraAModuleTypeMaxR_[modType] = r;
          mapExtraAModuleTypeMinZ_[modType] = z;
          mapExtraAModuleTypeMaxZ_[modType] = z;
        } else {
          if (mapExtraAModuleTypeMinR_[modType] > r)
            mapExtraAModuleTypeMinR_[modType] = r;
          if (mapExtraAModuleTypeMaxR_[modType] < r)
            mapExtraAModuleTypeMaxR_[modType] = r;
          if (mapExtraAModuleTypeMinZ_[modType] > z)
            mapExtraAModuleTypeMinZ_[modType] = z;
          if (mapExtraAModuleTypeMaxZ_[modType] < z)
            mapExtraAModuleTypeMaxZ_[modType] = z;
        }
      } else if (stub->barrel() && stub->layerId() == 2) {  // barrel layer 2
        if (mapExtraBModuleTypeMinR_.find(modType) == mapExtraBModuleTypeMinR_.end()) {
          mapExtraBModuleTypeMinR_[modType] = r;
          mapExtraBModuleTypeMaxR_[modType] = r;
          mapExtraBModuleTypeMinZ_[modType] = z;
          mapExtraBModuleTypeMaxZ_[modType] = z;
        } else {
          if (mapExtraBModuleTypeMinR_[modType] > r)
            mapExtraBModuleTypeMinR_[modType] = r;
          if (mapExtraBModuleTypeMaxR_[modType] < r)
            mapExtraBModuleTypeMaxR_[modType] = r;
          if (mapExtraBModuleTypeMinZ_[modType] > z)
            mapExtraBModuleTypeMinZ_[modType] = z;
          if (mapExtraBModuleTypeMaxZ_[modType] < z)
            mapExtraBModuleTypeMaxZ_[modType] = z;
        }
      } else if (!stub->barrel() && (stub->layerId() % 10 == 1 || stub->layerId() % 10 == 2)) {  // endcap wheel 1-2
        if (mapExtraCModuleTypeMinR_.find(modType) == mapExtraCModuleTypeMinR_.end()) {
          mapExtraCModuleTypeMinR_[modType] = r;
          mapExtraCModuleTypeMaxR_[modType] = r;
          mapExtraCModuleTypeMinZ_[modType] = z;
          mapExtraCModuleTypeMaxZ_[modType] = z;
        } else {
          if (mapExtraCModuleTypeMinR_[modType] > r)
            mapExtraCModuleTypeMinR_[modType] = r;
          if (mapExtraCModuleTypeMaxR_[modType] < r)
            mapExtraCModuleTypeMaxR_[modType] = r;
          if (mapExtraCModuleTypeMinZ_[modType] > z)
            mapExtraCModuleTypeMinZ_[modType] = z;
          if (mapExtraCModuleTypeMaxZ_[modType] < z)
            mapExtraCModuleTypeMaxZ_[modType] = z;
        }
      } else if (!stub->barrel() && (stub->layerId() % 10 == 3 || stub->layerId() % 10 == 4)) {  // endcap wheel 3-4
        if (mapExtraDModuleTypeMinR_.find(modType) == mapExtraDModuleTypeMinR_.end()) {
          mapExtraDModuleTypeMinR_[modType] = r;
          mapExtraDModuleTypeMaxR_[modType] = r;
          mapExtraDModuleTypeMinZ_[modType] = z;
          mapExtraDModuleTypeMaxZ_[modType] = z;
        } else {
          if (mapExtraDModuleTypeMinR_[modType] > r)
            mapExtraDModuleTypeMinR_[modType] = r;
          if (mapExtraDModuleTypeMaxR_[modType] < r)
            mapExtraDModuleTypeMaxR_[modType] = r;
          if (mapExtraDModuleTypeMinZ_[modType] > z)
            mapExtraDModuleTypeMinZ_[modType] = z;
          if (mapExtraDModuleTypeMaxZ_[modType] < z)
            mapExtraDModuleTypeMaxZ_[modType] = z;
        }
      } else {  // barrel layer 3-6 or endcap wheel 5.
        if (mapModuleTypeMinR_.find(modType) == mapModuleTypeMinR_.end()) {
          mapModuleTypeMinR_[modType] = r;
          mapModuleTypeMaxR_[modType] = r;
          mapModuleTypeMinZ_[modType] = z;
          mapModuleTypeMaxZ_[modType] = z;
        } else {
          if (mapModuleTypeMinR_[modType] > r)
            mapModuleTypeMinR_[modType] = r;
          if (mapModuleTypeMaxR_[modType] < r)
            mapModuleTypeMaxR_[modType] = r;
          if (mapModuleTypeMinZ_[modType] > z)
            mapModuleTypeMinZ_[modType] = z;
          if (mapModuleTypeMaxZ_[modType] < z)
            mapModuleTypeMaxZ_[modType] = z;
        }
      }
    }

    //=== Make denominator of tracking efficiency plots

    for (const TP& tp : vTPs) {
      if (tp.useForEff()) {  // Check TP is good for efficiency measurement.
        // Plot kinematics of all good TP.
        hisTPinvptForEff_->Fill(1. / tp.pt());
        hisTPetaForEff_->Fill(tp.eta());
        hisTPphiForEff_->Fill(tp.phi0());
        // Plot also production point of all good TP.
        hisTPd0ForEff_->Fill(std::abs(tp.d0()));
        hisTPz0ForEff_->Fill(std::abs(tp.z0()));

        if (tp.useForAlgEff()) {  // Check TP is good for algorithmic efficiency measurement.
          hisTPinvptForAlgEff_->Fill(1. / tp.pt());
          hisTPetaForAlgEff_->Fill(tp.eta());
          hisTPphiForAlgEff_->Fill(tp.phi0());
          // Plot also production point of all good TP.
          hisTPd0ForAlgEff_->Fill(std::abs(tp.d0()));
          hisTPz0ForAlgEff_->Fill(std::abs(tp.z0()));
        }
      }
    }
  }

void Histos::fillRphiHT ( const Array2D< std::unique_ptr< HTrphi >> &  mHtRphis )

protectedvirtual

Definition at line 445 of file Histos.cc.

Referenced by fill().

                                                                     {
    //--- Loop over (eta,phi) sectors, counting the number of stubs in the HT array of each.

    // Allow only one thread to run this function at a time (UNCOMMENT IF YOU ADD HISTOS HERE)
    //static std::mutex myMutex;
    //std::lock_guard<std::mutex> myGuard(myMutex);
  }

void Histos::fillRZfilters ( const Array2D< std::unique_ptr< Make3Dtracks >> &  mMake3Dtrks )

protectedvirtual

Definition at line 463 of file Histos.cc.

Referenced by fill().

                                                                                 {
    // Allow only one thread to run this function at a time (UNCOMMENT IF YOU ADD HISTOS HERE)
    //static std::mutex myMutex;
    //std::lock_guard<std::mutex> myGuard(myMutex);
  }

virtual void tmtt::Histos::fillTrackCands ( const InputData &  inputData, const Array2D< std::unique_ptr< Make3Dtracks >> &  mMake3Dtrks, const std::string &  tName  )

protectedvirtual

Referenced by fill().

virtual void tmtt::Histos::fillTrackCands ( const InputData &  inputData, const std::vector< L1track3D > &  tracks, const std::string &  tName  )

protectedvirtual

void Histos::fillTrackFitting ( const InputData &  inputData, const std::map< std::string, std::list< const L1fittedTrack * >> &  mapFinalTracks  )

protectedvirtual

Definition at line 1029 of file Histos.cc.

References funct::abs(), tmtt::TP::d0(), reco::deltaPhi(), tmtt::TP::eta(), tmtt::InputData::getTPs(), hisD0ResVsTrueEta_, hisD0ResVsTrueInvPt_, hisEtaResVsTrueEta_, hisEtaResVsTrueInvPt_, hisFitChi2DofRphiMatched_, hisFitChi2DofRphiUnmatched_, hisFitChi2DofRzMatched_, hisFitChi2DofRzUnmatched_, hisFitD0Matched_, hisFitD0Unmatched_, hisFitEtaMatched_, hisFitEtaUnmatched_, hisFitPhi0Matched_, hisFitPhi0Unmatched_, hisFitQinvPtMatched_, hisFitQinvPtUnmatched_, hisFitTPd0ForAlgEff_, hisFitTPd0ForEff_, hisFitTPetaForAlgEff_, hisFitTPetaForEff_, hisFitTPinvptForAlgEff_, hisFitTPinvptForEff_, hisFitTPphiForAlgEff_, hisFitTPphiForEff_, hisFitTPz0ForAlgEff_, hisFitTPz0ForEff_, hisFitZ0Matched_, hisFitZ0Unmatched_, hisKalmanChi2DofSkipLay0Matched_, hisKalmanChi2DofSkipLay0Unmatched_, hisKalmanChi2DofSkipLay1Matched_, hisKalmanChi2DofSkipLay1Unmatched_, hisKalmanChi2DofSkipLay2Matched_, hisKalmanChi2DofSkipLay2Unmatched_, hisKalmanNumUpdateCalls_, hisNumFitTrks_, hisNumFitTrksPerNon_, hisPerfFitTPetaForAlgEff_, hisPerfFitTPetaForEff_, hisPerfFitTPinvptForAlgEff_, hisPerfFitTPinvptForEff_, hisPhi0ResVsTrueEta_, hisPhi0ResVsTrueInvPt_, hisQoverPtResVsTrueEta_, hisQoverPtResVsTrueInvPt_, hisStubsPerFitTrack_, hisZ0ResVsTrueEta_, hisZ0ResVsTrueInvPt_, mutex, pileupDistInMC::num, tmtt::Settings::numPhiNonants(), numPhiSectors_, tmtt::TP::phi0(), profDupFitTrksVsEta_, profDupFitTrksVsInvPt_, profFitChi2DofRphiVsInvPtMatched_, profFitChi2DofRphiVsInvPtUnmatched_, profNumFitTracks_, profStubsOnFitTracks_, tmtt::TP::qOverPt(), settings_, cmsswSequenceInfo::tp, trackFitters_, tmtt::TP::useForAlgEff(), and tmtt::TP::z0().

Referenced by fill().

                                                                                               {
    // Allow only one thread to run this function at a time
    static std::mutex myMutex;
    std::lock_guard<std::mutex> myGuard(myMutex);

    const list<TP>& vTPs = inputData.getTPs();

    // Loop over all the fitting algorithms we are trying.
    for (const string& fitName : trackFitters_) {
      const list<const L1fittedTrack*>& fittedTracks = mapFinalTracks.at(fitName);  // Get fitted tracks.

      // Count tracks
      unsigned int nFitTracks = 0;
      unsigned int nFitsMatchingTP = 0;

      const unsigned int numPhiNonants = settings_->numPhiNonants();
      vector<unsigned int> nFitTracksPerNonant(numPhiNonants, 0);

      for (const L1fittedTrack* fitTrk : fittedTracks) {
        nFitTracks++;

        // Get matched truth particle, if any.
        const TP* tp = fitTrk->matchedTP();
        if (tp != nullptr)
          nFitsMatchingTP++;
        // Count fitted tracks per nonant.
        unsigned int iNonant = (numPhiSectors_ > 0) ? floor(fitTrk->iPhiSec() * numPhiNonants / (numPhiSectors_))
                                                    : 0;  // phi nonant number
        nFitTracksPerNonant[iNonant]++;
      }

      profNumFitTracks_[fitName]->Fill(1, nFitTracks);
      profNumFitTracks_[fitName]->Fill(2, nFitsMatchingTP);

      hisNumFitTrks_[fitName]->Fill(nFitTracks);
      for (const unsigned int& num : nFitTracksPerNonant) {
        hisNumFitTrksPerNon_[fitName]->Fill(num);
      }

      // Count stubs assigned to fitted tracks.
      unsigned int nTotStubs = 0;
      for (const L1fittedTrack* fitTrk : fittedTracks) {
        unsigned int nStubs = fitTrk->numStubs();
        hisStubsPerFitTrack_[fitName]->Fill(nStubs);
        nTotStubs += nStubs;
      }
      profStubsOnFitTracks_[fitName]->Fill(1., nTotStubs);

      // Note truth particles that are successfully fitted. And which give rise to duplicate tracks.

      map<const TP*, bool> tpRecoedMap;  // Note which truth particles were successfully fitted.
      map<const TP*, bool>
          tpPerfRecoedMap;  // Note which truth particles were successfully fitted with no incorrect hits.
      map<const TP*, unsigned int> tpRecoedDup;  // Note that this TP gave rise to duplicate tracks.
      for (const TP& tp : vTPs) {
        tpRecoedMap[&tp] = false;
        tpPerfRecoedMap[&tp] = false;
        unsigned int nMatch = 0;
        for (const L1fittedTrack* fitTrk : fittedTracks) {
          const TP* assocTP = fitTrk->matchedTP();  // Get the TP the fitted track matches to, if any.
          if (assocTP == &tp) {
            tpRecoedMap[&tp] = true;
            if (fitTrk->purity() == 1.)
              tpPerfRecoedMap[&tp] = true;
            nMatch++;
          }
        }
        tpRecoedDup[&tp] = nMatch;
      }

      // Count truth particles that are successfully fitted.

      unsigned int nFittedTPs = 0;
      unsigned int nFittedTPsForEff = 0;
      for (const TP& tp : vTPs) {
        if (tpRecoedMap[&tp]) {  // Was this truth particle successfully fitted?
          nFittedTPs++;
          if (tp.useForEff())
            nFittedTPsForEff++;
        }
      }

      profNumFitTracks_[fitName]->Fill(6, nFittedTPs);
      profNumFitTracks_[fitName]->Fill(7, nFittedTPsForEff);

      // Loop over fitted tracks again.

      for (const L1fittedTrack* fitTrk : fittedTracks) {
        // Info for specific track fit algorithms.
        unsigned int nSkippedLayers = 0;
        unsigned int numUpdateCalls = 0;
        if (fitName.find("KF") != string::npos) {
          fitTrk->infoKF(nSkippedLayers, numUpdateCalls);
          hisKalmanNumUpdateCalls_[fitName]->Fill(numUpdateCalls);
        }

        //--- Compare fitted tracks that match truth particles to those that don't.

        // Get matched truth particle, if any.
        const TP* tp = fitTrk->matchedTP();

        if (tp != nullptr) {
          hisFitQinvPtMatched_[fitName]->Fill(fitTrk->qOverPt());
          hisFitPhi0Matched_[fitName]->Fill(fitTrk->phi0());
          hisFitD0Matched_[fitName]->Fill(fitTrk->d0());
          hisFitZ0Matched_[fitName]->Fill(fitTrk->z0());
          hisFitEtaMatched_[fitName]->Fill(fitTrk->eta());

          // Only plot matched chi2 for tracks with no incorrect stubs.
          if (fitTrk->purity() == 1.) {
            hisFitChi2DofRphiMatched_[fitName]->Fill(fitTrk->chi2rphi() / fitTrk->numDOFrphi());
            hisFitChi2DofRzMatched_[fitName]->Fill(fitTrk->chi2rz() / fitTrk->numDOFrz());
            profFitChi2DofRphiVsInvPtMatched_[fitName]->Fill(std::abs(fitTrk->qOverPt()),
                                                             (fitTrk->chi2rphi() / fitTrk->numDOFrphi()));

            if (fitName.find("KF") != string::npos) {
              // No. of skipped layers on track during Kalman track fit.
              if (nSkippedLayers == 0) {
                hisKalmanChi2DofSkipLay0Matched_[fitName]->Fill(fitTrk->chi2dof());
              } else if (nSkippedLayers == 1) {
                hisKalmanChi2DofSkipLay1Matched_[fitName]->Fill(fitTrk->chi2dof());
              } else if (nSkippedLayers >= 2) {
                hisKalmanChi2DofSkipLay2Matched_[fitName]->Fill(fitTrk->chi2dof());
              }
            }
          }

        } else {
          hisFitQinvPtUnmatched_[fitName]->Fill(fitTrk->qOverPt());
          hisFitPhi0Unmatched_[fitName]->Fill(fitTrk->phi0());
          hisFitD0Unmatched_[fitName]->Fill(fitTrk->d0());
          hisFitZ0Unmatched_[fitName]->Fill(fitTrk->z0());
          hisFitEtaUnmatched_[fitName]->Fill(fitTrk->eta());

          hisFitChi2DofRphiUnmatched_[fitName]->Fill(fitTrk->chi2rphi() / fitTrk->numDOFrphi());
          hisFitChi2DofRzUnmatched_[fitName]->Fill(fitTrk->chi2rz() / fitTrk->numDOFrz());
          profFitChi2DofRphiVsInvPtUnmatched_[fitName]->Fill(std::abs(fitTrk->qOverPt()),
                                                             (fitTrk->chi2rphi() / fitTrk->numDOFrphi()));

          if (fitName.find("KF") != string::npos) {
            // No. of skipped layers on track during Kalman track fit.
            if (nSkippedLayers == 0) {
              hisKalmanChi2DofSkipLay0Unmatched_[fitName]->Fill(fitTrk->chi2dof());
            } else if (nSkippedLayers == 1) {
              hisKalmanChi2DofSkipLay1Unmatched_[fitName]->Fill(fitTrk->chi2dof());
            } else if (nSkippedLayers >= 2) {
              hisKalmanChi2DofSkipLay2Unmatched_[fitName]->Fill(fitTrk->chi2dof());
            }
          }
        }
      }

      // Study helix param resolution.

      for (const L1fittedTrack* fitTrk : fittedTracks) {
        const TP* tp = fitTrk->matchedTP();
        if (tp != nullptr) {
          // IRT
          if ((resPlotOpt_ && tp->useForAlgEff()) ||
              (not resPlotOpt_)) {  // Check TP is good for efficiency measurement (& also comes from signal event if requested)

            // Plot helix parameter resolution against eta or Pt.
            hisQoverPtResVsTrueEta_[fitName]->Fill(std::abs(tp->eta()), std::abs(fitTrk->qOverPt() - tp->qOverPt()));
            hisPhi0ResVsTrueEta_[fitName]->Fill(std::abs(tp->eta()),
                                                std::abs(reco::deltaPhi(fitTrk->phi0(), tp->phi0())));
            hisEtaResVsTrueEta_[fitName]->Fill(std::abs(tp->eta()), std::abs(fitTrk->eta() - tp->eta()));
            hisZ0ResVsTrueEta_[fitName]->Fill(std::abs(tp->eta()), std::abs(fitTrk->z0() - tp->z0()));
            hisD0ResVsTrueEta_[fitName]->Fill(std::abs(tp->eta()), std::abs(fitTrk->d0() - tp->d0()));

            hisQoverPtResVsTrueInvPt_[fitName]->Fill(std::abs(tp->qOverPt()),
                                                     std::abs(fitTrk->qOverPt() - tp->qOverPt()));
            hisPhi0ResVsTrueInvPt_[fitName]->Fill(std::abs(tp->qOverPt()),
                                                  std::abs(reco::deltaPhi(fitTrk->phi0(), tp->phi0())));
            hisEtaResVsTrueInvPt_[fitName]->Fill(std::abs(tp->qOverPt()), std::abs(fitTrk->eta() - tp->eta()));
            hisZ0ResVsTrueInvPt_[fitName]->Fill(std::abs(tp->qOverPt()), std::abs(fitTrk->z0() - tp->z0()));
            hisD0ResVsTrueInvPt_[fitName]->Fill(std::abs(tp->qOverPt()), std::abs(fitTrk->d0() - tp->d0()));
          }
        }
      }

      //=== Study duplicate tracks.

      for (const TP& tp : vTPs) {
        if (tpRecoedMap[&tp]) {  // Was this truth particle successfully fitted?
          profDupFitTrksVsEta_[fitName]->Fill(std::abs(tp.eta()), tpRecoedDup[&tp] - 1);
          profDupFitTrksVsInvPt_[fitName]->Fill(std::abs(tp.qOverPt()), tpRecoedDup[&tp] - 1);
        }
      }

      //=== Study tracking efficiency by looping over tracking particles.

      for (const TP& tp : vTPs) {
        if (tp.useForEff()) {  // Check TP is good for efficiency measurement.

          // If TP was reconstucted by HT, then plot its kinematics.
          if (tpRecoedMap[&tp]) {  // This truth particle was successfully fitted.
            hisFitTPinvptForEff_[fitName]->Fill(1. / tp.pt());
            hisFitTPetaForEff_[fitName]->Fill(tp.eta());
            hisFitTPphiForEff_[fitName]->Fill(tp.phi0());
            // Plot also production point of all good reconstructed TP.
            hisFitTPd0ForEff_[fitName]->Fill(std::abs(tp.d0()));
            hisFitTPz0ForEff_[fitName]->Fill(std::abs(tp.z0()));
            // Also plot efficiency to perfectly reconstruct the track (no fake hits)
            if (tpPerfRecoedMap[&tp]) {  // This truth particle was successfully fitted with no incorrect hits.
              hisPerfFitTPinvptForEff_[fitName]->Fill(1. / tp.pt());
              hisPerfFitTPetaForEff_[fitName]->Fill(tp.eta());
            }
            if (tp.useForAlgEff()) {  // Check TP is good for algorithmic efficiency measurement.
              hisFitTPinvptForAlgEff_[fitName]->Fill(1. / tp.pt());
              hisFitTPetaForAlgEff_[fitName]->Fill(tp.eta());
              hisFitTPphiForAlgEff_[fitName]->Fill(tp.phi0());
              // Plot also production point of all good reconstructed TP.
              hisFitTPd0ForAlgEff_[fitName]->Fill(std::abs(tp.d0()));
              hisFitTPz0ForAlgEff_[fitName]->Fill(std::abs(tp.z0()));
              // Also plot efficiency to perfectly reconstruct the track (no fake hits)
              if (tpPerfRecoedMap[&tp]) {
                hisPerfFitTPinvptForAlgEff_[fitName]->Fill(1. / tp.pt());
                hisPerfFitTPetaForAlgEff_[fitName]->Fill(tp.eta());
              }
            }
          }
        }
      }
    }
  }

void Histos::makeEfficiencyPlot ( TFileDirectory &  inputDir, TEfficiency *  outputEfficiency, TH1F *  pass, TH1F *  all, TString  name, TString  title  )

protectedvirtual

Definition at line 1464 of file Histos.cc.

References python.cmstools::all(), and TFileDirectory::make().

Referenced by plotTrackEffAfterFit(), and plotTrackEfficiency().

                                                                                                                   {
    outputEfficiency = inputDir.make<TEfficiency>(*pass, *all);
    outputEfficiency->SetName(name);
    outputEfficiency->SetTitle(title);
  }

virtual void tmtt::Histos::plotHybridDupRemovalEfficiency ( )

inlineprotectedvirtual

Definition at line 94 of file Histos.h.

Referenced by endJobAnalysis().

{};

TFileDirectory Histos::plotTrackEffAfterFit ( const std::string &  fitName )

protectedvirtual

Definition at line 1363 of file Histos.cc.

References fs_, hisFitTPd0ForAlgEff_, hisFitTPd0ForEff_, hisFitTPetaForAlgEff_, hisFitTPetaForEff_, hisFitTPinvptForAlgEff_, hisFitTPinvptForEff_, hisFitTPphiForAlgEff_, hisFitTPphiForEff_, hisFitTPz0ForAlgEff_, hisFitTPz0ForEff_, hisPerfFitTPetaForAlgEff_, hisPerfFitTPetaForEff_, hisPerfFitTPinvptForAlgEff_, hisPerfFitTPinvptForEff_, hisTPd0ForAlgEff_, hisTPd0ForEff_, hisTPetaForAlgEff_, hisTPetaForEff_, hisTPinvptForAlgEff_, hisTPinvptForEff_, hisTPphiForAlgEff_, hisTPphiForEff_, hisTPz0ForAlgEff_, hisTPz0ForEff_, makeEfficiencyPlot(), TFileService::mkdir(), alignCSCRings::s, teffAlgEffFitVsD0_, teffAlgEffFitVsEta_, teffAlgEffFitVsInvPt_, teffAlgEffFitVsPhi_, teffAlgEffFitVsZ0_, teffEffFitVsD0_, teffEffFitVsEta_, teffEffFitVsInvPt_, teffEffFitVsPhi_, teffEffFitVsZ0_, teffPerfAlgEffFitVsEta_, teffPerfAlgEffFitVsInvPt_, teffPerfEffFitVsEta_, and teffPerfEffFitVsInvPt_.

Referenced by endJobAnalysis().

                                                                   {
    // Define lambda function to facilitate adding "fitName" to directory & histogram names.
    auto addn = [fitName](const string& s) { return TString::Format("%s_%s", s.c_str(), fitName.c_str()); };

    TFileDirectory inputDir = fs_->mkdir(addn("Effi").Data());
    // Plot tracking efficiency
    makeEfficiencyPlot(inputDir,
                       teffEffFitVsInvPt_[fitName],
                       hisFitTPinvptForEff_[fitName],
                       hisTPinvptForEff_,
                       addn("EffFitVsInvPt"),
                       "; 1/Pt; Tracking efficiency");
    makeEfficiencyPlot(inputDir,
                       teffEffFitVsEta_[fitName],
                       hisFitTPetaForEff_[fitName],
                       hisTPetaForEff_,
                       addn("EffFitVsEta"),
                       "; #eta; Tracking efficiency");
    makeEfficiencyPlot(inputDir,
                       teffEffFitVsPhi_[fitName],
                       hisFitTPphiForEff_[fitName],
                       hisTPphiForEff_,
                       addn("EffFitVsPhi"),
                       "; #phi; Tracking efficiency");

    makeEfficiencyPlot(inputDir,
                       teffEffFitVsD0_[fitName],
                       hisFitTPd0ForEff_[fitName],
                       hisTPd0ForEff_,
                       addn("EffFitVsD0"),
                       "; d0 (cm); Tracking efficiency");
    makeEfficiencyPlot(inputDir,
                       teffEffFitVsZ0_[fitName],
                       hisFitTPz0ForEff_[fitName],
                       hisTPz0ForEff_,
                       addn("EffFitVsZ0"),
                       "; z0 (cm); Tracking efficiency");

    // Also plot efficiency to reconstruct track perfectly.
    makeEfficiencyPlot(inputDir,
                       teffPerfEffFitVsInvPt_[fitName],
                       hisPerfFitTPinvptForEff_[fitName],
                       hisTPinvptForEff_,
                       addn("PerfEffFitVsInvPt"),
                       "; 1/Pt; Tracking perfect efficiency");
    makeEfficiencyPlot(inputDir,
                       teffPerfEffFitVsEta_[fitName],
                       hisPerfFitTPetaForEff_[fitName],
                       hisTPetaForEff_,
                       addn("PerfEffFitVsEta"),
                       "; #eta; Tracking perfect efficiency");

    // Plot algorithmic tracking efficiency
    makeEfficiencyPlot(inputDir,
                       teffAlgEffFitVsInvPt_[fitName],
                       hisFitTPinvptForAlgEff_[fitName],
                       hisTPinvptForAlgEff_,
                       addn("AlgEffFitVsInvPt"),
                       "; 1/Pt; Tracking efficiency");
    makeEfficiencyPlot(inputDir,
                       teffAlgEffFitVsEta_[fitName],
                       hisFitTPetaForAlgEff_[fitName],
                       hisTPetaForAlgEff_,
                       addn("AlgEffFitVsEta"),
                       "; #eta; Tracking efficiency");
    makeEfficiencyPlot(inputDir,
                       teffAlgEffFitVsPhi_[fitName],
                       hisFitTPphiForAlgEff_[fitName],
                       hisTPphiForAlgEff_,
                       addn("AlgEffFitVsPhi"),
                       "; #phi; Tracking efficiency");

    makeEfficiencyPlot(inputDir,
                       teffAlgEffFitVsD0_[fitName],
                       hisFitTPd0ForAlgEff_[fitName],
                       hisTPd0ForAlgEff_,
                       addn("AlgEffFitVsD0"),
                       "; d0 (cm); Tracking efficiency");
    makeEfficiencyPlot(inputDir,
                       teffAlgEffFitVsZ0_[fitName],
                       hisFitTPz0ForAlgEff_[fitName],
                       hisTPz0ForAlgEff_,
                       addn("AlgEffFitVsZ0"),
                       "; z0 (cm); Tracking efficiency");

    // Also plot algorithmic efficiency to reconstruct track perfectly.
    makeEfficiencyPlot(inputDir,
                       teffPerfAlgEffFitVsInvPt_[fitName],
                       hisPerfFitTPinvptForAlgEff_[fitName],
                       hisTPinvptForAlgEff_,
                       addn("PerfAlgEffFitVsInvPt"),
                       "; 1/Pt; Tracking perfect efficiency");
    makeEfficiencyPlot(inputDir,
                       teffPerfAlgEffFitVsEta_[fitName],
                       hisPerfFitTPetaForAlgEff_[fitName],
                       hisTPetaForAlgEff_,
                       addn("Perf AlgEffFitVsEta"),
                       "; #eta; Tracking perfect efficiency");
    return inputDir;
  }

TFileDirectory Histos::plotTrackEfficiency ( const std::string &  tName )

protectedvirtual

Definition at line 1258 of file Histos.cc.

References fs_, hisPerfRecoTPetaForAlgEff_, hisPerfRecoTPetaForEff_, hisPerfRecoTPinvptForAlgEff_, hisPerfRecoTPinvptForEff_, hisRecoTPd0ForAlgEff_, hisRecoTPd0ForEff_, hisRecoTPetaForAlgEff_, hisRecoTPetaForEff_, hisRecoTPinvptForAlgEff_, hisRecoTPinvptForEff_, hisRecoTPphiForAlgEff_, hisRecoTPphiForEff_, hisRecoTPz0ForAlgEff_, hisRecoTPz0ForEff_, hisTPd0ForAlgEff_, hisTPd0ForEff_, hisTPetaForAlgEff_, hisTPetaForEff_, hisTPinvptForAlgEff_, hisTPinvptForEff_, hisTPphiForAlgEff_, hisTPphiForEff_, hisTPz0ForAlgEff_, hisTPz0ForEff_, makeEfficiencyPlot(), TFileService::mkdir(), alignCSCRings::s, teffAlgEffVsD0_, teffAlgEffVsEta_, teffAlgEffVsInvPt_, teffAlgEffVsPhi_, teffAlgEffVsZ0_, teffEffVsD0_, teffEffVsEta_, teffEffVsInvPt_, teffEffVsPhi_, teffEffVsZ0_, teffPerfAlgEffVsEta_, teffPerfAlgEffVsInvPt_, teffPerfEffVsEta_, and teffPerfEffVsInvPt_.

Referenced by endJobAnalysis().

                                                                {
    // Define lambda function to facilitate adding "tName" to directory & histogram names.
    auto addn = [tName](const string& s) { return TString::Format("%s_%s", s.c_str(), tName.c_str()); };

    TFileDirectory inputDir = fs_->mkdir(addn("Effi").Data());
    // Plot tracking efficiency
    makeEfficiencyPlot(inputDir,
                       teffEffVsInvPt_[tName],
                       hisRecoTPinvptForEff_[tName],
                       hisTPinvptForEff_,
                       addn("EffVsInvPt"),
                       "; 1/Pt; Tracking efficiency");
    makeEfficiencyPlot(inputDir,
                       teffEffVsEta_[tName],
                       hisRecoTPetaForEff_[tName],
                       hisTPetaForEff_,
                       addn("EffVsEta"),
                       "; #eta; Tracking efficiency");

    makeEfficiencyPlot(inputDir,
                       teffEffVsPhi_[tName],
                       hisRecoTPphiForEff_[tName],
                       hisTPphiForEff_,
                       addn("EffVsPhi"),
                       "; #phi; Tracking efficiency");

    makeEfficiencyPlot(inputDir,
                       teffEffVsD0_[tName],
                       hisRecoTPd0ForEff_[tName],
                       hisTPd0ForEff_,
                       addn("EffVsD0"),
                       "; d0 (cm); Tracking efficiency");
    makeEfficiencyPlot(inputDir,
                       teffEffVsZ0_[tName],
                       hisRecoTPz0ForEff_[tName],
                       hisTPz0ForEff_,
                       addn("EffVsZ0"),
                       "; z0 (cm); Tracking efficiency");

    // Also plot efficiency to reconstruct track perfectly.
    makeEfficiencyPlot(inputDir,
                       teffPerfEffVsInvPt_[tName],
                       hisPerfRecoTPinvptForEff_[tName],
                       hisTPinvptForEff_,
                       addn("PerfEffVsInvPt"),
                       "; 1/Pt; Tracking perfect efficiency");
    makeEfficiencyPlot(inputDir,
                       teffPerfEffVsEta_[tName],
                       hisPerfRecoTPetaForEff_[tName],
                       hisTPetaForEff_,
                       addn("PerfEffVsEta"),
                       "; #eta; Tracking perfect efficiency");

    // Plot algorithmic tracking efficiency
    makeEfficiencyPlot(inputDir,
                       teffAlgEffVsInvPt_[tName],
                       hisRecoTPinvptForAlgEff_[tName],
                       hisTPinvptForAlgEff_,
                       addn("AlgEffVsInvPt"),
                       "; 1/Pt; Tracking efficiency");
    makeEfficiencyPlot(inputDir,
                       teffAlgEffVsEta_[tName],
                       hisRecoTPetaForAlgEff_[tName],
                       hisTPetaForAlgEff_,
                       addn("AlgEffVsEta"),
                       "; #eta; Tracking efficiency");
    makeEfficiencyPlot(inputDir,
                       teffAlgEffVsPhi_[tName],
                       hisRecoTPphiForAlgEff_[tName],
                       hisTPphiForAlgEff_,
                       addn("AlgEffVsPhi"),
                       "; #phi; Tracking efficiency");

    makeEfficiencyPlot(inputDir,
                       teffAlgEffVsD0_[tName],
                       hisRecoTPd0ForAlgEff_[tName],
                       hisTPd0ForAlgEff_,
                       addn("AlgEffVsD0"),
                       "; d0 (cm); Tracking efficiency");
    makeEfficiencyPlot(inputDir,
                       teffAlgEffVsZ0_[tName],
                       hisRecoTPz0ForAlgEff_[tName],
                       hisTPz0ForAlgEff_,
                       addn("AlgEffVsZ0"),
                       "; z0 (cm); Tracking efficiency");

    // Also plot algorithmic efficiency to reconstruct track perfectly.
    makeEfficiencyPlot(inputDir,
                       teffPerfAlgEffVsInvPt_[tName],
                       hisPerfRecoTPinvptForAlgEff_[tName],
                       hisTPinvptForAlgEff_,
                       addn("PerfAlgEffVsInvPt"),
                       "; 1/Pt; Tracking perfect efficiency");
    makeEfficiencyPlot(inputDir,
                       teffPerfAlgEffVsEta_[tName],
                       hisPerfRecoTPetaForAlgEff_[tName],
                       hisTPetaForAlgEff_,
                       addn("PerfAlgEffVsEta"),
                       "; #eta; Tracking perfect efficiency");

    return inputDir;
  }

virtual void tmtt::Histos::plotTrackletSeedEfficiency ( )

inlineprotectedvirtual

Definition at line 92 of file Histos.h.

Referenced by endJobAnalysis().

{};

void Histos::printFitTrackPerformance ( const std::string &  fitName )

protectedvirtual

Definition at line 1519 of file Histos.cc.

References submitPVResolutionJobs::count, alignCSCRings::e, hisFitTPinvptForAlgEff_, hisPerfFitTPinvptForAlgEff_, hisTPinvptForAlgEff_, profNumFitTracks_, profStubsOnFitTracks_, tmtt::Settings::rzFilterName(), settings_, mathSSE::sqrt(), and useRZfilter_.

Referenced by endJobAnalysis().

                                                             {
    float numFitTracks = profNumFitTracks_[fitName]->GetBinContent(1);   // No. of track cands
    float numFitTracksErr = profNumFitTracks_[fitName]->GetBinError(1);  // No. of track cands uncertainty
    float numMatchedFitTracksIncDups =
        profNumFitTracks_[fitName]->GetBinContent(2);  // Ditto, counting only those matched to TP
    float numMatchedFitTracksExcDups = profNumFitTracks_[fitName]->GetBinContent(6);  // Ditto, but excluding duplicates
    float numFakeFitTracks = numFitTracks - numMatchedFitTracksIncDups;
    float numExtraDupFitTracks = numMatchedFitTracksIncDups - numMatchedFitTracksExcDups;
    float fracFakeFit = numFakeFitTracks / (numFitTracks + 1.0e-6);
    float fracDupFit = numExtraDupFitTracks / (numFitTracks + 1.0e-6);

    float numStubsOnFitTracks = profStubsOnFitTracks_[fitName]->GetBinContent(1);
    float meanStubsPerFitTrack =
        numStubsOnFitTracks / (numFitTracks + 1.0e-6);  //protection against demoninator equals zero.
    unsigned int numFitTPforAlg = hisFitTPinvptForAlgEff_[fitName]->GetEntries();
    // Histograms of input truth particles (e.g. hisTPinvptForAlgEff_), used for denominator of efficiencies, are identical,
    // irrespective of whether made after HT or after r-z track filter, so always use the former.
    unsigned int numTPforAlg = hisTPinvptForAlgEff_->GetEntries();
    unsigned int numPerfFitTPforAlg = hisPerfFitTPinvptForAlgEff_[fitName]->GetEntries();
    float fitEff = float(numFitTPforAlg) / (numTPforAlg + 1.0e-6);  //protection against demoninator equals zero.
    float fitEffErr = sqrt(fitEff * (1 - fitEff) / (numTPforAlg + 1.0e-6));  // uncertainty
    float fitPerfEff =
        float(numPerfFitTPforAlg) / (numTPforAlg + 1.0e-6);  //protection against demoninator equals zero.
    float fitPerfEffErr = sqrt(fitPerfEff * (1 - fitPerfEff) / (numTPforAlg + 1.0e-6));  // uncertainty

    // Does this fitter require r-z track filter to be run before it?
    bool useRZfilt = (std::count(useRZfilter_.begin(), useRZfilter_.end(), fitName) > 0);

    PrintL1trk() << "=========================================================================";
    PrintL1trk() << "                    TRACK FIT SUMMARY FOR: " << fitName;
    PrintL1trk() << "Number of fitted track candidates found per event = " << numFitTracks << " +- " << numFitTracksErr;
    PrintL1trk() << "                     with mean stubs/track = " << meanStubsPerFitTrack;
    PrintL1trk() << "Fraction of fitted tracks that are fake = " << fracFakeFit;
    PrintL1trk() << "Fraction of fitted tracks that are genuine, but extra duplicates = " << fracDupFit;
    PrintL1trk() << "Algorithmic fitting efficiency = " << numFitTPforAlg << "/" << numTPforAlg << " = " << fitEff
                 << " +- " << fitEffErr;
    PrintL1trk() << "Perfect algorithmic fitting efficiency = " << numPerfFitTPforAlg << "/" << numTPforAlg << " = "
                 << fitPerfEff << " +- " << fitPerfEffErr << " (no incorrect hits)";
    if (useRZfilt)
      PrintL1trk() << "(The above fitter used the '" << settings_->rzFilterName() << "' r-z track filter.)";
  }

virtual void tmtt::Histos::printHybridDupRemovalPerformance ( )

inlineprotectedvirtual

Definition at line 110 of file Histos.h.

Referenced by endJobAnalysis().

{};

virtual void tmtt::Histos::printTrackletSeedFindingPerformance ( )

inlineprotectedvirtual

Definition at line 107 of file Histos.h.

Referenced by endJobAnalysis().

{};

void Histos::printTrackPerformance ( const std::string &  tName )

protectedvirtual

Definition at line 1473 of file Histos.cc.

References alignCSCRings::e, hisPerfRecoTPinvptForAlgEff_, hisRecoTPinvptForAlgEff_, hisTPinvptForAlgEff_, profNumTrackCands_, profStubsOnTracks_, and mathSSE::sqrt().

Referenced by endJobAnalysis().

                                                        {
    float numTrackCands = profNumTrackCands_[tName]->GetBinContent(1);   // No. of track cands
    float numTrackCandsErr = profNumTrackCands_[tName]->GetBinError(1);  // No. of track cands uncertainty
    float numMatchedTrackCandsIncDups =
        profNumTrackCands_[tName]->GetBinContent(2);  // Ditto, counting only those matched to TP
    float numMatchedTrackCandsExcDups = profNumTrackCands_[tName]->GetBinContent(6);  // Ditto, but excluding duplicates
    float numFakeTracks = numTrackCands - numMatchedTrackCandsIncDups;
    float numExtraDupTracks = numMatchedTrackCandsIncDups - numMatchedTrackCandsExcDups;
    float fracFake = numFakeTracks / (numTrackCands + 1.0e-6);
    float fracDup = numExtraDupTracks / (numTrackCands + 1.0e-6);

    float numStubsOnTracks = profStubsOnTracks_[tName]->GetBinContent(1);
    float meanStubsPerTrack =
        numStubsOnTracks / (numTrackCands + 1.0e-6);  //protection against demoninator equals zero.
    unsigned int numRecoTPforAlg = hisRecoTPinvptForAlgEff_[tName]->GetEntries();
    // Histograms of input truth particles (e.g. hisTPinvptForAlgEff_), used for denominator of efficiencies, are identical,
    // irrespective of whether made after HT or after r-z track filter, so always use the former.
    unsigned int numTPforAlg = hisTPinvptForAlgEff_->GetEntries();
    unsigned int numPerfRecoTPforAlg = hisPerfRecoTPinvptForAlgEff_[tName]->GetEntries();
    float algEff = float(numRecoTPforAlg) / (numTPforAlg + 1.0e-6);  //protection against demoninator equals zero.
    float algEffErr = sqrt(algEff * (1 - algEff) / (numTPforAlg + 1.0e-6));  // uncertainty
    float algPerfEff =
        float(numPerfRecoTPforAlg) / (numTPforAlg + 1.0e-6);  //protection against demoninator equals zero.
    float algPerfEffErr = sqrt(algPerfEff * (1 - algPerfEff) / (numTPforAlg + 1.0e-6));  // uncertainty

    PrintL1trk() << "=========================================================================";
    if (tName == "HT") {
      PrintL1trk() << "               TRACK-FINDING SUMMARY AFTER HOUGH TRANSFORM             ";
    } else if (tName == "RZ") {
      PrintL1trk() << "               TRACK-FINDING SUMMARY AFTER R-Z TRACK FILTER            ";
    } else if (tName == "TRACKLET") {
      PrintL1trk() << "               TRACK-FINDING SUMMARY AFTER TRACKLET PATTERN RECO       ";
    }
    PrintL1trk() << "Number of track candidates found per event = " << numTrackCands << " +- " << numTrackCandsErr;
    PrintL1trk() << "                     with mean stubs/track = " << meanStubsPerTrack;
    PrintL1trk() << "Fraction of track cands that are fake = " << fracFake;
    PrintL1trk() << "Fraction of track cands that are genuine, but extra duplicates = " << fracDup;

    PrintL1trk() << std::fixed << std::setprecision(4) << "Algorithmic tracking efficiency = " << numRecoTPforAlg << "/"
                 << numTPforAlg << " = " << algEff << " +- " << algEffErr;
    PrintL1trk() << "Perfect algorithmic tracking efficiency = " << numPerfRecoTPforAlg << "/" << numTPforAlg << " = "
                 << algPerfEff << " +- " << algPerfEffErr << " (no incorrect hits)";
  }

void Histos::trackerGeometryAnalysis ( const std::list< TrackerModule > &  listTrackerModule )

virtual

Definition at line 1705 of file Histos.cc.

References funct::abs(), tmtt::Settings::bApprox_gradient(), tmtt::Settings::bApprox_intercept(), bApproxMistake_, enabled(), fs_, graphBVsZoverR_, mps_fire::i, TFileDirectory::make(), TFileService::mkdir(), HLT_FULL_cff::nSigma, AlCaHLTBitMon_ParallelJobs::p, settings_, and tmtt::Settings::useApproxB().

Referenced by tmtt::TMTrackProducer::globalEndJob().

                                                                                   {
    // Don't bother producing summary if user didn't request histograms via TFileService in their cfg.
    if (not this->enabled())
      return;

    map<float, float> dataForGraph;
    for (const TrackerModule& trackerModule : listTrackerModule) {
      if (trackerModule.tiltedBarrel()) {
        float paramB = trackerModule.paramB();
        float zOverR = std::abs(trackerModule.minZ()) / trackerModule.minR();
        dataForGraph[paramB] = zOverR;  // Only store each unique paramB once, to get better histo weights.
      }
    }

    const int nEntries = dataForGraph.size();
    vector<float> vecParamB(nEntries);
    vector<float> vecZoverR(nEntries);
    unsigned int i = 0;
    for (const auto& p : dataForGraph) {
      vecParamB[i] = p.first;
      vecZoverR[i] = p.second;
      i++;
    }

    PrintL1trk() << "=========================================================================";
    PrintL1trk() << "--- Fit to cfg params for FPGA-friendly approximation to B parameter in GP & KF ---";
    PrintL1trk() << "--- (used to allowed for tilted barrel modules)                                 ---";

    TFileDirectory inputDir = fs_->mkdir("InputDataB");
    graphBVsZoverR_ = inputDir.make<TGraph>(nEntries, &vecZoverR[0], &vecParamB[0]);
    graphBVsZoverR_->SetNameTitle("B vs module Z/R", "; Module Z/R; B");
    graphBVsZoverR_->Fit("pol1", "q");
    TF1* fittedFunction = graphBVsZoverR_->GetFunction("pol1");
    double gradient = fittedFunction->GetParameter(1);
    double intercept = fittedFunction->GetParameter(0);
    double gradient_err = fittedFunction->GetParError(1);
    double intercept_err = fittedFunction->GetParError(0);
    PrintL1trk() << "         BApprox_gradient (fitted)  = " << gradient << " +- " << gradient_err;
    PrintL1trk() << "         BApprox_intercept (fitted) = " << intercept << " +- " << intercept_err;
    PrintL1trk() << "=========================================================================";

    // Check fitted params consistent with those assumed in cfg file.
    if (settings_->useApproxB()) {
      double gradientDiff = std::abs(gradient - settings_->bApprox_gradient());
      double interceptDiff = std::abs(intercept - settings_->bApprox_intercept());
      constexpr unsigned int nSigma = 5;
      if (gradientDiff > nSigma * gradient_err ||
          interceptDiff > nSigma * intercept_err) {  // Uncertainty independent of number of events
        PrintL1trk() << "\n WARNING: fitted parameters inconsistent with those specified in cfg file:";
        PrintL1trk() << "         BApprox_gradient  (cfg) = " << settings_->bApprox_gradient();
        PrintL1trk() << "         BApprox_intercept (cfg) = " << settings_->bApprox_intercept();
        bApproxMistake_ = true;  // Note that problem has occurred.
      }
    }
  }

Member Data Documentation

bool tmtt::Histos::bApproxMistake_

protected

Definition at line 357 of file Histos.h.

Referenced by endJobAnalysis(), and trackerGeometryAnalysis().

float tmtt::Histos::chosenRofZ_

protected

Definition at line 123 of file Histos.h.

Referenced by Histos().

edm::Service<TFileService> tmtt::Histos::fs_

protected

Definition at line 110 of file Histos.h.

Referenced by available(), bookEtaPhiSectors(), bookInputData(), bookRphiHT(), bookRZfilters(), bookTrackCands(), bookTrackFitting(), plotTrackEffAfterFit(), plotTrackEfficiency(), and trackerGeometryAnalysis().

unsigned int tmtt::Histos::genMinStubLayers_

protected

Definition at line 117 of file Histos.h.

Referenced by Histos().

TGraph* tmtt::Histos::graphBVsZoverR_

protected

Definition at line 145 of file Histos.h.

Referenced by trackerGeometryAnalysis().

TH1F* tmtt::Histos::hisBendResStub_

protected

Definition at line 142 of file Histos.h.

Referenced by bookInputData(), and fillInputData().

TH1F* tmtt::Histos::hisBendStub_

protected

Definition at line 141 of file Histos.h.

Referenced by bookInputData(), and fillInputData().

std::map<std::string, TProfile*> tmtt::Histos::hisD0ResVsTrueEta_

protected

Definition at line 253 of file Histos.h.

Referenced by bookTrackFitting(), and fillTrackFitting().

std::map<std::string, TProfile*> tmtt::Histos::hisD0ResVsTrueInvPt_

protected

Definition at line 259 of file Histos.h.

Referenced by bookTrackFitting(), and fillTrackFitting().

std::map<std::string, TH1F*> tmtt::Histos::hisEta_

protected

Definition at line 171 of file Histos.h.

Referenced by bookTrackCands().

std::map<std::string, TH1F*> tmtt::Histos::hisEtaRes_

protected

Definition at line 177 of file Histos.h.

Referenced by bookTrackCands().

std::map<std::string, TProfile*> tmtt::Histos::hisEtaResVsTrueEta_

protected

Definition at line 251 of file Histos.h.

Referenced by bookTrackFitting(), and fillTrackFitting().

std::map<std::string, TProfile*> tmtt::Histos::hisEtaResVsTrueInvPt_

protected

Definition at line 257 of file Histos.h.

Referenced by bookTrackFitting(), and fillTrackFitting().

std::map<std::string, TH1F*> tmtt::Histos::hisFitChi2DofRphiMatched_

protected

Definition at line 241 of file Histos.h.

Referenced by bookTrackFitting(), and fillTrackFitting().

std::map<std::string, TH1F*> tmtt::Histos::hisFitChi2DofRphiUnmatched_

protected

Definition at line 245 of file Histos.h.

Referenced by bookTrackFitting(), and fillTrackFitting().

std::map<std::string, TH1F*> tmtt::Histos::hisFitChi2DofRzMatched_

protected

Definition at line 242 of file Histos.h.

Referenced by bookTrackFitting(), and fillTrackFitting().

std::map<std::string, TH1F*> tmtt::Histos::hisFitChi2DofRzUnmatched_

protected

Definition at line 246 of file Histos.h.

Referenced by bookTrackFitting(), and fillTrackFitting().

std::map<std::string, TH1F*> tmtt::Histos::hisFitD0Matched_

protected

Definition at line 223 of file Histos.h.

Referenced by bookTrackFitting(), and fillTrackFitting().

std::map<std::string, TH1F*> tmtt::Histos::hisFitD0Unmatched_

protected

Definition at line 229 of file Histos.h.

Referenced by bookTrackFitting(), and fillTrackFitting().

std::map<std::string, TH1F*> tmtt::Histos::hisFitEtaMatched_

protected

Definition at line 225 of file Histos.h.

Referenced by bookTrackFitting(), and fillTrackFitting().

std::map<std::string, TH1F*> tmtt::Histos::hisFitEtaUnmatched_

protected

Definition at line 231 of file Histos.h.

Referenced by bookTrackFitting(), and fillTrackFitting().

std::map<std::string, TH1F*> tmtt::Histos::hisFitPhi0Matched_

protected

Definition at line 222 of file Histos.h.

Referenced by bookTrackFitting(), and fillTrackFitting().

std::map<std::string, TH1F*> tmtt::Histos::hisFitPhi0Unmatched_

protected

Definition at line 228 of file Histos.h.

Referenced by bookTrackFitting(), and fillTrackFitting().

std::map<std::string, TH1F*> tmtt::Histos::hisFitQinvPtMatched_

protected

Definition at line 221 of file Histos.h.

Referenced by bookTrackFitting(), and fillTrackFitting().

std::map<std::string, TH1F*> tmtt::Histos::hisFitQinvPtUnmatched_

protected

Definition at line 227 of file Histos.h.

Referenced by bookTrackFitting(), and fillTrackFitting().

std::map<std::string, TH1F*> tmtt::Histos::hisFitTPd0ForAlgEff_

protected

Definition at line 275 of file Histos.h.

Referenced by bookTrackFitting(), fillTrackFitting(), and plotTrackEffAfterFit().

std::map<std::string, TH1F*> tmtt::Histos::hisFitTPd0ForEff_

protected

Definition at line 268 of file Histos.h.

Referenced by bookTrackFitting(), fillTrackFitting(), and plotTrackEffAfterFit().

std::map<std::string, TH1F*> tmtt::Histos::hisFitTPetaForAlgEff_

protected

Definition at line 273 of file Histos.h.

Referenced by bookTrackFitting(), fillTrackFitting(), and plotTrackEffAfterFit().

std::map<std::string, TH1F*> tmtt::Histos::hisFitTPetaForEff_

protected

Definition at line 266 of file Histos.h.

Referenced by bookTrackFitting(), fillTrackFitting(), and plotTrackEffAfterFit().

std::map<std::string, TH1F*> tmtt::Histos::hisFitTPinvptForAlgEff_

protected

Definition at line 272 of file Histos.h.

Referenced by bookTrackFitting(), fillTrackFitting(), plotTrackEffAfterFit(), and printFitTrackPerformance().

std::map<std::string, TH1F*> tmtt::Histos::hisFitTPinvptForEff_

protected

Definition at line 265 of file Histos.h.

Referenced by bookTrackFitting(), fillTrackFitting(), and plotTrackEffAfterFit().

std::map<std::string, TH1F*> tmtt::Histos::hisFitTPphiForAlgEff_

protected

Definition at line 274 of file Histos.h.

Referenced by bookTrackFitting(), fillTrackFitting(), and plotTrackEffAfterFit().

std::map<std::string, TH1F*> tmtt::Histos::hisFitTPphiForEff_

protected

Definition at line 267 of file Histos.h.

Referenced by bookTrackFitting(), fillTrackFitting(), and plotTrackEffAfterFit().

std::map<std::string, TH1F*> tmtt::Histos::hisFitTPz0ForAlgEff_

protected

Definition at line 276 of file Histos.h.

Referenced by bookTrackFitting(), fillTrackFitting(), and plotTrackEffAfterFit().

std::map<std::string, TH1F*> tmtt::Histos::hisFitTPz0ForEff_

protected

Definition at line 269 of file Histos.h.

Referenced by bookTrackFitting(), fillTrackFitting(), and plotTrackEffAfterFit().

std::map<std::string, TH1F*> tmtt::Histos::hisFitZ0Matched_

protected

Definition at line 224 of file Histos.h.

Referenced by bookTrackFitting(), and fillTrackFitting().

std::map<std::string, TH1F*> tmtt::Histos::hisFitZ0Unmatched_

protected

Definition at line 230 of file Histos.h.

Referenced by bookTrackFitting(), and fillTrackFitting().

std::map<std::string, TH1F*> tmtt::Histos::hisFracMatchStubsOnTracks_

protected

Definition at line 164 of file Histos.h.

Referenced by bookTrackCands().

std::map<std::string, TH1F*> tmtt::Histos::hisKalmanChi2DofSkipLay0Matched_

protected

Definition at line 234 of file Histos.h.

Referenced by bookTrackFitting(), and fillTrackFitting().

std::map<std::string, TH1F*> tmtt::Histos::hisKalmanChi2DofSkipLay0Unmatched_

protected

Definition at line 237 of file Histos.h.

Referenced by bookTrackFitting(), and fillTrackFitting().

std::map<std::string, TH1F*> tmtt::Histos::hisKalmanChi2DofSkipLay1Matched_

protected

Definition at line 235 of file Histos.h.

Referenced by bookTrackFitting(), and fillTrackFitting().

std::map<std::string, TH1F*> tmtt::Histos::hisKalmanChi2DofSkipLay1Unmatched_

protected

Definition at line 238 of file Histos.h.

Referenced by bookTrackFitting(), and fillTrackFitting().

std::map<std::string, TH1F*> tmtt::Histos::hisKalmanChi2DofSkipLay2Matched_

protected

Definition at line 236 of file Histos.h.

Referenced by bookTrackFitting(), and fillTrackFitting().

std::map<std::string, TH1F*> tmtt::Histos::hisKalmanChi2DofSkipLay2Unmatched_

protected

Definition at line 239 of file Histos.h.

Referenced by bookTrackFitting(), and fillTrackFitting().

std::map<std::string, TH1F*> tmtt::Histos::hisKalmanNumUpdateCalls_

protected

Definition at line 233 of file Histos.h.

Referenced by bookTrackFitting(), and fillTrackFitting().

std::map<std::string, TH1F*> tmtt::Histos::hisLayersPerTrack_

protected

Definition at line 160 of file Histos.h.

Referenced by bookTrackCands().

TH1F* tmtt::Histos::hisNumEtaSecsPerStub_

protected

Definition at line 148 of file Histos.h.

Referenced by bookEtaPhiSectors(), and fillEtaPhiSectors().

std::map<std::string, TH1F*> tmtt::Histos::hisNumFitTrks_

protected

Definition at line 215 of file Histos.h.

Referenced by bookTrackFitting(), and fillTrackFitting().

std::map<std::string, TH1F*> tmtt::Histos::hisNumFitTrksPerNon_

protected

Definition at line 216 of file Histos.h.

Referenced by bookTrackFitting(), and fillTrackFitting().

std::map<std::string, TH1F*> tmtt::Histos::hisNumFitTrksPerSect_

protected

Definition at line 217 of file Histos.h.

TH1F* tmtt::Histos::hisNumLayersPerTP_

protected

Definition at line 135 of file Histos.h.

Referenced by bookInputData(), and fillInputData().

TH1F* tmtt::Histos::hisNumPhiSecsPerStub_

protected

Definition at line 149 of file Histos.h.

Referenced by bookEtaPhiSectors(), and fillEtaPhiSectors().

TH1F* tmtt::Histos::hisNumPSLayersPerTP_

protected

Definition at line 136 of file Histos.h.

Referenced by bookInputData(), and fillInputData().

std::map<std::string, TH1F*> tmtt::Histos::hisNumStubsPerLink_

protected

Definition at line 162 of file Histos.h.

Referenced by bookTrackCands().

TH1F* tmtt::Histos::hisNumStubsPerSec_

protected

Definition at line 150 of file Histos.h.

Referenced by bookEtaPhiSectors(), and fillEtaPhiSectors().

std::map<std::string, TH1F*> tmtt::Histos::hisNumTracksVsQoverPt_

protected

Definition at line 155 of file Histos.h.

Referenced by bookTrackCands().

std::map<std::string, TH1F*> tmtt::Histos::hisNumTrksPerNon_

protected

Definition at line 156 of file Histos.h.

Referenced by bookTrackCands().

std::map<std::string, TH1F*> tmtt::Histos::hisPerfFitTPetaForAlgEff_

protected

Definition at line 278 of file Histos.h.

Referenced by bookTrackFitting(), fillTrackFitting(), and plotTrackEffAfterFit().

std::map<std::string, TH1F*> tmtt::Histos::hisPerfFitTPetaForEff_

protected

Definition at line 271 of file Histos.h.

Referenced by bookTrackFitting(), fillTrackFitting(), and plotTrackEffAfterFit().

std::map<std::string, TH1F*> tmtt::Histos::hisPerfFitTPinvptForAlgEff_

protected

Definition at line 277 of file Histos.h.

Referenced by bookTrackFitting(), fillTrackFitting(), plotTrackEffAfterFit(), and printFitTrackPerformance().

std::map<std::string, TH1F*> tmtt::Histos::hisPerfFitTPinvptForEff_

protected

Definition at line 270 of file Histos.h.

Referenced by bookTrackFitting(), fillTrackFitting(), and plotTrackEffAfterFit().

std::map<std::string, TH1F*> tmtt::Histos::hisPerfRecoTPetaForAlgEff_

protected

Definition at line 210 of file Histos.h.

Referenced by bookTrackCands(), and plotTrackEfficiency().

std::map<std::string, TH1F*> tmtt::Histos::hisPerfRecoTPetaForEff_

protected

Definition at line 201 of file Histos.h.

Referenced by bookTrackCands(), and plotTrackEfficiency().

std::map<std::string, TH1F*> tmtt::Histos::hisPerfRecoTPinvptForAlgEff_

protected

Definition at line 209 of file Histos.h.

Referenced by bookTrackCands(), plotTrackEfficiency(), and printTrackPerformance().

std::map<std::string, TH1F*> tmtt::Histos::hisPerfRecoTPinvptForEff_

protected

Definition at line 200 of file Histos.h.

Referenced by bookTrackCands(), and plotTrackEfficiency().

std::map<std::string, TH1F*> tmtt::Histos::hisPhi0_

protected

Definition at line 170 of file Histos.h.

Referenced by bookTrackCands().

std::map<std::string, TH1F*> tmtt::Histos::hisPhi0Res_

protected

Definition at line 176 of file Histos.h.

Referenced by bookTrackCands().

std::map<std::string, TProfile*> tmtt::Histos::hisPhi0ResVsTrueEta_

protected

Definition at line 250 of file Histos.h.

Referenced by bookTrackFitting(), and fillTrackFitting().

std::map<std::string, TProfile*> tmtt::Histos::hisPhi0ResVsTrueInvPt_

protected

Definition at line 256 of file Histos.h.

Referenced by bookTrackFitting(), and fillTrackFitting().

std::map<std::string, TH1F*> tmtt::Histos::hisQoverPt_

protected

Definition at line 169 of file Histos.h.

Referenced by bookTrackCands().

std::map<std::string, TH1F*> tmtt::Histos::hisQoverPtRes_

protected

Definition at line 175 of file Histos.h.

Referenced by bookTrackCands().

std::map<std::string, TProfile*> tmtt::Histos::hisQoverPtResVsTrueEta_

protected

Definition at line 249 of file Histos.h.

Referenced by bookTrackFitting(), and fillTrackFitting().

std::map<std::string, TProfile*> tmtt::Histos::hisQoverPtResVsTrueInvPt_

protected

Definition at line 255 of file Histos.h.

Referenced by bookTrackFitting(), and fillTrackFitting().

std::map<std::string, TH1F*> tmtt::Histos::hisRecoTPd0ForAlgEff_

protected

Definition at line 206 of file Histos.h.

Referenced by bookTrackCands(), and plotTrackEfficiency().

std::map<std::string, TH1F*> tmtt::Histos::hisRecoTPd0ForEff_

protected

Definition at line 197 of file Histos.h.

Referenced by bookTrackCands(), and plotTrackEfficiency().

std::map<std::string, TH1F*> tmtt::Histos::hisRecoTPetaForAlgEff_

protected

Definition at line 204 of file Histos.h.

Referenced by bookTrackCands(), and plotTrackEfficiency().

std::map<std::string, TH1F*> tmtt::Histos::hisRecoTPetaForEff_

protected

Definition at line 195 of file Histos.h.

Referenced by bookTrackCands(), and plotTrackEfficiency().

std::map<std::string, TH1F*> tmtt::Histos::hisRecoTPinvptForAlgEff_

protected

Definition at line 203 of file Histos.h.

Referenced by bookTrackCands(), plotTrackEfficiency(), and printTrackPerformance().

std::map<std::string, TH1F*> tmtt::Histos::hisRecoTPinvptForEff_

protected

Definition at line 194 of file Histos.h.

Referenced by bookTrackCands(), endJobAnalysis(), and plotTrackEfficiency().

std::map<std::string, TH1F*> tmtt::Histos::hisRecoTPphiForAlgEff_

protected

Definition at line 205 of file Histos.h.

Referenced by bookTrackCands(), and plotTrackEfficiency().

std::map<std::string, TH1F*> tmtt::Histos::hisRecoTPphiForEff_

protected

Definition at line 196 of file Histos.h.

Referenced by bookTrackCands(), and plotTrackEfficiency().

std::map<std::string, TH1F*> tmtt::Histos::hisRecoTPz0ForAlgEff_

protected

Definition at line 207 of file Histos.h.

Referenced by bookTrackCands(), and plotTrackEfficiency().

std::map<std::string, TH1F*> tmtt::Histos::hisRecoTPz0ForEff_

protected

Definition at line 198 of file Histos.h.

Referenced by bookTrackCands(), and plotTrackEfficiency().

TProfile* tmtt::Histos::hisStubIneffiVsEta_

protected

Definition at line 140 of file Histos.h.

Referenced by bookInputData(), and fillInputData().

TProfile* tmtt::Histos::hisStubIneffiVsInvPt_

protected

Definition at line 139 of file Histos.h.

Referenced by bookInputData(), and fillInputData().

TProfile* tmtt::Histos::hisStubKillFE_

protected

Definition at line 138 of file Histos.h.

Referenced by bookInputData(), and fillInputData().

std::map<std::string, TH1F*> tmtt::Histos::hisStubsOnTracksPerNon_

protected

Definition at line 158 of file Histos.h.

Referenced by bookTrackCands().

std::map<std::string, TH1F*> tmtt::Histos::hisStubsPerFitTrack_

protected

Definition at line 218 of file Histos.h.

Referenced by bookTrackFitting(), and fillTrackFitting().

std::map<std::string, TH1F*> tmtt::Histos::hisStubsPerTrack_

protected

Definition at line 159 of file Histos.h.

Referenced by bookTrackCands().

TH1F* tmtt::Histos::hisStubsVsEta_

protected

Definition at line 132 of file Histos.h.

Referenced by bookInputData(), and fillInputData().

TH1F* tmtt::Histos::hisStubsVsR_

protected

Definition at line 133 of file Histos.h.

Referenced by bookInputData(), and fillInputData().

TH1F* tmtt::Histos::hisTPd0ForAlgEff_

protected

Definition at line 190 of file Histos.h.

Referenced by bookInputData(), fillInputData(), plotTrackEffAfterFit(), and plotTrackEfficiency().

TH1F* tmtt::Histos::hisTPd0ForEff_

protected

Definition at line 184 of file Histos.h.

Referenced by bookInputData(), fillInputData(), plotTrackEffAfterFit(), and plotTrackEfficiency().

TH1F* tmtt::Histos::hisTPetaForAlgEff_

protected

Definition at line 188 of file Histos.h.

Referenced by bookInputData(), fillInputData(), plotTrackEffAfterFit(), and plotTrackEfficiency().

TH1F* tmtt::Histos::hisTPetaForEff_

protected

Definition at line 182 of file Histos.h.

Referenced by bookInputData(), fillInputData(), plotTrackEffAfterFit(), and plotTrackEfficiency().

TH1F* tmtt::Histos::hisTPinvptForAlgEff_

protected

Definition at line 187 of file Histos.h.

Referenced by bookInputData(), fillInputData(), plotTrackEffAfterFit(), plotTrackEfficiency(), printFitTrackPerformance(), and printTrackPerformance().

TH1F* tmtt::Histos::hisTPinvptForEff_

protected

Definition at line 181 of file Histos.h.

Referenced by bookInputData(), fillInputData(), plotTrackEffAfterFit(), and plotTrackEfficiency().

TH1F* tmtt::Histos::hisTPphiForAlgEff_

protected

Definition at line 189 of file Histos.h.

Referenced by bookInputData(), fillInputData(), plotTrackEffAfterFit(), and plotTrackEfficiency().

TH1F* tmtt::Histos::hisTPphiForEff_

protected

Definition at line 183 of file Histos.h.

Referenced by bookInputData(), fillInputData(), plotTrackEffAfterFit(), and plotTrackEfficiency().

TH1F* tmtt::Histos::hisTPz0ForAlgEff_

protected

Definition at line 191 of file Histos.h.

Referenced by bookInputData(), fillInputData(), plotTrackEffAfterFit(), and plotTrackEfficiency().

TH1F* tmtt::Histos::hisTPz0ForEff_

protected

Definition at line 185 of file Histos.h.

Referenced by bookInputData(), fillInputData(), plotTrackEffAfterFit(), and plotTrackEfficiency().

std::map<std::string, TH1F*> tmtt::Histos::hisZ0_

protected

Definition at line 172 of file Histos.h.

Referenced by bookTrackCands().

std::map<std::string, TH1F*> tmtt::Histos::hisZ0Res_

protected

Definition at line 178 of file Histos.h.

Referenced by bookTrackCands().

std::map<std::string, TProfile*> tmtt::Histos::hisZ0ResVsTrueEta_

protected

Definition at line 252 of file Histos.h.

Referenced by bookTrackFitting(), and fillTrackFitting().

std::map<std::string, TProfile*> tmtt::Histos::hisZ0ResVsTrueInvPt_

protected

Definition at line 258 of file Histos.h.

Referenced by bookTrackFitting(), and fillTrackFitting().

float tmtt::Histos::houghMinPt_

protected

Definition at line 120 of file Histos.h.

Referenced by Histos().

unsigned int tmtt::Histos::houghNbinsPhi_

protected

Definition at line 122 of file Histos.h.

Referenced by Histos().

unsigned int tmtt::Histos::houghNbinsPt_

protected

Definition at line 121 of file Histos.h.

Referenced by Histos().

std::map<unsigned int, float> tmtt::Histos::mapBarrelLayerMaxR_

protected

Definition at line 329 of file Histos.h.

Referenced by endJobAnalysis(), and fillInputData().

std::map<unsigned int, float> tmtt::Histos::mapBarrelLayerMinR_

protected

Definition at line 328 of file Histos.h.

Referenced by endJobAnalysis(), and fillInputData().

std::map<unsigned int, float> tmtt::Histos::mapEndcapWheelMaxZ_

protected

Definition at line 332 of file Histos.h.

Referenced by endJobAnalysis(), and fillInputData().

std::map<unsigned int, float> tmtt::Histos::mapEndcapWheelMinZ_

protected

Definition at line 331 of file Histos.h.

Referenced by endJobAnalysis(), and fillInputData().

std::map<unsigned int, float> tmtt::Histos::mapExtraAModuleTypeMaxR_

protected

Definition at line 341 of file Histos.h.

Referenced by endJobAnalysis(), and fillInputData().

std::map<unsigned int, float> tmtt::Histos::mapExtraAModuleTypeMaxZ_

protected

Definition at line 343 of file Histos.h.

Referenced by endJobAnalysis(), and fillInputData().

std::map<unsigned int, float> tmtt::Histos::mapExtraAModuleTypeMinR_

protected

Definition at line 340 of file Histos.h.

Referenced by endJobAnalysis(), and fillInputData().

std::map<unsigned int, float> tmtt::Histos::mapExtraAModuleTypeMinZ_

protected

Definition at line 342 of file Histos.h.

Referenced by endJobAnalysis(), and fillInputData().

std::map<unsigned int, float> tmtt::Histos::mapExtraBModuleTypeMaxR_

protected

Definition at line 345 of file Histos.h.

Referenced by endJobAnalysis(), and fillInputData().

std::map<unsigned int, float> tmtt::Histos::mapExtraBModuleTypeMaxZ_

protected

Definition at line 347 of file Histos.h.

Referenced by endJobAnalysis(), and fillInputData().

std::map<unsigned int, float> tmtt::Histos::mapExtraBModuleTypeMinR_

protected

Definition at line 344 of file Histos.h.

Referenced by endJobAnalysis(), and fillInputData().

std::map<unsigned int, float> tmtt::Histos::mapExtraBModuleTypeMinZ_

protected

Definition at line 346 of file Histos.h.

Referenced by endJobAnalysis(), and fillInputData().

std::map<unsigned int, float> tmtt::Histos::mapExtraCModuleTypeMaxR_

protected

Definition at line 349 of file Histos.h.

Referenced by endJobAnalysis(), and fillInputData().

std::map<unsigned int, float> tmtt::Histos::mapExtraCModuleTypeMaxZ_

protected

Definition at line 351 of file Histos.h.

Referenced by endJobAnalysis(), and fillInputData().

std::map<unsigned int, float> tmtt::Histos::mapExtraCModuleTypeMinR_

protected

Definition at line 348 of file Histos.h.

Referenced by endJobAnalysis(), and fillInputData().

std::map<unsigned int, float> tmtt::Histos::mapExtraCModuleTypeMinZ_

protected

Definition at line 350 of file Histos.h.

Referenced by endJobAnalysis(), and fillInputData().

std::map<unsigned int, float> tmtt::Histos::mapExtraDModuleTypeMaxR_

protected

Definition at line 353 of file Histos.h.

Referenced by endJobAnalysis(), and fillInputData().

std::map<unsigned int, float> tmtt::Histos::mapExtraDModuleTypeMaxZ_

protected

Definition at line 355 of file Histos.h.

Referenced by endJobAnalysis(), and fillInputData().

std::map<unsigned int, float> tmtt::Histos::mapExtraDModuleTypeMinR_

protected

Definition at line 352 of file Histos.h.

Referenced by endJobAnalysis(), and fillInputData().

std::map<unsigned int, float> tmtt::Histos::mapExtraDModuleTypeMinZ_

protected

Definition at line 354 of file Histos.h.

Referenced by endJobAnalysis(), and fillInputData().

std::map<unsigned int, float> tmtt::Histos::mapModuleTypeMaxR_

protected

Definition at line 336 of file Histos.h.

Referenced by endJobAnalysis(), and fillInputData().

std::map<unsigned int, float> tmtt::Histos::mapModuleTypeMaxZ_

protected

Definition at line 338 of file Histos.h.

Referenced by endJobAnalysis(), and fillInputData().

std::map<unsigned int, float> tmtt::Histos::mapModuleTypeMinR_

protected

Definition at line 335 of file Histos.h.

Referenced by endJobAnalysis(), and fillInputData().

std::map<unsigned int, float> tmtt::Histos::mapModuleTypeMinZ_

protected

Definition at line 337 of file Histos.h.

Referenced by endJobAnalysis(), and fillInputData().

unsigned int tmtt::Histos::numEtaRegions_

protected

Definition at line 119 of file Histos.h.

Referenced by bookTrackCands(), fillEtaPhiSectors(), and Histos().

std::map<std::string, unsigned int> tmtt::Histos::numFitAlgEff_

protected

Definition at line 320 of file Histos.h.

std::map<std::string, unsigned int> tmtt::Histos::numFitAlgEffPass_

protected

Definition at line 324 of file Histos.h.

std::map<std::string, unsigned int> tmtt::Histos::numFitPerfAlgEff_

protected

Definition at line 321 of file Histos.h.

std::map<std::string, unsigned int> tmtt::Histos::numFitPerfAlgEffPass_

protected

Definition at line 325 of file Histos.h.

unsigned int tmtt::Histos::numPhiSectors_

protected

Definition at line 118 of file Histos.h.

Referenced by fillEtaPhiSectors(), fillTrackFitting(), and Histos().

bool tmtt::Histos::oldSumW2opt_

protected

Definition at line 129 of file Histos.h.

Referenced by book(), and endJobAnalysis().

bool tmtt::Histos::printedGeomAnalysis_

protected

Definition at line 359 of file Histos.h.

std::map<std::string, TProfile*> tmtt::Histos::profDupFitTrksVsEta_

protected

Definition at line 261 of file Histos.h.

Referenced by bookTrackFitting(), and fillTrackFitting().

std::map<std::string, TProfile*> tmtt::Histos::profDupFitTrksVsInvPt_

protected

Definition at line 262 of file Histos.h.

Referenced by bookTrackFitting(), and fillTrackFitting().

std::map<std::string, TProfile*> tmtt::Histos::profDupTracksVsEta_

protected

Definition at line 165 of file Histos.h.

Referenced by bookTrackCands().

std::map<std::string, TProfile*> tmtt::Histos::profDupTracksVsInvPt_

protected

Definition at line 166 of file Histos.h.

Referenced by bookTrackCands().

std::map<std::string, TProfile*> tmtt::Histos::profFitChi2DofRphiVsInvPtMatched_

protected

Definition at line 243 of file Histos.h.

Referenced by bookTrackFitting(), and fillTrackFitting().

std::map<std::string, TProfile*> tmtt::Histos::profFitChi2DofRphiVsInvPtUnmatched_

protected

Definition at line 247 of file Histos.h.

Referenced by bookTrackFitting(), and fillTrackFitting().

std::map<std::string, TProfile*> tmtt::Histos::profMeanStubsPerLink_

protected

Definition at line 163 of file Histos.h.

Referenced by bookTrackCands().

std::map<std::string, TProfile*> tmtt::Histos::profNumFitTracks_

protected

Definition at line 214 of file Histos.h.

Referenced by bookTrackFitting(), fillTrackFitting(), and printFitTrackPerformance().

std::map<std::string, TProfile*> tmtt::Histos::profNumTrackCands_

protected

Definition at line 153 of file Histos.h.

Referenced by bookTrackCands(), and printTrackPerformance().

std::map<std::string, TProfile*> tmtt::Histos::profNumTracksVsEta_

protected

Definition at line 154 of file Histos.h.

Referenced by bookTrackCands().

std::map<std::string, TProfile*> tmtt::Histos::profStubsOnFitTracks_

protected

Definition at line 219 of file Histos.h.

Referenced by bookTrackFitting(), fillTrackFitting(), and printFitTrackPerformance().

std::map<std::string, TProfile*> tmtt::Histos::profStubsOnTracks_

protected

Definition at line 157 of file Histos.h.

Referenced by bookTrackCands(), and printTrackPerformance().

bool tmtt::Histos::ranRZfilter_

protected

Definition at line 126 of file Histos.h.

Referenced by book(), endJobAnalysis(), fill(), and Histos().

bool tmtt::Histos::resPlotOpt_

protected

Definition at line 127 of file Histos.h.

Referenced by Histos().

const Settings* tmtt::Histos::settings_

protected

Definition at line 116 of file Histos.h.

Referenced by enabled(), endJobAnalysis(), fillInputData(), fillTrackFitting(), printFitTrackPerformance(), and trackerGeometryAnalysis().

std::map<std::string, TEfficiency*> tmtt::Histos::teffAlgEffFitVsD0_

protected

Definition at line 313 of file Histos.h.

Referenced by plotTrackEffAfterFit().

std::map<std::string, TEfficiency*> tmtt::Histos::teffAlgEffFitVsEta_

protected

Definition at line 311 of file Histos.h.

Referenced by plotTrackEffAfterFit().

std::map<std::string, TEfficiency*> tmtt::Histos::teffAlgEffFitVsInvPt_

protected

Definition at line 310 of file Histos.h.

Referenced by plotTrackEffAfterFit().

std::map<std::string, TEfficiency*> tmtt::Histos::teffAlgEffFitVsPhi_

protected

Definition at line 312 of file Histos.h.

Referenced by plotTrackEffAfterFit().

std::map<std::string, TEfficiency*> tmtt::Histos::teffAlgEffFitVsZ0_

protected

Definition at line 314 of file Histos.h.

Referenced by plotTrackEffAfterFit().

std::map<std::string, TEfficiency*> tmtt::Histos::teffAlgEffVsD0_

protected

Definition at line 289 of file Histos.h.

Referenced by plotTrackEfficiency().

std::map<std::string, TEfficiency*> tmtt::Histos::teffAlgEffVsEta_

protected

Definition at line 293 of file Histos.h.

Referenced by plotTrackEfficiency().

std::map<std::string, TEfficiency*> tmtt::Histos::teffAlgEffVsInvPt_

protected

Definition at line 292 of file Histos.h.

Referenced by plotTrackEfficiency().

std::map<std::string, TEfficiency*> tmtt::Histos::teffAlgEffVsPhi_

protected

Definition at line 294 of file Histos.h.

Referenced by plotTrackEfficiency().

std::map<std::string, TEfficiency*> tmtt::Histos::teffAlgEffVsZ0_

protected

Definition at line 290 of file Histos.h.

Referenced by plotTrackEfficiency().

std::map<std::string, TEfficiency*> tmtt::Histos::teffEffFitVsD0_

protected

Definition at line 304 of file Histos.h.

Referenced by plotTrackEffAfterFit().

std::map<std::string, TEfficiency*> tmtt::Histos::teffEffFitVsEta_

protected

Definition at line 302 of file Histos.h.

Referenced by plotTrackEffAfterFit().

std::map<std::string, TEfficiency*> tmtt::Histos::teffEffFitVsInvPt_

protected

Definition at line 301 of file Histos.h.

Referenced by plotTrackEffAfterFit().

std::map<std::string, TEfficiency*> tmtt::Histos::teffEffFitVsPhi_

protected

Definition at line 303 of file Histos.h.

Referenced by plotTrackEffAfterFit().

std::map<std::string, TEfficiency*> tmtt::Histos::teffEffFitVsZ0_

protected

Definition at line 305 of file Histos.h.

Referenced by plotTrackEffAfterFit().

std::map<std::string, TEfficiency*> tmtt::Histos::teffEffVsD0_

protected

Definition at line 284 of file Histos.h.

Referenced by plotTrackEfficiency().

std::map<std::string, TEfficiency*> tmtt::Histos::teffEffVsEta_

protected

Definition at line 282 of file Histos.h.

Referenced by plotTrackEfficiency().

std::map<std::string, TEfficiency*> tmtt::Histos::teffEffVsInvPt_

protected

Definition at line 281 of file Histos.h.

Referenced by plotTrackEfficiency().

std::map<std::string, TEfficiency*> tmtt::Histos::teffEffVsPhi_

protected

Definition at line 283 of file Histos.h.

Referenced by plotTrackEfficiency().

std::map<std::string, TEfficiency*> tmtt::Histos::teffEffVsZ0_

protected

Definition at line 285 of file Histos.h.

Referenced by plotTrackEfficiency().

std::map<std::string, TEfficiency*> tmtt::Histos::teffPerfAlgEffFitVsEta_

protected

Definition at line 317 of file Histos.h.

Referenced by plotTrackEffAfterFit().

std::map<std::string, TEfficiency*> tmtt::Histos::teffPerfAlgEffFitVsInvPt_

protected

Definition at line 316 of file Histos.h.

Referenced by plotTrackEffAfterFit().

std::map<std::string, TEfficiency*> tmtt::Histos::teffPerfAlgEffVsEta_

protected

Definition at line 298 of file Histos.h.

Referenced by plotTrackEfficiency().

std::map<std::string, TEfficiency*> tmtt::Histos::teffPerfAlgEffVsInvPt_

protected

Definition at line 296 of file Histos.h.

Referenced by plotTrackEfficiency().

std::map<std::string, TEfficiency*> tmtt::Histos::teffPerfAlgEffVsPt_

protected

Definition at line 297 of file Histos.h.

std::map<std::string, TEfficiency*> tmtt::Histos::teffPerfEffFitVsEta_

protected

Definition at line 308 of file Histos.h.

Referenced by plotTrackEffAfterFit().

std::map<std::string, TEfficiency*> tmtt::Histos::teffPerfEffFitVsInvPt_

protected

Definition at line 307 of file Histos.h.

Referenced by plotTrackEffAfterFit().

std::map<std::string, TEfficiency*> tmtt::Histos::teffPerfEffVsEta_

protected

Definition at line 288 of file Histos.h.

Referenced by plotTrackEfficiency().

std::map<std::string, TEfficiency*> tmtt::Histos::teffPerfEffVsInvPt_

protected

Definition at line 287 of file Histos.h.

Referenced by plotTrackEfficiency().

std::vector<std::string> tmtt::Histos::trackFitters_

protected

Definition at line 124 of file Histos.h.

Referenced by bookTrackFitting(), endJobAnalysis(), fillTrackFitting(), and Histos().

std::vector<std::string> tmtt::Histos::useRZfilter_

protected

Definition at line 125 of file Histos.h.

Referenced by Histos(), and printFitTrackPerformance().