PatVertexAnalyzer Class Reference

Inheritance diagram for PatVertexAnalyzer:

Public Member Functions
void	analyze (const edm::Event &event, const edm::EventSetup &es) override
void	beginJob () override
	PatVertexAnalyzer (const edm::ParameterSet &params)
	constructor and destructor More...
	~PatVertexAnalyzer () override
Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzer ()
SerialTaskQueue *	globalLuminosityBlocksQueue ()
SerialTaskQueue *	globalRunsQueue ()
ModuleDescription const &	moduleDescription () const
std::string	workerType () const
	~EDAnalyzer () override
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)
virtual	~EDConsumerBase () noexcept(false)

Private Attributes
edm::EDGetTokenT< reco::GenParticleCollection >	genParticlesToken_
TH1 *	nTracks_
TH1 *	nVertices_
edm::EDGetTokenT< reco::VertexCollection >	srcToken_
TH1 *	x_
TH1 *	xDelta_
TH1 *	xErr_
TH1 *	xPull_
TH1 *	y_
TH1 *	yDelta_
TH1 *	yErr_
TH1 *	yPull_
TH1 *	z_
TH1 *	zDelta_
TH1 *	zErr_
TH1 *	zPull_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)
static bool	wantsGlobalLuminosityBlocks ()
static bool	wantsGlobalRuns ()
static bool	wantsStreamLuminosityBlocks ()
static bool	wantsStreamRuns ()
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Definition at line 26 of file PatVertexAnalyzer.cc.

Constructor & Destructor Documentation

PatVertexAnalyzer::PatVertexAnalyzer

(

const edm::ParameterSet &

params

)

constructor and destructor

Definition at line 48 of file PatVertexAnalyzer.cc.

    : srcToken_(consumes<reco::VertexCollection>(params.getParameter<edm::InputTag>("src"))),
      genParticlesToken_(consumes<reco::GenParticleCollection>(params.getParameter<edm::InputTag>("mc"))) {}

PatVertexAnalyzer::~PatVertexAnalyzer ( )

override

Definition at line 52 of file PatVertexAnalyzer.cc.

{}

Member Function Documentation

void PatVertexAnalyzer::analyze ( const edm::Event &  event, const edm::EventSetup &  es  )

override

Definition at line 75 of file PatVertexAnalyzer.cc.

References DEFINE_FWK_MODULE, genParticlesToken_, nTracks_, nVertices_, MetAnalyzer::pv(), srcToken_, reco::Vertex::tracksSize(), reco::Vertex::x(), x_, xDelta_, xErr_, reco::Vertex::xError(), xPull_, reco::Vertex::y(), y_, yDelta_, yErr_, reco::Vertex::yError(), yPull_, reco::Vertex::z(), z_, zDelta_, zErr_, reco::Vertex::zError(), and zPull_.

                                                                              {
  // handle to the primary vertex collection
  edm::Handle<reco::VertexCollection> pvHandle;
  event.getByToken(srcToken_, pvHandle);

  // handle to the generator particles (i.e. the MC truth)
  edm::Handle<reco::GenParticleCollection> genParticlesHandle;
  event.getByToken(genParticlesToken_, genParticlesHandle);

  // extract the position of the simulated vertex
  math::XYZPoint simPV = (*genParticlesHandle)[2].vertex();

  // the number of reconstructed primary vertices
  nVertices_->Fill(pvHandle->size());

  // if we have at least one, use the first (highest pt^2 sum)
  if (!pvHandle->empty()) {
    const reco::Vertex &pv = (*pvHandle)[0];

    nTracks_->Fill(pv.tracksSize());

    x_->Fill(pv.x());
    y_->Fill(pv.y());
    z_->Fill(pv.z());

    xErr_->Fill(pv.xError());
    yErr_->Fill(pv.yError());
    zErr_->Fill(pv.zError());

    xDelta_->Fill(pv.x() - simPV.X());
    yDelta_->Fill(pv.y() - simPV.Y());
    zDelta_->Fill(pv.z() - simPV.Z());

    xPull_->Fill((pv.x() - simPV.X()) / pv.xError());
    yPull_->Fill((pv.y() - simPV.Y()) / pv.yError());
    zPull_->Fill((pv.z() - simPV.Z()) / pv.zError());

    // we could access the tracks using the
    // pv.tracks_begin() ... pv.tracks_end() iterators
  }
}

void PatVertexAnalyzer::beginJob ( void  )

overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 54 of file PatVertexAnalyzer.cc.

References TFileService::make(), nTracks_, nVertices_, x_, xDelta_, xErr_, xPull_, y_, yDelta_, yErr_, yPull_, z_, zDelta_, zErr_, and zPull_.

                                 {
  // retrieve handle to auxiliary service
  //  used for storing histograms into ROOT file
  edm::Service<TFileService> fs;

  nVertices_ = fs->make<TH1F>("nVertices", "number of reconstructed primary vertices", 50, 0, 50);
  nTracks_ = fs->make<TH1F>("nTracks", "number of tracks at primary vertex", 100, 0, 300);
  x_ = fs->make<TH1F>("pvX", "primary vertex x", 100, -0.1, 0.1);
  y_ = fs->make<TH1F>("pvY", "primary vertex y", 100, -0.1, 0.1);
  z_ = fs->make<TH1F>("pvZ", "primary vertex z", 100, -30, 30);
  xErr_ = fs->make<TH1F>("pvErrorX", "primary vertex x error", 100, 0, 0.005);
  yErr_ = fs->make<TH1F>("pvErrorY", "primary vertex y error", 100, 0, 0.005);
  zErr_ = fs->make<TH1F>("pvErrorZ", "primary vertex z error", 100, 0, 0.01);
  xDelta_ = fs->make<TH1F>("pvDeltaX", "x shift wrt simulated vertex", 100, -0.01, 0.01);
  yDelta_ = fs->make<TH1F>("pvDeltaY", "y shift wrt simulated vertex", 100, -0.01, 0.01);
  zDelta_ = fs->make<TH1F>("pvDeltaZ", "z shift wrt simulated vertex", 100, -0.02, 0.02);
  xPull_ = fs->make<TH1F>("pvPullX", "primary vertex x pull", 100, -5, 5);
  yPull_ = fs->make<TH1F>("pvPullY", "primary vertex y pull", 100, -5, 5);
  zPull_ = fs->make<TH1F>("pvPullZ", "primary vertex z pull", 100, -5, 5);
}

Member Data Documentation

edm::EDGetTokenT<reco::GenParticleCollection> PatVertexAnalyzer::genParticlesToken_

private

Definition at line 39 of file PatVertexAnalyzer.cc.

Referenced by analyze().

TH1 * PatVertexAnalyzer::nTracks_

private

Definition at line 41 of file PatVertexAnalyzer.cc.

Referenced by analyze(), and beginJob().

TH1* PatVertexAnalyzer::nVertices_

private

Definition at line 41 of file PatVertexAnalyzer.cc.

Referenced by analyze(), and beginJob().

edm::EDGetTokenT<reco::VertexCollection> PatVertexAnalyzer::srcToken_

private

Definition at line 38 of file PatVertexAnalyzer.cc.

Referenced by analyze().

TH1* PatVertexAnalyzer::x_

private

Definition at line 42 of file PatVertexAnalyzer.cc.

Referenced by analyze(), and beginJob().

TH1* PatVertexAnalyzer::xDelta_

private

Definition at line 44 of file PatVertexAnalyzer.cc.

Referenced by analyze(), and beginJob().

TH1* PatVertexAnalyzer::xErr_

private

Definition at line 43 of file PatVertexAnalyzer.cc.

Referenced by analyze(), and beginJob().

TH1* PatVertexAnalyzer::xPull_

private

Definition at line 45 of file PatVertexAnalyzer.cc.

Referenced by analyze(), and beginJob().

TH1 * PatVertexAnalyzer::y_

private

Definition at line 42 of file PatVertexAnalyzer.cc.

Referenced by analyze(), and beginJob().

TH1 * PatVertexAnalyzer::yDelta_

private

Definition at line 44 of file PatVertexAnalyzer.cc.

Referenced by analyze(), and beginJob().

TH1 * PatVertexAnalyzer::yErr_

private

Definition at line 43 of file PatVertexAnalyzer.cc.

Referenced by analyze(), and beginJob().

TH1 * PatVertexAnalyzer::yPull_

private

Definition at line 45 of file PatVertexAnalyzer.cc.

Referenced by analyze(), and beginJob().

TH1 * PatVertexAnalyzer::z_

private

Definition at line 42 of file PatVertexAnalyzer.cc.

Referenced by analyze(), and beginJob().

TH1 * PatVertexAnalyzer::zDelta_

private

Definition at line 44 of file PatVertexAnalyzer.cc.

Referenced by analyze(), and beginJob().

TH1 * PatVertexAnalyzer::zErr_

private

Definition at line 43 of file PatVertexAnalyzer.cc.

Referenced by analyze(), and beginJob().

TH1 * PatVertexAnalyzer::zPull_

private

Definition at line 45 of file PatVertexAnalyzer.cc.

Referenced by analyze(), and beginJob().