ConversionSeedFinder Class Reference

#include <ConversionSeedFinder.h>

Inheritance diagram for ConversionSeedFinder:

Public Member Functions
void	clear ()
	ConversionSeedFinder ()
	ConversionSeedFinder (const edm::ParameterSet &config, edm::ConsumesCollector &iC)
const MeasurementTracker *	getMeasurementTracker () const
std::vector< const DetLayer * > const &	layerList () const
virtual void	makeSeeds (const edm::Handle< edm::View< reco::CaloCluster > > &allBc)=0
TrajectorySeedCollection &	seeds ()
virtual void	setCandidate (float e, GlobalPoint pos)
void	setEvent (const edm::Event &e)
void	setEventSetup (const edm::EventSetup &es)
	Initialize EventSetup objects at each event. More...
void	setMeasurementTracker (const MeasurementTracker *tracker) const
void	setNavigationSchool (const NavigationSchool *navigation)
virtual	~ConversionSeedFinder ()

Protected Member Functions
void	findLayers ()
void	findLayers (const FreeTrajectoryState &fts)
void	printLayer (int i) const
FreeTrajectoryState	trackStateFromClusters (int aCharge, const GlobalPoint &gpOrigine, PropagationDirection dir, float scaleFactor) const

Protected Attributes
edm::EDGetTokenT< reco::BeamSpot >	beamSpotToken_
PropagationDirection	dir_
edm::EDGetTokenT< MeasurementTrackerEvent >	measurementTrkToken_
float	theBCEnergy_
GlobalPoint	theBCPosition_
reco::BeamSpot	theBeamSpot_
edm::ESHandle< GeometricSearchTracker >	theGeomSearchTracker_
edm::ESGetToken< GeometricSearchTracker, TrackerRecoGeometryRecord >	theGeomSearchTrackerToken_
std::vector< const DetLayer * >	theLayerList_
const MeasurementTracker *	theMeasurementTracker_
edm::ESGetToken< MeasurementTracker, CkfComponentsRecord >	theMeasurementTrackerToken_
edm::ESHandle< MagneticField >	theMF_
edm::ESGetToken< MagneticField, IdealMagneticFieldRecord >	theMFToken_
const NavigationSchool *	theNavigationSchool_ = nullptr
const Propagator *	thePropagatorAlongMomentum_
edm::ESGetToken< Propagator, TrackingComponentsRecord >	thePropagatorAlongMomentumToken_
const Propagator *	thePropagatorOppositeToMomentum_
edm::ESGetToken< Propagator, TrackingComponentsRecord >	thePropagatorOppositeToMomentumToken_
reco::CaloCluster *	theSC_
float	theSCenergy_
GlobalPoint	theSCPosition_
TrajectorySeedCollection	theSeeds_
edm::Handle< MeasurementTrackerEvent >	theTrackerData_
const TrackingGeometry *	theTrackerGeom_
KFUpdator	theUpdator_

Detailed Description

Author

Nancy Marinelli, U. of Notre Dame, US

Definition at line 52 of file ConversionSeedFinder.h.

Constructor & Destructor Documentation

ConversionSeedFinder() [1/2]

ConversionSeedFinder::ConversionSeedFinder

(

)

ConversionSeedFinder() [2/2]

ConversionSeedFinder::ConversionSeedFinder	(	const edm::ParameterSet &	config,
		edm::ConsumesCollector &	iC
	)

Definition at line 16 of file ConversionSeedFinder.cc.

References beamSpotToken_, edm::ConsumesCollector::consumes(), edm::ConsumesCollector::esConsumes(), ProducerED_cfi::InputTag, LogDebug, HLT_2023v12_cff::measurementTrackerName, measurementTrkToken_, AlCaHLTBitMon_QueryRunRegistry::string, theGeomSearchTrackerToken_, theMeasurementTrackerToken_, theMFToken_, thePropagatorAlongMomentumToken_, and thePropagatorOppositeToMomentumToken_.

    :  //  conf_(config),
      theUpdator_() {
  measurementTrkToken_ = iC.consumes<MeasurementTrackerEvent>(
      edm::InputTag("MeasurementTrackerEvent"));  //hardcoded because the original was and no time to fix (sigh)
  beamSpotToken_ = iC.consumes<reco::BeamSpot>(
      edm::InputTag("offlineBeamSpot"));  //hardcoded because the original was and no time to fix (sigh)

  theMFToken_ = iC.esConsumes();
  theGeomSearchTrackerToken_ = iC.esConsumes();
  thePropagatorAlongMomentumToken_ = iC.esConsumes(edm::ESInputTag("", "alongMomElePropagator"));
  thePropagatorOppositeToMomentumToken_ = iC.esConsumes(edm::ESInputTag("", "oppositeToMomElePropagator"));
  LogDebug("ConversionSeedFinder") << " CTOR "
                                   << "\n";

  auto measurementTrackerName = config.getParameter<std::string>("MeasurementTrackerName");
  theMeasurementTrackerToken_ = iC.esConsumes(edm::ESInputTag("", measurementTrackerName));
}

~ConversionSeedFinder()

virtual ConversionSeedFinder::~ConversionSeedFinder ( )

inlinevirtual

Definition at line 57 of file ConversionSeedFinder.h.

{}

Member Function Documentation

clear()

void ConversionSeedFinder::clear ( void  )

inline

Definition at line 76 of file ConversionSeedFinder.h.

References theSeeds_.

Referenced by ConversionTrackCandidateProducer::produce().

{ theSeeds_.clear(); }

findLayers() [1/2]

void ConversionSeedFinder::findLayers ( )

protected

Definition at line 57 of file ConversionSeedFinder.cc.

References alongMomentum, ALCARECOTkAlJpsiMuMu_cff::charge, reco::BeamSpot::position(), theBeamSpot_, trackStateFromClusters(), and bphysicsOniaDQM_cfi::vertex.

Referenced by OutInConversionSeedFinder::makeSeeds(), and InOutConversionSeedFinder::makeSeeds().

                                      {
  int charge;
  //List the DetLayers crossed by a straight line from the centre of the
  //detector to the supercluster position
  //  GlobalPoint  vertex(0.,0.,0.);
  GlobalPoint vertex(theBeamSpot_.position().x(), theBeamSpot_.position().y(), theBeamSpot_.position().z());
  charge = -1;
  FreeTrajectoryState theStraightLineFTS = trackStateFromClusters(charge, vertex, alongMomentum, 1.);

  findLayers(theStraightLineFTS);
}

findLayers() [2/2]

void ConversionSeedFinder::findLayers ( const FreeTrajectoryState &  fts )

protected

Definition at line 108 of file ConversionSeedFinder.cc.

References LayerCollector::allLayers(), alongMomentum, getMeasurementTracker(), mps_fire::i, printLayer(), theLayerList_, theMF_, and theNavigationSchool_.

                                                                     {
  theLayerList_.clear();

  StraightLinePropagator prop(&(*theMF_), alongMomentum);

  LayerCollector collector(theNavigationSchool_, &prop, this->getMeasurementTracker(), 5., 5.);

  theLayerList_ = collector.allLayers(traj);

  for (unsigned int i = 0; i < theLayerList_.size(); ++i) {
    printLayer(i);
  }
}

getMeasurementTracker()

const MeasurementTracker* ConversionSeedFinder::getMeasurementTracker ( ) const

inline

Definition at line 69 of file ConversionSeedFinder.h.

References theMeasurementTracker_.

Referenced by OutInConversionSeedFinder::completeSeed(), InOutConversionSeedFinder::completeSeed(), findLayers(), InOutConversionSeedFinder::findSeeds(), setEvent(), and OutInConversionSeedFinder::startSeed().

{ return theMeasurementTracker_; }

layerList()

std::vector<const DetLayer*> const& ConversionSeedFinder::layerList ( ) const

inline

Definition at line 66 of file ConversionSeedFinder.h.

References theLayerList_.

Referenced by InOutConversionSeedFinder::completeSeed(), InOutConversionSeedFinder::fillClusterSeeds(), InOutConversionSeedFinder::findSeeds(), and OutInConversionSeedFinder::startSeed().

{ return theLayerList_; }

makeSeeds()

virtual void ConversionSeedFinder::makeSeeds ( const edm::Handle< edm::View< reco::CaloCluster > > &  allBc )

pure virtual

Implemented in InOutConversionSeedFinder, and OutInConversionSeedFinder.

printLayer()

void ConversionSeedFinder::printLayer ( int  i ) const

protected

Definition at line 122 of file ConversionSeedFinder.cc.

References GeomDetEnumerators::barrel, mps_fire::i, nano_mu_digi_cff::layer, and theLayerList_.

Referenced by InOutConversionSeedFinder::completeSeed(), InOutConversionSeedFinder::fillClusterSeeds(), and findLayers().

                                                 {
  const DetLayer* layer = theLayerList_[i];
  if (layer->location() == GeomDetEnumerators::barrel) {
    //    const BarrelDetLayer * barrelLayer = dynamic_cast<const BarrelDetLayer*>(layer);
    //float r = barrelLayer->specificSurface().radius();
    //    std::cout   <<  " barrel layer radius " << r << " " << barrelLayer->specificSurface().bounds().length()/2. << "\n";

  } else {
    //    const ForwardDetLayer * forwardLayer = dynamic_cast<const ForwardDetLayer*>(layer);
    // float z =  fabs(forwardLayer->surface().position().z());
    //    std::cout   << " forward layer position " << z << " " << forwardLayer->specificSurface().innerRadius() << " " << forwardLayer->specificSurface().outerRadius() << "\n";
  }
}

seeds()

TrajectorySeedCollection& ConversionSeedFinder::seeds ( )

inline

Definition at line 61 of file ConversionSeedFinder.h.

References theSeeds_.

Referenced by ConversionTrackCandidateProducer::buildCollections().

{ return theSeeds_; }

setCandidate()

virtual void ConversionSeedFinder::setCandidate ( float  e, GlobalPoint  pos  )

inlinevirtual

Definition at line 62 of file ConversionSeedFinder.h.

References MillePedeFileConverter_cfg::e, theSCenergy_, and theSCPosition_.

Referenced by ConversionTrackCandidateProducer::buildCollections().

                                                      {
    theSCenergy_ = e;
    theSCPosition_ = pos;
  }

setEvent()

void ConversionSeedFinder::setEvent

(

const edm::Event &

e

)

Definition at line 35 of file ConversionSeedFinder.cc.

References beamSpotToken_, MeasurementTracker::geomTracker(), edm::Event::getByToken(), getMeasurementTracker(), measurementTrkToken_, theBeamSpot_, theTrackerData_, and theTrackerGeom_.

Referenced by ConversionTrackCandidateProducer::produce().

                                                       {
  theTrackerGeom_ = this->getMeasurementTracker()->geomTracker();

  //get the BeamSpot
  edm::Handle<reco::BeamSpot> recoBeamSpotHandle;
  evt.getByToken(beamSpotToken_, recoBeamSpotHandle);
  theBeamSpot_ = *recoBeamSpotHandle;

  evt.getByToken(measurementTrkToken_, theTrackerData_);
}

setEventSetup()

void ConversionSeedFinder::setEventSetup

(

const edm::EventSetup &

es

)

Initialize EventSetup objects at each event.

Definition at line 46 of file ConversionSeedFinder.cc.

References edm::EventSetup::getData(), edm::EventSetup::getHandle(), theGeomSearchTracker_, theGeomSearchTrackerToken_, theMeasurementTracker_, theMeasurementTrackerToken_, theMF_, theMFToken_, thePropagatorAlongMomentum_, thePropagatorAlongMomentumToken_, thePropagatorOppositeToMomentum_, and thePropagatorOppositeToMomentumToken_.

Referenced by ConversionTrackCandidateProducer::setEventSetup().

                                                                {
  theGeomSearchTracker_ = es.getHandle(theGeomSearchTrackerToken_);
  theMF_ = es.getHandle(theMFToken_);

  theMeasurementTracker_ = &es.getData(theMeasurementTrackerToken_);

  thePropagatorAlongMomentum_ = &es.getData(thePropagatorAlongMomentumToken_);

  thePropagatorOppositeToMomentum_ = &es.getData(thePropagatorOppositeToMomentumToken_);
}

setMeasurementTracker()

void ConversionSeedFinder::setMeasurementTracker ( const MeasurementTracker *  tracker ) const

inline

Definition at line 68 of file ConversionSeedFinder.h.

{ ; }

setNavigationSchool()

void ConversionSeedFinder::setNavigationSchool ( const NavigationSchool *  navigation )

inline

Definition at line 73 of file ConversionSeedFinder.h.

References theNavigationSchool_.

Referenced by ConversionTrackCandidateProducer::beginRun().

{ theNavigationSchool_ = navigation; }

trackStateFromClusters()

FreeTrajectoryState ConversionSeedFinder::trackStateFromClusters ( int  aCharge, const GlobalPoint &  gpOrigine, PropagationDirection  dir, float  scaleFactor  ) const

protected

Definition at line 69 of file ConversionSeedFinder.cc.

References alongMomentum, ALCARECOTkAlJpsiMuMu_cff::charge, DeadROC_duringRun::dir, visualization-live-secondInstance_cfg::m, PV3DBase< T, PVType, FrameType >::perp(), scaleFactor, mathSSE::sqrt(), theMF_, theSCenergy_, theSCPosition_, PV3DBase< T, PVType, FrameType >::theta(), Vector3DBase< T, FrameTag >::unit(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by findLayers(), and OutInConversionSeedFinder::startSeed().

                                                                                          {
  double caloEnergy = theSCenergy_ * scaleFactor;

  GlobalVector radiusCalo = theSCPosition_ - theOrigin;

  GlobalVector momentumWithoutCurvature = radiusCalo.unit() * caloEnergy;

  GlobalTrajectoryParameters gtp;
  if (dir == alongMomentum) {
    gtp = GlobalTrajectoryParameters(theOrigin, momentumWithoutCurvature, charge, &(*theMF_));
  } else {
    gtp = GlobalTrajectoryParameters(theSCPosition_, momentumWithoutCurvature, charge, &(*theMF_));
  }

  // now create error matrix
  // dpos = 4mm/sqrt(E), dtheta = move vertex by 1sigma
  float dpos = 0.4 / sqrt(theSCenergy_);
  dpos *= 2.;
  float dphi = dpos / theSCPosition_.perp();
  //  float dp = 0.03 * sqrt(theCaloEnergy);
  //  float dp = theCaloEnergy / sqrt(12.); // for fun
  float theta1 = theSCPosition_.theta();
  float theta2 = atan2(double(theSCPosition_.perp()), theSCPosition_.z() - 5.5);
  float dtheta = theta1 - theta2;
  AlgebraicSymMatrix55 m;
  m[0][0] = 1.;
  m[1][1] = dpos * dpos;
  m[2][2] = dpos * dpos;
  m[3][3] = dphi * dphi;
  m[4][4] = dtheta * dtheta;

  FreeTrajectoryState fts(gtp, CurvilinearTrajectoryError(m));

  return fts;
}

Member Data Documentation

beamSpotToken_

edm::EDGetTokenT<reco::BeamSpot> ConversionSeedFinder::beamSpotToken_

protected

Definition at line 106 of file ConversionSeedFinder.h.

Referenced by ConversionSeedFinder(), and setEvent().

dir_

PropagationDirection ConversionSeedFinder::dir_

protected

Definition at line 110 of file ConversionSeedFinder.h.

measurementTrkToken_

edm::EDGetTokenT<MeasurementTrackerEvent> ConversionSeedFinder::measurementTrkToken_

protected

Definition at line 107 of file ConversionSeedFinder.h.

Referenced by ConversionSeedFinder(), and setEvent().

theBCEnergy_

float ConversionSeedFinder::theBCEnergy_

protected

Definition at line 117 of file ConversionSeedFinder.h.

Referenced by OutInConversionSeedFinder::makeSeeds(), and OutInConversionSeedFinder::makeTrackState().

theBCPosition_

GlobalPoint ConversionSeedFinder::theBCPosition_

protected

Definition at line 116 of file ConversionSeedFinder.h.

Referenced by OutInConversionSeedFinder::makeEstimator(), OutInConversionSeedFinder::makeSeeds(), and OutInConversionSeedFinder::makeTrackState().

theBeamSpot_

reco::BeamSpot ConversionSeedFinder::theBeamSpot_

protected

Definition at line 122 of file ConversionSeedFinder.h.

Referenced by findLayers(), OutInConversionSeedFinder::makeTrackState(), and setEvent().

theGeomSearchTracker_

edm::ESHandle<GeometricSearchTracker> ConversionSeedFinder::theGeomSearchTracker_

protected

Definition at line 98 of file ConversionSeedFinder.h.

Referenced by setEventSetup().

theGeomSearchTrackerToken_

edm::ESGetToken<GeometricSearchTracker, TrackerRecoGeometryRecord> ConversionSeedFinder::theGeomSearchTrackerToken_

protected

Definition at line 101 of file ConversionSeedFinder.h.

Referenced by ConversionSeedFinder(), and setEventSetup().

theLayerList_

std::vector<const DetLayer*> ConversionSeedFinder::theLayerList_

protected

Definition at line 114 of file ConversionSeedFinder.h.

Referenced by OutInConversionSeedFinder::completeSeed(), findLayers(), layerList(), and printLayer().

theMeasurementTracker_

const MeasurementTracker* ConversionSeedFinder::theMeasurementTracker_

protected

Definition at line 93 of file ConversionSeedFinder.h.

Referenced by getMeasurementTracker(), and setEventSetup().

theMeasurementTrackerToken_

edm::ESGetToken<MeasurementTracker, CkfComponentsRecord> ConversionSeedFinder::theMeasurementTrackerToken_

protected

Definition at line 102 of file ConversionSeedFinder.h.

Referenced by ConversionSeedFinder(), and setEventSetup().

theMF_

edm::ESHandle<MagneticField> ConversionSeedFinder::theMF_

protected

Definition at line 97 of file ConversionSeedFinder.h.

Referenced by OutInConversionSeedFinder::createSeedFTS(), findLayers(), OutInConversionSeedFinder::makeTrackState(), setEventSetup(), and trackStateFromClusters().

theMFToken_

edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> ConversionSeedFinder::theMFToken_

protected

Definition at line 100 of file ConversionSeedFinder.h.

Referenced by ConversionSeedFinder(), and setEventSetup().

theNavigationSchool_

const NavigationSchool* ConversionSeedFinder::theNavigationSchool_ = nullptr

protected

Definition at line 95 of file ConversionSeedFinder.h.

Referenced by findLayers(), and setNavigationSchool().

thePropagatorAlongMomentum_

const Propagator* ConversionSeedFinder::thePropagatorAlongMomentum_

protected

Definition at line 119 of file ConversionSeedFinder.h.

Referenced by InOutConversionSeedFinder::createSeed(), InOutConversionSeedFinder::findSeeds(), setEventSetup(), and OutInConversionSeedFinder::startSeed().

thePropagatorAlongMomentumToken_

edm::ESGetToken<Propagator, TrackingComponentsRecord> ConversionSeedFinder::thePropagatorAlongMomentumToken_

protected

Definition at line 103 of file ConversionSeedFinder.h.

Referenced by ConversionSeedFinder(), and setEventSetup().

thePropagatorOppositeToMomentum_

const Propagator* ConversionSeedFinder::thePropagatorOppositeToMomentum_

protected

Definition at line 120 of file ConversionSeedFinder.h.

Referenced by OutInConversionSeedFinder::createSeed(), InOutConversionSeedFinder::fillClusterSeeds(), setEventSetup(), and OutInConversionSeedFinder::startSeed().

thePropagatorOppositeToMomentumToken_

edm::ESGetToken<Propagator, TrackingComponentsRecord> ConversionSeedFinder::thePropagatorOppositeToMomentumToken_

protected

Definition at line 104 of file ConversionSeedFinder.h.

Referenced by ConversionSeedFinder(), and setEventSetup().

theSC_

reco::CaloCluster* ConversionSeedFinder::theSC_

protected

Definition at line 111 of file ConversionSeedFinder.h.

theSCenergy_

float ConversionSeedFinder::theSCenergy_

protected

Definition at line 112 of file ConversionSeedFinder.h.

Referenced by OutInConversionSeedFinder::createSeedFTS(), setCandidate(), and trackStateFromClusters().

theSCPosition_

GlobalPoint ConversionSeedFinder::theSCPosition_

protected

Definition at line 91 of file ConversionSeedFinder.h.

Referenced by OutInConversionSeedFinder::createSeedFTS(), OutInConversionSeedFinder::fixPointRadius(), OutInConversionSeedFinder::makeSeeds(), setCandidate(), and trackStateFromClusters().

theSeeds_

TrajectorySeedCollection ConversionSeedFinder::theSeeds_

protected

Definition at line 90 of file ConversionSeedFinder.h.

Referenced by clear(), OutInConversionSeedFinder::createSeed(), InOutConversionSeedFinder::createSeed(), OutInConversionSeedFinder::makeSeeds(), InOutConversionSeedFinder::makeSeeds(), and seeds().

theTrackerData_

edm::Handle<MeasurementTrackerEvent> ConversionSeedFinder::theTrackerData_

protected

Definition at line 123 of file ConversionSeedFinder.h.

Referenced by OutInConversionSeedFinder::completeSeed(), InOutConversionSeedFinder::completeSeed(), InOutConversionSeedFinder::findSeeds(), setEvent(), and OutInConversionSeedFinder::startSeed().

theTrackerGeom_

const TrackingGeometry* ConversionSeedFinder::theTrackerGeom_

protected

Definition at line 94 of file ConversionSeedFinder.h.

Referenced by InOutConversionSeedFinder::fillClusterSeeds(), and setEvent().

theUpdator_

KFUpdator ConversionSeedFinder::theUpdator_

protected

Definition at line 109 of file ConversionSeedFinder.h.

Referenced by OutInConversionSeedFinder::createSeed(), and InOutConversionSeedFinder::createSeed().