Inheritance diagram for L1TrackFastJetProducer:

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
	L1TrackFastJetProducer (const edm::ParameterSet &)

	~L1TrackFastJetProducer () 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 (edm::ConfigurationDescriptions &descriptions)

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

Private Attributes
double	coneSize_

double	deltaZ0Cut_

bool	displaced_

bool	doTightChi2_

edm::EDGetTokenT< TkPrimaryVertexCollection >	pvToken_

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

double	trkBendChi2Max_

float	trkChi2dofMax_

float	trkChi2dofTightChi2_

float	trkEtaMax_

int	trkNPSStubMin_

int	trkNStubMin_

float	trkPtMin_

float	trkPtTightChi2_

float	trkZMax_

edm::ESGetToken< TrackerTopology, TrackerTopologyRcd >	tTopoToken_

Detailed Description

Definition at line 51 of file L1TrackFastJetProducer.cc.

Member Typedef Documentation

◆ L1TTTrackCollectionType

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

Definition at line 54 of file L1TrackFastJetProducer.cc.

◆ L1TTTrackType

typedef TTTrack<Ref_Phase2TrackerDigi_> L1TrackFastJetProducer::L1TTTrackType

Definition at line 53 of file L1TrackFastJetProducer.cc.

Constructor & Destructor Documentation

◆ L1TrackFastJetProducer()

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

explicit

Definition at line 87 of file L1TrackFastJetProducer.cc.

     : trackToken_(consumes<std::vector<TTTrack<Ref_Phase2TrackerDigi_> > >(
           iConfig.getParameter<edm::InputTag>("L1TrackInputTag"))),
       pvToken_(consumes<TkPrimaryVertexCollection>(iConfig.getParameter<edm::InputTag>("L1PrimaryVertexTag"))),
       tTopoToken_(esConsumes<TrackerTopology, TrackerTopologyRcd>(edm::ESInputTag("", ""))) {
   trkZMax_ = (float)iConfig.getParameter<double>("trk_zMax");
   trkChi2dofMax_ = (float)iConfig.getParameter<double>("trk_chi2dofMax");
   trkBendChi2Max_ = iConfig.getParameter<double>("trk_bendChi2Max");
   trkPtMin_ = (float)iConfig.getParameter<double>("trk_ptMin");
   trkEtaMax_ = (float)iConfig.getParameter<double>("trk_etaMax");
   trkNStubMin_ = (int)iConfig.getParameter<int>("trk_nStubMin");
   trkNPSStubMin_ = (int)iConfig.getParameter<int>("trk_nPSStubMin");
   deltaZ0Cut_ = (float)iConfig.getParameter<double>("deltaZ0Cut");
   coneSize_ = (float)iConfig.getParameter<double>("coneSize");
   doTightChi2_ = iConfig.getParameter<bool>("doTightChi2");
   trkPtTightChi2_ = (float)iConfig.getParameter<double>("trk_ptTightChi2");
   trkChi2dofTightChi2_ = (float)iConfig.getParameter<double>("trk_chi2dofTightChi2");
   displaced_ = iConfig.getParameter<bool>("displaced");
   if (displaced_)
     produces<TkJetCollection>("L1TrackFastJetsExtended");
   else
     produces<TkJetCollection>("L1TrackFastJets");
 }

References coneSize_, deltaZ0Cut_, displaced_, doTightChi2_, dqmMemoryStats::float, edm::ParameterSet::getParameter(), createfilelist::int, trkBendChi2Max_, trkChi2dofMax_, trkChi2dofTightChi2_, trkEtaMax_, trkNPSStubMin_, trkNStubMin_, trkPtMin_, trkPtTightChi2_, and trkZMax_.

◆ ~L1TrackFastJetProducer()

L1TrackFastJetProducer::~L1TrackFastJetProducer ( )

override

Definition at line 112 of file L1TrackFastJetProducer.cc.

112 {}

Member Function Documentation

◆ fillDescriptions()

void L1TrackFastJetProducer::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 220 of file L1TrackFastJetProducer.cc.

                                                                                         {
   //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
   edm::ParameterSetDescription desc;
   desc.setUnknown();
   descriptions.addDefault(desc);
 }

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

◆ produce()

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

overrideprivate

Definition at line 115 of file L1TrackFastJetProducer.cc.

                                                                                   {
   std::unique_ptr<TkJetCollection> L1TrackFastJets(new TkJetCollection);
  
   // L1 tracks
   edm::Handle<std::vector<TTTrack<Ref_Phase2TrackerDigi_> > > TTTrackHandle;
   iEvent.getByToken(trackToken_, TTTrackHandle);
   std::vector<TTTrack<Ref_Phase2TrackerDigi_> >::const_iterator iterL1Track;
  
   // Tracker Topology
   const TrackerTopology& tTopo = iSetup.getData(tTopoToken_);
  
   edm::Handle<TkPrimaryVertexCollection> TkPrimaryVertexHandle;
   iEvent.getByToken(pvToken_, TkPrimaryVertexHandle);
  
   fastjet::JetDefinition jet_def(fastjet::antikt_algorithm, coneSize_);
   std::vector<fastjet::PseudoJet> JetInputs;
  
   float recoVtx = TkPrimaryVertexHandle->begin()->zvertex();
   unsigned int this_l1track = 0;
   for (iterL1Track = TTTrackHandle->begin(); iterL1Track != TTTrackHandle->end(); iterL1Track++) {
     this_l1track++;
     float trk_pt = iterL1Track->momentum().perp();
     float trk_z0 = iterL1Track->z0();
     float trk_chi2dof = iterL1Track->chi2Red();
     float trk_bendchi2 = iterL1Track->stubPtConsistency();
     std::vector<edm::Ref<edmNew::DetSetVector<TTStub<Ref_Phase2TrackerDigi_> >, TTStub<Ref_Phase2TrackerDigi_> > >
         theStubs = iterL1Track->getStubRefs();
     int trk_nstub = (int)theStubs.size();
  
     if (fabs(trk_z0) > trkZMax_)
       continue;
     if (fabs(iterL1Track->momentum().eta()) > trkEtaMax_)
       continue;
     if (trk_pt < trkPtMin_)
       continue;
     if (trk_nstub < trkNStubMin_)
       continue;
     if (trk_chi2dof > trkChi2dofMax_)
       continue;
     if (trk_bendchi2 > trkBendChi2Max_)
       continue;
     if (doTightChi2_ && (trk_pt > trkPtTightChi2_ && trk_chi2dof > trkChi2dofTightChi2_))
       continue;
  
     int trk_nPS = 0;
     for (int istub = 0; istub < trk_nstub; istub++) {
       DetId detId(theStubs.at(istub)->getDetId());
       bool tmp_isPS = false;
       if (detId.det() == DetId::Detector::Tracker) {
         if (detId.subdetId() == StripSubdetector::TOB && tTopo.tobLayer(detId) <= 3)
           tmp_isPS = true;
         else if (detId.subdetId() == StripSubdetector::TID && tTopo.tidRing(detId) <= 9)
           tmp_isPS = true;
       }
       if (tmp_isPS)
         trk_nPS++;
     }
     if (trk_nPS < trkNPSStubMin_)
       continue;
     if (fabs(recoVtx - trk_z0) > deltaZ0Cut_)
       continue;
  
     fastjet::PseudoJet psuedoJet(iterL1Track->momentum().x(),
                                  iterL1Track->momentum().y(),
                                  iterL1Track->momentum().z(),
                                  iterL1Track->momentum().mag());
     JetInputs.push_back(psuedoJet);                     // input tracks for clustering
     JetInputs.back().set_user_index(this_l1track - 1);  // save track index in the collection
   }                                                     // end loop over tracks
  
   fastjet::ClusterSequence cs(JetInputs, jet_def);  // define the output jet collection
   std::vector<fastjet::PseudoJet> JetOutputs =
       fastjet::sorted_by_pt(cs.inclusive_jets(0));  // output jet collection, pT-ordered
  
   for (unsigned int ijet = 0; ijet < JetOutputs.size(); ++ijet) {
     math::XYZTLorentzVector jetP4(
         JetOutputs[ijet].px(), JetOutputs[ijet].py(), JetOutputs[ijet].pz(), JetOutputs[ijet].modp());
     float sumpt = 0;
     float avgZ = 0;
     std::vector<edm::Ptr<L1TTTrackType> > L1TrackPtrs;
     std::vector<fastjet::PseudoJet> fjConstituents = fastjet::sorted_by_pt(cs.constituents(JetOutputs[ijet]));
  
     for (unsigned int i = 0; i < fjConstituents.size(); ++i) {
       auto index = fjConstituents[i].user_index();
       edm::Ptr<L1TTTrackType> trkPtr(TTTrackHandle, index);
       L1TrackPtrs.push_back(trkPtr);  // L1Tracks in the jet
       sumpt = sumpt + trkPtr->momentum().perp();
       avgZ = avgZ + trkPtr->momentum().perp() * trkPtr->z0();
     }
     avgZ = avgZ / sumpt;
     edm::Ref<JetBxCollection> jetRef;
     TkJet trkJet(jetP4, jetRef, L1TrackPtrs, avgZ);
     L1TrackFastJets->push_back(trkJet);
   }  //end loop over Jet Outputs
  
   if (displaced_)
     iEvent.put(std::move(L1TrackFastJets), "L1TrackFastJetsExtended");
   else
     iEvent.put(std::move(L1TrackFastJets), "L1TrackFastJets");
 }

References coneSize_, fwrapper::cs, deltaZ0Cut_, displaced_, doTightChi2_, edm::EventSetup::getData(), mps_fire::i, iEvent, createfilelist::int, fwrapper::jet_def, L1TrackFastJetProducer_cfi::L1TrackFastJets, TTTrack< T >::momentum(), eostools::move(), PV3DBase< T, PVType, FrameType >::perp(), pvToken_, multPhiCorr_741_25nsDY_cfi::px, multPhiCorr_741_25nsDY_cfi::py, StripSubdetector::TID, TrackerTopology::tidRing(), StripSubdetector::TOB, TrackerTopology::tobLayer(), align::Tracker, trackToken_, trkBendChi2Max_, trkChi2dofMax_, trkChi2dofTightChi2_, trkEtaMax_, trkNPSStubMin_, trkNStubMin_, trkPtMin_, trkPtTightChi2_, trkZMax_, tTopoToken_, and TTTrack< T >::z0().

Member Data Documentation

◆ coneSize_

double L1TrackFastJetProducer::coneSize_

private

Definition at line 75 of file L1TrackFastJetProducer.cc.

Referenced by L1TrackFastJetProducer(), and produce().

◆ deltaZ0Cut_

double L1TrackFastJetProducer::deltaZ0Cut_

private

Definition at line 74 of file L1TrackFastJetProducer.cc.

Referenced by L1TrackFastJetProducer(), and produce().

◆ displaced_

bool L1TrackFastJetProducer::displaced_

private

Definition at line 79 of file L1TrackFastJetProducer.cc.

Referenced by L1TrackFastJetProducer(), and produce().

◆ doTightChi2_

bool L1TrackFastJetProducer::doTightChi2_

private

Definition at line 76 of file L1TrackFastJetProducer.cc.

Referenced by L1TrackFastJetProducer(), and produce().

◆ pvToken_

edm::EDGetTokenT<TkPrimaryVertexCollection> L1TrackFastJetProducer::pvToken_

private

Definition at line 82 of file L1TrackFastJetProducer.cc.

Referenced by produce().

◆ trackToken_

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

private

Definition at line 81 of file L1TrackFastJetProducer.cc.

Referenced by produce().

◆ trkBendChi2Max_

double L1TrackFastJetProducer::trkBendChi2Max_

private

Definition at line 69 of file L1TrackFastJetProducer.cc.

Referenced by L1TrackFastJetProducer(), and produce().

◆ trkChi2dofMax_

float L1TrackFastJetProducer::trkChi2dofMax_

private

Definition at line 68 of file L1TrackFastJetProducer.cc.

Referenced by L1TrackFastJetProducer(), and produce().

◆ trkChi2dofTightChi2_

float L1TrackFastJetProducer::trkChi2dofTightChi2_

private

Definition at line 78 of file L1TrackFastJetProducer.cc.

Referenced by L1TrackFastJetProducer(), and produce().

◆ trkEtaMax_

float L1TrackFastJetProducer::trkEtaMax_

private

Definition at line 71 of file L1TrackFastJetProducer.cc.

Referenced by L1TrackFastJetProducer(), and produce().

◆ trkNPSStubMin_

int L1TrackFastJetProducer::trkNPSStubMin_

private

Definition at line 73 of file L1TrackFastJetProducer.cc.

Referenced by L1TrackFastJetProducer(), and produce().

◆ trkNStubMin_

int L1TrackFastJetProducer::trkNStubMin_

private

Definition at line 72 of file L1TrackFastJetProducer.cc.

Referenced by L1TrackFastJetProducer(), and produce().

◆ trkPtMin_

float L1TrackFastJetProducer::trkPtMin_

private

Definition at line 70 of file L1TrackFastJetProducer.cc.

Referenced by L1TrackFastJetProducer(), and produce().

◆ trkPtTightChi2_

float L1TrackFastJetProducer::trkPtTightChi2_

private

Definition at line 77 of file L1TrackFastJetProducer.cc.

Referenced by L1TrackFastJetProducer(), and produce().

◆ trkZMax_

float L1TrackFastJetProducer::trkZMax_

private

Definition at line 67 of file L1TrackFastJetProducer.cc.

Referenced by L1TrackFastJetProducer(), and produce().

◆ tTopoToken_

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

private

Definition at line 83 of file L1TrackFastJetProducer.cc.

Referenced by produce().

Public Types

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ L1TTTrackCollectionType

◆ L1TTTrackType

Constructor & Destructor Documentation

◆ L1TrackFastJetProducer()

◆ ~L1TrackFastJetProducer()

Member Function Documentation

◆ fillDescriptions()

◆ produce()

Member Data Documentation

◆ coneSize_

◆ deltaZ0Cut_

◆ displaced_

◆ doTightChi2_

◆ pvToken_

◆ trackToken_

◆ trkBendChi2Max_

◆ trkChi2dofMax_

◆ trkChi2dofTightChi2_

◆ trkEtaMax_

◆ trkNPSStubMin_

◆ trkNStubMin_

◆ trkPtMin_

◆ trkPtTightChi2_

◆ trkZMax_

◆ tTopoToken_