|
76 | 76 |
|
77 | 77 | #include <math.h> // FIXME: Replace M_PI |
78 | 78 |
|
79 | | -using namespace std; |
80 | 79 | using namespace o2; |
81 | 80 | using namespace o2::aod; |
82 | 81 | using namespace o2::mch; |
@@ -169,9 +168,9 @@ struct VarTrack { |
169 | 168 | }; |
170 | 169 |
|
171 | 170 | struct VarClusters { |
172 | | - vector<vector<float>> posClusters; // (x,y,z) |
173 | | - vector<vector<float>> errorClusters; // (ex,ey) |
174 | | - vector<int> DEIDs; |
| 171 | + std::vector<std::vector<float>> posClusters; // (x,y,z) |
| 172 | + std::vector<std::vector<float>> errorClusters; // (ex,ey) |
| 173 | + std::vector<int> DEIDs; |
175 | 174 | }; |
176 | 175 |
|
177 | 176 | struct muonQa { |
@@ -1399,7 +1398,7 @@ struct muonQa { |
1399 | 1398 | { |
1400 | 1399 | int removable = 0; |
1401 | 1400 | auto clustersSliced = mchcls.sliceBy(perMuon, muon.globalIndex()); // Slice clusters by muon id |
1402 | | - vector<vector<float>> posClusters; |
| 1401 | + std::vector<std::vector<float>> posClusters; |
1403 | 1402 |
|
1404 | 1403 | int clIndex = -1; |
1405 | 1404 | // Get re-aligned clusters associated to current track |
@@ -1511,9 +1510,9 @@ struct muonQa { |
1511 | 1510 | fgValues.phi = propmuon.getPhi(); |
1512 | 1511 |
|
1513 | 1512 | fgValues.p = propmuon.getP(); |
1514 | | - fgValues.px = propmuon.getP() * sin(M_PI / 2 - atan(propmuon.getTgl())) * cos(propmuon.getPhi()); |
1515 | | - fgValues.py = propmuon.getP() * sin(M_PI / 2 - atan(propmuon.getTgl())) * sin(propmuon.getPhi()); |
1516 | | - fgValues.pz = propmuon.getP() * cos(M_PI / 2 - atan(propmuon.getTgl())); |
| 1513 | + fgValues.px = propmuon.getP() * std::sin(M_PI / 2 - std::atan(propmuon.getTgl())) * std::cos(propmuon.getPhi()); |
| 1514 | + fgValues.py = propmuon.getP() * std::sin(M_PI / 2 - std::atan(propmuon.getTgl())) * std::sin(propmuon.getPhi()); |
| 1515 | + fgValues.pz = propmuon.getP() * std::cos(M_PI / 2 - std::atan(propmuon.getTgl())); |
1517 | 1516 | } |
1518 | 1517 |
|
1519 | 1518 | if (endPoint == kToDCA) { |
@@ -1583,8 +1582,8 @@ struct muonQa { |
1583 | 1582 | double pMCH = fgValuesMCH.p; |
1584 | 1583 | int sign = fgValuesMCH.sign; |
1585 | 1584 |
|
1586 | | - double px = pMCH * sin(M_PI / 2 - atan(muon.tgl())) * cos(muon.phi()); |
1587 | | - double py = pMCH * sin(M_PI / 2 - atan(muon.tgl())) * sin(muon.phi()); |
| 1585 | + double px = pMCH * std::sin(M_PI / 2 - std::atan(muon.tgl())) * std::cos(muon.phi()); |
| 1586 | + double py = pMCH * std::sin(M_PI / 2 - std::atan(muon.tgl())) * std::sin(muon.phi()); |
1588 | 1587 | // double pz = pMCH * cos(M_PI / 2 - atan(mft.tgl())); |
1589 | 1588 | double pt = std::sqrt(std::pow(px, 2) + std::pow(py, 2)); |
1590 | 1589 |
|
@@ -1654,9 +1653,9 @@ struct muonQa { |
1654 | 1653 | fgValues.phi = propmuon.getPhi(); |
1655 | 1654 |
|
1656 | 1655 | fgValues.p = propmuon.getP(); |
1657 | | - fgValues.px = propmuon.getP() * sin(M_PI / 2 - atan(propmuon.getTgl())) * cos(propmuon.getPhi()); |
1658 | | - fgValues.py = propmuon.getP() * sin(M_PI / 2 - atan(propmuon.getTgl())) * sin(propmuon.getPhi()); |
1659 | | - fgValues.pz = propmuon.getP() * cos(M_PI / 2 - atan(propmuon.getTgl())); |
| 1656 | + fgValues.px = propmuon.getP() * std::sin(M_PI / 2 - std::atan(propmuon.getTgl())) * std::cos(propmuon.getPhi()); |
| 1657 | + fgValues.py = propmuon.getP() * std::sin(M_PI / 2 - std::atan(propmuon.getTgl())) * std::sin(propmuon.getPhi()); |
| 1658 | + fgValues.pz = propmuon.getP() * std::cos(M_PI / 2 - std::atan(propmuon.getTgl())); |
1660 | 1659 | } |
1661 | 1660 |
|
1662 | 1661 | if (endPoint == kToDCA) { |
|
0 commit comments