16 #include "CLHEP/Units/SystemOfUnits.h"
38 std::cout <<
"FATAL TRT_BarrelDriftTimeData::TRT_BarrelDriftTimeData digversion < 11 ("
39 << digversion <<
") is no longer supported. The job will die now :(" << std::endl;
40 throw std::runtime_error(
"TRT_BarrelDriftTimeData: unsupported digversion");
330 std::cout <<
"FATAL TRT_BarrelDriftTimeData::TRT_BarrelDriftTimeData strawGas ("
331 <<
m_strawGas <<
") must be 0(Xe), 1(Kr) or 2(Ar). The job will die now :(" << std::endl;
332 throw std::runtime_error(
"TRT_BarrelDriftTimeData: unsupported gas type");
354 unsigned int indexHigh = nTabulated-1;
355 for (
unsigned int testIndex = 0; testIndex < nTabulated-1; testIndex++) {
357 indexHigh = testIndex;
361 unsigned int indexLow = 0;
363 std::cout <<
"TRT_BarrelDriftTimeData::DriftTimeAtNoField::Warning"
364 <<
" - distance is lower than first tabulated distance" << std::endl;
368 indexLow = indexHigh - 1;
374 double lowDistWeight = 1.0-(dist - distLow)/(distHigh-distLow);
375 double highDistWeight = 1.0-(distHigh - dist)/(distHigh-distLow);
391 unsigned int indexHigh = nTabulated-1;
392 for (
unsigned int testIndex = 0; testIndex < nTabulated-1; testIndex++) {
394 indexHigh = testIndex;
398 unsigned int indexLow = 0;
400 std::cout <<
"TRTDT::Warning - distance is lower than first tabulated distance" << std::endl;
404 indexLow = indexHigh - 1;
410 double lowDistWeight = 1.0-(dist - distLow)/(distHigh-distLow);
411 double highDistWeight = 1.0-(distHigh - dist)/(distHigh-distLow);