#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

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)

	~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 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::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::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 9 of file PixelVertexProducerMedian.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 35 of file PixelVertexProducerMedian.cc.

Referenced by ComparePairs::operator()().

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

PixelVertexProducerMedian::~PixelVertexProducerMedian ( )

override

Definition at line 42 of file PixelVertexProducerMedian.cc.

References produce().

43 {

44 }

Member Function Documentation

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

override

Definition at line 48 of file PixelVertexProducerMedian.cc.

References connectstrParser::f1, edm::Event::getByLabel(), trackerHits::histo, mps_fire::i, LogTrace, eostools::move(), edm::Handle< T >::product(), EnergyCorrector::pt, edm::Event::put(), jetUpdater_cfi::sort, AlCaHLTBitMon_QueryRunRegistry::string, HiIsolationCommonParameters_cff::track, findElectronsInSiStrips_cfi::trackCollection, l1t::tracks, tracks_, and electrons_cff::vertices.

Referenced by ~PixelVertexProducerMedian().

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

Member Data Documentation

edm::ParameterSet PixelVertexProducerMedian::theConfig

private

Definition at line 17 of file PixelVertexProducerMedian.h.

double PixelVertexProducerMedian::thePtMin

private

Definition at line 18 of file PixelVertexProducerMedian.h.

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation