#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_

std::vector< const DetLayer * >	theLayerList_

const MeasurementTracker *	theMeasurementTracker_

std::string	theMeasurementTrackerName_

edm::ESHandle< MagneticField >	theMF_

const NavigationSchool *	theNavigationSchool_ = nullptr

const Propagator *	thePropagatorAlongMomentum_

const Propagator *	thePropagatorOppositeToMomentum_

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 47 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.

     :  //  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)
  
   LogDebug("ConversionSeedFinder") << " CTOR "
                                    << "\n";
  
   theMeasurementTrackerName_ = config.getParameter<std::string>("MeasurementTrackerName");
 }

References beamSpotToken_, edm::ConsumesCollector::consumes(), HLT_2018_cff::InputTag, LogDebug, measurementTrkToken_, AlCaHLTBitMon_QueryRunRegistry::string, and theMeasurementTrackerName_.

◆ ~ConversionSeedFinder()

virtual ConversionSeedFinder::~ConversionSeedFinder ( )

inlinevirtual

Definition at line 52 of file ConversionSeedFinder.h.

52 {}

Member Function Documentation

◆ clear()

void ConversionSeedFinder::clear ( void )

inline

Definition at line 71 of file ConversionSeedFinder.h.

71 { theSeeds_.clear(); }

References theSeeds_.

◆ findLayers() [1/2]

void ConversionSeedFinder::findLayers ( )

protected

Definition at line 58 of file ConversionSeedFinder.cc.

                                       {
   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);
 }

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

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

◆ findLayers() [2/2]

void ConversionSeedFinder::findLayers ( const FreeTrajectoryState & fts )

protected

Definition at line 109 of file ConversionSeedFinder.cc.

                                                                      {
   theLayerList_.clear();
  
   StraightLinePropagator prop(&(*theMF_), alongMomentum);
  
   StartingLayerFinder starter(&prop, this->getMeasurementTracker());
  
   LayerCollector collector(theNavigationSchool_, &prop, &starter, 5., 5.);
  
   theLayerList_ = collector.allLayers(traj);
  
   for (unsigned int i = 0; i < theLayerList_.size(); ++i) {
     printLayer(i);
   }
 }

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

◆ getMeasurementTracker()

const MeasurementTracker* ConversionSeedFinder::getMeasurementTracker ( ) const

inline

Definition at line 64 of file ConversionSeedFinder.h.

64 { return theMeasurementTracker_; }

References theMeasurementTracker_.

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

◆ layerList()

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

inline

Definition at line 61 of file ConversionSeedFinder.h.

61 { return theLayerList_; }

References theLayerList_.

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

◆ 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 125 of file ConversionSeedFinder.cc.

                                                  {
   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";
   }
 }

References GeomDetEnumerators::barrel, mps_fire::i, DetLayer::location(), and theLayerList_.

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

◆ seeds()

TrajectorySeedCollection& ConversionSeedFinder::seeds ( )

inline

Definition at line 56 of file ConversionSeedFinder.h.

56 { return theSeeds_; }

References theSeeds_.

◆ setCandidate()

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

inlinevirtual

Definition at line 57 of file ConversionSeedFinder.h.

                                                       {
     theSCenergy_ = e;
     theSCPosition_ = pos;
   }

References MillePedeFileConverter_cfg::e, theSCenergy_, and theSCPosition_.

◆ setEvent()

void ConversionSeedFinder::setEvent ( const edm::Event & e )

Definition at line 30 of file ConversionSeedFinder.cc.

                                                        {
   theTrackerGeom_ = this->getMeasurementTracker()->geomTracker();
  
   //get the BeamSpot
   edm::Handle<reco::BeamSpot> recoBeamSpotHandle;
   evt.getByToken(beamSpotToken_, recoBeamSpotHandle);
   theBeamSpot_ = *recoBeamSpotHandle;
  
   evt.getByToken(measurementTrkToken_, theTrackerData_);
 }

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

◆ setEventSetup()

void ConversionSeedFinder::setEventSetup ( const edm::EventSetup & es )

Initialize EventSetup objects at each event.

Definition at line 41 of file ConversionSeedFinder.cc.

                                                                 {
   es.get<TrackerRecoGeometryRecord>().get(theGeomSearchTracker_);
   es.get<IdealMagneticFieldRecord>().get(theMF_);
  
   edm::ESHandle<MeasurementTracker> measurementTrackerHandle;
   es.get<CkfComponentsRecord>().get(theMeasurementTrackerName_, measurementTrackerHandle);
   theMeasurementTracker_ = measurementTrackerHandle.product();
  
   edm::ESHandle<Propagator> propagatorAlongMomHandle;
   es.get<TrackingComponentsRecord>().get("alongMomElePropagator", propagatorAlongMomHandle);
   thePropagatorAlongMomentum_ = &(*propagatorAlongMomHandle);
  
   edm::ESHandle<Propagator> propagatorOppoToMomHandle;
   es.get<TrackingComponentsRecord>().get("oppositeToMomElePropagator", propagatorOppoToMomHandle);
   thePropagatorOppositeToMomentum_ = &(*propagatorOppoToMomHandle);
 }

References edm::EventSetup::get(), get, edm::ESHandle< T >::product(), theGeomSearchTracker_, theMeasurementTracker_, theMeasurementTrackerName_, theMF_, thePropagatorAlongMomentum_, and thePropagatorOppositeToMomentum_.

◆ setMeasurementTracker()

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

inline

Definition at line 63 of file ConversionSeedFinder.h.

63 { ; }

◆ setNavigationSchool()

void ConversionSeedFinder::setNavigationSchool ( const NavigationSchool * navigation )

inline

Definition at line 68 of file ConversionSeedFinder.h.

68 { theNavigationSchool_ = navigation; }

References theNavigationSchool_.

◆ trackStateFromClusters()

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

protected

Definition at line 70 of file ConversionSeedFinder.cc.

                                                                                           {
   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;
 }

References alongMomentum, ALCARECOTkAlJpsiMuMu_cff::charge, DeadROC_duringRun::dir, visualization-live-secondInstance_cfg::m, PV3DBase< T, PVType, FrameType >::perp(), multiplicitycorr_cfi::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().

Member Data Documentation

◆ beamSpotToken_

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

protected

Definition at line 96 of file ConversionSeedFinder.h.

Referenced by ConversionSeedFinder(), and setEvent().

◆ dir_

PropagationDirection ConversionSeedFinder::dir_

protected

Definition at line 100 of file ConversionSeedFinder.h.

◆ measurementTrkToken_

edm::EDGetTokenT<MeasurementTrackerEvent> ConversionSeedFinder::measurementTrkToken_

protected

Definition at line 97 of file ConversionSeedFinder.h.

Referenced by ConversionSeedFinder(), and setEvent().

◆ theBCEnergy_

float ConversionSeedFinder::theBCEnergy_

protected

Definition at line 107 of file ConversionSeedFinder.h.

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

◆ theBCPosition_

GlobalPoint ConversionSeedFinder::theBCPosition_

protected

Definition at line 106 of file ConversionSeedFinder.h.

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

◆ theBeamSpot_

reco::BeamSpot ConversionSeedFinder::theBeamSpot_

protected

Definition at line 112 of file ConversionSeedFinder.h.

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

◆ theGeomSearchTracker_

edm::ESHandle<GeometricSearchTracker> ConversionSeedFinder::theGeomSearchTracker_

protected

Definition at line 94 of file ConversionSeedFinder.h.

Referenced by setEventSetup().

◆ theLayerList_

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

protected

Definition at line 104 of file ConversionSeedFinder.h.

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

◆ theMeasurementTracker_

const MeasurementTracker* ConversionSeedFinder::theMeasurementTracker_

protected

Definition at line 89 of file ConversionSeedFinder.h.

Referenced by getMeasurementTracker(), and setEventSetup().

◆ theMeasurementTrackerName_

std::string ConversionSeedFinder::theMeasurementTrackerName_

protected

Definition at line 88 of file ConversionSeedFinder.h.

Referenced by ConversionSeedFinder(), and setEventSetup().

◆ theMF_

edm::ESHandle<MagneticField> ConversionSeedFinder::theMF_

protected

Definition at line 93 of file ConversionSeedFinder.h.

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

◆ theNavigationSchool_

const NavigationSchool* ConversionSeedFinder::theNavigationSchool_ = nullptr

protected

Definition at line 91 of file ConversionSeedFinder.h.

Referenced by findLayers(), and setNavigationSchool().

◆ thePropagatorAlongMomentum_

const Propagator* ConversionSeedFinder::thePropagatorAlongMomentum_

protected

Definition at line 109 of file ConversionSeedFinder.h.

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

◆ thePropagatorOppositeToMomentum_

const Propagator* ConversionSeedFinder::thePropagatorOppositeToMomentum_

protected

Definition at line 110 of file ConversionSeedFinder.h.

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

◆ theSC_

reco::CaloCluster* ConversionSeedFinder::theSC_

protected

Definition at line 101 of file ConversionSeedFinder.h.

◆ theSCenergy_

float ConversionSeedFinder::theSCenergy_

protected

Definition at line 102 of file ConversionSeedFinder.h.

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

◆ theSCPosition_

GlobalPoint ConversionSeedFinder::theSCPosition_

protected

Definition at line 86 of file ConversionSeedFinder.h.

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

◆ theSeeds_

TrajectorySeedCollection ConversionSeedFinder::theSeeds_

protected

Definition at line 85 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 113 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 90 of file ConversionSeedFinder.h.

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

◆ theUpdator_

KFUpdator ConversionSeedFinder::theUpdator_

protected

Definition at line 99 of file ConversionSeedFinder.h.

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

Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ ConversionSeedFinder() [1/2]

◆ ConversionSeedFinder() [2/2]

◆ ~ConversionSeedFinder()

Member Function Documentation

◆ clear()

◆ findLayers() [1/2]

◆ findLayers() [2/2]

◆ getMeasurementTracker()

◆ layerList()

◆ makeSeeds()

◆ printLayer()

◆ seeds()

◆ setCandidate()

◆ setEvent()

◆ setEventSetup()

◆ setMeasurementTracker()

◆ setNavigationSchool()

◆ trackStateFromClusters()

Member Data Documentation

◆ beamSpotToken_

◆ dir_

◆ measurementTrkToken_

◆ theBCEnergy_

◆ theBCPosition_

◆ theBeamSpot_

◆ theGeomSearchTracker_

◆ theLayerList_

◆ theMeasurementTracker_

◆ theMeasurementTrackerName_

◆ theMF_

◆ theNavigationSchool_

◆ thePropagatorAlongMomentum_

◆ thePropagatorOppositeToMomentum_

◆ theSC_

◆ theSCenergy_

◆ theSCPosition_

◆ theSeeds_

◆ theTrackerData_

◆ theTrackerGeom_

◆ theUpdator_