HLTmumutkVtxProducer Class Reference

#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 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::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 43 of file HLTmumutkVtxProducer.h.

Constructor & Destructor Documentation

HLTmumutkVtxProducer()

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

explicit

Definition at line 31 of file HLTmumutkVtxProducer.cc.

    : 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")),
      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 237 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 54 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 221 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 231 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 72 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:
  edm::ESHandle<TransientTrackBuilder> theB;
  iSetup.get<TransientTrackRecord>().get("TransientTrackBuilder", theB);
 
  //get the beamspot position
  edm::Handle<reco::BeamSpot> recoBeamSpotHandle;
  iEvent.getByToken(beamSpotToken_, recoBeamSpotHandle);
 
  //get the b field
  ESHandle<MagneticField> bFieldHandle;
  iSetup.get<IdealMagneticFieldRecord>().get(mfName_, bFieldHandle);
  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::EventSetup::get(), get, edm::Ref< C, T, F >::get(), trigger::TriggerRefsCollections::getObjects(), iEvent, initialFreeState(), TransientVertex::isValid(), LogDebug, maxEta_, maxInvMass_, mfName_, minD0Significance_, minInvMass_, minPt_, eostools::move(), muCandToken_, overlap(), AlCaHLTBitMon_ParallelJobs::p, p1, p2, p3, previousCandToken_, edm::ESHandle< T >::product(), Measurement1D::significance(), mathSSE::sqrt(), thirdTrackMass_, 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 72 of file HLTmumutkVtxProducer.h.

beamSpotToken_

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

private

Definition at line 73 of file HLTmumutkVtxProducer.h.

Referenced by produce().

maxEta_

const double HLTmumutkVtxProducer::maxEta_

private

Definition at line 65 of file HLTmumutkVtxProducer.h.

Referenced by produce().

maxInvMass_

const double HLTmumutkVtxProducer::maxInvMass_

private

Definition at line 68 of file HLTmumutkVtxProducer.h.

Referenced by produce().

mfName_

const std::string HLTmumutkVtxProducer::mfName_

private

Definition at line 62 of file HLTmumutkVtxProducer.h.

Referenced by produce().

minD0Significance_

const double HLTmumutkVtxProducer::minD0Significance_

private

Definition at line 69 of file HLTmumutkVtxProducer.h.

Referenced by produce().

minInvMass_

const double HLTmumutkVtxProducer::minInvMass_

private

Definition at line 67 of file HLTmumutkVtxProducer.h.

Referenced by produce().

minPt_

const double HLTmumutkVtxProducer::minPt_

private

Definition at line 66 of file HLTmumutkVtxProducer.h.

Referenced by produce().

muCandTag_

const edm::InputTag HLTmumutkVtxProducer::muCandTag_

private

Definition at line 55 of file HLTmumutkVtxProducer.h.

muCandToken_

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

private

Definition at line 56 of file HLTmumutkVtxProducer.h.

Referenced by produce().

overlapDR_

const double HLTmumutkVtxProducer::overlapDR_

private

Definition at line 70 of file HLTmumutkVtxProducer.h.

Referenced by overlap().

previousCandTag_

const edm::InputTag HLTmumutkVtxProducer::previousCandTag_

private

Definition at line 59 of file HLTmumutkVtxProducer.h.

previousCandToken_

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

private

Definition at line 60 of file HLTmumutkVtxProducer.h.

Referenced by produce().

thirdTrackMass_

const double HLTmumutkVtxProducer::thirdTrackMass_

private

Definition at line 64 of file HLTmumutkVtxProducer.h.

Referenced by produce().

trkCandTag_

const edm::InputTag HLTmumutkVtxProducer::trkCandTag_

private

Definition at line 57 of file HLTmumutkVtxProducer.h.

trkCandToken_

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

private

Definition at line 58 of file HLTmumutkVtxProducer.h.

Referenced by produce().