#include <HLTmumutkVtxProducer.h>

Inheritance diagram for HLTmumutkVtxProducer:

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

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

	~HLTmumutkVtxProducer ()

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducer ()

Public Member Functions inherited from edm::ProducerBase
	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

virtual	~ProducerBase ()

Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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

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

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

Private Attributes
edm::InputTag	beamSpotTag_

edm::EDGetTokenT< reco::BeamSpot >	beamSpotToken_

const double	maxEta_

const double	maxInvMass_

const std::string	mfName_

const double	minD0Significance_

const double	minInvMass_

const double	minPt_

edm::InputTag	muCandTag_

edm::EDGetTokenT < reco::RecoChargedCandidateCollection >	muCandToken_

const double	overlapDR_

edm::InputTag	previousCandTag_

edm::EDGetTokenT < trigger::TriggerFilterObjectWithRefs >	previousCandToken_

const double	thirdTrackMass_

edm::InputTag	trkCandTag_

edm::EDGetTokenT < reco::RecoChargedCandidateCollection >	trkCandToken_

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Protected Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

HLT Filter for b to (mumu) + X

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

Definition at line 42 of file HLTmumutkVtxProducer.h.

Constructor & Destructor Documentation

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 ( )

Definition at line 53 of file HLTmumutkVtxProducer.cc.

53 {}

Member Function Documentation

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

private

Definition at line 225 of file HLTmumutkVtxProducer.cc.

References i, and convertSQLiteXML::ok.

Referenced by produce().

                                                                                                                {
   bool ok=false;
   for (unsigned int i=0; i<refVect.size(); i++) {
     if ( refVect[i]->get<TrackRef>() == trackref ) {
       ok=true;
       break;
     }
   }
   return ok;
 }

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

static

Definition at line 56 of file HLTmumutkVtxProducer.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), 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);
 }

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

staticprivate

Definition at line 209 of file HLTmumutkVtxProducer.cc.

References reco::TrackBase::charge(), reco::TrackBase::covariance(), 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);
 }

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

private

Definition at line 220 of file HLTmumutkVtxProducer.cc.

References deltaR(), and overlapDR_.

Referenced by produce().

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

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

virtual

Implements edm::EDProducer.

Definition at line 74 of file HLTmumutkVtxProducer.cc.

References funct::abs(), beamSpotToken_, checkPreviousCand(), edm::EventSetup::get(), edm::Ref< C, T, F >::get(), edm::Event::getByToken(), initialFreeState(), TransientVertex::isValid(), LogDebug, maxEta_, maxInvMass_, mfName_, minD0Significance_, minInvMass_, minPt_, muCandToken_, overlap(), AlCaHLTBitMon_ParallelJobs::p, p1, p2, p3, previousCandToken_, edm::ESHandle< class >::product(), edm::Event::put(), Measurement1D::significance(), mathSSE::sqrt(), thirdTrackMass_, TrajectoryStateClosestToBeamLine::transverseImpactParameter(), trigger::TriggerMuon, trkCandToken_, KalmanVertexFitter::vertex(), and GoodVertex_cfg::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:
   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);
 
   auto_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->size() < 1 )   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(vertexCollection);
 }

Member Data Documentation

edm::InputTag HLTmumutkVtxProducer::beamSpotTag_

private

Definition at line 72 of file HLTmumutkVtxProducer.h.

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

private

Definition at line 73 of file HLTmumutkVtxProducer.h.

Referenced by produce().

const double HLTmumutkVtxProducer::maxEta_

private

Definition at line 65 of file HLTmumutkVtxProducer.h.

Referenced by produce().

const double HLTmumutkVtxProducer::maxInvMass_

private

Definition at line 68 of file HLTmumutkVtxProducer.h.

Referenced by produce().

const std::string HLTmumutkVtxProducer::mfName_

private

Definition at line 62 of file HLTmumutkVtxProducer.h.

Referenced by produce().

const double HLTmumutkVtxProducer::minD0Significance_

private

Definition at line 69 of file HLTmumutkVtxProducer.h.

Referenced by produce().

const double HLTmumutkVtxProducer::minInvMass_

private

Definition at line 67 of file HLTmumutkVtxProducer.h.

Referenced by produce().

const double HLTmumutkVtxProducer::minPt_

private

Definition at line 66 of file HLTmumutkVtxProducer.h.

Referenced by produce().

edm::InputTag HLTmumutkVtxProducer::muCandTag_

private

Definition at line 55 of file HLTmumutkVtxProducer.h.

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

private

Definition at line 56 of file HLTmumutkVtxProducer.h.

Referenced by produce().

const double HLTmumutkVtxProducer::overlapDR_

private

Definition at line 70 of file HLTmumutkVtxProducer.h.

Referenced by overlap().

edm::InputTag HLTmumutkVtxProducer::previousCandTag_

private

Definition at line 59 of file HLTmumutkVtxProducer.h.

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

private

Definition at line 60 of file HLTmumutkVtxProducer.h.

Referenced by produce().

const double HLTmumutkVtxProducer::thirdTrackMass_

private

Definition at line 64 of file HLTmumutkVtxProducer.h.

Referenced by produce().

edm::InputTag HLTmumutkVtxProducer::trkCandTag_

private

Definition at line 57 of file HLTmumutkVtxProducer.h.

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

private

Definition at line 58 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

Member Function Documentation

Member Data Documentation