11 #include <analysis/ContinuumSuppression/CleoCones.h>
21 const std::vector<TVector3>& p3_cms_roe,
22 const TVector3& thrustB,
23 bool calc_CleoCones_with_all,
24 bool calc_CleoCones_with_roe
33 if (calc_CleoCones_with_all ==
true) {
34 for (
int i = 1; i <= 9; i++) {
35 float momentum_flow_all = 0;
36 for (
auto& iter0 : p3_cms_all) {
50 float angle = ((180 * (thrustB.Angle(iter0))) / M_PI);
51 if (((((i - 1) * 10) <= angle) && (angle < (i * 10))) || (((180 - (i * 10)) < angle) && (angle <= (180 - ((i - 1) * 10))))) {
52 momentum_flow_all += iter0.Mag();
55 if ((i == 9) && (angle == 90)) {
56 momentum_flow_all += iter0.Mag();
66 if (calc_CleoCones_with_roe ==
true) {
67 for (
int i = 1; i <= 9; i++) {
68 float momentum_flow_roe = 0;
69 for (
auto& iter1 : p3_cms_roe) {
70 float angle = ((180 * (thrustB.Angle(iter1))) / M_PI);
71 if (((((i - 1) * 10) <= angle) && (angle < (i * 10))) || (((180 - (i * 10)) < angle) && (angle <= (180 - ((i - 1) * 10))))) {
72 momentum_flow_roe += iter1.Mag();
74 if ((i == 9) && (angle == 90)) {
75 momentum_flow_roe += iter1.Mag();