| 1 |
pam-fi |
1.1 |
/* |
| 2 |
|
|
* TrkDedxHCut.cpp |
| 3 |
|
|
* |
| 4 |
pam-fi |
1.2 |
* Created on: 15-mar-2010 |
| 5 |
|
|
* Author: S. Ricciarini |
| 6 |
pam-fi |
1.1 |
*/ |
| 7 |
|
|
|
| 8 |
|
|
/*! @file TrkDedxHCut.cpp The TrkDedxHCut class implementation file */ |
| 9 |
|
|
|
| 10 |
|
|
#include "TrkDedxHCut.h" |
| 11 |
|
|
|
| 12 |
|
|
const Double_t TrkDedxHCut::_xLow[] = { 0.126414, 0.140142, 0.152714, 0.165854, 0.178357, 0.194244, 0.211669, 0.242859, |
| 13 |
|
|
0.264646, 0.293385, 0.314258, 0.360565, 0.392911, 0.450807, 0.508423, 0.603735, 0.680897, 0.754838, 0.851312, |
| 14 |
|
|
1.00032, 1.36201, 1.37, 1000. }; // _nLow=23 elements |
| 15 |
|
|
|
| 16 |
|
|
const Double_t TrkDedxHCut::_yLow[] = { 18.9221, 16.0752, 14.2214, 12.6324, 11.2421, 9.71928, 8.1965, 6.93856, 6.14407, |
| 17 |
|
|
5.15095, 4.6875, 3.8268, 3.49576, 2.89989, 2.56886, 2.17161, 1.90678, 1.70816, 1.50953, 1.27754, 0.715042, 0, 0 }; // _nLow=23 elements |
| 18 |
|
|
|
| 19 |
pam-fi |
1.2 |
Int_t TrkDedxHCut::Check(PamLevel2 *event) { |
| 20 |
pam-fi |
1.1 |
|
| 21 |
|
|
TrkTrack *trkTrack = event->GetTrack(0)->GetTrkTrack(); |
| 22 |
|
|
|
| 23 |
|
|
Float_t rigMod = trkTrack->GetRigidity(); |
| 24 |
|
|
Float_t dedx = trkTrack->GetDEDX(); |
| 25 |
|
|
|
| 26 |
pam-fi |
1.2 |
if (rigMod < 0.126414) |
| 27 |
pam-fi |
1.1 |
return 0; // a proton below 126.414 MV is rejected |
| 28 |
|
|
|
| 29 |
pam-fi |
1.2 |
// lower limit taken from N. De Simone (interpolation of _xLow) |
| 30 |
|
|
|
| 31 |
|
|
for (Int_t i = 0; i < _nLow - 1; i++) { |
| 32 |
pam-fi |
1.1 |
if (_xLow[i] <= rigMod && rigMod < _xLow[i + 1]) { |
| 33 |
pam-fi |
1.2 |
Double_t slope = (_yLow[i + 1] - _yLow[i]) / (_xLow[i + 1] - _xLow[i]); |
| 34 |
|
|
Double_t y = _yLow[i] + slope * (rigMod - _xLow[i]); |
| 35 |
pam-fi |
1.1 |
|
| 36 |
|
|
if (dedx < y) |
| 37 |
|
|
return 0; |
| 38 |
|
|
} |
| 39 |
|
|
} |
| 40 |
|
|
|
| 41 |
pam-fi |
1.2 |
// higher limit taken from N. Mori lower limit for Helium |
| 42 |
pam-fi |
1.1 |
|
| 43 |
pam-fi |
1.2 |
// Compute dE/dx for X and Y separately |
| 44 |
|
|
Float_t totReleaseX = 0.; |
| 45 |
|
|
Float_t totReleaseY = 0.; |
| 46 |
|
|
Int_t nX = 0, nY = 0; |
| 47 |
|
|
Float_t dEdxView = 0; |
| 48 |
|
|
|
| 49 |
|
|
for (Int_t ip = 0; ip < 6; ip++) { |
| 50 |
|
|
// X view |
| 51 |
|
|
dEdxView = trkTrack->GetDEDX(ip, 0); |
| 52 |
|
|
if (dEdxView > 0 && trkTrack->XGood(ip)) { |
| 53 |
|
|
totReleaseX += dEdxView; |
| 54 |
|
|
nX++; |
| 55 |
|
|
} |
| 56 |
|
|
// Y view |
| 57 |
|
|
dEdxView = trkTrack->GetDEDX(ip, 1); |
| 58 |
|
|
if (dEdxView > 0 && trkTrack->YGood(ip)) { |
| 59 |
|
|
totReleaseY += dEdxView; |
| 60 |
|
|
nY++; |
| 61 |
pam-fi |
1.1 |
} |
| 62 |
|
|
} |
| 63 |
|
|
|
| 64 |
pam-fi |
1.2 |
Float_t dEdxX = totReleaseX / nX; // nX assumed > 0 |
| 65 |
|
|
Float_t dEdxY = totReleaseY / nY; // nY assumed > 0 |
| 66 |
|
|
|
| 67 |
|
|
Float_t denHigh = pow(rigMod, 1.8); // The power of the denominator is the same for X and Y (with current calibration) |
| 68 |
|
|
if (dEdxX > 3.7 + 4.6 / denHigh) { |
| 69 |
|
|
return 0; |
| 70 |
|
|
} |
| 71 |
|
|
if (dEdxY > 3.3 + 4.9 / denHigh) { |
| 72 |
|
|
return 0; |
| 73 |
pam-fi |
1.1 |
} |
| 74 |
pam-fi |
1.2 |
|
| 75 |
pam-fi |
1.1 |
return CUTOK; |
| 76 |
|
|
|
| 77 |
|
|
} |
Parent Directory
| 
Revision Log