L1TrackFastJetProducer Class Reference

Inheritance diagram for L1TrackFastJetProducer:

Public Types
typedef std::vector< L1TTTrackType >	L1TTTrackCollectionType
typedef edm::RefVector< L1TTTrackCollectionType >	L1TTTrackRefCollectionType
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 bool	displaced_
const EDGetTokenT< L1TTTrackRefCollectionType >	trackToken_

Detailed Description

Definition at line 54 of file L1TrackFastJetProducer.cc.

Member Typedef Documentation

L1TTTrackCollectionType

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

Definition at line 57 of file L1TrackFastJetProducer.cc.

L1TTTrackRefCollectionType

typedef edm::RefVector<L1TTTrackCollectionType> L1TrackFastJetProducer::L1TTTrackRefCollectionType

Definition at line 58 of file L1TrackFastJetProducer.cc.

L1TTTrackType

typedef TTTrack<Ref_Phase2TrackerDigi_> L1TrackFastJetProducer::L1TTTrackType

Definition at line 56 of file L1TrackFastJetProducer.cc.

Constructor & Destructor Documentation

L1TrackFastJetProducer()

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

explicit

Definition at line 76 of file L1TrackFastJetProducer.cc.

References displaced_.

    : coneSize_((float)iConfig.getParameter<double>("coneSize")),
      displaced_(iConfig.getParameter<bool>("displaced")),
      trackToken_(consumes<L1TTTrackRefCollectionType>(iConfig.getParameter<InputTag>("L1TrackInputTag"))) {
  if (displaced_)
    produces<TkJetCollection>("L1TrackFastJetsExtended");
  else
    produces<TkJetCollection>("L1TrackFastJets");
}

~L1TrackFastJetProducer()

L1TrackFastJetProducer::~L1TrackFastJetProducer ( )

override

Definition at line 87 of file L1TrackFastJetProducer.cc.

{}

Member Function Documentation

fillDescriptions()

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

static

Definition at line 142 of file L1TrackFastJetProducer.cc.

References submitPVResolutionJobs::desc, and ProducerED_cfi::InputTag.

                                                                                        {
  //The following says we do not know what parameters are allowed so do no validation
  edm::ParameterSetDescription desc;
  desc.add<edm::InputTag>(
      "L1PVertexInputTag",
      edm::InputTag("l1tTrackVertexAssociationProducerForJets", "Level1TTTracksSelectedAssociated"));
  desc.add<double>("coneSize", 0.5);
  desc.add<bool>("displaced", false);
}

produce()

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

overrideprivate

Definition at line 90 of file L1TrackFastJetProducer.cc.

References coneSize_, callgraph::cs, displaced_, mps_fire::i, iEvent, PV3DBase< T, PVType, FrameType >::mag(), TTTrack< T >::momentum(), eostools::move(), PV3DBase< T, PVType, FrameType >::perp(), multPhiCorr_741_25nsDY_cfi::px, multPhiCorr_741_25nsDY_cfi::py, edm::RefVector< C, T, F >::size(), trackToken_, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), PV3DBase< T, PVType, FrameType >::z(), and TTTrack< T >::z0().

                                                                                  {
  std::unique_ptr<TkJetCollection> L1TrackFastJets(new TkJetCollection);

  // L1 tracks
  edm::Handle<L1TTTrackRefCollectionType> TTTrackHandle;
  iEvent.getByToken(trackToken_, TTTrackHandle);

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

  for (unsigned int this_l1track = 0; this_l1track < TTTrackHandle->size(); this_l1track++) {
    edm::Ptr<L1TTTrackType> iterL1Track(TTTrackHandle, this_l1track);

    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);  // 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 69 of file L1TrackFastJetProducer.cc.

Referenced by produce().

displaced_

const bool L1TrackFastJetProducer::displaced_

private

Definition at line 70 of file L1TrackFastJetProducer.cc.

Referenced by L1TrackFastJetProducer(), and produce().

trackToken_

const EDGetTokenT<L1TTTrackRefCollectionType> L1TrackFastJetProducer::trackToken_

private

Definition at line 72 of file L1TrackFastJetProducer.cc.

Referenced by produce().