Inheritance diagram for l1tpf::PFTrackProducerFromL1Tracks:

Public Member Functions
	PFTrackProducerFromL1Tracks (const edm::ParameterSet &)

	~PFTrackProducerFromL1Tracks () 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

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

Private Attributes
edm::ESGetToken< MagneticField, IdealMagneticFieldRecord >	BFieldTag_

edm::ESWatcher< IdealMagneticFieldRecord >	BFieldWatcher_

float	fBz_

int	nParam_

std::vector< StringCutObjectSelector< l1t::PFTrack::L1TTTrackType > >	qualityBitSetters_

bool	redigitizeTrackWord_

l1tpf::ParametricResolution	resolCalo_

l1tpf::ParametricResolution	resolTrk_

edm::EDGetTokenT< std::vector< l1t::PFTrack::L1TTTrackType > >	TrackTag_

Additional Inherited Members
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

Detailed Description

Definition at line 18 of file PFTrackProducerFromL1Tracks.cc.

Constructor & Destructor Documentation

◆ PFTrackProducerFromL1Tracks()

l1tpf::PFTrackProducerFromL1Tracks::PFTrackProducerFromL1Tracks ( const edm::ParameterSet & iConfig )

explicit

Definition at line 38 of file PFTrackProducerFromL1Tracks.cc.

     : TrackTag_(consumes<std::vector<l1t::PFTrack::L1TTTrackType>>(iConfig.getParameter<edm::InputTag>("L1TrackTag"))),
       BFieldTag_{esConsumes<MagneticField, IdealMagneticFieldRecord>()},
       nParam_(iConfig.getParameter<unsigned int>("nParam")),
       resolCalo_(iConfig.getParameter<edm::ParameterSet>("resolCalo")),
       resolTrk_(iConfig.getParameter<edm::ParameterSet>("resolTrack")),
       redigitizeTrackWord_(iConfig.getParameter<bool>("redigitizeTrackWord")) {
   for (const auto &cut : iConfig.getParameter<std::vector<std::string>>("qualityBits")) {
     qualityBitSetters_.emplace_back(cut);
   }
   produces<l1t::PFTrackCollection>();
 }

◆ ~PFTrackProducerFromL1Tracks()

l1tpf::PFTrackProducerFromL1Tracks::~PFTrackProducerFromL1Tracks ( )

inlineoverride

Definition at line 21 of file PFTrackProducerFromL1Tracks.cc.

21 {}

Member Function Documentation

◆ produce()

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

overrideprivate

Definition at line 51 of file PFTrackProducerFromL1Tracks.cc.

                                                                                             {
   if (BFieldWatcher_.check(iSetup)) {
     fBz_ = iSetup.getData(BFieldTag_).inTesla(GlobalPoint(0, 0, 0)).z();
   }
 
   std::unique_ptr<l1t::PFTrackCollection> out(new l1t::PFTrackCollection());
 
   // https://github.com/skinnari/cmssw/blob/80c19f1b721325c3a02ee0482f72fb974a4c3bf7/L1Trigger/TrackFindingTracklet/test/L1TrackNtupleMaker.cc
   edm::Handle<std::vector<l1t::PFTrack::L1TTTrackType>> htracks;
   iEvent.getByToken(TrackTag_, htracks);
   const auto &tracks = *htracks;
 
   for (unsigned int i = 0, n = tracks.size(); i < n; ++i) {
     const auto &tk = tracks[i];
 
     float pt = tk.momentum().perp();
     float eta = tk.momentum().eta();
     float phi = tk.momentum().phi();
     float z0 = tk.POCA().z();  //cm
     int charge = tk.rInv() > 0 ? +1 : -1;
 
     reco::Candidate::PolarLorentzVector p4p(pt, eta, phi, 0.137);  // pion mass
     reco::Particle::LorentzVector p4(p4p.X(), p4p.Y(), p4p.Z(), p4p.E());
     reco::Particle::Point vtx(tk.POCA().x(), tk.POCA().y(), z0);
 
     auto caloetaphi = l1tpf::propagateToCalo(p4, math::XYZTLorentzVector(0., 0., z0, 0.), charge, fBz_);
 
     float trkErr = resolTrk_(pt, std::abs(eta));
     float caloErr = resolCalo_(pt, std::abs(eta));
     int quality = 0;
     for (int i = 0, n = qualityBitSetters_.size(), bit = 1; i < n; ++i, bit <<= 1) {
       if (qualityBitSetters_[i](tk))
         quality += bit;
     }
     out->emplace_back(charge,
                       p4,
                       vtx,
                       l1t::PFTrack::TrackRef(htracks, i),
                       nParam_,
                       caloetaphi.first,
                       caloetaphi.second,
                       trkErr,
                       caloErr,
                       quality);
 
     if (redigitizeTrackWord_) {
       // Force re-digitization if an old TTrack object is read from an EDM file, and update the quaility bit for now
       l1t::PFTrack::L1TTTrackType trackCopy = tk;
       trackCopy.setTrackWordBits();  // important
       TTTrack_TrackWord &tw = out->back().trackWord();
       tw.setTrackWord(trackCopy.getValidWord(),
                       trackCopy.getRinvWord(),
                       trackCopy.getPhiWord(),
                       trackCopy.getTanlWord(),
                       trackCopy.getZ0Word(),
                       trackCopy.getD0Word(),
                       trackCopy.getChi2RPhiWord(),
                       trackCopy.getChi2RZWord(),
                       trackCopy.getBendChi2Word(),
                       trackCopy.getHitPatternWord(),
                       trackCopy.getMVAQualityWord(),
                       ap_uint<TTTrack_TrackWord::kMVAOtherSize>(quality));
     }
   }
   iEvent.put(std::move(out));
 }