#include <HLTmumutkVtxProducer.h>

Inheritance diagram for HLTmumutkVtxProducer:

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

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

	~HLTmumutkVtxProducer () override=default

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
bool	checkPreviousCand (const reco::TrackRef &trackref, const std::vector< reco::RecoChargedCandidateRef > &ref2) const

bool	overlap (const reco::TrackRef &trackref1, const reco::TrackRef &trackref2)

Static Private Member Functions
static FreeTrajectoryState	initialFreeState (const reco::Track &, const MagneticField *)

Private Attributes
const edm::InputTag	beamSpotTag_

const edm::EDGetTokenT< reco::BeamSpot >	beamSpotToken_

const edm::ESGetToken< MagneticField, IdealMagneticFieldRecord >	idealMagneticFieldRecordToken_

const double	maxEta_

const double	maxInvMass_

const std::string	mfName_

const double	minD0Significance_

const double	minInvMass_

const double	minPt_

const edm::InputTag	muCandTag_

const edm::EDGetTokenT< reco::RecoChargedCandidateCollection >	muCandToken_

const double	overlapDR2_

const edm::InputTag	previousCandTag_

const edm::EDGetTokenT< trigger::TriggerFilterObjectWithRefs >	previousCandToken_

const double	thirdTrackMass_

const edm::ESGetToken< TransientTrackBuilder, TransientTrackRecord >	transientTrackRecordToken_

const edm::InputTag	trkCandTag_

const edm::EDGetTokenT< reco::RecoChargedCandidateCollection >	trkCandToken_

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

HLT Filter for b to (mumu) + X

Implementation: <Notes on="" implementation>="">

Definition at line 46 of file HLTmumutkVtxProducer.h.

Constructor & Destructor Documentation

◆ HLTmumutkVtxProducer()

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

explicit

Definition at line 25 of file HLTmumutkVtxProducer.cc.

     : transientTrackRecordToken_(esConsumes(edm::ESInputTag("", "TransientTrackBuilder"))),
       muCandTag_(iConfig.getParameter<edm::InputTag>("MuCand")),
       muCandToken_(consumes<reco::RecoChargedCandidateCollection>(muCandTag_)),
       trkCandTag_(iConfig.getParameter<edm::InputTag>("TrackCand")),
       trkCandToken_(consumes<reco::RecoChargedCandidateCollection>(trkCandTag_)),
       previousCandTag_(iConfig.getParameter<edm::InputTag>("PreviousCandTag")),
       previousCandToken_(consumes<trigger::TriggerFilterObjectWithRefs>(previousCandTag_)),
       mfName_(iConfig.getParameter<std::string>("SimpleMagneticField")),
       idealMagneticFieldRecordToken_(esConsumes(edm::ESInputTag("", mfName_))),
       thirdTrackMass_(iConfig.getParameter<double>("ThirdTrackMass")),
       maxEta_(iConfig.getParameter<double>("MaxEta")),
       minPt_(iConfig.getParameter<double>("MinPt")),
       minInvMass_(iConfig.getParameter<double>("MinInvMass")),
       maxInvMass_(iConfig.getParameter<double>("MaxInvMass")),
       minD0Significance_(iConfig.getParameter<double>("MinD0Significance")),
       // minimum delta-R^2 threshold (with sign) for non-overlapping tracks
       overlapDR2_(iConfig.getParameter<double>("OverlapDR") * std::abs(iConfig.getParameter<double>("OverlapDR"))),
       beamSpotTag_(iConfig.getParameter<edm::InputTag>("BeamSpotTag")),
       beamSpotToken_(consumes<reco::BeamSpot>(beamSpotTag_)) {
   produces<VertexCollection>();
 }

◆ ~HLTmumutkVtxProducer()

HLTmumutkVtxProducer::~HLTmumutkVtxProducer ( )

overridedefault

Member Function Documentation

◆ checkPreviousCand()

bool HLTmumutkVtxProducer::checkPreviousCand	(	const reco::TrackRef &	trackref,
		const std::vector< reco::RecoChargedCandidateRef > &	ref2
	)		const

private

Definition at line 226 of file HLTmumutkVtxProducer.cc.

References mps_fire::i, and convertSQLiteXML::ok.

Referenced by produce().

                                                                                                    {
   bool ok = false;
   for (auto& i : refVect) {
     if (i->get<TrackRef>() == trackref) {
       ok = true;
       break;
     }
   }
   return ok;
 }

◆ fillDescriptions()

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

static

Definition at line 48 of file HLTmumutkVtxProducer.cc.

References edm::ConfigurationDescriptions::add(), submitPVResolutionJobs::desc, ProducerED_cfi::InputTag, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                       {
   edm::ParameterSetDescription desc;
   desc.add<edm::InputTag>("MuCand", edm::InputTag("hltMuTracks"));
   desc.add<edm::InputTag>("TrackCand", edm::InputTag("hltMumukAllConeTracks"));
   desc.add<edm::InputTag>("PreviousCandTag", edm::InputTag("hltDisplacedmumuFilterDoubleMu4Jpsi"));
   desc.add<std::string>("SimpleMagneticField", "");
   desc.add<double>("ThirdTrackMass", 0.493677);
   desc.add<double>("MaxEta", 2.5);
   desc.add<double>("MinPt", 3.0);
   desc.add<double>("MinInvMass", 0.0);
   desc.add<double>("MaxInvMass", 99999.);
   desc.add<double>("MinD0Significance", 0.0);
   desc.add<double>("OverlapDR", 1.44e-4);
   desc.add<edm::InputTag>("BeamSpotTag", edm::InputTag("hltOfflineBeamSpot"));
   descriptions.add("HLTmumutkVtxProducer", desc);
 }

◆ initialFreeState()

FreeTrajectoryState HLTmumutkVtxProducer::initialFreeState	(	const reco::Track &	tk,
		const MagneticField *	field
	)

staticprivate

Definition at line 212 of file HLTmumutkVtxProducer.cc.

References reco::TrackBase::charge(), reco::TrackBase::covariance(), submitPVResolutionJobs::err, reco::TrackBase::momentum(), and reco::TrackBase::vertex().

Referenced by produce().

                                                                                                           {
   Basic3DVector<float> pos(tk.vertex());
   GlobalPoint gpos(pos);
   Basic3DVector<float> mom(tk.momentum());
   GlobalVector gmom(mom);
   GlobalTrajectoryParameters par(gpos, gmom, tk.charge(), field);
   CurvilinearTrajectoryError err(tk.covariance());
   return FreeTrajectoryState(par, err);
 }

◆ overlap()

bool HLTmumutkVtxProducer::overlap	(	const reco::TrackRef &	trackref1,
		const reco::TrackRef &	trackref2
	)

private

Definition at line 222 of file HLTmumutkVtxProducer.cc.

References reco::deltaR2(), and overlapDR2_.

Referenced by produce().

                                                                                        {
   return (reco::deltaR2(trackref1->eta(), trackref1->phi(), trackref2->eta(), trackref2->phi()) < overlapDR2_);
 }

◆ produce()

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

override

Definition at line 65 of file HLTmumutkVtxProducer.cc.

References funct::abs(), beamSpotToken_, checkPreviousCand(), StorageManager_cfg::e1, SiPixelPhase1Clusters_cfi::e3, edm::Ref< C, T, F >::get(), edm::EventSetup::getHandle(), trigger::TriggerRefsCollections::getObjects(), idealMagneticFieldRecordToken_, iEvent, initialFreeState(), TransientVertex::isValid(), LogDebug, maxEta_, maxInvMass_, minD0Significance_, minInvMass_, minPt_, eostools::move(), muCandToken_, overlap(), AlCaHLTBitMon_ParallelJobs::p, LaserDQM_cfg::p1, SiStripOfflineCRack_cfg::p2, chargedHadronTrackResolutionFilter_cfi::p3, previousCandToken_, Measurement1D::significance(), mathSSE::sqrt(), thirdTrackMass_, transientTrackRecordToken_, TrajectoryStateClosestToBeamLine::transverseImpactParameter(), trigger::TriggerMuon, trkCandToken_, bphysicsOniaDQM_cfi::vertex, KalmanVertexFitter::vertex(), and spclusmultinvestigator_cfi::vertexCollection.

                                                                                 {
   const double MuMass(0.106);
   const double MuMass2(MuMass * MuMass);
   const double thirdTrackMass2(thirdTrackMass_ * thirdTrackMass_);
 
   // get hold of muon trks
   Handle<RecoChargedCandidateCollection> mucands;
   iEvent.getByToken(muCandToken_, mucands);
 
   // get the transient track builder
   auto const& theB = iSetup.getHandle(transientTrackRecordToken_);
 
   //get the beamspot position
   edm::Handle<reco::BeamSpot> recoBeamSpotHandle;
   iEvent.getByToken(beamSpotToken_, recoBeamSpotHandle);
 
   //get the b field
   auto const& bFieldHandle = iSetup.getHandle(idealMagneticFieldRecordToken_);
   const MagneticField* magField = bFieldHandle.product();
   TSCBLBuilderNoMaterial blsBuilder;
 
   // get track candidates around displaced muons
   Handle<RecoChargedCandidateCollection> trkcands;
   iEvent.getByToken(trkCandToken_, trkcands);
 
   unique_ptr<VertexCollection> vertexCollection(new VertexCollection());
 
   // Ref to Candidate object to be recorded in filter object
   RecoChargedCandidateRef refMu1;
   RecoChargedCandidateRef refMu2;
   RecoChargedCandidateRef refTrk;
 
   double e1, e2, e3;
   Particle::LorentzVector p, p1, p2, p3;
 
   if (mucands->size() < 2)
     return;
   if (trkcands->empty())
     return;
 
   RecoChargedCandidateCollection::const_iterator mucand1;
   RecoChargedCandidateCollection::const_iterator mucand2;
   RecoChargedCandidateCollection::const_iterator trkcand;
 
   // get the objects passing the previous filter
   Handle<TriggerFilterObjectWithRefs> previousCands;
   iEvent.getByToken(previousCandToken_, previousCands);
 
   vector<RecoChargedCandidateRef> vPrevCands;
   previousCands->getObjects(TriggerMuon, vPrevCands);
 
   for (mucand1 = mucands->begin(); mucand1 != mucands->end(); ++mucand1) {
     TrackRef trk1 = mucand1->get<TrackRef>();
     LogDebug("HLTmumutkVtxProducer") << " 1st muon: q*pt= " << trk1->charge() * trk1->pt() << ", eta= " << trk1->eta()
                                      << ", hits= " << trk1->numberOfValidHits();
 
     //first check if this muon passed the previous filter
     if (!checkPreviousCand(trk1, vPrevCands))
       continue;
 
     // eta and pt cut
     if (fabs(trk1->eta()) > maxEta_)
       continue;
     if (trk1->pt() < minPt_)
       continue;
 
     mucand2 = mucand1;
     ++mucand2;
     for (; mucand2 != mucands->end(); mucand2++) {
       TrackRef trk2 = mucand2->get<TrackRef>();
 
       LogDebug("HLTDisplacedMumukFilter") << " 2nd muon: q*pt= " << trk2->charge() * trk2->pt()
                                           << ", eta= " << trk2->eta() << ", hits= " << trk2->numberOfValidHits();
 
       //first check if this muon passed the previous filter
       if (!checkPreviousCand(trk2, vPrevCands))
         continue;
       // eta and pt cut
       if (fabs(trk2->eta()) > maxEta_)
         continue;
       if (trk2->pt() < minPt_)
         continue;
 
       //loop on track collection
       for (trkcand = trkcands->begin(); trkcand != trkcands->end(); ++trkcand) {
         TrackRef trk3 = trkcand->get<TrackRef>();
         if (overlap(trk1, trk3))
           continue;
         if (overlap(trk2, trk3))
           continue;
 
         LogDebug("HLTDisplacedMumukFilter") << " 3rd track: q*pt= " << trk3->charge() * trk3->pt()
                                             << ", eta= " << trk3->eta() << ", hits= " << trk3->numberOfValidHits();
 
         // eta and pt cut
         if (fabs(trk3->eta()) > maxEta_)
           continue;
         if (trk3->pt() < minPt_)
           continue;
 
         // Combined system
         e1 = sqrt(trk1->momentum().Mag2() + MuMass2);
         e2 = sqrt(trk2->momentum().Mag2() + MuMass2);
         e3 = sqrt(trk3->momentum().Mag2() + thirdTrackMass2);
 
         p1 = Particle::LorentzVector(trk1->px(), trk1->py(), trk1->pz(), e1);
         p2 = Particle::LorentzVector(trk2->px(), trk2->py(), trk2->pz(), e2);
         p3 = Particle::LorentzVector(trk3->px(), trk3->py(), trk3->pz(), e3);
 
         p = p1 + p2 + p3;
 
         //invariant mass cut
         double invmass = abs(p.mass());
         LogDebug("HLTDisplacedMumukFilter") << " Invmass= " << invmass;
         if (invmass < minInvMass_)
           continue;
         if (invmass > maxInvMass_)
           continue;
 
         // do the vertex fit
         vector<TransientTrack> t_tks;
         t_tks.push_back((*theB).build(&trk1));
         t_tks.push_back((*theB).build(&trk2));
         t_tks.push_back((*theB).build(&trk3));
         if (t_tks.size() != 3)
           continue;
 
         FreeTrajectoryState InitialFTS = initialFreeState(*trk3, magField);
         TrajectoryStateClosestToBeamLine tscb(blsBuilder(InitialFTS, *recoBeamSpotHandle));
         double d0sig = tscb.transverseImpactParameter().significance();
         if (d0sig < minD0Significance_)
           continue;
 
         KalmanVertexFitter kvf;
         TransientVertex tv = kvf.vertex(t_tks);
         if (!tv.isValid())
           continue;
         Vertex vertex = tv;
 
         // put vertex in the event
         vertexCollection->push_back(vertex);
       }
     }
   }
   iEvent.put(std::move(vertexCollection));
 }

Member Data Documentation

◆ beamSpotTag_

const edm::InputTag HLTmumutkVtxProducer::beamSpotTag_

private

Definition at line 80 of file HLTmumutkVtxProducer.h.

◆ beamSpotToken_

const edm::EDGetTokenT<reco::BeamSpot> HLTmumutkVtxProducer::beamSpotToken_

private

Definition at line 81 of file HLTmumutkVtxProducer.h.

Referenced by produce().

◆ idealMagneticFieldRecordToken_

const edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> HLTmumutkVtxProducer::idealMagneticFieldRecordToken_

private

Definition at line 70 of file HLTmumutkVtxProducer.h.

Referenced by produce().

◆ maxEta_

const double HLTmumutkVtxProducer::maxEta_

private

Definition at line 73 of file HLTmumutkVtxProducer.h.

Referenced by produce().

◆ maxInvMass_

const double HLTmumutkVtxProducer::maxInvMass_

private

Definition at line 76 of file HLTmumutkVtxProducer.h.

Referenced by produce().

◆ mfName_

const std::string HLTmumutkVtxProducer::mfName_

private

Definition at line 69 of file HLTmumutkVtxProducer.h.

◆ minD0Significance_

const double HLTmumutkVtxProducer::minD0Significance_

private

Definition at line 77 of file HLTmumutkVtxProducer.h.

Referenced by produce().

◆ minInvMass_

const double HLTmumutkVtxProducer::minInvMass_

private

Definition at line 75 of file HLTmumutkVtxProducer.h.

Referenced by produce().

◆ minPt_

const double HLTmumutkVtxProducer::minPt_

private

Definition at line 74 of file HLTmumutkVtxProducer.h.

Referenced by produce().

◆ muCandTag_

const edm::InputTag HLTmumutkVtxProducer::muCandTag_

private

Definition at line 62 of file HLTmumutkVtxProducer.h.

◆ muCandToken_

const edm::EDGetTokenT<reco::RecoChargedCandidateCollection> HLTmumutkVtxProducer::muCandToken_

private

Definition at line 63 of file HLTmumutkVtxProducer.h.

Referenced by produce().

◆ overlapDR2_

const double HLTmumutkVtxProducer::overlapDR2_

private

Definition at line 78 of file HLTmumutkVtxProducer.h.

Referenced by overlap().

◆ previousCandTag_

const edm::InputTag HLTmumutkVtxProducer::previousCandTag_

private

Definition at line 66 of file HLTmumutkVtxProducer.h.

◆ previousCandToken_

const edm::EDGetTokenT<trigger::TriggerFilterObjectWithRefs> HLTmumutkVtxProducer::previousCandToken_

private

Definition at line 67 of file HLTmumutkVtxProducer.h.

Referenced by produce().

◆ thirdTrackMass_

const double HLTmumutkVtxProducer::thirdTrackMass_

private

Definition at line 72 of file HLTmumutkVtxProducer.h.

Referenced by produce().

◆ transientTrackRecordToken_

const edm::ESGetToken<TransientTrackBuilder, TransientTrackRecord> HLTmumutkVtxProducer::transientTrackRecordToken_

private

Definition at line 60 of file HLTmumutkVtxProducer.h.

Referenced by produce().

◆ trkCandTag_

const edm::InputTag HLTmumutkVtxProducer::trkCandTag_

private

Definition at line 64 of file HLTmumutkVtxProducer.h.

◆ trkCandToken_

const edm::EDGetTokenT<reco::RecoChargedCandidateCollection> HLTmumutkVtxProducer::trkCandToken_

private

Definition at line 65 of file HLTmumutkVtxProducer.h.

Referenced by produce().

Public Member Functions

Static Public Member Functions

Private Member Functions

Static Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ HLTmumutkVtxProducer()

◆ ~HLTmumutkVtxProducer()

Member Function Documentation

◆ checkPreviousCand()

◆ fillDescriptions()

◆ initialFreeState()

◆ overlap()

◆ produce()

Member Data Documentation

◆ beamSpotTag_

◆ beamSpotToken_

◆ idealMagneticFieldRecordToken_

◆ maxEta_

◆ maxInvMass_

◆ mfName_

◆ minD0Significance_

◆ minInvMass_

◆ minPt_

◆ muCandTag_

◆ muCandToken_

◆ overlapDR2_

◆ previousCandTag_

◆ previousCandToken_

◆ thirdTrackMass_

◆ transientTrackRecordToken_

◆ trkCandTag_

◆ trkCandToken_