#include <PixelVertexProducerClusters.h>

Inheritance diagram for PixelVertexProducerClusters:

Public Member Functions
int	getContainedHits (const std::vector< VertexHit > &hits, float z0, float &chi)

	PixelVertexProducerClusters (const edm::ParameterSet &ps)

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

	~PixelVertexProducerClusters ()

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
std::vector< ConsumesInfo >	consumesInfo () const

	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

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) 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	beginRun (edm::Run const &run, edm::EventSetup const &es) override

Private Attributes
edm::ParameterSet	theConfig

const TrackerGeometry *	theTracker

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 PixelVertexProducerClusters.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 48 of file PixelVertexProducerClusters.cc.

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

PixelVertexProducerClusters::~PixelVertexProducerClusters ( )

Definition at line 56 of file PixelVertexProducerClusters.cc.

57 {

58 }

Member Function Documentation

void PixelVertexProducerClusters::beginRun	(	edm::Run const &	run,
		edm::EventSetup const &	es
	)

overrideprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 62 of file PixelVertexProducerClusters.cc.

References edm::EventSetup::get(), and edm::ESHandle< class >::product().

 {
   // Get tracker geometry
   edm::ESHandle<TrackerGeometry> trackerHandle;
   es.get<TrackerDigiGeometryRecord>().get(trackerHandle);
   theTracker = trackerHandle.product();
 }

int PixelVertexProducerClusters::getContainedHits	(	const std::vector< VertexHit > &	hits,
		float	z0,
		float &	chi
	)

Definition at line 72 of file PixelVertexProducerClusters.cc.

References gen::n, AlCaHLTBitMon_ParallelJobs::p, and hit::z.

 {
   int n = 0;
   chi = 0.;
 
   for(vector<VertexHit>::const_iterator hit = hits.begin();
                                         hit!= hits.end(); hit++)
   {
     // Predicted cluster width in y direction
     float p = 2 * fabs(hit->z - z0)/hit->r + 0.5; // FIXME
 
     if(fabs(p - hit->w) <= 1.)
     { 
       chi += fabs(p - hit->w);
       n++;
     }
   }
 
   return n;
 }

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

overridevirtual

Implements edm::EDProducer.

Definition at line 95 of file PixelVertexProducerClusters.cc.

References edmNew::DetSetVector< T >::data(), alignCSCRings::e, edm::EventSetup::get(), edm::Event::getByLabel(), PixelTopology::isItEdgePixelInX(), PixelTopology::isItEdgePixelInY(), LogTrace, PixelSubdetector::PixelBarrel, edm::Handle< T >::product(), edm::Event::put(), VertexHit::r, edmNew::DetSetVector< T >::size(), PixelGeomDetUnit::specificTopology(), HLT_25ns14e33_v1_cff::vertices, VertexHit::w, and VertexHit::z.

 {
   //Retrieve tracker topology from geometry
   edm::ESHandle<TrackerTopology> tTopo;
   es.get<IdealGeometryRecord>().get(tTopo);
 
 
   // Get pixel hit collections
   edm::Handle<SiPixelRecHitCollection> pixelColl;
   ev.getByLabel("siPixelRecHits",      pixelColl);
 
   const SiPixelRecHitCollection* thePixelHits = pixelColl.product();
 
   std::auto_ptr<reco::VertexCollection> vertices(new reco::VertexCollection);
 
   if(thePixelHits->size() > 0)
   {
     vector<VertexHit> hits;
 
     for(SiPixelRecHitCollection::DataContainer::const_iterator
            recHit = thePixelHits->data().begin(),
            recHitEnd = thePixelHits->data().end();
            recHit != recHitEnd; ++recHit)
     {
       if(recHit->isValid())
       {
 //      if(!(recHit->isOnEdge() || recHit->hasBadPixels()))
         DetId id = recHit->geographicalId();
         const PixelGeomDetUnit* pgdu =
           dynamic_cast<const PixelGeomDetUnit*>(theTracker->idToDetUnit(id));
         const PixelTopology* theTopol = ( &(pgdu->specificTopology()) );
         vector<SiPixelCluster::Pixel> pixels = recHit->cluster()->pixels();
 
         bool pixelOnEdge = false;
         for(vector<SiPixelCluster::Pixel>::const_iterator
                pixel = pixels.begin(); pixel!= pixels.end(); pixel++)
         {
           int pos_x = (int)pixel->x;
           int pos_y = (int)pixel->y;
 
           if(theTopol->isItEdgePixelInX(pos_x) ||
               theTopol->isItEdgePixelInY(pos_y))
           { pixelOnEdge = true; break; }
         }
     
         if(!pixelOnEdge)
         if(id.subdetId() == int(PixelSubdetector::PixelBarrel))  
         {
           
  
           LocalPoint lpos = LocalPoint(recHit->localPosition().x(),
                                        recHit->localPosition().y(),
                                        recHit->localPosition().z());
 
           GlobalPoint gpos = theTracker->idToDet(id)->toGlobal(lpos);
 
           VertexHit hit;
           hit.z = gpos.z(); 
           hit.r = gpos.perp(); 
           hit.w = recHit->cluster()->sizeY();
 
           hits.push_back(hit);
         }
       }
     }
 
     int nhits; int nhits_max = 0;
     float chi; float chi_max = 1e+9;
 
     float zest = 0;
  
     for(float z0 = -15.9; z0 <= 15.95; z0 += 0.1)
     {
       nhits = getContainedHits(hits, z0, chi);
 
       if(nhits > 0)
       {
         if(nhits >  nhits_max)
         { chi_max = 1e+9; nhits_max = nhits; }
 
         if(nhits >= nhits_max)
         if(chi < chi_max)
         { chi_max = chi; zest = z0; }
       }
     }
 
     LogTrace("MinBiasTracking")
       << "  [vertex position] estimated = " << zest 
       << " | pixel barrel hits = " << thePixelHits->size();
 
     reco::Vertex::Error err;
     err(2,2) = 0.6 * 0.6;
     reco::Vertex ver(reco::Vertex::Point(0,0,zest), err, 0, 1, 1);
     vertices->push_back(ver);
   }
 
   ev.put(vertices);
 }

Member Data Documentation

edm::ParameterSet PixelVertexProducerClusters::theConfig

private

Definition at line 23 of file PixelVertexProducerClusters.h.

const TrackerGeometry* PixelVertexProducerClusters::theTracker

private

Definition at line 25 of file PixelVertexProducerClusters.h.

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation