|
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
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
|
#include "DataFormats/TrackerRecHit2D/interface/SiTrackerMultiRecHit.h"
using namespace std;
using namespace edm;
SiTrackerMultiRecHit::SiTrackerMultiRecHit(const LocalPoint& pos,
const LocalError& err,
GeomDet const& idet,
const std::vector<std::pair<const TrackingRecHit*, float> >& aHitMap,
double annealing)
: BaseTrackerRecHit(pos, err, idet, trackerHitRTTI::multi) {
for (std::vector<std::pair<const TrackingRecHit*, float> >::const_iterator ihit = aHitMap.begin();
ihit != aHitMap.end();
ihit++) {
theHits.push_back(ihit->first->clone());
theWeights.push_back(ihit->second);
}
annealing_ = annealing;
}
bool SiTrackerMultiRecHit::sharesInput(const TrackingRecHit* other, SharedInputType what) const {
if (geographicalId() != other->geographicalId() && what == all)
return false;
vector<const TrackingRecHit*> otherhits = other->recHits();
if (what == all) {
if (theHits.size() != other->recHits().size())
return false;
for (vector<const TrackingRecHit*>::iterator otherhit = otherhits.begin(); otherhit != otherhits.end();
++otherhit) {
bool found = false;
for (OwnVector<TrackingRecHit>::const_iterator hit = theHits.begin(); hit != theHits.end(); ++hit) {
if ((hit)->sharesInput(*otherhit, all)) {
found = true;
break;
}
}
if (found == false) {
return false;
}
}
return true;
} else {
for (OwnVector<TrackingRecHit>::const_iterator hit = theHits.begin(); hit != theHits.end(); ++hit) {
if (!otherhits.empty()) {
for (vector<const TrackingRecHit*>::iterator otherhit = otherhits.begin(); otherhit != otherhits.end();
++otherhit) {
if ((hit)->sharesInput(*otherhit, some))
return true;
}
} else { //otherwise it should be a single rechit
if ((hit)->sharesInput(other, some))
return true;
}
}
return false;
}
}
vector<const TrackingRecHit*> SiTrackerMultiRecHit::recHits() const {
vector<const TrackingRecHit*> myhits;
for (edm::OwnVector<TrackingRecHit>::const_iterator ihit = theHits.begin(); ihit != theHits.end(); ihit++) {
myhits.push_back(&*ihit);
}
return myhits;
}
vector<TrackingRecHit*> SiTrackerMultiRecHit::recHits() { return theHits.data(); }
int SiTrackerMultiRecHit::dimension() const {
//supposing all the hits inside of a MRH have the same id == same type
int randomComponent = 0;
if (theHits[randomComponent].dimension() == 1) {
return 1;
} else if (theHits[randomComponent].dimension() == 2) {
return 2;
} else {
return 0;
}
}
void SiTrackerMultiRecHit::getKfComponents(KfComponentsHolder& holder) const {
if (dimension() == 1)
getKfComponents1D(holder);
if (dimension() == 2)
getKfComponents2D(holder);
}
|