Skip to content
Snippets Groups Projects
Commit b7d4b877 authored by Miha Muskinja's avatar Miha Muskinja Committed by Miha Muskinja
Browse files

replace sincos

parent 6801e614
7 merge requests!69091Fix correlated smearing bug in JER in JetUncertainties in 22.0,!58791DataQualityConfigurations: Modify L1Calo config for web display,!51674Fixing hotSpotInHIST for Run3 HIST,!50012RecExConfig: Adjust log message levels from GetRunNumber and GetLBNumber,!46784MuonCondInterface: Enable thread-safety checking.,!46776Updated LArMonitoring config file for WD to match new files produced using MT,!40107Vectorized sincos calculation in LArWheelCalculator
...@@ -48,7 +48,11 @@ namespace LArWheelCalculator_Impl ...@@ -48,7 +48,11 @@ namespace LArWheelCalculator_Impl
const double cos_a = scalpha.cs, sin_a = scalpha.sn; const double cos_a = scalpha.cs, sin_a = scalpha.sn;
#else // parameterized sine #else // parameterized sine
double cos_a, sin_a; double cos_a, sin_a;
lwc()->parameterized_sincos(P.y(), sin_a, cos_a); lwc()->m_vsincos_par.eval(P.y(), sin_a, cos_a);
// std::cout << sin_a << " " << cos_a << std::endl;
// double cos_a_old, sin_a_old;
// lwc()->parameterized_sincos(P.y(), sin_a_old, cos_a_old);
// std::cout << "old " << sin_a_old << " " << cos_a_old << std::endl;
#endif #endif
// determination of the nearest quarter-wave number // determination of the nearest quarter-wave number
int nqwave = (z < 0.) ? 0 : int(z / lwc()->m_QuarterWaveLength); int nqwave = (z < 0.) ? 0 : int(z / lwc()->m_QuarterWaveLength);
...@@ -160,7 +164,11 @@ namespace LArWheelCalculator_Impl ...@@ -160,7 +164,11 @@ namespace LArWheelCalculator_Impl
double cos_a = scalpha.cs, sin_a = scalpha.sn; double cos_a = scalpha.cs, sin_a = scalpha.sn;
#else // parameterized sine #else // parameterized sine
double cos_a, sin_a; double cos_a, sin_a;
lwc()->parameterized_sincos(P.y(), sin_a, cos_a); lwc()->m_vsincos_par.eval(P.y(), sin_a, cos_a);
// std::cout << sin_a << " " << cos_a << std::endl;
// double cos_a_old, sin_a_old;
// lwc()->parameterized_sincos(P.y(), sin_a_old, cos_a_old);
// std::cout << "old " << sin_a_old << " " << cos_a_old << std::endl;
#endif #endif
bool sqw = false; bool sqw = false;
...@@ -251,7 +259,11 @@ namespace LArWheelCalculator_Impl ...@@ -251,7 +259,11 @@ namespace LArWheelCalculator_Impl
const double cos_a = scalpha.cs, sin_a = scalpha.sn; const double cos_a = scalpha.cs, sin_a = scalpha.sn;
#else // parameterized sine #else // parameterized sine
double cos_a, sin_a; double cos_a, sin_a;
lwc()->parameterized_sincos(P.y(), sin_a, cos_a); lwc()->m_vsincos_par.eval(P.y(), sin_a, cos_a);
// std::cout << sin_a << " " << cos_a << std::endl;
// double cos_a_old, sin_a_old;
// lwc()->parameterized_sincos(P.y(), sin_a_old, cos_a_old);
// std::cout << "old " << sin_a_old << " " << cos_a_old << std::endl;
#endif #endif
int nqwave; int nqwave;
...@@ -339,7 +351,11 @@ namespace LArWheelCalculator_Impl ...@@ -339,7 +351,11 @@ namespace LArWheelCalculator_Impl
double cos_a = scalpha.cs, sin_a = scalpha.sn; double cos_a = scalpha.cs, sin_a = scalpha.sn;
#else // parameterized sine #else // parameterized sine
double cos_a, sin_a; double cos_a, sin_a;
lwc()->parameterized_sincos(P.y(), sin_a, cos_a); lwc()->m_vsincos_par.eval(P.y(), sin_a, cos_a);
// std::cout << sin_a << " " << cos_a << std::endl;
// double cos_a_old, sin_a_old;
// lwc()->parameterized_sincos(P.y(), sin_a_old, cos_a_old);
// std::cout << "old " << sin_a_old << " " << cos_a_old << std::endl;
#endif #endif
bool sqw = false; bool sqw = false;
...@@ -434,7 +450,11 @@ namespace LArWheelCalculator_Impl ...@@ -434,7 +450,11 @@ namespace LArWheelCalculator_Impl
// parameterized sine // parameterized sine
#else #else
double cos_a, sin_a; double cos_a, sin_a;
lwc()->parameterized_sincos(P.y(), sin_a, cos_a); lwc()->m_vsincos_par.eval(P.y(), sin_a, cos_a);
// std::cout << sin_a << " " << cos_a << std::endl;
// double cos_a_old, sin_a_old;
// lwc()->parameterized_sincos(P.y(), sin_a_old, cos_a_old);
// std::cout << "old " << sin_a_old << " " << cos_a_old << std::endl;
#endif #endif
// determination of the nearest quarter-wave number // determination of the nearest quarter-wave number
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment