PixelVertexProducerMedian Class Reference

#include <PixelVertexProducerMedian.h>

Inheritance diagram for PixelVertexProducerMedian:

Public Member Functions
	PixelVertexProducerMedian (const edm::ParameterSet &ps)
void	produce (edm::Event &ev, const edm::EventSetup &es) override
	~PixelVertexProducerMedian () override
Public Member Functions inherited from edm::EDProducer
	EDProducer ()
SerialTaskQueue *	globalLuminosityBlocksQueue ()
SerialTaskQueue *	globalRunsQueue ()
ModuleDescription const &	moduleDescription () const
	~EDProducer () override
Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const
	ProducerBase ()
std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const
std::vector< bool > const &	recordProvenanceList () const
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...
void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)
TypeLabelList const &	typeLabelList () const
	used by the fwk to register the list of products of this module More...
	~ProducerBase () noexcept(false) 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 &&)=default
	EDConsumerBase (EDConsumerBase const &)=delete
ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const
std::vector< ESProxyIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const
std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (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 &	operator= (EDConsumerBase &&)=default
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
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::ParameterSet	theConfig
double	thePtMin

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType
Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const *, const char *, edm::ProductResolverIndex > >
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
static bool	wantsGlobalLuminosityBlocks ()
static bool	wantsGlobalRuns ()
static bool	wantsInputProcessBlocks ()
static bool	wantsProcessBlocks ()
static bool	wantsStreamLuminosityBlocks ()
static bool	wantsStreamRuns ()
Protected Member Functions inherited from edm::ProducerBase
template<class ProductType >
BranchAliasSetterT< ProductType >	produces ()
	declare what type of product will make and with which optional label More...
template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces ()
template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces ()
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)
template<BranchType B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)
template<Transition B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)
template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)
template<class ProductType >
BranchAliasSetterT< ProductType >	produces (std::string instanceName)
template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)
ProducesCollector	producesCollector ()
Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept
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<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag) noexcept
template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

Detailed Description

Definition at line 12 of file PixelVertexProducerMedian.h.

Constructor & Destructor Documentation

PixelVertexProducerMedian()

PixelVertexProducerMedian::PixelVertexProducerMedian ( const edm::ParameterSet &  ps )

explicit

Definition at line 29 of file PixelVertexProducerMedian.cc.

                                                                              : theConfig(ps) {
  produces<reco::VertexCollection>();
}

~PixelVertexProducerMedian()

PixelVertexProducerMedian::~PixelVertexProducerMedian ( )

override

Definition at line 34 of file PixelVertexProducerMedian.cc.

{}

Member Function Documentation

produce()

void PixelVertexProducerMedian::produce ( edm::Event &  ev, const edm::EventSetup &  es  )

overridevirtual

Implements edm::EDProducer.

Definition at line 37 of file PixelVertexProducerMedian.cc.

                                                                             {
  // Get pixel tracks
  edm::Handle<reco::TrackCollection> trackCollection;
  std::string trackCollectionName = theConfig.getParameter<std::string>("TrackCollection");
  ev.getByLabel(trackCollectionName, trackCollection);
  const reco::TrackCollection tracks_ = *(trackCollection.product());
 
  thePtMin = theConfig.getParameter<double>("PtMin");
 
  // Select tracks
  std::vector<const reco::Track*> tracks;
  for (unsigned int i = 0; i < tracks_.size(); i++) {
    if (tracks_[i].pt() > thePtMin) {
      reco::TrackRef recTrack(trackCollection, i);
      tracks.push_back(&(*recTrack));
    }
  }
 
  LogTrace("MinBiasTracking") << " [VertexProducer] selected tracks: " << tracks.size() << " (out of " << tracks_.size()
                              << ")";
 
  auto vertices = std::make_unique<reco::VertexCollection>();
 
  if (!tracks.empty()) {
    // Sort along vertex z position
    std::sort(tracks.begin(), tracks.end(), ComparePairs());
 
    // Median
    float med;
    if (tracks.size() % 2 == 0)
      med = (tracks[tracks.size() / 2 - 1]->vz() + tracks[tracks.size() / 2]->vz()) / 2;
    else
      med = tracks[tracks.size() / 2]->vz();
 
    LogTrace("MinBiasTracking") << "  [vertex position] median    = " << med << " cm";
 
    if (tracks.size() > 10) {
      // Binning around med, halfWidth
      int nBin = 100;
      float halfWidth = 0.1;  // cm
 
      // Most probable
      TH1F histo("histo", "histo", nBin, -halfWidth, halfWidth);
 
      for (std::vector<const reco::Track*>::const_iterator track = tracks.begin(); track != tracks.end(); track++)
        if (fabs((*track)->vz() - med) < halfWidth)
          histo.Fill((*track)->vz() - med);
 
      LogTrace("MinBiasTracking") << "  [vertex position] most prob = "
                                  << med + histo.GetBinCenter(histo.GetMaximumBin()) << " cm";
 
      // Fit above max/2
      histo.Sumw2();
 
      TF1 f1("f1", "[0]*exp(-0.5 * ((x-[1])/[2])^2) + [3]");
      f1.SetParameters(10., 0., 0.01, 1.);
 
      histo.Fit(&f1, "QN");
 
      LogTrace("MinBiasTracking") << "  [vertex position] fitted    = " << med + f1.GetParameter(1) << " +- "
                                  << f1.GetParError(1) << " cm";
 
      // Store
      reco::Vertex::Error err;
      err(2, 2) = f1.GetParError(1) * f1.GetParError(1);
      reco::Vertex ver(reco::Vertex::Point(0, 0, med + f1.GetParameter(1)), err, 0, 1, 1);
      vertices->push_back(ver);
    } else {
      // Store
      reco::Vertex::Error err;
      err(2, 2) = 0.1 * 0.1;
      reco::Vertex ver(reco::Vertex::Point(0, 0, med), err, 0, 1, 1);
      vertices->push_back(ver);
    }
  }
  ev.put(std::move(vertices));
}

References submitPVResolutionJobs::err, ev, DeadROC_duringRun::f1, edm::ParameterSet::getParameter(), timingPdfMaker::histo, mps_fire::i, LogTrace, eostools::move(), PFElectronDQMAnalyzer_cfi::nBin, DiDispStaMuonMonitor_cfi::pt, AlCaHLTBitMon_QueryRunRegistry::string, tools::TF1, theConfig, thePtMin, HLT_FULL_cff::track, duplicaterechits_cfi::trackCollection, PDWG_EXOHSCP_cff::tracks, and pwdgSkimBPark_cfi::vertices.

Member Data Documentation

theConfig

edm::ParameterSet PixelVertexProducerMedian::theConfig

private

Definition at line 19 of file PixelVertexProducerMedian.h.

Referenced by produce().

thePtMin

double PixelVertexProducerMedian::thePtMin

private

Definition at line 20 of file PixelVertexProducerMedian.h.

Referenced by produce().