#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

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	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>
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)

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.

34 {}

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));
 }