|
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
|
#ifndef DATAFORMATS_TRACKCANDIDATE_TRACKCANDIDATE_H
#define DATAFORMATS_TRACKCANDIDATE_TRACKCANDIDATE_H
#include "DataFormats/Common/interface/RefToBase.h"
#include "DataFormats/Common/interface/OwnVector.h"
#include "DataFormats/TrackCandidate/interface/TrajectoryStopReasons.h"
#include "DataFormats/TrackingRecHit/interface/TrackingRecHit.h"
#include "DataFormats/TrajectorySeed/interface/TrajectorySeed.h"
#include "FWCore/Utilities/interface/Range.h"
#include <utility>
/** A track candidate is
- a TSOS or equivalent (here a PTrajectoryStateOnDet)
- a vector of rechits (here via the OwnVector interface)
- a TrajectorySeed (to be confirmed as matching the final track)
- a reference to the TrajectorySeed in the origianl collection
of seeds. Often this collection is not saved on disk and
therefore the reference may be invalid.
only the second is compulsory,the other three can be empty / not present
**/
class TrackCandidate {
public:
typedef edm::OwnVector<TrackingRecHit> RecHitContainer;
TrackCandidate()
: rh_(), seed_(), state_(), seedRef_(), nLoops_(0), stopReason_((uint8_t)StopReason::UNINITIALIZED) {}
explicit TrackCandidate(RecHitContainer& rh)
: rh_(), seed_(), state_(), seedRef_(), nLoops_(0), stopReason_((uint8_t)StopReason::UNINITIALIZED) {
rh_.swap(rh);
}
TrackCandidate(RecHitContainer& rh,
TrajectorySeed const& s,
PTrajectoryStateOnDet const& st,
signed char nLoops = 0,
uint8_t stopReason = (uint8_t)StopReason::UNINITIALIZED)
: rh_(), seed_(s), state_(st), seedRef_(), nLoops_(nLoops), stopReason_(stopReason) {
rh_.swap(rh);
}
TrackCandidate(RecHitContainer& rh,
TrajectorySeed const& s,
PTrajectoryStateOnDet const& st,
const edm::RefToBase<TrajectorySeed>& seedRef,
signed char nLoops = 0,
uint8_t stopReason = (uint8_t)StopReason::UNINITIALIZED)
: rh_(), seed_(s), state_(st), seedRef_(seedRef), nLoops_(nLoops), stopReason_(stopReason) {
rh_.swap(rh);
}
PTrajectoryStateOnDet const& trajectoryStateOnDet() const { return state_; }
edm::Range<RecHitContainer::const_iterator> recHits() const { return {rh_.begin(), rh_.end()}; }
auto nRecHits() const { return rh_.size(); }
TrajectorySeed const& seed() const { return seed_; }
bool isLooper() const { return (nLoops_ > 0); }
signed char nLoops() const { return nLoops_; }
uint8_t stopReason() const { return stopReason_; }
void setNLoops(signed char value) { nLoops_ = value; }
void setStopReason(uint8_t value) { stopReason_ = value; }
/** return the edm::reference to the trajectory seed in the original
* seeds collection. If the collection has been dropped from the
* Event, the reference may be invalid. Its validity should be tested,
* before the reference is actually used.
*/
edm::RefToBase<TrajectorySeed> seedRef() const { return seedRef_; }
void setSeedRef(edm::RefToBase<TrajectorySeed>& seedRef) { seedRef_ = seedRef; }
private:
RecHitContainer rh_;
TrajectorySeed seed_;
PTrajectoryStateOnDet state_;
edm::RefToBase<TrajectorySeed> seedRef_;
signed char nLoops_;
uint8_t stopReason_;
};
#endif
|