#include <PixelVertexFinding/interface/PixelVertexProducer.h>

Inheritance diagram for PixelVertexProducer:

Public Member Functions
	PixelVertexProducer (const edm::ParameterSet &)

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

	~PixelVertexProducer ()

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 Attributes
DivisiveVertexFinder *	dvf_

bool	method2

double	ptMin_

edm::EDGetTokenT< reco::BeamSpot >	token_BeamSpot

edm::EDGetTokenT < reco::TrackCollection >	token_Tracks

edm::InputTag	trackCollName

int	verbose_

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

Description: This produces 1D (z only) primary vertexes using only pixel information.

Implementation: This producer can use either the Divisive Primary Vertex Finder or the Histogramming Primary Vertex Finder (currently not implemented). It relies on the PixelTripletProducer and PixelTrackProducer having already been run upstream. This is code ported from ORCA originally written by S Cucciarelli, M Konecki, D Kotlinski.

Definition at line 35 of file PixelVertexProducer.h.

Constructor & Destructor Documentation

PixelVertexProducer::PixelVertexProducer ( const edm::ParameterSet & conf )

explicit

Definition at line 13 of file PixelVertexProducer.cc.

References dvf_, edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), method2, ptMin_, AlCaHLTBitMon_QueryRunRegistry::string, token_BeamSpot, token_Tracks, trackCollName, and verbose_.

   : verbose_(0), dvf_(0), ptMin_(1.0)
 {
   // Register my product
   produces<reco::VertexCollection>();
 
   // Setup shop
   verbose_           = conf.getParameter<int>("Verbosity"); // 0 silent, 1 chatty, 2 loud
   std::string finder = conf.getParameter<std::string>("Finder"); // DivisiveVertexFinder
   bool useError      = conf.getParameter<bool>("UseError"); // true
   bool wtAverage     = conf.getParameter<bool>("WtAverage"); // true
   double zOffset     = conf.getParameter<double>("ZOffset"); // 5.0 sigma
   double zSeparation = conf.getParameter<double>("ZSeparation"); // 0.05 cm
   int ntrkMin        = conf.getParameter<int>("NTrkMin"); // 3
   // Tracking requirements before sending a track to be considered for vtx
   ptMin_ = conf.getParameter<double>("PtMin"); // 1.0 GeV
   trackCollName = conf.getParameter<edm::InputTag>("TrackCollection");
   token_Tracks = consumes<reco::TrackCollection>(trackCollName);
   token_BeamSpot = consumes<reco::BeamSpot>(conf.getParameter<edm::InputTag>("beamSpot"));
   method2 = conf.getParameter<bool>("Method2");
 
   double track_pt_min   = ptMin_;
   double track_pt_max   = 10.;
   double track_chi2_max = 9999999.;
   double track_prob_min = -1.;
 
   if ( conf.exists("PVcomparer") ) {
     edm::ParameterSet PVcomparerPSet = conf.getParameter<edm::ParameterSet>("PVcomparer");
     track_pt_min   = PVcomparerPSet.getParameter<double>("track_pt_min");    
     if (track_pt_min != ptMin_) {
       if (track_pt_min < ptMin_)
     edm::LogWarning("PixelVertexProducer") << "minimum track pT setting differs between PixelVertexProducer (" << ptMin_ << ") and PVcomparer (" << track_pt_min << ") [PVcomparer considers tracks w/ lower threshold than PixelVertexProducer does] !!!";
       else
     edm::LogInfo("PixelVertexProducer") << "minimum track pT setting differs between PixelVertexProducer (" << ptMin_ << ") and PVcomparer (" << track_pt_min << ") !!!";
     }
     track_pt_max   = PVcomparerPSet.getParameter<double>("track_pt_max");
     track_chi2_max = PVcomparerPSet.getParameter<double>("track_chi2_max");
     track_prob_min = PVcomparerPSet.getParameter<double>("track_prob_min");
   }
 
   if (finder == "DivisiveVertexFinder") {
     if (verbose_ > 0) edm::LogInfo("PixelVertexProducer") << ": Using the DivisiveVertexFinder\n";
     dvf_ = new DivisiveVertexFinder(track_pt_min,track_pt_max,track_chi2_max,track_prob_min,zOffset, ntrkMin, useError, zSeparation, wtAverage, verbose_);
   }
   else { // Finder not supported, or you made a mistake in your request
     // throw an exception once I figure out how CMSSW does this
   }
 }

PixelVertexProducer::~PixelVertexProducer ( )

Definition at line 63 of file PixelVertexProducer.cc.

References dvf_.

                                           {
   delete dvf_;
 }

Member Function Documentation

void PixelVertexProducer::produce	(	edm::Event &	e,
		const edm::EventSetup &	es
	)

virtual

Implements edm::EDProducer.

Definition at line 67 of file PixelVertexProducer.cc.

References reco::Vertex::add(), dvf_, reco::BeamSpot::dxdz(), reco::BeamSpot::dydz(), relativeConstraints::error, DivisiveVertexFinder::findVertexes(), DivisiveVertexFinder::findVertexesAlt(), edm::Event::getByToken(), i, edm::HandleBase::isValid(), method2, convertSQLiteXML::ok, reco::BeamSpot::position(), edm::Handle< T >::product(), RecoTauCleanerPlugins::pt, ptMin_, edm::RefVector< C, T, F >::push_back(), edm::Event::put(), reco::BeamSpot::rotatedCovariance3D(), edm::RefVector< C, T, F >::size(), mathSSE::sqrt(), token_BeamSpot, token_Tracks, trackCollName, testEve_cfg::tracks, findQualityFiles::v, verbose_, x, reco::BeamSpot::x0(), detailsBasic3DVector::y, reco::BeamSpot::y0(), detailsBasic3DVector::z, and reco::BeamSpot::z0().

                                                                       {
 
   // First fish the pixel tracks out of the event
   edm::Handle<reco::TrackCollection> trackCollection;
   e.getByToken(token_Tracks,trackCollection);
   const reco::TrackCollection tracks = *(trackCollection.product());
   if (verbose_ > 0) edm::LogInfo("PixelVertexProducer") << ": Found " << tracks.size() << " tracks in TrackCollection called " << trackCollName << "\n";
   
 
   // Second, make a collection of pointers to the tracks we want for the vertex finder
   reco::TrackRefVector trks;
   for (unsigned int i=0; i<tracks.size(); i++) {
     if (tracks[i].pt() > ptMin_)     
       trks.push_back( reco::TrackRef(trackCollection, i) );
   }
   if (verbose_ > 0) edm::LogInfo("PixelVertexProducer") << ": Selected " << trks.size() << " of these tracks for vertexing\n";
 
   edm::Handle<reco::BeamSpot> bsHandle;
   e.getByToken(token_BeamSpot,bsHandle);
   math::XYZPoint myPoint(0.,0.,0.);
   if (bsHandle.isValid()) myPoint = math::XYZPoint(bsHandle->x0(),bsHandle->y0(), 0. ); //FIXME: fix last coordinate with vertex.z() at same time
 
   // Third, ship these tracks off to be vertexed
   std::auto_ptr<reco::VertexCollection> vertexes(new reco::VertexCollection);
   bool ok;
   if (method2) {
     ok = dvf_->findVertexesAlt(trks,       // input
                  *vertexes,myPoint); // output
     if (verbose_ > 0) edm::LogInfo("PixelVertexProducer") << "Method2 returned status of " << ok;
   }
   else {
     ok = dvf_->findVertexes(trks,       // input
                 *vertexes); // output
     if (verbose_ > 0) edm::LogInfo("PixelVertexProducer") << "Method1 returned status of " << ok;
   }
 
   if (verbose_ > 0) {
     edm::LogInfo("PixelVertexProducer") << ": Found " << vertexes->size() << " vertexes\n";
     for (unsigned int i=0; i<vertexes->size(); ++i) {
       edm::LogInfo("PixelVertexProducer") << "Vertex number " << i << " has " << (*vertexes)[i].tracksSize() << " tracks with a position of " << (*vertexes)[i].z() << " +- " << std::sqrt( (*vertexes)[i].covariance(2,2) );
     }
   }
 
 
   if(bsHandle.isValid())
     {
       const reco::BeamSpot & bs = *bsHandle;
       
       for (unsigned int i=0; i<vertexes->size(); ++i) {
     double z=(*vertexes)[i].z();
     double x=bs.x0()+bs.dxdz()*(z-bs.z0());
     double y=bs.y0()+bs.dydz()*(z-bs.z0()); 
     reco::Vertex v( reco::Vertex::Point(x,y,z), (*vertexes)[i].error(),(*vertexes)[i].chi2() , (*vertexes)[i].ndof() , (*vertexes)[i].tracksSize());
     //Copy also the tracks 
     for (std::vector<reco::TrackBaseRef >::const_iterator it = (*vertexes)[i].tracks_begin();
          it !=(*vertexes)[i].tracks_end(); it++) {
       v.add( *it );
     }
     (*vertexes)[i]=v;
     
       }
     }
   else
     {
       edm::LogWarning("PixelVertexProducer") << "No beamspot found. Using returning vertexes with (0,0,Z) ";
     } 
   
   if(!vertexes->size() && bsHandle.isValid()){
     
     const reco::BeamSpot & bs = *bsHandle;
       
       GlobalError bse(bs.rotatedCovariance3D());
       if ( (bse.cxx() <= 0.) ||
        (bse.cyy() <= 0.) ||
        (bse.czz() <= 0.) ) {
     AlgebraicSymMatrix33 we;
     we(0,0)=10000;
     we(1,1)=10000;
     we(2,2)=10000;
     vertexes->push_back(reco::Vertex(bs.position(), we,0.,0.,0));
     
     edm::LogInfo("PixelVertexProducer") <<"No vertices found. Beamspot with invalid errors " << bse.matrix() << std::endl
                            << "Will put Vertex derived from dummy-fake BeamSpot into Event.\n"
                            << (*vertexes)[0].x() << "\n"
                            << (*vertexes)[0].y() << "\n"
                            << (*vertexes)[0].z() << "\n";
       } else {
     vertexes->push_back(reco::Vertex(bs.position(),
                      bs.rotatedCovariance3D(),0.,0.,0));
     
     edm::LogInfo("PixelVertexProducer") << "No vertices found. Will put Vertex derived from BeamSpot into Event:\n"
                            << (*vertexes)[0].x() << "\n"
                            << (*vertexes)[0].y() << "\n"
                            << (*vertexes)[0].z() << "\n";
       }
   }
       
   else if(!vertexes->size() && !bsHandle.isValid())
     {
       edm::LogWarning("PixelVertexProducer") << "No beamspot and no vertex found. No vertex returned.";
     }
   
   e.put(vertexes);
   
 }