12 |
extern "C" { |
extern "C" { |
13 |
void dotrack_(int*, double*, double*, double*, double*, int*); |
void dotrack_(int*, double*, double*, double*, double*, int*); |
14 |
void dotrack2_(int*, double*, double*, double*, double*,double*, double*, double*,int*); |
void dotrack2_(int*, double*, double*, double*, double*,double*, double*, double*,int*); |
15 |
void mini2_(int*,int*,int*); |
void mini2_(int*,int*,int*); |
16 |
void guess_(); |
void guess_(); |
17 |
|
void gufld_(float*, float*); |
18 |
|
float risxeta2_(float *); |
19 |
|
float risxeta3_(float *); |
20 |
|
float risxeta4_(float *); |
21 |
|
float risyeta2_(float *); |
22 |
} |
} |
23 |
|
|
24 |
//-------------------------------------- |
//-------------------------------------- |
52 |
ayv[ip] = 0; |
ayv[ip] = 0; |
53 |
dedx_x[ip] = 0; |
dedx_x[ip] = 0; |
54 |
dedx_y[ip] = 0; |
dedx_y[ip] = 0; |
55 |
|
multmaxx[ip] = 0; |
56 |
|
multmaxy[ip] = 0; |
57 |
|
seedx[ip] = 0; |
58 |
|
seedy[ip] = 0; |
59 |
|
xpu[ip] = 0; |
60 |
|
ypu[ip] = 0; |
61 |
|
|
62 |
}; |
}; |
63 |
// clx = 0; |
|
64 |
// cly = 0; |
// TrkParams::SetTrackingMode(); |
65 |
// clx = new TRefArray(6,0); //forse causa memory leak??? |
// TrkParams::SetPrecisionFactor(); |
66 |
// cly = new TRefArray(6,0); //forse causa memory leak??? |
// TrkParams::SetStepMin(); |
67 |
// clx = TRefArray(6,0); |
TrkParams::SetMiniDefault(); |
|
// cly = TRefArray(6,0); |
|
|
|
|
|
TrkParams::SetTrackingMode(); |
|
|
TrkParams::SetPrecisionFactor(); |
|
|
TrkParams::SetStepMin(); |
|
68 |
TrkParams::SetPFA(); |
TrkParams::SetPFA(); |
69 |
|
|
70 |
|
int ngf = TrkParams::nGF; |
71 |
|
for(int i=0; i<ngf; i++){ |
72 |
|
xGF[i] = 0.; |
73 |
|
yGF[i] = 0.; |
74 |
|
} |
75 |
|
|
76 |
|
|
77 |
}; |
}; |
78 |
//-------------------------------------- |
//-------------------------------------- |
79 |
// |
// |
105 |
ayv[ip] = t.ayv[ip]; |
ayv[ip] = t.ayv[ip]; |
106 |
dedx_x[ip] = t.dedx_x[ip]; |
dedx_x[ip] = t.dedx_x[ip]; |
107 |
dedx_y[ip] = t.dedx_y[ip]; |
dedx_y[ip] = t.dedx_y[ip]; |
108 |
|
multmaxx[ip] = t.multmaxx[ip]; |
109 |
|
multmaxy[ip] = t.multmaxy[ip]; |
110 |
|
seedx[ip] = t.seedx[ip]; |
111 |
|
seedy[ip] = t.seedy[ip]; |
112 |
|
xpu[ip] = t.xpu[ip]; |
113 |
|
ypu[ip] = t.ypu[ip]; |
114 |
}; |
}; |
115 |
// clx = 0; |
|
116 |
// cly = 0; |
// TrkParams::SetTrackingMode(); |
117 |
// if(t.clx)clx = new TRefArray(*(t.clx)); |
// TrkParams::SetPrecisionFactor(); |
118 |
// if(t.cly)cly = new TRefArray(*(t.cly)); |
// TrkParams::SetStepMin(); |
119 |
// clx = TRefArray(t.clx); |
TrkParams::SetMiniDefault(); |
|
// cly = TRefArray(t.cly); |
|
|
|
|
|
TrkParams::SetTrackingMode(); |
|
|
TrkParams::SetPrecisionFactor(); |
|
|
TrkParams::SetStepMin(); |
|
120 |
TrkParams::SetPFA(); |
TrkParams::SetPFA(); |
121 |
|
|
122 |
|
int ngf = TrkParams::nGF; |
123 |
|
for(int i=0; i<ngf; i++){ |
124 |
|
xGF[i] = t.xGF[i]; |
125 |
|
yGF[i] = t.yGF[i]; |
126 |
|
} |
127 |
}; |
}; |
128 |
//-------------------------------------- |
//-------------------------------------- |
129 |
// |
// |
156 |
t.ayv[ip] = ayv[ip]; |
t.ayv[ip] = ayv[ip]; |
157 |
t.dedx_x[ip] = dedx_x[ip]; |
t.dedx_x[ip] = dedx_x[ip]; |
158 |
t.dedx_y[ip] = dedx_y[ip]; |
t.dedx_y[ip] = dedx_y[ip]; |
159 |
|
t.multmaxx[ip] = multmaxx[ip]; |
160 |
|
t.multmaxy[ip] = multmaxy[ip]; |
161 |
|
t.seedx[ip] = seedx[ip]; |
162 |
|
t.seedy[ip] = seedy[ip]; |
163 |
|
t.xpu[ip] = xpu[ip]; |
164 |
|
t.ypu[ip] = ypu[ip]; |
165 |
|
|
166 |
}; |
}; |
167 |
|
int ngf = TrkParams::nGF; |
168 |
|
for(int i=0; i<ngf; i++){ |
169 |
|
t.xGF[i] = xGF[i]; |
170 |
|
t.yGF[i] = yGF[i]; |
171 |
|
} |
172 |
|
|
|
// t.clx = TRefArray(clx); |
|
|
// t.cly = TRefArray(cly); |
|
173 |
|
|
174 |
}; |
}; |
175 |
//-------------------------------------- |
//-------------------------------------- |
182 |
* @param t pointer to an object of the class Trajectory, |
* @param t pointer to an object of the class Trajectory, |
183 |
* which z coordinates should be previously initialized by calling the proper constructor ( Trajectory::Trajectory(int n, float* zin) ). |
* which z coordinates should be previously initialized by calling the proper constructor ( Trajectory::Trajectory(int n, float* zin) ). |
184 |
* @return error flag. |
* @return error flag. |
185 |
|
* |
186 |
|
* >>> OBSOLETE !!! use TrkTrack::DoTrack2(Trajectory* t) instead |
187 |
|
* |
188 |
*/ |
*/ |
189 |
int TrkTrack::DoTrack(Trajectory* t){ |
int TrkTrack::DoTrack(Trajectory* t){ |
190 |
|
|
191 |
|
cout << " int TrkTrack::DoTrack(Trajectory* t) --->> OBSOLETE !!! "<<endl; |
192 |
|
cout << " use int TrkTrack::DoTrack2(Trajectory* t)"<<endl; |
193 |
|
|
194 |
double *dxout = new double[t->npoint]; |
double *dxout = new double[t->npoint]; |
195 |
double *dyout = new double[t->npoint]; |
double *dyout = new double[t->npoint]; |
196 |
double *dzin = new double[t->npoint]; |
double *dzin = new double[t->npoint]; |
322 |
return dedx; |
return dedx; |
323 |
}; |
}; |
324 |
/** |
/** |
325 |
* Returns 1 if the cluster on a tracker view includes bad strips. |
* Returns 1 if the cluster on a tracker view includes bad strips |
326 |
|
* (at least one bad strip among the four strip used by p.f.a.) |
327 |
* @param ip plane (0-5) |
* @param ip plane (0-5) |
328 |
* @param iv view (0=x 1=y) |
* @param iv view (0=x 1=y) |
329 |
*/ |
*/ |
396 |
} |
} |
397 |
return (last_plane-first_plane+1); |
return (last_plane-first_plane+1); |
398 |
} |
} |
399 |
|
/** |
400 |
|
* Returns the reduced chi-square of track x-projection |
401 |
|
*/ |
402 |
|
Float_t TrkTrack::GetChi2X(){ |
403 |
|
float chiq=0; |
404 |
|
for(int ip=0; ip<6; ip++)if(XGood(ip))chiq+= pow((xv[ip]-xm[ip])/resx[ip],2.); |
405 |
|
if(GetNX()>3)chiq=chiq/(GetNX()-3); |
406 |
|
else chiq=0; |
407 |
|
if(chiq==0)cout << " Float_t TrkTrack::GetChi2X() -- WARNING -- value not defined "<<chiq<<endl; |
408 |
|
return chiq; |
409 |
|
} |
410 |
|
/** |
411 |
|
* Returns the reduced chi-square of track y-projection |
412 |
|
*/ |
413 |
|
Float_t TrkTrack::GetChi2Y(){ |
414 |
|
float chiq=0; |
415 |
|
for(int ip=0; ip<6; ip++)if(YGood(ip))chiq+= pow((yv[ip]-ym[ip])/resy[ip],2.); |
416 |
|
if(GetNY()>2)chiq=chiq/(GetNY()-2); |
417 |
|
else chiq=0; |
418 |
|
if(chiq==0)cout << " Float_t TrkTrack::GetChi2Y() -- WARNING -- value not defined "<<chiq<<endl; |
419 |
|
return chiq; |
420 |
|
} |
421 |
|
/** |
422 |
|
* Returns the logarythm of the likeliwood-function of track x-projection |
423 |
|
*/ |
424 |
|
Float_t TrkTrack::GetLnLX(){ |
425 |
|
float lnl=0; |
426 |
|
for(int ip=0; ip<6; ip++) |
427 |
|
if( XGood(ip) && tailx[ip]!=0 ) |
428 |
|
lnl += (tailx[ip]+1.) * log( (tailx[ip]*pow(resx[ip],2.) + pow(xv[ip]-xm[ip],2.)) / (tailx[ip]*pow(resx[ip],2)) ); |
429 |
|
if(GetNX()>3)lnl=lnl/(GetNX()-3); |
430 |
|
else lnl=0; |
431 |
|
if(lnl==0){ |
432 |
|
cout << " Float_t TrkTrack::GetLnLX() -- WARNING -- value not defined "<<lnl<<endl; |
433 |
|
Dump(); |
434 |
|
} |
435 |
|
return lnl; |
436 |
|
|
437 |
|
} |
438 |
|
/** |
439 |
|
* Returns the logarythm of the likeliwood-function of track y-projection |
440 |
|
*/ |
441 |
|
Float_t TrkTrack::GetLnLY(){ |
442 |
|
float lnl=0; |
443 |
|
for(int ip=0; ip<6; ip++) |
444 |
|
if( YGood(ip) && taily[ip]!=0 ) |
445 |
|
lnl += (taily[ip]+1.) * log( (taily[ip]*pow(resy[ip],2.) + pow(yv[ip]-ym[ip],2.)) / (taily[ip]*pow(resy[ip],2)) ); |
446 |
|
if(GetNY()>2)lnl=lnl/(GetNY()-2); |
447 |
|
else lnl=0; |
448 |
|
if(lnl==0){ |
449 |
|
cout << " Float_t TrkTrack::GetLnLY() -- WARNING -- value not defined "<<lnl<<endl; |
450 |
|
Dump(); |
451 |
|
} |
452 |
|
return lnl; |
453 |
|
|
454 |
|
} |
455 |
|
/** |
456 |
|
* Returns the effective angle, relative to the sensor, on each plane. |
457 |
|
* @param ip plane (0-5) |
458 |
|
* @param iv view (0=x 1=y) |
459 |
|
*/ |
460 |
|
Float_t TrkTrack::GetEffectiveAngle(int ip, int iv){ |
461 |
|
|
462 |
|
if(ip<0 || ip>5){ |
463 |
|
cout << "Float_t TrkTrack::GetEffectiveAngle(int "<<ip<<", int "<<iv<<") ==> wrong input"<<endl; |
464 |
|
return 0.; |
465 |
|
} |
466 |
|
|
467 |
|
float v[3]={xv[ip],yv[ip],zv[ip]}; |
468 |
|
//----------------------------------------- |
469 |
|
// effective angle (relative to the sensor) |
470 |
|
//----------------------------------------- |
471 |
|
float axv_geo = axv[ip]; |
472 |
|
float muhall_h = 297.61; //cm**2/Vs |
473 |
|
float BY = TrkParams::GetBY(v); |
474 |
|
float axv_eff = 0; |
475 |
|
if(ip==5) axv_geo = -1*axv_geo; |
476 |
|
if(ip==5) BY = -1*BY; |
477 |
|
axv_eff = 180.*atan( tan(axv_geo*acos(-1.)/180.) + muhall_h * BY * 0.0001)/acos(-1.); |
478 |
|
//----------------------------------------- |
479 |
|
// effective angle (relative to the sensor) |
480 |
|
//----------------------------------------- |
481 |
|
float ayv_geo = ayv[ip]; |
482 |
|
float muhall_e = 1258.18; //cm**2/Vs |
483 |
|
float BX = TrkParams::GetBX(v); |
484 |
|
float ayv_eff = 0; |
485 |
|
ayv_eff = 180.*atan( tan(ayv_geo*acos(-1.)/180.) + muhall_e * BX * 0.0001)/acos(-1.); |
486 |
|
|
487 |
|
if (iv==0)return axv_eff; |
488 |
|
else if(iv==1)return ayv_eff; |
489 |
|
else{ |
490 |
|
cout << "Float_t TrkTrack::GetEffectiveAngle(int "<<ip<<", int "<<iv<<") ==> wrong input"<<endl; |
491 |
|
return 0.; |
492 |
|
} |
493 |
|
|
494 |
|
}; |
495 |
|
|
496 |
//-------------------------------------- |
//-------------------------------------- |
497 |
// |
// |
498 |
// |
// |
514 |
cout << endl << "zv : "; for(int i=0; i<6; i++)cout << zv[i] << " "; |
cout << endl << "zv : "; for(int i=0; i<6; i++)cout << zv[i] << " "; |
515 |
cout << endl << "resx : "; for(int i=0; i<6; i++)cout << resx[i] << " "; |
cout << endl << "resx : "; for(int i=0; i<6; i++)cout << resx[i] << " "; |
516 |
cout << endl << "resy : "; for(int i=0; i<6; i++)cout << resy[i] << " "; |
cout << endl << "resy : "; for(int i=0; i<6; i++)cout << resy[i] << " "; |
517 |
|
cout << endl << "tailx : "; for(int i=0; i<6; i++)cout << tailx[i] << " "; |
518 |
|
cout << endl << "taily : "; for(int i=0; i<6; i++)cout << taily[i] << " "; |
519 |
cout << endl << "coval : "; for(int i=0; i<5; i++)cout << coval[0][i]<<" "; |
cout << endl << "coval : "; for(int i=0; i<5; i++)cout << coval[0][i]<<" "; |
520 |
cout << endl << " "; for(int i=0; i<5; i++)cout << coval[1][i]<<" "; |
cout << endl << " "; for(int i=0; i<5; i++)cout << coval[1][i]<<" "; |
521 |
cout << endl << " "; for(int i=0; i<5; i++)cout << coval[2][i]<<" "; |
cout << endl << " "; for(int i=0; i<5; i++)cout << coval[2][i]<<" "; |
523 |
cout << endl << " "; for(int i=0; i<5; i++)cout << coval[4][i]<<" "; |
cout << endl << " "; for(int i=0; i<5; i++)cout << coval[4][i]<<" "; |
524 |
cout << endl << "dedx_x : "; for(int i=0; i<6; i++)cout << dedx_x[i] << " "; |
cout << endl << "dedx_x : "; for(int i=0; i<6; i++)cout << dedx_x[i] << " "; |
525 |
cout << endl << "dedx_y : "; for(int i=0; i<6; i++)cout << dedx_y[i] << " "; |
cout << endl << "dedx_y : "; for(int i=0; i<6; i++)cout << dedx_y[i] << " "; |
526 |
|
cout << endl << "maxs x : "; for(int i=0; i<6; i++)cout << GetClusterX_MaxStrip(i) << " "; |
527 |
|
cout << endl << "maxs y : "; for(int i=0; i<6; i++)cout << GetClusterY_MaxStrip(i) << " "; |
528 |
|
cout << endl << "mult x : "; for(int i=0; i<6; i++)cout << GetClusterX_Multiplicity(i) << " "; |
529 |
|
cout << endl << "mult y : "; for(int i=0; i<6; i++)cout << GetClusterY_Multiplicity(i) << " "; |
530 |
|
cout << endl << "seed x : "; for(int i=0; i<6; i++)cout << GetClusterX_Seed(i) << " "; |
531 |
|
cout << endl << "seed y : "; for(int i=0; i<6; i++)cout << GetClusterY_Seed(i) << " "; |
532 |
|
cout << endl << "xpu : "; for(int i=0; i<6; i++)cout << xpu[i] << " "; |
533 |
|
cout << endl << "ypu : "; for(int i=0; i<6; i++)cout << ypu[i] << " "; |
534 |
|
|
535 |
cout << endl; |
cout << endl; |
536 |
} |
} |
537 |
/** |
/** |
550 |
for(int i=0; i<6; i++) resy[i]=*ry++; |
for(int i=0; i<6; i++) resy[i]=*ry++; |
551 |
} |
} |
552 |
/** |
/** |
553 |
* Set the TrkTrack good measurement. |
* Set the TrkTrack tails position resolution |
554 |
|
*/ |
555 |
|
void TrkTrack::SetTail(double *tx, double *ty, double factor){ |
556 |
|
for(int i=0; i<6; i++) tailx[i]=factor*(*tx++); |
557 |
|
for(int i=0; i<6; i++) taily[i]=factor*(*ty++); |
558 |
|
} |
559 |
|
/** |
560 |
|
* Set the TrkTrack Student parameter (resx,resy,tailx,taily) |
561 |
|
* from previous gausian fit |
562 |
|
*@param flag =0 standard, =1 with noise correction |
563 |
|
*/ |
564 |
|
void TrkTrack::SetStudentParam(int flag){ |
565 |
|
float sx[11]={0.000128242, |
566 |
|
0.000136942, |
567 |
|
0.000162718, |
568 |
|
0.000202644, |
569 |
|
0.00025597, |
570 |
|
0.000317456, |
571 |
|
0.000349048, |
572 |
|
0.000384638, |
573 |
|
0.000457295, |
574 |
|
0.000512319, |
575 |
|
0.000538573}; |
576 |
|
float tx[11]={1.79402, |
577 |
|
2.04876, |
578 |
|
2.88376, |
579 |
|
3.3, |
580 |
|
3.14084, |
581 |
|
4.07686, |
582 |
|
4.44736, |
583 |
|
3.5179, |
584 |
|
3.38697, |
585 |
|
3.45739, |
586 |
|
3.18627}; |
587 |
|
float sy[11]={0.000483075, |
588 |
|
0.000466925, |
589 |
|
0.000431658, |
590 |
|
0.000428317, |
591 |
|
0.000433854, |
592 |
|
0.000444044, |
593 |
|
0.000482098, |
594 |
|
0.000537579, |
595 |
|
0.000636279, |
596 |
|
0.000741998, |
597 |
|
0.000864261}; |
598 |
|
float ty[11]={0.997032, |
599 |
|
1.11147, |
600 |
|
1.18526, |
601 |
|
1.61404, |
602 |
|
2.21908, |
603 |
|
3.08959, |
604 |
|
4.48833, |
605 |
|
4.42687, |
606 |
|
4.65253, |
607 |
|
4.52043, |
608 |
|
4.29926}; |
609 |
|
int index; |
610 |
|
float fact; |
611 |
|
for(int i=0; i<6; i++) { |
612 |
|
index = int((fabs(axv[i])+1.)/2.); |
613 |
|
if(index>10) index=10; |
614 |
|
tailx[i]=tx[index]; |
615 |
|
if(flag==1) { |
616 |
|
if(fabs(axv[i])<=10.) fact = resx[i]/risxeta2_(&(axv[i])); |
617 |
|
if(fabs(axv[i])>10.&&fabs(axv[i])<=15.) fact = resx[i]/risxeta3_(&(axv[i])); |
618 |
|
if(fabs(axv[i])>15.) fact = resx[i]/risxeta4_(&(axv[i])); |
619 |
|
} else fact = 1.; |
620 |
|
resx[i] = sx[index]*fact; |
621 |
|
} |
622 |
|
for(int i=0; i<6; i++) { |
623 |
|
index = int((fabs(ayv[i])+1.)/2.); |
624 |
|
if(index>10) index=10; |
625 |
|
taily[i]=ty[index]; |
626 |
|
if(flag==1) fact = resy[i]/risyeta2_(&(ayv[i])); |
627 |
|
else fact = 1.; |
628 |
|
resy[i] = sy[index]*fact; |
629 |
|
} |
630 |
|
} |
631 |
|
/** |
632 |
|
* Set the TrkTrack good measurement |
633 |
*/ |
*/ |
634 |
void TrkTrack::SetGood(int *xg, int *yg){ |
void TrkTrack::SetGood(int *xg, int *yg){ |
635 |
|
|
647 |
// path_.error = 0; |
// path_.error = 0; |
648 |
// readb_(); |
// readb_(); |
649 |
|
|
650 |
TrkParams::SetTrackingMode(); |
// TrkParams::SetTrackingMode(); |
651 |
TrkParams::SetPrecisionFactor(); |
// TrkParams::SetPrecisionFactor(); |
652 |
TrkParams::SetStepMin(); |
// TrkParams::SetStepMin(); |
653 |
|
TrkParams::SetMiniDefault(); |
654 |
|
|
655 |
TrkParams::Set(path,1); |
TrkParams::Set(path,1); |
656 |
TrkParams::Load(1); |
TrkParams::Load(1); |
692 |
|
|
693 |
track.resx[i]=resx[i]; |
track.resx[i]=resx[i]; |
694 |
track.resy[i]=resy[i]; |
track.resy[i]=resy[i]; |
695 |
|
track.tailx[i]=tailx[i]; |
696 |
|
track.taily[i]=taily[i]; |
697 |
} |
} |
698 |
|
|
699 |
for(int i=0; i<5; i++) track.al[i]=al[i]; |
for(int i=0; i<5; i++) track.al[i]=al[i]; |
749 |
* |
* |
750 |
* @see TrkParams::SetPFA(int) |
* @see TrkParams::SetPFA(int) |
751 |
*/ |
*/ |
752 |
void TrkTrack::EvaluateClusterPositions(){ |
Bool_t TrkTrack::EvaluateClusterPositions(){ |
753 |
|
|
754 |
// cout << "void TrkTrack::GetClusterPositions() "<<endl; |
// cout << "void TrkTrack::GetClusterositions() "<<endl; |
755 |
|
|
756 |
TrkParams::Load( ); |
TrkParams::Load( ); |
757 |
if( !TrkParams::IsLoaded() )return; |
if( !TrkParams::IsLoaded() )return false; |
758 |
|
|
759 |
for(int ip=0; ip<6; ip++){ |
for(int ip=0; ip<6; ip++){ |
760 |
// cout << ip<<" ** "<<xm[ip]<<" / "<<ym[ip]<<endl;; |
// cout << ip<<" ** "<<xm[ip]<<" / "<<ym[ip]<<endl;; |
773 |
float bfy = 10*TrkParams::GetBY(v);//Tesla |
float bfy = 10*TrkParams::GetBY(v);//Tesla |
774 |
int ipp=ip+1; |
int ipp=ip+1; |
775 |
xyzpam_(&ipp,&icx,&icy,&ladder,&sensor,&ax,&ay,&bfx,&bfy); |
xyzpam_(&ipp,&icx,&icy,&ladder,&sensor,&ax,&ay,&bfx,&bfy); |
776 |
|
if(icx<0 || icy<0)return false; |
777 |
} |
} |
778 |
|
return true; |
779 |
} |
} |
780 |
/** |
/** |
781 |
* \brief Tracking method. It calls F77 mini routine. |
* \brief Tracking method. It calls F77 mini routine. |
799 |
* @see EvaluateClusterPositions() |
* @see EvaluateClusterPositions() |
800 |
* |
* |
801 |
* The fitting procedure can be varied by changing the tracking mode, |
* The fitting procedure can be varied by changing the tracking mode, |
802 |
* the fit-precision factor and the minimum number of step. |
* the fit-precision factor, the minimum number of step, etc. |
803 |
* @see SetTrackingMode(int) |
* @see SetTrackingMode(int) |
804 |
* @see SetPrecisionFactor(double) |
* @see SetPrecisionFactor(double) |
805 |
* @see SetStepMin(int) |
* @see SetStepMin(int) |
806 |
|
* @see SetDeltaB(int,double) |
807 |
*/ |
*/ |
808 |
void TrkTrack::Fit(double pfixed, int& fail, int iprint, int froml1){ |
void TrkTrack::Fit(double pfixed, int& fail, int iprint, int froml1){ |
809 |
|
|
814 |
|
|
815 |
extern cMini2track track_; |
extern cMini2track track_; |
816 |
fail = 0; |
fail = 0; |
817 |
|
|
818 |
FillMiniStruct(track_); |
FillMiniStruct(track_); |
819 |
|
|
820 |
if(froml1!=0)EvaluateClusterPositions(); |
if(froml1!=0){ |
821 |
|
if( !EvaluateClusterPositions() ){ |
822 |
|
cout << "void TrkTrack::Fit("<<pfixed<<","<<fail<<","<<iprint<<","<<froml1<<") --- ERROR evaluating cluster positions "<<endl; |
823 |
|
FillMiniStruct(track_) ; |
824 |
|
fail = 1; |
825 |
|
return; |
826 |
|
} |
827 |
|
}else{ |
828 |
|
FillMiniStruct(track_); |
829 |
|
} |
830 |
|
|
831 |
// if fit variables have been reset, evaluate the initial guess |
// if fit variables have been reset, evaluate the initial guess |
832 |
if(al[0]==-9999.&&al[1]==-9999.&&al[2]==-9999.&&al[3]==-9999.&&al[4]==-9999.)guess_(); |
if(al[0]==-9999.&&al[1]==-9999.&&al[2]==-9999.&&al[3]==-9999.&&al[4]==-9999.)guess_(); |
899 |
track_.fact = fact; |
track_.fact = fact; |
900 |
} |
} |
901 |
/** |
/** |
902 |
* Set the factor scale for tracking precision |
* Set the minimum number of steps for tracking precision |
903 |
*/ |
*/ |
904 |
void TrkTrack::SetStepMin(int istepmin){ |
void TrkTrack::SetStepMin(int istepmin){ |
905 |
extern cMini2track track_; |
extern cMini2track track_; |
906 |
track_.istepmin = istepmin; |
track_.istepmin = istepmin; |
907 |
} |
} |
908 |
|
/** |
909 |
|
* Set deltaB parameters (id=0,1). By default they are set to zero. |
910 |
|
*/ |
911 |
|
void TrkTrack::SetDeltaB(int id, double db){ |
912 |
|
if(id!=0 && id!=1)cout << "void TrkTrack::SetDeltaB(int id,double db) -- wrong input parameters: "<<id<<" "<<db<<endl; |
913 |
|
TrkParams::SetDeltaB(id,db); |
914 |
|
} |
915 |
|
|
916 |
|
/** |
917 |
|
* Returns true if the track is inside the magnet cavity. |
918 |
|
* @param toll Tolerance around the nominal volume (toll>0 define an inner fiducial volume) |
919 |
|
*/ |
920 |
|
Bool_t TrkTrack::IsInsideCavity(float toll){ |
921 |
|
|
922 |
|
// float xmagntop, ymagntop, xmagnbottom, ymagnbottom; |
923 |
|
// xmagntop = xv[0] + (ZMAGNHIGH-zv[0])*tan(acos(-1.0)*axv[0]/180.); |
924 |
|
// ymagntop = yv[0] + (ZMAGNHIGH-zv[0])*tan(acos(-1.0)*ayv[0]/180.); |
925 |
|
// xmagnbottom = xv[5] + (ZMAGNLOW-zv[5])*tan(acos(-1.0)*axv[5]/180.); |
926 |
|
// ymagnbottom = yv[5] + (ZMAGNLOW-zv[5])*tan(acos(-1.0)*ayv[5]/180.); |
927 |
|
// if( xmagntop>XMAGNLOW && xmagntop<XMAGNHIGH && |
928 |
|
// ymagntop>YMAGNLOW && ymagntop<YMAGNHIGH && |
929 |
|
// xmagnbottom>XMAGNLOW && xmagnbottom<XMAGNHIGH && |
930 |
|
// ymagnbottom>YMAGNLOW && ymagnbottom<YMAGNHIGH ) return(true); |
931 |
|
// else return(false); |
932 |
|
|
933 |
|
int ngf = TrkParams::nGF; |
934 |
|
for(int i=0; i<ngf; i++){ |
935 |
|
// |
936 |
|
// cout << endl << TrkParams::GF_element[i]; |
937 |
|
if( |
938 |
|
TrkParams::GF_element[i].CompareTo("CUF") && |
939 |
|
TrkParams::GF_element[i].CompareTo("T2") && |
940 |
|
TrkParams::GF_element[i].CompareTo("T3") && |
941 |
|
TrkParams::GF_element[i].CompareTo("T4") && |
942 |
|
TrkParams::GF_element[i].CompareTo("T5") && |
943 |
|
TrkParams::GF_element[i].CompareTo("CLF") && |
944 |
|
true)continue; |
945 |
|
// apply condition only within the cavity |
946 |
|
// cout << " -- "<<xGF[i]<<" "<<yGF[i]; |
947 |
|
if( |
948 |
|
xGF[i] <= TrkParams::xGF_min[i] + toll || |
949 |
|
xGF[i] >= TrkParams::xGF_max[i] - toll || |
950 |
|
yGF[i] <= TrkParams::yGF_min[i] + toll || |
951 |
|
yGF[i] >= TrkParams::yGF_max[i] - toll || |
952 |
|
false){ |
953 |
|
|
954 |
|
return false; |
955 |
|
} |
956 |
|
} |
957 |
|
return true; |
958 |
|
|
959 |
|
|
960 |
|
} |
961 |
|
/** |
962 |
|
* Returns true if the track is inside the nominal acceptance, which is defined |
963 |
|
* by the intersection among magnet cavity, silicon-plane sensitive area and |
964 |
|
* ToF sensitive area (nominal values from the official document used to |
965 |
|
* calculate the geometrical factor) |
966 |
|
*/ |
967 |
|
Bool_t TrkTrack::IsInsideAcceptance(){ |
968 |
|
|
969 |
|
int ngf = TrkParams::nGF; |
970 |
|
for(int i=0; i<ngf; i++){ |
971 |
|
if( |
972 |
|
xGF[i] <= TrkParams::xGF_min[i] || |
973 |
|
xGF[i] >= TrkParams::xGF_max[i] || |
974 |
|
yGF[i] <= TrkParams::yGF_min[i] || |
975 |
|
yGF[i] >= TrkParams::yGF_max[i] || |
976 |
|
false)return false; |
977 |
|
} |
978 |
|
return true; |
979 |
|
|
980 |
|
} |
981 |
/** |
/** |
982 |
* Method to retrieve ID (0,1,...) of x-cluster (if any) associated to this track. |
* Method to retrieve ID (0,1,...) of x-cluster (if any) associated to this track. |
983 |
* If no cluster is associated, ID=-1. |
* If no cluster is associated, ID=-1. |
994 |
Int_t TrkTrack::GetClusterY_ID(int ip){ |
Int_t TrkTrack::GetClusterY_ID(int ip){ |
995 |
return ((Int_t)fabs(ygood[ip]))%10000000-1; |
return ((Int_t)fabs(ygood[ip]))%10000000-1; |
996 |
}; |
}; |
997 |
|
|
998 |
/** |
/** |
999 |
* Method to retrieve the ladder (0-4, increasing x) traversed by the track on this plane. |
* Method to retrieve the ladder (0-4, increasing x) traversed by the track on this plane. |
1000 |
* If no ladder is traversed (dead area) the metod retuns -1. |
* If no ladder is traversed (dead area) the metod retuns -1. |
1043 |
if(bad)ygood[ip]=-ygood[ip]; |
if(bad)ygood[ip]=-ygood[ip]; |
1044 |
}; |
}; |
1045 |
|
|
1046 |
|
/** |
1047 |
|
* \brief Average X |
1048 |
|
* Average value of <xv>, evaluated from the first to the last hit x view. |
1049 |
|
*/ |
1050 |
|
Float_t TrkTrack::GetXav(){ |
1051 |
|
|
1052 |
|
int first_plane = -1; |
1053 |
|
int last_plane = -1; |
1054 |
|
for(Int_t ip=0; ip<6; ip++){ |
1055 |
|
if( XGood(ip) && first_plane == -1 )first_plane = ip; |
1056 |
|
if( XGood(ip) && first_plane != -1 )last_plane = ip; |
1057 |
|
} |
1058 |
|
if( first_plane == -1 || last_plane == -1){ |
1059 |
|
return -100; |
1060 |
|
} |
1061 |
|
if( last_plane-first_plane+1 ==0 )return -100; |
1062 |
|
|
1063 |
|
Float_t av = 0; |
1064 |
|
for(int ip=first_plane; ip<=last_plane; ip++)av+=xv[ip]; |
1065 |
|
|
1066 |
|
return (av/(last_plane-first_plane+1)); |
1067 |
|
} |
1068 |
|
/** |
1069 |
|
* \brief Average Y |
1070 |
|
* Average value of <yv>, evaluated from the first to the last hit x view. |
1071 |
|
*/ |
1072 |
|
Float_t TrkTrack::GetYav(){ |
1073 |
|
|
1074 |
|
int first_plane = -1; |
1075 |
|
int last_plane = -1; |
1076 |
|
for(Int_t ip=0; ip<6; ip++){ |
1077 |
|
if( XGood(ip) && first_plane == -1 )first_plane = ip; |
1078 |
|
if( XGood(ip) && first_plane != -1 )last_plane = ip; |
1079 |
|
} |
1080 |
|
if( first_plane == -1 || last_plane == -1){ |
1081 |
|
return -100; |
1082 |
|
} |
1083 |
|
if( last_plane-first_plane+1 ==0 )return -100; |
1084 |
|
|
1085 |
|
Float_t av = 0; |
1086 |
|
for(int ip=first_plane; ip<=last_plane; ip++)av+=yv[ip]; |
1087 |
|
|
1088 |
|
return (av/(last_plane-first_plane+1)); |
1089 |
|
} |
1090 |
|
/** |
1091 |
|
* \brief Average Z |
1092 |
|
* Average value of <zv>, evaluated from the first to the last hit x view. |
1093 |
|
*/ |
1094 |
|
Float_t TrkTrack::GetZav(){ |
1095 |
|
|
1096 |
|
int first_plane = -1; |
1097 |
|
int last_plane = -1; |
1098 |
|
for(Int_t ip=0; ip<6; ip++){ |
1099 |
|
if( XGood(ip) && first_plane == -1 )first_plane = ip; |
1100 |
|
if( XGood(ip) && first_plane != -1 )last_plane = ip; |
1101 |
|
} |
1102 |
|
if( first_plane == -1 || last_plane == -1){ |
1103 |
|
return -100; |
1104 |
|
} |
1105 |
|
if( last_plane-first_plane+1 ==0 )return -100; |
1106 |
|
|
1107 |
|
Float_t av = 0; |
1108 |
|
for(int ip=first_plane; ip<=last_plane; ip++)av+=zv[ip]; |
1109 |
|
|
1110 |
|
return (av/(last_plane-first_plane+1)); |
1111 |
|
} |
1112 |
|
|
1113 |
|
/** |
1114 |
|
* \brief Number of column traversed |
1115 |
|
*/ |
1116 |
|
Int_t TrkTrack::GetNColumns(){ |
1117 |
|
int sensors[] = {0,0,0,0,0,0}; |
1118 |
|
for(int ip=0; ip<6; ip++){ |
1119 |
|
int sensorid = GetLadder(ip)+3*GetSensor(ip); |
1120 |
|
if(XGood(ip)||YGood(ip)) |
1121 |
|
if(sensorid>=0 && sensorid<6)sensors[sensorid]=1; |
1122 |
|
} |
1123 |
|
int nsensors=0; |
1124 |
|
for(int is=0; is<6; is++)nsensors += sensors[is]; |
1125 |
|
return nsensors; |
1126 |
|
}; |
1127 |
|
/** |
1128 |
|
* \brief Give the maximum energy release |
1129 |
|
*/ |
1130 |
|
Float_t TrkTrack::GetDEDX_max(int ip, int iv){ |
1131 |
|
Float_t max=0; |
1132 |
|
int pfrom = 0; |
1133 |
|
int pto = 6; |
1134 |
|
int vfrom = 0; |
1135 |
|
int vto = 2; |
1136 |
|
if(ip>=0&&ip<6){ |
1137 |
|
pfrom = ip; |
1138 |
|
pto = ip+1; |
1139 |
|
} |
1140 |
|
if(iv>=0&&iv<2){ |
1141 |
|
vfrom = iv; |
1142 |
|
vto = iv+1; |
1143 |
|
} |
1144 |
|
for(int i=pfrom; i<pto; i++) |
1145 |
|
for(int j=0; j<vto; j++) |
1146 |
|
if(GetDEDX(i,j)>max)max=GetDEDX(i,j); |
1147 |
|
|
1148 |
|
return max; |
1149 |
|
|
1150 |
|
}; |
1151 |
|
|
1152 |
|
/** |
1153 |
|
* \brief Give the minimum energy release |
1154 |
|
*/ |
1155 |
|
Float_t TrkTrack::GetDEDX_min(int ip, int iv){ |
1156 |
|
Float_t min=100000000; |
1157 |
|
int pfrom = 0; |
1158 |
|
int pto = 6; |
1159 |
|
int vfrom = 0; |
1160 |
|
int vto = 2; |
1161 |
|
if(ip>=0&&ip<6){ |
1162 |
|
pfrom = ip; |
1163 |
|
pto = ip+1; |
1164 |
|
} |
1165 |
|
if(iv>=0&&iv<2){ |
1166 |
|
vfrom = iv; |
1167 |
|
vto = iv+1; |
1168 |
|
} |
1169 |
|
for(int i=pfrom; i<pto; i++) |
1170 |
|
for(int j=0; j<vto; j++) |
1171 |
|
if(GetDEDX(i,j)<min)min=GetDEDX(i,j); |
1172 |
|
|
1173 |
|
return min; |
1174 |
|
|
1175 |
|
}; |
1176 |
|
|
1177 |
|
/** |
1178 |
|
* \brief Give the maximum spatial residual release |
1179 |
|
*/ |
1180 |
|
Float_t TrkTrack::GetResidual_max(int ip, int iv){ |
1181 |
|
Float_t max=0; |
1182 |
|
int pfrom = 0; |
1183 |
|
int pto = 6; |
1184 |
|
int vfrom = 0; |
1185 |
|
int vto = 2; |
1186 |
|
if(ip>=0&&ip<6){ |
1187 |
|
pfrom = ip; |
1188 |
|
pto = ip+1; |
1189 |
|
} |
1190 |
|
if(iv>=0&&iv<2){ |
1191 |
|
vfrom = iv; |
1192 |
|
vto = iv+1; |
1193 |
|
} |
1194 |
|
for(int i=pfrom; i<pto; i++){ |
1195 |
|
for(int j=0; j<vto; j++){ |
1196 |
|
if(j==0 && XGood(i) && fabs(xm[i]-xv[i])>fabs(max))max=xv[i]-xm[i]; |
1197 |
|
if(j==1 && YGood(i) && fabs(ym[i]-yv[i])>fabs(max))max=yv[i]-ym[i]; |
1198 |
|
} |
1199 |
|
} |
1200 |
|
return max; |
1201 |
|
|
1202 |
|
}; |
1203 |
|
|
1204 |
|
|
1205 |
|
/** |
1206 |
|
* \brief Give the maximum multiplicity on the x view |
1207 |
|
*/ |
1208 |
|
Int_t TrkTrack::GetClusterX_Multiplicity_max(){ |
1209 |
|
int max=0; |
1210 |
|
for(int ip=0; ip<6; ip++) |
1211 |
|
if(GetClusterX_Multiplicity(ip)>max)max=GetClusterX_Multiplicity(ip); |
1212 |
|
return max; |
1213 |
|
}; |
1214 |
|
/** |
1215 |
|
* \brief Give the minimum multiplicity on the x view |
1216 |
|
*/ |
1217 |
|
Int_t TrkTrack::GetClusterX_Multiplicity_min(){ |
1218 |
|
int min=50; |
1219 |
|
for(int ip=0; ip<6; ip++) |
1220 |
|
if(GetClusterX_Multiplicity(ip)<min)min=GetClusterX_Multiplicity(ip); |
1221 |
|
return min; |
1222 |
|
}; |
1223 |
|
/** |
1224 |
|
* \brief Give the maximum multiplicity on the x view |
1225 |
|
*/ |
1226 |
|
Int_t TrkTrack::GetClusterY_Multiplicity_max(){ |
1227 |
|
int max=0; |
1228 |
|
for(int ip=0; ip<6; ip++) |
1229 |
|
if(GetClusterY_Multiplicity(ip)>max)max=GetClusterY_Multiplicity(ip); |
1230 |
|
return max; |
1231 |
|
}; |
1232 |
|
/** |
1233 |
|
* \brief Give the minimum multiplicity on the x view |
1234 |
|
*/ |
1235 |
|
Int_t TrkTrack::GetClusterY_Multiplicity_min(){ |
1236 |
|
int min=50; |
1237 |
|
for(int ip=0; ip<6; ip++) |
1238 |
|
if(GetClusterY_Multiplicity(ip)<min)min=GetClusterY_Multiplicity(ip); |
1239 |
|
return min; |
1240 |
|
}; |
1241 |
|
|
1242 |
|
/** |
1243 |
|
* \brief Give the minimum seed on the x view |
1244 |
|
*/ |
1245 |
|
Float_t TrkTrack::GetClusterX_Seed_min(){ |
1246 |
|
Float_t min=100000; |
1247 |
|
for(int ip=0; ip<6; ip++) |
1248 |
|
if(XGood(ip) && GetClusterX_Seed(ip)<min)min=GetClusterX_Seed(ip); |
1249 |
|
return min; |
1250 |
|
}; |
1251 |
|
/** |
1252 |
|
* \brief Give the minimum seed on the x view |
1253 |
|
*/ |
1254 |
|
Float_t TrkTrack::GetClusterY_Seed_min(){ |
1255 |
|
Float_t min=100000; |
1256 |
|
for(int ip=0; ip<6; ip++) |
1257 |
|
if(YGood(ip) && GetClusterY_Seed(ip)<min)min=GetClusterY_Seed(ip); |
1258 |
|
return min; |
1259 |
|
}; |
1260 |
|
|
1261 |
|
|
1262 |
//-------------------------------------- |
//-------------------------------------- |
1263 |
// |
// |
1264 |
// |
// |
1292 |
dedx_y[ip] = 0; |
dedx_y[ip] = 0; |
1293 |
|
|
1294 |
}; |
}; |
1295 |
|
int ngf = TrkParams::nGF; |
1296 |
|
for(int i=0; i<ngf; i++){ |
1297 |
|
xGF[i] = 0.; |
1298 |
|
yGF[i] = 0.; |
1299 |
|
} |
1300 |
// if(clx)clx->Clear(); |
// if(clx)clx->Clear(); |
1301 |
// if(cly)cly->Clear(); |
// if(cly)cly->Clear(); |
1302 |
// clx.Clear(); |
// clx.Clear(); |
1323 |
//-------------------------------------- |
//-------------------------------------- |
1324 |
TrkSinglet::TrkSinglet(){ |
TrkSinglet::TrkSinglet(){ |
1325 |
// cout << "TrkSinglet::TrkSinglet() " << GetUniqueID()<<endl; |
// cout << "TrkSinglet::TrkSinglet() " << GetUniqueID()<<endl; |
1326 |
plane = 0; |
// plane = 0; |
1327 |
coord[0] = 0; |
// coord[0] = 0; |
1328 |
coord[1] = 0; |
// coord[1] = 0; |
1329 |
sgnl = 0; |
// sgnl = 0; |
1330 |
|
// multmax = 0; |
1331 |
// cls = 0; |
// cls = 0; |
1332 |
|
Clear(); |
1333 |
}; |
}; |
1334 |
//-------------------------------------- |
//-------------------------------------- |
1335 |
// |
// |
1341 |
coord[0] = s.coord[0]; |
coord[0] = s.coord[0]; |
1342 |
coord[1] = s.coord[1]; |
coord[1] = s.coord[1]; |
1343 |
sgnl = s.sgnl; |
sgnl = s.sgnl; |
1344 |
|
multmax = s.multmax; |
1345 |
// cls = 0;//<<<<pointer |
// cls = 0;//<<<<pointer |
1346 |
// cls = TRef(s.cls); |
// cls = TRef(s.cls); |
1347 |
}; |
}; |
1352 |
void TrkSinglet::Dump(){ |
void TrkSinglet::Dump(){ |
1353 |
int i=0; |
int i=0; |
1354 |
cout << endl << "========== Singlet " ; |
cout << endl << "========== Singlet " ; |
1355 |
cout << endl << "plane : " << plane; |
cout << endl << "plane : " << plane; |
1356 |
cout << endl << "coord[2] : "; while( i<2 && cout << coord[i] << " ") i++; |
cout << endl << "coord[2] : "; while( i<2 && cout << coord[i] << " ") i++; |
1357 |
cout << endl << "sgnl : " << sgnl; |
cout << endl << "sgnl : " << sgnl; |
1358 |
|
cout << endl << "max.strip : "; |
1359 |
|
cout << endl << "multiplicity : "; |
1360 |
} |
} |
1361 |
//-------------------------------------- |
//-------------------------------------- |
1362 |
// |
// |
1369 |
coord[0]=-999; |
coord[0]=-999; |
1370 |
coord[1]=-999; |
coord[1]=-999; |
1371 |
sgnl=0; |
sgnl=0; |
1372 |
|
multmax = 0; |
1373 |
|
|
1374 |
} |
} |
1375 |
//-------------------------------------- |
//-------------------------------------- |
1405 |
|
|
1406 |
// |
// |
1407 |
cout << endl << endl << "=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-"; |
cout << endl << endl << "=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-"; |
1408 |
cout << endl << "good : "; for(int i=0; i<12; i++) cout << good[i]<<" "; |
cout << endl << "good : "; for(int i=0; i<12; i++) cout << hex <<" 0x"<< good[i]<<dec; |
1409 |
cout << endl << "ntrk() : " << this->ntrk() ; |
cout << endl << "ntrk() : " << ntrk() ; |
1410 |
cout << endl << "nclsx() : " << this->nclsx(); |
cout << endl << "nclsx() : " << nclsx(); |
1411 |
cout << endl << "nclsy() : " << this->nclsy(); |
cout << endl << "nclsy() : " << nclsy(); |
1412 |
if(Track){ |
if(Track){ |
1413 |
TClonesArray &t = *Track; |
TClonesArray &t = *Track; |
1414 |
for(int i=0; i<ntrk(); i++) ((TrkTrack *)t[i])->Dump(); |
for(int i=0; i<ntrk(); i++) ((TrkTrack *)t[i])->Dump(); |
1415 |
} |
} |
1416 |
if(SingletX){ |
// if(SingletX){ |
1417 |
TClonesArray &sx = *SingletX; |
// TClonesArray &sx = *SingletX; |
1418 |
for(int i=0; i<nclsx(); i++) ((TrkSinglet *)sx[i])->Dump(); |
// for(int i=0; i<nclsx(); i++) ((TrkSinglet *)sx[i])->Dump(); |
1419 |
} |
// } |
1420 |
if(SingletY){ |
// if(SingletY){ |
1421 |
TClonesArray &sy = *SingletY; |
// TClonesArray &sy = *SingletY; |
1422 |
for(int i=0; i<nclsy(); i++) ((TrkSinglet *)sy[i])->Dump(); |
// for(int i=0; i<nclsy(); i++) ((TrkSinglet *)sy[i])->Dump(); |
1423 |
} |
// } |
1424 |
|
cout << endl; |
1425 |
} |
} |
1426 |
|
/** |
1427 |
|
* \brief Dump processing status |
1428 |
|
*/ |
1429 |
|
void TrkLevel2::StatusDump(int view){ |
1430 |
|
cout << "DSP n. "<<view+1<<" status: "<<hex<<good[view]<<endl; |
1431 |
|
}; |
1432 |
|
/** |
1433 |
|
* \brief Check event status |
1434 |
|
* |
1435 |
|
* Check the event status, according to a flag-mask given as input. |
1436 |
|
* Return true if the view passes the check. |
1437 |
|
* |
1438 |
|
* @param view View number (0-11) |
1439 |
|
* @param flagmask Mask of flags to check (eg. flagmask=0x111 no missing packet, |
1440 |
|
* no crc error, no software alarm) |
1441 |
|
* |
1442 |
|
* @see TrkLevel2 class definition to know how the status flag is defined |
1443 |
|
* |
1444 |
|
*/ |
1445 |
|
Bool_t TrkLevel2::StatusCheck(int view, int flagmask){ |
1446 |
|
|
1447 |
|
if( view<0 || view >= 12)return false; |
1448 |
|
return !(good[view]&flagmask); |
1449 |
|
|
1450 |
|
}; |
1451 |
|
|
1452 |
|
|
1453 |
//-------------------------------------- |
//-------------------------------------- |
1454 |
// |
// |
1455 |
// |
// |
1480 |
* forced (see TrkLevel2::GetVKMask(int,int)) or |
* forced (see TrkLevel2::GetVKMask(int,int)) or |
1481 |
* for this event only (TrkLevel2::GetVKFlag(int,int)). |
* for this event only (TrkLevel2::GetVKFlag(int,int)). |
1482 |
* @param iv Tracker view (0-11) |
* @param iv Tracker view (0-11) |
1483 |
* @param ivk Viking-chip number (0-23) |
* @param ivk Viking-chip number (0-23) |
1484 |
*/ |
*/ |
1485 |
Bool_t TrkLevel2::IsMaskedVK(int iv, int ivk){ |
Bool_t TrkLevel2::IsMaskedVK(int iv, int ivk){ |
1486 |
return !(GetVKMask(iv,ivk)&&GetVKFlag(iv,ivk) ); |
return !(GetVKMask(iv,ivk)&&GetVKFlag(iv,ivk) ); |
1565 |
t_track->ayv[ip] = l2->ayv_nt[i][ip]; |
t_track->ayv[ip] = l2->ayv_nt[i][ip]; |
1566 |
t_track->dedx_x[ip] = l2->dedx_x[i][ip]; |
t_track->dedx_x[ip] = l2->dedx_x[i][ip]; |
1567 |
t_track->dedx_y[ip] = l2->dedx_y[i][ip]; |
t_track->dedx_y[ip] = l2->dedx_y[i][ip]; |
1568 |
|
t_track->multmaxx[ip] = l2->multmaxx[i][ip]; |
1569 |
|
t_track->multmaxy[ip] = l2->multmaxy[i][ip]; |
1570 |
|
t_track->seedx[ip] = l2->seedx[i][ip]; |
1571 |
|
t_track->seedy[ip] = l2->seedy[i][ip]; |
1572 |
|
t_track->xpu[ip] = l2->xpu[i][ip]; |
1573 |
|
t_track->ypu[ip] = l2->ypu[i][ip]; |
1574 |
//----------------------------------------------------- |
//----------------------------------------------------- |
1575 |
//----------------------------------------------------- |
//----------------------------------------------------- |
1576 |
//----------------------------------------------------- |
//----------------------------------------------------- |
1577 |
//----------------------------------------------------- |
//----------------------------------------------------- |
1578 |
}; |
}; |
1579 |
|
// ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ |
1580 |
|
// evaluated coordinates (to define GF) |
1581 |
|
// ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ |
1582 |
|
int ngf = TrkParams::nGF; |
1583 |
|
float *zgf = TrkParams::zGF; |
1584 |
|
Trajectory tgf = Trajectory(ngf,zgf); |
1585 |
|
tgf.DoTrack2(t_track->al);//<<<< integrate the trajectory |
1586 |
|
for(int ip=0; ip<ngf; ip++){ |
1587 |
|
t_track->xGF[ip] = tgf.x[ip]; |
1588 |
|
t_track->yGF[ip] = tgf.y[ip]; |
1589 |
|
} |
1590 |
|
|
1591 |
// if(t_track->IsSaturated())t_track->Dump(); |
// if(t_track->IsSaturated())t_track->Dump(); |
1592 |
new(t[i]) TrkTrack(*t_track); |
new(t[i]) TrkTrack(*t_track); |
1593 |
t_track->Clear(); |
t_track->Clear(); |
1594 |
}; |
};//end loop over track |
1595 |
|
|
1596 |
// ---------------- |
// ---------------- |
1597 |
// *** SINGLETS *** |
// *** SINGLETS *** |
1603 |
t_singlet->coord[0] = l2->xs[i][0]; |
t_singlet->coord[0] = l2->xs[i][0]; |
1604 |
t_singlet->coord[1] = l2->xs[i][1]; |
t_singlet->coord[1] = l2->xs[i][1]; |
1605 |
t_singlet->sgnl = l2->signlxs[i]; |
t_singlet->sgnl = l2->signlxs[i]; |
1606 |
|
t_singlet->multmax = l2->multmaxsx[i]; |
1607 |
|
if(l2->sxbad[i]>0) t_singlet->multmax = -1*t_singlet->multmax; |
1608 |
//----------------------------------------------------- |
//----------------------------------------------------- |
1609 |
// if(l1) t_singlet->cls = l1->GetCluster(l2->clsx[i]-1); |
// if(l1) t_singlet->cls = l1->GetCluster(l2->clsx[i]-1); |
1610 |
//----------------------------------------------------- |
//----------------------------------------------------- |
1618 |
t_singlet->coord[0] = l2->ys[i][0]; |
t_singlet->coord[0] = l2->ys[i][0]; |
1619 |
t_singlet->coord[1] = l2->ys[i][1]; |
t_singlet->coord[1] = l2->ys[i][1]; |
1620 |
t_singlet->sgnl = l2->signlys[i]; |
t_singlet->sgnl = l2->signlys[i]; |
1621 |
|
t_singlet->multmax = l2->multmaxsy[i]; |
1622 |
|
if(l2->sybad[i]>0) t_singlet->multmax = -1*t_singlet->multmax; |
1623 |
//----------------------------------------------------- |
//----------------------------------------------------- |
1624 |
// if(l1) t_singlet->cls = l1->GetCluster(l2->clsy[i]-1); |
// if(l1) t_singlet->cls = l1->GetCluster(l2->clsy[i]-1); |
1625 |
//----------------------------------------------------- |
//----------------------------------------------------- |
1626 |
new(sy[i]) TrkSinglet(*t_singlet); |
new(sy[i]) TrkSinglet(*t_singlet); |
1627 |
t_singlet->Clear(); |
t_singlet->Clear(); |
1628 |
}; |
}; |
1629 |
|
|
1630 |
|
|
1631 |
|
|
1632 |
delete t_track; |
delete t_track; |
1633 |
delete t_singlet; |
delete t_singlet; |
1803 |
TrkTrack *TrkLevel2::GetStoredTrack(int is){ |
TrkTrack *TrkLevel2::GetStoredTrack(int is){ |
1804 |
|
|
1805 |
if(is >= this->ntrk()){ |
if(is >= this->ntrk()){ |
1806 |
cout << "** TrkLevel2 ** Track "<< is << "doen not exits! " << endl; |
cout << "TrkTrack *TrkLevel2::GetStoredTrack(int) >> Track "<< is << "doen not exits! " << endl; |
1807 |
cout << " Stored tracks ntrk() = "<< this->ntrk() << endl; |
cout << "Stored tracks ntrk() = "<< this->ntrk() << endl; |
1808 |
return 0; |
return 0; |
1809 |
} |
} |
1810 |
if(!Track){ |
if(!Track){ |
1825 |
TrkSinglet *TrkLevel2::GetSingletX(int is){ |
TrkSinglet *TrkLevel2::GetSingletX(int is){ |
1826 |
|
|
1827 |
if(is >= this->nclsx()){ |
if(is >= this->nclsx()){ |
1828 |
cout << "** TrkLevel2 ** Singlet "<< is << "doen not exits! " << endl; |
cout << "TrkSinglet *TrkLevel2::GetSingletX(int) >> Singlet "<< is << "doen not exits! " << endl; |
1829 |
cout << " Stored x-singlets nclsx() = "<< this->nclsx() << endl; |
cout << "Stored x-singlets nclsx() = "<< this->nclsx() << endl; |
1830 |
return 0; |
return 0; |
1831 |
} |
} |
1832 |
if(!SingletX)return 0; |
if(!SingletX)return 0; |
1845 |
TrkSinglet *TrkLevel2::GetSingletY(int is){ |
TrkSinglet *TrkLevel2::GetSingletY(int is){ |
1846 |
|
|
1847 |
if(is >= this->nclsy()){ |
if(is >= this->nclsy()){ |
1848 |
cout << "** TrkLevel2 ** Singlet "<< is << "doen not exits! " << endl; |
cout << "TrkSinglet *TrkLevel2::GetSingletY(int) >> Singlet "<< is << "doen not exits! " << endl; |
1849 |
cout << " Stored y-singlets nclsy() = "<< this->nclsx() << endl; |
cout << "Stored y-singlets nclsx() = "<< this->nclsx() << endl; |
1850 |
return 0; |
return 0; |
1851 |
} |
} |
1852 |
if(!SingletY)return 0; |
if(!SingletY)return 0; |
1866 |
TrkTrack *TrkLevel2::GetTrack(int it){ |
TrkTrack *TrkLevel2::GetTrack(int it){ |
1867 |
|
|
1868 |
if(it >= this->GetNTracks()){ |
if(it >= this->GetNTracks()){ |
1869 |
cout << "** TrkLevel2 ** Track "<< it << "does not exits! " << endl; |
cout << "TrkTrack *TrkLevel2::GetTrack(int) >> Track "<< it << "does not exits! " << endl; |
1870 |
cout << " Physical tracks GetNTracks() = "<< this->ntrk() << endl; |
cout << "Physical tracks GetNTracks() = "<< this->ntrk() << endl; |
1871 |
return 0; |
return 0; |
1872 |
} |
} |
1873 |
|
|
1904 |
TrkTrack *TrkLevel2::GetTrackImage(int it){ |
TrkTrack *TrkLevel2::GetTrackImage(int it){ |
1905 |
|
|
1906 |
if(it >= this->GetNTracks()){ |
if(it >= this->GetNTracks()){ |
1907 |
cout << "** TrkLevel2 ** Track "<< it << "does not exits! " << endl; |
cout << "TrkTrack *TrkLevel2::GetTrackImage(int) >> Track "<< it << "does not exits! " << endl; |
1908 |
cout << " Physical tracks GetNTracks() = "<< this->ntrk() << endl; |
cout << "Physical tracks GetNTracks() = "<< this->ntrk() << endl; |
1909 |
return 0; |
return 0; |
1910 |
} |
} |
1911 |
|
|
1914 |
TrkTrack *track = (TrkTrack*)sorted->At(it); |
TrkTrack *track = (TrkTrack*)sorted->At(it); |
1915 |
|
|
1916 |
if(!track->HasImage()){ |
if(!track->HasImage()){ |
1917 |
cout << "** TrkLevel2 ** Track "<< it << "does not have image! " << endl; |
cout << "TrkTrack *TrkLevel2::GetTrackImage(int) >> Track "<< it << "does not have image! " << endl; |
1918 |
return 0; |
return 0; |
1919 |
} |
} |
1920 |
if(!Track)return 0; |
if(!Track)return 0; |
1941 |
// path_.error = 0; |
// path_.error = 0; |
1942 |
// readb_(); |
// readb_(); |
1943 |
|
|
1944 |
TrkParams::SetTrackingMode(); |
// TrkParams::SetTrackingMode(); |
1945 |
TrkParams::SetPrecisionFactor(); |
// TrkParams::SetPrecisionFactor(); |
1946 |
TrkParams::SetStepMin(); |
// TrkParams::SetStepMin(); |
1947 |
|
TrkParams::SetMiniDefault(); |
1948 |
|
|
1949 |
TrkParams::Set(path,1); |
TrkParams::Set(path,1); |
1950 |
TrkParams::Load(1); |
TrkParams::Load(1); |