HIProtoTrackFilter Class Reference

#include <HIProtoTrackFilter.h>

Inheritance diagram for HIProtoTrackFilter:

Public Member Functions
	HIProtoTrackFilter (const edm::ParameterSet &ps, const edm::EventSetup &es)
	HIProtoTrackFilter (const edm::ParameterSet &ps)
virtual bool	operator() (const reco::Track *, const PixelTrackFilter::Hits &hits) const
virtual void	update (edm::Event &ev)
virtual	~HIProtoTrackFilter ()
Public Member Functions inherited from PixelTrackFilter
virtual bool	operator() (const reco::Track *) const
virtual	~PixelTrackFilter ()

Private Attributes
bool	doVariablePtMin
const reco::BeamSpot *	theBeamSpot
edm::InputTag	theBeamSpotTag
double	theChi2Max
double	thePtMin
edm::InputTag	theSiPixelRecHits
double	theTIPMax
double	theVariablePtMin

Additional Inherited Members
Public Types inherited from PixelTrackFilter
typedef std::vector< const TrackingRecHit * >	Hits

Detailed Description

Definition at line 11 of file HIProtoTrackFilter.h.

Constructor & Destructor Documentation

HIProtoTrackFilter::HIProtoTrackFilter	(	const edm::ParameterSet &	ps,
		const edm::EventSetup &	es
	)

Definition at line 23 of file HIProtoTrackFilter.cc.

                                                                                          :
theTIPMax( ps.getParameter<double>("tipMax") ),
theChi2Max( ps.getParameter<double>("chi2") ),
thePtMin( ps.getParameter<double>("ptMin") ),
doVariablePtMin( ps.getParameter<bool>("doVariablePtMin") ),
theBeamSpotTag( ps.getParameter<InputTag>("beamSpot")),
theSiPixelRecHits( ps.getParameter<InputTag>("siPixelRecHits")),
theBeamSpot(0),
theVariablePtMin(0)
{ 
}

HIProtoTrackFilter::HIProtoTrackFilter

(

const edm::ParameterSet &

ps

)

Definition at line 36 of file HIProtoTrackFilter.cc.

                                                                 :
theTIPMax( ps.getParameter<double>("tipMax") ),
theChi2Max( ps.getParameter<double>("chi2") ),
thePtMin( ps.getParameter<double>("ptMin") ),
doVariablePtMin( ps.getParameter<bool>("doVariablePtMin") ),
theBeamSpotTag( ps.getParameter<InputTag>("beamSpot")),
theSiPixelRecHits( ps.getParameter<InputTag>("siPixelRecHits")),
theBeamSpot(0),
theVariablePtMin(0)
{ 
}

HIProtoTrackFilter::~HIProtoTrackFilter ( )

virtual

Definition at line 49 of file HIProtoTrackFilter.cc.

{ }

Member Function Documentation

bool HIProtoTrackFilter::operator() ( const reco::Track *  track, const PixelTrackFilter::Hits &  hits  ) const

virtual

Reimplemented from PixelTrackFilter.

Definition at line 53 of file HIProtoTrackFilter.cc.

References reco::TrackBase::chi2(), doVariablePtMin, reco::TrackBase::dxy(), reco::BeamSpot::position(), reco::TrackBase::pt(), theBeamSpot, theChi2Max, thePtMin, theTIPMax, and theVariablePtMin.

{

  if (!track) return false; 

  float minpt = thePtMin;
  if(doVariablePtMin) minpt = theVariablePtMin;

  if (track->chi2() > theChi2Max || track->pt() < minpt) return false; 
  
  math::XYZPoint vtxPoint(0.0,0.0,0.0);
  
  if(theBeamSpot)
    vtxPoint = theBeamSpot->position();
  
  double d0=0.0;
  d0 = -1.*track->dxy(vtxPoint);
  
  if (theTIPMax>0 && fabs(d0)>theTIPMax) return false;
  
  return true;
}

void HIProtoTrackFilter::update ( edm::Event &  ev )

virtual

Reimplemented from PixelTrackFilter.

Definition at line 77 of file HIProtoTrackFilter.cc.

References edmNew::copyDetSetRange(), edm::Event::getByLabel(), LogTrace, TrackerLayerIdAccessor::pixelBarrelLayer(), edm::Handle< T >::product(), theBeamSpot, theBeamSpotTag, thePtMin, theSiPixelRecHits, theVariablePtMin, reco::BeamSpot::x0(), reco::BeamSpot::y0(), and reco::BeamSpot::z0().

Referenced by progressbar.ProgressBar::__next__(), relval_steps.Matrix::__setitem__(), relval_steps.Steps::__setitem__(), Vispa.Gui.VispaWidget.VispaWidget::autosize(), Vispa.Views.LineDecayView.LineDecayContainer::createObject(), Vispa.Views.LineDecayView.LineDecayContainer::deselectAllObjects(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::deselectAllWidgets(), Vispa.Gui.VispaWidget.VispaWidget::enableAutosizing(), progressbar.ProgressBar::finish(), Vispa.Gui.MenuWidget.MenuWidget::leaveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseMoveEvent(), Vispa.Gui.MenuWidget.MenuWidget::mouseMoveEvent(), Vispa.Views.LineDecayView.LineDecayContainer::mouseMoveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseReleaseEvent(), Vispa.Views.LineDecayView.LineDecayContainer::objectMoved(), relval_steps.Steps::overwrite(), Vispa.Views.LineDecayView.LineDecayContainer::removeObject(), Vispa.Gui.ConnectableWidget.ConnectableWidget::removePorts(), Vispa.Gui.FindDialog.FindDialog::reset(), Vispa.Gui.PortConnection.PointToPointConnection::select(), Vispa.Gui.VispaWidget.VispaWidget::select(), Vispa.Views.LineDecayView.LineDecayContainer::select(), Vispa.Gui.VispaWidget.VispaWidget::setText(), Vispa.Gui.VispaWidget.VispaWidget::setTitle(), Vispa.Gui.ZoomableWidget.ZoomableWidget::setZoom(), Vispa.Views.LineDecayView.LineDecayContainer::setZoom(), and Vispa.Gui.PortConnection.PointToPointConnection::updateConnection().

{
  
  // Get the beam spot
  edm::Handle<reco::BeamSpot> bsHandle;
  ev.getByLabel( theBeamSpotTag, bsHandle);
  theBeamSpot = bsHandle.product();
  
  if(theBeamSpot) {
    LogInfo("HeavyIonVertexing") 
      << "[HIProtoTrackFilter] Proto track selection based on beamspot"
      << "\n   (x,y,z) = (" << theBeamSpot->x0() << "," << theBeamSpot->y0() << "," << theBeamSpot->z0() << ")";
  } else {
    LogError("HeavyIonVertexing") // this can be made a warning when operator() is fixed
      << "No beamspot found with tag '" << theBeamSpotTag << "'";
  }

  // Estimate multiplicity
  edm::Handle<SiPixelRecHitCollection> recHitColl;
  ev.getByLabel(theSiPixelRecHits, recHitColl);
  
  vector<const TrackingRecHit*> theChosenHits;   
  TrackerLayerIdAccessor acc;    
  edmNew::copyDetSetRange(*recHitColl,theChosenHits,acc.pixelBarrelLayer(1));    
  float estMult = theChosenHits.size();
  
  theVariablePtMin=thePtMin;

  // parameterize ptMin such that a roughly constant number of selected prototracks passed are to vertexing
  float varPtCutoff = 1500; //cutoff for variable ptMin
  if(estMult < varPtCutoff) {
    theVariablePtMin = 0.075;
    if(estMult > 0) theVariablePtMin = (13. - (varPtCutoff/estMult) )/12.; 
    if(theVariablePtMin<0.075) theVariablePtMin = 0.075; // don't lower the cut past 75 MeV
  }
  
  LogTrace("heavyIonHLTVertexing")<<"   [HIProtoTrackFilter: theVariablePtMin: " << theVariablePtMin << "]";
  
  
  return;
  
}

Member Data Documentation

bool HIProtoTrackFilter::doVariablePtMin

private

Definition at line 21 of file HIProtoTrackFilter.h.

Referenced by operator()().

const reco::BeamSpot* HIProtoTrackFilter::theBeamSpot

private

Definition at line 24 of file HIProtoTrackFilter.h.

Referenced by operator()(), and update().

edm::InputTag HIProtoTrackFilter::theBeamSpotTag

private

Definition at line 22 of file HIProtoTrackFilter.h.

Referenced by update().

double HIProtoTrackFilter::theChi2Max

private

Definition at line 20 of file HIProtoTrackFilter.h.

Referenced by operator()().

double HIProtoTrackFilter::thePtMin

private

Definition at line 20 of file HIProtoTrackFilter.h.

Referenced by operator()(), and update().

edm::InputTag HIProtoTrackFilter::theSiPixelRecHits

private

Definition at line 23 of file HIProtoTrackFilter.h.

Referenced by update().

double HIProtoTrackFilter::theTIPMax

private

Definition at line 19 of file HIProtoTrackFilter.h.

Referenced by operator()().

double HIProtoTrackFilter::theVariablePtMin

private

Definition at line 25 of file HIProtoTrackFilter.h.

Referenced by operator()(), and update().