#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 19 of file JetMuonHitsIDHelper.h.

Constructor & Destructor Documentation

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

inline

Definition at line 23 of file JetMuonHitsIDHelper.h.

23 {}

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 25 of file JetMuonHitsIDHelper.h.

25 {}

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 49 of file JetMuonHitsIDHelper.cc.

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

Referenced by JetIDProducer::produce().

 {
 
   // 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_++;
       }
     }
   }
 }