Kamikaze reweighting compilation fails
Dear all,
I am trying to use Kamikaze reweighting to reweight my NLO aTGC sample. I am using MG5_aMC_v2_3_3.
The NLO aTGC model is by Celine Degrande: EWdim6NLO. My process is NLO WZ with 0 and 1 jets:
set group_subprocesses False
import model EWdim6NLO-
generate p p > w+ l+ l- QED=3 [QCD] @0
add process p p > w- l+ l- QED=3 [QCD] @1
I successfully generate events, but I get an error while compiling reweighting step. I run generate_events with reweighting ON:
1 Perturbative order of the calculation: order=NLO
2 Fixed order (no event generation and no MC@[N]LO matching): fixed_order=OFF
3 Shower the generated events: shower=OFF
4 Decay particles with the MadSpin module: madspin=OFF
5 Add weights to the events based on changing model parameters: reweight=ON
And my output:
=======
....
INFO: Updating the number of unweighted events per channel
Intermediate results:
Random seed: 34
Total cross-section: 2.299e+00 +- 6.5e-03 pb
Total abs(cross-section): 2.720e+00 +- 6.6e-03 pb
INFO: Generating events
...
INFO: Doing reweight
INFO: Idle: 0, Running: 0, Completed: 23 [ current time: 06h16 ]
INFO: Collecting events
INFO:
Summary:
Process p p > w+ l+ l- QED=3 [QCD] @0 ; p p > w- l+ l- QED=3 [QCD] @1
Run at p-p collider (6500.0 + 6500.0 GeV)
Total cross-section: 2.286e+00 +- 6.5e-03 pb
Number of events generated: 1000
Parton shower to be used: PYTHIA8
Fraction of negative weights: 0.08
Total running time : 6m 48s
INFO: The /nfs_scratch/
INFO: Events generated
reweight -from_cards
INFO: Extracting the banner ...
WARNING: Remember that the reweighting is performed at Leading Order. NLO precision is not guarantee.
...
INFO: process: p p > w+ l+ l-
INFO: options: QED=3
INFO: Running Reweight
launch
INFO: detected model: EWdim6NLO-
INFO: generating the square matrix element for reweighting
INFO: generate p p > w+ l+ l- QED=3 @0 --no_warning=
INFO: Done 16.53
set Dim6 1 6
INFO: modify param_card information BLOCK dim6 with id (1,) set to 6.0
set Dim6 2 0
INFO: modify param_card information BLOCK dim6 with id (2,) set to 0.0
set Dim6 3 0
INFO: modify param_card information BLOCK dim6 with id (3,) set to 0.0
Do you want to edit a card (press enter to bypass editing)?
1 / param : param_card.dat
9 / plot : plot_card.dat
you can also
- enter the path to a valid card or banner.
- use the 'set' command to modify a parameter directly.
The set option works only for param_card and run_card.
Type 'help set' for more information on this command.
- call an external program (ASperGE/
Type 'help' for the list of available command
[0, done, 1, param, 9, plot, enter path]
The answer to the previous question is not set in your input file
Use 0 value
INFO: starts to compute weight for events with the following modification to the param_card:
INFO: set param_card dim6 1 6.0 # orig: 4.0
set param_card dim6 2 0.0 # orig: 3.0
set param_card dim6 3 0.0 # orig: 15.0
INFO: Event nb 0 current time: 06h16
Command "launch " interrupted with error:
MadGraph5Error: A compilation Error occurs when trying to compile /nfs_scratch/
The compilation fails with the following output message:
touch __init__.py
f2py -L../../lib/ -ldhelas -lmodel -c matrix.f -m matrix2py --fcompiler=
Unknown vendor: "gfortran"
running build
running config_cc
unifing config_cc, config, build_clib, build_ext, build commands --compiler options
running config_fc
unifing config_fc, config, build_clib, build_ext, build commands --fcompiler options
running build_src
build_src
building extension "matrix2py" sources
f2py options: []
f2py:> /tmp/tmpWGXxcc/
creating /tmp/tmpWGXxcc
creating /tmp/tmpWGXxcc/
Reading fortran codes...
Reading file 'matrix.f' (format:fix,strict)
Line #85 in matrix.f:" DATA (NHEL(I, 1),I=1,5) / 1,-1,-1, 1,-1/"
analyzeline: implied-DO list "(nhel(i, 1),i=1,5)" is not supported. Skipping.
Line #86 in matrix.f:" DATA (NHEL(I, 2),I=1,5) / 1,-1,-1, 1, 1/"
analyzeline: implied-DO list "(nhel(i, 2),i=1,5)" is not supported. Skipping.
Line #87 in matrix.f:" DATA (NHEL(I, 3),I=1,5) / 1,-1,-1,-1,-1/"
analyzeline: implied-DO list "(nhel(i, 3),i=1,5)" is not supported. Skipping.
Line #88 in matrix.f:" DATA (NHEL(I, 4),I=1,5) / 1,-1,-1,-1, 1/"
analyzeline: implied-DO list "(nhel(i, 4),i=1,5)" is not supported. Skipping.
Line #89 in matrix.f:" DATA (NHEL(I, 5),I=1,5) / 1,-1, 0, 1,-1/"
analyzeline: implied-DO list "(nhel(i, 5),i=1,5)" is not supported. Skipping.
Line #90 in matrix.f:" DATA (NHEL(I, 6),I=1,5) / 1,-1, 0, 1, 1/"
analyzeline: implied-DO list "(nhel(i, 6),i=1,5)" is not supported. Skipping.
Line #91 in matrix.f:" DATA (NHEL(I, 7),I=1,5) / 1,-1, 0,-1,-1/"
analyzeline: implied-DO list "(nhel(i, 7),i=1,5)" is not supported. Skipping.
Line #92 in matrix.f:" DATA (NHEL(I, 8),I=1,5) / 1,-1, 0,-1, 1/"
analyzeline: implied-DO list "(nhel(i, 8),i=1,5)" is not supported. Skipping.
Line #93 in matrix.f:" DATA (NHEL(I, 9),I=1,5) / 1,-1, 1, 1,-1/"
analyzeline: implied-DO list "(nhel(i, 9),i=1,5)" is not supported. Skipping.
Line #94 in matrix.f:" DATA (NHEL(I, 10),I=1,5) / 1,-1, 1, 1, 1/"
analyzeline: implied-DO list "(nhel(i, 10),i=1,5)" is not supported. Skipping.
Line #95 in matrix.f:" DATA (NHEL(I, 11),I=1,5) / 1,-1, 1,-1,-1/"
analyzeline: implied-DO list "(nhel(i, 11),i=1,5)" is not supported. Skipping.
Line #96 in matrix.f:" DATA (NHEL(I, 12),I=1,5) / 1,-1, 1,-1, 1/"
analyzeline: implied-DO list "(nhel(i, 12),i=1,5)" is not supported. Skipping.
Line #97 in matrix.f:" DATA (NHEL(I, 13),I=1,5) / 1, 1,-1, 1,-1/"
analyzeline: implied-DO list "(nhel(i, 13),i=1,5)" is not supported. Skipping.
Line #98 in matrix.f:" DATA (NHEL(I, 14),I=1,5) / 1, 1,-1, 1, 1/"
analyzeline: implied-DO list "(nhel(i, 14),i=1,5)" is not supported. Skipping.
Line #99 in matrix.f:" DATA (NHEL(I, 15),I=1,5) / 1, 1,-1,-1,-1/"
analyzeline: implied-DO list "(nhel(i, 15),i=1,5)" is not supported. Skipping.
Line #100 in matrix.f:" DATA (NHEL(I, 16),I=1,5) / 1, 1,-1,-1, 1/"
analyzeline: implied-DO list "(nhel(i, 16),i=1,5)" is not supported. Skipping.
Line #101 in matrix.f:" DATA (NHEL(I, 17),I=1,5) / 1, 1, 0, 1,-1/"
analyzeline: implied-DO list "(nhel(i, 17),i=1,5)" is not supported. Skipping.
Line #102 in matrix.f:" DATA (NHEL(I, 18),I=1,5) / 1, 1, 0, 1, 1/"
analyzeline: implied-DO list "(nhel(i, 18),i=1,5)" is not supported. Skipping.
Line #103 in matrix.f:" DATA (NHEL(I, 19),I=1,5) / 1, 1, 0,-1,-1/"
analyzeline: implied-DO list "(nhel(i, 19),i=1,5)" is not supported. Skipping.
Line #104 in matrix.f:" DATA (NHEL(I, 20),I=1,5) / 1, 1, 0,-1, 1/"
analyzeline: implied-DO list "(nhel(i, 20),i=1,5)" is not supported. Skipping.
Line #105 in matrix.f:" DATA (NHEL(I, 21),I=1,5) / 1, 1, 1, 1,-1/"
analyzeline: implied-DO list "(nhel(i, 21),i=1,5)" is not supported. Skipping.
Line #106 in matrix.f:" DATA (NHEL(I, 22),I=1,5) / 1, 1, 1, 1, 1/"
analyzeline: implied-DO list "(nhel(i, 22),i=1,5)" is not supported. Skipping.
Line #107 in matrix.f:" DATA (NHEL(I, 23),I=1,5) / 1, 1, 1,-1,-1/"
analyzeline: implied-DO list "(nhel(i, 23),i=1,5)" is not supported. Skipping.
Line #108 in matrix.f:" DATA (NHEL(I, 24),I=1,5) / 1, 1, 1,-1, 1/"
analyzeline: implied-DO list "(nhel(i, 24),i=1,5)" is not supported. Skipping.
Line #109 in matrix.f:" DATA (NHEL(I, 25),I=1,5) /-1,-1,-1, 1,-1/"
analyzeline: implied-DO list "(nhel(i, 25),i=1,5)" is not supported. Skipping.
Line #110 in matrix.f:" DATA (NHEL(I, 26),I=1,5) /-1,-1,-1, 1, 1/"
analyzeline: implied-DO list "(nhel(i, 26),i=1,5)" is not supported. Skipping.
Line #111 in matrix.f:" DATA (NHEL(I, 27),I=1,5) /-1,-1,-1,-1,-1/"
analyzeline: implied-DO list "(nhel(i, 27),i=1,5)" is not supported. Skipping.
Line #112 in matrix.f:" DATA (NHEL(I, 28),I=1,5) /-1,-1,-1,-1, 1/"
analyzeline: implied-DO list "(nhel(i, 28),i=1,5)" is not supported. Skipping.
Line #113 in matrix.f:" DATA (NHEL(I, 29),I=1,5) /-1,-1, 0, 1,-1/"
analyzeline: implied-DO list "(nhel(i, 29),i=1,5)" is not supported. Skipping.
Line #114 in matrix.f:" DATA (NHEL(I, 30),I=1,5) /-1,-1, 0, 1, 1/"
analyzeline: implied-DO list "(nhel(i, 30),i=1,5)" is not supported. Skipping.
Line #115 in matrix.f:" DATA (NHEL(I, 31),I=1,5) /-1,-1, 0,-1,-1/"
analyzeline: implied-DO list "(nhel(i, 31),i=1,5)" is not supported. Skipping.
Line #116 in matrix.f:" DATA (NHEL(I, 32),I=1,5) /-1,-1, 0,-1, 1/"
analyzeline: implied-DO list "(nhel(i, 32),i=1,5)" is not supported. Skipping.
Line #117 in matrix.f:" DATA (NHEL(I, 33),I=1,5) /-1,-1, 1, 1,-1/"
analyzeline: implied-DO list "(nhel(i, 33),i=1,5)" is not supported. Skipping.
Line #118 in matrix.f:" DATA (NHEL(I, 34),I=1,5) /-1,-1, 1, 1, 1/"
analyzeline: implied-DO list "(nhel(i, 34),i=1,5)" is not supported. Skipping.
Line #119 in matrix.f:" DATA (NHEL(I, 35),I=1,5) /-1,-1, 1,-1,-1/"
analyzeline: implied-DO list "(nhel(i, 35),i=1,5)" is not supported. Skipping.
Line #120 in matrix.f:" DATA (NHEL(I, 36),I=1,5) /-1,-1, 1,-1, 1/"
analyzeline: implied-DO list "(nhel(i, 36),i=1,5)" is not supported. Skipping.
Line #121 in matrix.f:" DATA (NHEL(I, 37),I=1,5) /-1, 1,-1, 1,-1/"
analyzeline: implied-DO list "(nhel(i, 37),i=1,5)" is not supported. Skipping.
Line #122 in matrix.f:" DATA (NHEL(I, 38),I=1,5) /-1, 1,-1, 1, 1/"
analyzeline: implied-DO list "(nhel(i, 38),i=1,5)" is not supported. Skipping.
Line #123 in matrix.f:" DATA (NHEL(I, 39),I=1,5) /-1, 1,-1,-1,-1/"
analyzeline: implied-DO list "(nhel(i, 39),i=1,5)" is not supported. Skipping.
Line #124 in matrix.f:" DATA (NHEL(I, 40),I=1,5) /-1, 1,-1,-1, 1/"
analyzeline: implied-DO list "(nhel(i, 40),i=1,5)" is not supported. Skipping.
Line #125 in matrix.f:" DATA (NHEL(I, 41),I=1,5) /-1, 1, 0, 1,-1/"
analyzeline: implied-DO list "(nhel(i, 41),i=1,5)" is not supported. Skipping.
Line #126 in matrix.f:" DATA (NHEL(I, 42),I=1,5) /-1, 1, 0, 1, 1/"
analyzeline: implied-DO list "(nhel(i, 42),i=1,5)" is not supported. Skipping.
Line #127 in matrix.f:" DATA (NHEL(I, 43),I=1,5) /-1, 1, 0,-1,-1/"
analyzeline: implied-DO list "(nhel(i, 43),i=1,5)" is not supported. Skipping.
Line #128 in matrix.f:" DATA (NHEL(I, 44),I=1,5) /-1, 1, 0,-1, 1/"
analyzeline: implied-DO list "(nhel(i, 44),i=1,5)" is not supported. Skipping.
Line #129 in matrix.f:" DATA (NHEL(I, 45),I=1,5) /-1, 1, 1, 1,-1/"
analyzeline: implied-DO list "(nhel(i, 45),i=1,5)" is not supported. Skipping.
Line #130 in matrix.f:" DATA (NHEL(I, 46),I=1,5) /-1, 1, 1, 1, 1/"
analyzeline: implied-DO list "(nhel(i, 46),i=1,5)" is not supported. Skipping.
Line #131 in matrix.f:" DATA (NHEL(I, 47),I=1,5) /-1, 1, 1,-1,-1/"
analyzeline: implied-DO list "(nhel(i, 47),i=1,5)" is not supported. Skipping.
Line #132 in matrix.f:" DATA (NHEL(I, 48),I=1,5) /-1, 1, 1,-1, 1/"
analyzeline: implied-DO list "(nhel(i, 48),i=1,5)" is not supported. Skipping.
Reading file 'coupl.inc' (format:fix)
Line #216 in matrix.f:" DATA (CF(I, 1),I= 1, 1) / 3/"
analyzeline: implied-DO list "(cf(i, 1),i= 1, 1)" is not supported. Skipping.
Reading file 'coupl.inc' (format:fix)
Line #299 in matrix.f:" CALL SETPARA(PATH) !first call to setup the paramaters "
analyzeline: No name/args pattern found for line.
Post-
Block: matrix2py
Block: smatrixhel
Block: smatrix
Block: matrix
Block: get_me
Block: initialise
Post-processing (stage 2)...
Building modules...
Building module "matrix2py"...
Constructing wrapper function "smatrixhel"...
ans = smatrixhel(p,hel)
Constructing wrapper function "smatrix"...
ans = smatrix(p)
Creating wrapper for Fortran function "matrix"
Constructing wrapper function "matrix"...
matrix = matrix(p,nhel,ic)
Constructing wrapper function "get_me"...
ans = get_me(
Constructing wrapper function "initialise"...
Constructing COMMON block support for "heluserchoice"...
userhel
Constructing COMMON block support for "widths"...
Constructing COMMON block support for "rscale"...
mu_r
Constructing COMMON block support for "weak"...
gal
Constructing COMMON block support for "couplings"...
Constructing COMMON block support for "strong"...
g
Constructing COMMON block support for "masses"...
Wrote C/API module "matrix2py" to file "/tmp/tmpWGXxcc
Fortran 77 wrappers are saved to "/tmp/tmpWGXxcc
adding '/tmp/tmpWGXxcc
adding '/tmp/tmpWGXxcc
copying /cvmfs/
copying /cvmfs/
adding '/tmp/tmpWGXxcc
build_src: building npy-pkg config files
running build_ext
customize UnixCCompiler
customize UnixCCompiler using build_ext
customize Gnu95FCompiler
Found executable /cvmfs/
customize Gnu95FCompiler using build_ext
building 'matrix2py' extension
compiling C sources
C compiler: gcc -pthread -fno-strict-
creating /tmp/tmpWGXxcc/tmp
creating /tmp/tmpWGXxcc/
creating /tmp/tmpWGXxcc/
compile options: '-I/tmp/
gcc: /tmp/tmpWGXxcc/
gcc: /tmp/tmpWGXxcc/
compiling Fortran sources
Fortran f77 compiler: /cvmfs/
Fortran f90 compiler: /cvmfs/
Fortran fix compiler: /cvmfs/
compile options: '-I/tmp/
gfortran:f77: matrix.f
matrix.
MATRIX = MATRIX+
Warning: Possible change of value in conversion from COMPLEX(8) to REAL(8) at (1)
gfortran:f77: /tmp/tmpWGXxcc/
/cvmfs/
/cvmfs/
collect2: error: ld returned 1 exit status
/cvmfs/
collect2: error: ld returned 1 exit status
error: Command "/cvmfs/
make: *** [matrix2py.so] Error 1
<======
What is the solution to this compilation error in reweighting?
Also, could you please confirm that I am trying to run Kamikaze reweighting in a correct way?
I do not see the "keep_rwgt_info" parameter in my cards... Can you please point me to where this is set?
Thank you in advance!
Greetings,
Senka
Question information
- Language:
- English Edit question
- Status:
- Answered
- Assignee:
- No assignee Edit question
- Last query:
- Last reply:
Can you help with this problem?
Provide an answer of your own, or ask senka duric for more information if necessary.