@@ -259,30 +259,30 @@ struct flowQC {
259259 float centrality = getCentrality (collision);
260260
261261 // EP method
262+ float QmodFT0A_EP = collision.qFT0A ();
262263 float psiFT0A_EP = collision.psiFT0A ();
263- float QxFT0A_EP = std::cos (2 * psiFT0A_EP);
264- float QyFT0A_EP = std::sin (2 * psiFT0A_EP);
265- float QmodFT0A_EP = std::hypot (QxFT0A_EP, QyFT0A_EP);
264+ float QxFT0A_EP = QmodFT0A_EP * std::cos (2 * psiFT0A_EP);
265+ float QyFT0A_EP = QmodFT0A_EP * std::sin (2 * psiFT0A_EP);
266266
267+ float QmodFT0C_EP = collision.qFT0C ();
267268 float psiFT0C_EP = collision.psiFT0C ();
268- float QxFT0C_EP = std::cos (2 * psiFT0C_EP);
269- float QyFT0C_EP = std::sin (2 * psiFT0C_EP);
270- float QmodFT0C_EP = std::hypot (QxFT0C_EP, QyFT0C_EP);
269+ float QxFT0C_EP = QmodFT0C_EP * std::cos (2 * psiFT0C_EP);
270+ float QyFT0C_EP = QmodFT0C_EP * std::sin (2 * psiFT0C_EP);
271271
272+ float QmodTPCl_EP = collision.qTPCL ();
272273 float psiTPCl_EP = collision.psiTPCL ();
273- float QxTPCl_EP = std::cos (2 * psiTPCl_EP);
274- float QyTPCl_EP = std::sin (2 * psiTPCl_EP);
275- float QmodTPCl_EP = std::hypot (QxTPCl_EP, QyTPCl_EP);
274+ float QxTPCl_EP = QmodTPCl_EP * std::cos (2 * psiTPCl_EP);
275+ float QyTPCl_EP = QmodTPCl_EP * std::sin (2 * psiTPCl_EP);
276276
277+ float QmodTPCr_EP = collision.qTPCR ();
277278 float psiTPCr_EP = collision.psiTPCR ();
278- float QxTPCr_EP = std::cos (2 * psiTPCr_EP);
279- float QyTPCr_EP = std::sin (2 * psiTPCr_EP);
280- float QmodTPCr_EP = std::hypot (QxTPCr_EP, QyTPCr_EP);
279+ float QxTPCr_EP = QmodTPCr_EP * std::cos (2 * psiTPCr_EP);
280+ float QyTPCr_EP = QmodTPCr_EP * std::sin (2 * psiTPCr_EP);
281281
282+ float QmodTPC_EP = collision.qTPC ();
282283 float psiTPC_EP = collision.psiTPC ();
283- float QxTPC_EP = std::cos (2 * psiTPC_EP);
284- float QyTPC_EP = std::sin (2 * psiTPC_EP);
285- float QmodTPC_EP = std::hypot (QxTPC_EP, QyTPC_EP);
284+ float QxTPC_EP = QmodTPC_EP * std::cos (2 * psiTPC_EP);
285+ float QyTPC_EP = QmodTPC_EP * std::sin (2 * psiTPC_EP);
286286
287287 // Qvec method
288288 float QxFT0A_Qvec = collision.qvecFT0ARe ();
0 commit comments