11 |
long int RefitTrackAction::_currEvent = -1; |
long int RefitTrackAction::_currEvent = -1; |
12 |
TString RefitTrackAction::_currTrkParamsFile = ""; |
TString RefitTrackAction::_currTrkParamsFile = ""; |
13 |
|
|
14 |
|
float RefitTrackAction::_xmBackup[6], RefitTrackAction::_ymBackup[6]; |
15 |
|
int RefitTrackAction::_layerXBackup[6], RefitTrackAction::_layerYBackup[6]; |
16 |
|
|
17 |
void RefitTrackAction::OnGood(PamLevel2 *event) { |
void RefitTrackAction::OnGood(PamLevel2 *event) { |
18 |
|
|
19 |
// Do something only if at least one track exists |
// Do something only if at least one track exists |
59 |
} |
} |
60 |
|
|
61 |
// Apply the gaussian spread |
// Apply the gaussian spread |
62 |
|
|
63 |
|
cout << GetName() << " " << trkTrack->xm[0]; |
64 |
|
cout.flush(); |
65 |
for (Int_t ip = 0; ip < 6; ip++) { |
for (Int_t ip = 0; ip < 6; ip++) { |
66 |
// gaussian distributed random number Gauss(mean,sigma) (in cm) |
// gaussian distributed random number Gauss(mean,sigma) (in cm) |
67 |
trkTrack->xm[ip] += _r3->Gaus(0., _spreadX); |
trkTrack->xm[ip] += _r3->Gaus(0., _spreadX); |
68 |
trkTrack->ym[ip] += _r3->Gaus(0., _spreadY); |
trkTrack->ym[ip] += _r3->Gaus(0., _spreadY); |
69 |
} |
} |
70 |
|
cout << GetName() << " " << trkTrack->xm[0] << endl; |
71 |
// Mask the layers |
// Mask the layers |
72 |
if (_maskedViews) { |
if (_maskedViews) { |
73 |
if ((_maskedViews & T1X) == T1X) |
if ((_maskedViews & T1X) == T1X) |