|
18 | 18 | #include "PWGEM/Dilepton/Utils/EMFwdTrack.h" |
19 | 19 | #include "PWGEM/Dilepton/Utils/EMTrack.h" |
20 | 20 | #include "PWGEM/Dilepton/Utils/EMTrackUtilities.h" |
21 | | -// #include "PWGEM/Dilepton/Utils/EventHistograms.h" |
22 | 21 | #include "PWGEM/Dilepton/Utils/EventMixingHandler.h" |
23 | 22 | #include "PWGEM/Dilepton/Utils/PairUtilities.h" |
24 | 23 |
|
@@ -111,8 +110,8 @@ struct DileptonPolarization { |
111 | 110 |
|
112 | 111 | struct : ConfigurableGroup { |
113 | 112 | std::string prefix = "dileptoncut_group"; |
114 | | - Configurable<float> cfg_min_mass{"cfg_min_mass", 0.0, "min mass"}; |
115 | | - Configurable<float> cfg_max_mass{"cfg_max_mass", 1e+10, "max mass"}; |
| 113 | + Configurable<float> cfg_min_pair_mass{"cfg_min_pair_mass", 0.0, "min pair mass"}; |
| 114 | + Configurable<float> cfg_max_pair_mass{"cfg_max_pair_mass", 1e+10, "max pair mass"}; |
116 | 115 | Configurable<float> cfg_min_pair_pt{"cfg_min_pair_pt", 0.0, "min pair pT"}; |
117 | 116 | Configurable<float> cfg_max_pair_pt{"cfg_max_pair_pt", 1e+10, "max pair pT"}; |
118 | 117 | Configurable<float> cfg_min_pair_y{"cfg_min_pair_y", -0.9, "min pair rapidity"}; |
@@ -464,6 +463,14 @@ struct DileptonPolarization { |
464 | 463 | ROOT::Math::PtEtaPhiMVector v2(dilepton.pt2(), dilepton.eta2(), dilepton.phi2(), leptonM2); |
465 | 464 | ROOT::Math::PtEtaPhiMVector v12 = v1 + v2; |
466 | 465 |
|
| 466 | + if (v12.M() < dileptoncuts.cfg_min_pair_mass || dileptoncuts.cfg_max_pair_mass < v12.M()) { |
| 467 | + return false; |
| 468 | + } |
| 469 | + |
| 470 | + if (v12.Pt() < dileptoncuts.cfg_min_pair_pt || dileptoncuts.cfg_max_pair_pt < v12.Pt()) { |
| 471 | + return false; |
| 472 | + } |
| 473 | + |
467 | 474 | if (v12.Rapidity() < dileptoncuts.cfg_min_pair_y || dileptoncuts.cfg_max_pair_y < v12.Rapidity()) { |
468 | 475 | return false; |
469 | 476 | } |
|
0 commit comments