Inheritance diagram for L1TrackJetProducer:

Public Types
typedef 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
virtual EtaPhiBin *	L1_cluster (EtaPhiBin *phislice)

	L1TrackJetProducer (const ParameterSet &)

void	L2_cluster (vector< Ptr< L1TTTrackType > > L1TrkPtrs_, vector< int > ttrk_, vector< int > tdtrk_, vector< int > ttdtrk_, MaxZBin &mzb)

bool	trackQualityCuts (float trk_pt, int trk_nstub, float trk_chi2, float trk_bendchi2, float trk_d0)

	~L1TrackJetProducer () 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

Static Public Member Functions
static void	fillDescriptions (ConfigurationDescriptions &descriptions)

Private Member Functions
void	beginStream (StreamID) override

void	endStream () override

void	produce (Event &, const EventSetup &) override

Private Attributes
float	d0CutNStubs4_

float	d0CutNStubs5_

bool	displaced_

int	etaBins_

float	etaStep_

int	highpTJetMinTrackMultiplicity_

vector< Ptr< L1TTTrackType > >	L1TrkPtrs_

int	lowpTJetMinTrackMultiplicity_

double	minJetEtHighPt_

double	minJetEtLowPt_

double	minTrkJetpT_

float	nStubs4DisplacedBendLoose_

float	nStubs4DisplacedBendTight_

float	nStubs4DisplacedChi2Loose_

float	nStubs4DisplacedChi2Tight_

float	nStubs5DisplacedBendLoose_

float	nStubs5DisplacedBendTight_

float	nStubs5DisplacedChi2Loose_

float	nStubs5DisplacedChi2Tight_

int	phiBins_

float	phiStep_

vector< int >	tdtrk_

const EDGetTokenT< vector < TTTrack < Ref_Phase2TrackerDigi_ > > >	trackToken_

float	trkBendChi2Max_

float	trkChi2dofMax_

float	trkEtaMax_

int	trkNPSStubMin_

float	trkPtMax_

float	trkPtMin_

float	trkZMax_

vector< int >	ttdtrk_

edm::ESGetToken < TrackerTopology, TrackerTopologyRcd >	tTopoToken_

vector< int >	ttrk_

vector< int >	zBinCount_

int	zBins_

float	zStep_

Detailed Description

Definition at line 43 of file L1TrackJetProducer.cc.

Member Typedef Documentation

typedef vector<L1TTTrackType> L1TrackJetProducer::L1TTTrackCollectionType

Definition at line 48 of file L1TrackJetProducer.cc.

typedef TTTrack<Ref_Phase2TrackerDigi_> L1TrackJetProducer::L1TTTrackType

Definition at line 47 of file L1TrackJetProducer.cc.

Constructor & Destructor Documentation

L1TrackJetProducer::L1TrackJetProducer ( const ParameterSet & iConfig )

explicit

Definition at line 102 of file L1TrackJetProducer.cc.

References d0CutNStubs4_, d0CutNStubs5_, displaced_, etaBins_, etaStep_, edm::ParameterSet::getParameter(), highpTJetMinTrackMultiplicity_, lowpTJetMinTrackMultiplicity_, M_PI, minJetEtHighPt_, minJetEtLowPt_, minTrkJetpT_, nStubs4DisplacedBendLoose_, nStubs4DisplacedBendTight_, nStubs4DisplacedChi2Loose_, nStubs4DisplacedChi2Tight_, nStubs5DisplacedBendLoose_, nStubs5DisplacedBendTight_, nStubs5DisplacedChi2Loose_, nStubs5DisplacedChi2Tight_, phiBins_, phiStep_, trkBendChi2Max_, trkChi2dofMax_, trkEtaMax_, trkNPSStubMin_, trkPtMax_, trkPtMin_, trkZMax_, zBins_, and zStep_.

     : trackToken_(
           consumes<vector<TTTrack<Ref_Phase2TrackerDigi_> > >(iConfig.getParameter<InputTag>("L1TrackInputTag"))),
       tTopoToken_(esConsumes<TrackerTopology, TrackerTopologyRcd>(edm::ESInputTag("", ""))) {
   trkZMax_ = (float)iConfig.getParameter<double>("trk_zMax");
   trkPtMax_ = (float)iConfig.getParameter<double>("trk_ptMax");
   trkPtMin_ = (float)iConfig.getParameter<double>("trk_ptMin");
   trkEtaMax_ = (float)iConfig.getParameter<double>("trk_etaMax");
   trkChi2dofMax_ = (float)iConfig.getParameter<double>("trk_chi2dofMax");
   trkBendChi2Max_ = (float)iConfig.getParameter<double>("trk_bendChi2Max");
   trkNPSStubMin_ = (int)iConfig.getParameter<int>("trk_nPSStubMin");
   minTrkJetpT_ = iConfig.getParameter<double>("minTrkJetpT");
   minJetEtLowPt_ = iConfig.getParameter<double>("minJetEtLowPt");
   minJetEtHighPt_ = iConfig.getParameter<double>("minJetEtHighPt");
   etaBins_ = (int)iConfig.getParameter<int>("etaBins");
   phiBins_ = (int)iConfig.getParameter<int>("phiBins");
   zBins_ = (int)iConfig.getParameter<int>("zBins");
   d0CutNStubs4_ = (float)iConfig.getParameter<double>("d0_cutNStubs4");
   d0CutNStubs5_ = (float)iConfig.getParameter<double>("d0_cutNStubs5");
   lowpTJetMinTrackMultiplicity_ = (int)iConfig.getParameter<int>("lowpTJetMinTrackMultiplicity");
   highpTJetMinTrackMultiplicity_ = (int)iConfig.getParameter<int>("highpTJetMinTrackMultiplicity");
   displaced_ = iConfig.getParameter<bool>("displaced");
   nStubs4DisplacedChi2Loose_ = (float)iConfig.getParameter<double>("nStubs4DisplacedChi2_Loose");
   nStubs5DisplacedChi2Loose_ = (float)iConfig.getParameter<double>("nStubs5DisplacedChi2_Loose");
   nStubs4DisplacedBendLoose_ = (float)iConfig.getParameter<double>("nStubs4Displacedbend_Loose");
   nStubs5DisplacedBendLoose_ = (float)iConfig.getParameter<double>("nStubs5Displacedbend_Loose");
   nStubs4DisplacedChi2Tight_ = (float)iConfig.getParameter<double>("nStubs4DisplacedChi2_Tight");
   nStubs5DisplacedChi2Tight_ = (float)iConfig.getParameter<double>("nStubs5DisplacedChi2_Tight");
   nStubs4DisplacedBendTight_ = (float)iConfig.getParameter<double>("nStubs4Displacedbend_Tight");
   nStubs5DisplacedBendTight_ = (float)iConfig.getParameter<double>("nStubs5Displacedbend_Tight");
 
   zStep_ = 2.0 * trkZMax_ / zBins_;
   etaStep_ = 2.0 * trkEtaMax_ / etaBins_;  //etaStep is the width of an etabin
   phiStep_ = 2 * M_PI / phiBins_;          
 
   if (displaced_)
     produces<TkJetCollection>("L1TrackJetsExtended");
   else
     produces<TkJetCollection>("L1TrackJets");
 }

L1TrackJetProducer::~L1TrackJetProducer ( )

override

Definition at line 143 of file L1TrackJetProducer.cc.

143 {}

Member Function Documentation

void L1TrackJetProducer::beginStream ( StreamID )

overrideprivate

Definition at line 642 of file L1TrackJetProducer.cc.

642 {}

void L1TrackJetProducer::endStream ( )

overrideprivate

Definition at line 644 of file L1TrackJetProducer.cc.

644 {}

void L1TrackJetProducer::fillDescriptions ( ConfigurationDescriptions & descriptions )

static

Definition at line 662 of file L1TrackJetProducer.cc.

References edm::ConfigurationDescriptions::addDefault(), submitPVResolutionJobs::desc, and edm::ParameterSetDescription::setUnknown().

                                                                                  {
   //The following says we do not know what parameters are allowed so do no validation
   // Please change this to state exactly what you do use, even if it is no parameters
   ParameterSetDescription desc;
   desc.setUnknown();
   descriptions.addDefault(desc);
 }

EtaPhiBin * L1TrackJetProducer::L1_cluster ( EtaPhiBin * phislice )

virtual

Definition at line 554 of file L1TrackJetProducer.cc.

References HLT_FULL_cff::clusters, EtaPhiBin::eta, PVValHelper::eta, etaBins_, etaStep_, mps_fire::i, visualization-live-secondInstance_cfg::m, DeadROCCounter::nclust, EtaPhiBin::numtdtrks, EtaPhiBin::numtracks, EtaPhiBin::numttdtrks, EtaPhiBin::numttrks, EtaPhiBin::pTtot, and EtaPhiBin::used.

Referenced by L2_cluster().

                                                              {
   EtaPhiBin *clusters = new EtaPhiBin[etaBins_ / 2];
   if (clusters == nullptr)
     edm::LogWarning("L1TrackJetProducer") << "Clusters memory not assigned!\n";
 
   // Find eta-phibin with maxpT, make center of cluster, add neighbors if not already used
   float my_pt, left_pt, right_pt, right2pt;
   int nclust = 0;
   right2pt = 0;
   for (int etabin = 0; etabin < etaBins_; ++etabin) {
     // assign values for my pT and neighbors' pT
     if (phislice[etabin].used)
       continue;
     my_pt = phislice[etabin].pTtot;
     if (etabin > 0 && !phislice[etabin - 1].used) {
       left_pt = phislice[etabin - 1].pTtot;
     } else
       left_pt = 0;
     if (etabin < etaBins_ - 1 && !phislice[etabin + 1].used) {
       right_pt = phislice[etabin + 1].pTtot;
       if (etabin < etaBins_ - 2 && !phislice[etabin + 2].used) {
         right2pt = phislice[etabin + 2].pTtot;
       } else
         right2pt = 0;
     } else
       right_pt = 0;
 
     // if I'm not a cluster, move on
     if (my_pt < left_pt || my_pt <= right_pt) {
       // if unused pT in the left neighbor, spit it out as a cluster
       if (left_pt > 0) {
         clusters[nclust] = phislice[etabin - 1];
         phislice[etabin - 1].used = true;
         nclust++;
       }
       continue;
     }
 
     // I guess I'm a cluster-- should I use my right neighbor?
     // Note: left neighbor will definitely be used because if it
     //       didn't belong to me it would have been used already
     clusters[nclust] = phislice[etabin];
     phislice[etabin].used = true;
     if (left_pt > 0) {
       if (clusters != nullptr) {
         clusters[nclust].pTtot += left_pt;
         clusters[nclust].numtracks += phislice[etabin - 1].numtracks;
         clusters[nclust].numttrks += phislice[etabin - 1].numttrks;
         clusters[nclust].numtdtrks += phislice[etabin - 1].numtdtrks;
         clusters[nclust].numttdtrks += phislice[etabin - 1].numttdtrks;
       }
     }
     if (my_pt >= right2pt && right_pt > 0) {
       if (clusters != nullptr) {
         clusters[nclust].pTtot += right_pt;
         clusters[nclust].numtracks += phislice[etabin + 1].numtracks;
         clusters[nclust].numttrks += phislice[etabin + 1].numttrks;
         clusters[nclust].numtdtrks += phislice[etabin + 1].numtdtrks;
         clusters[nclust].numttdtrks += phislice[etabin + 1].numttdtrks;
         phislice[etabin + 1].used = true;
       }
     }
     nclust++;
   }  // for each etabin
 
   // Now merge clusters, if necessary
   for (int m = 0; m < nclust - 1; ++m) {
     if (fabs(clusters[m + 1].eta - clusters[m].eta) < 1.5 * etaStep_) {
       if (clusters[m + 1].pTtot > clusters[m].pTtot) {
         clusters[m].eta = clusters[m + 1].eta;
       }
       clusters[m].pTtot += clusters[m + 1].pTtot;
       clusters[m].numtracks += clusters[m + 1].numtracks;  // Previous version didn't add tracks when merging
       clusters[m].numttrks += clusters[m + 1].numttrks;
       clusters[m].numtdtrks += clusters[m + 1].numtdtrks;
       clusters[m].numttdtrks += clusters[m + 1].numttdtrks;
       for (int m1 = m + 1; m1 < nclust - 1; ++m1)
         clusters[m1] = clusters[m1 + 1];
       nclust--;
       m = -1;
     }  // end if clusters neighbor in eta
   }    // end for (m) loop
 
   for (int i = nclust; i < etaBins_ / 2; ++i)  // zero out remaining unused clusters
     clusters[i].pTtot = 0;
   return clusters;
 }

void L1TrackJetProducer::L2_cluster	(	vector< Ptr< L1TTTrackType > >	L1TrkPtrs_,
		vector< int >	ttrk_,
		vector< int >	tdtrk_,
		vector< int >	ttdtrk_,
		MaxZBin &	mzb
	)

Definition at line 271 of file L1TrackJetProducer.cc.

References MaxZBin::clusters, EtaPhiBin::eta, PVValHelper::eta, etaBins_, etaStep_, highpTJetMinTrackMultiplicity_, MaxZBin::ht, mps_fire::i, dqmiolumiharvest::j, isotrackApplyRegressor::k, L1_cluster(), L1TrkPtrs_, lowpTJetMinTrackMultiplicity_, visualization-live-secondInstance_cfg::m, M_PI, minJetEtHighPt_, minJetEtLowPt_, minTrkJetpT_, dqmiodumpmetadata::n, MaxZBin::nclust, DeadROCCounter::nclust, EtaPhiBin::numtdtrks, EtaPhiBin::numtracks, EtaPhiBin::numttdtrks, EtaPhiBin::numttrks, phi, EtaPhiBin::phi, phiBins_, phimax, phimin, phiStep_, EtaPhiBin::pTtot, trkEtaMax_, trkPtMax_, trkZMax_, EtaPhiBin::used, z, MaxZBin::zbincenter, zBinCount_, zBins_, MaxZBin::znum, and zStep_.

Referenced by produce().

                                                                                                                       {
   const int nz = zBins_;
   MaxZBin all_zBins[nz];
   for (int z = 0; z < nz; ++z)
     all_zBins[z] = MaxZBin{0, 0, 0, nullptr, 0};
 
   float zmin = -1.0 * trkZMax_;
   float zmax = zmin + 2 * zStep_;
 
   EtaPhiBin epbins[phiBins_][etaBins_];  // create grid of phiBins
   float phi = -1.0 * M_PI;
   float eta;
   float etamin, etamax, phimin, phimax;
   for (int i = 0; i < phiBins_; ++i) {
     eta = -1.0 * trkEtaMax_;
     for (int j = 0; j < etaBins_; ++j) {
       phimin = phi;
       phimax = phi + phiStep_;
       etamin = eta;
       eta = eta + etaStep_;
       etamax = eta;
       epbins[i][j].phi = (phimin + phimax) / 2.0;
       epbins[i][j].eta = (etamin + etamax) / 2.0;
     }  // for each etabin
     phi = phi + phiStep_;
   }  // for each phibin (finished creating epbins)
   mzb = all_zBins[0];
   int ntracks = L1TrkPtrs_.size();
   // uninitalized arrays
   EtaPhiBin *L1clusters[phiBins_];
   EtaPhiBin L2cluster[ntracks];
 
   for (int zbin = 0; zbin < zBins_ - 1; ++zbin) {
     for (int i = 0; i < phiBins_; ++i) {  //First initialize pT, numtracks, used to 0 (or false)
       for (int j = 0; j < etaBins_; ++j) {
         epbins[i][j].pTtot = 0;
         epbins[i][j].used = false;
         epbins[i][j].numtracks = 0;
         epbins[i][j].numttrks = 0;
         epbins[i][j].numtdtrks = 0;
         epbins[i][j].numttdtrks = 0;
       }  //for each etabin
       L1clusters[i] = epbins[i];
     }  //for each phibin
 
     for (unsigned int k = 0; k < L1TrkPtrs_.size(); ++k) {
       float trkpt = L1TrkPtrs_[k]->momentum().perp();
       float trketa = L1TrkPtrs_[k]->momentum().eta();
       float trkphi = L1TrkPtrs_[k]->momentum().phi();
       float trkZ = L1TrkPtrs_[k]->z0();
 
       for (int i = 0; i < phiBins_; ++i) {
         for (int j = 0; j < etaBins_; ++j) {
           L2cluster[k] = epbins[i][j];
           if ((zmin <= trkZ && zmax >= trkZ) &&
               ((epbins[i][j].eta - etaStep_ / 2.0 <= trketa && epbins[i][j].eta + etaStep_ / 2.0 >= trketa) &&
                epbins[i][j].phi - phiStep_ / 2.0 <= trkphi && epbins[i][j].phi + phiStep_ / 2.0 >= trkphi &&
                (zBinCount_[k] != 2))) {
             zBinCount_.at(k) = zBinCount_.at(k) + 1;
             if (trkpt < trkPtMax_)
               epbins[i][j].pTtot += trkpt;
             else
               epbins[i][j].pTtot += trkPtMax_;
             epbins[i][j].numttrks += ttrk_[k];
             epbins[i][j].numtdtrks += tdtrk_[k];
             epbins[i][j].numttdtrks += ttdtrk_[k];
             ++epbins[i][j].numtracks;
           }  // if right bin
         }    // for each phibin: j loop
       }      // for each phibin: i loop
     }        // end loop over tracks
 
     for (int phislice = 0; phislice < phiBins_; ++phislice) {
       L1clusters[phislice] = L1_cluster(epbins[phislice]);
       if (L1clusters[phislice] != nullptr) {
         for (int ind = 0; L1clusters[phislice][ind].pTtot != 0; ++ind) {
           L1clusters[phislice][ind].used = false;
         }
       }
     }
 
     //Create clusters array to hold output cluster data for Layer2; can't have more clusters than tracks.
     //Find eta-phibin with maxpT, make center of cluster, add neighbors if not already used.
     float hipT = 0;
     int nclust = 0;
     int phibin = 0;
     int imax = -1;
     int index1;  //index of clusters array for each phislice
     float E1 = 0;
     float E0 = 0;
     float E2 = 0;
     int trx1, trx2;
     int ttrk1, ttrk2;
     int tdtrk1, tdtrk2;
     int ttdtrk1, ttdtrk2;
     int used1, used2, used3, used4;
 
     for (phibin = 0; phibin < phiBins_; ++phibin) {  //Find eta-phibin with highest pT
       while (true) {
         hipT = 0;
         for (index1 = 0; L1clusters[phibin][index1].pTtot > 0; ++index1) {
           if (!L1clusters[phibin][index1].used && L1clusters[phibin][index1].pTtot >= hipT) {
             hipT = L1clusters[phibin][index1].pTtot;
             imax = index1;
           }
         }  // for each index within the phibin
 
         if (hipT == 0)
           break;    // If highest pT is 0, all bins are used
         E0 = hipT;  // E0 is pT of first phibin of the cluster
         E1 = 0;
         E2 = 0;
         trx1 = 0;
         trx2 = 0;
         ttrk1 = 0;
         ttrk2 = 0;
         tdtrk1 = 0;
         tdtrk2 = 0;
         ttdtrk1 = 0;
         ttdtrk2 = 0;
         L2cluster[nclust] = L1clusters[phibin][imax];
         L1clusters[phibin][imax].used = true;
         // Add pT of upper neighbor
         // E1 is pT of the middle phibin (should be highest pT)
         if (phibin != phiBins_ - 1) {
           used1 = -1;
           used2 = -1;
           for (index1 = 0; L1clusters[phibin + 1][index1].pTtot != 0; ++index1) {
             if (L1clusters[phibin + 1][index1].used)
               continue;
             if (fabs(L1clusters[phibin + 1][index1].eta - L1clusters[phibin][imax].eta) <= 1.5 * etaStep_) {
               E1 += L1clusters[phibin + 1][index1].pTtot;
               trx1 += L1clusters[phibin + 1][index1].numtracks;
               ttrk1 += L1clusters[phibin + 1][index1].numttrks;
               tdtrk1 += L1clusters[phibin + 1][index1].numtdtrks;
               ttdtrk1 += L1clusters[phibin + 1][index1].numttdtrks;
               if (used1 < 0)
                 used1 = index1;
               else
                 used2 = index1;
             }  // if cluster is within one phibin
           }    // for each cluster in above phibin
 
           if (E1 < E0) {  // if E1 isn't higher, E0 and E1 are their own cluster
             L2cluster[nclust].pTtot += E1;
             L2cluster[nclust].numtracks += trx1;
             L2cluster[nclust].numttrks += ttrk1;
             L2cluster[nclust].numtdtrks += tdtrk1;
             L2cluster[nclust].numttdtrks += ttdtrk1;
             if (used1 >= 0)
               L1clusters[phibin + 1][used1].used = true;
             if (used2 >= 0)
               L1clusters[phibin + 1][used2].used = true;
             nclust++;
             continue;
           }
 
           if (phibin != phiBins_ - 2) {  // E2 will be the pT of the third phibin (should be lower than E1)
             used3 = -1;
             used4 = -1;
             for (index1 = 0; L1clusters[phibin + 2][index1].pTtot != 0; ++index1) {
               if (L1clusters[phibin + 2][index1].used)
                 continue;
               if (fabs(L1clusters[phibin + 2][index1].eta - L1clusters[phibin][imax].eta) <= 1.5 * etaStep_) {
                 E2 += L1clusters[phibin + 2][index1].pTtot;
                 trx2 += L1clusters[phibin + 2][index1].numtracks;
                 ttrk2 += L1clusters[phibin + 2][index1].numttrks;
                 tdtrk2 += L1clusters[phibin + 2][index1].numtdtrks;
                 ttdtrk2 += L1clusters[phibin + 2][index1].numttdtrks;
                 if (used3 < 0)
                   used3 = index1;
                 else
                   used4 = index1;
               }
             }
             // if indeed E2 < E1, add E1 and E2 to E0, they're all a cluster together
             // otherwise, E0 is its own cluster
             if (E2 < E1) {
               L2cluster[nclust].pTtot += E1 + E2;
               L2cluster[nclust].numtracks += trx1 + trx2;
               L2cluster[nclust].numttrks += ttrk1 + ttrk2;
               L2cluster[nclust].numtdtrks += tdtrk1 + tdtrk2;
               L2cluster[nclust].numttdtrks += ttdtrk1 + ttdtrk2;
               L2cluster[nclust].phi = L1clusters[phibin + 1][used1].phi;
               if (used1 >= 0)
                 L1clusters[phibin + 1][used1].used = true;
               if (used2 >= 0)
                 L1clusters[phibin + 1][used2].used = true;
               if (used3 >= 0)
                 L1clusters[phibin + 2][used3].used = true;
               if (used4 >= 0)
                 L1clusters[phibin + 2][used4].used = true;
             }
             nclust++;
             continue;
           }  // end Not phiBins-2
           else {
             L2cluster[nclust].pTtot += E1;
             L2cluster[nclust].numtracks += trx1;
             L2cluster[nclust].numttrks += ttrk1;
             L2cluster[nclust].numtdtrks += tdtrk1;
             L2cluster[nclust].numttdtrks += ttdtrk1;
             L2cluster[nclust].phi = L1clusters[phibin + 1][used1].phi;
             if (used1 >= 0)
               L1clusters[phibin + 1][used1].used = true;
             if (used2 >= 0)
               L1clusters[phibin + 1][used2].used = true;
             nclust++;
             continue;
           }
         }       //End not last phibin(23)
         else {  //if it is phibin 23
           L1clusters[phibin][imax].used = true;
           nclust++;
         }
       }  // while hipT not 0
     }    // for each phibin
 
     for (phibin = 0; phibin < phiBins_; ++phibin)
       delete[] L1clusters[phibin];
 
     // Now merge clusters, if necessary
     for (int m = 0; m < nclust - 1; ++m) {
       for (int n = m + 1; n < nclust; ++n)
         if (L2cluster[n].eta == L2cluster[m].eta && (fabs(L2cluster[n].phi - L2cluster[m].phi) < 1.5 * phiStep_ ||
                                                      fabs(L2cluster[n].phi - L2cluster[m].phi) > 6.0)) {
           if (L2cluster[n].pTtot > L2cluster[m].pTtot)
             L2cluster[m].phi = L2cluster[n].phi;
           L2cluster[m].pTtot += L2cluster[n].pTtot;
           L2cluster[m].numtracks += L2cluster[n].numtracks;
           L2cluster[m].numttrks += L2cluster[n].numttrks;
           L2cluster[m].numtdtrks += L2cluster[n].numtdtrks;
           L2cluster[m].numttdtrks += L2cluster[n].numttdtrks;
           for (int m1 = n; m1 < nclust - 1; ++m1)
             L2cluster[m1] = L2cluster[m1 + 1];
           nclust--;
           m = -1;
           break;  //?????
         }         // end if clusters neighbor in eta
     }             // end for (m) loop
 
     // sum up all pTs in this zbin to find ht
     float ht = 0;
     for (int k = 0; k < nclust; ++k) {
       if (L2cluster[k].pTtot > minJetEtLowPt_ && L2cluster[k].numtracks < lowpTJetMinTrackMultiplicity_)
         continue;
       if (L2cluster[k].pTtot > minJetEtHighPt_ && L2cluster[k].numtracks < highpTJetMinTrackMultiplicity_)
         continue;
       if (L2cluster[k].pTtot > minTrkJetpT_)
         ht += L2cluster[k].pTtot;
     }
 
     // if ht is larger than previous max, this is the new vertex zbin
     all_zBins[zbin].znum = zbin;
     all_zBins[zbin].clusters = new EtaPhiBin[nclust];
     all_zBins[zbin].nclust = nclust;
     all_zBins[zbin].zbincenter = (zmin + zmax) / 2.0;
     for (int k = 0; k < nclust; ++k) {
       all_zBins[zbin].clusters[k].phi = L2cluster[k].phi;
       all_zBins[zbin].clusters[k].eta = L2cluster[k].eta;
       all_zBins[zbin].clusters[k].pTtot = L2cluster[k].pTtot;
       all_zBins[zbin].clusters[k].numtracks = L2cluster[k].numtracks;
       all_zBins[zbin].clusters[k].numttrks = L2cluster[k].numttrks;
       all_zBins[zbin].clusters[k].numtdtrks = L2cluster[k].numtdtrks;
       all_zBins[zbin].clusters[k].numttdtrks = L2cluster[k].numttdtrks;
     }
     all_zBins[zbin].ht = ht;
     if (ht >= mzb.ht) {
       mzb = all_zBins[zbin];
       mzb.zbincenter = (zmin + zmax) / 2.0;
     }
     // Prepare for next zbin!
     zmin = zmin + zStep_;
     zmax = zmax + zStep_;
   }  // for each zbin
   for (int zbin = 0; zbin < zBins_ - 1; ++zbin) {
     if (zbin == mzb.znum)
       continue;
     delete[] all_zBins[zbin].clusters;
   }
 }

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

overrideprivate

Definition at line 145 of file L1TrackJetProducer.cc.

References MaxZBin::clusters, funct::cos(), d0CutNStubs4_, d0CutNStubs5_, displaced_, EtaPhiBin::eta, etaStep_, edm::Event::getByToken(), edm::EventSetup::getData(), dqmiolumiharvest::j, reco::btau::jetEta, reco::btau::jetPhi, reco::btau::jetPt, L1TrkPtrs_, L2_cluster(), eostools::move(), MaxZBin::nclust, nStubs4DisplacedBendTight_, nStubs4DisplacedChi2Tight_, nStubs5DisplacedBendTight_, nStubs5DisplacedChi2Tight_, EtaPhiBin::numtdtrks, EtaPhiBin::numtracks, EtaPhiBin::numttdtrks, EtaPhiBin::numttrks, EtaPhiBin::phi, phiStep_, EtaPhiBin::pTtot, edm::Event::put(), funct::sin(), submitPVValidationJobs::t, tdtrk_, StripSubdetector::TID, TrackerTopology::tidRing(), StripSubdetector::TOB, TrackerTopology::tobLayer(), align::Tracker, trackQualityCuts(), trackToken_, trkEtaMax_, trkNPSStubMin_, trkPtMin_, trkZMax_, ttdtrk_, tTopoToken_, ttrk_, MaxZBin::zbincenter, zBinCount_, and zStep_.

                                                                         {
   unique_ptr<TkJetCollection> L1L1TrackJetProducer(new TkJetCollection);
 
   // For TTStubs
   const TrackerTopology &tTopo = iSetup.getData(tTopoToken_);
 
   edm::Handle<vector<TTTrack<Ref_Phase2TrackerDigi_> > > TTTrackHandle;
   iEvent.getByToken(trackToken_, TTTrackHandle);
   vector<TTTrack<Ref_Phase2TrackerDigi_> >::const_iterator iterL1Track;
 
   L1TrkPtrs_.clear();
   zBinCount_.clear();
   ttrk_.clear();
   tdtrk_.clear();
   ttdtrk_.clear();
 
   unsigned int this_l1track = 0;
   for (iterL1Track = TTTrackHandle->begin(); iterL1Track != TTTrackHandle->end(); iterL1Track++) {
     edm::Ptr<L1TTTrackType> trkPtr(TTTrackHandle, this_l1track);
     this_l1track++;
     float trk_pt = trkPtr->momentum().perp();
     int trk_nstubs = (int)trkPtr->getStubRefs().size();
     float trk_chi2dof = trkPtr->chi2Red();
     float trk_d0 = trkPtr->d0();
     float trk_bendchi2 = trkPtr->stubPtConsistency();
 
     int trk_nPS = 0;
     for (int istub = 0; istub < trk_nstubs; istub++) {  // loop over the stubs
       DetId detId(trkPtr->getStubRefs().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))
           trk_nPS++;
       }
     }
 
     if (trk_nPS < trkNPSStubMin_)
       continue;
     if (!trackQualityCuts(trk_pt, trk_nstubs, trk_chi2dof, trk_bendchi2, fabs(trk_d0)))
       continue;
     if (fabs(iterL1Track->z0()) > trkZMax_)
       continue;
     if (fabs(iterL1Track->momentum().eta()) > trkEtaMax_)
       continue;
     if (trk_pt < trkPtMin_)
       continue;
     L1TrkPtrs_.push_back(trkPtr);
     zBinCount_.push_back(0);
 
     if ((fabs(trk_d0) > d0CutNStubs5_ && trk_nstubs >= 5) || (trk_nstubs == 4 && fabs(trk_d0) > d0CutNStubs4_))
       tdtrk_.push_back(1);
     else
       tdtrk_.push_back(0);  //displaced track
     if ((trk_nstubs >= 5 && trk_chi2dof < nStubs5DisplacedChi2Tight_ && trk_bendchi2 < nStubs5DisplacedBendTight_) ||
         (trk_nstubs == 4 && trk_chi2dof < nStubs4DisplacedChi2Tight_ && trk_bendchi2 < nStubs4DisplacedBendTight_))
       ttrk_.push_back(1);
     else
       ttrk_.push_back(0);
     if ((trk_nstubs >= 5 && trk_chi2dof < nStubs5DisplacedChi2Tight_ && trk_bendchi2 < nStubs5DisplacedBendTight_ &&
          fabs(trk_d0) > d0CutNStubs5_) ||
         (trk_nstubs == 4 && trk_chi2dof < nStubs4DisplacedChi2Tight_ && trk_bendchi2 < nStubs4DisplacedBendTight_ &&
          fabs(trk_d0) > d0CutNStubs4_))
       ttdtrk_.push_back(1);
     else
       ttdtrk_.push_back(0);
   }
 
   if (!L1TrkPtrs_.empty()) {
     MaxZBin mzb;
 
     L2_cluster(L1TrkPtrs_, ttrk_, tdtrk_, ttdtrk_, mzb);
     edm::Ref<JetBxCollection> jetRef;  //null, no Calo Jet Ref
     vector<Ptr<L1TTTrackType> > L1TrackAssocJet;
     if (mzb.clusters != nullptr) {
       for (int j = 0; j < mzb.nclust; ++j) {
         //FILL Two Layer Jets for Jet Collection
         if (mzb.clusters[j].pTtot <= trkPtMin_)
           continue;  //protects against reading bad memory
         if (mzb.clusters[j].numtracks < 1)
           continue;
         if (mzb.clusters[j].numtracks > 5000)
           continue;
         float jetEta = mzb.clusters[j].eta;
         float jetPhi = mzb.clusters[j].phi;
         float jetPt = mzb.clusters[j].pTtot;
         float jetPx = jetPt * cos(jetPhi);
         float jetPy = jetPt * sin(jetPhi);
         float jetPz = jetPt * sinh(jetEta);
         float jetP = jetPt * cosh(jetEta);
         int totalTighttrk_ = mzb.clusters[j].numttrks;
         int totalDisptrk = mzb.clusters[j].numtdtrks;
         int totalTightDisptrk = mzb.clusters[j].numttdtrks;
 
         math::XYZTLorentzVector jetP4(jetPx, jetPy, jetPz, jetP);
         L1TrackAssocJet.clear();
         for (unsigned int t = 0; t < L1TrkPtrs_.size(); ++t) {
           if (L1TrackAssocJet.size() == (unsigned int)mzb.clusters[j].numtracks)
             break;
           float deta = L1TrkPtrs_[t]->momentum().eta() - jetEta;
           float dphi = L1TrkPtrs_[t]->momentum().phi() - jetPhi;
           float dZ = fabs(mzb.zbincenter - L1TrkPtrs_[t]->z0());
           if (dZ < zStep_ && fabs(deta) < etaStep_ * 2.0 && fabs(dphi) < phiStep_ * 2.0) {
             L1TrackAssocJet.push_back(L1TrkPtrs_[t]);
           }
         }
         TkJet trkJet(jetP4,
                      L1TrackAssocJet,
                      mzb.zbincenter,
                      mzb.clusters[j].numtracks,
                      totalTighttrk_,
                      totalDisptrk,
                      totalTightDisptrk);
         //trkJet.setDispCounters(DispCounters);
         if (!L1TrackAssocJet.empty())
           L1L1TrackJetProducer->push_back(trkJet);
       }
     }
     //free(mzb.clusters);
     if (displaced_)
       iEvent.put(std::move(L1L1TrackJetProducer), "L1TrackJetsExtended");
     else
       iEvent.put(std::move(L1L1TrackJetProducer), "L1TrackJets");
     delete[] mzb.clusters;
   }
 }

bool L1TrackJetProducer::trackQualityCuts	(	float	trk_pt,
		int	trk_nstub,
		float	trk_chi2,
		float	trk_bendchi2,
		float	trk_d0
	)

Definition at line 646 of file L1TrackJetProducer.cc.

References d0CutNStubs4_, displaced_, nStubs4DisplacedBendLoose_, nStubs4DisplacedBendTight_, nStubs4DisplacedChi2Loose_, nStubs4DisplacedChi2Tight_, nStubs5DisplacedBendLoose_, and trkBendChi2Max_.

Referenced by produce().

                                                                                    {
   bool PassQuality = false;
   if (trk_bendchi2 < trkBendChi2Max_ && trk_chi2 < trkChi2dofMax_ && trk_nstub >= 4 && !displaced_)
     PassQuality = true;
   if (displaced_ && trk_bendchi2 < nStubs4DisplacedBendTight_ && trk_chi2 < nStubs4DisplacedChi2Tight_ &&
       trk_nstub == 4 && trk_d0 <= d0CutNStubs4_)
     PassQuality = true;
   if (displaced_ && trk_bendchi2 < nStubs4DisplacedBendLoose_ && trk_chi2 < nStubs4DisplacedChi2Loose_ &&
       trk_nstub == 4 && trk_d0 > d0CutNStubs4_)
     PassQuality = true;
   if (displaced_ && trk_bendchi2 < nStubs5DisplacedBendLoose_ && trk_chi2 < nStubs5DisplacedChi2Loose_ && trk_nstub > 4)
     PassQuality = true;
   return PassQuality;
 }