9 |
|
|
10 |
#include "GeoFieldCut.h" |
#include "GeoFieldCut.h" |
11 |
|
|
12 |
int GeoFieldCut::Check(PamLevel2 * event) { |
int GeoFieldCut::Check(PamLevel2 *event) { |
13 |
|
|
14 |
if (event->GetOrbitalInfo()->Babs > _bAbs) // Babs = absolute value of geomagnetic field (Gauss) |
Float_t bAbs = event->GetOrbitalInfo()->Babs; |
15 |
return CUTOK; |
if (bAbs < _bAbsMin) // Babs = absolute value of geomagnetic field (Gauss) |
|
else |
|
16 |
return 0; |
return 0; |
17 |
|
if (_bAbsMax > _bAbsMin && bAbs > _bAbsMax) |
18 |
|
return 0; |
19 |
|
|
20 |
|
//cout << bAbs << endl; |
21 |
|
return CUTOK; |
22 |
} |
} |