CMS 3D CMS Logo

All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
SeedGeneratorForCRack.cc
Go to the documentation of this file.
6  const SiStripRecHit2DCollection& collrphi,
7  const SiStripMatchedRecHit2DCollection& collmatched,
8  const edm::EventSetup& iSetup) {
10  iSetup.get<TrackerDigiGeometryRecord>().get(tracker);
13  theUpdator = new KFUpdator();
14 
15  // get the transient builder
16 
18 
19  iSetup.get<TransientRecHitRecord>().get(builderName, theBuilder);
20  TTTRHBuilder = theBuilder.product();
21  CosmicLayerPairs cosmiclayers(geometry);
22  cosmiclayers.init(collstereo, collrphi, collmatched, iSetup);
23  thePairGenerator = new CosmicHitPairGenerator(cosmiclayers, iSetup);
24  HitPairs.clear();
26  LogDebug("CosmicSeedFinder") << "Initialized with " << HitPairs.size() << " hit pairs" << std::endl;
27 }
28 
30  : conf_(conf),
31  region(conf.getParameter<double>("ptMin"),
32  conf.getParameter<double>("originRadius"),
33  conf.getParameter<double>("originHalfLength"),
34  conf.getParameter<double>("originZPosition")) {
35  seedpt = conf_.getParameter<double>("SeedPt");
36  builderName = conf_.getParameter<std::string>("TTRHBuilder");
37  geometry = conf_.getUntrackedParameter<std::string>("GeometricStructure", "STANDARD");
38  multipleScatteringFactor = conf_.getUntrackedParameter<double>("multipleScatteringFactor", 1.0);
39  seedMomentum = conf_.getUntrackedParameter<double>("SeedMomentum", 1);
40 }
41 
43  seeds(output, iSetup, region);
44  delete thePairGenerator;
45  delete thePropagatorAl;
46  delete thePropagatorOp;
47  delete theUpdator;
48 }
49 
51  const edm::EventSetup& iSetup,
52  const TrackingRegion& region) {
53  for (unsigned int is = 0; is < HitPairs.size(); is++) {
54  GlobalPoint inner = tracker->idToDet((*(HitPairs[is].inner())).geographicalId())
55  ->surface()
56  .toGlobal((*(HitPairs[is].inner())).localPosition());
57  GlobalPoint outer = tracker->idToDet((*(HitPairs[is].outer())).geographicalId())
58  ->surface()
59  .toGlobal((*(HitPairs[is].outer())).localPosition());
60 
61  LogDebug("CosmicSeedFinder") << "inner point of the seed " << inner << " outer point of the seed " << outer;
62  SeedingHitSet::ConstRecHitPointer inrhit = HitPairs[is].inner();
63  SeedingHitSet::ConstRecHitPointer outrhit = HitPairs[is].outer();
64 
66  hits.push_back(HitPairs[is].outer()->hit()->clone());
67 
68  for (int i = 0; i < 2; i++) {
69  //FIRST STATE IS CALCULATED CONSIDERING THAT THE CHARGE CAN BE POSITIVE OR NEGATIVE
70  int predsign = (2 * i) - 1;
71  if ((outer.y() - inner.y()) > 0) {
72  GlobalVector momentum = GlobalVector(inner - outer);
73  momentum = momentum.unit() * seedMomentum;
74  GlobalTrajectoryParameters Gtp(inner, momentum, predsign, &(*magfield));
75  AlgebraicSymMatrix55 errMatrix = ROOT::Math::SMatrixIdentity();
76  TSOS innerState = TSOS(Gtp,
77  CurvilinearTrajectoryError(errMatrix),
78  tracker->idToDet((HitPairs[is].inner()->hit())->geographicalId())->surface());
79  const TSOS innerUpdated = theUpdator->update(innerState, *inrhit);
80  //Cosmic Seed update inner...
81  LogDebug("CosmicSeedFinder") << " FirstTSOS " << innerUpdated;
82 
83  //First propagation
84  const TSOS outerState = thePropagatorOp->propagate(
85  innerUpdated, tracker->idToDet((*(HitPairs[is].outer())).geographicalId())->surface());
86  if (outerState.isValid()) {
87  LogDebug("CosmicSeedFinder") << "outerState " << outerState;
88  TSOS outerUpdated = theUpdator->update(outerState, *outrhit);
89  //fudge factor for multiple scattering
91  if (outerUpdated.isValid()) {
92  LogDebug("CosmicSeedFinder") << "outerUpdated " << outerUpdated;
93 
95  outerUpdated, (*(HitPairs[is].outer())).geographicalId().rawId());
96  output.push_back(TrajectorySeed(PTraj, hits, alongMomentum));
97 
98  } else
99  edm::LogWarning("CosmicSeedFinder") << " SeedForCosmics first update failed ";
100  } else
101  edm::LogWarning("CosmicSeedFinder") << " SeedForCosmics first propagation failed ";
102 
103  } else {
104  GlobalVector momentum = GlobalVector(outer - inner);
105  momentum = momentum.unit() * seedMomentum;
106  GlobalTrajectoryParameters Gtp(inner, momentum, predsign, &(*magfield));
107  AlgebraicSymMatrix55 errMatrix = ROOT::Math::SMatrixIdentity();
108  TSOS innerState = TSOS(Gtp,
109  CurvilinearTrajectoryError(errMatrix),
110  tracker->idToDet((HitPairs[is].inner()->hit())->geographicalId())->surface());
111  const TSOS innerUpdated = theUpdator->update(innerState, *inrhit);
112  LogDebug("CosmicSeedFinder") << " FirstTSOS " << innerState;
113 
114  //First propagation
115  const TSOS outerState = thePropagatorOp->propagate(
116  innerUpdated, tracker->idToDet((*(HitPairs[is].outer())).geographicalId())->surface());
117  if (outerState.isValid()) {
118  LogDebug("CosmicSeedFinder") << "outerState " << outerState;
119  TSOS outerUpdated = theUpdator->update(outerState, *outrhit);
120  //fudge factor for multiple scattering
122  if (outerUpdated.isValid()) {
123  LogDebug("CosmicSeedFinder") << "outerUpdated " << outerUpdated;
125  outerUpdated, (*(HitPairs[is].outer())).geographicalId().rawId());
126  output.push_back(TrajectorySeed(PTraj, hits, oppositeToMomentum));
127 
128  } else
129  edm::LogWarning("CosmicSeedFinder") << " SeedForCosmics first update failed ";
130  } else
131  edm::LogWarning("CosmicSeedFinder") << " SeedForCosmics first propagation failed ";
132  }
133  }
134  }
135 }
#define LogDebug(id)
GlobalPoint toGlobal(const Point2DBase< Scalar, LocalTag > lp) const
Definition: Surface.h:81
T getParameter(std::string const &) const
T getUntrackedParameter(std::string const &, T const &) const
const TransientTrackingRecHitBuilder * TTTRHBuilder
void seeds(TrajectorySeedCollection &output, const edm::EventSetup &c, const TrackingRegion &region)
PTrajectoryStateOnDet persistentState(const TrajectoryStateOnSurface &ts, unsigned int detid)
T y() const
Definition: PV3DBase.h:60
SeedGeneratorForCRack(const edm::ParameterSet &conf)
TrajectoryStateOnSurface TSOS
const Plane & surface() const
The nominal surface of the GeomDet.
Definition: GeomDet.h:37
void hitPairs(const TrackingRegion &reg, OrderedHitPairs &prs, const edm::EventSetup &iSetup)
form base class
void push_back(D *&d)
Definition: OwnVector.h:326
TrajectoryStateOnSurface update(const TrajectoryStateOnSurface &, const TrackingRecHit &) const override
Definition: KFUpdator.cc:177
BaseTrackerRecHit const * ConstRecHitPointer
Definition: SeedingHitSet.h:10
std::vector< TrajectorySeed > TrajectorySeedCollection
PropagatorWithMaterial * thePropagatorAl
unsigned int size() const override
edm::ESHandle< MagneticField > magfield
void init(const SiStripRecHit2DCollection &collstereo, const SiStripRecHit2DCollection &collrphi, const SiStripMatchedRecHit2DCollection &collmatched, const edm::EventSetup &iSetup)
Vector3DBase unit() const
Definition: Vector3DBase.h:54
PropagatorWithMaterial * thePropagatorOp
ROOT::Math::SMatrix< double, 5, 5, ROOT::Math::MatRepSym< double, 5 > > AlgebraicSymMatrix55
TEveGeoShape * clone(const TEveElement *element, TEveElement *parent)
Definition: eve_macros.cc:135
TrajectoryStateOnSurface propagate(STA const &state, SUR const &surface) const
Definition: Propagator.h:50
void init(const SiStripRecHit2DCollection &collstereo, const SiStripRecHit2DCollection &collrphi, const SiStripMatchedRecHit2DCollection &collmatched, const edm::EventSetup &c)
GlobalTrackingRegion region
T get() const
Definition: EventSetup.h:73
void run(TrajectorySeedCollection &, const edm::EventSetup &c)
const TrackerGeomDet * idToDet(DetId) const override
CosmicHitPairGenerator * thePairGenerator
T const * product() const
Definition: ESHandle.h:86
Global3DVector GlobalVector
Definition: GlobalVector.h:10
edm::ESHandle< TrackerGeometry > tracker