#include <VertexBeamspotOrigins.h>

Public Types
enum	OperationMode { OperationMode::BEAM_SPOT_FIXED, OperationMode::BEAM_SPOT_SIGMA, OperationMode::VERTICES_FIXED, OperationMode::VERTICES_SIGMA }

using	Origins = std::vector< std::pair< GlobalPoint, float > >

Public Member Functions
Origins	origins (const edm::Event &iEvent) const

	VertexBeamspotOrigins (const edm::ParameterSet &regPSet, edm::ConsumesCollector &&iC)

	VertexBeamspotOrigins (const edm::ParameterSet &regPSet, edm::ConsumesCollector &iC)

	~VertexBeamspotOrigins ()=default

Static Public Member Functions
static void	fillDescriptions (edm::ParameterSetDescription &desc, const std::string &defaultBeamSpot="offlineBeamSpot", const std::string &defaultVertex="firstStepPrimaryVertices", int defaultMaxVertices=-1)

Private Attributes
int	m_maxNVertices

float	m_nSigmaZBeamSpot

float	m_nSigmaZVertex

OperationMode	m_operationMode

float	m_zErrorBeamSpot

float	m_zErrorVertex

edm::EDGetTokenT< reco::BeamSpot >	token_beamSpot

edm::EDGetTokenT< reco::VertexCollection >	token_vertex

Detailed Description

Definition at line 15 of file VertexBeamspotOrigins.h.

Member Typedef Documentation

using VertexBeamspotOrigins::Origins = std::vector<std::pair<GlobalPoint, float> >

Definition at line 17 of file VertexBeamspotOrigins.h.

Member Enumeration Documentation

enum VertexBeamspotOrigins::OperationMode

strong

Enumerator
BEAM_SPOT_FIXED
BEAM_SPOT_SIGMA
VERTICES_FIXED
VERTICES_SIGMA

Definition at line 18 of file VertexBeamspotOrigins.h.

18 { BEAM_SPOT_FIXED, BEAM_SPOT_SIGMA, VERTICES_FIXED, VERTICES_SIGMA };

Constructor & Destructor Documentation

VertexBeamspotOrigins::VertexBeamspotOrigins	(	const edm::ParameterSet &	regPSet,
		edm::ConsumesCollector &&	iC
	)

inline

Definition at line 20 of file VertexBeamspotOrigins.h.

References fillDescriptions(), iEvent, origins(), AlCaHLTBitMon_QueryRunRegistry::string, and ~VertexBeamspotOrigins().

21 : VertexBeamspotOrigins(regPSet, iC) {}

VertexBeamspotOrigins::VertexBeamspotOrigins

VertexBeamspotOrigins(const edm::ParameterSet &regPSet, edm::ConsumesCollector &&iC)

Definition: VertexBeamspotOrigins.h:20

VertexBeamspotOrigins::VertexBeamspotOrigins	(	const edm::ParameterSet &	regPSet,
		edm::ConsumesCollector &	iC
	)

Definition at line 3 of file VertexBeamspotOrigins.cc.

References BEAM_SPOT_FIXED, BEAM_SPOT_SIGMA, edm::ConsumesCollector::consumes(), Exception, edm::ParameterSet::getParameter(), m_maxNVertices, m_nSigmaZBeamSpot, m_nSigmaZVertex, m_operationMode, m_zErrorBeamSpot, m_zErrorVertex, AlCaHLTBitMon_QueryRunRegistry::string, token_beamSpot, token_vertex, VERTICES_FIXED, and VERTICES_SIGMA.

                                                                                                      {
   // operation mode
   std::string operationModeString = regPSet.getParameter<std::string>("operationMode");
   if (operationModeString == "BeamSpotFixed")
     m_operationMode = OperationMode::BEAM_SPOT_FIXED;
   else if (operationModeString == "BeamSpotSigma")
     m_operationMode = OperationMode::BEAM_SPOT_SIGMA;
   else if (operationModeString == "VerticesFixed")
     m_operationMode = OperationMode::VERTICES_FIXED;
   else if (operationModeString == "VerticesSigma")
     m_operationMode = OperationMode::VERTICES_SIGMA;
   else
     throw cms::Exception("Configuration") << "Unknown operation mode string: " << operationModeString;
 
   token_beamSpot = iC.consumes<reco::BeamSpot>(regPSet.getParameter<edm::InputTag>("beamSpot"));
   m_maxNVertices = 1;
   if (m_operationMode == OperationMode::VERTICES_FIXED || m_operationMode == OperationMode::VERTICES_SIGMA) {
     token_vertex = iC.consumes<reco::VertexCollection>(regPSet.getParameter<edm::InputTag>("vertexCollection"));
     m_maxNVertices = regPSet.getParameter<int>("maxNVertices");
   }
 
   // mode-dependent z-halflength of tracking regions
   m_zErrorBeamSpot = regPSet.getParameter<double>("zErrorBeamSpot");
   if (m_operationMode == OperationMode::VERTICES_SIGMA)
     m_nSigmaZVertex = regPSet.getParameter<double>("nSigmaZVertex");
   if (m_operationMode == OperationMode::VERTICES_FIXED)
     m_zErrorVertex = regPSet.getParameter<double>("zErrorVertex");
   m_nSigmaZBeamSpot = -1.;
   if (m_operationMode == OperationMode::BEAM_SPOT_SIGMA) {
     m_nSigmaZBeamSpot = regPSet.getParameter<double>("nSigmaZBeamSpot");
     if (m_nSigmaZBeamSpot < 0.)
       throw cms::Exception("Configuration") << "nSigmaZBeamSpot must be positive for BeamSpotSigma mode!";
   }
 }

VertexBeamspotOrigins::~VertexBeamspotOrigins ( )

default

Referenced by VertexBeamspotOrigins().

Member Function Documentation

void VertexBeamspotOrigins::fillDescriptions	(	edm::ParameterSetDescription &	desc,
		const std::string &	defaultBeamSpot = `"offlineBeamSpot"`,
		const std::string &	defaultVertex = `"firstStepPrimaryVertices"`,
		int	defaultMaxVertices = `-1`
	)

static

Definition at line 38 of file VertexBeamspotOrigins.cc.

References edm::ParameterSetDescription::add(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by PixelInactiveAreaTrackingRegionsSeedingLayersProducer::fillDescriptions(), AreaSeededTrackingRegionsProducer::fillDescriptions(), CandidatePointSeededTrackingRegionsProducer::fillDescriptions(), and VertexBeamspotOrigins().

                                                                      {
   desc.add<std::string>("operationMode", "BeamSpotFixed");
   desc.add<edm::InputTag>("beamSpot", defaultBeamSpot);
   desc.add<edm::InputTag>("vertexCollection", defaultVertex);
   desc.add<int>("maxNVertices", defaultMaxVertices);
 
   desc.add<double>("nSigmaZBeamSpot", 4.);
   desc.add<double>("zErrorBeamSpot", 24.2);
   desc.add<double>("nSigmaZVertex", 3.);
   desc.add<double>("zErrorVertex", 0.2);
 }

VertexBeamspotOrigins::Origins VertexBeamspotOrigins::origins ( const edm::Event & iEvent ) const

Definition at line 53 of file VertexBeamspotOrigins.cc.

References BEAM_SPOT_FIXED, BEAM_SPOT_SIGMA, edm::Event::getByToken(), edm::HandleBase::isValid(), m_maxNVertices, m_nSigmaZBeamSpot, m_nSigmaZVertex, m_operationMode, m_zErrorBeamSpot, m_zErrorVertex, runTheMatrix::ret, reco::BeamSpot::sigmaZ(), token_beamSpot, token_vertex, findQualityFiles::v, pwdgSkimBPark_cfi::vertices, VERTICES_FIXED, VERTICES_SIGMA, reco::BeamSpot::x0(), reco::BeamSpot::y0(), reco::BeamSpot::z0(), and reco::BeamSpot::z0Error().

Referenced by PixelInactiveAreaTrackingRegionsSeedingLayersProducer::produce(), AreaSeededTrackingRegionsProducer::regions(), CandidatePointSeededTrackingRegionsProducer::regions(), and VertexBeamspotOrigins().

                                                                                         {
   Origins ret;
 
   // always need the beam spot (as a fall back strategy for vertex modes)
   edm::Handle<reco::BeamSpot> bs;
   iEvent.getByToken(token_beamSpot, bs);
   if (!bs.isValid())
     return ret;
 
   // this is a default origin for all modes
   GlobalPoint default_origin(bs->x0(), bs->y0(), bs->z0());
 
   // fill the origins and halfLengths depending on the mode
   if (m_operationMode == OperationMode::BEAM_SPOT_FIXED || m_operationMode == OperationMode::BEAM_SPOT_SIGMA) {
     ret.emplace_back(
         default_origin,
         (m_operationMode == OperationMode::BEAM_SPOT_FIXED) ? m_zErrorBeamSpot : m_nSigmaZBeamSpot * bs->sigmaZ());
   } else if (m_operationMode == OperationMode::VERTICES_FIXED || m_operationMode == OperationMode::VERTICES_SIGMA) {
     edm::Handle<reco::VertexCollection> vertices;
     iEvent.getByToken(token_vertex, vertices);
     int n_vert = 0;
     for (const auto& v : *vertices) {
       if (v.isFake() || !v.isValid())
         continue;
 
       ret.emplace_back(
           GlobalPoint(v.x(), v.y(), v.z()),
           (m_operationMode == OperationMode::VERTICES_FIXED) ? m_zErrorVertex : m_nSigmaZVertex * v.zError());
       ++n_vert;
       if (m_maxNVertices >= 0 && n_vert >= m_maxNVertices) {
         break;
       }
     }
     // no-vertex fall-back case:
     if (ret.empty()) {
       ret.emplace_back(default_origin, (m_nSigmaZBeamSpot > 0.) ? m_nSigmaZBeamSpot * bs->z0Error() : m_zErrorBeamSpot);
     }
   }
 
   return ret;
 }