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
const double	coneSize_

const double	deltaZ0Cut_

const bool	displaced_

const bool	doTightChi2_

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

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

bool	selectTrkMatchGenLoose_

bool	selectTrkMatchGenOrPU_

bool	selectTrkMatchGenTight_

edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	tGeomToken_

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

const double	trkBendChi2Max_

const float	trkChi2dofMax_

const float	trkChi2dofTightChi2_

const float	trkEtaMax_

const int	trkNPSStubMin_

const int	trkNStubMin_

const float	trkPtMin_

const float	trkPtTightChi2_

const float	trkZMax_

edm::ESGetToken< TrackerTopology, TrackerTopologyRcd >	tTopoToken_

Detailed Description

Definition at line 56 of file L1FastTrackingJetProducer.cc.

Member Typedef Documentation

◆ L1TTTrackCollectionType

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

Definition at line 59 of file L1FastTrackingJetProducer.cc.

◆ L1TTTrackType

typedef TTTrack<Ref_Phase2TrackerDigi_> L1FastTrackingJetProducer::L1TTTrackType

Definition at line 58 of file L1FastTrackingJetProducer.cc.

Constructor & Destructor Documentation

◆ L1FastTrackingJetProducer()

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

explicit

Definition at line 97 of file L1FastTrackingJetProducer.cc.

References displaced_.

     : 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")),
       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("", ""))),
       tGeomToken_(esConsumes<TrackerGeometry, TrackerDigiGeometryRecord>(edm::ESInputTag("", ""))) {
   if (displaced_)
     produces<TkJetCollection>("L1FastTrackingJetsExtended");
   else
     produces<TkJetCollection>("L1FastTrackingJets");
 }

◆ ~L1FastTrackingJetProducer()

L1FastTrackingJetProducer::~L1FastTrackingJetProducer ( )

override

Definition at line 128 of file L1FastTrackingJetProducer.cc.

128 {}

Member Function Documentation

◆ beginJob()

void L1FastTrackingJetProducer::beginJob ( void )

privatevirtual

Definition at line 251 of file L1FastTrackingJetProducer.cc.

251 {}

◆ endJob()

void L1FastTrackingJetProducer::endJob ( void )

privatevirtual

Definition at line 253 of file L1FastTrackingJetProducer.cc.

253 {}

◆ fillDescriptions()

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

static

Definition at line 255 of file L1FastTrackingJetProducer.cc.

References edm::ConfigurationDescriptions::add(), submitPVResolutionJobs::desc, HLT_2022v12_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);
   }*/
 }

◆ produce()

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

overrideprivate

Definition at line 131 of file L1FastTrackingJetProducer.cc.

References funct::abs(), coneSize_, callgraph::cs, deltaZ0Cut_, displaced_, doTightChi2_, genToken_, edm::EventSetup::getData(), mps_fire::i, iEvent, createfilelist::int, L1FastTrackingJetProducer_cfi::L1FastTrackingJets, TTTrack< T >::momentum(), eostools::move(), PV3DBase< T, PVType, FrameType >::perp(), pvToken_, multPhiCorr_741_25nsDY_cfi::px, multPhiCorr_741_25nsDY_cfi::py, selectTrkMatchGenLoose_, selectTrkMatchGenOrPU_, selectTrkMatchGenTight_, 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().

                                                                                      {
   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_);
   //const TrackerGeometry& tGeom = iSetup.getData(tGeomToken_);
 
   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");
 }