Skip to content
Closed
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
6 changes: 3 additions & 3 deletions ALICE3/TableProducer/OTF/onTheFlyRICHPID.cxx
Original file line number Diff line number Diff line change
Expand Up @@ -226,19 +226,19 @@
l_aerogel_z[i_central_mirror] = std::sqrt(1.0 + m_val * m_val) * R_min * square_size_barrel_cylinder / (std::sqrt(1.0 + m_val * m_val) * R_max - m_val * square_size_barrel_cylinder);
T_r_plus_g[i_central_mirror] = R_max - R_min;
float t = std::tan(std::atan(m_val) + std::atan(square_size_barrel_cylinder / (2.0 * R_max * std::sqrt(1.0 + m_val * m_val) - square_size_barrel_cylinder * m_val)));
theta_max[i_central_mirror] = M_PI / 2.0 - std::atan(t);

Check failure on line 229 in ALICE3/TableProducer/OTF/onTheFlyRICHPID.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[external-pi]

Use the PI constant (and its multiples and fractions) defined in o2::constants::math.
theta_min[i_central_mirror] = M_PI / 2.0 + std::atan(t);

Check failure on line 230 in ALICE3/TableProducer/OTF/onTheFlyRICHPID.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[external-pi]

Use the PI constant (and its multiples and fractions) defined in o2::constants::math.
mProjectiveLengthInner = R_min * t;
aerogel_rindex[i_central_mirror] = bRichRefractiveIndexSector[0];
for (int i = i_central_mirror + 1; i < number_of_sectors_in_z; i++) {
float par_a = t;
float par_b = 2.0 * R_max / square_size_z;
m_val = (std::sqrt(par_a * par_a * par_b * par_b + par_b * par_b - 1.0) + par_a * par_b * par_b) / (par_b * par_b - 1.0);
theta_min[i] = M_PI / 2.0 - std::atan(t);

Check failure on line 237 in ALICE3/TableProducer/OTF/onTheFlyRICHPID.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[external-pi]

Use the PI constant (and its multiples and fractions) defined in o2::constants::math.
theta_max[2 * i_central_mirror - i] = M_PI / 2.0 + std::atan(t);

Check failure on line 238 in ALICE3/TableProducer/OTF/onTheFlyRICHPID.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[external-pi]

Use the PI constant (and its multiples and fractions) defined in o2::constants::math.
t = std::tan(std::atan(m_val) + std::atan(square_size_z / (2.0 * R_max * std::sqrt(1.0 + m_val * m_val) - square_size_z * m_val)));
theta_max[i] = M_PI / 2.0 - std::atan(t);

Check failure on line 240 in ALICE3/TableProducer/OTF/onTheFlyRICHPID.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[external-pi]

Use the PI constant (and its multiples and fractions) defined in o2::constants::math.
theta_min[2 * i_central_mirror - i] = M_PI / 2.0 + std::atan(t);

Check failure on line 241 in ALICE3/TableProducer/OTF/onTheFlyRICHPID.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[external-pi]

Use the PI constant (and its multiples and fractions) defined in o2::constants::math.
// Forward sectors
theta_bi[i] = std::atan(m_val);
R0_tilt[i] = R_max - square_size_z / 2.0 * std::sin(std::atan(m_val));
Expand Down Expand Up @@ -269,11 +269,11 @@
float par_a = t;
float par_b = 2.0 * R_max / square_size_z;
m_val = (std::sqrt(par_a * par_a * par_b * par_b + par_b * par_b - 1.0) + par_a * par_b * par_b) / (par_b * par_b - 1.0);
theta_min[i] = M_PI / 2.0 - std::atan(t);

Check failure on line 272 in ALICE3/TableProducer/OTF/onTheFlyRICHPID.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[external-pi]

Use the PI constant (and its multiples and fractions) defined in o2::constants::math.
theta_max[2 * i_central_mirror - i - 1] = M_PI / 2.0 + std::atan(t);

Check failure on line 273 in ALICE3/TableProducer/OTF/onTheFlyRICHPID.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[external-pi]

Use the PI constant (and its multiples and fractions) defined in o2::constants::math.
t = std::tan(std::atan(m_val) + std::atan(square_size_z / (2.0 * R_max * std::sqrt(1.0 + m_val * m_val) - square_size_z * m_val)));
theta_max[i] = M_PI / 2.0 - std::atan(t);

Check failure on line 275 in ALICE3/TableProducer/OTF/onTheFlyRICHPID.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[external-pi]

Use the PI constant (and its multiples and fractions) defined in o2::constants::math.
theta_min[2 * i_central_mirror - i - 1] = M_PI / 2.0 + std::atan(t);

Check failure on line 276 in ALICE3/TableProducer/OTF/onTheFlyRICHPID.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[external-pi]

Use the PI constant (and its multiples and fractions) defined in o2::constants::math.
// Forward sectors
theta_bi[i] = std::atan(m_val);
R0_tilt[i] = R_max - square_size_z / 2.0 * std::sin(std::atan(m_val));
Expand Down Expand Up @@ -308,11 +308,11 @@
gap_thickness[i] = (det_centers[i] - rad_centers[i]).Mag() - bRichRadiatorThickness / 2.0;
}
// DEBUG
//std::cout << std::endl << std::endl;
//for (int i = 0; i < number_of_sectors_in_z; i++) {
// std::cout << std::endl << std::endl;
// for (int i = 0; i < number_of_sectors_in_z; i++) {
// std::cout << "Sector " << i << ": Gap = " << gap_thickness[i] << " cm, Index aerogel = " << aerogel_rindex[i] << ", Gap thickness = " << gap_thickness[i] << " cm" << std::endl;
//}
//std::cout << std::endl << std::endl;
// std::cout << std::endl << std::endl;
}

void init(o2::framework::InitContext& /*initContext*/)
Expand Down
Loading