#include <HLTmumutkVtxProducer.h>

Inheritance diagram for HLTmumutkVtxProducer:

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

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

	~HLTmumutkVtxProducer () override

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

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	overlapDR_

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<>
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Detailed Description

HLT Filter for b to (mumu) + X

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

Definition at line 47 of file HLTmumutkVtxProducer.h.

Constructor & Destructor Documentation

◆ HLTmumutkVtxProducer()

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

explicit

Definition at line 26 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")),
       overlapDR_(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 232 of file HLTmumutkVtxProducer.cc.

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

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

Referenced by produce().

◆ fillDescriptions()

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

static

Definition at line 51 of file HLTmumutkVtxProducer.cc.

                                                                                       {
   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);
 }

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

◆ initialFreeState()

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

staticprivate

Definition at line 216 of file HLTmumutkVtxProducer.cc.

                                                                                                           {
   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);
 }

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

Referenced by produce().

◆ overlap()

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

private

Definition at line 226 of file HLTmumutkVtxProducer.cc.

                                                                                        {
   if (deltaR(trackref1->eta(), trackref1->phi(), trackref2->eta(), trackref2->phi()) < overlapDR_)
     return true;
   return false;
 }

References PbPb_ZMuSkimMuonDPG_cff::deltaR, and overlapDR_.

Referenced by produce().

◆ produce()

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

override

Definition at line 69 of file HLTmumutkVtxProducer.cc.

                                                                                 {
   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));
 }

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, p1, p2, p3, previousCandToken_, Measurement1D::significance(), mathSSE::sqrt(), thirdTrackMass_, transientTrackRecordToken_, TrajectoryStateClosestToBeamLine::transverseImpactParameter(), trigger::TriggerMuon, trkCandToken_, bphysicsOniaDQM_cfi::vertex, KalmanVertexFitter::vertex(), and spclusmultinvestigator_cfi::vertexCollection.

Member Data Documentation

◆ beamSpotTag_

const edm::InputTag HLTmumutkVtxProducer::beamSpotTag_

private

Definition at line 79 of file HLTmumutkVtxProducer.h.

◆ beamSpotToken_

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

private

Definition at line 80 of file HLTmumutkVtxProducer.h.

Referenced by produce().

◆ idealMagneticFieldRecordToken_

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

private

Definition at line 69 of file HLTmumutkVtxProducer.h.

Referenced by produce().

◆ maxEta_

const double HLTmumutkVtxProducer::maxEta_

private

Definition at line 72 of file HLTmumutkVtxProducer.h.

Referenced by produce().

◆ maxInvMass_

const double HLTmumutkVtxProducer::maxInvMass_

private

Definition at line 75 of file HLTmumutkVtxProducer.h.

Referenced by produce().

◆ mfName_

const std::string HLTmumutkVtxProducer::mfName_

private

Definition at line 68 of file HLTmumutkVtxProducer.h.

◆ minD0Significance_

const double HLTmumutkVtxProducer::minD0Significance_

private

Definition at line 76 of file HLTmumutkVtxProducer.h.

Referenced by produce().

◆ minInvMass_

const double HLTmumutkVtxProducer::minInvMass_

private

Definition at line 74 of file HLTmumutkVtxProducer.h.

Referenced by produce().

◆ minPt_

const double HLTmumutkVtxProducer::minPt_

private

Definition at line 73 of file HLTmumutkVtxProducer.h.

Referenced by produce().

◆ muCandTag_

const edm::InputTag HLTmumutkVtxProducer::muCandTag_

private

Definition at line 61 of file HLTmumutkVtxProducer.h.

◆ muCandToken_

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

private

Definition at line 62 of file HLTmumutkVtxProducer.h.

Referenced by produce().

◆ overlapDR_

const double HLTmumutkVtxProducer::overlapDR_

private

Definition at line 77 of file HLTmumutkVtxProducer.h.

Referenced by overlap().

◆ previousCandTag_

const edm::InputTag HLTmumutkVtxProducer::previousCandTag_

private

Definition at line 65 of file HLTmumutkVtxProducer.h.

◆ previousCandToken_

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

private

Definition at line 66 of file HLTmumutkVtxProducer.h.

Referenced by produce().

◆ thirdTrackMass_

const double HLTmumutkVtxProducer::thirdTrackMass_

private

Definition at line 71 of file HLTmumutkVtxProducer.h.

Referenced by produce().

◆ transientTrackRecordToken_

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

private

Definition at line 59 of file HLTmumutkVtxProducer.h.

Referenced by produce().

◆ trkCandTag_

const edm::InputTag HLTmumutkVtxProducer::trkCandTag_

private

Definition at line 63 of file HLTmumutkVtxProducer.h.

◆ trkCandToken_

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

private

Definition at line 64 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_

◆ overlapDR_

◆ previousCandTag_

◆ previousCandToken_

◆ thirdTrackMass_

◆ transientTrackRecordToken_

◆ trkCandTag_

◆ trkCandToken_