16 #include "CLHEP/Units/SystemOfUnits.h"
324 std::cout <<
"FATAL TRT_BarrelDriftTimeData::TRT_BarrelDriftTimeData strawGas ("
325 <<
m_strawGas <<
") must be 0(Xe), 1(Kr) or 2(Ar). The job will die now :(" << std::endl;
326 throw std::runtime_error(
"TRT_BarrelDriftTimeData: unsupported gas type");
348 unsigned int indexHigh = nTabulated-1;
349 for (
unsigned int testIndex = 0; testIndex < nTabulated-1; testIndex++) {
351 indexHigh = testIndex;
355 unsigned int indexLow = 0;
357 std::cout <<
"TRT_BarrelDriftTimeData::DriftTimeAtNoField::Warning"
358 <<
" - distance is lower than first tabulated distance" << std::endl;
362 indexLow = indexHigh - 1;
368 double lowDistWeight = 1.0-(dist - distLow)/(distHigh-distLow);
369 double highDistWeight = 1.0-(distHigh - dist)/(distHigh-distLow);
385 unsigned int indexHigh = nTabulated-1;
386 for (
unsigned int testIndex = 0; testIndex < nTabulated-1; testIndex++) {
388 indexHigh = testIndex;
392 unsigned int indexLow = 0;
394 std::cout <<
"TRTDT::Warning - distance is lower than first tabulated distance" << std::endl;
398 indexLow = indexHigh - 1;
404 double lowDistWeight = 1.0-(dist - distLow)/(distHigh-distLow);
405 double highDistWeight = 1.0-(distHigh - dist)/(distHigh-distLow);