#include <SeedForPhotonConversion1Leg.h>

Public Member Functions
	SeedForPhotonConversion1Leg (const edm::ParameterSet &cfg)

const TrajectorySeed *	trajectorySeed (TrajectorySeedCollection &seedCollection, const SeedingHitSet &hits, const GlobalPoint &vertex, const GlobalVector &vertexBounds, float ptmin, const edm::EventSetup &es, float cotTheta, std::stringstream &ss)

	~SeedForPhotonConversion1Leg ()

Static Public Attributes
static const int	cotTheta_Max =99999

Protected Member Functions
const TrajectorySeed *	buildSeed (TrajectorySeedCollection &seedCollection, const SeedingHitSet &hits, const FreeTrajectoryState &fts, const edm::EventSetup &es) const

bool	checkHit (const TrajectoryStateOnSurface &, const SeedingHitSet::ConstRecHitPointer &hit, const edm::EventSetup &es) const

CurvilinearTrajectoryError	initialError (const GlobalVector &vertexBounds, float ptMin, float sinTheta) const

GlobalTrajectoryParameters	initialKinematic (const SeedingHitSet &hits, const GlobalPoint &vertexPos, const edm::EventSetup &es, const float cotTheta) const

SeedingHitSet::RecHitPointer	refitHit (SeedingHitSet::ConstRecHitPointer hit, const TrajectoryStateOnSurface &state, const TkClonerImpl &cloner) const

Protected Attributes
PrintRecoObjects	po

std::stringstream *	pss

double	theBOFFMomentum

std::string	thePropagatorLabel

std::string	TTRHBuilder

Detailed Description

Definition at line 22 of file SeedForPhotonConversion1Leg.h.

Constructor & Destructor Documentation

SeedForPhotonConversion1Leg::SeedForPhotonConversion1Leg ( const edm::ParameterSet & cfg )

inline

Definition at line 26 of file SeedForPhotonConversion1Leg.h.

                                                            :
     thePropagatorLabel (cfg.getParameter<std::string>("propagator"))
     ,theBOFFMomentum   (cfg.getParameter<double>("SeedMomentumForBOFF"))
     ,TTRHBuilder       (cfg.getParameter<std::string>("TTRHBuilder"))
       {}

SeedForPhotonConversion1Leg::~SeedForPhotonConversion1Leg ( )

inline

Definition at line 33 of file SeedForPhotonConversion1Leg.h.

References hfClusterShapes_cfi::hits, and ptmin.

33 {}

Member Function Documentation

const TrajectorySeed * SeedForPhotonConversion1Leg::buildSeed	(	TrajectorySeedCollection &	seedCollection,
		const SeedingHitSet &	hits,
		const FreeTrajectoryState &	fts,
		const edm::EventSetup &	es
	)		const

protected

Definition at line 134 of file SeedForPhotonConversion1Leg.cc.

References alongMomentum, checkHit(), TrackingRecHit::geographicalId(), edm::EventSetup::get(), TrajectoryStateOnSurface::globalMomentum(), TrackerGeometry::idToDet(), TrajectoryStateOnSurface::isValid(), TrackingRecHit::localPosition(), PV3DBase< T, PVType, FrameType >::perp(), trajectoryStateTransform::persistentState(), po, PrintRecoObjects::print(), edm::ESHandle< T >::product(), Propagator::propagate(), PhotonConversionTrajectorySeedProducerFromQuadruplets_cfi::propagator, edm::OwnVector< T, P >::push_back(), DetId::rawId(), refitHit(), SeedingHitSet::size(), thePropagatorLabel, trackingTruthProducer_cfi::tracker, TTRHBuilder, KFUpdator::update(), TrackInfoProducer_cfi::updatedState, and gsfElectronCkfTrackCandidateMaker_cff::updator.

Referenced by trajectorySeed().

 {
   // FIXME all this stuff shoould go in an initialized...
 
   // get tracker
   edm::ESHandle<TrackerGeometry> tracker;
   es.get<TrackerDigiGeometryRecord>().get(tracker);
   
   // get propagator
   edm::ESHandle<Propagator>  propagatorHandle;
   es.get<TrackingComponentsRecord>().get(thePropagatorLabel, propagatorHandle);
   const Propagator*  propagator = &(*propagatorHandle);
   
   edm::ESHandle<TransientTrackingRecHitBuilder> builderH;
   es.get<TransientRecHitRecord>().get(TTRHBuilder, builderH);
   auto builder = (TkTransientTrackingRecHitBuilder const *)(builderH.product());
   auto cloner = (*builder).cloner();
 
   // get updator
   KFUpdator  updator;
   
 
 
   // Now update initial state track using information from seed hits.
   
   TrajectoryStateOnSurface updatedState;
   edm::OwnVector<TrackingRecHit> seedHits;
   
   const TrackingRecHit* hit = nullptr;
   for ( unsigned int iHit = 0; iHit < hits.size() && iHit<1; iHit++) {
     hit = hits[iHit];
     TrajectoryStateOnSurface state = (iHit==0) ? 
       propagator->propagate(fts,tracker->idToDet(hit->geographicalId())->surface())
       : propagator->propagate(updatedState, tracker->idToDet(hit->geographicalId())->surface());
     if (!state.isValid()) return nullptr;
     
     SeedingHitSet::ConstRecHitPointer tth = hits[iHit]; 
     
     std::unique_ptr<BaseTrackerRecHit> newtth(refitHit( tth, state, cloner));
 
     
     if (!checkHit(state,&*newtth,es)) return nullptr;
 
     updatedState =  updator.update(state, *newtth);
     if (!updatedState.isValid()) return nullptr;
     
     seedHits.push_back(newtth.release());
 #ifdef mydebug_seed
     uint32_t detid = hit->geographicalId().rawId();
     (*pss) << "\n[SeedForPhotonConversion1Leg] hit " << iHit;
     po.print(*pss, detid);
     (*pss) << " "  << detid << "\t lp " << hit->localPosition()
        << " tth " << tth->localPosition() << " newtth " << newtth->localPosition() << " state " << state.globalMomentum().perp();
 #endif
   } 
   
   if(!hit) return nullptr;
 
   PTrajectoryStateOnDet const & PTraj =
       trajectoryStateTransform::persistentState(updatedState, hit->geographicalId().rawId());
   
   seedCollection.push_back( TrajectorySeed(PTraj,seedHits,alongMomentum));
   return &seedCollection.back();
 }

bool SeedForPhotonConversion1Leg::checkHit	(	const TrajectoryStateOnSurface &	,
		const SeedingHitSet::ConstRecHitPointer &	hit,
		const edm::EventSetup &	es
	)		const

inlineprotected

Definition at line 47 of file SeedForPhotonConversion1Leg.h.

References ALCARECOTkAlBeamHalo_cff::ptMin.

Referenced by buildSeed().

50 { return true; }

CurvilinearTrajectoryError SeedForPhotonConversion1Leg::initialError	(	const GlobalVector &	vertexBounds,
		float	ptMin,
		float	sinTheta
	)		const

protected

Definition at line 105 of file SeedForPhotonConversion1Leg.cc.

References patCaloMETCorrections_cff::C, GlobalErrorBase< T, ErrorWeightType >::cxx(), GlobalErrorBase< T, ErrorWeightType >::czz(), SiStripPI::max, funct::sqr(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by initialKinematic(), and trajectorySeed().

 {
   // Set initial uncertainty on track parameters, using only P.V. constraint and no hit
   // information.
   GlobalError vertexErr( sqr(vertexBounds.x()), 0, 
              sqr(vertexBounds.y()), 0, 0,
              sqr(vertexBounds.z())
              );
   
  
   AlgebraicSymMatrix55 C = ROOT::Math::SMatrixIdentity();
 
 // FIXME: minC00. Prevent apriori uncertainty in 1/P from being too small, 
 // to avoid instabilities.
 // N.B. This parameter needs optimising ...
   float sin2th = sqr(sinTheta);
   float minC00 = 1.0;
   C[0][0] = std::max(sin2th/sqr(ptMin), minC00);
   float zErr = vertexErr.czz();
   float transverseErr = vertexErr.cxx(); // assume equal cxx cyy
   C[3][3] = transverseErr;
   C[4][4] = zErr*sin2th + transverseErr*(1-sin2th);
 
   return CurvilinearTrajectoryError(C);
 }

GlobalTrajectoryParameters SeedForPhotonConversion1Leg::initialKinematic	(	const SeedingHitSet &	hits,
		const GlobalPoint &	vertexPos,
		const edm::EventSetup &	es,
		const float	cotTheta
	)		const

protected

Definition at line 46 of file SeedForPhotonConversion1Leg.cc.

References GlobalTrajectoryParameters::charge(), cotTheta_Max, edm::EventSetup::get(), initialError(), GlobalTrajectoryParameters::magneticField(), GlobalTrajectoryParameters::momentum(), MagneticField::nominalValue(), PV3DBase< T, PVType, FrameType >::perp(), po, GlobalTrajectoryParameters::position(), PrintRecoObjects::print(), theBOFFMomentum, GlobalTrajectoryParameters::transverseCurvature(), Vector3DBase< T, FrameTag >::unit(), PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().

Referenced by trajectorySeed().

 {
   GlobalTrajectoryParameters kine;
 
   SeedingHitSet::ConstRecHitPointer tth1 = hits[0];
   SeedingHitSet::ConstRecHitPointer tth2 = hits[1];
 
    // FIXME optimize: move outside loop
     edm::ESHandle<MagneticField> bfield;
     es.get<IdealMagneticFieldRecord>().get(bfield);
     float nomField = bfield->nominalValue();
 
   FastHelix helix(tth2->globalPosition(), tth1->globalPosition(), vertexPos, nomField, &*bfield, vertexPos);
   kine = helix.stateAtVertex();
 
   //force the pz/pt equal to the measured one
   if(fabs(cotTheta)<cotTheta_Max)
     kine = GlobalTrajectoryParameters(kine.position(),
                       GlobalVector(kine.momentum().x(),kine.momentum().y(),kine.momentum().perp()*cotTheta),
                       kine.charge(),
                       & kine.magneticField()
                       );
   else
     kine = GlobalTrajectoryParameters(kine.position(),
                       GlobalVector(kine.momentum().x(),kine.momentum().y(),kine.momentum().perp()*cotTheta_Max),
                       kine.charge(),
                       & kine.magneticField()
                       );
 
 #ifdef mydebug_seed
   uint32_t detid;
   (*pss) << "[SeedForPhotonConversion1Leg] initialKinematic tth1 " ;
   detid=tth1->geographicalId().rawId();
   po.print(*pss, detid );
   (*pss) << " \t " << detid << " " << tth1->localPosition()  << " " << tth1->globalPosition()    ;
   detid= tth2->geographicalId().rawId();
   (*pss) << " \n\t tth2 ";
   po.print(*pss, detid );
   (*pss) << " \t " << detid << " " << tth2->localPosition()  << " " << tth2->globalPosition()  
      << "\nhelix momentum " << kine.momentum() << " pt " << kine.momentum().perp() << " radius " << 1/kine.transverseCurvature(); 
 #endif
 
   bool isBOFF =(0==nomField);;
   if (isBOFF && (theBOFFMomentum > 0)) {
     kine = GlobalTrajectoryParameters(kine.position(),
                               kine.momentum().unit() * theBOFFMomentum,
                               kine.charge(),
                               &*bfield);
   }
   return kine;
 }

SeedingHitSet::RecHitPointer SeedForPhotonConversion1Leg::refitHit	(	SeedingHitSet::ConstRecHitPointer	hit,
		const TrajectoryStateOnSurface &	state,
		const TkClonerImpl &	cloner
	)		const

protected

Definition at line 203 of file SeedForPhotonConversion1Leg.cc.

Referenced by buildSeed().

 {
   //const TransientTrackingRecHit* a= hit.get();
   //return const_cast<TransientTrackingRecHit*> (a);
   //This was modified otherwise the rechit will have just the local x component and local y=0
   // To understand how to modify for pixels
 
   //const TSiStripRecHit2DLocalPos* b = dynamic_cast<const TSiStripRecHit2DLocalPos*>(a);
   //return const_cast<TSiStripRecHit2DLocalPos*>(b);
   return (SeedingHitSet::RecHitPointer)(cloner(*hit,state));
 }

const TrajectorySeed * SeedForPhotonConversion1Leg::trajectorySeed	(	TrajectorySeedCollection &	seedCollection,
		const SeedingHitSet &	hits,
		const GlobalPoint &	vertex,
		const GlobalVector &	vertexBounds,
		float	ptmin,
		const edm::EventSetup &	es,
		float	cotTheta,
		std::stringstream &	ss
	)

Definition at line 24 of file SeedForPhotonConversion1Leg.cc.

References buildSeed(), relativeConstraints::error, initialError(), initialKinematic(), GlobalTrajectoryParameters::momentum(), funct::sin(), SeedingHitSet::size(), and PV3DBase< T, PVType, FrameType >::theta().

 {
   pss = &ss;
   if ( hits.size() < 2) return nullptr;
 
   GlobalTrajectoryParameters kine = initialKinematic(hits, vertex, es, cotTheta);
   float sinTheta = sin(kine.momentum().theta());
 
   CurvilinearTrajectoryError error = initialError(vertexBounds, ptmin,  sinTheta);
   FreeTrajectoryState fts(kine, error);
 
   return buildSeed(seedCollection,hits,fts,es); 
 }