L1FastTrackingJetProducer Class Reference

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

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

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

~L1FastTrackingJetProducer()

L1FastTrackingJetProducer::~L1FastTrackingJetProducer ( )

override

Definition at line 126 of file L1FastTrackingJetProducer.cc.

{}

Member Function Documentation

beginJob()

void L1FastTrackingJetProducer::beginJob ( void  )

privatevirtual

Definition at line 248 of file L1FastTrackingJetProducer.cc.

{}

endJob()

void L1FastTrackingJetProducer::endJob ( void  )

privatevirtual

Definition at line 250 of file L1FastTrackingJetProducer.cc.

{}

fillDescriptions()

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

static

Definition at line 252 of file L1FastTrackingJetProducer.cc.

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

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

produce()

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

overrideprivate

Definition at line 129 of file L1FastTrackingJetProducer.cc.

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

References funct::abs(), coneSize_, fwrapper::cs, deltaZ0Cut_, displaced_, doTightChi2_, genToken_, edm::EventSetup::getData(), mps_fire::i, iEvent, createfilelist::int, fwrapper::jet_def, 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().

Member Data Documentation

coneSize_

double L1FastTrackingJetProducer::coneSize_

private

Definition at line 80 of file L1FastTrackingJetProducer.cc.

Referenced by L1FastTrackingJetProducer(), and produce().

deltaZ0Cut_

double L1FastTrackingJetProducer::deltaZ0Cut_

private

Definition at line 79 of file L1FastTrackingJetProducer.cc.

Referenced by L1FastTrackingJetProducer(), and produce().

displaced_

bool L1FastTrackingJetProducer::displaced_

private

Definition at line 84 of file L1FastTrackingJetProducer.cc.

Referenced by L1FastTrackingJetProducer(), and produce().

doTightChi2_

bool L1FastTrackingJetProducer::doTightChi2_

private

Definition at line 81 of file L1FastTrackingJetProducer.cc.

Referenced by L1FastTrackingJetProducer(), and produce().

genToken_

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

private

Definition at line 91 of file L1FastTrackingJetProducer.cc.

Referenced by produce().

pvToken_

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

private

Definition at line 90 of file L1FastTrackingJetProducer.cc.

Referenced by produce().

selectTrkMatchGenLoose_

bool L1FastTrackingJetProducer::selectTrkMatchGenLoose_

private

Definition at line 86 of file L1FastTrackingJetProducer.cc.

Referenced by L1FastTrackingJetProducer(), and produce().

selectTrkMatchGenOrPU_

bool L1FastTrackingJetProducer::selectTrkMatchGenOrPU_

private

Definition at line 87 of file L1FastTrackingJetProducer.cc.

Referenced by L1FastTrackingJetProducer(), and produce().

selectTrkMatchGenTight_

bool L1FastTrackingJetProducer::selectTrkMatchGenTight_

private

Definition at line 85 of file L1FastTrackingJetProducer.cc.

Referenced by L1FastTrackingJetProducer(), and produce().

trackToken_

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

private

Definition at line 89 of file L1FastTrackingJetProducer.cc.

Referenced by produce().

trkBendChi2Max_

double L1FastTrackingJetProducer::trkBendChi2Max_

private

Definition at line 74 of file L1FastTrackingJetProducer.cc.

Referenced by L1FastTrackingJetProducer(), and produce().

trkChi2dofMax_

float L1FastTrackingJetProducer::trkChi2dofMax_

private

Definition at line 73 of file L1FastTrackingJetProducer.cc.

Referenced by L1FastTrackingJetProducer(), and produce().

trkChi2dofTightChi2_

float L1FastTrackingJetProducer::trkChi2dofTightChi2_

private

Definition at line 83 of file L1FastTrackingJetProducer.cc.

Referenced by L1FastTrackingJetProducer(), and produce().

trkEtaMax_

float L1FastTrackingJetProducer::trkEtaMax_

private

Definition at line 76 of file L1FastTrackingJetProducer.cc.

Referenced by L1FastTrackingJetProducer(), and produce().

trkNPSStubMin_

int L1FastTrackingJetProducer::trkNPSStubMin_

private

Definition at line 78 of file L1FastTrackingJetProducer.cc.

Referenced by L1FastTrackingJetProducer(), and produce().

trkNStubMin_

int L1FastTrackingJetProducer::trkNStubMin_

private

Definition at line 77 of file L1FastTrackingJetProducer.cc.

Referenced by L1FastTrackingJetProducer(), and produce().

trkPtMin_

float L1FastTrackingJetProducer::trkPtMin_

private

Definition at line 75 of file L1FastTrackingJetProducer.cc.

Referenced by L1FastTrackingJetProducer(), and produce().

trkPtTightChi2_

float L1FastTrackingJetProducer::trkPtTightChi2_

private

Definition at line 82 of file L1FastTrackingJetProducer.cc.

Referenced by L1FastTrackingJetProducer(), and produce().

trkZMax_

float L1FastTrackingJetProducer::trkZMax_

private

Definition at line 72 of file L1FastTrackingJetProducer.cc.

Referenced by L1FastTrackingJetProducer(), and produce().

tTopoToken_

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

private

Definition at line 92 of file L1FastTrackingJetProducer.cc.

Referenced by produce().