#include <TrackingRegionsFromBeamSpotAndL2Tau.h>

Inheritance diagram for TrackingRegionsFromBeamSpotAndL2Tau:

Public Member Functions
virtual std::vector < TrackingRegion * >	regions (const edm::Event &e, const edm::EventSetup &es) const

	TrackingRegionsFromBeamSpotAndL2Tau (const edm::ParameterSet &conf, edm::ConsumesCollector &&iC)

virtual	~TrackingRegionsFromBeamSpotAndL2Tau ()

Public Member Functions inherited from TrackingRegionProducer
virtual	~TrackingRegionProducer ()

Private Attributes
float	m_deltaEta

float	m_deltaPhi

float	m_jetMaxEta

int	m_jetMaxN

float	m_jetMinPt

float	m_originHalfLength

float	m_originRadius

bool	m_precise

float	m_ptMin

bool	m_searchOpt

RectangularEtaPhiTrackingRegion::UseMeasurementTracker	m_whereToUseMeasurementTracker

edm::EDGetTokenT< reco::BeamSpot >	token_beamSpot

edm::EDGetTokenT < reco::CandidateView >	token_jet

edm::EDGetTokenT < MeasurementTrackerEvent >	token_measurementTracker

Detailed Description

class TrackingRegionsFromBeamSpotAndL2Tau plugin for creating eta-phi TrackingRegions in directions of L2 taus

Definition at line 25 of file TrackingRegionsFromBeamSpotAndL2Tau.h.

Constructor & Destructor Documentation

TrackingRegionsFromBeamSpotAndL2Tau::TrackingRegionsFromBeamSpotAndL2Tau	(	const edm::ParameterSet &	conf,
		edm::ConsumesCollector &&	iC
	)

inlineexplicit

Definition at line 29 of file TrackingRegionsFromBeamSpotAndL2Tau.h.

References edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), m_deltaEta, m_deltaPhi, m_jetMaxEta, m_jetMaxN, m_jetMinPt, m_originHalfLength, m_originRadius, m_precise, m_ptMin, m_searchOpt, m_whereToUseMeasurementTracker, AlCaHLTBitMon_QueryRunRegistry::string, token_beamSpot, token_jet, and token_measurementTracker.

   {
     edm::LogInfo ("TrackingRegionsFromBeamSpotAndL2Tau") << "Enter the TrackingRegionsFromBeamSpotAndL2Tau";
 
     edm::ParameterSet regionPSet = conf.getParameter<edm::ParameterSet>("RegionPSet");
 
     m_ptMin            = regionPSet.getParameter<double>("ptMin");
     m_originRadius     = regionPSet.getParameter<double>("originRadius");
     m_originHalfLength = regionPSet.getParameter<double>("originHalfLength");
     m_deltaEta         = regionPSet.getParameter<double>("deltaEta");
     m_deltaPhi         = regionPSet.getParameter<double>("deltaPhi");
     token_jet          = iC.consumes<reco::CandidateView>(regionPSet.getParameter<edm::InputTag>("JetSrc"));
     m_jetMinPt         = regionPSet.getParameter<double>("JetMinPt");
     m_jetMaxEta        = regionPSet.getParameter<double>("JetMaxEta");
     m_jetMaxN          = regionPSet.getParameter<int>("JetMaxN");
     token_beamSpot     = iC.consumes<reco::BeamSpot>(regionPSet.getParameter<edm::InputTag>("beamSpot"));
     m_precise          = regionPSet.getParameter<bool>("precise");
 
     if (regionPSet.exists("searchOpt")) m_searchOpt = regionPSet.getParameter<bool>("searchOpt");
     else                                m_searchOpt = false;
 
     m_whereToUseMeasurementTracker = RectangularEtaPhiTrackingRegion::UseMeasurementTracker::kForSiStrips;
     if (regionPSet.exists("measurementTrackerName"))
     {
       // FIXME: when next time altering the configuration of this
       // class, please change the types of the following parameters:
       // - whereToUseMeasurementTracker to at least int32 or to a string
       //   corresponding to the UseMeasurementTracker enumeration
       // - measurementTrackerName to InputTag
       if (regionPSet.exists("whereToUseMeasurementTracker"))
         m_whereToUseMeasurementTracker = RectangularEtaPhiTrackingRegion::doubleToUseMeasurementTracker(regionPSet.getParameter<double>("whereToUseMeasurementTracker"));
       if(m_whereToUseMeasurementTracker != RectangularEtaPhiTrackingRegion::UseMeasurementTracker::kNever)
         token_measurementTracker = iC.consumes<MeasurementTrackerEvent>(regionPSet.getParameter<std::string>("measurementTrackerName"));
     }
   }

virtual TrackingRegionsFromBeamSpotAndL2Tau::~TrackingRegionsFromBeamSpotAndL2Tau ( )

inlinevirtual

Definition at line 66 of file TrackingRegionsFromBeamSpotAndL2Tau.h.

66 {}

Member Function Documentation

virtual std::vector<TrackingRegion* > TrackingRegionsFromBeamSpotAndL2Tau::regions	(	const edm::Event &	e,
		const edm::EventSetup &	es
	)		const

inlinevirtual

Implements TrackingRegionProducer.

Definition at line 69 of file TrackingRegionsFromBeamSpotAndL2Tau.h.

References funct::abs(), reco::Candidate::eta(), edm::Event::getByToken(), i, edm::EDGetTokenT< T >::isUninitialized(), edm::HandleBase::isValid(), metsig::jet, m_deltaEta, m_deltaPhi, m_jetMaxEta, m_jetMaxN, m_jetMinPt, m_originHalfLength, m_originRadius, m_precise, m_ptMin, m_searchOpt, m_whereToUseMeasurementTracker, reco::Candidate::momentum(), edm::Handle< T >::product(), reco::Candidate::pt(), query::result, token_beamSpot, token_jet, token_measurementTracker, reco::BeamSpot::x0(), reco::BeamSpot::y0(), and reco::BeamSpot::z0().

   {
     std::vector<TrackingRegion* > result;
 
     // use beam spot to pick up the origin
     edm::Handle<reco::BeamSpot> bsHandle;
     e.getByToken( token_beamSpot, bsHandle);
     if(!bsHandle.isValid()) return result;
     const reco::BeamSpot & bs = *bsHandle;
     GlobalPoint origin(bs.x0(), bs.y0(), bs.z0());
 
     // pick up the candidate objects of interest
     edm::Handle< reco::CandidateView > objects;
     e.getByToken( token_jet, objects );
     size_t n_objects = objects->size();
     if (n_objects == 0) return result;
 
     const MeasurementTrackerEvent *measurementTracker = nullptr;
     if(!token_measurementTracker.isUninitialized()) {
       edm::Handle<MeasurementTrackerEvent> hmte;
       e.getByToken(token_measurementTracker, hmte);
       measurementTracker = hmte.product();
     }
 
     // create maximum JetMaxN tracking regions in directions of 
     // highest pt jets that are above threshold and are within allowed eta
     // (we expect that jet collection was sorted in decreasing pt order)
     int n_regions = 0;
     for (size_t i =0; i < n_objects && n_regions < m_jetMaxN; ++i)
     {
       const reco::Candidate & jet = (*objects)[i];
       if ( jet.pt() < m_jetMinPt || std::abs(jet.eta()) > m_jetMaxEta ) continue;
       
       GlobalVector direction(jet.momentum().x(), jet.momentum().y(), jet.momentum().z());
 
       RectangularEtaPhiTrackingRegion* etaphiRegion = new RectangularEtaPhiTrackingRegion(
           direction,
           origin,
           m_ptMin,
           m_originRadius,
           m_originHalfLength,
           m_deltaEta,
           m_deltaPhi,
           m_whereToUseMeasurementTracker,
           m_precise,
           measurementTracker,
           m_searchOpt
       );
       result.push_back(etaphiRegion);
       ++n_regions;
     }
     //std::cout<<"nregions = "<<n_regions<<std::endl;
     return result;
   }