#include <HIBestVertexProducer.h>

Inheritance diagram for HIBestVertexProducer:

Public Member Functions
	HIBestVertexProducer (const edm::ParameterSet &ps)

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

	~HIBestVertexProducer ()

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducer ()

Public Member Functions inherited from edm::ProducerBase
	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

virtual	~ProducerBase ()

Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Member Functions
void	beginJob ()

Private Attributes
edm::EDGetTokenT < reco::VertexCollection >	theAdaptiveVertexCollection

edm::EDGetTokenT< reco::BeamSpot >	theBeamSpotTag

edm::ParameterSet	theConfig

edm::EDGetTokenT < reco::VertexCollection >	theMedianVertexCollection

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

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 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 12 of file HIBestVertexProducer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 17 of file HIBestVertexProducer.cc.

                             : theConfig(ps),
   theBeamSpotTag(consumes<reco::BeamSpot>(ps.getParameter<edm::InputTag>("beamSpotLabel"))),
   theMedianVertexCollection(consumes<reco::VertexCollection>(ps.getParameter<edm::InputTag>("medianVertexCollection"))),
   theAdaptiveVertexCollection(consumes<reco::VertexCollection>(ps.getParameter<edm::InputTag>("adaptiveVertexCollection")))
 {
   produces<reco::VertexCollection>();
 }

HIBestVertexProducer::~HIBestVertexProducer ( )

Definition at line 27 of file HIBestVertexProducer.cc.

28 {

29 }

Member Function Documentation

void HIBestVertexProducer::beginJob ( void )

privatevirtual

Reimplemented from edm::EDProducer.

Definition at line 32 of file HIBestVertexProducer.cc.

33 {

34 }

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

virtual

Implements edm::EDProducer.

Definition at line 38 of file HIBestVertexProducer.cc.

References SiPixelRawToDigiRegional_cfi::beamSpot, reco::BeamSpot::BeamWidthX(), reco::BeamSpot::BeamWidthY(), edm::Event::getByToken(), edm::HandleBase::isValid(), reco::BeamSpot::position(), funct::pow(), edm::Handle< T >::product(), edm::Event::put(), reco::BeamSpot::sigmaZ(), reco::BeamSpot::x0(), and reco::BeamSpot::y0().

 {
   
   // 1. use best adaptive vertex preferentially
   // 2. use median vertex in case where adaptive algorithm failed
   // 3. use beamspot if netither vertexing method succeeds
 
   // New vertex collection
   std::auto_ptr<reco::VertexCollection> newVertexCollection(new reco::VertexCollection);
 
   //** Get precise adaptive vertex **/
   edm::Handle<reco::VertexCollection> vc1;
   ev.getByToken(theAdaptiveVertexCollection, vc1);
   const reco::VertexCollection *vertices1 = vc1.product();
 
   if(vertices1->size()==0)
     LogError("HeavyIonVertexing") << "adaptive vertex collection is empty!" << endl;
 
   if(vertices1->begin()->zError()<3) { 
   
     reco::VertexCollection::const_iterator vertex1 = vertices1->begin();
     newVertexCollection->push_back(*vertex1);
 
     LogInfo("HeavyIonVertexing") << "adaptive vertex:\n vz = (" 
                  << vertex1->x() << ", " << vertex1->y() << ", " << vertex1->z() << ")" 
                  << "\n error = ("
                  << vertex1->xError() << ", " << vertex1->yError() << ", " 
                  << vertex1->zError() << ")" << endl;
   } else {
     
     //** Get fast median vertex **/
     edm::Handle<reco::VertexCollection> vc2;
     ev.getByToken(theMedianVertexCollection, vc2);
     const reco::VertexCollection * vertices2 = vc2.product();
     
     //** Get beam spot position and error **/
     reco::BeamSpot beamSpot;
     edm::Handle<reco::BeamSpot> beamSpotHandle;
     ev.getByToken(theBeamSpotTag, beamSpotHandle);
 
     if( beamSpotHandle.isValid() ) 
       beamSpot = *beamSpotHandle;
     else
       LogError("HeavyIonVertexing") << "no beamspot found "  << endl;
 
     if(vertices2->size() > 0) { 
       
       reco::VertexCollection::const_iterator vertex2 = vertices2->begin();
       reco::Vertex::Error err;
       err(0,0)=pow(beamSpot.BeamWidthX(),2);
       err(1,1)=pow(beamSpot.BeamWidthY(),2);
       err(2,2)=pow(vertex2->zError(),2);
       reco::Vertex newVertex(reco::Vertex::Point(beamSpot.x0(),beamSpot.y0(),vertex2->z()),
                  err, 0, 1, 1);
       newVertexCollection->push_back(newVertex);  
 
       LogInfo("HeavyIonVertexing") << "median vertex + beamspot: \n position = (" 
                    << newVertex.x() << ", " << newVertex.y() << ", " << newVertex.z() << ")" 
                    << "\n error = ("
                    << newVertex.xError() << ", " << newVertex.yError() << ", " 
                    << newVertex.zError() << ")" << endl;
       
     } else { 
       
       reco::Vertex::Error err;
       err(0,0)=pow(beamSpot.BeamWidthX(),2);
       err(1,1)=pow(beamSpot.BeamWidthY(),2);
       err(2,2)=pow(beamSpot.sigmaZ(),2);
       reco::Vertex newVertex(beamSpot.position(), 
                  err, 0, 0, 1);
       newVertexCollection->push_back(newVertex);  
 
       LogInfo("HeavyIonVertexing") << "beam spot: \n position = (" 
                    << newVertex.x() << ", " << newVertex.y() << ", " << newVertex.z() << ")" 
                    << "\n error = ("
                    << newVertex.xError() << ", " << newVertex.yError() << ", " 
                    << newVertex.zError() << ")" << endl;
       
     }
     
   }
   
   // put new vertex collection into event
   ev.put(newVertexCollection);
   
 }

Member Data Documentation

edm::EDGetTokenT<reco::VertexCollection> HIBestVertexProducer::theAdaptiveVertexCollection

private

Definition at line 24 of file HIBestVertexProducer.h.

edm::EDGetTokenT<reco::BeamSpot> HIBestVertexProducer::theBeamSpotTag

private

Definition at line 22 of file HIBestVertexProducer.h.

edm::ParameterSet HIBestVertexProducer::theConfig

private

Definition at line 21 of file HIBestVertexProducer.h.

edm::EDGetTokenT<reco::VertexCollection> HIBestVertexProducer::theMedianVertexCollection

private

Definition at line 23 of file HIBestVertexProducer.h.

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation