#include <ConversionSeedFinder.h>

Inheritance diagram for ConversionSeedFinder:

Public Member Functions
void	clear ()

	ConversionSeedFinder ()

	ConversionSeedFinder (const edm::ParameterSet &config)

const MeasurementTracker *	getMeasurementTracker () const

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

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

TrajectorySeedCollection &	seeds ()

virtual void	setCandidate (float e, GlobalPoint pos) const

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

virtual	~ConversionSeedFinder ()

Protected Member Functions
void	findLayers () const

void	findLayers (const FreeTrajectoryState &fts) const

void	printLayer (int i) const

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

Protected Attributes
edm::ParameterSet	conf_

PropagationDirection	dir_

float	theBCEnergy_

GlobalPoint	theBCPosition_

reco::BeamSpot	theBeamSpot_

edm::ESHandle < GeometricSearchTracker >	theGeomSearchTracker_

std::vector< const DetLayer * >	theLayerList_

const MeasurementTracker *	theMeasurementTracker_

edm::ESHandle< MagneticField >	theMF_

const Propagator *	thePropagatorAlongMomentum_

const Propagator *	thePropagatorOppositeToMomentum_

reco::CaloCluster *	theSC_

float	theSCenergy_

GlobalPoint	theSCPosition_

TrajectorySeedCollection	theSeeds_

const TrackingGeometry *	theTrackerGeom_

KFUpdator	theUpdator_

Detailed Description

Id:: ConversionSeedFinder.h,v 1.11 2008/11/19 16:47:44 nancy Exp

Date:: 2008/11/19 16:47:44

Revision:: 1.11

Author: Nancy Marinelli, U. of Notre Dame, US

Definition at line 50 of file ConversionSeedFinder.h.

Constructor & Destructor Documentation

ConversionSeedFinder::ConversionSeedFinder ( )

ConversionSeedFinder::ConversionSeedFinder ( const edm::ParameterSet & config )

Definition at line 18 of file ConversionSeedFinder.cc.

References LogDebug.

                                                                         : 
   conf_(config),
   theUpdator_()
 {
 
   LogDebug("ConversionSeedFinder")  << " CTOR " << "\n";
 
 
 }

virtual ConversionSeedFinder::~ConversionSeedFinder ( )

inlinevirtual

Definition at line 58 of file ConversionSeedFinder.h.

58 {}

Member Function Documentation

void ConversionSeedFinder::clear ( void )

inline

Definition at line 77 of file ConversionSeedFinder.h.

References theSeeds_.

Referenced by python.Vispa.Views.BoxDecayView.BoxDecayView::closeEvent(), ConversionTrackCandidateProducer::produce(), python.Vispa.Views.LineDecayView.LineDecayView::setDataObjects(), and python.Vispa.Views.BoxDecayView.BoxDecayView::updateContent().

                {
     theSeeds_.clear();
   }

void ConversionSeedFinder::findLayers ( ) const

protected

Definition at line 68 of file ConversionSeedFinder.cc.

References alongMomentum, DeDxDiscriminatorTools::charge(), reco::BeamSpot::position(), theBeamSpot_, and trackStateFromClusters().

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

void ConversionSeedFinder::findLayers ( const FreeTrajectoryState & fts ) const

protected

Definition at line 127 of file ConversionSeedFinder.cc.

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

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

const MeasurementTracker* ConversionSeedFinder::getMeasurementTracker ( ) const

inline

Definition at line 71 of file ConversionSeedFinder.h.

References theMeasurementTracker_.

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

71 {return theMeasurementTracker_;}

ConversionSeedFinder::theMeasurementTracker_

const MeasurementTracker * theMeasurementTracker_

Definition: ConversionSeedFinder.h:101

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

inline

Definition at line 67 of file ConversionSeedFinder.h.

References theLayerList_.

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

67 { return theLayerList_;}

ConversionSeedFinder::theLayerList_

std::vector< const DetLayer * > theLayerList_

Definition: ConversionSeedFinder.h:115

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

pure virtual

Implemented in InOutConversionSeedFinder, and OutInConversionSeedFinder.

void ConversionSeedFinder::printLayer ( int i ) const

protected

Definition at line 151 of file ConversionSeedFinder.cc.

References GeomDetEnumerators::barrel, i, DetLayer::location(), 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";
   }
 }

TrajectorySeedCollection& ConversionSeedFinder::seeds ( )

inline

Definition at line 65 of file ConversionSeedFinder.h.

References theSeeds_.

Referenced by ConversionTrackCandidateProducer::buildCollections().

65 { return theSeeds_;}

ConversionSeedFinder::theSeeds_

TrajectorySeedCollection theSeeds_

Definition: ConversionSeedFinder.h:97

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

inlinevirtual

Definition at line 66 of file ConversionSeedFinder.h.

References pos, theSCenergy_, and theSCPosition_.

Referenced by ConversionTrackCandidateProducer::buildCollections().

66 { theSCenergy_=e; theSCPosition_= pos; }

ConversionSeedFinder::theSCPosition_

GlobalPoint theSCPosition_

Definition: ConversionSeedFinder.h:98

pos

Definition: Histograms.cc:19

ConversionSeedFinder::theSCenergy_

float theSCenergy_

Definition: ConversionSeedFinder.h:112

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

Definition at line 30 of file ConversionSeedFinder.cc.

References MeasurementTracker::geomTracker(), edm::Event::getByLabel(), getMeasurementTracker(), theBeamSpot_, theMeasurementTracker_, theTrackerGeom_, and MeasurementTracker::update().

Referenced by ConversionTrackCandidateProducer::produce().

                                                           {
  
   theMeasurementTracker_->update(evt);
   theTrackerGeom_= this->getMeasurementTracker()->geomTracker();
 
   //get the BeamSpot
   edm::Handle<reco::BeamSpot> recoBeamSpotHandle;
   evt.getByLabel("offlineBeamSpot",recoBeamSpotHandle);
   theBeamSpot_ = *recoBeamSpotHandle;
 
 
 
 }

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

Initialize EventSetup objects at each event.

Definition at line 44 of file ConversionSeedFinder.cc.

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

Referenced by ConversionTrackCandidateProducer::setEventSetup().

                                                                    {
   es.get<TrackerRecoGeometryRecord>().get( theGeomSearchTracker_);
   es.get<IdealMagneticFieldRecord>().get( theMF_ );
 
 
   edm::ESHandle<MeasurementTracker> measurementTrackerHandle;
   es.get<CkfComponentsRecord>().get(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);
 
 
 
 
 }

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

inline

Definition at line 70 of file ConversionSeedFinder.h.

70 { ; }

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

protected

Definition at line 84 of file ConversionSeedFinder.cc.

References alongMomentum, m, PV3DBase< T, PVType, FrameType >::perp(), 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;
   AlgebraicSymMatrix  m(5,1) ;
   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 ;
 
 
 }