#include <HLTDisplacedmumuVtxProducer.h>

Inheritance diagram for HLTDisplacedmumuVtxProducer:

Public Member Functions
virtual void	beginJob ()

virtual void	endJob ()

	HLTDisplacedmumuVtxProducer (const edm::ParameterSet &)

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

	~HLTDisplacedmumuVtxProducer ()

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

virtual	~EDProducer ()

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

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

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

virtual	~ProducerBase ()

Private Attributes
int	chargeOpt_

double	maxEta_

double	maxInvMass_

double	minInvMass_

double	minPt_

double	minPtPair_

edm::InputTag	src_

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

typedef WorkerT< EDProducer >	WorkerType

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

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDProducer
CurrentProcessingContext const *	currentContext () const

Protected Member Functions inherited from edm::ProducerBase
template<class TProducer , class TMethod >
void	callWhenNewProductsRegistered (TProducer *iProd, TMethod iMethod)

Detailed Description

Definition at line 24 of file HLTDisplacedmumuVtxProducer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 33 of file HLTDisplacedmumuVtxProducer.cc.

                                                                                       : 
     src_ (iConfig.getParameter<edm::InputTag>("Src")),
     maxEta_ (iConfig.getParameter<double>("MaxEta")),
     minPt_ (iConfig.getParameter<double>("MinPt")),
     minPtPair_ (iConfig.getParameter<double>("MinPtPair")),
     minInvMass_ (iConfig.getParameter<double>("MinInvMass")),
     maxInvMass_ (iConfig.getParameter<double>("MaxInvMass")),
     chargeOpt_ (iConfig.getParameter<int>("ChargeOpt"))
 {
     produces<VertexCollection>();
 }

HLTDisplacedmumuVtxProducer::~HLTDisplacedmumuVtxProducer ( )

Definition at line 46 of file HLTDisplacedmumuVtxProducer.cc.

 {
 
 }

Member Function Documentation

void HLTDisplacedmumuVtxProducer::beginJob ( void )

virtual

Reimplemented from edm::EDProducer.

Definition at line 53 of file HLTDisplacedmumuVtxProducer.cc.

 {
 
 }

void HLTDisplacedmumuVtxProducer::endJob ( void )

virtual

Reimplemented from edm::EDProducer.

Definition at line 59 of file HLTDisplacedmumuVtxProducer.cc.

 {
     
 }

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

virtual

Implements edm::EDProducer.

Definition at line 65 of file HLTDisplacedmumuVtxProducer.cc.

References abs, chargeOpt_, edm::EventSetup::get(), edm::Ref< C, T, F >::get(), edm::Event::getByLabel(), TransientVertex::isValid(), LogDebug, maxEta_, maxInvMass_, minInvMass_, minPt_, minPtPair_, L1TEmulatorMonitor_cff::p, p1, p2, edm::Event::put(), mathSSE::sqrt(), src_, KalmanVertexFitter::vertex(), and GoodVertex_cfg::vertexCollection.

 {
     double const MuMass = 0.106;
     double const MuMass2 = MuMass*MuMass;
     
     
     // get hold of muon trks
     Handle<RecoChargedCandidateCollection> mucands;
     iEvent.getByLabel (src_,mucands);
     
     //get the transient track builder:
     edm::ESHandle<TransientTrackBuilder> theB;
     iSetup.get<TransientTrackRecord>().get("TransientTrackBuilder",theB);
 
     std::auto_ptr<VertexCollection> vertexCollection(new VertexCollection());
 
     // look at all mucands,  check cuts and make vertices
     double e1,e2;
     Particle::LorentzVector p,p1,p2;
     
     RecoChargedCandidateCollection::const_iterator cand1;
     RecoChargedCandidateCollection::const_iterator cand2;
     for (cand1=mucands->begin(); cand1!=mucands->end(); cand1++) {
            TrackRef tk1 = cand1->get<TrackRef>();
           
            LogDebug("HLTDisplacedMumuFilter") << " 1st muon in loop: q*pt= " << tk1->charge()*tk1->pt() << ", eta= " << tk1->eta() << ", hits= " << tk1->numberOfValidHits();
             
            // cuts
            if (fabs(tk1->eta())>maxEta_) continue;
            if (tk1->pt() < minPt_) continue;
           
            cand2 = cand1; cand2++;
            for (; cand2!=mucands->end(); cand2++) {
                  TrackRef tk2 = cand2->get<TrackRef>();
          
              // eta cut
              LogDebug("HLTMuonDimuonFilter") << " 2nd muon in loop: q*pt= " << tk2->charge()*tk2->pt() << ", eta= " << tk2->eta() << ", hits= " << tk2->numberOfValidHits() << ", d0= " << tk2->d0();
              
              // cuts
              if (fabs(tk2->eta())>maxEta_) continue;
              if (tk2->pt() < minPt_) continue;
              
              // opposite sign or same sign
              if (chargeOpt_<0) {
                if (tk1->charge()*tk2->charge()>0) continue;
              } else if (chargeOpt_>0) {
                if (tk1->charge()*tk2->charge()<0) continue;
              }
              
              // Combined dimuon system
              e1 = sqrt(tk1->momentum().Mag2()+MuMass2);
              e2 = sqrt(tk2->momentum().Mag2()+MuMass2);
              p1 = Particle::LorentzVector(tk1->px(),tk1->py(),tk1->pz(),e1);
              p2 = Particle::LorentzVector(tk2->px(),tk2->py(),tk2->pz(),e2);
              p = p1+p2;
              
              
              if (p.pt()<minPtPair_) continue;
              
              double invmass = abs(p.mass());
              LogDebug("HLTDisplacedMumuFilter") << " ... 1-2 invmass= " << invmass;
              
              if (invmass<minInvMass_) continue;
              if (invmass>maxInvMass_) continue;
              
              // do the vertex fit
              vector<TransientTrack> t_tks;
              TransientTrack ttkp1 = (*theB).build(&tk1);
              TransientTrack ttkp2 = (*theB).build(&tk2);
              t_tks.push_back(ttkp1);
              t_tks.push_back(ttkp2);
              
         
              if (t_tks.size()!=2) 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);
 }