L1TrackFastJetProducer Class Reference

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
const double	coneSize_
const double	deltaZ0Cut_
const bool	displaced_
const bool	doTightChi2_
edm::EDGetTokenT< VertexCollection >	pvToken_
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 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.

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")),
      trackToken_(consumes<std::vector<TTTrack<Ref_Phase2TrackerDigi_> > >(
          iConfig.getParameter<edm::InputTag>("L1TrackInputTag"))),
      pvToken_(consumes<VertexCollection>(iConfig.getParameter<edm::InputTag>("L1PrimaryVertexTag"))),
      tTopoToken_(esConsumes<TrackerTopology, TrackerTopologyRcd>(edm::ESInputTag("", ""))) {
  if (displaced_)
    produces<TkJetCollection>("L1TrackFastJetsExtended");
  else
    produces<TkJetCollection>("L1TrackFastJets");
}

~L1TrackFastJetProducer()

L1TrackFastJetProducer::~L1TrackFastJetProducer ( )

override

Definition at line 112 of file L1TrackFastJetProducer.cc.

{}

Member Function Documentation

fillDescriptions()

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

static

Definition at line 220 of file L1TrackFastJetProducer.cc.

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

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

produce()

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

overrideprivate

Definition at line 115 of file L1TrackFastJetProducer.cc.

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

                                                                                  {
  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<l1t::VertexCollection> L1PrimaryVertexHandle;
  iEvent.getByToken(pvToken_, L1PrimaryVertexHandle);

  fastjet::JetDefinition jet_def(fastjet::antikt_algorithm, coneSize_);
  std::vector<fastjet::PseudoJet> JetInputs;

  float recoVtx = L1PrimaryVertexHandle->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;

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

Member Data Documentation

coneSize_

const double L1TrackFastJetProducer::coneSize_

private

Definition at line 75 of file L1TrackFastJetProducer.cc.

Referenced by produce().

deltaZ0Cut_

const double L1TrackFastJetProducer::deltaZ0Cut_

private

Definition at line 74 of file L1TrackFastJetProducer.cc.

Referenced by produce().

displaced_

const bool L1TrackFastJetProducer::displaced_

private

Definition at line 79 of file L1TrackFastJetProducer.cc.

Referenced by L1TrackFastJetProducer(), and produce().

doTightChi2_

const bool L1TrackFastJetProducer::doTightChi2_

private

Definition at line 76 of file L1TrackFastJetProducer.cc.

Referenced by produce().

pvToken_

edm::EDGetTokenT<VertexCollection> 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_

const double L1TrackFastJetProducer::trkBendChi2Max_

private

Definition at line 69 of file L1TrackFastJetProducer.cc.

Referenced by produce().

trkChi2dofMax_

const float L1TrackFastJetProducer::trkChi2dofMax_

private

Definition at line 68 of file L1TrackFastJetProducer.cc.

Referenced by produce().

trkChi2dofTightChi2_

const float L1TrackFastJetProducer::trkChi2dofTightChi2_

private

Definition at line 78 of file L1TrackFastJetProducer.cc.

Referenced by produce().

trkEtaMax_

const float L1TrackFastJetProducer::trkEtaMax_

private

Definition at line 71 of file L1TrackFastJetProducer.cc.

Referenced by produce().

trkNPSStubMin_

const int L1TrackFastJetProducer::trkNPSStubMin_

private

Definition at line 73 of file L1TrackFastJetProducer.cc.

Referenced by produce().

trkNStubMin_

const int L1TrackFastJetProducer::trkNStubMin_

private

Definition at line 72 of file L1TrackFastJetProducer.cc.

Referenced by produce().

trkPtMin_

const float L1TrackFastJetProducer::trkPtMin_

private

Definition at line 70 of file L1TrackFastJetProducer.cc.

Referenced by produce().

trkPtTightChi2_

const float L1TrackFastJetProducer::trkPtTightChi2_

private

Definition at line 77 of file L1TrackFastJetProducer.cc.

Referenced by produce().

trkZMax_

const float L1TrackFastJetProducer::trkZMax_

private

Definition at line 67 of file L1TrackFastJetProducer.cc.

Referenced by produce().

tTopoToken_

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

private

Definition at line 83 of file L1TrackFastJetProducer.cc.

Referenced by produce().