#include <L1TPhase2OuterTrackerTkMET.h>

Inheritance diagram for L1TPhase2OuterTrackerTkMET:

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

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

	L1TPhase2OuterTrackerTkMET (const edm::ParameterSet &)

	~L1TPhase2OuterTrackerTkMET () override

Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) final

void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final

void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final

void	beginStream (edm::StreamID id) final

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

	DQMEDAnalyzer ()

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

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

virtual bool	getCanSaveByLumi ()

Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
	EDProducer ()=default

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

Public Attributes
MonitorElement *	TkMET_PU = nullptr

MonitorElement *	TkMET_QualityCuts = nullptr

MonitorElement *	Track_BendChi2 = nullptr

MonitorElement *	Track_Chi2Dof = nullptr

MonitorElement *	Track_Eta = nullptr

MonitorElement *	Track_N = nullptr

MonitorElement *	Track_NPSstubs = nullptr

MonitorElement *	Track_NStubs = nullptr

MonitorElement *	Track_Pt = nullptr

MonitorElement *	Track_VtxZ = nullptr

Private Attributes
float	bendchi2Max

float	chi2dofMax

edm::ParameterSet	conf_

float	DeltaZ

int	HighPtTracks

const edm::ESGetToken < TrackerTopology, TrackerTopologyRcd >	m_topoToken

float	maxEta

float	maxPt

float	maxZ0

float	minPt

int	nStubsmin

int	nStubsPSmin

edm::EDGetTokenT < l1t::TkPrimaryVertexCollection >	pvToken

std::string	topFolderName_

edm::EDGetTokenT< std::vector < TTTrack < Ref_Phase2TrackerDigi_ > > >	ttTrackToken_

Additional Inherited Members
Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore

typedef dqm::reco::MonitorElement	MonitorElement

Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
using	CacheTypes = CacheContexts< T...>

using	GlobalCache = typename CacheTypes::GlobalCache

using	HasAbility = AbilityChecker< T...>

using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache

using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache

using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >

using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache

using	RunCache = typename CacheTypes::RunCache

using	RunContext = RunContextT< RunCache, GlobalCache >

using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Static Public Member Functions inherited from DQMEDAnalyzer
static void	globalEndJob (DQMEDAnalyzerGlobalCache const *)

static void	globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)

static void	globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)

static std::unique_ptr < DQMEDAnalyzerGlobalCache >	initializeGlobalCache (edm::ParameterSet const &)

Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

unsigned int	streamId_

Detailed Description

Definition at line 22 of file L1TPhase2OuterTrackerTkMET.h.

Constructor & Destructor Documentation

L1TPhase2OuterTrackerTkMET::L1TPhase2OuterTrackerTkMET ( const edm::ParameterSet & iConfig )

explicit

Definition at line 39 of file L1TPhase2OuterTrackerTkMET.cc.

References bendchi2Max, chi2dofMax, conf_, DeltaZ, edm::ParameterSet::getParameter(), HighPtTracks, maxEta, maxPt, maxZ0, minPt, nStubsmin, nStubsPSmin, pvToken, AlCaHLTBitMon_QueryRunRegistry::string, topFolderName_, and ttTrackToken_.

     : conf_(iConfig), m_topoToken(esConsumes()) {
   topFolderName_ = conf_.getParameter<std::string>("TopFolderName");
   ttTrackToken_ =
       consumes<std::vector<TTTrack<Ref_Phase2TrackerDigi_> > >(conf_.getParameter<edm::InputTag>("TTTracksTag"));
   pvToken = consumes<l1t::TkPrimaryVertexCollection>(conf_.getParameter<edm::InputTag>("L1VertexInputTag"));
 
   maxZ0 = conf_.getParameter<double>("maxZ0");
   DeltaZ = conf_.getParameter<double>("DeltaZ");
   chi2dofMax = conf_.getParameter<double>("chi2dofMax");
   bendchi2Max = conf_.getParameter<double>("bendchi2Max");
   minPt = conf_.getParameter<double>("minPt");
   nStubsmin = iConfig.getParameter<int>("nStubsmin");
   nStubsPSmin = iConfig.getParameter<int>("nStubsPSmin");
   maxPt = conf_.getParameter<double>("maxPt");
   maxEta = conf_.getParameter<double>("maxEta");
   HighPtTracks = iConfig.getParameter<int>("HighPtTracks");
 }

L1TPhase2OuterTrackerTkMET::~L1TPhase2OuterTrackerTkMET ( )

override

Definition at line 58 of file L1TPhase2OuterTrackerTkMET.cc.

                                                         {
   // do anything here that needs to be done at desctruction time
   // (e.g. close files, deallocate resources etc.)
 }

Member Function Documentation

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

overridevirtual

Make the pointer

Reimplemented from DQMEDAnalyzer.

Definition at line 66 of file L1TPhase2OuterTrackerTkMET.cc.

References bendchi2Max, chi2dofMax, funct::cos(), DeltaZ, PVValHelper::eta, dqm::impl::MonitorElement::Fill(), edm::Event::getByToken(), edm::EventSetup::getData(), HighPtTracks, edm::HandleBase::isValid(), m_topoToken, maxEta, maxPt, maxZ0, minPt, nStubsmin, nStubsPSmin, phi, DiDispStaMuonMonitor_cfi::pt, pvToken, funct::sin(), mathSSE::sqrt(), StripSubdetector::TID, TrackerTopology::tidRing(), TkMET_PU, TkMET_QualityCuts, StripSubdetector::TOB, TrackerTopology::tobLayer(), Track_BendChi2, Track_Chi2Dof, Track_Eta, Track_N, Track_NPSstubs, Track_NStubs, Track_Pt, Track_VtxZ, align::Tracker, and ttTrackToken_.

                                                                                             {
   // L1 Primaries
   edm::Handle<l1t::TkPrimaryVertexCollection> L1VertexHandle;
   iEvent.getByToken(pvToken, L1VertexHandle);
 
   edm::Handle<std::vector<TTTrack<Ref_Phase2TrackerDigi_> > > TTTrackHandle;
   iEvent.getByToken(ttTrackToken_, TTTrackHandle);
 
   // for PS stubs
   // Tracker Topology
   const TrackerTopology* const tTopo = &iSetup.getData(m_topoToken);
 
   // Adding protection
   if (!TTTrackHandle.isValid()) {
     edm::LogWarning("L1TPhase2OuterTrackerTkMET") << "cant find tracks" << std::endl;
     return;
   }
   if (!L1VertexHandle.isValid()) {
     edm::LogWarning("L1TPhase2OuterTrackerTkMET") << "cant find vertex" << std::endl;
     return;
   }
   float sumPx = 0;
   float sumPy = 0;
   float etTot = 0;
   double sumPx_PU = 0;
   double sumPy_PU = 0;
   double etTot_PU = 0;
   int nTracks_counter = 0;
 
   float zVTX = L1VertexHandle->begin()->zvertex();
   unsigned int tkCnt = 0;
   for (const auto& trackIter : *TTTrackHandle) {
     edm::Ptr<TTTrack<Ref_Phase2TrackerDigi_> > tempTrackPtr(TTTrackHandle, tkCnt++);  
     float pt = tempTrackPtr->momentum().perp();
     float phi = tempTrackPtr->momentum().phi();
     float eta = tempTrackPtr->momentum().eta();
     std::vector<edm::Ref<edmNew::DetSetVector<TTStub<Ref_Phase2TrackerDigi_> >, TTStub<Ref_Phase2TrackerDigi_> > >
         theStubs = trackIter.getStubRefs();
     int nstubs = (int)theStubs.size();
 
     float chi2dof = tempTrackPtr->chi2Red();
     float bendchi2 = tempTrackPtr->stubPtConsistency();
     float z0 = tempTrackPtr->z0();
 
     if (pt < minPt || fabs(z0) > maxZ0 || fabs(eta) > maxEta || chi2dof > chi2dofMax || bendchi2 > bendchi2Max)
       continue;
     if (maxPt > 0 && pt > maxPt) {
       if (HighPtTracks == 0)
         continue;  // ignore these very high PT tracks: truncate
       if (HighPtTracks == 1)
         pt = maxPt;  // saturate
     }
 
     int nPS = 0.;  // number of stubs in PS modules
     // loop over the stubs
     for (unsigned int istub = 0; istub < (unsigned int)theStubs.size(); istub++) {
       DetId detId(theStubs.at(istub)->getDetId());
       if (detId.det() == DetId::Detector::Tracker) {
         if ((detId.subdetId() == StripSubdetector::TOB && tTopo->tobLayer(detId) <= 3) ||
             (detId.subdetId() == StripSubdetector::TID && tTopo->tidRing(detId) <= 9))
           nPS++;
       }
     }
 
     if (nstubs < nStubsmin || nPS < nStubsPSmin)
       continue;
 
     // construct deltaZ cut to be based on track eta
     if (fabs(eta) >= 0 && fabs(eta) < 0.7)
       DeltaZ = 0.4;
     else if (fabs(eta) >= 0.7 && fabs(eta) < 1.0)
       DeltaZ = 0.6;
     else if (fabs(eta) >= 1.0 && fabs(eta) < 1.2)
       DeltaZ = 0.76;
     else if (fabs(eta) >= 1.2 && fabs(eta) < 1.6)
       DeltaZ = 1.0;
     else if (fabs(eta) >= 1.6 && fabs(eta) < 2.0)
       DeltaZ = 1.7;
     else if (fabs(eta) >= 2.0 && fabs(eta) <= 2.4)
       DeltaZ = 2.2;
 
     if (fabs(z0 - zVTX) <= DeltaZ) {
       nTracks_counter++;
       Track_Pt->Fill(pt);
       Track_NStubs->Fill(nstubs);
       Track_NPSstubs->Fill(nPS);
       Track_Eta->Fill(eta);
       Track_VtxZ->Fill(z0);
       Track_Chi2Dof->Fill(chi2dof);
       Track_BendChi2->Fill(bendchi2);
 
       sumPx += pt * cos(phi);
       sumPy += pt * sin(phi);
       etTot += pt;
     } else {  // PU sums
       sumPx_PU += pt * cos(phi);
       sumPy_PU += pt * sin(phi);
       etTot_PU += pt;
     }
   }  // end loop over tracks
 
   Track_N->Fill(nTracks_counter);
   float et = sqrt(sumPx * sumPx + sumPy * sumPy);
   double etmiss_PU = sqrt(sumPx_PU * sumPx_PU + sumPy_PU * sumPy_PU);
 
   math::XYZTLorentzVector missingEt(-sumPx, -sumPy, 0, et);
 
   TkMET_QualityCuts->Fill(missingEt.Pt());
   TkMET_PU->Fill(etmiss_PU);
 
 }  // end of method

void L1TPhase2OuterTrackerTkMET::bookHistograms	(	DQMStore::IBooker &	iBooker,
		edm::Run const &	run,
		edm::EventSetup const &	es
	)

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 180 of file L1TPhase2OuterTrackerTkMET.cc.

References dqm::implementation::IBooker::book1D(), conf_, edm::ParameterSet::getParameter(), dqm::impl::MonitorElement::setAxisTitle(), dqm::implementation::NavigatorBase::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, TkMET_PU, TkMET_QualityCuts, topFolderName_, Track_BendChi2, Track_Chi2Dof, Track_Eta, Track_N, Track_NPSstubs, Track_NStubs, Track_Pt, and Track_VtxZ.

                                                                          {
   std::string HistoName;
   iBooker.setCurrentFolder(topFolderName_ + "/TkMET_Tracks/");
 
   // Num of L1Tracks in tkMET selection
   HistoName = "Track_N";
   edm::ParameterSet psTrack_N = conf_.getParameter<edm::ParameterSet>("TH1_NTracks");
   Track_N = iBooker.book1D(HistoName,
                            HistoName,
                            psTrack_N.getParameter<int32_t>("Nbinsx"),
                            psTrack_N.getParameter<double>("xmin"),
                            psTrack_N.getParameter<double>("xmax"));
   Track_N->setAxisTitle("# L1 Tracks", 1);
   Track_N->setAxisTitle("# Events", 2);
 
   //Pt of the tracks
   edm::ParameterSet psTrack_Pt = conf_.getParameter<edm::ParameterSet>("TH1_Track_Pt");
   HistoName = "Track_Pt";
   Track_Pt = iBooker.book1D(HistoName,
                             HistoName,
                             psTrack_Pt.getParameter<int32_t>("Nbinsx"),
                             psTrack_Pt.getParameter<double>("xmin"),
                             psTrack_Pt.getParameter<double>("xmax"));
   Track_Pt->setAxisTitle("p_{T} [GeV]", 1);
   Track_Pt->setAxisTitle("# L1 Tracks", 2);
 
   //Eta
   edm::ParameterSet psTrack_Eta = conf_.getParameter<edm::ParameterSet>("TH1_Track_Eta");
   HistoName = "Track_Eta";
   Track_Eta = iBooker.book1D(HistoName,
                              HistoName,
                              psTrack_Eta.getParameter<int32_t>("Nbinsx"),
                              psTrack_Eta.getParameter<double>("xmin"),
                              psTrack_Eta.getParameter<double>("xmax"));
   Track_Eta->setAxisTitle("#eta", 1);
   Track_Eta->setAxisTitle("# L1 Tracks", 2);
 
   //VtxZ
   edm::ParameterSet psTrack_VtxZ = conf_.getParameter<edm::ParameterSet>("TH1_Track_VtxZ");
   HistoName = "Track_VtxZ";
   Track_VtxZ = iBooker.book1D(HistoName,
                               HistoName,
                               psTrack_VtxZ.getParameter<int32_t>("Nbinsx"),
                               psTrack_VtxZ.getParameter<double>("xmin"),
                               psTrack_VtxZ.getParameter<double>("xmax"));
   Track_VtxZ->setAxisTitle("L1 Track vertex position z [cm]", 1);
   Track_VtxZ->setAxisTitle("# L1 Tracks", 2);
 
   //chi2dof
   edm::ParameterSet psTrack_Chi2Dof = conf_.getParameter<edm::ParameterSet>("TH1_Track_Chi2Dof");
   HistoName = "Track_Chi2Dof";
   Track_Chi2Dof = iBooker.book1D(HistoName,
                                  HistoName,
                                  psTrack_Chi2Dof.getParameter<int32_t>("Nbinsx"),
                                  psTrack_Chi2Dof.getParameter<double>("xmin"),
                                  psTrack_Chi2Dof.getParameter<double>("xmax"));
   Track_Chi2Dof->setAxisTitle("L1 Track #chi^{2}/D.O.F.", 1);
   Track_Chi2Dof->setAxisTitle("# L1 Tracks", 2);
 
   //bend chi2
   edm::ParameterSet psTrack_BendChi2 = conf_.getParameter<edm::ParameterSet>("TH1_Track_BendChi2");
   HistoName = "Track_BendChi2";
   Track_BendChi2 = iBooker.book1D(HistoName,
                                   HistoName,
                                   psTrack_BendChi2.getParameter<int32_t>("Nbinsx"),
                                   psTrack_BendChi2.getParameter<double>("xmin"),
                                   psTrack_BendChi2.getParameter<double>("xmax"));
   Track_BendChi2->setAxisTitle("L1 Track Bend #chi^{2}", 1);
   Track_BendChi2->setAxisTitle("# L1 Tracks", 2);
 
   //nstubs
   edm::ParameterSet psTrack_NStubs = conf_.getParameter<edm::ParameterSet>("TH1_Track_NStubs");
   HistoName = "Track_NStubs";
   Track_NStubs = iBooker.book1D(HistoName,
                                 HistoName,
                                 psTrack_NStubs.getParameter<int32_t>("Nbinsx"),
                                 psTrack_NStubs.getParameter<double>("xmin"),
                                 psTrack_NStubs.getParameter<double>("xmax"));
   Track_NStubs->setAxisTitle("# L1 Stubs", 1);
   Track_NStubs->setAxisTitle("# L1 Tracks", 2);
 
   //nPSstubs
   edm::ParameterSet psTrack_NPSstubs = conf_.getParameter<edm::ParameterSet>("TH1_Track_NPSstubs");
   HistoName = "Track_NPSstubs";
   Track_NPSstubs = iBooker.book1D(HistoName,
                                   HistoName,
                                   psTrack_NPSstubs.getParameter<int32_t>("Nbinsx"),
                                   psTrack_NPSstubs.getParameter<double>("xmin"),
                                   psTrack_NPSstubs.getParameter<double>("xmax"));
   Track_NPSstubs->setAxisTitle("# PS Stubs", 1);
   Track_NPSstubs->setAxisTitle("# L1 Tracks", 2);
 
   iBooker.setCurrentFolder(topFolderName_ + "/TkMET");
   //loose tkMET
   edm::ParameterSet psTrack_TkMET = conf_.getParameter<edm::ParameterSet>("TH1_Track_TkMET");
   HistoName = "TkMET_QualityCuts";
   TkMET_QualityCuts = iBooker.book1D(HistoName,
                                      HistoName,
                                      psTrack_TkMET.getParameter<int32_t>("Nbinsx"),
                                      psTrack_TkMET.getParameter<double>("xmin"),
                                      psTrack_TkMET.getParameter<double>("xmax"));
   TkMET_QualityCuts->setAxisTitle("L1 Track MET [GeV]", 1);
   TkMET_QualityCuts->setAxisTitle("# Events", 2);
 
   //tkMET -- PU only
   HistoName = "TkMET_PU";
   TkMET_PU = iBooker.book1D(HistoName,
                             HistoName,
                             psTrack_TkMET.getParameter<int32_t>("Nbinsx"),
                             psTrack_TkMET.getParameter<double>("xmin"),
                             psTrack_TkMET.getParameter<double>("xmax"));
   TkMET_PU->setAxisTitle("L1 Track MET (PU only) [GeV]", 1);
   TkMET_PU->setAxisTitle("# Events", 2);
 }  //end of method