Inheritance diagram for ShiftedPFCandidateProducerByMatchedObject:

Classes
struct	objectEntryType

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

	~ShiftedPFCandidateProducerByMatchedObject () 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
double	dR2match_Object_

double	dR2match_PFCandidate_

double	dRmatch_Object_

double	dRmatch_PFCandidate_

std::vector< objectEntryType >	objects_

edm::EDGetTokenT < reco::PFCandidateCollection >	srcPFCandidates_

edm::EDGetTokenT< edm::View < reco::Candidate > >	srcShiftedObjects_

edm::EDGetTokenT< edm::View < reco::Candidate > >	srcUnshiftedObjects_

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

Vary energy of PFCandidates coinciding in eta-phi with selected electrons/muons/tau-jets/jets by electron/muon/tau-jet/jet energy uncertainty.

Author: Christian Veelken, LLR

Definition at line 28 of file ShiftedPFCandidateProducerByMatchedObject.cc.

Constructor & Destructor Documentation

ShiftedPFCandidateProducerByMatchedObject::ShiftedPFCandidateProducerByMatchedObject ( const edm::ParameterSet & cfg )

explicit

Definition at line 71 of file ShiftedPFCandidateProducerByMatchedObject.cc.

References dR2match_Object_, dR2match_PFCandidate_, dRmatch_Object_, dRmatch_PFCandidate_, edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), srcPFCandidates_, srcShiftedObjects_, and srcUnshiftedObjects_.

                                                                                                                {
   srcPFCandidates_ = consumes<reco::PFCandidateCollection>(cfg.getParameter<edm::InputTag>("srcPFCandidates"));
   srcUnshiftedObjects_ = consumes<edm::View<reco::Candidate> >(cfg.getParameter<edm::InputTag>("srcUnshiftedObjects"));
   srcShiftedObjects_ = consumes<edm::View<reco::Candidate> >(cfg.getParameter<edm::InputTag>("srcShiftedObjects"));
 
   dRmatch_PFCandidate_ = cfg.getParameter<double>("dRmatch_PFCandidate");
   dR2match_PFCandidate_ = dRmatch_PFCandidate_ * dRmatch_PFCandidate_;
   dRmatch_Object_ = cfg.exists("dRmatch_Object") ? cfg.getParameter<double>("dRmatch_Object") : 0.1;
   dR2match_Object_ = dRmatch_Object_ * dRmatch_Object_;
   produces<reco::PFCandidateCollection>();
 }

ShiftedPFCandidateProducerByMatchedObject::~ShiftedPFCandidateProducerByMatchedObject ( )

override

Definition at line 83 of file ShiftedPFCandidateProducerByMatchedObject.cc.

                                                                                       {
   // nothing to be done yet...
 }

Member Function Documentation

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

overrideprivate

Definition at line 87 of file ShiftedPFCandidateProducerByMatchedObject.cc.

References reco::deltaR2(), dR2match_Object_, dR2match_PFCandidate_, dRDefault, edm::Event::getByToken(), ResonanceBuilder::mass, eostools::move(), objects_, edm::Event::put(), reco::LeafCandidate::setP4(), edm::shift, mathSSE::sqrt(), srcPFCandidates_, srcShiftedObjects_, and srcUnshiftedObjects_.

                                                                                               {
   edm::Handle<reco::PFCandidateCollection> originalPFCandidates;
   evt.getByToken(srcPFCandidates_, originalPFCandidates);
 
   typedef edm::View<reco::Candidate> CandidateView;
 
   edm::Handle<CandidateView> unshiftedObjects;
   evt.getByToken(srcUnshiftedObjects_, unshiftedObjects);
 
   edm::Handle<CandidateView> shiftedObjects;
   evt.getByToken(srcShiftedObjects_, shiftedObjects);
 
   objects_.clear();
 
   CandidateView::const_iterator shiftedObjectP4_matched;
   bool isMatched_Object = false;
   double dR2bestMatch_Object = dRDefault;
   for (CandidateView::const_iterator unshiftedObject = unshiftedObjects->begin();
        unshiftedObject != unshiftedObjects->end();
        ++unshiftedObject) {
     isMatched_Object = false;
     dR2bestMatch_Object = dRDefault;
 
     for (CandidateView::const_iterator shiftedObject = shiftedObjects->begin(); shiftedObject != shiftedObjects->end();
          ++shiftedObject) {
       double dR2 = deltaR2(unshiftedObject->p4(), shiftedObject->p4());
       if (dR2 < dR2match_Object_ && dR2 < dR2bestMatch_Object) {
         shiftedObjectP4_matched = shiftedObject;
         isMatched_Object = true;
         dR2bestMatch_Object = dR2;
       }
     }
     if (isMatched_Object) {
       objects_.push_back(
           objectEntryType(shiftedObjectP4_matched->p4(), unshiftedObject->p4(), sqrt(dR2bestMatch_Object)));
     }
   }
 
   auto shiftedPFCandidates = std::make_unique<reco::PFCandidateCollection>();
 
   for (reco::PFCandidateCollection::const_iterator originalPFCandidate = originalPFCandidates->begin();
        originalPFCandidate != originalPFCandidates->end();
        ++originalPFCandidate) {
     double shift = 0.;
     bool applyShift = false;
     double dR2bestMatch_PFCandidate = dRDefault;
     for (std::vector<objectEntryType>::const_iterator object = objects_.begin(); object != objects_.end(); ++object) {
       if (!object->isValidMatch_)
         continue;
       double dR2 = deltaR2(originalPFCandidate->p4(), object->unshiftedObjectP4_);
       if (dR2 < dR2match_PFCandidate_ && dR2 < dR2bestMatch_PFCandidate) {
         shift = object->shift_;
         applyShift = true;
         dR2bestMatch_PFCandidate = dR2;
       }
     }
 
     reco::Candidate::LorentzVector shiftedPFCandidateP4 = originalPFCandidate->p4();
     if (applyShift) {
       double shiftedPx = (1. + shift) * originalPFCandidate->px();
       double shiftedPy = (1. + shift) * originalPFCandidate->py();
       double shiftedPz = (1. + shift) * originalPFCandidate->pz();
       double mass = originalPFCandidate->mass();
       double shiftedEn = sqrt(shiftedPx * shiftedPx + shiftedPy * shiftedPy + shiftedPz * shiftedPz + mass * mass);
       shiftedPFCandidateP4.SetPxPyPzE(shiftedPx, shiftedPy, shiftedPz, shiftedEn);
     }
 
     reco::PFCandidate shiftedPFCandidate(*originalPFCandidate);
     shiftedPFCandidate.setP4(shiftedPFCandidateP4);
 
     shiftedPFCandidates->push_back(shiftedPFCandidate);
   }
 
   evt.put(std::move(shiftedPFCandidates));
 }