Inheritance diagram for Run3ScoutingParticleToRecoPFCandidateProducer:

Public Member Functions
void	beginStream (edm::StreamID) override

void	clearVars ()

reco::PFCandidate	createPFCand (Run3ScoutingParticle scoutingparticle)

void	createPFCandidates (edm::Handle< std::vector< Run3ScoutingParticle >> scoutingparticleHandle, std::unique_ptr< reco::PFCandidateCollection > &pfcands)

void	createPFCandidatesSK (edm::Handle< std::vector< Run3ScoutingParticle >> scoutingparticleHandle, std::unique_ptr< reco::PFCandidateCollection > &pfcands)

void	endStream () override

void	produce (edm::Event &iEvent, edm::EventSetup const &setup) override

	Run3ScoutingParticleToRecoPFCandidateProducer (const edm::ParameterSet &)

	~Run3ScoutingParticleToRecoPFCandidateProducer () 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 Attributes
std::vector< float >	dxy_

std::vector< float >	dxysig_

std::vector< float >	dz_

std::vector< float >	dzsig_

const edm::EDGetTokenT< std::vector< Run3ScoutingParticle > >	input_scoutingparticle_token_

std::vector< int >	lostInnerHits_

std::vector< float >	normchi2_

const edm::ESGetToken< HepPDT::ParticleDataTable, edm::DefaultRecord >	particletable_token_

const HepPDT::ParticleDataTable *	pdTable_

std::vector< int >	quality_

std::vector< float >	trkEta_

std::vector< float >	trkPhi_

std::vector< float >	trkPt_

bool	use_CHS_

bool	use_softKiller_

std::vector< int8_t >	vertexIndex_

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 19 of file Run3ScoutingParticleToRecoPFCandidateProducer.cc.

Constructor & Destructor Documentation

◆ Run3ScoutingParticleToRecoPFCandidateProducer()

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

explicit

Definition at line 59 of file Run3ScoutingParticleToRecoPFCandidateProducer.cc.

     : input_scoutingparticle_token_(consumes(iConfig.getParameter<edm::InputTag>("scoutingparticle"))),
       particletable_token_(esConsumes<HepPDT::ParticleDataTable, edm::DefaultRecord>()),
       use_softKiller_(iConfig.getParameter<bool>("softKiller")),
       use_CHS_(iConfig.getParameter<bool>("CHS")) {
   //register products
   produces<reco::PFCandidateCollection>();
   produces<edm::ValueMap<int>>("vertexIndex");
   produces<edm::ValueMap<float>>("normchi2");
   produces<edm::ValueMap<float>>("dz");
   produces<edm::ValueMap<float>>("dxy");
   produces<edm::ValueMap<float>>("dzsig");
   produces<edm::ValueMap<float>>("dxysig");
   produces<edm::ValueMap<int>>("lostInnerHits");
   produces<edm::ValueMap<int>>("quality");
   produces<edm::ValueMap<float>>("trkPt");
   produces<edm::ValueMap<float>>("trkEta");
   produces<edm::ValueMap<float>>("trkPhi");
 }

◆ ~Run3ScoutingParticleToRecoPFCandidateProducer()

Run3ScoutingParticleToRecoPFCandidateProducer::~Run3ScoutingParticleToRecoPFCandidateProducer ( )

overridedefault

Member Function Documentation

◆ beginStream()

void Run3ScoutingParticleToRecoPFCandidateProducer::beginStream ( edm::StreamID )

inlineoverride

Definition at line 25 of file Run3ScoutingParticleToRecoPFCandidateProducer.cc.

25 {}

◆ clearVars()

void Run3ScoutingParticleToRecoPFCandidateProducer::clearVars ( )

Definition at line 262 of file Run3ScoutingParticleToRecoPFCandidateProducer.cc.

References dxy_, dxysig_, dz_, dzsig_, lostInnerHits_, normchi2_, quality_, trkEta_, trkPhi_, trkPt_, and vertexIndex_.

Referenced by produce().

                                                               {
   vertexIndex_.clear();
   normchi2_.clear();
   dz_.clear();
   dxy_.clear();
   dzsig_.clear();
   dxysig_.clear();
   lostInnerHits_.clear();
   quality_.clear();
   trkPt_.clear();
   trkEta_.clear();
   trkPhi_.clear();
 }

◆ createPFCand()

reco::PFCandidate Run3ScoutingParticleToRecoPFCandidateProducer::createPFCand ( Run3ScoutingParticle scoutingparticle )

Definition at line 82 of file Run3ScoutingParticleToRecoPFCandidateProducer.cc.

References funct::cos(), dxy_, dxysig_, dz_, dzsig_, hcalRecHitTable_cff::energy, LogDebug, lostInnerHits_, visualization-live-secondInstance_cfg::m, normchi2_, or, AlCaHLTBitMon_ParallelJobs::p, LHEGenericFilter_cfi::ParticleID, pdTable_, pfDeepBoostedJetPreprocessParams_cfi::pfcand, pfElectronTranslator_cfi::PFCandidate, multPhiCorr_741_25nsDY_cfi::px, multPhiCorr_741_25nsDY_cfi::py, submitPVResolutionJobs::q, quality_, HLT_2024v13_cff::relativeTrackVars, run3scouting_cff::scoutingparticle, funct::sin(), mathSSE::sqrt(), trkEta_, trkPhi_, trkPt_, and vertexIndex_.

Referenced by createPFCandidates(), and createPFCandidatesSK().

                                                                                                                  {
   auto m = pdTable_->particle(HepPDT::ParticleID(scoutingparticle.pdgId())) != nullptr
                ? pdTable_->particle(HepPDT::ParticleID(scoutingparticle.pdgId()))->mass()
                : -99.f;
   auto q = pdTable_->particle(HepPDT::ParticleID(scoutingparticle.pdgId())) != nullptr
                ? pdTable_->particle(HepPDT::ParticleID(scoutingparticle.pdgId()))->charge()
                : -99.f;
   if (m < -90 or q < -90) {
     LogDebug("createPFCand") << "<Run3ScoutingParticleToRecoPFCandidateProducer::createPFCand>:" << std::endl
                              << "Unrecognisable pdgId - skipping particle" << std::endl;
     return reco::PFCandidate();
   }
 
   float px = scoutingparticle.pt() * cos(scoutingparticle.phi());
   float py = scoutingparticle.pt() * sin(scoutingparticle.phi());
   float pz = scoutingparticle.pt() * sinh(scoutingparticle.eta());
   float p = scoutingparticle.pt() * cosh(scoutingparticle.eta());
   float energy = std::sqrt(p * p + m * m);
   reco::Particle::LorentzVector p4(px, py, pz, energy);
 
   static const reco::PFCandidate dummy;
   auto pfcand = reco::PFCandidate(q, p4, dummy.translatePdgIdToType(scoutingparticle.pdgId()));
 
   bool relativeTrackVars = scoutingparticle.relative_trk_vars();
   vertexIndex_.push_back(scoutingparticle.vertex());
   normchi2_.push_back(scoutingparticle.normchi2());
   dz_.push_back(scoutingparticle.dz());
   dxy_.push_back(scoutingparticle.dxy());
   dzsig_.push_back(scoutingparticle.dzsig());
   dxysig_.push_back(scoutingparticle.dxysig());
   lostInnerHits_.push_back(scoutingparticle.lostInnerHits());
   quality_.push_back(scoutingparticle.quality());
   trkPt_.push_back(relativeTrackVars ? scoutingparticle.trk_pt() + scoutingparticle.pt() : scoutingparticle.trk_pt());
   trkEta_.push_back(relativeTrackVars ? scoutingparticle.trk_eta() + scoutingparticle.eta()
                                       : scoutingparticle.trk_eta());
   trkPhi_.push_back(relativeTrackVars ? scoutingparticle.trk_phi() + scoutingparticle.phi()
                                       : scoutingparticle.trk_phi());
 
   return pfcand;
 }

◆ createPFCandidates()

void Run3ScoutingParticleToRecoPFCandidateProducer::createPFCandidates	(	edm::Handle< std::vector< Run3ScoutingParticle >>	scoutingparticleHandle,
		std::unique_ptr< reco::PFCandidateCollection > &	pfcands
	)

Definition at line 123 of file Run3ScoutingParticleToRecoPFCandidateProducer.cc.

References createPFCand(), pfDeepBoostedJetPreprocessParams_cfi::pfcand, HLT_2024v13_cff::pfcands, run3scouting_cff::scoutingparticle, and use_CHS_.

Referenced by produce().

                                                        {
   for (unsigned int icand = 0; icand < scoutingparticleHandle->size(); ++icand) {
     auto &scoutingparticle = (*scoutingparticleHandle)[icand];
 
     if (use_CHS_ and scoutingparticle.vertex() > 0)
       continue;
 
     auto pfcand = createPFCand(scoutingparticle);
     if (pfcand.energy() != 0)
       pfcands->push_back(pfcand);
   }
 }

◆ createPFCandidatesSK()

void Run3ScoutingParticleToRecoPFCandidateProducer::createPFCandidatesSK	(	edm::Handle< std::vector< Run3ScoutingParticle >>	scoutingparticleHandle,
		std::unique_ptr< reco::PFCandidateCollection > &	pfcands
	)

Definition at line 138 of file Run3ScoutingParticleToRecoPFCandidateProducer.cc.

References createPFCand(), LogDebug, visualization-live-secondInstance_cfg::m, LHEGenericFilter_cfi::ParticleID, pdTable_, pfDeepBoostedJetPreprocessParams_cfi::pfcand, HLT_2024v13_cff::pfcands, and run3scouting_cff::scoutingparticle.

Referenced by produce().

                                                        {
   std::vector<fastjet::PseudoJet> fj;
 
   for (auto iter = scoutingparticleHandle->begin(),
             ibegin = scoutingparticleHandle->begin(),
             iend = scoutingparticleHandle->end();
        iter != iend;
        ++iter) {
     auto m = pdTable_->particle(HepPDT::ParticleID(iter->pdgId())) != nullptr
                  ? pdTable_->particle(HepPDT::ParticleID(iter->pdgId()))->mass()
                  : -99.f;
     if (m < -90) {
       LogDebug("createPFCandidatesSK") << "<Run3ScoutingParticleToRecoPFCandidateProducer::createPFCandidatesSK>:"
                                        << std::endl
                                        << "Unrecognisable pdgId - skipping particle" << std::endl;
       continue;
     }
     math::PtEtaPhiMLorentzVector p4(iter->pt(), iter->eta(), iter->phi(), m);
     fj.push_back(fastjet::PseudoJet(p4.px(), p4.py(), p4.pz(), p4.energy()));
     fj.back().set_user_index(iter - ibegin);
   }
 
   fastjet::contrib::SoftKiller soft_killer(5, 0.4);
   std::vector<fastjet::PseudoJet> soft_killed_particles = soft_killer(fj);
 
   for (auto &particle : soft_killed_particles) {
     const Run3ScoutingParticle scoutingparticle = scoutingparticleHandle->at(particle.user_index());
     auto pfcand = createPFCand(scoutingparticle);
     if (pfcand.energy() != 0)
       pfcands->push_back(pfcand);
   }
 }

◆ endStream()

void Run3ScoutingParticleToRecoPFCandidateProducer::endStream ( )

inlineoverride

Definition at line 27 of file Run3ScoutingParticleToRecoPFCandidateProducer.cc.

27 {}

◆ fillDescriptions()

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

static

Definition at line 277 of file Run3ScoutingParticleToRecoPFCandidateProducer.cc.

References edm::ConfigurationDescriptions::addWithDefaultLabel(), submitPVResolutionJobs::desc, and ProducerED_cfi::InputTag.

                                                                                                                {
   edm::ParameterSetDescription desc;
   desc.add<edm::InputTag>("scoutingparticle", edm::InputTag("hltScoutingPFPacker"));
   desc.add<bool>("softKiller", false);
   desc.add<bool>("CHS", false);
   descriptions.addWithDefaultLabel(desc);
 }

◆ produce()

void Run3ScoutingParticleToRecoPFCandidateProducer::produce	(	edm::Event &	iEvent,
		edm::EventSetup const &	setup
	)

override

Definition at line 174 of file Run3ScoutingParticleToRecoPFCandidateProducer.cc.

References clearVars(), createPFCandidates(), createPFCandidatesSK(), dxy_, dxysig_, dz_, dzsig_, edm::helper::Filler< Map >::fill(), iEvent, input_scoutingparticle_token_, edm::helper::Filler< Map >::insert(), lostInnerHits_, eostools::move(), normchi2_, particletable_token_, pdTable_, HLT_2024v13_cff::pfcands, quality_, singleTopDQM_cfi::setup, trkEta_, trkPhi_, trkPt_, use_softKiller_, and vertexIndex_.

                                                                                                         {
   using namespace edm;
 
   auto pdt = setup.getHandle(particletable_token_);
   pdTable_ = pdt.product();
 
   Handle<std::vector<Run3ScoutingParticle>> scoutingparticleHandle;
   iEvent.getByToken(input_scoutingparticle_token_, scoutingparticleHandle);
 
   auto pfcands = std::make_unique<reco::PFCandidateCollection>();
 
   if (use_softKiller_) {
     createPFCandidatesSK(scoutingparticleHandle, pfcands);
   } else {
     createPFCandidates(scoutingparticleHandle, pfcands);
   }
 
   edm::OrphanHandle<reco::PFCandidateCollection> oh = iEvent.put(std::move(pfcands));
 
   std::unique_ptr<edm::ValueMap<int>> vertexIndex_VM(new edm::ValueMap<int>());
   edm::ValueMap<int>::Filler filler_vertexIndex(*vertexIndex_VM);
   filler_vertexIndex.insert(oh, vertexIndex_.begin(), vertexIndex_.end());
   filler_vertexIndex.fill();
   iEvent.put(std::move(vertexIndex_VM), "vertexIndex");
 
   std::unique_ptr<edm::ValueMap<float>> normchi2_VM(new edm::ValueMap<float>());
   edm::ValueMap<float>::Filler filler_normchi2(*normchi2_VM);
   filler_normchi2.insert(oh, normchi2_.begin(), normchi2_.end());
   filler_normchi2.fill();
   iEvent.put(std::move(normchi2_VM), "normchi2");
 
   std::unique_ptr<edm::ValueMap<float>> dz_VM(new edm::ValueMap<float>());
   edm::ValueMap<float>::Filler filler_dz(*dz_VM);
   filler_dz.insert(oh, dz_.begin(), dz_.end());
   filler_dz.fill();
   iEvent.put(std::move(dz_VM), "dz");
 
   std::unique_ptr<edm::ValueMap<float>> dxy_VM(new edm::ValueMap<float>());
   edm::ValueMap<float>::Filler filler_dxy(*dxy_VM);
   filler_dxy.insert(oh, dxy_.begin(), dxy_.end());
   filler_dxy.fill();
   iEvent.put(std::move(dxy_VM), "dxy");
 
   std::unique_ptr<edm::ValueMap<float>> dzsig_VM(new edm::ValueMap<float>());
   edm::ValueMap<float>::Filler filler_dzsig(*dzsig_VM);
   filler_dzsig.insert(oh, dzsig_.begin(), dzsig_.end());
   filler_dzsig.fill();
   iEvent.put(std::move(dzsig_VM), "dzsig");
 
   std::unique_ptr<edm::ValueMap<float>> dxysig_VM(new edm::ValueMap<float>());
   edm::ValueMap<float>::Filler filler_dxysig(*dxysig_VM);
   filler_dxysig.insert(oh, dxysig_.begin(), dxysig_.end());
   filler_dxysig.fill();
   iEvent.put(std::move(dxysig_VM), "dxysig");
 
   std::unique_ptr<edm::ValueMap<int>> lostInnerHits_VM(new edm::ValueMap<int>());
   edm::ValueMap<int>::Filler filler_lostInnerHits(*lostInnerHits_VM);
   filler_lostInnerHits.insert(oh, lostInnerHits_.begin(), lostInnerHits_.end());
   filler_lostInnerHits.fill();
   iEvent.put(std::move(lostInnerHits_VM), "lostInnerHits");
 
   std::unique_ptr<edm::ValueMap<int>> quality_VM(new edm::ValueMap<int>());
   edm::ValueMap<int>::Filler filler_quality(*quality_VM);
   filler_quality.insert(oh, quality_.begin(), quality_.end());
   filler_quality.fill();
   iEvent.put(std::move(quality_VM), "quality");
 
   std::unique_ptr<edm::ValueMap<float>> trkPt_VM(new edm::ValueMap<float>());
   edm::ValueMap<float>::Filler filler_trkPt(*trkPt_VM);
   filler_trkPt.insert(oh, trkPt_.begin(), trkPt_.end());
   filler_trkPt.fill();
   iEvent.put(std::move(trkPt_VM), "trkPt");
 
   std::unique_ptr<edm::ValueMap<float>> trkEta_VM(new edm::ValueMap<float>());
   edm::ValueMap<float>::Filler filler_trkEta(*trkEta_VM);
   filler_trkEta.insert(oh, trkEta_.begin(), trkEta_.end());
   filler_trkEta.fill();
   iEvent.put(std::move(trkEta_VM), "trkEta");
 
   std::unique_ptr<edm::ValueMap<float>> trkPhi_VM(new edm::ValueMap<float>());
   edm::ValueMap<float>::Filler filler_trkPhi(*trkPhi_VM);
   filler_trkPhi.insert(oh, trkPhi_.begin(), trkPhi_.end());
   filler_trkPhi.fill();
   iEvent.put(std::move(trkPhi_VM), "trkPhi");
 
   clearVars();
 }