BOSS 8.0.0
BESIII Offline Software System
Loading...
Searching...
No Matches
MdcTrack Class Reference

#include <MdcTrack.h>

Public Member Functions

 MdcTrack (TrkRecoTrk *aTrack)
 MdcTrack (int nsupers, const TrkExchangePar &par, double chisq, TrkContext &, double trackT0)
 ~MdcTrack ()
bool operator== (const MdcTrack &tk) const
TrkRecoTrk & track ()
const TrkRecoTrk & track () const
void setTrack (TrkRecoTrk *trk)
int hasCurled () const
const MdcLayer * firstLayer () const
const MdcLayer * lastLayer () const
void setHasCurled (bool c=true)
void setFirstLayer (const MdcLayer *l)
void setLastLayer (const MdcLayer *l)
int projectToR (double radius, BesAngle &phiIntersect, int lCurl=0) const
int projectToR (double radius, BesAngle &phiIntersect, double &arcLength, int lCurl=0) const
void storeTrack (int trackId, RecMdcTrackCol *trackList, RecMdcHitCol *hitList, int tkStat)
 MdcTrack (TrkRecoTrk *aTrack)
 MdcTrack (int nsupers, const TrkExchangePar &par, double chisq, TrkContext &, double trackT0)
 ~MdcTrack ()
bool operator== (const MdcTrack &tk) const
TrkRecoTrk & track ()
const TrkRecoTrk & track () const
void setTrack (TrkRecoTrk *trk)
int hasCurled () const
const MdcLayer * firstLayer () const
const MdcLayer * lastLayer () const
void setHasCurled (bool c=true)
void setFirstLayer (const MdcLayer *l)
void setLastLayer (const MdcLayer *l)
int projectToR (double radius, BesAngle &phiIntersect, int lCurl=0) const
int projectToR (double radius, BesAngle &phiIntersect, double &arcLength, int lCurl=0) const
void storeTrack (int trackId, RecMdcTrackCol *trackList, RecMdcHitCol *hitList, int tkStat)
 MdcTrack (TrkRecoTrk *aTrack)
 MdcTrack (int nsupers, const TrkExchangePar &par, double chisq, TrkContext &, double trackT0)
 ~MdcTrack ()
bool operator== (const MdcTrack &tk) const
TrkRecoTrk & track ()
const TrkRecoTrk & track () const
void setTrack (TrkRecoTrk *trk)
int hasCurled () const
const MdcLayer * firstLayer () const
const MdcLayer * lastLayer () const
void setHasCurled (bool c=true)
void setFirstLayer (const MdcLayer *l)
void setLastLayer (const MdcLayer *l)
int projectToR (double radius, BesAngle &phiIntersect, int lCurl=0) const
int projectToR (double radius, BesAngle &phiIntersect, double &arcLength, int lCurl=0) const
void storeTrack (int trackId, RecMdcTrackCol *trackList, RecMdcHitCol *hitList, int tkStat)

Detailed Description

Definition at line 19 of file InstallArea/x86_64-el9-gcc13-dbg/include/MdcTrkRecon/MdcTrack.h.

Constructor & Destructor Documentation

◆ MdcTrack() [1/6]

MdcTrack::MdcTrack ( TrkRecoTrk * aTrack)

Definition at line 41 of file MdcTrack.cxx.

41 {
42 //************************************************************************/
43 _theTrack = aTrack;
44 _firstLayer = _lastLayer = 0;
45 _haveCurled = 0;
46}

Referenced by operator==(), and storeTrack().

◆ MdcTrack() [2/6]

MdcTrack::MdcTrack ( int nsupers,
const TrkExchangePar & par,
double chisq,
TrkContext & context,
double trackT0 )

Definition at line 49 of file MdcTrack.cxx.

50 {
51 //************************************************************************/
52 TrkCircleMaker maker;
53 _theTrack = maker.makeTrack( par, chisq, context, trackT0 );
54 _firstLayer = _lastLayer = 0;
55 _haveCurled = 0;
56}
TrkCircleMaker
TrkSimpleMaker< TrkCircleRep > TrkCircleMaker
Definition InstallArea/x86_64-el9-gcc13-dbg/include/TrkFitter/TrkCircleMaker.h:14
TrkSimpleMaker::makeTrack
TrkRecoTrk * makeTrack(const TrkExchangePar &helix, const double chi2, const TrkContext &, double trackT0) const

◆ ~MdcTrack() [1/3]

MdcTrack::~MdcTrack ( )

Definition at line 59 of file MdcTrack.cxx.

59 {
60 //************************************************************************/
61 delete _theTrack;
62}

◆ MdcTrack() [3/6]

MdcTrack::MdcTrack ( TrkRecoTrk * aTrack)

◆ MdcTrack() [4/6]

MdcTrack::MdcTrack ( int nsupers,
const TrkExchangePar & par,
double chisq,
TrkContext & ,
double trackT0 )

◆ ~MdcTrack() [2/3]

MdcTrack::~MdcTrack ( )

◆ MdcTrack() [5/6]

MdcTrack::MdcTrack ( TrkRecoTrk * aTrack)

◆ MdcTrack() [6/6]

MdcTrack::MdcTrack ( int nsupers,
const TrkExchangePar & par,
double chisq,
TrkContext & ,
double trackT0 )

◆ ~MdcTrack() [3/3]

MdcTrack::~MdcTrack ( )

Member Function Documentation

◆ firstLayer() [1/3]

const MdcLayer * MdcTrack::firstLayer ( ) const
inline

Definition at line 31 of file InstallArea/x86_64-el9-gcc13-dbg/include/MdcTrkRecon/MdcTrack.h.

31{ return _firstLayer; }

Referenced by MdcTrackList::pickHits(), and MdcSegGrouper::transferHits().

◆ firstLayer() [2/3]

const MdcLayer * MdcTrack::firstLayer ( ) const
inline

Definition at line 31 of file InstallArea/x86_64-el9-gcc13-opt/include/MdcTrkRecon/MdcTrack.h.

31{ return _firstLayer; }

◆ firstLayer() [3/3]

const MdcLayer * MdcTrack::firstLayer ( ) const
inline

Definition at line 31 of file Reconstruction/MdcPatRec/MdcTrkRecon/include/MdcTrkRecon/MdcTrack.h.

31{ return _firstLayer; }

◆ hasCurled() [1/3]

int MdcTrack::hasCurled ( ) const
inline

Definition at line 30 of file InstallArea/x86_64-el9-gcc13-dbg/include/MdcTrkRecon/MdcTrack.h.

30{ return _haveCurled; }

Referenced by MdcTrackList::pickHits(), and MdcSegGrouper::transferHits().

◆ hasCurled() [2/3]

int MdcTrack::hasCurled ( ) const
inline

Definition at line 30 of file InstallArea/x86_64-el9-gcc13-opt/include/MdcTrkRecon/MdcTrack.h.

30{ return _haveCurled; }

◆ hasCurled() [3/3]

int MdcTrack::hasCurled ( ) const
inline

Definition at line 30 of file Reconstruction/MdcPatRec/MdcTrkRecon/include/MdcTrkRecon/MdcTrack.h.

30{ return _haveCurled; }

◆ lastLayer() [1/3]

const MdcLayer * MdcTrack::lastLayer ( ) const
inline

Definition at line 32 of file InstallArea/x86_64-el9-gcc13-dbg/include/MdcTrkRecon/MdcTrack.h.

32{ return _lastLayer; }

Referenced by MdcTrackList::pickHits(), and MdcSegGrouper::transferHits().

◆ lastLayer() [2/3]

const MdcLayer * MdcTrack::lastLayer ( ) const
inline

Definition at line 32 of file InstallArea/x86_64-el9-gcc13-opt/include/MdcTrkRecon/MdcTrack.h.

32{ return _lastLayer; }

◆ lastLayer() [3/3]

const MdcLayer * MdcTrack::lastLayer ( ) const
inline

Definition at line 32 of file Reconstruction/MdcPatRec/MdcTrkRecon/include/MdcTrkRecon/MdcTrack.h.

32{ return _lastLayer; }

◆ operator==() [1/3]

bool MdcTrack::operator== ( const MdcTrack & tk) const

Definition at line 138 of file MdcTrack.cxx.

138 {
139 //----------------------------------------------------------------
140
141 return ( track() == tk.track() );
142}

◆ operator==() [2/3]

bool MdcTrack::operator== ( const MdcTrack & tk) const

◆ operator==() [3/3]

bool MdcTrack::operator== ( const MdcTrack & tk) const

◆ projectToR() [1/6]

int MdcTrack::projectToR ( double radius,
BesAngle & phiIntersect,
double & arcLength,
int lCurl = 0 ) const

Definition at line 97 of file MdcTrack.cxx.

98 {
99 //************************************************************************/
100
101 // return -1 if no intersection
102 // Only valid for projecting track outwards from point of closest approach
103 // Returns arclength from POCA
104
105 const TrkFit* tkFit = track().fitResult();
106 if ( tkFit == 0 ) return -1;
107 TrkExchangePar par = tkFit->helix( 0.0 );
108 double d0 = par.d0();
109 double omega = par.omega();
110 double phi0 = par.phi0();
111 double r2d2 = radius * radius - d0 * d0;
112 if ( r2d2 < 0 ) return -1;
113 double rinv = 1. / radius;
114 double k2dinv = 1. / ( 1 + omega * d0 );
115 if ( k2dinv < 0.0 ) return -1;
116 double arg;
117
118 if ( !lCurl )
119 {
120 arg = d0 * rinv + 0.5 * omega * rinv * r2d2 * k2dinv;
121 if ( fabs( arg ) > 1.0 ) return -1;
122 phiIntersect.setRad( phi0 + asin( arg ) );
123 }
124 else
125 {
126 arg = -d0 * rinv - 0.5 * omega * rinv * r2d2 * k2dinv;
127 if ( fabs( arg ) > 1.0 ) return -1;
128 phiIntersect.setRad( phi0 + Constants::pi + asin( arg ) );
129 }
130
131 arg = 0.5 * omega * sqrt( r2d2 * k2dinv );
132 arcLength = 2. * asin( arg ) / omega;
133
134 return 0;
135}
arg
double arg(const EvtComplex &c)
Definition EvtComplex.hh:195
TrkFit::helix
virtual TrkExchangePar helix(double fltL) const =0
TrkRecoTrk::fitResult
const TrkFit * fitResult() const
Definition TrkRecoTrk.cxx:312

◆ projectToR() [2/6]

int MdcTrack::projectToR ( double radius,
BesAngle & phiIntersect,
double & arcLength,
int lCurl = 0 ) const

◆ projectToR() [3/6]

int MdcTrack::projectToR ( double radius,
BesAngle & phiIntersect,
double & arcLength,
int lCurl = 0 ) const

◆ projectToR() [4/6]

int MdcTrack::projectToR ( double radius,
BesAngle & phiIntersect,
int lCurl = 0 ) const

Definition at line 64 of file MdcTrack.cxx.

64 {
65 //************************************************************************/
66 // note that these are currently circle-only routines
67
68 // return -1 if no intersection
69
70 const TrkFit* tkFit = track().fitResult();
71 if ( tkFit == 0 ) return -1;
72 TrkExchangePar par = tkFit->helix( 0.0 );
73 double d0 = par.d0();
74 double omega = par.omega();
75 double phi0 = par.phi0();
76 double r2d2 = radius * radius - d0 * d0;
77 if ( r2d2 < 0 ) return -1;
78 double rinv = 1. / radius;
79 double k2dinv = 1. / ( 1 + omega * d0 );
80 if ( !lCurl )
81 {
82 double arg = d0 * rinv + 0.5 * omega * rinv * r2d2 * k2dinv;
83 if ( fabs( arg ) > 1.0 ) return -1;
84 phiIntersect.setRad( phi0 + asin( arg ) );
85 }
86 else
87 {
88 double arg = -d0 * rinv - 0.5 * omega * rinv * r2d2 * k2dinv;
89 if ( fabs( arg ) > 1.0 ) return -1;
90 phiIntersect.setRad( phi0 + Constants::pi + asin( arg ) );
91 }
92
93 return 0;
94}

Referenced by MdcTrackList::pickHits().

◆ projectToR() [5/6]

int MdcTrack::projectToR ( double radius,
BesAngle & phiIntersect,
int lCurl = 0 ) const

◆ projectToR() [6/6]

int MdcTrack::projectToR ( double radius,
BesAngle & phiIntersect,
int lCurl = 0 ) const

◆ setFirstLayer() [1/3]

void MdcTrack::setFirstLayer ( const MdcLayer * l)
inline

Definition at line 35 of file InstallArea/x86_64-el9-gcc13-dbg/include/MdcTrkRecon/MdcTrack.h.

35{ _firstLayer = l; }

Referenced by MdcTrackList::pickHits(), and MdcSegGrouper::transferHits().

◆ setFirstLayer() [2/3]

void MdcTrack::setFirstLayer ( const MdcLayer * l)
inline

Definition at line 35 of file InstallArea/x86_64-el9-gcc13-opt/include/MdcTrkRecon/MdcTrack.h.

35{ _firstLayer = l; }

◆ setFirstLayer() [3/3]

void MdcTrack::setFirstLayer ( const MdcLayer * l)
inline

Definition at line 35 of file Reconstruction/MdcPatRec/MdcTrkRecon/include/MdcTrkRecon/MdcTrack.h.

35{ _firstLayer = l; }

◆ setHasCurled() [1/3]

void MdcTrack::setHasCurled ( bool c = true)
inline

Definition at line 34 of file InstallArea/x86_64-el9-gcc13-dbg/include/MdcTrkRecon/MdcTrack.h.

34{ _haveCurled = c; }

Referenced by MdcTrackList::pickHits().

◆ setHasCurled() [2/3]

void MdcTrack::setHasCurled ( bool c = true)
inline

Definition at line 34 of file InstallArea/x86_64-el9-gcc13-opt/include/MdcTrkRecon/MdcTrack.h.

34{ _haveCurled = c; }

◆ setHasCurled() [3/3]

void MdcTrack::setHasCurled ( bool c = true)
inline

Definition at line 34 of file Reconstruction/MdcPatRec/MdcTrkRecon/include/MdcTrkRecon/MdcTrack.h.

34{ _haveCurled = c; }

◆ setLastLayer() [1/3]

void MdcTrack::setLastLayer ( const MdcLayer * l)
inline

Definition at line 36 of file InstallArea/x86_64-el9-gcc13-dbg/include/MdcTrkRecon/MdcTrack.h.

36{ _lastLayer = l; }

Referenced by MdcTrackList::pickHits(), and MdcSegGrouper::transferHits().

◆ setLastLayer() [2/3]

void MdcTrack::setLastLayer ( const MdcLayer * l)
inline

Definition at line 36 of file InstallArea/x86_64-el9-gcc13-opt/include/MdcTrkRecon/MdcTrack.h.

36{ _lastLayer = l; }

◆ setLastLayer() [3/3]

void MdcTrack::setLastLayer ( const MdcLayer * l)
inline

Definition at line 36 of file Reconstruction/MdcPatRec/MdcTrkRecon/include/MdcTrkRecon/MdcTrack.h.

36{ _lastLayer = l; }

◆ setTrack() [1/3]

void MdcTrack::setTrack ( TrkRecoTrk * trk)
inline

Definition at line 29 of file InstallArea/x86_64-el9-gcc13-dbg/include/MdcTrkRecon/MdcTrack.h.

29{ _theTrack = trk; }

Referenced by MdcTrackListBase::arbitrateHits().

◆ setTrack() [2/3]

void MdcTrack::setTrack ( TrkRecoTrk * trk)
inline

Definition at line 29 of file InstallArea/x86_64-el9-gcc13-opt/include/MdcTrkRecon/MdcTrack.h.

29{ _theTrack = trk; }

◆ setTrack() [3/3]

void MdcTrack::setTrack ( TrkRecoTrk * trk)
inline

Definition at line 29 of file Reconstruction/MdcPatRec/MdcTrkRecon/include/MdcTrkRecon/MdcTrack.h.

29{ _theTrack = trk; }

◆ storeTrack() [1/3]

void MdcTrack::storeTrack ( int trackId,
RecMdcTrackCol * trackList,
RecMdcHitCol * hitList,
int tkStat )

Definition at line 145 of file MdcTrack.cxx.

146 {
147 // check the fit succeeded
148 // if (status().failure()) { return; }
149
150 // get the results of the fit to this track
151 const TrkFit* fitresult = track().fitResult();
152 // check the fit worked
153 if ( fitresult == 0 ) return;
154
155 // new a Rec. Track MdcTrack
156 RecMdcTrack* recMdcTrack = new RecMdcTrack();
157 const TrkHitList* aList = track().hits();
158 // some track info
159 const BField& theField = track().bField();
160 double Bz = theField.bFieldZ();
161 // Fit info
162 int hitId = 0;
163 int nHits = 0;
164 int nSt = 0;
165 // int nAct=0; int nSeg=0;
166 // int maxlay = 0; int minlay = 43; int seg[11];/* 11 slayer */ int segme;
167 //****************************
168 //* active flag open this
169 //****************************
170 //* if (allHit>0){
171 //* nHits = aList->nHit();//add inActive hit also
172 //* }else{
173 //* nHits = fitresult->nMdc();// = nActive()
174 //* }
175 //****************************
176 // for 5.1.0 use all hits
177 nHits = aList->nHit();
178 //****************************
179 // use active only
180 // nHits = fitresult->nMdc();// = nActive()
181 //****************************
182
183 int q = fitresult->charge();
184 double chisq = fitresult->chisq();
185 int nDof = fitresult->nDof();
186 // track parameters at closest approach to beamline
187 double fltLenPoca = 0.0;
188 TrkExchangePar helix = fitresult->helix( fltLenPoca );
189 // std::cout<< __FILE__ << " phi0 " << helix.phi0() << " omega "
190 // <<helix.omega()<<std::endl;
191 double phi0 = helix.phi0();
192 double tanDip = helix.tanDip();
193 // double z0 = helix.z0();
194 double d0 = helix.d0();
195 // momenta and position
196 // CLHEP::Hep3Vector mom = fitresult->momentum(fltLenPoca);
197 HepPoint3D poca = fitresult->position( fltLenPoca );
198
199 double helixPar[5];
200 // dro =-d0
201 helixPar[0] = -d0;
202 // phi0 = phi0 - pi/2 [0,2pi)
203 double tphi0 = phi0 - Constants::pi / 2.;
204 if ( tphi0 < 0. ) tphi0 += Constants::pi * 2.;
205 helixPar[1] = tphi0;
206 // kappa = q/pxy
207 double pxy = fitresult->pt();
208 if ( pxy == 0. ) helixPar[2] = 9999.;
209 else helixPar[2] = q / fabs( pxy );
210 if ( pxy > 9999. ) helixPar[2] = 0.00001;
211 // dz = z0
212 helixPar[3] = helix.z0();
213 // tanl
214 helixPar[4] = tanDip;
215 // error
216 // V(Y)=S * V(X) * ST , mS = S , mVy = V(Y) , helix.covariance() = V(X)
217 HepSymMatrix mS( helix.params().num_row(), 0 );
218 mS[0][0] = -1.; // dr0=-d0
219 mS[1][1] = 1.;
220 mS[2][2] = -333.567 / Bz; // pxy -> cpa
221 mS[3][3] = 1.;
222 mS[4][4] = 1.;
223 HepSymMatrix mVy = helix.covariance().similarity( mS );
224 double errorMat[15];
225 int k = 0;
226 for ( int ie = 0; ie < 5; ie++ )
227 {
228 for ( int je = ie; je < 5; je++ )
229 {
230 errorMat[k] = mVy[ie][je];
231 k++;
232 }
233 }
234 double p, px, py, pz;
235 px = pxy * ( -sin( helixPar[1] ) );
236 py = pxy * cos( helixPar[1] );
237 pz = pxy * helixPar[4];
238 p = sqrt( pxy * pxy + pz * pz );
239 // theta, phi
240 double theta = acos( pz / p );
241 double phi = atan2( py, px );
242 recMdcTrack->setTrackId( trackId );
243 recMdcTrack->setHelix( helixPar );
244 recMdcTrack->setCharge( q );
245 recMdcTrack->setPxy( pxy );
246 recMdcTrack->setPx( px );
247 recMdcTrack->setPy( py );
248 recMdcTrack->setPz( pz );
249 recMdcTrack->setP( p );
250 recMdcTrack->setTheta( theta );
251 recMdcTrack->setPhi( phi );
252 recMdcTrack->setPoca( poca );
253 recMdcTrack->setX( poca.x() ); // poca
254 recMdcTrack->setY( poca.y() );
255 recMdcTrack->setZ( poca.z() );
256 recMdcTrack->setR( sqrt( poca.x() * poca.x() + poca.y() * poca.y() ) );
257 HepPoint3D pivot( 0., 0., 0. );
258 recMdcTrack->setPivot( pivot );
259 recMdcTrack->setVX0( 0. ); // pivot
260 recMdcTrack->setVY0( 0. );
261 recMdcTrack->setVZ0( 0. );
262 recMdcTrack->setError( errorMat );
263 recMdcTrack->setError( mVy );
264 recMdcTrack->setChi2( chisq );
265 recMdcTrack->setNdof( nDof );
266 recMdcTrack->setStat( tkStat );
267 recMdcTrack->setNhits( nHits );
268 //-----------hit list-------------
269 HitRefVec hitRefVec;
270 vector<int> vecHits;
271 // for fiterm
272 int maxLayerId = -1;
273 int minLayerId = 43;
274 double fiTerm = 999.;
275 const MdcRecoHitOnTrack* fiTermHot = NULL;
276 TrkHitList::hot_iterator hot = aList->begin();
277 for ( ; hot != aList->end(); hot++ )
278 {
279 const MdcRecoHitOnTrack* recoHot;
280 recoHot = dynamic_cast<const MdcRecoHitOnTrack*>( &( *hot ) );
281 // if (!recoHot->mdcHit()) return;
282 RecMdcHit* recMdcHit = new RecMdcHit;
283 // id
284 recMdcHit->setId( hitId );
285 //---------BESIII----------
286 // phiWire - phiHit <0 flagLR=0 left driftleft<0 ;
287 // phiWire - phiHit >0 flagLR=1 right driftright>0;
288 // flag = 2 ambig;
289 //-----Babar wireAmbig----
290 //-1 -> right, 0 -> don't know, +1 -> left
291 //+1 phiWire-phiHit<0 Left,
292 //-1 phiWire-phiHit>0 Right,
293 // 0 don't know
294 // ambig() w.r.t track direction
295 // wireAmbig() w.r.t. wire location
296 double hotWireAmbig = recoHot->wireAmbig();
297 double driftDist = fabs( recoHot->drift() );
298 double sigma = recoHot->hitRms();
299 double doca = fabs( recoHot->dcaToWire() );
300 int flagLR = 2;
301 if ( hotWireAmbig == 1 )
302 {
303 flagLR = 0; // left driftDist <0
304 doca *= -1.; // 2012-07-19
305 }
306 else if ( hotWireAmbig == -1 )
307 {
308 flagLR = 1; // right driftDist >0
309 }
310 else if ( hotWireAmbig == 0 )
311 {
312 flagLR = 2; // don't know
313 }
314 recMdcHit->setFlagLR( flagLR );
315 recMdcHit->setDriftDistLeft( ( -1 ) * driftDist );
316 recMdcHit->setErrDriftDistLeft( sigma );
317 recMdcHit->setDriftDistRight( driftDist );
318 recMdcHit->setErrDriftDistRight( sigma );
319 // Mdc Id
320 Identifier mdcId = recoHot->mdcHit()->digi()->identify();
321 recMdcHit->setMdcId( mdcId );
322 // corrected ADC
323
324 // contribution to chisquare
325 // fill chisq
326 double res = 999., rese = 999.;
327 if ( recoHot->resid( res, rese, false ) ) {}
328 else {}
329 double deltaChi = 0;
330 recoHot->getFitStuff( deltaChi ); // yzhang open 2010-09-20
331 recMdcHit->setChisqAdd( deltaChi * deltaChi );
332 // set tracks containing this hit
333 recMdcHit->setTrkId( trackId );
334 // doca
335 recMdcHit->setDoca( doca ); // doca sign left<0
336 // entrance angle
337
338 recMdcHit->setEntra( recoHot->entranceAngle() );
339 // z of hit
340 HepPoint3D pos;
341 Hep3Vector dir;
342 double fltLen = recoHot->fltLen();
343 recMdcHit->setFltLen( fltLen );
344 recoHot->getParentRep()->traj().getInfo( fltLen, pos, dir );
345 recMdcHit->setZhit( pos.z() );
346 double tof = recoHot->getParentRep()->arrivalTime( fltLen );
347 recMdcHit->setTdc( recoHot->mdcHit()->tdcIndex() );
348 recMdcHit->setAdc( recoHot->mdcHit()->adcIndex() );
349 // drift time
350 recMdcHit->setDriftT( recoHot->mdcHit()->driftTime( tof, pos.z() ) );
351 // for fiterm
352 int layerId = recoHot->mdcHit()->layernumber();
353 int wireId = recoHot->mdcHit()->wirenumber();
354 // std::cout<<" MdcTrack::store() ("<<layerId<<","<<wireId<<") fltLen "<<fltLen<<"
355 // wireAmbig "<<hotWireAmbig<<" y "<<pos.y()<<std::endl;
356 //<<recoHot->entranceAngle()<<std::endl;
357 if ( layerId >= maxLayerId )
358 {
359 maxLayerId = layerId;
360 fiTermHot = recoHot;
361 }
362 if ( layerId < minLayerId ) { minLayerId = layerId; }
363 // status flag
364 if ( recoHot->isActive() ) { recMdcHit->setStat( 1 ); }
365 else { recMdcHit->setStat( 0 ); }
366 // for 5.1.0 use all hits
367 if ( recoHot->layer()->view() ) { ++nSt; }
368 hitList->push_back( recMdcHit );
369 SmartRef<RecMdcHit> refHit( recMdcHit );
370 hitRefVec.push_back( refHit );
371 vecHits.push_back( mdcId.get_value() );
372 ++hitId;
373 }
374 // fi terminal angle
375 if ( fiTermHot != NULL )
376 { fiTerm = ( 1. / sqrt( 1. + tanDip * tanDip ) ) * fiTermHot->fltLen() * helix.omega(); }
377 recMdcTrack->setFiTerm( fiTerm );
378 // number of stereo hits contained
379 recMdcTrack->setNster( nSt );
380 recMdcTrack->setFirstLayer( maxLayerId );
381 recMdcTrack->setLastLayer( minLayerId );
382 recMdcTrack->setVecHits( hitRefVec );
383 trackList->push_back( recMdcTrack );
384} // end of storeTrack
HepPoint3D
HepGeom::Point3D< double > HepPoint3D
Definition CosmicGenerator.h:36
sin
double sin(const BesAngle a)
Definition InstallArea/x86_64-el9-gcc13-dbg/include/MdcGeom/BesAngle.h:185
cos
double cos(const BesAngle a)
Definition InstallArea/x86_64-el9-gcc13-dbg/include/MdcGeom/BesAngle.h:187
HitRefVec
SmartRefVector< RecMdcHit > HitRefVec
Definition InstallArea/x86_64-el9-gcc13-dbg/include/MdcRecEvent/RecMdcTrack.h:22
q
****INTEGER imax DOUBLE PRECISION m_pi *DOUBLE PRECISION m_amfin DOUBLE PRECISION m_Chfin DOUBLE PRECISION m_Xenph DOUBLE PRECISION m_sinw2 DOUBLE PRECISION m_GFermi DOUBLE PRECISION m_MfinMin DOUBLE PRECISION m_ta2 INTEGER m_out INTEGER m_KeyFSR INTEGER m_KeyQCD *COMMON c_Semalib $ !copy of input $ !CMS energy $ !beam mass $ !final mass $ !beam charge $ !final charge $ !smallest final mass $ !Z mass $ !Z width $ !EW mixing angle $ !Gmu Fermi $ alphaQED at q
Definition KKsem.h:33
DstMdcTrack::setError
void setError(double err[15])
Definition DstMdcTrack.cxx:112
DstMdcTrack::setHelix
void setHelix(double helix[5])
Definition DstMdcTrack.cxx:96
DstMdcTrack::setPoca
void setPoca(double poca[3])
Definition DstMdcTrack.cxx:104
MdcHitOnTrack::wireAmbig
int wireAmbig() const
Definition MdcHitOnTrack.cxx:359
MdcHitOnTrack::dcaToWire
double dcaToWire() const
Definition MdcHitOnTrack.cxx:112
MdcHitOnTrack::entranceAngle
double entranceAngle() const
Definition MdcHitOnTrack.cxx:158
MdcHitOnTrack::layer
const MdcLayer * layer() const
Definition MdcHitOnTrack.cxx:330
MdcHit::driftTime
double driftTime(double tof, double z) const
Definition MdcHit.cxx:142
MdcRecoHitOnTrack::mdcHit
const MdcHit * mdcHit() const
Definition MdcRecoHitOnTrack.cxx:29
RawData::identify
virtual Identifier identify() const
Definition RawData.cxx:15
Trajectory::getInfo
virtual void getInfo(double fltLen, HepPoint3D &pos, Hep3Vector &direction) const =0
TrkAbsFit::pt
virtual double pt(double fltL=0.) const =0
TrkAbsFit::chisq
virtual double chisq() const =0
TrkAbsFit::charge
virtual int charge() const =0
TrkAbsFit::nDof
virtual int nDof() const =0
TrkAbsFit::traj
virtual const TrkDifTraj & traj() const =0
TrkAbsFit::position
virtual HepPoint3D position(double fltL) const =0
TrkHitOnTrk::resid
double resid(bool exclude=false) const
Definition TrkHitOnTrk.cxx:153
TrkHitOnTrk::getFitStuff
TrkErrCode getFitStuff(HepVector &derivs, double &deltaChi) const
Definition TrkHitOnTrk.cxx:201
TrkRep::arrivalTime
virtual double arrivalTime(double fltL) const
Definition TrkRep.cxx:158

Referenced by MdcMergeDups::store(), MdcTrackListBase::store(), and MdcxCosmicSewer::store().

◆ storeTrack() [2/3]

void MdcTrack::storeTrack ( int trackId,
RecMdcTrackCol * trackList,
RecMdcHitCol * hitList,
int tkStat )

◆ storeTrack() [3/3]

void MdcTrack::storeTrack ( int trackId,
RecMdcTrackCol * trackList,
RecMdcHitCol * hitList,
int tkStat )

◆ track() [1/6]

TrkRecoTrk & MdcTrack::track ( )
inline

Definition at line 27 of file InstallArea/x86_64-el9-gcc13-dbg/include/MdcTrkRecon/MdcTrack.h.

27{ return *_theTrack; }

Referenced by MdcTrackListBase::arbitrateHits(), MdcSegGrouper::combineSegs(), MdcTrackList::createFromSegs(), MdcTrackListCsmc::createFromSegs(), MdcTrackList::dropMultiHotInLayer(), MdcTrackList::dumpAxCombine(), MdcTrackList::dumpAxFill(), MdcTrackList::dumpCircle(), MdcTrackList::dumpHelix(), MdcTrackList::dumpStCombine(), MdcTrkRecon::fillTrackList(), MdcSegGrouperSt::fillWithSegs(), MdcTrackList::finishCircle(), MdcTrackList::finishHelix(), operator==(), MdcTrackList::pickHits(), MdcTrackListBase::plot(), projectToR(), projectToR(), MdcSegGrouperSt::storePar(), storeTrack(), and MdcSegGrouper::transferHits().

◆ track() [2/6]

TrkRecoTrk & MdcTrack::track ( )
inline

Definition at line 27 of file InstallArea/x86_64-el9-gcc13-opt/include/MdcTrkRecon/MdcTrack.h.

27{ return *_theTrack; }

◆ track() [3/6]

TrkRecoTrk & MdcTrack::track ( )
inline

Definition at line 27 of file Reconstruction/MdcPatRec/MdcTrkRecon/include/MdcTrkRecon/MdcTrack.h.

27{ return *_theTrack; }

◆ track() [4/6]

const TrkRecoTrk & MdcTrack::track ( ) const
inline

Definition at line 28 of file InstallArea/x86_64-el9-gcc13-dbg/include/MdcTrkRecon/MdcTrack.h.

28{ return *_theTrack; }

◆ track() [5/6]

const TrkRecoTrk & MdcTrack::track ( ) const
inline

Definition at line 28 of file InstallArea/x86_64-el9-gcc13-opt/include/MdcTrkRecon/MdcTrack.h.

28{ return *_theTrack; }

◆ track() [6/6]

const TrkRecoTrk & MdcTrack::track ( ) const
inline

Definition at line 28 of file Reconstruction/MdcPatRec/MdcTrkRecon/include/MdcTrkRecon/MdcTrack.h.

28{ return *_theTrack; }
The documentation for this class was generated from the following files: