Inheritance diagram for PFTICLProducer:

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

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

	~PFTICLProducer () 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
const edm::EDGetTokenT < reco::MuonCollection >	muons_

std::unique_ptr< PFMuonAlgo >	pfmu_

const edm::EDGetTokenT < edm::ValueMap< float > >	srcTrackTime_

const edm::EDGetTokenT < edm::ValueMap< float > >	srcTrackTimeError_

const edm::EDGetTokenT < edm::ValueMap< float > >	srcTrackTimeQuality_

const edm::EDGetTokenT < edm::View< TICLCandidate > >	ticl_candidates_

const float	timingQualityThreshold_

const bool	useMTDTiming_

const bool	useTimingAverage_

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 PFTICLProducer.cc.

Constructor & Destructor Documentation

PFTICLProducer::PFTICLProducer ( const edm::ParameterSet & conf )

Definition at line 43 of file PFTICLProducer.cc.

     : useMTDTiming_(conf.getParameter<bool>("useMTDTiming")),
       useTimingAverage_(conf.getParameter<bool>("useTimingAverage")),
       timingQualityThreshold_(conf.getParameter<double>("timingQualityThreshold")),
       ticl_candidates_(consumes<edm::View<TICLCandidate>>(conf.getParameter<edm::InputTag>("ticlCandidateSrc"))),
       srcTrackTime_(consumes<edm::ValueMap<float>>(conf.getParameter<edm::InputTag>("trackTimeValueMap"))),
       srcTrackTimeError_(consumes<edm::ValueMap<float>>(conf.getParameter<edm::InputTag>("trackTimeErrorMap"))),
       srcTrackTimeQuality_(consumes<edm::ValueMap<float>>(conf.getParameter<edm::InputTag>("trackTimeQualityMap"))),
       muons_(consumes<reco::MuonCollection>(conf.getParameter<edm::InputTag>("muonSrc"))),
       pfmu_(std::make_unique<PFMuonAlgo>(conf.getParameterSet("pfMuonAlgoParameters"),
                                          false)) {  // postMuonCleaning = false
   produces<reco::PFCandidateCollection>();
 }

PFTICLProducer::~PFTICLProducer ( )

inlineoverride

Definition at line 21 of file PFTICLProducer.cc.

21 {}

Member Function Documentation

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

static

Definition at line 57 of file PFTICLProducer.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), submitPVResolutionJobs::desc, PFMuonAlgo::fillPSetDescription(), and HLT_FULL_cff::InputTag.

                                                                                 {
   edm::ParameterSetDescription desc;
   desc.add<edm::InputTag>("ticlCandidateSrc", edm::InputTag("ticlTrackstersMerge"));
   desc.add<edm::InputTag>("trackTimeValueMap", edm::InputTag("tofPID:t0"));
   desc.add<edm::InputTag>("trackTimeErrorMap", edm::InputTag("tofPID:sigmat0"));
   desc.add<edm::InputTag>("trackTimeQualityMap", edm::InputTag("mtdTrackQualityMVA:mtdQualMVA"));
   desc.add<double>("timingQualityThreshold", 0.5);
   desc.add<bool>("useMTDTiming", true);
   desc.add<bool>("useTimingAverage", false);
   // For PFMuonAlgo
   desc.add<edm::InputTag>("muonSrc", edm::InputTag("muons1stStep"));
   edm::ParameterSetDescription psd_PFMuonAlgo;
   PFMuonAlgo::fillPSetDescription(psd_PFMuonAlgo);
   desc.add<edm::ParameterSetDescription>("pfMuonAlgoParameters", psd_PFMuonAlgo);
   //
   descriptions.add("pfTICLProducer", desc);
 }

void PFTICLProducer::produce	(	edm::Event &	evt,
		const edm::EventSetup &	es
	)

override

Definition at line 75 of file PFTICLProducer.cc.

References funct::abs(), HLT_FULL_cff::candidates, RecoTauCleanerPlugins::charge, reco::PFCandidate::e, reco::PFCandidate::gamma, edm::Event::getByToken(), edm::Event::getHandle(), reco::PFCandidate::h, reco::PFCandidate::h0, eostools::move(), reco::PFCandidate::mu, PFMuonAlgo::muAssocToTrack(), patZpeak::muons, muons_, pfmu_, funct::pow(), edm::Event::productGetter(), edm::Event::put(), mathSSE::sqrt(), srcTrackTime_, srcTrackTimeError_, srcTrackTimeQuality_, submitPVValidationJobs::t, ticl_candidates_, timingQualityThreshold_, useMTDTiming_, useTimingAverage_, and reco::PFCandidate::X.

                                                                    {
   //get TICLCandidates
   edm::Handle<edm::View<TICLCandidate>> ticl_cand_h;
   evt.getByToken(ticl_candidates_, ticl_cand_h);
   const auto ticl_candidates = *ticl_cand_h;
   edm::Handle<edm::ValueMap<float>> trackTimeH, trackTimeErrH, trackTimeQualH;
   evt.getByToken(srcTrackTime_, trackTimeH);
   evt.getByToken(srcTrackTimeError_, trackTimeErrH);
   evt.getByToken(srcTrackTimeQuality_, trackTimeQualH);
   const auto muonH = evt.getHandle(muons_);
   const auto muons = *muonH;
 
   auto candidates = std::make_unique<reco::PFCandidateCollection>();
 
   for (const auto& ticl_cand : ticl_candidates) {
     const auto abs_pdg_id = std::abs(ticl_cand.pdgId());
     const auto charge = ticl_cand.charge();
     const auto& four_mom = ticl_cand.p4();
     double ecal_energy = 0.;
 
     for (const auto& t : ticl_cand.tracksters()) {
       double ecal_energy_fraction = t->raw_em_pt() / t->raw_pt();
       ecal_energy += t->raw_energy() * ecal_energy_fraction;
     }
     double hcal_energy = ticl_cand.rawEnergy() - ecal_energy;
     // fix for floating point rounding could go slightly below 0
     hcal_energy = hcal_energy < 0 ? 0 : hcal_energy;
 
     reco::PFCandidate::ParticleType part_type;
     switch (abs_pdg_id) {
       case 11:
         part_type = reco::PFCandidate::e;
         break;
       case 13:
         part_type = reco::PFCandidate::mu;
         break;
       case 22:
         part_type = reco::PFCandidate::gamma;
         break;
       case 130:
         part_type = reco::PFCandidate::h0;
         break;
       case 211:
         part_type = reco::PFCandidate::h;
         break;
       // default also handles neutral pions (111) for the time being (not yet foreseen in PFCandidate)
       default:
         part_type = reco::PFCandidate::X;
     }
 
     candidates->emplace_back(charge, four_mom, part_type);
 
     auto& candidate = candidates->back();
     candidate.setEcalEnergy(ecal_energy, ecal_energy);
     candidate.setHcalEnergy(hcal_energy, hcal_energy);
     if (candidate.charge()) {  // otherwise PFCandidate throws
       // Construct edm::Ref from edm::Ptr. As of now, assumes type to be reco::Track. To be extended (either via
       // dynamic type checking or configuration) if additional track types are needed.
       reco::TrackRef trackref(ticl_cand.trackPtr().id(), int(ticl_cand.trackPtr().key()), &evt.productGetter());
       candidate.setTrackRef(trackref);
       // Utilize PFMuonAlgo
       const int muId = PFMuonAlgo::muAssocToTrack(trackref, muons);
       if (muId != -1) {
         const reco::MuonRef muonref = reco::MuonRef(muonH, muId);
         const bool allowLoose = (part_type == reco::PFCandidate::mu);
         // Redefine pfmuon candidate kinematics and add muonref
         pfmu_->reconstructMuon(candidate, muonref, allowLoose);
       }
     }
 
     // HGCAL timing as default values
     auto time = ticl_cand.time();
     auto timeE = ticl_cand.timeError();
 
     if (useMTDTiming_ and candidate.charge()) {
       // Ignore HGCAL timing until it will be TOF corrected
       time = -99.;
       timeE = -1.;
       // Check MTD timing availability
       const bool assocQuality = (*trackTimeQualH)[candidate.trackRef()] > timingQualityThreshold_;
       if (assocQuality) {
         const auto timeHGC = time;
         const auto timeEHGC = timeE;
         const auto timeMTD = (*trackTimeH)[candidate.trackRef()];
         const auto timeEMTD = (*trackTimeErrH)[candidate.trackRef()];
 
         if (useTimingAverage_ && (timeEMTD > 0 && timeEHGC > 0)) {
           // Compute weighted average between HGCAL and MTD timing
           const auto invTimeESqHGC = pow(timeEHGC, -2);
           const auto invTimeESqMTD = pow(timeEMTD, -2);
           timeE = (invTimeESqHGC * invTimeESqMTD) / (invTimeESqHGC + invTimeESqMTD);
           time = (timeHGC * invTimeESqHGC + timeMTD * invTimeESqMTD) * timeE;
           timeE = sqrt(timeE);
         } else if (timeEMTD > 0) {  // Ignore HGCal timing until it will be TOF corrected
           time = timeMTD;
           timeE = timeEMTD;
         }
       }
     }
     candidate.setTime(time, timeE);
   }
 
   evt.put(std::move(candidates));
 }

Member Data Documentation

const edm::EDGetTokenT<reco::MuonCollection> PFTICLProducer::muons_

private

Definition at line 36 of file PFTICLProducer.cc.

Referenced by produce().

std::unique_ptr<PFMuonAlgo> PFTICLProducer::pfmu_

private

Definition at line 38 of file PFTICLProducer.cc.

Referenced by produce().

const edm::EDGetTokenT<edm::ValueMap<float> > PFTICLProducer::srcTrackTime_

private

Definition at line 35 of file PFTICLProducer.cc.

Referenced by produce().

const edm::EDGetTokenT<edm::ValueMap<float> > PFTICLProducer::srcTrackTimeError_

private

Definition at line 35 of file PFTICLProducer.cc.

Referenced by produce().

const edm::EDGetTokenT<edm::ValueMap<float> > PFTICLProducer::srcTrackTimeQuality_

private

Definition at line 35 of file PFTICLProducer.cc.

Referenced by produce().

const edm::EDGetTokenT<edm::View<TICLCandidate> > PFTICLProducer::ticl_candidates_

private

Definition at line 34 of file PFTICLProducer.cc.

Referenced by produce().

const float PFTICLProducer::timingQualityThreshold_

private

Definition at line 31 of file PFTICLProducer.cc.

Referenced by produce().

const bool PFTICLProducer::useMTDTiming_

private

Definition at line 29 of file PFTICLProducer.cc.

Referenced by produce().

const bool PFTICLProducer::useTimingAverage_

private

Definition at line 30 of file PFTICLProducer.cc.

Referenced by produce().

Public Member Functions

Static Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation