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 52 of file L1TrackFastJetProducer.cc.

Member Typedef Documentation

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

Definition at line 55 of file L1TrackFastJetProducer.cc.

typedef TTTrack<Ref_Phase2TrackerDigi_> L1TrackFastJetProducer::L1TTTrackType

Definition at line 54 of file L1TrackFastJetProducer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 88 of file L1TrackFastJetProducer.cc.

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

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

L1TrackFastJetProducer::~L1TrackFastJetProducer ( )

override

Definition at line 113 of file L1TrackFastJetProducer.cc.

113 {}

Member Function Documentation

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

static

Definition at line 221 of file L1TrackFastJetProducer.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
   edm::ParameterSetDescription desc;
   desc.setUnknown();
   descriptions.addDefault(desc);
 }

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

overrideprivate

Definition at line 116 of file L1TrackFastJetProducer.cc.

References coneSize_, fwrapper::cs, deltaZ0Cut_, displaced_, doTightChi2_, edm::Event::getByToken(), edm::EventSetup::getData(), mps_fire::i, fwrapper::jet_def, L1TrackFastJetProducer_cfi::L1TrackFastJets, eostools::move(), edm::Event::put(), pvToken_, StripSubdetector::TID, TrackerTopology::tidRing(), StripSubdetector::TOB, TrackerTopology::tobLayer(), align::Tracker, trackToken_, trkBendChi2Max_, trkChi2dofMax_, trkChi2dofTightChi2_, trkEtaMax_, trkNPSStubMin_, trkNStubMin_, trkPtMin_, trkPtTightChi2_, trkZMax_, and tTopoToken_.

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