Inheritance diagram for L1TrackerEtMissProducer:

Public Types
typedef std::vector< L1TTTrackType >	L1TTTrackCollectionType

typedef TTTrack< Ref_Phase2TrackerDigi_ >	L1TTTrackType

typedef VertexCollection	L1VertexCollectionType

typedef Vertex	L1VertexType

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
virtual void	beginJob ()

virtual void	endJob ()

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

Private Attributes
const float	bendChi2Max_

const float	Chi2RphidofMax_

const float	Chi2RzdofMax_

bool	debug_

float	deltaZ_

const bool	displaced_

vector< double >	etaRegions_

const int	highPtTracks_

std::string	L1ExtendedMetCollectionName

std::string	L1MetCollectionName

const float	maxEta_

const float	maxPt_

const float	maxZ0_

const float	minPt_

const int	nPSStubsMin_

const int	nStubsmin_

const edm::EDGetTokenT< VertexCollection >	pvToken_

const edm::EDGetTokenT< L1TTTrackCollectionType >	trackToken_

edm::ESGetToken< TrackerTopology, TrackerTopologyRcd >	tTopoToken_

vector< double >	z0Thresholds_

Detailed Description

Definition at line 33 of file L1TrackerEtMissProducer.cc.

Member Typedef Documentation

◆ L1TTTrackCollectionType

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

Definition at line 36 of file L1TrackerEtMissProducer.cc.

◆ L1TTTrackType

typedef TTTrack<Ref_Phase2TrackerDigi_> L1TrackerEtMissProducer::L1TTTrackType

Definition at line 35 of file L1TrackerEtMissProducer.cc.

◆ L1VertexCollectionType

typedef VertexCollection L1TrackerEtMissProducer::L1VertexCollectionType

Definition at line 39 of file L1TrackerEtMissProducer.cc.

◆ L1VertexType

typedef Vertex L1TrackerEtMissProducer::L1VertexType

Definition at line 38 of file L1TrackerEtMissProducer.cc.

Constructor & Destructor Documentation

◆ L1TrackerEtMissProducer()

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

explicit

Definition at line 75 of file L1TrackerEtMissProducer.cc.

References displaced_, edm::ParameterSet::getParameter(), L1ExtendedMetCollectionName, L1MetCollectionName, and AlCaHLTBitMon_QueryRunRegistry::string.

     : maxZ0_((float)iConfig.getParameter<double>("maxZ0")),
       deltaZ_((float)iConfig.getParameter<double>("deltaZ")),
       Chi2RphidofMax_((float)iConfig.getParameter<double>("chi2rphidofMax")),
       Chi2RzdofMax_((float)iConfig.getParameter<double>("chi2rzdofMax")),
       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")),
       z0Thresholds_(iConfig.getParameter<std::vector<double>>("z0Thresholds")),
       etaRegions_(iConfig.getParameter<std::vector<double>>("etaRegions")),
       debug_(iConfig.getParameter<bool>("debug")),
       pvToken_(consumes<L1VertexCollectionType>(iConfig.getParameter<edm::InputTag>("L1VertexInputTag"))),
       trackToken_(consumes<L1TTTrackCollectionType>(iConfig.getParameter<edm::InputTag>("L1TrackInputTag"))),
       tTopoToken_(esConsumes<TrackerTopology, TrackerTopologyRcd>(edm::ESInputTag("", ""))) {
   L1MetCollectionName = (std::string)iConfig.getParameter<std::string>("L1MetCollectionName");
 
   if (displaced_) {
     L1ExtendedMetCollectionName = (std::string)iConfig.getParameter<std::string>("L1MetExtendedCollectionName");
     produces<TkEtMissCollection>(L1ExtendedMetCollectionName);
   } else
     produces<TkEtMissCollection>(L1MetCollectionName);
 }

◆ ~L1TrackerEtMissProducer()

L1TrackerEtMissProducer::~L1TrackerEtMissProducer ( )

override

Definition at line 103 of file L1TrackerEtMissProducer.cc.

103 {}

Member Function Documentation

◆ beginJob()

void L1TrackerEtMissProducer::beginJob ( void )

privatevirtual

Definition at line 249 of file L1TrackerEtMissProducer.cc.

249 {}

◆ endJob()

void L1TrackerEtMissProducer::endJob ( void )

privatevirtual

Definition at line 251 of file L1TrackerEtMissProducer.cc.

251 {}

◆ produce()

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

overrideprivate

Definition at line 105 of file L1TrackerEtMissProducer.cc.

References funct::abs(), bendChi2Max_, Chi2RphidofMax_, Chi2RzdofMax_, funct::cos(), debug_, deltaZ_, displaced_, EgHLTOffHistBins_cfi::et, PVValHelper::eta, etaRegions_, edm::EventSetup::getData(), highPtTracks_, iEvent, createfilelist::int, edm::HandleBase::isValid(), L1ExtendedMetCollectionName, L1MetCollectionName, maxEta_, maxPt_, maxZ0_, minPt_, eostools::move(), nPSStubsMin_, nStubsmin_, phi, DiDispStaMuonMonitor_cfi::pt, pvToken_, funct::sin(), mathSSE::sqrt(), StripSubdetector::TID, TrackerTopology::tidRing(), StripSubdetector::TOB, TrackerTopology::tobLayer(), align::Tracker, trackToken_, tTopoToken_, HLTMuonOfflineAnalyzer_cfi::z0, and z0Thresholds_.

                                                                                    {
   using namespace edm;
 
   std::unique_ptr<TkEtMissCollection> METCollection(new TkEtMissCollection);
 
   // Tracker Topology
   const TrackerTopology& tTopo = iSetup.getData(tTopoToken_);
 
   edm::Handle<L1VertexCollectionType> L1VertexHandle;
   iEvent.getByToken(pvToken_, L1VertexHandle);
 
   edm::Handle<L1TTTrackCollectionType> L1TTTrackHandle;
   iEvent.getByToken(trackToken_, L1TTTrackHandle);
   L1TTTrackCollectionType::const_iterator trackIter;
 
   if (!L1VertexHandle.isValid()) {
     LogError("L1TrackerEtMissProducer") << "\nWarning: VertexCollection 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()->z0();
 
   int numtracks = 0;
   int numqualitytracks = 0;
   int numassoctracks = 0;
 
   for (trackIter = L1TTTrackHandle->begin(); trackIter != L1TTTrackHandle->end(); ++trackIter) {
     numtracks++;
     float pt = trackIter->momentum().perp();
     float phi = trackIter->momentum().phi();
     float eta = trackIter->momentum().eta();
     float chi2rphidof = trackIter->chi2XYRed();
     float chi2rzdof = trackIter->chi2ZRed();
     float bendChi2 = trackIter->stubPtConsistency();
     float z0 = trackIter->z0();
     //unsigned int Sector = trackIter->phiSector();
     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 (std::abs(z0) > maxZ0_)
       continue;
     if (std::abs(eta) > maxEta_)
       continue;
     if (chi2rphidof > Chi2RphidofMax_)
       continue;
     if (chi2rzdof > Chi2RzdofMax_)
       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;
 
     numqualitytracks++;
 
     if (!displaced_) {  // if displaced, deltaZ = 3.0 cm, very loose
       // construct deltaZ cut to be based on track eta
       for (unsigned int reg = 0; reg < etaRegions_.size() - 1; reg++) {
         if (std::abs(eta) >= etaRegions_[reg] && std::abs(eta) < etaRegions_[reg + 1]) {
           deltaZ_ = z0Thresholds_[reg];
           break;
         }
       }
       if (std::abs(eta) >= etaRegions_[etaRegions_.size() - 1]) {
         deltaZ_ = z0Thresholds_[etaRegions_.size() - 1];
         break;
       }
     }
 
     if (std::abs(z0 - zVTX) <= deltaZ_) {
       numassoctracks++;
       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 etphi = atan2(sumPy, sumPx);
 
   math::XYZTLorentzVector missingEt(-sumPx, -sumPy, 0, et);
 
   if (debug_) {
     edm::LogVerbatim("L1TrackerEtMissProducer") << "====Global Pt===="
                                                 << "\n"
                                                 << "Px: " << sumPx << "| Py: " << sumPy << "\n"
                                                 << "====MET==="
                                                 << "\n"
                                                 << "MET: " << et << "| Phi: " << etphi << "\n"
 
                                                 << "# Intial Tracks: " << numtracks << "\n"
                                                 << "# Tracks after Quality Cuts: " << numqualitytracks << "\n"
                                                 << "# Tracks Associated to Vertex: " << numassoctracks << "\n"
                                                 << "========================================================"
                                                 << "\n";
   }
 
   int ibx = 0;
   METCollection->push_back(TkEtMiss(missingEt, TkEtMiss::kMET, etphi, numassoctracks, ibx));
 
   if (displaced_)
     iEvent.put(std::move(METCollection), L1ExtendedMetCollectionName);
   else
     iEvent.put(std::move(METCollection), L1MetCollectionName);
 }  // end producer

Member Data Documentation

◆ bendChi2Max_

const float L1TrackerEtMissProducer::bendChi2Max_

private

Definition at line 54 of file L1TrackerEtMissProducer.cc.

Referenced by produce().

◆ Chi2RphidofMax_

const float L1TrackerEtMissProducer::Chi2RphidofMax_

private

Definition at line 52 of file L1TrackerEtMissProducer.cc.

Referenced by produce().

◆ Chi2RzdofMax_

const float L1TrackerEtMissProducer::Chi2RzdofMax_

private

Definition at line 53 of file L1TrackerEtMissProducer.cc.

Referenced by produce().

◆ debug_

bool L1TrackerEtMissProducer::debug_

private

Definition at line 64 of file L1TrackerEtMissProducer.cc.

Referenced by produce().

◆ deltaZ_

float L1TrackerEtMissProducer::deltaZ_

private

Definition at line 51 of file L1TrackerEtMissProducer.cc.

Referenced by produce().

◆ displaced_

const bool L1TrackerEtMissProducer::displaced_

private

Definition at line 61 of file L1TrackerEtMissProducer.cc.

Referenced by L1TrackerEtMissProducer(), and produce().

◆ etaRegions_

vector<double> L1TrackerEtMissProducer::etaRegions_

private

Definition at line 63 of file L1TrackerEtMissProducer.cc.

Referenced by produce().

◆ highPtTracks_

const int L1TrackerEtMissProducer::highPtTracks_

private

Definition at line 60 of file L1TrackerEtMissProducer.cc.

Referenced by produce().

◆ L1ExtendedMetCollectionName

std::string L1TrackerEtMissProducer::L1ExtendedMetCollectionName

private

Definition at line 67 of file L1TrackerEtMissProducer.cc.

Referenced by L1TrackerEtMissProducer(), and produce().

◆ L1MetCollectionName

std::string L1TrackerEtMissProducer::L1MetCollectionName

private

Definition at line 66 of file L1TrackerEtMissProducer.cc.

Referenced by L1TrackerEtMissProducer(), and produce().

◆ maxEta_

const float L1TrackerEtMissProducer::maxEta_

private

Definition at line 59 of file L1TrackerEtMissProducer.cc.

Referenced by produce().

◆ maxPt_

const float L1TrackerEtMissProducer::maxPt_

private

Definition at line 58 of file L1TrackerEtMissProducer.cc.

Referenced by produce().

◆ maxZ0_

const float L1TrackerEtMissProducer::maxZ0_

private

Definition at line 50 of file L1TrackerEtMissProducer.cc.

Referenced by produce().

◆ minPt_

const float L1TrackerEtMissProducer::minPt_

private

Definition at line 55 of file L1TrackerEtMissProducer.cc.

Referenced by produce().

◆ nPSStubsMin_

const int L1TrackerEtMissProducer::nPSStubsMin_

private

Definition at line 57 of file L1TrackerEtMissProducer.cc.

Referenced by produce().

◆ nStubsmin_

const int L1TrackerEtMissProducer::nStubsmin_

private

Definition at line 56 of file L1TrackerEtMissProducer.cc.

Referenced by produce().

◆ pvToken_

const edm::EDGetTokenT<VertexCollection> L1TrackerEtMissProducer::pvToken_

private

Definition at line 69 of file L1TrackerEtMissProducer.cc.

Referenced by produce().

◆ trackToken_

const edm::EDGetTokenT<L1TTTrackCollectionType> L1TrackerEtMissProducer::trackToken_

private

Definition at line 70 of file L1TrackerEtMissProducer.cc.

Referenced by produce().

◆ tTopoToken_

edm::ESGetToken<TrackerTopology, TrackerTopologyRcd> L1TrackerEtMissProducer::tTopoToken_

private

Definition at line 71 of file L1TrackerEtMissProducer.cc.

Referenced by produce().

◆ z0Thresholds_

vector<double> L1TrackerEtMissProducer::z0Thresholds_

private

Definition at line 62 of file L1TrackerEtMissProducer.cc.

Referenced by produce().

Public Types

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ L1TTTrackCollectionType

◆ L1TTTrackType

◆ L1VertexCollectionType

◆ L1VertexType

Constructor & Destructor Documentation

◆ L1TrackerEtMissProducer()

◆ ~L1TrackerEtMissProducer()

Member Function Documentation

◆ beginJob()

◆ endJob()

◆ produce()

Member Data Documentation

◆ bendChi2Max_

◆ Chi2RphidofMax_

◆ Chi2RzdofMax_

◆ debug_

◆ deltaZ_

◆ displaced_

◆ etaRegions_

◆ highPtTracks_

◆ L1ExtendedMetCollectionName

◆ L1MetCollectionName

◆ maxEta_

◆ maxPt_

◆ maxZ0_

◆ minPt_

◆ nPSStubsMin_

◆ nStubsmin_

◆ pvToken_

◆ trackToken_

◆ tTopoToken_

◆ z0Thresholds_