| 68 |
float rigThreshold = _thresholdCoeff * event->GetOrbitalInfo()->GetCutoffSVL(); |
float rigThreshold = _thresholdCoeff * event->GetOrbitalInfo()->GetCutoffSVL(); |
| 69 |
float rig = event->GetTrack(0)->GetTrkTrack()->GetRigidity(); |
float rig = event->GetTrack(0)->GetTrkTrack()->GetRigidity(); |
| 70 |
|
|
| 71 |
_rootHisto.Fill(rigThreshold, rig); // X,Y in the ROOT histogram |
//_rootHisto.Fill(rigThreshold, rig); // X,Y in the ROOT histogram |
| 72 |
|
|
| 73 |
// Event rigidity bin identification. |
// Event rigidity bin identification. |
| 74 |
|
|
| 84 |
|
|
| 85 |
// Threshold rigidity bin identification. |
// Threshold rigidity bin identification. |
| 86 |
|
|
| 87 |
if (rigThreshold < 0 || rigThreshold >= _bins[(int) _textHisto.GetNRows()]) { |
if (rigThreshold <= 0 || rigThreshold >= _bins[(int) _textHisto.GetNRows()]) { |
| 88 |
return; |
return; |
| 89 |
} |
} |
| 90 |
|
|
| 91 |
if (0 <= rigThreshold && rigThreshold < _bins[0]) |
if (0 < rigThreshold && rigThreshold < _bins[0]) |
| 92 |
_zeroCutoffBins[i]++; |
_zeroCutoffBins[i]++; |
| 93 |
else { |
else { |
| 94 |
int j = 1; |
int j = 1; |
| 124 |
outTextFile.close(); |
outTextFile.close(); |
| 125 |
|
|
| 126 |
// Write the ROOT file |
// Write the ROOT file |
| 127 |
TFile outRootFile((_outFileBase + "-" + GetName() + ".root"), "RECREATE"); |
/*TFile outRootFile((_outFileBase + "-" + GetName() + ".root"), "RECREATE"); |
| 128 |
outRootFile.cd(); |
outRootFile.cd(); |
| 129 |
_rootHisto.Write(); |
_rootHisto.Write(); |
| 130 |
outRootFile.Close(); |
outRootFile.Close(); |
| 131 |
|
*/ |
| 132 |
} |
} |
| 133 |
} |
} |