#include <TrackTypeMonitor.h>

Inheritance diagram for TrackTypeMonitor:

Public Member Functions
	TrackTypeMonitor (const edm::ParameterSet &)

Public Member Functions inherited from DQMEDAnalyzer
virtual void	beginRun (edm::Run const &, edm::EventSetup const &) final

virtual void	beginStream (edm::StreamID id) final

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

	DQMEDAnalyzer (void)

virtual void	endLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, dqmDetails::NoCache *) const final

virtual void	endRunSummary (edm::Run const &, edm::EventSetup const &, dqmDetails::NoCache *) const final

uint32_t	streamId () const

Public Member Functions inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
	EDAnalyzer ()=default

Public Member Functions inherited from edm::stream::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDAnalyzerBase ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Protected Member Functions
void	analyze (edm::Event const &iEvent, edm::EventSetup const &iSetup) override

void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &)

void	endLuminosityBlock (edm::LuminosityBlock const &lumiSeg, edm::EventSetup const &eSetup) override

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Private Member Functions
void	fillHistograms (const reco::Track &track, int indx)

Private Attributes
const edm::InputTag	bsTag_

const edm::EDGetTokenT < reco::BeamSpot >	bsToken_

MonitorElement *	dphiH_

MonitorElement *	drH_

const edm::InputTag	electronTag_

const edm::EDGetTokenT < reco::GsfElectronCollection >	electronToken_

std::string	folderName_

MonitorElement *	hcounterH_

unsigned long long	m_cacheID_

std::string	moduleName_

const edm::InputTag	muonTag_

const edm::EDGetTokenT < reco::MuonCollection >	muonToken_

std::vector< MonitorElement * >	nTracksHList_

edm::ParameterSet	parameters_

std::vector< MonitorElement * >	trackChi2bynDOFHList_

std::vector< MonitorElement * >	trackdzHList_

std::vector< MonitorElement * >	trackEtaHList_

std::vector< MonitorElement * >	trackPhiHList_

std::vector< MonitorElement * >	trackPHList_

std::vector< MonitorElement * >	trackPterrHList_

std::vector< MonitorElement * >	trackPtHList_

std::vector< MonitorElement * >	trackqOverpHList_

const std::string	trackQuality_

const edm::InputTag	trackTag_

const edm::EDGetTokenT < reco::TrackCollection >	trackToken_

bool	verbose_

const edm::InputTag	vertexTag_

const edm::EDGetTokenT < reco::VertexCollection >	vertexToken_

Additional Inherited Members
Public Types inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
typedef CacheContexts< T...>	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T...>	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT < LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Public Types inherited from edm::stream::EDAnalyzerBase
typedef EDAnalyzerAdaptorBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from DQMEDAnalyzer
static std::shared_ptr < dqmDetails::NoCache >	globalBeginLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *)

static std::shared_ptr < dqmDetails::NoCache >	globalBeginRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *)

static void	globalEndLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const , dqmDetails::NoCache )

static void	globalEndRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const , dqmDetails::NoCache )

Static Public Member Functions inherited from edm::stream::EDAnalyzerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Detailed Description

Definition at line 26 of file TrackTypeMonitor.h.

Constructor & Destructor Documentation

TrackTypeMonitor::TrackTypeMonitor ( const edm::ParameterSet & ps )

Definition at line 37 of file TrackTypeMonitor.cc.

References nTracksHList_, trackChi2bynDOFHList_, trackEtaHList_, trackPhiHList_, trackPHList_, trackPterrHList_, trackPtHList_, and trackqOverpHList_.

                                                            : 
   parameters_(ps),
   moduleName_(parameters_.getUntrackedParameter<std::string>("ModuleName", "TrackTypeMonitor")),
   folderName_(parameters_.getUntrackedParameter<std::string>("FolderName", "highPurityTracks")),
   verbose_(parameters_.getUntrackedParameter<bool>("verbose", false)),
   muonTag_(ps.getUntrackedParameter<edm::InputTag>("muonInputTag", edm::InputTag("muons"))),
   electronTag_(ps.getUntrackedParameter<edm::InputTag>("electronInputTag", edm::InputTag("gedGsfElectrons"))),
   trackTag_(parameters_.getUntrackedParameter<edm::InputTag>("trackInputTag", edm::InputTag("generalTracks"))),
   bsTag_(parameters_.getUntrackedParameter<edm::InputTag>("offlineBeamSpot", edm::InputTag("offlineBeamSpot"))),
   vertexTag_(parameters_.getUntrackedParameter<edm::InputTag>("vertexTag", edm::InputTag("offlinePrimaryVertices"))),
   electronToken_(consumes<reco::GsfElectronCollection>(electronTag_)),
   muonToken_(consumes<reco::MuonCollection>(muonTag_)),
   trackToken_(consumes<reco::TrackCollection>(trackTag_)),
   bsToken_(consumes<reco::BeamSpot>(bsTag_)),
   vertexToken_(consumes<reco::VertexCollection>(vertexTag_)),
   trackQuality_(parameters_.getUntrackedParameter<std::string>("trackQuality", "highPurity"))
 {
   trackEtaHList_.clear();
   trackPhiHList_.clear();
   trackPHList_.clear();
   trackPtHList_.clear();
   trackPterrHList_.clear();
   trackqOverpHList_.clear();
   trackChi2bynDOFHList_.clear();
   nTracksHList_.clear();
 }

Member Function Documentation

void TrackTypeMonitor::analyze	(	edm::Event const &	iEvent,
		edm::EventSetup const &	iSetup
	)

overrideprotectedvirtual

Implements edm::stream::EDAnalyzerBase.

Definition at line 234 of file TrackTypeMonitor.cc.

References SiPixelRawToDigiRegional_cfi::beamSpot, bsToken_, reco::TrackBase::chi2(), reco::TrackBase::d0(), SiPixelRawToDigiRegional_cfi::deltaPhi, deltaR(), dphiH_, drH_, reco::TrackBase::dxy(), reco::TrackBase::dz(), electronToken_, reco::TrackBase::eta(), MonitorElement::Fill(), fillHistograms(), edm::Event::getByToken(), hcounterH_, reco::TrackBase::hitPattern(), edm::Ref< C, T, F >::isNonnull(), edm::HandleBase::isValid(), bookConverter::max, muonToken_, ndof, reco::TrackBase::ndof(), nTracksHList_, reco::HitPattern::numberOfValidPixelHits(), reco::HitPattern::numberOfValidStripHits(), getGTfromDQMFile::obj, reco::TrackBase::phi(), reco::Vertex::position(), HLT_25ns14e33_v1_cff::quality, reco::TrackBase::qualityByName(), python.multivaluedict::sort(), reco::GsfElectron::PflowIsolationVariables::sumChargedHadronPt, reco::MuonPFIsolation::sumChargedParticlePt, reco::MuonPFIsolation::sumNeutralHadronEt, reco::GsfElectron::PflowIsolationVariables::sumNeutralHadronEt, reco::MuonPFIsolation::sumPhotonEt, reco::GsfElectron::PflowIsolationVariables::sumPhotonEt, reco::MuonPFIsolation::sumPUPt, reco::GsfElectron::PflowIsolationVariables::sumPUPt, trackQuality_, testEve_cfg::tracks, trackToken_, verbose_, and vertexToken_.

                                                                                   {
 
   // read the beam spot
   edm::Handle<reco::BeamSpot> beamSpot;
   iEvent.getByToken(bsToken_, beamSpot);
 
   std::vector<const reco::Track*> isoTrkList;
 
   // muons
   edm::Handle<reco::MuonCollection> muonColl;
   iEvent.getByToken(muonToken_, muonColl);
 
   if (muonColl.isValid()) {
     for (auto const& muo : *muonColl) {
       if (muo.isGlobalMuon() && muo.isPFMuon() && std::fabs(muo.eta()) <= 2.1 && muo.pt() > 1) {
     reco::TrackRef gtrkref = muo.globalTrack();
     if (!gtrkref.isNonnull()) continue;     
     const reco::Track* gtk = &(*gtrkref);
     double chi2 = gtk->chi2();
     double ndof = gtk->ndof();
     double chbyndof = (ndof > 0) ? chi2/ndof : 0;
     
     const reco::HitPattern& hitp = gtk->hitPattern();
     int nPixelHits = hitp.numberOfValidPixelHits();
     int nStripHits = hitp.numberOfValidStripHits();
     
     reco::TrackRef itrkref = muo.innerTrack();  // tracker segment only
     if (!itrkref.isNonnull()) continue;
     const reco::Track* tk = &(*itrkref);
     double trkd0 = tk->d0();
     double trkdz = tk->dz();
     if (beamSpot.isValid()) {
       trkd0 = -(tk->dxy(beamSpot->position()));
       trkdz = tk->dz(beamSpot->position());
     }
    
     // Hits/section in the muon chamber
     int nChambers = muo.numberOfChambers();
     int nMatches = muo.numberOfMatches();
     int nMatchedStations = muo.numberOfMatchedStations();
      
     // PF Isolation      
     const reco::MuonPFIsolation& pfIso04 = muo.pfIsolationR04();
     double absiso = pfIso04.sumChargedParticlePt + std::max(0.0, pfIso04.sumNeutralHadronEt + pfIso04.sumPhotonEt - 0.5 * pfIso04.sumPUPt);
     if (chbyndof < 10
         &&  std::fabs(trkd0) < 0.02
         &&  std::fabs(trkdz) < 20
         &&  nPixelHits > 1
         &&  nStripHits > 8
         &&  nChambers > 2
         &&  nMatches > 2 
         &&  nMatchedStations > 2
         &&  absiso/muo.pt() < 0.3)
       {
         isoTrkList.push_back(gtk);
         fillHistograms(*gtk, 0);
       }
       }
     }
   }
   
   // electrons
   edm::Handle<reco::GsfElectronCollection> electronColl;
   iEvent.getByToken(electronToken_, electronColl);
 
   if (electronColl.isValid()) {
     for (auto const& ele : *electronColl) {
       if (!ele.ecalDriven()) continue;
       if (ele.pt() < 5) continue;
       hcounterH_->Fill(0);
 
       double hOverE = ele.hadronicOverEm();
       double sigmaee = ele.sigmaIetaIeta();
       double deltaPhiIn = ele.deltaPhiSuperClusterTrackAtVtx();
       double deltaEtaIn = ele.deltaEtaSuperClusterTrackAtVtx();
 
       if (ele.isEB()) {
     if (std::fabs(deltaPhiIn) >= .15
         && std::fabs(deltaEtaIn) >= .007
         && hOverE >= .12
         && sigmaee >= .01) continue;
       }
       else if (ele.isEE()) {
     if (std::fabs(deltaPhiIn) >= .10
         && std::fabs(deltaEtaIn) >= .009
         && hOverE >= .10
         && sigmaee >= .03) continue;
       }
       hcounterH_->Fill(1);
 
       reco::GsfTrackRef gsftrk = ele.gsfTrack();
       if (!gsftrk.isNonnull()) continue;
       const reco::GsfTrack* trk = &(*gsftrk);     
       double trkd0 = trk->d0();
       double trkdz = trk->dz();
       if (beamSpot.isValid()) {
     trkd0 = -(trk->dxy(beamSpot->position()));
     trkdz = trk->dz(beamSpot->position());
       }
       double chi2 = trk->chi2();
       double ndof = trk->ndof();
       double chbyndof = (ndof > 0) ? chi2/ndof : 0;
       if (chbyndof >= 10) continue;
       hcounterH_->Fill(2);
 
       if (std::fabs(trkd0) >= 0.02 || std::fabs(trkdz) >= 20) continue;
       hcounterH_->Fill(3);
 
       const reco::HitPattern& hitp = trk->hitPattern();
       int nPixelHits = hitp.numberOfValidPixelHits();
       if (nPixelHits < 1) continue;
       hcounterH_->Fill(4);
 
       int nStripHits = hitp.numberOfValidStripHits();
       if (nStripHits < 8) continue;
       hcounterH_->Fill(5);
 
       reco::GsfElectron::PflowIsolationVariables pfIso = ele.pfIsolationVariables();
       float absiso = pfIso.sumChargedHadronPt + std::max(0.0, pfIso.sumNeutralHadronEt + pfIso.sumPhotonEt - 0.5 * pfIso.sumPUPt);
       float eiso = absiso/ele.pt();
       if (eiso > 0.2) continue;
       hcounterH_->Fill(6);
     
       isoTrkList.push_back(trk);
       fillHistograms(*trk, 0);
     }
   }
   nTracksHList_.at(0)->Fill(isoTrkList.size());
 
   // Read track collection
   edm::Handle<reco::TrackCollection> tracks;
   iEvent.getByToken(trackToken_, tracks);
   
   // Read vertex collection
   edm::Handle<reco::VertexCollection> vertexColl;
   iEvent.getByToken(vertexToken_, vertexColl);
 
   if (tracks.isValid()) {
     std::vector<const reco::Track*> nisoTrkList;
     const reco::Vertex& vit = vertexColl->front(); // Highest sumPt vertex
     reco::Track::TrackQuality quality = reco::Track::qualityByName(trackQuality_); 
     
     int nmatch = 0;
     int nTracks = 0, nMBTracks = 0;
     for (auto const& track: *tracks) {
       if (!track.quality(quality)) continue;
       ++nTracks;
       fillHistograms(track, 3);
 
       // now classify primary and underlying tracks using vertex information
       double dxy = track.dxy(vit.position());
       double dz  = track.dz(vit.position());
       if (std::fabs(dxy) < 0.02 && std::fabs(dz) < 20) { // primary tracks
     // remove tracks associated to the identified particles
         double drmin = 999;
     for (auto const& tk : isoTrkList) {
           double dr = deltaR(track.eta(), track.phi(), tk->eta(), tk->phi());
           if (dr < drmin) drmin = dr;
     }
     if (drmin < 0.01) {
           if (verbose_) {
         edm::LogInfo("TrackTypeMonitor") << " Match: " << ++nmatch << " drmin = " << drmin;
         edm::LogInfo("TrackTypeMonitor") << " Track Pt: " << track.pt() << " Eta: " << track.eta() << " Phi: " << track.phi();
         for (auto const& isotk : isoTrkList) {
           edm::LogInfo("TrackTypeMonitor") << " Lepton Pt: " << isotk->pt() << " Eta: " << isotk->eta() << " Phi: " << isotk->phi();
         }
           }
           continue;
         }
         
     fillHistograms(track, 1);
     nisoTrkList.push_back(&track);
       } 
       else {
         ++nMBTracks;
     fillHistograms(track, 2); //non-primary tracks
       }     
     }
     nTracksHList_.at(1)->Fill(nisoTrkList.size());
     nTracksHList_.at(2)->Fill(nMBTracks);
     nTracksHList_.at(3)->Fill(nTracks);
     
     std::sort(nisoTrkList.begin(), nisoTrkList.end(), PtComparator<reco::Track>());
     const reco::Track* isoTrk = isoTrkList.at(0);
     for (auto const& obj : nisoTrkList) {
       if (obj->pt() > 5) {
     double dphi = deltaPhi(isoTrk->phi(), obj->phi());
     dphiH_->Fill(std::fabs(dphi));
 
     double dr = deltaR(isoTrk->eta(), isoTrk->phi(), obj->eta(), obj->phi());
     drH_->Fill(dr);
       }
     }
   }  
 }

void TrackTypeMonitor::bookHistograms	(	DQMStore::IBooker &	iBook,
		edm::Run const &	iRun,
		edm::EventSetup const &	iSetup
	)

protectedvirtual

Implements DQMEDAnalyzer.

Definition at line 63 of file TrackTypeMonitor.cc.

References DQMStore::IBooker::book1D(), dphiH_, drH_, folderName_, edm::ParameterSet::getParameter(), hcounterH_, moduleName_, nTracksHList_, parameters_, DQMStore::IBooker::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, trackChi2bynDOFHList_, trackdzHList_, trackEtaHList_, trackPhiHList_, trackPHList_, trackPterrHList_, trackPtHList_, and trackqOverpHList_.

                                                                                                              {
   edm::ParameterSet TrackEtaHistoPar = parameters_.getParameter<edm::ParameterSet>("TrackEtaPar");
   edm::ParameterSet TrackPhiHistoPar = parameters_.getParameter<edm::ParameterSet>("TrackPhiPar");
   edm::ParameterSet TrackPHistoPar = parameters_.getParameter<edm::ParameterSet>("TrackPPar");
   edm::ParameterSet TrackPtHistoPar = parameters_.getParameter<edm::ParameterSet>("TrackPtPar");
   edm::ParameterSet TrackPterrHistoPar = parameters_.getParameter<edm::ParameterSet>("TrackPterrPar");
   edm::ParameterSet TrackqOverpHistoPar = parameters_.getParameter<edm::ParameterSet>("TrackqOverpPar");
   edm::ParameterSet TrackdzHistoPar = parameters_.getParameter<edm::ParameterSet>("TrackdzPar");
   edm::ParameterSet TrackChi2bynDOFHistoPar = parameters_.getParameter<edm::ParameterSet>("TrackChi2bynDOFPar");
   edm::ParameterSet nTracksHistoPar = parameters_.getParameter<edm::ParameterSet>("nTracksPar");
 
   std::string currentFolder = moduleName_ + "/" + folderName_ ;
   iBook.setCurrentFolder(currentFolder.c_str());
 
   trackEtaHList_.push_back(iBook.book1D("trackEtaIso", "Isolated Track Eta", 
                     TrackEtaHistoPar.getParameter<int32_t>("Xbins"), 
                     TrackEtaHistoPar.getParameter<double>("Xmin"), 
                     TrackEtaHistoPar.getParameter<double>("Xmax")));
   trackEtaHList_.push_back(iBook.book1D("trackEtaNoIso", "NonIsolated Track Eta", 
                     TrackEtaHistoPar.getParameter<int32_t>("Xbins"), 
                     TrackEtaHistoPar.getParameter<double>("Xmin"), 
                     TrackEtaHistoPar.getParameter<double>("Xmax")));
   trackEtaHList_.push_back(iBook.book1D("trackEtaUL", "Underlying Track Eta", 
                     TrackEtaHistoPar.getParameter<int32_t>("Xbins"), 
                     TrackEtaHistoPar.getParameter<double>("Xmin"), 
                     TrackEtaHistoPar.getParameter<double>("Xmax")));
   trackEtaHList_.push_back(iBook.book1D("trackEtaALL", "All Track Eta", 
                     TrackEtaHistoPar.getParameter<int32_t>("Xbins"),
                                         TrackEtaHistoPar.getParameter<double>("Xmin"),
                                         TrackEtaHistoPar.getParameter<double>("Xmax")));
 
   trackPhiHList_.push_back(iBook.book1D("trackPhiIso", "Isolated Track Phi", 
                     TrackPhiHistoPar.getParameter<int32_t>("Xbins"),
                                         TrackPhiHistoPar.getParameter<double>("Xmin"),
                                         TrackPhiHistoPar.getParameter<double>("Xmax")));
   trackPhiHList_.push_back(iBook.book1D("trackPhiNonIso", "NonIsolated Track Phi", 
                     TrackPhiHistoPar.getParameter<int32_t>("Xbins"),
                     TrackPhiHistoPar.getParameter<double>("Xmin"),
                     TrackPhiHistoPar.getParameter<double>("Xmax")));
   trackPhiHList_.push_back(iBook.book1D("trackPhiUL", "Underlying Track Phi", 
                     TrackPhiHistoPar.getParameter<int32_t>("Xbins"),
                     TrackPhiHistoPar.getParameter<double>("Xmin"),
                     TrackPhiHistoPar.getParameter<double>("Xmax")));
   trackPhiHList_.push_back(iBook.book1D("trackPhiALL", "All Track Phi", 
                     TrackPhiHistoPar.getParameter<int32_t>("Xbins"),
                                         TrackPhiHistoPar.getParameter<double>("Xmin"),
                                         TrackPhiHistoPar.getParameter<double>("Xmax")));
 
   trackPHList_.push_back(iBook.book1D("trackPIso", "Isolated Track P", 
                       TrackPHistoPar.getParameter<int32_t>("Xbins"),
                       TrackPHistoPar.getParameter<double>("Xmin"),
                       TrackPHistoPar.getParameter<double>("Xmax")));
   trackPHList_.push_back(iBook.book1D("trackPNonIso", "NonIsolated Track P", 
                       TrackPHistoPar.getParameter<int32_t>("Xbins"),
                       TrackPHistoPar.getParameter<double>("Xmin"),
                       TrackPHistoPar.getParameter<double>("Xmax")));
   trackPHList_.push_back(iBook.book1D("trackPUL", "Underlying Track P", 
                       TrackPHistoPar.getParameter<int32_t>("Xbins"),
                       TrackPHistoPar.getParameter<double>("Xmin"),
                       TrackPHistoPar.getParameter<double>("Xmax")));
   trackPHList_.push_back(iBook.book1D("trackPALL", "All Track P", 
                       TrackPHistoPar.getParameter<int32_t>("Xbins"),
                                       TrackPHistoPar.getParameter<double>("Xmin"),
                                       TrackPHistoPar.getParameter<double>("Xmax")));
 
   trackPtHList_.push_back(iBook.book1D("trackPtIsolated", "Isolated Track Pt", 
                        TrackPtHistoPar.getParameter<int32_t>("Xbins"),
                       TrackPtHistoPar.getParameter<double>("Xmin"),
                       TrackPtHistoPar.getParameter<double>("Xmax")));
   trackPtHList_.push_back(iBook.book1D("trackPtNonIsolated", "NonIsolated Track Pt", 
                        TrackPtHistoPar.getParameter<int32_t>("Xbins"),
                                       TrackPtHistoPar.getParameter<double>("Xmin"),
                                       TrackPtHistoPar.getParameter<double>("Xmax")));
   trackPtHList_.push_back(iBook.book1D("trackPtUL", "Underlying Track Pt", 
                        TrackPtHistoPar.getParameter<int32_t>("Xbins"),
                                       TrackPtHistoPar.getParameter<double>("Xmin"),
                                       TrackPtHistoPar.getParameter<double>("Xmax")));
   trackPtHList_.push_back(iBook.book1D("trackPtALL", "All Track Pt", 
                        TrackPtHistoPar.getParameter<int32_t>("Xbins"),
                        TrackPtHistoPar.getParameter<double>("Xmin"),
                        TrackPtHistoPar.getParameter<double>("Xmax")));
   
   trackPterrHList_.push_back(iBook.book1D("trackPterrIsolated", "Isolated Track Pterr", 
                       TrackPterrHistoPar.getParameter<int32_t>("Xbins"),
                       TrackPterrHistoPar.getParameter<double>("Xmin"),
                       TrackPterrHistoPar.getParameter<double>("Xmax")));
   trackPterrHList_.push_back(iBook.book1D("trackPterrNonIsolated", "NonIsolated Track Pterr", 
                       TrackPterrHistoPar.getParameter<int32_t>("Xbins"),
                       TrackPterrHistoPar.getParameter<double>("Xmin"),
                       TrackPterrHistoPar.getParameter<double>("Xmax")));
   trackPterrHList_.push_back(iBook.book1D("trackPterrUL", "Underlying Track Pterr", 
                       TrackPterrHistoPar.getParameter<int32_t>("Xbins"),
                       TrackPterrHistoPar.getParameter<double>("Xmin"),
                       TrackPterrHistoPar.getParameter<double>("Xmax")));
   trackPterrHList_.push_back(iBook.book1D("trackPterrALL", "All Track Pterr", 
                       TrackPterrHistoPar.getParameter<int32_t>("Xbins"),
                       TrackPterrHistoPar.getParameter<double>("Xmin"),
                       TrackPterrHistoPar.getParameter<double>("Xmax")));
  
   trackqOverpHList_.push_back(iBook.book1D("trackqOverpIsolated", "Isolated Track qOverp", 
                        TrackqOverpHistoPar.getParameter<int32_t>("Xbins"),
                        TrackqOverpHistoPar.getParameter<double>("Xmin"),
                        TrackqOverpHistoPar.getParameter<double>("Xmax")));
   trackqOverpHList_.push_back(iBook.book1D("trackqOverpNonIsolated", "NonIsolated Track qOverp", 
                        TrackqOverpHistoPar.getParameter<int32_t>("Xbins"),
                        TrackqOverpHistoPar.getParameter<double>("Xmin"),
                        TrackqOverpHistoPar.getParameter<double>("Xmax")));
   trackqOverpHList_.push_back(iBook.book1D("trackqOverpUL", "Underlying Track qOverp", 
                        TrackqOverpHistoPar.getParameter<int32_t>("Xbins"),
                        TrackqOverpHistoPar.getParameter<double>("Xmin"),
                        TrackqOverpHistoPar.getParameter<double>("Xmax")));
   trackqOverpHList_.push_back(iBook.book1D("trackqOverpALL", "All Track qOverp", 
                        TrackqOverpHistoPar.getParameter<int32_t>("Xbins"),
                        TrackqOverpHistoPar.getParameter<double>("Xmin"),
                        TrackqOverpHistoPar.getParameter<double>("Xmax")));
 
   trackdzHList_.push_back(iBook.book1D("trackdzIsolated", "Isolated Track dz", 
                        TrackdzHistoPar.getParameter<int32_t>("Xbins"),
                        TrackdzHistoPar.getParameter<double>("Xmin"),
                        TrackdzHistoPar.getParameter<double>("Xmax")));
   trackdzHList_.push_back(iBook.book1D("trackdzNonIsolated", "NonIsolated Track dz", 
                        TrackdzHistoPar.getParameter<int32_t>("Xbins"),
                        TrackdzHistoPar.getParameter<double>("Xmin"),
                        TrackdzHistoPar.getParameter<double>("Xmax")));
   trackdzHList_.push_back(iBook.book1D("trackdzUL", "Underlying Track dz", 
                        TrackdzHistoPar.getParameter<int32_t>("Xbins"),
                        TrackdzHistoPar.getParameter<double>("Xmin"),
                        TrackdzHistoPar.getParameter<double>("Xmax")));
   trackdzHList_.push_back(iBook.book1D("trackdzALL", "All Track dz", 
                        TrackdzHistoPar.getParameter<int32_t>("Xbins"),
                        TrackdzHistoPar.getParameter<double>("Xmin"),
                        TrackdzHistoPar.getParameter<double>("Xmax")));
 
   trackChi2bynDOFHList_.push_back(iBook.book1D("trackChi2bynDOFIsolated", "Isolated Track Chi2bynDOF", 
                            TrackChi2bynDOFHistoPar.getParameter<int32_t>("Xbins"),
                            TrackChi2bynDOFHistoPar.getParameter<double>("Xmin"),
                            TrackChi2bynDOFHistoPar.getParameter<double>("Xmax")));
   trackChi2bynDOFHList_.push_back(iBook.book1D("trackChi2bynDOFNonIsolated", "NonIsolated Track Chi2bynDOF", 
                            TrackChi2bynDOFHistoPar.getParameter<int32_t>("Xbins"),
                            TrackChi2bynDOFHistoPar.getParameter<double>("Xmin"),
                            TrackChi2bynDOFHistoPar.getParameter<double>("Xmax")));
   trackChi2bynDOFHList_.push_back(iBook.book1D("trackChi2bynDOFUL", "Underlying Track Chi2bynDOF", 
                            TrackChi2bynDOFHistoPar.getParameter<int32_t>("Xbins"),
                            TrackChi2bynDOFHistoPar.getParameter<double>("Xmin"),
                            TrackChi2bynDOFHistoPar.getParameter<double>("Xmax")));
   trackChi2bynDOFHList_.push_back(iBook.book1D("trackChi2bynDOFAll", "All Track Chi2bynDOF", 
                            TrackChi2bynDOFHistoPar.getParameter<int32_t>("Xbins"),
                            TrackChi2bynDOFHistoPar.getParameter<double>("Xmin"),
                            TrackChi2bynDOFHistoPar.getParameter<double>("Xmax")));
 
   nTracksHList_.push_back(iBook.book1D("nTracksIsolated", "Isolated Track nTracks", 
                        nTracksHistoPar.getParameter<int32_t>("Xbins"),
                        nTracksHistoPar.getParameter<double>("Xmin"),
                        nTracksHistoPar.getParameter<double>("Xmax")));
   nTracksHList_.push_back(iBook.book1D("nTracksNonIsolated", "NonIsolated Track nTracks", 
                        nTracksHistoPar.getParameter<int32_t>("Xbins"),
                        nTracksHistoPar.getParameter<double>("Xmin"),
                        nTracksHistoPar.getParameter<double>("Xmax")));
   nTracksHList_.push_back(iBook.book1D("nTracksUL", "Underlying Track nTracks", 
                        nTracksHistoPar.getParameter<int32_t>("Xbins"),
                        nTracksHistoPar.getParameter<double>("Xmin"),
                        nTracksHistoPar.getParameter<double>("Xmax")));
   nTracksHList_.push_back(iBook.book1D("nTracksAll", "All Track nTracks", 
                        nTracksHistoPar.getParameter<int32_t>("Xbins"),
                        nTracksHistoPar.getParameter<double>("Xmin"),
                        nTracksHistoPar.getParameter<double>("Xmax")));
 
   hcounterH_ = iBook.book1D("hcounter", "hcounter", 7, -0.5, 6.5);
   dphiH_ = iBook.book1D("dphi", "dphi", 100, 0, 7);
   drH_ = iBook.book1D("dr", "dr", 100, 0, 6);
 }

void TrackTypeMonitor::endLuminosityBlock	(	edm::LuminosityBlock const &	lumiSeg,
		edm::EventSetup const &	eSetup
	)

overrideprotectedvirtual

Reimplemented from edm::stream::EDAnalyzerBase.

Definition at line 456 of file TrackTypeMonitor.cc.

456 {

457 }

void TrackTypeMonitor::fillHistograms	(	const reco::Track &	track,
		int	indx
	)

private

Definition at line 429 of file TrackTypeMonitor.cc.

References reco::TrackBase::chi2(), reco::TrackBase::dz(), reco::TrackBase::eta(), ndof, reco::TrackBase::ndof(), reco::TrackBase::p(), reco::TrackBase::phi(), reco::TrackBase::pt(), reco::TrackBase::ptError(), reco::TrackBase::qoverp(), trackChi2bynDOFHList_, trackdzHList_, trackEtaHList_, trackPhiHList_, trackPHList_, trackPterrHList_, trackPtHList_, and trackqOverpHList_.

Referenced by analyze().

                                                                       {
   if (indx >= 0 && indx < static_cast<int>(trackEtaHList_.size())) 
     trackEtaHList_.at(indx)->Fill(track.eta());
   
   if (indx >= 0 && indx < static_cast<int>(trackPhiHList_.size())) 
     trackPhiHList_.at(indx)->Fill(track.phi());
 
   if (indx >= 0 && indx < static_cast<int>(trackPHList_.size())) 
     trackPHList_.at(indx)->Fill(track.p());
 
   if (indx >= 0 && indx < static_cast<int>(trackPtHList_.size()))
     trackPtHList_.at(indx)->Fill(track.pt());
 
   if (indx >= 0 || indx < static_cast<int>(trackPterrHList_.size()))
     trackPterrHList_.at(indx)->Fill(track.ptError());
 
   if (indx >= 0 || indx < static_cast<int>(trackqOverpHList_.size()))
     trackqOverpHList_.at(indx)->Fill(track.qoverp());
 
   if (indx >= 0 || indx < static_cast<int>(trackqOverpHList_.size()))
     trackdzHList_.at(indx)->Fill(track.dz());
 
   double chi2 = track.chi2();
   double ndof = track.ndof();
   double chbyndof = (ndof > 0) ? chi2/ndof : 0;
   trackChi2bynDOFHList_.at(indx)->Fill(chbyndof);  
 }

Member Data Documentation

const edm::InputTag TrackTypeMonitor::bsTag_

private

Definition at line 49 of file TrackTypeMonitor.h.

const edm::EDGetTokenT<reco::BeamSpot> TrackTypeMonitor::bsToken_

private

Definition at line 55 of file TrackTypeMonitor.h.

Referenced by analyze().

MonitorElement* TrackTypeMonitor::dphiH_

private

Definition at line 71 of file TrackTypeMonitor.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* TrackTypeMonitor::drH_

private

Definition at line 72 of file TrackTypeMonitor.h.

Referenced by analyze(), and bookHistograms().

const edm::InputTag TrackTypeMonitor::electronTag_

private

Definition at line 47 of file TrackTypeMonitor.h.

const edm::EDGetTokenT<reco::GsfElectronCollection> TrackTypeMonitor::electronToken_

private

Definition at line 52 of file TrackTypeMonitor.h.

Referenced by analyze().

std::string TrackTypeMonitor::folderName_

private

Definition at line 43 of file TrackTypeMonitor.h.

Referenced by bookHistograms().

MonitorElement* TrackTypeMonitor::hcounterH_

private

Definition at line 70 of file TrackTypeMonitor.h.

Referenced by analyze(), and bookHistograms().

unsigned long long TrackTypeMonitor::m_cacheID_

private

Definition at line 74 of file TrackTypeMonitor.h.

std::string TrackTypeMonitor::moduleName_

private

Definition at line 42 of file TrackTypeMonitor.h.

Referenced by bookHistograms().

const edm::InputTag TrackTypeMonitor::muonTag_

private

Definition at line 46 of file TrackTypeMonitor.h.

const edm::EDGetTokenT<reco::MuonCollection> TrackTypeMonitor::muonToken_

private

Definition at line 53 of file TrackTypeMonitor.h.

Referenced by analyze().

std::vector<MonitorElement*> TrackTypeMonitor::nTracksHList_

private

Definition at line 67 of file TrackTypeMonitor.h.

Referenced by analyze(), bookHistograms(), and TrackTypeMonitor().

edm::ParameterSet TrackTypeMonitor::parameters_

private

Definition at line 40 of file TrackTypeMonitor.h.

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

std::vector<MonitorElement*> TrackTypeMonitor::trackChi2bynDOFHList_

private

Definition at line 66 of file TrackTypeMonitor.h.

Referenced by bookHistograms(), fillHistograms(), and TrackTypeMonitor().

std::vector<MonitorElement*> TrackTypeMonitor::trackdzHList_

private

Definition at line 68 of file TrackTypeMonitor.h.

Referenced by bookHistograms(), and fillHistograms().

std::vector<MonitorElement*> TrackTypeMonitor::trackEtaHList_

private

Definition at line 60 of file TrackTypeMonitor.h.

Referenced by bookHistograms(), fillHistograms(), and TrackTypeMonitor().

std::vector<MonitorElement*> TrackTypeMonitor::trackPhiHList_

private

Definition at line 61 of file TrackTypeMonitor.h.

Referenced by bookHistograms(), fillHistograms(), and TrackTypeMonitor().

std::vector<MonitorElement*> TrackTypeMonitor::trackPHList_

private

Definition at line 62 of file TrackTypeMonitor.h.

Referenced by bookHistograms(), fillHistograms(), and TrackTypeMonitor().

std::vector<MonitorElement*> TrackTypeMonitor::trackPterrHList_

private

Definition at line 64 of file TrackTypeMonitor.h.

Referenced by bookHistograms(), fillHistograms(), and TrackTypeMonitor().

std::vector<MonitorElement*> TrackTypeMonitor::trackPtHList_

private

Definition at line 63 of file TrackTypeMonitor.h.

Referenced by bookHistograms(), fillHistograms(), and TrackTypeMonitor().

std::vector<MonitorElement*> TrackTypeMonitor::trackqOverpHList_

private

Definition at line 65 of file TrackTypeMonitor.h.

Referenced by bookHistograms(), fillHistograms(), and TrackTypeMonitor().

const std::string TrackTypeMonitor::trackQuality_

private

Definition at line 58 of file TrackTypeMonitor.h.

Referenced by analyze().

const edm::InputTag TrackTypeMonitor::trackTag_

private

Definition at line 48 of file TrackTypeMonitor.h.

const edm::EDGetTokenT<reco::TrackCollection> TrackTypeMonitor::trackToken_

private

Definition at line 54 of file TrackTypeMonitor.h.

Referenced by analyze().

bool TrackTypeMonitor::verbose_

private

Definition at line 44 of file TrackTypeMonitor.h.

Referenced by analyze().

const edm::InputTag TrackTypeMonitor::vertexTag_

private

Definition at line 50 of file TrackTypeMonitor.h.

const edm::EDGetTokenT<reco::VertexCollection> TrackTypeMonitor::vertexToken_

private

Definition at line 56 of file TrackTypeMonitor.h.

Referenced by analyze().

Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation