Main function of the track finding by the cellular automaton.
92{
95
97
98 const std::map<ITrackType, CDCHitVector>& mcTracksByMCParticleIdx =
100
101 std::size_t nAxialTracks = mcTracksByMCParticleIdx.size();
102 outputAxialTracks.reserve(outputAxialTracks.size() + nAxialTracks);
103
104 for (const std::pair<ITrackType, CDCHitVector> mcTracksAndMCParticleIdx : mcTracksByMCParticleIdx) {
105
106 const CDCHitVector& mcTrack = mcTracksAndMCParticleIdx.second;
107
108 outputAxialTracks.push_back(CDCTrack());
109 CDCTrack& axialTrack = outputAxialTracks.back();
110 bool reachedOuterMostLayer = false;
111 for (const CDCHit* ptrHit : mcTrack) {
112
114
115
116
117 if (ptrHit->getISuperLayer() == 8 and ptrHit->getILayer() == 5) {
118 reachedOuterMostLayer = true;
119 }
120 if (reachedOuterMostLayer and ptrHit->getILayer() != 5) {
121 break;
122 }
123 }
124
125 const CDCWireHit* wireHit = simHitLookUp.
getWireHit(ptrHit, inputWireHits);
126 if (not wireHit) continue;
127
129 if (not recoHit2D.
isAxial())
continue;
130
132 CDCRecoHit3D recoHit3D(recoHit2D.
getRLWireHit(), {tmp.X(), tmp.Y(), 0}, NAN);
133 axialTrack.push_back(recoHit3D);
134 }
135
136
137 if (axialTrack.size() < 5) outputAxialTracks.pop_back();
138 }
139
140 CDC::RealisticTDCCountTranslator tdcCountTranslator;
141 for (CDCTrack& track : outputAxialTracks) {
142 for (CDCRecoHit3D& recoHit3D : track) {
143 ROOT::Math::XYVector recoPos2D = recoHit3D.getRecoPos2D();
144 ROOT::Math::XYVector flightDirection = recoHit3D.getFlightDirection2D();
145 double alpha = ROOT::Math::VectorUtil::DeltaPhi(recoPos2D, flightDirection);
146
147 const CDCWire& wire = recoHit3D.getWire();
148 const bool rl = recoHit3D.getRLInfo() == ERightLeft::c_Right;
149
150 double driftLength = std::fabs(recoHit3D.getSignedRecoDriftLength());
152
153
154 const double beta = 1;
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173 driftLength = recoHit3D.getWireHit().getRefDriftLength();
174
175 } else {
176
179 rl,
181 alpha);
182
183 driftLength += gRandom->Gaus(0, std::sqrt(driftLengthVariance));
184 }
185 bool snapRecoPos = true;
186 recoHit3D.setRecoDriftLength(driftLength, snapRecoPos);
187 }
188 }
189
191 CDCKarimakiFitter fitter;
192 for (CDCTrack& track : outputAxialTracks) {
193 CDCTrajectory2D trajectory2D = fitter.fit(track);
196 }
197 } else {
199 for (CDCTrack& track : outputAxialTracks) {
203 }
204 }
205
207 for (CDCTrack& track : outputAxialTracks) {
209 }
210 }
211}
const WireID & getWireID() const
Getter for the wire id.
double getRefZ() const
Getter for the wire reference z coordinate Gives the wire's reference z coordinate.
double getDriftLengthResolution(double driftLength, const WireID &wireID=WireID(), bool leftRight=false, double z=0, double alpha=0, double=static_cast< double >(TMath::Pi()/2.)) override
Get position resolution^2 corresponding to the drift length from getDriftLength of this class.
bool m_param_reconstructedDriftLength
Parameter : Setup the drift length as it can be estimated from two dimensional information.
bool m_param_reconstructedPositions
Parameter : Switch to reconstruct the positions in the tracks imitating the Legendre finder.
bool m_param_fit
Parameter : Fit the track instead of forwarding the MC truth information.
bool m_param_useOnlyBeforeTOP
Parameter : Cut tracks after the last layer of the CDC has been reached, assuming the tracks left the...
TrackingUtilities::CDCTrajectory3D getTrajectory3D(const ACDCHitCollection *ptrHits) const
Returns the trajectory of the collection of hits.
static const CDCMCTrackLookUp & getInstance()
Getter for the singletone instance.
static const CDCMCTrackStore & getInstance()
Getter for the singletone instance.
std::vector< const CDCHit * > CDCHitVector
Type for an ordered sequence of pointers to the CDCHit.
const std::map< ITrackType, Belle2::TrackFindingCDC::CDCMCTrackStore::CDCHitVector > & getMCTracksByMCParticleIdx() const
Getter for the stored Monte Carlo tracks ordered by their Monte Carlo Id.
static const CDCSimHitLookUp & getInstance()
Getter for the singletone instance.
TrackingUtilities::CDCRecoHit2D getClosestPrimaryRecoHit2D(const CDCHit *ptrHit, const std::vector< TrackingUtilities::CDCWireHit > &wireHits) const
Construct an TrackingUtilities::CDCRecoHit2D from the closest primary CDCSimHit information related t...
const TrackingUtilities::CDCWireHit * getWireHit(const CDCHit *ptrHit, const std::vector< TrackingUtilities::CDCWireHit > &wireHits) const
Retrieve the wire hit the given CDCHit form the given wire hits.
static const FlightTimeEstimator & instance(std::unique_ptr< FlightTimeEstimator > replacement=nullptr)
Getter for the instance.
virtual double getFlightTime2D(const ROOT::Math::XYVector &, double, double=1) const
Default estimator for the flight time.
const CDCRLWireHit & getRLWireHit() const
Getter for the oriented wire hit associated with the reconstructed hit.
bool isAxial() const
Indicator if the underlying wire is axial.
ROOT::Math::XYVector getRecoPos2D() const
Getter for the position in the reference plane.
double setLocalOrigin(const ROOT::Math::XYVector &localOrigin)
Setter for the origin of the local coordinate system.
CDCTrajectory2D getTrajectory2D() const
Getter for the two dimensional trajectory.
static CDCTrajectorySZ basicAssumption()
Constructs a basic assumption, what the z0 start position and the sz slope are, including some broad ...
static void normalizeTrack(TrackingUtilities::CDCTrack &track)
Refit and resort the track. Unmask all hits.
1.13.2