Inheritance diagram for L1FastTrackingJetProducer:

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

	~L1FastTrackingJetProducer () 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
virtual void	beginJob ()

virtual void	endJob ()

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

Private Attributes
double	coneSize_

double	deltaZ0Cut_

bool	displaced_

bool	doTightChi2_

const edm::EDGetTokenT < TTTrackAssociationMap < Ref_Phase2TrackerDigi_ > >	genToken_

edm::EDGetTokenT< std::vector < l1t::Vertex > >	pvToken_

bool	selectTrkMatchGenLoose_

bool	selectTrkMatchGenOrPU_

bool	selectTrkMatchGenTight_

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 57 of file L1FastTrackingJetProducer.cc.

Member Typedef Documentation

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

Definition at line 60 of file L1FastTrackingJetProducer.cc.

typedef TTTrack<Ref_Phase2TrackerDigi_> L1FastTrackingJetProducer::L1TTTrackType

Definition at line 59 of file L1FastTrackingJetProducer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 97 of file L1FastTrackingJetProducer.cc.

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

     : trackToken_(consumes<std::vector<TTTrack<Ref_Phase2TrackerDigi_>>>(
           iConfig.getParameter<edm::InputTag>("L1TrackInputTag"))),
       pvToken_(consumes<std::vector<l1t::Vertex>>(iConfig.getParameter<edm::InputTag>("L1PrimaryVertexTag"))),
       genToken_(
           consumes<TTTrackAssociationMap<Ref_Phase2TrackerDigi_>>(iConfig.getParameter<edm::InputTag>("GenInfo"))),
       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");
   selectTrkMatchGenTight_ = iConfig.getParameter<bool>("selectTrkMatchGenTight");
   selectTrkMatchGenLoose_ = iConfig.getParameter<bool>("selectTrkMatchGenLoose");
   selectTrkMatchGenOrPU_ = iConfig.getParameter<bool>("selectTrkMatchGenOrPU");
   if (displaced_)
     produces<TkJetCollection>("L1FastTrackingJetsExtended");
   else
     produces<TkJetCollection>("L1FastTrackingJets");
 }

L1FastTrackingJetProducer::~L1FastTrackingJetProducer ( )

override

Definition at line 127 of file L1FastTrackingJetProducer.cc.

127 {}

Member Function Documentation

void L1FastTrackingJetProducer::beginJob ( void )

privatevirtual

Definition at line 249 of file L1FastTrackingJetProducer.cc.

249 {}

void L1FastTrackingJetProducer::endJob ( void )

privatevirtual

Definition at line 251 of file L1FastTrackingJetProducer.cc.

251 {}

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

static

Definition at line 253 of file L1FastTrackingJetProducer.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), submitPVResolutionJobs::desc, HLT_FULL_cff::InputTag, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                            {
   {
     // L1FastTrackingJets
     edm::ParameterSetDescription desc;
     desc.add<edm::InputTag>("L1TrackInputTag", edm::InputTag("TTTracksFromTrackletEmulation", "Level1TTTracks"));
     desc.add<std::string>("L1PrimaryVertexTag", "l1vertices");
     desc.add<edm::InputTag>("GenInfo", edm::InputTag("TTTrackAssociatorFromPixelDigis", "Level1TTTracks"));
     desc.add<double>("trk_zMax", 15.0);
     desc.add<double>("trk_chi2dofMax", 10.0);
     desc.add<double>("trk_bendChi2Max", 2.2);
     desc.add<double>("trk_ptMin", 2.0);
     desc.add<double>("trk_etaMax", 2.5);
     desc.add<int>("trk_nStubMin", 4);
     desc.add<int>("trk_nPSStubMin", -1);
     desc.add<double>("deltaZ0Cut", 0.5);
     desc.add<bool>("doTightChi2", true);
     desc.add<double>("trk_ptTightChi2", 20.0);
     desc.add<double>("trk_chi2dofTightChi2", 5.0);
     desc.add<double>("coneSize", 0.4);
     desc.add<bool>("displaced", false);
     desc.add<bool>("selectTrkMatchGenTight", true);
     desc.add<bool>("selectTrkMatchGenLoose", false);
     desc.add<bool>("selectTrkMatchGenOrPU", false);
     descriptions.add("L1FastTrackingJets", desc);
     // or use the following to generate the label from the module's C++ type
     //descriptions.addWithDefaultLabel(desc);
   }
   /*{                                                                                                                           
     // L1FastTrackingJetsExtended                                                                                             
     desc.add<double>("trk_bendChi2Max", 2.4);
     desc.add<double>("trk_ptMin", 3.0);
     desc.add<double>("trk_etaMax", 2.5);
     desc.add<int>("trk_nStubMin", 4);
     desc.add<int>("trk_nPSStubMin", -1);
     desc.add<double>("deltaZ0Cut", 3.0);
     desc.add<bool>("doTightChi2", true);
     desc.add<double>("trk_ptTightChi2", 20.0);
     desc.add<double>("trk_chi2dofTightChi2", 5.0);
     desc.add<double>("coneSize", 0.4);
     desc.add<bool>("displaced", true);
     desc.add<bool>("selectTrkMatchGenTight", true);
     desc.add<bool>("selectTrkMatchGenLoose", false);
     desc.add<bool>("selectTrkMatchGenOrPU", false);
     descriptions.add("L1FastTrackingJetsExtended", desc);
     // or use the following to generate the label from the module's C++ type
     //descriptions.addWithDefaultLabel(desc);
   }*/
 }

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

overrideprivate

Definition at line 130 of file L1FastTrackingJetProducer.cc.

References funct::abs(), coneSize_, fwrapper::cs, deltaZ0Cut_, displaced_, doTightChi2_, genToken_, edm::Event::getByToken(), edm::EventSetup::getData(), mps_fire::i, edm::EventBase::isRealData(), fwrapper::jet_def, L1FastTrackingJetProducer_cfi::L1FastTrackingJets, eostools::move(), edm::Event::put(), pvToken_, selectTrkMatchGenLoose_, selectTrkMatchGenOrPU_, selectTrkMatchGenTight_, 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> L1FastTrackingJets(new TkJetCollection);
 
   // L1 tracks
   edm::Handle<std::vector<TTTrack<Ref_Phase2TrackerDigi_>>> TTTrackHandle;
   iEvent.getByToken(trackToken_, TTTrackHandle);
   std::vector<TTTrack<Ref_Phase2TrackerDigi_>>::const_iterator iterL1Track;
 
   // Gen
   edm::Handle<TTTrackAssociationMap<Ref_Phase2TrackerDigi_>> MCTrkAssociation;
   if (!iEvent.isRealData())
     iEvent.getByToken(genToken_, MCTrkAssociation);
 
   // Tracker Topology
   const TrackerTopology& tTopo = iSetup.getData(tTopoToken_);
 
   edm::Handle<std::vector<l1t::Vertex>> L1VertexHandle;
   iEvent.getByToken(pvToken_, L1VertexHandle);
 
   fastjet::JetDefinition jet_def(fastjet::antikt_algorithm, coneSize_);
   std::vector<fastjet::PseudoJet> JetInputs;
 
   float recoVtx = L1VertexHandle->begin()->z0();
   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 (std::abs(trk_z0) > trkZMax_)
       continue;
     if (std::abs(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 (std::abs(recoVtx - trk_z0) > deltaZ0Cut_)
       continue;
     if (!iEvent.isRealData()) {
       edm::Ptr<TTTrack<Ref_Phase2TrackerDigi_>> trk_ptr(TTTrackHandle, this_l1track);
 
       if (!(MCTrkAssociation->isGenuine(trk_ptr)) && selectTrkMatchGenTight_)
         continue;
       if (!(MCTrkAssociation->isLooselyGenuine(trk_ptr) || MCTrkAssociation->isGenuine(trk_ptr)) &&
           selectTrkMatchGenLoose_)
         continue;
       if (!(MCTrkAssociation->isLooselyGenuine(trk_ptr) || MCTrkAssociation->isGenuine(trk_ptr) ||
             MCTrkAssociation->isCombinatoric(trk_ptr)) &&
           selectTrkMatchGenOrPU_)
         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);
     L1FastTrackingJets->push_back(trkJet);
   }  //end loop over Jet Outputs
 
   if (displaced_)
     iEvent.put(std::move(L1FastTrackingJets), "L1FastTrackingJetsExtended");
   else
     iEvent.put(std::move(L1FastTrackingJets), "L1FastTrackingJets");
 }