Definition at line 145 of file L1TrackJetProducer.cc.
References funct::abs(), TTTrack< T >::chi2Red(), MaxZBin::clusters, filterCSVwithJSON::copy, funct::cos(), TTTrack< T >::d0(), d0CutNStubs4_, d0CutNStubs5_, displaced_, dzPVTrk_, EtaPhiBin::eta, PVValHelper::eta, PV3DBase< T, PVType, FrameType >::eta(), egammaIdentification::eta_max, egammaIdentification::eta_min, etaBins_, etaStep_, edm::EventSetup::getData(), TTTrack< T >::getStubRefs(), highpTJetMinTrackMultiplicity_, highpTJetThreshold_, MaxZBin::ht, mps_fire::i, iEvent, createfilelist::int, dqmiolumiharvest::j, reco::btau::jetEta, reco::btau::jetPhi, reco::btau::jetPt, dqmdumpme::k, L1_clustering(), L1TrkPtrs_, L2_clustering(), lowpTJetMinTrackMultiplicity_, lowpTJetThreshold_, M_PI, minTrkJetpT_, TTTrack< T >::momentum(), eostools::move(), MaxZBin::nclust, nDisplacedTracks_, EtaPhiBin::numtdtrks, EtaPhiBin::numtracks, PV3DBase< T, PVType, FrameType >::perp(), phi, EtaPhiBin::phi, phiBins_, phiStep_, EtaPhiBin::pTtot, PVtxToken_, funct::sin(), TTTrack< T >::stubPtConsistency(), tdtrk_, StripSubdetector::TID, TrackerTopology::tidRing(), StripSubdetector::TOB, TrackerTopology::tobLayer(), align::Tracker, EtaPhiBin::trackidx, trackQualityCuts(), trackToken_, trkEtaMax_, trkNPSStubMin_, trkPtMax_, trkPtMin_, trkZMax_, tTopoToken_, HLTMuonOfflineAnalyzer_cfi::z0, TTTrack< T >::z0(), MaxZBin::zbincenter, zBins_, SiStripMonitorCluster_cfi::zmax, SiStripMonitorCluster_cfi::zmin, MaxZBin::znum, and zStep_.
156 float PVz = (PVtx->at(0)).
z0();
161 for (
unsigned int this_l1track = 0; this_l1track < TTTrackHandle->size(); this_l1track++) {
163 float trk_pt = trkPtr->momentum().perp();
164 int trk_nstubs = (
int)trkPtr->getStubRefs().size();
165 float trk_chi2dof = trkPtr->chi2Red();
166 float trk_d0 = trkPtr->d0();
167 float trk_bendchi2 = trkPtr->stubPtConsistency();
168 float trk_z0 = trkPtr->z0();
171 for (
int istub = 0; istub < trk_nstubs; istub++) {
172 DetId detId(trkPtr->getStubRefs().at(istub)->getDetId());
219 std::vector<float> zmins, zmaxs;
220 for (
int zbin = 0; zbin <
zBins_; zbin++) {
226 std::vector<std::vector<EtaPhiBin>> L1clusters;
228 std::vector<EtaPhiBin> L2clusters;
230 for (
unsigned int zbin = 0; zbin < zmins.size(); ++zbin) {
232 float zmin = zmins[zbin];
233 float zmax = zmaxs[zbin];
248 if (zbin == 0 &&
zmin == trkZ)
260 if ((trketa <
eta_min) || (trketa >
eta_max) || (trkphi < phi_min) || (trkphi > phi_max))
275 for (
int phibin = 0; phibin <
phiBins_; ++phibin) {
284 for (
unsigned int k = 0;
k < L2clusters.size(); ++
k) {
290 sum_pt += L2clusters[
k].pTtot;
299 mzb.
nclust = L2clusters.size();
303 vector<Ptr<L1TTTrackType>> L1TrackAssocJet;
304 for (
unsigned int j = 0;
j < mzb.
clusters.size(); ++
j) {
312 int totalDisptrk = mzb.
clusters[
j].numtdtrks;
316 L1TrackAssocJet.clear();
317 for (
unsigned int itrk = 0; itrk < mzb.
clusters[
j].trackidx.size(); itrk++)
322 if (!L1TrackAssocJet.empty())
323 L1L1TrackJetProducer->push_back(trkJet);
bool trackQualityCuts(int trk_nstub, float trk_chi2, float trk_bendchi2)
const int highpTJetMinTrackMultiplicity_
unsigned int tobLayer(const DetId &id) const
std::vector< EtaPhiBin > L2_clustering(std::vector< std::vector< EtaPhiBin >> &L1clusters, int phiBins_, float phiStep_, float etaStep_)
const float d0CutNStubs4_
const float lowpTJetThreshold_
T const & getData(const ESGetToken< T, R > &iToken) const noexcept(false)
Sin< T >::type sin(const T &t)
edm::ESGetToken< TrackerTopology, TrackerTopologyRcd > tTopoToken_
std::vector< unsigned int > trackidx
XYZTLorentzVectorD XYZTLorentzVector
Lorentz vector with cylindrical internal representation using pseudorapidity.
const float d0CutNStubs5_
const float highpTJetThreshold_
vector< Ptr< L1TTTrackType > > L1TrkPtrs_
std::vector< EtaPhiBin > L1_clustering(EtaPhiBin *phislice, int etaBins_, float etaStep_)
Cos< T >::type cos(const T &t)
Abs< T >::type abs(const T &t)
static constexpr auto TOB
const edm::EDGetTokenT< std::vector< l1t::Vertex > > PVtxToken_
const int lowpTJetMinTrackMultiplicity_
const int nDisplacedTracks_
const EDGetTokenT< vector< TTTrack< Ref_Phase2TrackerDigi_ > > > trackToken_
unsigned int tidRing(const DetId &id) const
std::vector< EtaPhiBin > clusters
std::vector< TkJet > TkJetCollection
static constexpr auto TID
const double minTrkJetpT_