|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
|
#include "Alignment/ReferenceTrajectories/interface/DualBzeroReferenceTrajectory.h"
#include "TrackingTools/TrajectoryState/interface/TrajectoryStateOnSurface.h"
#include "DataFormats/CLHEP/interface/AlgebraicObjects.h"
#include "DataFormats/TrajectoryState/interface/LocalTrajectoryParameters.h"
#include "Alignment/ReferenceTrajectories/interface/BzeroReferenceTrajectory.h"
DualBzeroReferenceTrajectory::DualBzeroReferenceTrajectory(const TrajectoryStateOnSurface& tsos,
const ConstRecHitContainer& forwardRecHits,
const ConstRecHitContainer& backwardRecHits,
const MagneticField* magField,
const reco::BeamSpot& beamSpot,
const ReferenceTrajectoryBase::Config& config)
: DualReferenceTrajectory(tsos.localParameters().mixedFormatVector().kSize - 1,
numberOfUsedRecHits(forwardRecHits) + numberOfUsedRecHits(backwardRecHits) - 1,
config),
theMomentumEstimate(config.momentumEstimate) {
theValidityFlag = construct(tsos, forwardRecHits, backwardRecHits, magField, beamSpot);
}
ReferenceTrajectory* DualBzeroReferenceTrajectory::construct(const TrajectoryStateOnSurface& referenceTsos,
const ConstRecHitContainer& recHits,
double mass,
MaterialEffects materialEffects,
const PropagationDirection propDir,
const MagneticField* magField,
bool useBeamSpot,
const reco::BeamSpot& beamSpot) const {
if (materialEffects >= breakPoints)
throw cms::Exception("BadConfig") << "[DualBzeroReferenceTrajectory::construct] Wrong MaterialEffects: "
<< materialEffects;
ReferenceTrajectoryBase::Config config(materialEffects, propDir, mass, theMomentumEstimate);
config.useBeamSpot = useBeamSpot;
config.hitsAreReverse = false;
return new BzeroReferenceTrajectory(referenceTsos, recHits, magField, beamSpot, config);
}
AlgebraicVector DualBzeroReferenceTrajectory::extractParameters(const TrajectoryStateOnSurface& referenceTsos) const {
AlgebraicVector param = asHepVector<5>(referenceTsos.localParameters().mixedFormatVector());
return param.sub(2, 5);
}
|