L1TrackerEtMissProducer Class Reference

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

L1TTTrackCollectionType

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

Definition at line 33 of file L1TrackerEtMissProducer.cc.

L1TTTrackType

typedef TTTrack<Ref_Phase2TrackerDigi_> L1TrackerEtMissProducer::L1TTTrackType

Definition at line 32 of file L1TrackerEtMissProducer.cc.

Constructor & Destructor Documentation

L1TrackerEtMissProducer()

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

explicit

Definition at line 62 of file L1TrackerEtMissProducer.cc.

    : 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");
}

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

~L1TrackerEtMissProducer()

L1TrackerEtMissProducer::~L1TrackerEtMissProducer ( )

override

Definition at line 85 of file L1TrackerEtMissProducer.cc.

{}

Member Function Documentation

produce()

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

overrideprivate

Definition at line 87 of file L1TrackerEtMissProducer.cc.

                                                                                   {
  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

References bendChi2Max_, chi2dofMax_, funct::cos(), deltaZ_, displaced_, EgHLTOffHistBins_cfi::et, PVValHelper::eta, edm::EventSetup::getData(), highPtTracks_, iEvent, createfilelist::int, edm::HandleBase::isValid(), 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_, and HLTMuonOfflineAnalyzer_cfi::z0.

Member Data Documentation

bendChi2Max_

float L1TrackerEtMissProducer::bendChi2Max_

private

Definition at line 48 of file L1TrackerEtMissProducer.cc.

Referenced by L1TrackerEtMissProducer(), and produce().

chi2dofMax_

float L1TrackerEtMissProducer::chi2dofMax_

private

Definition at line 47 of file L1TrackerEtMissProducer.cc.

Referenced by L1TrackerEtMissProducer(), and produce().

deltaZ_

float L1TrackerEtMissProducer::deltaZ_

private

Definition at line 45 of file L1TrackerEtMissProducer.cc.

Referenced by L1TrackerEtMissProducer(), and produce().

displaced_

bool L1TrackerEtMissProducer::displaced_

private

Definition at line 54 of file L1TrackerEtMissProducer.cc.

Referenced by L1TrackerEtMissProducer(), and produce().

highPtTracks_

int L1TrackerEtMissProducer::highPtTracks_

private

Definition at line 53 of file L1TrackerEtMissProducer.cc.

Referenced by L1TrackerEtMissProducer(), and produce().

maxEta_

float L1TrackerEtMissProducer::maxEta_

private

Definition at line 46 of file L1TrackerEtMissProducer.cc.

Referenced by L1TrackerEtMissProducer(), and produce().

maxPt_

float L1TrackerEtMissProducer::maxPt_

private

Definition at line 52 of file L1TrackerEtMissProducer.cc.

Referenced by L1TrackerEtMissProducer(), and produce().

maxZ0_

float L1TrackerEtMissProducer::maxZ0_

private

Definition at line 44 of file L1TrackerEtMissProducer.cc.

Referenced by L1TrackerEtMissProducer(), and produce().

minPt_

float L1TrackerEtMissProducer::minPt_

private

Definition at line 49 of file L1TrackerEtMissProducer.cc.

Referenced by L1TrackerEtMissProducer(), and produce().

nPSStubsMin_

int L1TrackerEtMissProducer::nPSStubsMin_

private

Definition at line 51 of file L1TrackerEtMissProducer.cc.

Referenced by L1TrackerEtMissProducer(), and produce().

nStubsmin_

int L1TrackerEtMissProducer::nStubsmin_

private

Definition at line 50 of file L1TrackerEtMissProducer.cc.

Referenced by L1TrackerEtMissProducer(), and produce().

pvToken_

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

private

Definition at line 56 of file L1TrackerEtMissProducer.cc.

Referenced by produce().

trackToken_

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

private

Definition at line 57 of file L1TrackerEtMissProducer.cc.

Referenced by produce().

tTopoToken_

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

private

Definition at line 58 of file L1TrackerEtMissProducer.cc.

Referenced by produce().