std::vector<TruthVariables> TruthLevel::Multiplets(Event const& event, PreCuts const& pre_cuts, TMVA::Reader const&)
{
INFO0;
return Jets(event, pre_cuts, [&](Particle & jet) {
// if (Problematic(jet, pre_cuts)) throw boca::Problematic();
return jet;
});
}
int TruthLevel::Train(Event const& event, PreCuts const& pre_cuts, Tag)
{
INFO0;
return SaveEntries(Jets(event, pre_cuts, [&](Particle & jet) {
// if (Problematic(jet, pre_cuts, tag)) throw boca::Problematic();
return jet;
}));
}
bool Jets(Cluster_Amplitude *ampl,int mode)
{
DEBUG_FUNC("mode = "<<mode);
msg_Debugging()<<*ampl<<"\n";
PHASIC::Jet_Finder *jf(ampl->JF<PHASIC::Jet_Finder>());
double q2cut(jf->Ycut()*sqr(rpa->gen.Ecms()));
double q2min(std::numeric_limits<double>::max());
size_t imin(0), jmin(0), kmin(0);
Flavour mofl;
for (size_t i(0);i<ampl->Legs().size();++i) {
Cluster_Leg *li(ampl->Leg(i));
for (size_t j(Max(i+1,ampl->NIn()));j<ampl->Legs().size();++j) {
Cluster_Leg *lj(ampl->Leg(j));
if (j<ampl->NIn()) continue;
for (size_t k(0);k<ampl->Legs().size();++k) {
if (k==i || k==j) continue;
Cluster_Leg *lk(ampl->Leg(k));
if (i<ampl->NIn() && k>=ampl->NIn()) continue;
if (lk->Flav().Strong() &&
li->Flav().Strong() && lj->Flav().Strong()) {
if (i<ampl->NIn()) li->SetMom(-li->Mom());
if (k<ampl->NIn()) lk->SetMom(-lk->Mom());
double q2ijk(pT2pythia(ampl,*li,*lj,*lk,i<ampl->NIn()?-1:1));
msg_Debugging()<<"Q_{"<<ID(li->Id())<<ID(lj->Id())
<<","<<ID(lk->Id())<<"} = "<<sqrt(q2ijk)<<"\n";
if (i<ampl->NIn()) li->SetMom(-li->Mom());
if (k<ampl->NIn()) lk->SetMom(-lk->Mom());
if (mode==0) {
if (q2ijk<q2cut) return false;
}
else {
if (q2ijk<q2min) {
q2min=q2ijk;
mofl=Flavour(kf_gluon);
if (li->Flav().IsGluon()) mofl=lj->Flav();
if (lj->Flav().IsGluon()) mofl=li->Flav();
imin=i;
jmin=j;
kmin=k;
}
}
}
}
}
}
if (mode!=0 && imin!=jmin) {
Vec4D_Vector p=Combine(*ampl,imin,jmin,kmin,mofl);
if (p.empty()) {
msg_Error()<<METHOD<<"(): Combine failed. Use R configuration."<<std::endl;
return Jets(ampl,0);
}
Cluster_Amplitude *bampl(Cluster_Amplitude::New());
bampl->SetProc(ampl->Proc<void>());
bampl->SetMS(ampl->MS());
bampl->SetNIn(ampl->NIn());
bampl->SetJF(ampl->JF<void>());
for (int i(0), j(0);i<ampl->Legs().size();++i) {
if (i==jmin) continue;
if (i==imin) {
bampl->CreateLeg(p[j],mofl,ampl->Leg(i)->Col());
bampl->Legs().back()->SetId(ampl->Leg(imin)->Id()|ampl->Leg(jmin)->Id());
bampl->Legs().back()->SetK(ampl->Leg(kmin)->Id());
}
else {
bampl->CreateLeg(p[j],ampl->Leg(i)->Flav(),ampl->Leg(i)->Col());
}
++j;
}
bool res=Jets(bampl,0);
bampl->Delete();
return res;
}
msg_Debugging()<<"--- Jet veto ---\n";
return true;
}
