Skip to content

Fix assert in pptt code generation #340

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
valassi merged 6 commits into from
Jan 24, 2022
Merged

Fix assert in pptt code generation #340

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
6 commits
Select commit Hold shift + click to select a range
ace6b6a
[pptt] add pptt to code generation - I imediately hit Olivier's bug #337
valassi Jan 20, 2022
517d1fd
[pptt] add uu~>tt~ and check that this is what fails in pp to tt (#337)
valassi Jan 24, 2022
2876844
[pptt] BUG FIX #337, remove assert(False) - this fixes uutt but not y…
valassi Jan 24, 2022
e9dbaeb
[pptt] (bug fix) remove a debug printout that I had forgotten enabled…
valassi Jan 24, 2022
38ac549
[pptt] minor bug fix in syncManu script (was incorrectly including li…
valassi Jan 24, 2022
bd2b9ab
[pptt] **COMPLETE PPTT ** regenerate auto and resync manu for all 3 p…
valassi Jan 24, 2022
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
1 change: 0 additions & 1 deletion ...cpp/CODEGEN/PLUGIN/CUDACPP_SA_OUTPUT/madgraph/iolibs/template_files/gpu/EventStatistics.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -139,7 +139,6 @@ namespace mg5amcCpu
<< s.tag << "StdDevWeight = ( " << s.stdWG()
<< std::string(16, ' ') << " )" << std::endl
<< std::defaultfloat; // default format: affects all floats
out << s.tag << "SqsMatrixElemValue = " << s.sqsMEdiff << std::endl;
}
};

Expand Down
3 changes: 2 additions & 1 deletion epochX/cudacpp/CODEGEN/PLUGIN/CUDACPP_SA_OUTPUT/model_handling.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -1010,6 +1010,7 @@ def write_process_cc_file(self, writer):
return out

# AV - replace the export_cpp.OneProcessExporterGPU method (improve formatting? actually keep all defaults!)
# [NB this is used in uu~>tt~ but not in gg>tt~ or e+e->mu+mu-, see issue #337]
@staticmethod
def coeff(ff_number, frac, is_imaginary, Nc_power, Nc_value=3):
"""Returns a nicely formatted string for the coefficients in JAMP lines"""
Expand All @@ -1024,7 +1025,7 @@ def coeff(ff_number, frac, is_imaginary, Nc_power, Nc_value=3):
return '-cxtype(0,1)*' # AV keep default (this is not used in eemumu - should use cI eventually)
else:
return '-' # AV keep default (eg jamp_sv[0] += -amp_sv[0])
assert(False)
###assert(False) # [this had been inserted to check if coeff is used at all, it is used in uu~>tt~, see #337]
res_str = '%+i.' % total_coeff.numerator
if total_coeff.denominator != 1:
# Check if total_coeff is an integer
Expand Down
9 changes: 6 additions & 3 deletions epochX/cudacpp/CODEGEN/generateAndCompare.sh
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -14,12 +14,15 @@ function codeGenAndDiff()
gg_tt)
cmd="generate g g > t t~"
;;
###gg_ttg)
### cmd="generate g g > t t~ g"
### ;;
gg_ttgg)
cmd="generate g g > t t~ g g"
;;
pp_tt)
cmd="generate p p > t t~"
;;
uu_tt)
cmd="generate u u~ > t t~"
;;
*)
echo -e "\nWARNING! Skipping unknown process '$proc'"
return
Expand Down
2 changes: 1 addition & 1 deletion epochX/cudacpp/CODEGEN/syncManu.sh
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -38,7 +38,7 @@ cd $(dirname $0)
for proc in ${processes}; do
cmdfile=$(mktemp)
echo "------------------------------------------------------------------"
./diffCode.sh ../${proc}.auto ../${proc} -r | grep diff | awk '{print "cp -dp", $(NF-1), $NF}' > ${cmdfile}
./diffCode.sh ../${proc}.auto ../${proc} -r | grep ^diff | awk '{print "cp -dp", $(NF-1), $NF}' > ${cmdfile}
cat ${cmdfile}
source ${cmdfile}
echo -e "\nPENDING DIFFERENCES (before copying to manual any new files added in auto):"
Expand Down
9 changes: 5 additions & 4 deletions epochX/cudacpp/ee_mumu.auto/CODEGEN_cudacpp_ee_mumu_log.txt
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,4 +1,5 @@
Running MG5 in debug mode
('WARNING: loading of madgraph too slow!!!', 0.8647480010986328)
************************************************************
* *
* W E L C O M E to *
Expand Down Expand Up @@ -50,7 +51,7 @@ generate e+ e- > mu+ mu-
No model currently active, so we import the Standard Model
INFO: load particles
INFO: load vertices
DEBUG: model prefixing takes 0.006859779357910156 
DEBUG: model prefixing takes 0.006945133209228516 
INFO: Restrict model sm with file models/sm/restrict_default.dat .
INFO: Change particles name to pass to MG5 convention
DEBUG: True [misc.py at line 2192] 
Expand Down Expand Up @@ -120,6 +121,6 @@ INFO: /data/avalassi/GPU2020/MG5aMC/2.7.0_gpu/CODEGEN_cudacpp_ee_mumu/src/. and
DEBUG: Entering PLUGIN_ProcessExporter.finalize [output.py at line 178] 
quit

real 0m4.038s
user 0m0.930s
sys 0m0.132s
real 0m5.407s
user 0m0.987s
sys 0m0.140s
1 change: 0 additions & 1 deletion epochX/cudacpp/ee_mumu.auto/SubProcesses/EventStatistics.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -139,7 +139,6 @@ namespace mg5amcCpu
<< s.tag << "StdDevWeight = ( " << s.stdWG()
<< std::string(16, ' ') << " )" << std::endl
<< std::defaultfloat; // default format: affects all floats
out << s.tag << "SqsMatrixElemValue = " << s.sqsMEdiff << std::endl;
}
};

Expand Down
1 change: 0 additions & 1 deletion epochX/cudacpp/ee_mumu/SubProcesses/EventStatistics.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -139,7 +139,6 @@ namespace mg5amcCpu
<< s.tag << "StdDevWeight = ( " << s.stdWG()
<< std::string(16, ' ') << " )" << std::endl
<< std::defaultfloat; // default format: affects all floats
out << s.tag << "SqsMatrixElemValue = " << s.sqsMEdiff << std::endl;
}
};

Expand Down
21 changes: 11 additions & 10 deletions epochX/cudacpp/gg_tt.auto/CODEGEN_cudacpp_gg_tt_log.txt
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,4 +1,5 @@
Running MG5 in debug mode
('WARNING: loading of madgraph too slow!!!', 1.7634224891662598)
************************************************************
* *
* W E L C O M E to *
Expand Down Expand Up @@ -50,7 +51,7 @@ generate g g > t t~
No model currently active, so we import the Standard Model
INFO: load particles
INFO: load vertices
DEBUG: model prefixing takes 0.006892681121826172 
DEBUG: model prefixing takes 0.007020235061645508 
INFO: Restrict model sm with file models/sm/restrict_default.dat .
INFO: Change particles name to pass to MG5 convention
DEBUG: True [misc.py at line 2192] 
Expand All @@ -67,7 +68,7 @@ INFO: Please specify coupling orders to bypass this step.
INFO: Trying coupling order WEIGHTED<=2: WEIGTHED IS QCD+2*QED
INFO: Trying process: g g > t t~ WEIGHTED<=2 @1
INFO: Process has 3 diagrams
1 processes with 3 diagrams generated in 0.010 s
1 processes with 3 diagrams generated in 0.077 s
Total: 1 processes with 3 diagrams
output standalone_cudacpp CODEGEN_cudacpp_gg_tt
Output will be done with PLUGIN: CUDACPP_SA_OUTPUT
Expand All @@ -87,16 +88,16 @@ FileWriter <class 'PLUGIN.CUDACPP_SA_OUTPUT.model_handling.PLUGIN_FileWriter'> f
DEBUG: Entering PLUGIN_OneProcessExporter.write_process_h_file [model_handling.py at line 983] 
FileWriter <class 'PLUGIN.CUDACPP_SA_OUTPUT.model_handling.PLUGIN_FileWriter'> for /data/avalassi/GPU2020/MG5aMC/2.7.0_gpu/CODEGEN_cudacpp_gg_tt/SubProcesses/P1_Sigma_sm_gg_ttx/./CPPProcess.cc
DEBUG: Entering PLUGIN_OneProcessExporter.write_process_cc_file [model_handling.py at line 1005] 
DEBUG: call =  vxxxxx( momenta,m_pars->%s, cHel[ihel][%d],%+d, w_sv[%d], %d ); [model_handling.py at line 1289] 
DEBUG: ('ZERO', 0, -1, 0, 0) [model_handling.py at line 1290] 
DEBUG: call =  vxxxxx( momenta,m_pars->%s, cHel[ihel][%d],%+d, w_sv[%d], %d ); [model_handling.py at line 1289] 
DEBUG: ('ZERO', 1, -1, 1, 1) [model_handling.py at line 1290] 
DEBUG: call =  vxxxxx( momenta,m_pars->%s, cHel[ihel][%d],%+d, w_sv[%d], %d ); [model_handling.py at line 1290] 
DEBUG: ('ZERO', 0, -1, 0, 0) [model_handling.py at line 1291] 
DEBUG: call =  vxxxxx( momenta,m_pars->%s, cHel[ihel][%d],%+d, w_sv[%d], %d ); [model_handling.py at line 1290] 
DEBUG: ('ZERO', 1, -1, 1, 1) [model_handling.py at line 1291] 
DEBUG: only one Matrix-element supported? 
INFO: Created files CPPProcess.h and CPPProcess.cc in directory /data/avalassi/GPU2020/MG5aMC/2.7.0_gpu/CODEGEN_cudacpp_gg_tt/SubProcesses/P1_Sigma_sm_gg_ttx/.
DEBUG: Entering PLUGIN_OneProcessExporter.edit_check_sa [model_handling.py at line 948] 
DEBUG: Entering PLUGIN_OneProcessExporter.edit_mgonGPU [model_handling.py at line 963] 
DEBUG: Entering PLUGIN_OneProcessExporter.edit_processidfile [model_handling.py at line 970] 
Generated helas calls for 1 subprocesses (3 diagrams) in 0.014 s
Generated helas calls for 1 subprocesses (3 diagrams) in 0.084 s
DEBUG: Entering PLUGIN_ProcessExporter.convert_model (create the model) [output.py at line 169] 
ALOHA: aloha creates routines (starting by VVV1 with options P0)
<class 'aloha.create_aloha.AbstractRoutine'> VVV1
Expand All @@ -117,6 +118,6 @@ INFO: /data/avalassi/GPU2020/MG5aMC/2.7.0_gpu/CODEGEN_cudacpp_gg_tt/src/. and /d
DEBUG: Entering PLUGIN_ProcessExporter.finalize [output.py at line 178] 
quit

real 0m4.133s
user 0m0.759s
sys 0m0.169s
real 0m13.913s
user 0m0.894s
sys 0m0.139s
1 change: 0 additions & 1 deletion epochX/cudacpp/gg_tt.auto/SubProcesses/EventStatistics.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -139,7 +139,6 @@ namespace mg5amcCpu
<< s.tag << "StdDevWeight = ( " << s.stdWG()
<< std::string(16, ' ') << " )" << std::endl
<< std::defaultfloat; // default format: affects all floats
out << s.tag << "SqsMatrixElemValue = " << s.sqsMEdiff << std::endl;
}
};

Expand Down
1 change: 0 additions & 1 deletion epochX/cudacpp/gg_tt/SubProcesses/EventStatistics.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -139,7 +139,6 @@ namespace mg5amcCpu
<< s.tag << "StdDevWeight = ( " << s.stdWG()
<< std::string(16, ' ') << " )" << std::endl
<< std::defaultfloat; // default format: affects all floats
out << s.tag << "SqsMatrixElemValue = " << s.sqsMEdiff << std::endl;
}
};

Expand Down
29 changes: 15 additions & 14 deletions epochX/cudacpp/gg_ttgg.auto/CODEGEN_cudacpp_gg_ttgg_log.txt
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,4 +1,5 @@
Running MG5 in debug mode
('WARNING: loading of madgraph too slow!!!', 1.7154383659362793)
************************************************************
* *
* W E L C O M E to *
Expand Down Expand Up @@ -50,7 +51,7 @@ generate g g > t t~ g g
No model currently active, so we import the Standard Model
INFO: load particles
INFO: load vertices
DEBUG: model prefixing takes 0.006979465484619141 
DEBUG: model prefixing takes 0.01733565330505371 
INFO: Restrict model sm with file models/sm/restrict_default.dat .
INFO: Change particles name to pass to MG5 convention
DEBUG: True [misc.py at line 2192] 
Expand All @@ -67,7 +68,7 @@ INFO: Please specify coupling orders to bypass this step.
INFO: Trying coupling order WEIGHTED<=4: WEIGTHED IS QCD+2*QED
INFO: Trying process: g g > t t~ g g WEIGHTED<=4 @1
INFO: Process has 123 diagrams
1 processes with 123 diagrams generated in 0.209 s
1 processes with 123 diagrams generated in 0.586 s
Total: 1 processes with 123 diagrams
output standalone_cudacpp CODEGEN_cudacpp_gg_ttgg
Output will be done with PLUGIN: CUDACPP_SA_OUTPUT
Expand All @@ -87,20 +88,20 @@ FileWriter <class 'PLUGIN.CUDACPP_SA_OUTPUT.model_handling.PLUGIN_FileWriter'> f
DEBUG: Entering PLUGIN_OneProcessExporter.write_process_h_file [model_handling.py at line 983] 
FileWriter <class 'PLUGIN.CUDACPP_SA_OUTPUT.model_handling.PLUGIN_FileWriter'> for /data/avalassi/GPU2020/MG5aMC/2.7.0_gpu/CODEGEN_cudacpp_gg_ttgg/SubProcesses/P1_Sigma_sm_gg_ttxgg/./CPPProcess.cc
DEBUG: Entering PLUGIN_OneProcessExporter.write_process_cc_file [model_handling.py at line 1005] 
DEBUG: call =  vxxxxx( momenta,m_pars->%s, cHel[ihel][%d],%+d, w_sv[%d], %d ); [model_handling.py at line 1289] 
DEBUG: ('ZERO', 0, -1, 0, 0) [model_handling.py at line 1290] 
DEBUG: call =  vxxxxx( momenta,m_pars->%s, cHel[ihel][%d],%+d, w_sv[%d], %d ); [model_handling.py at line 1289] 
DEBUG: ('ZERO', 1, -1, 1, 1) [model_handling.py at line 1290] 
DEBUG: call =  vxxxxx( momenta,m_pars->%s, cHel[ihel][%d],%+d, w_sv[%d], %d ); [model_handling.py at line 1289] 
DEBUG: ('ZERO', 4, 1, 4, 4) [model_handling.py at line 1290] 
DEBUG: call =  vxxxxx( momenta,m_pars->%s, cHel[ihel][%d],%+d, w_sv[%d], %d ); [model_handling.py at line 1289] 
DEBUG: ('ZERO', 5, 1, 5, 5) [model_handling.py at line 1290] 
DEBUG: call =  vxxxxx( momenta,m_pars->%s, cHel[ihel][%d],%+d, w_sv[%d], %d ); [model_handling.py at line 1290] 
DEBUG: ('ZERO', 0, -1, 0, 0) [model_handling.py at line 1291] 
DEBUG: call =  vxxxxx( momenta,m_pars->%s, cHel[ihel][%d],%+d, w_sv[%d], %d ); [model_handling.py at line 1290] 
DEBUG: ('ZERO', 1, -1, 1, 1) [model_handling.py at line 1291] 
DEBUG: call =  vxxxxx( momenta,m_pars->%s, cHel[ihel][%d],%+d, w_sv[%d], %d ); [model_handling.py at line 1290] 
DEBUG: ('ZERO', 4, 1, 4, 4) [model_handling.py at line 1291] 
DEBUG: call =  vxxxxx( momenta,m_pars->%s, cHel[ihel][%d],%+d, w_sv[%d], %d ); [model_handling.py at line 1290] 
DEBUG: ('ZERO', 5, 1, 5, 5) [model_handling.py at line 1291] 
DEBUG: only one Matrix-element supported? 
INFO: Created files CPPProcess.h and CPPProcess.cc in directory /data/avalassi/GPU2020/MG5aMC/2.7.0_gpu/CODEGEN_cudacpp_gg_ttgg/SubProcesses/P1_Sigma_sm_gg_ttxgg/.
DEBUG: Entering PLUGIN_OneProcessExporter.edit_check_sa [model_handling.py at line 948] 
DEBUG: Entering PLUGIN_OneProcessExporter.edit_mgonGPU [model_handling.py at line 963] 
DEBUG: Entering PLUGIN_OneProcessExporter.edit_processidfile [model_handling.py at line 970] 
Generated helas calls for 1 subprocesses (123 diagrams) in 0.562 s
Generated helas calls for 1 subprocesses (123 diagrams) in 1.149 s
DEBUG: Entering PLUGIN_ProcessExporter.convert_model (create the model) [output.py at line 169] 
ALOHA: aloha creates routines (starting by VVV1)
<class 'aloha.create_aloha.AbstractRoutine'> VVV1
Expand Down Expand Up @@ -137,6 +138,6 @@ INFO: /data/avalassi/GPU2020/MG5aMC/2.7.0_gpu/CODEGEN_cudacpp_gg_ttgg/src/. and
DEBUG: Entering PLUGIN_ProcessExporter.finalize [output.py at line 178] 
quit

real 0m12.617s
user 0m2.010s
sys 0m0.131s
real 0m10.858s
user 0m2.123s
sys 0m0.153s
1 change: 0 additions & 1 deletion epochX/cudacpp/gg_ttgg.auto/SubProcesses/EventStatistics.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -139,7 +139,6 @@ namespace mg5amcCpu
<< s.tag << "StdDevWeight = ( " << s.stdWG()
<< std::string(16, ' ') << " )" << std::endl
<< std::defaultfloat; // default format: affects all floats
out << s.tag << "SqsMatrixElemValue = " << s.sqsMEdiff << std::endl;
}
};

Expand Down
1 change: 0 additions & 1 deletion epochX/cudacpp/gg_ttgg/SubProcesses/EventStatistics.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -139,7 +139,6 @@ namespace mg5amcCpu
<< s.tag << "StdDevWeight = ( " << s.stdWG()
<< std::string(16, ' ') << " )" << std::endl
<< std::defaultfloat; // default format: affects all floats
out << s.tag << "SqsMatrixElemValue = " << s.sqsMEdiff << std::endl;
}
};

Expand Down