#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 20 of file SeedForPhotonConversion1Leg.h.

Constructor & Destructor Documentation

◆ SeedForPhotonConversion1Leg()

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

inline

Definition at line 24 of file SeedForPhotonConversion1Leg.h.

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

◆ ~SeedForPhotonConversion1Leg()

SeedForPhotonConversion1Leg::~SeedForPhotonConversion1Leg ( )

inline

Definition at line 30 of file SeedForPhotonConversion1Leg.h.

30 {}

Member Function Documentation

◆ buildSeed()

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

protected

Definition at line 124 of file SeedForPhotonConversion1Leg.cc.

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

References alongMomentum, checkHit(), edm::EventSetup::get(), get, TrajectoryStateOnSurface::globalMomentum(), hfClusterShapes_cfi::hits, TrajectoryStateOnSurface::isValid(), PV3DBase< T, PVType, FrameType >::perp(), trajectoryStateTransform::persistentState(), po, PrintRecoObjects::print(), edm::ESHandle< T >::product(), TrackCandidateProducer_cfi::propagator, edm::OwnVector< T, P >::push_back(), refitHit(), ElectronSeedTrackRefFix_cfi::seedCollection, thePropagatorLabel, PbPb_ZMuSkimMuonDPG_cff::tracker, TTRHBuilder, TrackInfoProducer_cfi::updatedState, and HLT_2018_cff::updator.

Referenced by trajectorySeed().

◆ checkHit()

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

inlineprotected

Definition at line 42 of file SeedForPhotonConversion1Leg.h.

                                                {
     return true;
   }

Referenced by buildSeed().

◆ initialError()

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

protected

Definition at line 101 of file SeedForPhotonConversion1Leg.cc.

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

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

Referenced by trajectorySeed().

◆ initialKinematic()

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

protected

Definition at line 48 of file SeedForPhotonConversion1Leg.cc.

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

References GlobalTrajectoryParameters::charge(), cotTheta_Max, edm::EventSetup::get(), get, hfClusterShapes_cfi::hits, 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().

◆ refitHit()

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

protected

Definition at line 192 of file SeedForPhotonConversion1Leg.cc.

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

Referenced by buildSeed().

◆ trajectorySeed()

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 27 of file SeedForPhotonConversion1Leg.cc.

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