#include <JetMuonHitsIDHelper.h>

Public Member Functions
void	calculate (const edm::Event &event, const edm::EventSetup &isetup, const reco::Jet &jet, const int iDbg=0)

void	fillDescription (edm::ParameterSetDescription &iDesc)

void	initValues ()

	JetMuonHitsIDHelper ()

	JetMuonHitsIDHelper (edm::ParameterSet const &pset, edm::ConsumesCollector &&iC)

int	numberOfHits1RPC () const

int	numberOfHits2RPC () const

int	numberOfHits3RPC () const

int	numberOfHits4RPC () const

int	numberOfHitsRPC () const

	~JetMuonHitsIDHelper ()

Private Attributes
edm::EDGetTokenT< RPCRecHitCollection >	input_rpchits_token_

bool	isRECO_

int	numberOfHits1RPC_

int	numberOfHits2RPC_

int	numberOfHits3RPC_

int	numberOfHits4RPC_

int	numberOfHitsRPC_

edm::InputTag	rpcRecHits_

Detailed Description

Definition at line 18 of file JetMuonHitsIDHelper.h.

Constructor & Destructor Documentation

reco::helper::JetMuonHitsIDHelper::JetMuonHitsIDHelper ( )

inline

Definition at line 21 of file JetMuonHitsIDHelper.h.

References muonDTDigis_cfi::pset.

21 {}

reco::helper::JetMuonHitsIDHelper::JetMuonHitsIDHelper	(	edm::ParameterSet const &	pset,
		edm::ConsumesCollector &&	iC
	)

Definition at line 32 of file JetMuonHitsIDHelper.cc.

References edm::ParameterSet::getParameter().

                                                                                                          {
   isRECO_ = true;  // This will be "true" initially, then if the product isn't found, set to false once
   numberOfHits1RPC_ = 0;
   numberOfHits2RPC_ = 0;
   numberOfHits3RPC_ = 0;
   numberOfHits4RPC_ = 0;
   numberOfHitsRPC_ = 0;
   rpcRecHits_ = pset.getParameter<edm::InputTag>("rpcRecHits");
 
   input_rpchits_token_ = iC.consumes<RPCRecHitCollection>(rpcRecHits_);
 }

reco::helper::JetMuonHitsIDHelper::~JetMuonHitsIDHelper ( )

inline

Definition at line 23 of file JetMuonHitsIDHelper.h.

References calculate(), fillDescription(), initValues(), and metsig::jet.

23 {}

Member Function Documentation

void reco::helper::JetMuonHitsIDHelper::calculate	(	const edm::Event &	event,
		const edm::EventSetup &	isetup,
		const reco::Jet &	jet,
		const int	iDbg = `0`
	)

Definition at line 44 of file JetMuonHitsIDHelper.cc.

References reco::deltaR(), reco::LeafCandidate::eta(), TrackingRecHit::geographicalId(), edm::EventSetup::get(), runTauDisplay::gp, GlobalTrackingGeometry::idToDet(), edm::HandleBase::isValid(), RPCRecHit::localPosition(), reco::LeafCandidate::phi(), RPCDetId::station(), and GeomDet::toGlobal().

Referenced by JetIDProducer::produce(), and ~JetMuonHitsIDHelper().

                                                                   {
   // initialize
   numberOfHits1RPC_ = 0;
   numberOfHits2RPC_ = 0;
   numberOfHits3RPC_ = 0;
   numberOfHits4RPC_ = 0;
   numberOfHitsRPC_ = 0;
 
   if (isRECO_) {  // This will be "true" initially, then if the product isn't found, set to false once
 
     // Get tracking geometry
     edm::ESHandle<GlobalTrackingGeometry> trackingGeometry;
     iSetup.get<GlobalTrackingGeometryRecord>().get(trackingGeometry);
 
     //####READ RPC RecHits Collection########
     //#In config: RpcRecHits     = cms.InputTag("rpcRecHits")
     edm::Handle<RPCRecHitCollection> rpcRecHits_handle;
     event.getByToken(input_rpchits_token_, rpcRecHits_handle);
 
     if (!rpcRecHits_handle.isValid()) {
       // don't throw exception if not running on RECO
       edm::LogWarning("DataNotAvailable") << "JetMuonHitsIDHelper will not be run at all, this is not a RECO file.";
       isRECO_ = false;
       return;
     }
 
     //####calculate rpc  variables for each jet########
 
     for (RPCRecHitCollection::const_iterator itRPC = rpcRecHits_handle->begin(), itRPCEnd = rpcRecHits_handle->end();
          itRPC != itRPCEnd;
          ++itRPC) {
       RPCRecHit const& hit = *itRPC;
       DetId detid = hit.geographicalId();
       LocalPoint lp = hit.localPosition();
       const GeomDet* gd = trackingGeometry->idToDet(detid);
       GlobalPoint gp = gd->toGlobal(lp);
       double dR2 = reco::deltaR(jet.eta(), jet.phi(), static_cast<double>(gp.eta()), static_cast<double>(gp.phi()));
       if (dR2 < 0.5) {
         RPCDetId rpcChamberId = (RPCDetId)detid;
         numberOfHitsRPC_++;
         if (rpcChamberId.station() == 1)
           numberOfHits1RPC_++;
         if (rpcChamberId.station() == 2)
           numberOfHits2RPC_++;
         if (rpcChamberId.station() == 3)
           numberOfHits3RPC_++;
         if (rpcChamberId.station() == 4)
           numberOfHits4RPC_++;
       }
     }
   }
 }