PWGJE: give a maximum-pT cut on jet constituents

[joonsukbae]

Hi Nima,
Existing track pT max cut can't allow jet rejection after jets finding w/o track pT cut both in jetfinderQA nor jetfinder tasks, thus I added this selection in jet acceptance task.

[nzardosh]

Hi Joonsuk, i think the function is cleaner like this:

template <typename T, typename U>
bool isAcceptedJet(U const& jet)
{

if (jetAreaFractionMin > -98.0) {
  if (jet.area() < jetAreaFractionMin * M_PI * (jet.r() / 100.0) * (jet.r() / 100.0)) {
    return false;
  }
}
bool checkConstituentPt = true;
bool checkConstituentMinPt = (leadingConstituentPtMin > -98.0);
bool checkConstituentMaxPt = (leadingConstituentPtMax < 9998.0);
if (!checkConstituentMinPt && !checkConstituentMaxPt) {
  checkConstituentPt = false;
}

if (checkConstituentPt){
  bool isMinLeadingConstituent = !checkConstituentMinPt;
  bool isMaxLeadingConstituent = true;

for (const auto& constituent : jet.template tracks_as<T>()) {
  double pt = constituent.pt();

  if (checkConstituentMinPt && pt >= leadingConstituentPtMin) {
    isMinLeadingConstituent = true;
  }
  if (checkConstituentMaxPt && pt > leadingConstituentPtMax) {
    isMaxLeadingConstituent = false;
  }
}
return isMinLeadingConstituent && isMaxLeadingConstituent;
}

return true;

}

can you please check it works and change accoridngly?

[joonsukbae]

Sure, thanks for the repair. I'll check it then fix it:)
Btw, would the lack of an early exit condition significantly impact performance? In the case of the MB MC dataset I’m using, the jets have an average of around 3 constituents each.

[nzardosh]

I think the problem with the early condition was that if we want to add a condition not related to the pT later then it would be skipped

[joonsukbae]

Please ignore this fetch at 21:23. Because my O2 installation failed in local, I pulled upstream.

[joonsukbae]

I checked your solution and it works very well in local