CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
List of all members | Public Member Functions | Static Public Attributes | Protected Member Functions | Protected Attributes
SeedForPhotonConversion1Leg Class Reference

#include <SeedForPhotonConversion1Leg.h>

Public Member Functions

 SeedForPhotonConversion1Leg (const edm::ParameterSet &cfg)
 
 SeedForPhotonConversion1Leg (const std::string &propagator="PropagatorWithMaterial", double seedMomentumForBOFF=-5.0)
 
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
 

Protected Attributes

TkClonerImpl cloner
 
PrintRecoObjects po
 
std::stringstream * pss
 
double theBOFFMomentum
 
std::string thePropagatorLabel
 

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.

26  :
27  thePropagatorLabel(cfg.getParameter<std::string>("propagator")),
28  theBOFFMomentum(cfg.existsAs<double>("SeedMomentumForBOFF") ? cfg.getParameter<double>("SeedMomentumForBOFF") : 5.0)
29  {}
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
edm::ParameterSet::existsAs
bool existsAs(std::string const &parameterName, bool trackiness=true) const
checks if a parameter exists as a given type
Definition: ParameterSet.h:184
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
SeedForPhotonConversion1Leg::SeedForPhotonConversion1Leg ( const std::string &  propagator = "PropagatorWithMaterial",
double  seedMomentumForBOFF = -5.0 
)
inline

Definition at line 31 of file SeedForPhotonConversion1Leg.h.

33  : thePropagatorLabel(propagator), theBOFFMomentum(seedMomentumForBOFF) { }
SeedForPhotonConversion1Leg::~SeedForPhotonConversion1Leg ( )
inline

Definition at line 36 of file SeedForPhotonConversion1Leg.h.

36 {}

Member Function Documentation

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

Definition at line 132 of file SeedForPhotonConversion1Leg.cc.

References alongMomentum, checkHit(), cloner, cond::rpcobgas::detid, TrackingRecHit::geographicalId(), edm::EventSetup::get(), TrajectoryStateOnSurface::globalMomentum(), TrajectoryStateOnSurface::isValid(), TrackingRecHit::localPosition(), PV3DBase< T, PVType, FrameType >::perp(), trajectoryStateTransform::persistentState(), po, PrintRecoObjects::print(), edm::ESHandle< class >::product(), Propagator::propagate(), edm::OwnVector< T, P >::push_back(), DetId::rawId(), refitHit(), SeedingHitSet::size(), thePropagatorLabel, patCandidatesForDimuonsSequences_cff::tracker, and MixedLayerPairs_cfi::TTRHBuilder.

Referenced by trajectorySeed().

137 {
138  // get tracker
139  edm::ESHandle<TrackerGeometry> tracker;
140  es.get<TrackerDigiGeometryRecord>().get(tracker);
141 
142  // get propagator
143  edm::ESHandle<Propagator> propagatorHandle;
144  es.get<TrackingComponentsRecord>().get(thePropagatorLabel, propagatorHandle);
145  const Propagator* propagator = &(*propagatorHandle);
146 
147  // get cloner (FIXME: add to config)
148  try {
149  auto TTRHBuilder = "WithTrackAngle";
150  edm::ESHandle<TransientTrackingRecHitBuilder> builderH;
151  es.get<TransientRecHitRecord>().get(TTRHBuilder, builderH);
152  auto builder = (TkTransientTrackingRecHitBuilder const *)(builderH.product());
153  cloner = (*builder).cloner();
154  } catch(...) {
155  auto TTRHBuilder = "hltESPTTRHBWithTrackAngle";
156  edm::ESHandle<TransientTrackingRecHitBuilder> builderH;
157  es.get<TransientRecHitRecord>().get(TTRHBuilder, builderH);
158  auto builder = (TkTransientTrackingRecHitBuilder const *)(builderH.product());
159  cloner = (*builder).cloner();
160  }
161 
162  // get updator
163  KFUpdator updator;
164 
165  // Now update initial state track using information from seed hits.
166 
167  TrajectoryStateOnSurface updatedState;
168  edm::OwnVector<TrackingRecHit> seedHits;
169 
170  const TrackingRecHit* hit = 0;
171  for ( unsigned int iHit = 0; iHit < hits.size() && iHit<1; iHit++) {
172  hit = hits[iHit];
173  TrajectoryStateOnSurface state = (iHit==0) ?
174  propagator->propagate(fts,tracker->idToDet(hit->geographicalId())->surface())
175  : propagator->propagate(updatedState, tracker->idToDet(hit->geographicalId())->surface());
176  if (!state.isValid()) return 0;
177 
178  SeedingHitSet::ConstRecHitPointer tth = hits[iHit];
179 
180  std::unique_ptr<BaseTrackerRecHit> newtth(refitHit( tth, state));
181 
182 
183  if (!checkHit(state,&*newtth,es)) return 0;
184 
185  updatedState = updator.update(state, *newtth);
186  if (!updatedState.isValid()) return 0;
187 
188  seedHits.push_back(newtth.release());
189 #ifdef mydebug_seed
190  uint32_t detid = hit->geographicalId().rawId();
191  (*pss) << "\n[SeedForPhotonConversion1Leg] hit " << iHit;
192  po.print(*pss, detid);
193  (*pss) << " " << detid << "\t lp " << hit->localPosition()
194  << " tth " << tth->localPosition() << " newtth " << newtth->localPosition() << " state " << state.globalMomentum().perp();
195 #endif
196  }
197 
198 
199  PTrajectoryStateOnDet const & PTraj =
200  trajectoryStateTransform::persistentState(updatedState, hit->geographicalId().rawId());
201 
202  seedCollection.push_back( TrajectorySeed(PTraj,seedHits,alongMomentum));
203  return &seedCollection.back();
204 }
pss
std::pair< ALIstring, ALIstring > pss
Definition: Fit.h:27
PV3DBase::perp
T perp() const
Definition: PV3DBase.h:72
trajectoryStateTransform::persistentState
PTrajectoryStateOnDet persistentState(const TrajectoryStateOnSurface &ts, unsigned int detid)
Definition: TrajectoryStateTransform.cc:16
MixedLayerPairs_cfi.TTRHBuilder
tuple TTRHBuilder
Definition: MixedLayerPairs_cfi.py:34
DetId::rawId
uint32_t rawId() const
get the raw id
Definition: DetId.h:43
edm::OwnVector::push_back
void push_back(D *&d)
Definition: OwnVector.h:274
SeedingHitSet::ConstRecHitPointer
BaseTrackerRecHit const * ConstRecHitPointer
Definition: SeedingHitSet.h:11
patCandidatesForDimuonsSequences_cff.tracker
tuple tracker
Definition: patCandidatesForDimuonsSequences_cff.py:69
PrintRecoObjects::print
void print(std::stringstream &ss, const SiStripCluster &clus)
Definition: PrintRecoObjects.cc:9
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:55
edm::ESHandle::product
T const * product() const
Definition: ESHandle.h:62
Propagator::propagate
TrajectoryStateOnSurface propagate(STA const &state, SUR const &surface) const
Definition: Propagator.h:56
TrajectoryStateOnSurface::globalMomentum
GlobalVector globalMomentum() const
Definition: TrajectoryStateOnSurface.h:92
SeedingHitSet::size
unsigned int size() const
Definition: SeedingHitSet.h:44
SeedForPhotonConversion1Leg::refitHit
SeedingHitSet::RecHitPointer refitHit(SeedingHitSet::ConstRecHitPointer hit, const TrajectoryStateOnSurface &state) const
Definition: SeedForPhotonConversion1Leg.cc:206
TrackingRecHit::geographicalId
DetId geographicalId() const
Definition: TrackingRecHit.h:120
SeedForPhotonConversion1Leg::checkHit
bool checkHit(const TrajectoryStateOnSurface &, const SeedingHitSet::ConstRecHitPointer &hit, const edm::EventSetup &es) const
Definition: SeedForPhotonConversion1Leg.h:50
TrackingRecHit::localPosition
virtual LocalPoint localPosition() const =0
bool SeedForPhotonConversion1Leg::checkHit ( const TrajectoryStateOnSurface &  ,
const SeedingHitSet::ConstRecHitPointer &  hit,
const edm::EventSetup &  es 
) const
inlineprotected

Definition at line 50 of file SeedForPhotonConversion1Leg.h.

Referenced by buildSeed().

53  { return true; }
CurvilinearTrajectoryError SeedForPhotonConversion1Leg::initialError ( const GlobalVector &  vertexBounds,
float  ptMin,
float  sinTheta 
) const
protected

Definition at line 103 of file SeedForPhotonConversion1Leg.cc.

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

Referenced by trajectorySeed().

107 {
108  // Set initial uncertainty on track parameters, using only P.V. constraint and no hit
109  // information.
110  GlobalError vertexErr( sqr(vertexBounds.x()), 0,
111  sqr(vertexBounds.y()), 0, 0,
112  sqr(vertexBounds.z())
113  );
114 
115 
116  AlgebraicSymMatrix55 C = ROOT::Math::SMatrixIdentity();
117 
118 // FIXME: minC00. Prevent apriori uncertainty in 1/P from being too small,
119 // to avoid instabilities.
120 // N.B. This parameter needs optimising ...
121  float sin2th = sqr(sinTheta);
122  float minC00 = 1.0;
123  C[0][0] = std::max(sin2th/sqr(ptMin), minC00);
124  float zErr = vertexErr.czz();
125  float transverseErr = vertexErr.cxx(); // assume equal cxx cyy
126  C[3][3] = transverseErr;
127  C[4][4] = zErr*sin2th + transverseErr*(1-sin2th);
128 
129  return CurvilinearTrajectoryError(C);
130 }
PV3DBase::y
T y() const
Definition: PV3DBase.h:63
AlgebraicSymMatrix55
ROOT::Math::SMatrix< double, 5, 5, ROOT::Math::MatRepSym< double, 5 > > AlgebraicSymMatrix55
Definition: AlgebraicROOTObjects.h:31
max
const T & max(const T &a, const T &b)
Definition: MaterialBudgetTrackerHistos.cc:4
PV3DBase::z
T z() const
Definition: PV3DBase.h:64
PtMinSelector_cfg.ptMin
tuple ptMin
Definition: PtMinSelector_cfg.py:44
funct::sqr
Square< F >::type sqr(const F &f)
Definition: Square.h:13
PV3DBase::x
T x() const
Definition: PV3DBase.h:62
GlobalTrajectoryParameters SeedForPhotonConversion1Leg::initialKinematic ( const SeedingHitSet &  hits,
const GlobalPoint &  vertexPos,
const edm::EventSetup &  es,
const float  cotTheta 
) const
protected

Definition at line 44 of file SeedForPhotonConversion1Leg.cc.

References GlobalTrajectoryParameters::charge(), cotTheta_Max, cond::rpcobgas::detid, edm::EventSetup::get(), GlobalTrajectoryParameters::magneticField(), GlobalTrajectoryParameters::momentum(), 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().

49 {
50  GlobalTrajectoryParameters kine;
51 
52  SeedingHitSet::ConstRecHitPointer tth1 = hits[0];
53  SeedingHitSet::ConstRecHitPointer tth2 = hits[1];
54 
55  // FIXME optimize: move outside loop
56  edm::ESHandle<MagneticField> bfield;
57  es.get<IdealMagneticFieldRecord>().get(bfield);
58  float nomField = bfield->nominalValue();
59 
60  FastHelix helix(tth2->globalPosition(), tth1->globalPosition(), vertexPos, nomField, &*bfield, vertexPos);
61  kine = helix.stateAtVertex();
62 
63  //force the pz/pt equal to the measured one
64  if(fabs(cotTheta)<cotTheta_Max)
65  kine = GlobalTrajectoryParameters(kine.position(),
66  GlobalVector(kine.momentum().x(),kine.momentum().y(),kine.momentum().perp()*cotTheta),
67  kine.charge(),
68  & kine.magneticField()
69  );
70  else
71  kine = GlobalTrajectoryParameters(kine.position(),
72  GlobalVector(kine.momentum().x(),kine.momentum().y(),kine.momentum().perp()*cotTheta_Max),
73  kine.charge(),
74  & kine.magneticField()
75  );
76 
77 #ifdef mydebug_seed
78  uint32_t detid;
79  (*pss) << "[SeedForPhotonConversion1Leg] initialKinematic tth1 " ;
80  detid=tth1->geographicalId().rawId();
81  po.print(*pss, detid );
82  (*pss) << " \t " << detid << " " << tth1->localPosition() << " " << tth1->globalPosition() ;
83  detid= tth2->geographicalId().rawId();
84  (*pss) << " \n\t tth2 ";
85  po.print(*pss, detid );
86  (*pss) << " \t " << detid << " " << tth2->localPosition() << " " << tth2->globalPosition()
87  << "\nhelix momentum " << kine.momentum() << " pt " << kine.momentum().perp() << " radius " << 1/kine.transverseCurvature();
88 #endif
89 
90  bool isBOFF =(0==nomField);;
91  if (isBOFF && (theBOFFMomentum > 0)) {
92  kine = GlobalTrajectoryParameters(kine.position(),
93  kine.momentum().unit() * theBOFFMomentum,
94  kine.charge(),
95  &*bfield);
96  }
97  return kine;
98 }
pss
std::pair< ALIstring, ALIstring > pss
Definition: Fit.h:27
PV3DBase::perp
T perp() const
Definition: PV3DBase.h:72
PV3DBase::y
T y() const
Definition: PV3DBase.h:63
SeedingHitSet::ConstRecHitPointer
BaseTrackerRecHit const * ConstRecHitPointer
Definition: SeedingHitSet.h:11
PrintRecoObjects::print
void print(std::stringstream &ss, const SiStripCluster &clus)
Definition: PrintRecoObjects.cc:9
Vector3DBase::unit
Vector3DBase unit() const
Definition: Vector3DBase.h:57
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:55
GlobalTrajectoryParameters::magneticField
const MagneticField & magneticField() const
Definition: GlobalTrajectoryParameters.h:130
PV3DBase::x
T x() const
Definition: PV3DBase.h:62
GlobalVector
Global3DVector GlobalVector
Definition: GlobalVector.h:10
SeedingHitSet::RecHitPointer SeedForPhotonConversion1Leg::refitHit ( SeedingHitSet::ConstRecHitPointer  hit,
const TrajectoryStateOnSurface &  state 
) const
protected

Definition at line 206 of file SeedForPhotonConversion1Leg.cc.

References cloner.

Referenced by buildSeed().

209 {
210  //const TransientTrackingRecHit* a= hit.get();
211  //return const_cast<TransientTrackingRecHit*> (a);
212  //This was modified otherwise the rechit will have just the local x component and local y=0
213  // To understand how to modify for pixels
214 
215  //const TSiStripRecHit2DLocalPos* b = dynamic_cast<const TSiStripRecHit2DLocalPos*>(a);
216  //return const_cast<TSiStripRecHit2DLocalPos*>(b);
217  return (SeedingHitSet::RecHitPointer)(cloner(*hit,state));
218 }
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 22 of file SeedForPhotonConversion1Leg.cc.

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

30 {
31  pss = &ss;
32  if ( hits.size() < 2) return 0;
33 
34  GlobalTrajectoryParameters kine = initialKinematic(hits, vertex, es, cotTheta);
35  float sinTheta = sin(kine.momentum().theta());
36 
37  CurvilinearTrajectoryError error = initialError(vertexBounds, ptmin, sinTheta);
38  FreeTrajectoryState fts(kine, error);
39 
40  return buildSeed(seedCollection,hits,fts,es);
41 }
pss
std::pair< ALIstring, ALIstring > pss
Definition: Fit.h:27
SeedForPhotonConversion1Leg::buildSeed
const TrajectorySeed * buildSeed(TrajectorySeedCollection &seedCollection, const SeedingHitSet &hits, const FreeTrajectoryState &fts, const edm::EventSetup &es) const
Definition: SeedForPhotonConversion1Leg.cc:132
SeedForPhotonConversion1Leg::initialKinematic
GlobalTrajectoryParameters initialKinematic(const SeedingHitSet &hits, const GlobalPoint &vertexPos, const edm::EventSetup &es, const float cotTheta) const
Definition: SeedForPhotonConversion1Leg.cc:44
funct::sin
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
PV3DBase::theta
Geom::Theta< T > theta() const
Definition: PV3DBase.h:75
ptmin
double ptmin
Definition: HydjetWrapper.h:85
SeedingHitSet::size
unsigned int size() const
Definition: SeedingHitSet.h:44
SeedForPhotonConversion1Leg::initialError
CurvilinearTrajectoryError initialError(const GlobalVector &vertexBounds, float ptMin, float sinTheta) const
Definition: SeedForPhotonConversion1Leg.cc:103

Member Data Documentation

TkClonerImpl SeedForPhotonConversion1Leg::cloner
mutableprotected

Definition at line 80 of file SeedForPhotonConversion1Leg.h.

Referenced by buildSeed(), and refitHit().

const int SeedForPhotonConversion1Leg::cotTheta_Max =99999
static

Definition at line 24 of file SeedForPhotonConversion1Leg.h.

Referenced by initialKinematic().

PrintRecoObjects SeedForPhotonConversion1Leg::po
protected

Definition at line 83 of file SeedForPhotonConversion1Leg.h.

Referenced by buildSeed(), and initialKinematic().

std::stringstream* SeedForPhotonConversion1Leg::pss
protected

Definition at line 82 of file SeedForPhotonConversion1Leg.h.

double SeedForPhotonConversion1Leg::theBOFFMomentum
protected

Definition at line 77 of file SeedForPhotonConversion1Leg.h.

Referenced by initialKinematic().

std::string SeedForPhotonConversion1Leg::thePropagatorLabel
protected

Definition at line 76 of file SeedForPhotonConversion1Leg.h.

Referenced by buildSeed().

Generated for CMSSW Reference Manual by doxygen 1.8.5