Skip to content
This repository has been archived by the owner on Jan 13, 2021. It is now read-only.

Commit

Permalink
Merge pull request #36 from usarica/master
Browse files Browse the repository at this point in the history 
Patch HJJ ME per-parton scenario
  • Loading branch information
@usarica
usarica authored Feb 28, 2018
2 parents a4ab4fa + 5346247 commit 6a537d6
Show file tree
Hide file tree
Showing 6 changed files with 781 additions and 765 deletions.
56 changes: 36 additions & 20 deletions MELA/src/TUtil.cc
View file
Original file line number Diff line number Diff line change
Expand Up @@ -4764,7 +4764,7 @@ double TUtil::HJJMatEl(
){
__modhiggsjj_MOD_evalamp_sbfh_unsymm_sa_select_exact(p4, &isel, &jsel, &ssel, &rsel, &msq_tmp);
MatElsq[jsel+5][isel+5] += msq_tmp; // Assign only those that match gen. info, if present at all.
if (verbosity >= TVar::DEBUG_VERBOSE) MELAout << "Channel (isel, jsel, rsel, ssel)=" << isel << ", " << jsel << ", " << rsel << ", " << ssel << '\t' << msq_tmp << endl;
if (verbosity >= TVar::DEBUG_VERBOSE) MELAout << "Channel (isel, jsel, rsel, ssel)=" << isel << ", " << jsel << ", " << ssel << ", " << rsel << '\t' << msq_tmp << endl;
}
}
else{ // gg->gg
Expand Down Expand Up @@ -4797,7 +4797,7 @@ double TUtil::HJJMatEl(
){
__modhiggsjj_MOD_evalamp_sbfh_unsymm_sa_select_exact(p4, &isel, &jsel, &ssel, &rsel, &msq_tmp);
MatElsq[jsel+5][isel+5] += msq_tmp; // Assign only those that match gen. info, if present at all.
if (verbosity >= TVar::DEBUG_VERBOSE) MELAout << "Channel (isel, jsel, rsel, ssel)=" << isel << ", " << jsel << ", " << rsel << ", " << ssel << '\t' << msq_tmp << endl;
if (verbosity >= TVar::DEBUG_VERBOSE) MELAout << "Channel (isel, jsel, rsel, ssel)=" << isel << ", " << jsel << ", " << ssel << ", " << rsel << '\t' << msq_tmp << endl;
}
}
else if ((isel>0 && jsel<0) || (isel<0 && jsel>0)){ // qQb/qbQ->?
Expand All @@ -4818,22 +4818,30 @@ double TUtil::HJJMatEl(
else{ rsel=2; ssel=-2; }
// The amplitude is aready multiplied by nf-1, so no need to calculate everything (nf-1) times.
if (
(partonIsUnknown[2] || MYIDUP_tmp[2]==rsel)
&&
(partonIsUnknown[3] || MYIDUP_tmp[3]==ssel)
(partonIsUnknown[2] && partonIsUnknown[3])
||
(partonIsUnknown[2] && std::abs(MYIDUP_tmp[3])!=std::abs(isel) && MYIDUP_tmp[3]<0)
||
(std::abs(MYIDUP_tmp[2])!=std::abs(isel) && MYIDUP_tmp[2]>0 && partonIsUnknown[3])
||
(std::abs(MYIDUP_tmp[2])!=std::abs(isel) && MYIDUP_tmp[2]==-MYIDUP_tmp[3] && MYIDUP_tmp[2]>0)
){
__modhiggsjj_MOD_evalamp_sbfh_unsymm_sa_select_exact(p4, &isel, &jsel, &rsel, &ssel, &msq_tmp);
MatElsq[jsel+5][isel+5] += msq_tmp; // Assign only those that match gen. info, if present at all.
if (verbosity >= TVar::DEBUG_VERBOSE) MELAout << "Channel (isel, jsel, rsel, ssel)=" << isel << ", " << jsel << ", " << rsel << ", " << ssel << '\t' << msq_tmp << endl;
}
if (
(partonIsUnknown[2] || MYIDUP_tmp[2]==ssel)
&&
(partonIsUnknown[3] || MYIDUP_tmp[3]==rsel)
(partonIsUnknown[2] && partonIsUnknown[3])
||
(partonIsUnknown[2] && std::abs(MYIDUP_tmp[3])!=std::abs(isel) && MYIDUP_tmp[3]>0)
||
(std::abs(MYIDUP_tmp[2])!=std::abs(isel) && MYIDUP_tmp[2]<0 && partonIsUnknown[3])
||
(std::abs(MYIDUP_tmp[2])!=std::abs(isel) && MYIDUP_tmp[2]==-MYIDUP_tmp[3] && MYIDUP_tmp[2]<0)
){
__modhiggsjj_MOD_evalamp_sbfh_unsymm_sa_select_exact(p4, &isel, &jsel, &ssel, &rsel, &msq_tmp);
MatElsq[jsel+5][isel+5] += msq_tmp; // Assign only those that match gen. info, if present at all.
if (verbosity >= TVar::DEBUG_VERBOSE) MELAout << "Channel (isel, jsel, rsel, ssel)=" << isel << ", " << jsel << ", " << rsel << ", " << ssel << '\t' << msq_tmp << endl;
if (verbosity >= TVar::DEBUG_VERBOSE) MELAout << "Channel (isel, jsel, rsel, ssel)=" << isel << ", " << jsel << ", " << ssel << ", " << rsel << '\t' << msq_tmp << endl;
}
}
else{ // qQb/qbQ->qQb/qbQ
Expand Down Expand Up @@ -4876,7 +4884,7 @@ double TUtil::HJJMatEl(
){
__modhiggsjj_MOD_evalamp_sbfh_unsymm_sa_select_exact(p4, &isel, &jsel, &ssel, &rsel, &msq_tmp);
MatElsq[jsel+5][isel+5] += msq_tmp; // Assign only those that match gen. info, if present at all.
if (verbosity >= TVar::DEBUG_VERBOSE) MELAout << "Channel (isel, jsel, rsel, ssel)=" << isel << ", " << jsel << ", " << rsel << ", " << ssel << '\t' << msq_tmp << endl;
if (verbosity >= TVar::DEBUG_VERBOSE) MELAout << "Channel (isel, jsel, rsel, ssel)=" << isel << ", " << jsel << ", " << ssel << ", " << rsel << '\t' << msq_tmp << endl;
}
}
} // End unswapped isel>=jsel cases
Expand Down Expand Up @@ -4913,7 +4921,7 @@ double TUtil::HJJMatEl(
){
__modhiggsjj_MOD_evalamp_sbfh_unsymm_sa_select_exact(p4, &isel, &jsel, &ssel, &rsel, &msq_tmp);
MatElsq[jsel+5][isel+5] += msq_tmp; // Assign only those that match gen. info, if present at all.
if (verbosity >= TVar::DEBUG_VERBOSE) MELAout << "Channel (isel, jsel, rsel, ssel)=" << isel << ", " << jsel << ", " << rsel << ", " << ssel << '\t' << msq_tmp << endl;
if (verbosity >= TVar::DEBUG_VERBOSE) MELAout << "Channel (isel, jsel, rsel, ssel)=" << isel << ", " << jsel << ", " << ssel << ", " << rsel << '\t' << msq_tmp << endl;
}
}
else if ((isel>0 && jsel<0) || (isel<0 && jsel>0)){ // qQb/qbQ->?
Expand All @@ -4934,22 +4942,30 @@ double TUtil::HJJMatEl(
else{ rsel=2; ssel=-2; }
// The amplitude is aready multiplied by nf-1, so no need to calculate everything (nf-1) times.
if (
(partonIsUnknown[2] || MYIDUP_tmp[2]==rsel)
&&
(partonIsUnknown[3] || MYIDUP_tmp[3]==ssel)
(partonIsUnknown[2] && partonIsUnknown[3])
||
(partonIsUnknown[2] && std::abs(MYIDUP_tmp[3])!=std::abs(isel) && MYIDUP_tmp[3]<0)
||
(std::abs(MYIDUP_tmp[2])!=std::abs(isel) && MYIDUP_tmp[2]>0 && partonIsUnknown[3])
||
(std::abs(MYIDUP_tmp[2])!=std::abs(isel) && MYIDUP_tmp[2]==-MYIDUP_tmp[3] && MYIDUP_tmp[2]>0)
){
__modhiggsjj_MOD_evalamp_sbfh_unsymm_sa_select_exact(p4, &isel, &jsel, &rsel, &ssel, &msq_tmp);
MatElsq[jsel+5][isel+5] += msq_tmp; // Assign only those that match gen. info, if present at all.
if (verbosity >= TVar::DEBUG_VERBOSE) MELAout << "Channel (isel, jsel, rsel, ssel)=" << isel << ", " << jsel << ", " << rsel << ", " << ssel << '\t' << msq_tmp << endl;
}
if (
(partonIsUnknown[2] || MYIDUP_tmp[2]==ssel)
&&
(partonIsUnknown[3] || MYIDUP_tmp[3]==rsel)
(partonIsUnknown[2] && partonIsUnknown[3])
||
(partonIsUnknown[2] && std::abs(MYIDUP_tmp[3])!=std::abs(isel) && MYIDUP_tmp[3]>0)
||
(std::abs(MYIDUP_tmp[2])!=std::abs(isel) && MYIDUP_tmp[2]<0 && partonIsUnknown[3])
||
(std::abs(MYIDUP_tmp[2])!=std::abs(isel) && MYIDUP_tmp[2]==-MYIDUP_tmp[3] && MYIDUP_tmp[2]<0)
){
__modhiggsjj_MOD_evalamp_sbfh_unsymm_sa_select_exact(p4, &isel, &jsel, &ssel, &rsel, &msq_tmp);
MatElsq[jsel+5][isel+5] += msq_tmp; // Assign only those that match gen. info, if present at all.
if (verbosity >= TVar::DEBUG_VERBOSE) MELAout << "Channel (isel, jsel, rsel, ssel)=" << isel << ", " << jsel << ", " << rsel << ", " << ssel << '\t' << msq_tmp << endl;
if (verbosity >= TVar::DEBUG_VERBOSE) MELAout << "Channel (isel, jsel, rsel, ssel)=" << isel << ", " << jsel << ", " << ssel << ", " << rsel << '\t' << msq_tmp << endl;
}
}
else{ // qQb/qbQ->qQb/qbQ
Expand All @@ -4969,7 +4985,7 @@ double TUtil::HJJMatEl(
){
__modhiggsjj_MOD_evalamp_sbfh_unsymm_sa_select_exact(p4, &isel, &jsel, &ssel, &rsel, &msq_tmp);
MatElsq[jsel+5][isel+5] += msq_tmp; // Assign only those that match gen. info, if present at all.
if (verbosity >= TVar::DEBUG_VERBOSE) MELAout << "Channel (isel, jsel, rsel, ssel)=" << isel << ", " << jsel << ", " << rsel << ", " << ssel << '\t' << msq_tmp << endl;
if (verbosity >= TVar::DEBUG_VERBOSE) MELAout << "Channel (isel, jsel, rsel, ssel)=" << isel << ", " << jsel << ", " << ssel << ", " << rsel << '\t' << msq_tmp << endl;
}
}
}
Expand All @@ -4992,7 +5008,7 @@ double TUtil::HJJMatEl(
){
__modhiggsjj_MOD_evalamp_sbfh_unsymm_sa_select_exact(p4, &isel, &jsel, &ssel, &rsel, &msq_tmp);
MatElsq[jsel+5][isel+5] += msq_tmp; // Assign only those that match gen. info, if present at all.
if (verbosity >= TVar::DEBUG_VERBOSE) MELAout << "Channel (isel, jsel, rsel, ssel)=" << isel << ", " << jsel << ", " << rsel << ", " << ssel << '\t' << msq_tmp << endl;
if (verbosity >= TVar::DEBUG_VERBOSE) MELAout << "Channel (isel, jsel, rsel, ssel)=" << isel << ", " << jsel << ", " << ssel << ", " << rsel << '\t' << msq_tmp << endl;
}
}
} // End swapped isel<jsel cases
Expand Down
Loading

0 comments on commit 6a537d6

Please sign in to comment.