Inheritance diagram for L1TrackerEtMissProducer:

Public Types
typedef std::vector < L1TTTrackType >	L1TTTrackCollectionType

typedef TTTrack < Ref_Phase2TrackerDigi_ >	L1TTTrackType

Public Types inherited from edm::stream::EDProducer<>
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

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

	~L1TrackerEtMissProducer () override

Public Member Functions inherited from edm::stream::EDProducer<>
	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

Private Member Functions
void	produce (edm::Event &, const edm::EventSetup &) override

Private Attributes
float	bendChi2Max_

float	chi2dofMax_

float	deltaZ_

bool	displaced_

int	highPtTracks_

float	maxEta_

float	maxPt_

float	maxZ0_

float	minPt_

int	nPSStubsMin_

int	nStubsmin_

const edm::EDGetTokenT < TkPrimaryVertexCollection >	pvToken_

const edm::EDGetTokenT < std::vector< TTTrack < Ref_Phase2TrackerDigi_ > > >	trackToken_

edm::ESGetToken < TrackerTopology, TrackerTopologyRcd >	tTopoToken_

Detailed Description

Definition at line 30 of file L1TrackerEtMissProducer.cc.

Member Typedef Documentation

typedef std::vector<L1TTTrackType> L1TrackerEtMissProducer::L1TTTrackCollectionType

Definition at line 33 of file L1TrackerEtMissProducer.cc.

typedef TTTrack<Ref_Phase2TrackerDigi_> L1TrackerEtMissProducer::L1TTTrackType

Definition at line 32 of file L1TrackerEtMissProducer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 62 of file L1TrackerEtMissProducer.cc.

References bendChi2Max_, chi2dofMax_, deltaZ_, displaced_, edm::ParameterSet::getParameter(), highPtTracks_, maxEta_, maxPt_, maxZ0_, minPt_, nPSStubsMin_, and nStubsmin_.

     : pvToken_(consumes<TkPrimaryVertexCollection>(iConfig.getParameter<edm::InputTag>("L1VertexInputTag"))),
       trackToken_(consumes<std::vector<TTTrack<Ref_Phase2TrackerDigi_> > >(
           iConfig.getParameter<edm::InputTag>("L1TrackInputTag"))),
       tTopoToken_(esConsumes<TrackerTopology, TrackerTopologyRcd>(edm::ESInputTag("", ""))) {
   maxZ0_ = (float)iConfig.getParameter<double>("maxZ0");
   deltaZ_ = (float)iConfig.getParameter<double>("deltaZ");
   chi2dofMax_ = (float)iConfig.getParameter<double>("chi2dofMax");
   bendChi2Max_ = (float)iConfig.getParameter<double>("bendChi2Max");
   minPt_ = (float)iConfig.getParameter<double>("minPt");
   nStubsmin_ = iConfig.getParameter<int>("nStubsmin");
   nPSStubsMin_ = iConfig.getParameter<int>("nPSStubsMin");
   maxPt_ = (float)iConfig.getParameter<double>("maxPt");
   maxEta_ = (float)iConfig.getParameter<double>("maxEta");
   highPtTracks_ = iConfig.getParameter<int>("highPtTracks");
   displaced_ = iConfig.getParameter<bool>("displaced");
 
   if (displaced_)
     produces<TkEtMissCollection>("L1TrackerEtMissExtended");
   else
     produces<TkEtMissCollection>("L1TrackerEtMiss");
 }

L1TrackerEtMissProducer::~L1TrackerEtMissProducer ( )

override

Definition at line 85 of file L1TrackerEtMissProducer.cc.

85 {}

Member Function Documentation

void L1TrackerEtMissProducer::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overrideprivate

Definition at line 87 of file L1TrackerEtMissProducer.cc.

References bendChi2Max_, chi2dofMax_, funct::cos(), deltaZ_, displaced_, PVValHelper::eta, edm::Event::getByToken(), edm::EventSetup::getData(), highPtTracks_, maxEta_, maxPt_, maxZ0_, minPt_, eostools::move(), nPSStubsMin_, nStubsmin_, phi, DiDispStaMuonMonitor_cfi::pt, edm::Event::put(), pvToken_, funct::sin(), mathSSE::sqrt(), StripSubdetector::TID, StripSubdetector::TOB, align::Tracker, trackToken_, and tTopoToken_.

                                                                                    {
   using namespace edm;
 
   std::unique_ptr<TkEtMissCollection> METCollection(new TkEtMissCollection);
 
   // Tracker Topology
   const TrackerTopology& tTopo = iSetup.getData(tTopoToken_);
 
   edm::Handle<TkPrimaryVertexCollection> L1VertexHandle;
   iEvent.getByToken(pvToken_, L1VertexHandle);
 
   edm::Handle<L1TTTrackCollectionType> L1TTTrackHandle;
   iEvent.getByToken(trackToken_, L1TTTrackHandle);
   L1TTTrackCollectionType::const_iterator trackIter;
 
   if (!L1VertexHandle.isValid()) {
     LogError("L1TrackerEtMissProducer") << "\nWarning: TkPrimaryVertexCollection not found in the event. Exit\n";
     return;
   }
 
   if (!L1TTTrackHandle.isValid()) {
     LogError("L1TrackerEtMissProducer") << "\nWarning: L1TTTrackCollection not found in the event. Exit\n";
     return;
   }
 
   float sumPx = 0;
   float sumPy = 0;
   float etTot = 0;
   double sumPx_PU = 0;
   double sumPy_PU = 0;
   double etTot_PU = 0;
   float zVTX = L1VertexHandle->begin()->zvertex();
 
   for (trackIter = L1TTTrackHandle->begin(); trackIter != L1TTTrackHandle->end(); ++trackIter) {
     float pt = trackIter->momentum().perp();
     float phi = trackIter->momentum().phi();
     float eta = trackIter->momentum().eta();
     float chi2dof = trackIter->chi2Red();
     float bendChi2 = trackIter->stubPtConsistency();
     float z0 = trackIter->z0();
     std::vector<edm::Ref<edmNew::DetSetVector<TTStub<Ref_Phase2TrackerDigi_> >, TTStub<Ref_Phase2TrackerDigi_> > >
         theStubs = trackIter->getStubRefs();
     int nstubs = (int)theStubs.size();
 
     if (pt < minPt_)
       continue;
     if (fabs(z0) > maxZ0_)
       continue;
     if (fabs(eta) > maxEta_)
       continue;
     if (chi2dof > chi2dofMax_)
       continue;
     if (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_)
       continue;
     if (nPS < nPSStubsMin_)
       continue;
 
     if (!displaced_) {  // if displaced, deltaZ = 3.0 cm, very loose
       // 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_) {
       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
 
   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);
   int ibx = 0;
   METCollection->push_back(TkEtMiss(missingEt, TkEtMiss::kMET, etTot, etmiss_PU, etTot_PU, ibx));
 
   if (displaced_)
     iEvent.put(std::move(METCollection), "L1TrackerEtMissExtended");
   else
     iEvent.put(std::move(METCollection), "L1TrackerEtMiss");
 }  // end producer