#include <TauAnalysis/EmbeddingProducer/plugins/MuMuForEmbeddingSelector.cc>

Inheritance diagram for MuMuForEmbeddingSelector:

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

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
edm::EDGetTokenT< reco::BeamSpot >	theBeamSpotLabel_

edm::EDGetTokenT< edm::View< pat::MET > >	theMETLabel_

edm::EDGetTokenT< edm::View< pat::MET > >	thePuppiMETLabel_

edm::EDGetTokenT< reco::VertexCollection >	theVertexLabel_

bool	use_zmass = false

edm::EDGetTokenT< edm::View< reco::CompositeCandidate > >	ZmumuCandidates_

Static Private Attributes
static constexpr double	zmass = 91.1876

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

Description: [one line class summary]

Implementation: [Notes on implementation]

Definition at line 46 of file MuMuForEmbeddingSelector.cc.

Constructor & Destructor Documentation

◆ MuMuForEmbeddingSelector()

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

explicit

Definition at line 76 of file MuMuForEmbeddingSelector.cc.

References edm::ParameterSet::getParameter(), theBeamSpotLabel_, theMETLabel_, thePuppiMETLabel_, theVertexLabel_, and use_zmass.

     : ZmumuCandidates_(consumes<edm::View<reco::CompositeCandidate>>(
           iConfig.getParameter<edm::InputTag>("ZmumuCandidatesCollection"))) {
   use_zmass = iConfig.getParameter<bool>("use_zmass");
   produces<edm::RefVector<pat::MuonCollection>>();
   produces<float>("oldMass");
   produces<float>("newMass");
   produces<float>("nPairCandidates");
   produces<bool>("isMediumLeadingMuon");
   produces<bool>("isTightLeadingMuon");
   produces<bool>("isMediumTrailingMuon");
   produces<bool>("isTightTrailingMuon");
   produces<float>("initialMETEt");
   produces<float>("initialMETphi");
   produces<float>("initialPuppiMETEt");
   produces<float>("initialPuppiMETphi");
   theVertexLabel_ = consumes<reco::VertexCollection>(iConfig.getParameter<edm::InputTag>("inputTagVertex"));
   theBeamSpotLabel_ = consumes<reco::BeamSpot>(iConfig.getParameter<edm::InputTag>("inputTagBeamSpot"));
   theMETLabel_ = consumes<edm::View<pat::MET>>(iConfig.getParameter<edm::InputTag>("Met"));
   thePuppiMETLabel_ = consumes<edm::View<pat::MET>>(iConfig.getParameter<edm::InputTag>("PuppiMet"));
   // now do what ever other initialization is needed
 }

Member Function Documentation

◆ fillDescriptions()

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

static

Definition at line 183 of file MuMuForEmbeddingSelector.cc.

References edm::ConfigurationDescriptions::addDefault(), and submitPVResolutionJobs::desc.

                                                                                           {
   // The following says we do not know what parameters are allowed so do no validation
   //  Please change this to state exactly what you do use, even if it is no parameters
   edm::ParameterSetDescription desc;
   desc.setUnknown();
   descriptions.addDefault(desc);
 }

◆ produce()

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

overrideprivate

Definition at line 104 of file MuMuForEmbeddingSelector.cc.

References funct::abs(), pwdgSkimBPark_cfi::beamSpot, edm::RefToBase< T >::castTo(), reco::CompositeCandidate::daughter(), iEvent, muon::isMediumMuon(), muon::isTightMuon(), reco::LeafCandidate::mass(), reco::Candidate::masterClone(), BTaggingMonitor_cfi::met, eostools::move(), BeamMonitor_cff::primaryVertex, dumpMFGeometry_cfg::prod, theBeamSpotLabel_, theMETLabel_, thePuppiMETLabel_, theVertexLabel_, use_zmass, bphysicsOniaDQM_cfi::vertex, L1BJetProducer_cff::vtx, zmass, SelectingProcedure_cff::ZmumuCandidates, and ZmumuCandidates_.

                                                                                     {
   using namespace edm;
   edm::Handle<edm::View<reco::CompositeCandidate>> ZmumuCandidatesHandle;
   iEvent.getByToken(ZmumuCandidates_, ZmumuCandidatesHandle);
   edm::View<reco::CompositeCandidate> ZmumuCandidates = *ZmumuCandidatesHandle;
   const reco::CompositeCandidate *chosenZCand = nullptr;
   const reco::CompositeCandidate *chosenZCand_zmass = nullptr;
   const reco::CompositeCandidate *chosenZCand_largest = nullptr;
   double massDifference = 9999;
   edm::Handle<reco::BeamSpot> beamSpot;
   iEvent.getByToken(theBeamSpotLabel_, beamSpot);
   edm::Handle<reco::VertexCollection> vertex;
   iEvent.getByToken(theVertexLabel_, vertex);
   edm::Handle<edm::View<pat::MET>> met;
   iEvent.getByToken(theMETLabel_, met);
   edm::Handle<edm::View<pat::MET>> puppimet;
   iEvent.getByToken(thePuppiMETLabel_, puppimet);
   // get primary vertex
   reco::Vertex::Point posVtx;
   reco::Vertex::Error errVtx;
   for (const auto &vtx : *vertex) {
     if (vtx.isValid() && !vtx.isFake()) {
       posVtx = vtx.position();
       errVtx = vtx.error();
       break;
     }
   }
   reco::Vertex primaryVertex(posVtx, errVtx);
 
   for (edm::View<reco::CompositeCandidate>::const_iterator iZCand = ZmumuCandidates.begin();
        iZCand != ZmumuCandidates.end();
        ++iZCand) {
     if (std::abs(zmass - iZCand->mass()) < massDifference) {
       massDifference = std::abs(zmass - iZCand->mass());
       chosenZCand_zmass = &(*iZCand);
     }
   }
   for (edm::View<reco::CompositeCandidate>::const_iterator iZCand = ZmumuCandidates.begin();
        iZCand != ZmumuCandidates.end();
        ++iZCand) {
     if (chosenZCand_largest == nullptr) {
       chosenZCand_largest = &(*iZCand);
     } else {
       if (iZCand->mass() > chosenZCand_largest->mass()) {
         chosenZCand_largest = &(*iZCand);
       }
     }
   }
   if (use_zmass) {
     chosenZCand = chosenZCand_zmass;
   } else {
     chosenZCand = chosenZCand_largest;
   }
 
   std::unique_ptr<edm::RefVector<pat::MuonCollection>> prod(new edm::RefVector<pat::MuonCollection>());
   prod->reserve(2);
   prod->push_back(chosenZCand->daughter(0)->masterClone().castTo<pat::MuonRef>());
   prod->push_back(chosenZCand->daughter(1)->masterClone().castTo<pat::MuonRef>());
   iEvent.put(std::move(prod));
   iEvent.put(std::make_unique<float>(chosenZCand_zmass->mass()), "oldMass");
   iEvent.put(std::make_unique<float>(chosenZCand_largest->mass()), "newMass");
   iEvent.put(std::make_unique<float>(ZmumuCandidates.size()), "nPairCandidates");
   iEvent.put(std::make_unique<bool>(chosenZCand->daughter(0)->masterClone().castTo<pat::MuonRef>()->isMediumMuon()),
              "isMediumLeadingMuon");
   iEvent.put(std::make_unique<bool>(
                  chosenZCand->daughter(0)->masterClone().castTo<pat::MuonRef>()->isTightMuon(primaryVertex)),
              "isTightLeadingMuon");
   iEvent.put(std::make_unique<bool>(chosenZCand->daughter(1)->masterClone().castTo<pat::MuonRef>()->isMediumMuon()),
              "isMediumTrailingMuon");
   iEvent.put(std::make_unique<bool>(
                  chosenZCand->daughter(1)->masterClone().castTo<pat::MuonRef>()->isTightMuon(primaryVertex)),
              "isTightTrailingMuon");
   iEvent.put(std::make_unique<float>(met->at(0).et()), "initialMETEt");
   iEvent.put(std::make_unique<float>(met->at(0).phi()), "initialMETphi");
   iEvent.put(std::make_unique<float>(puppimet->at(0).et()), "initialPuppiMETEt");
   iEvent.put(std::make_unique<float>(puppimet->at(0).phi()), "initialPuppiMETphi");
 }