aMCatNLOError : An error occurred during reweight.
Dear experts,
I'm generating
import model loop_sm-ckm
generate p p > t b~ j $$ w+ w- [QCD]
with a modified setscales.f to account for the recommended scale in the 4FS tchannel.
###
###
The crucial part to pick the kinematics of the b-quark in the setcales.f looks like:
do i=nincoming+
enddo
if (nexternal.eq.5) then
tmp = dot(pp(
ptb = (pt(pp(
tmpmass = sqrt(tmp)
elseif (nexternal.eq.6) then
tmp = dot(pp(
ptb = (pt(pp(
tmpmass = sqrt(tmp)
###
###
Now, I get the following error after trying to generate 250000 events:
INFO: Idle: 0, Running: 4, Completed: 64 [ 4h 19m ]
INFO: Idle: 0, Running: 3, Completed: 65 [ 4h 19m ]
INFO: Idle: 0, Running: 2, Completed: 66 [ 4h 19m ]
INFO: Idle: 0, Running: 1, Completed: 67 [ 4h 19m ]
INFO: Idle: 0, Running: 0, Completed: 68 [ 4h 21m ]
INFO: Doing reweight
INFO: Idle: 48, Running: 20, Completed: 0 [ current time: 17h26 ]
INFO: Idle: 47, Running: 20, Completed: 1 [ 8.3s ]
INFO: Idle: 46, Running: 20, Completed: 2 [ 8.4s ]
...
INFO: Idle: 0, Running: 10, Completed: 58 [ 12m 19s ]
INFO: Idle: 0, Running: 8, Completed: 60 [ 12m 49s ]
INFO: Idle: 0, Running: 4, Completed: 64 [ 13m 24s ]
INFO: Idle: 0, Running: 0, Completed: 68 [ 14m 4s ]
Error detected in "launch "
write debug file /afs/cern.
If you need help with this issue please contact us on https:/
aMCatNLOError : An error occurred during reweight. Check the'reweight_
quit
####
#### The reweight_
####
==== POWHEG (OR POSSIBLY HERWIG) WILL USE LHAPDFv6 ====
LHAPDF 6.1.4 loading /cvmfs/
NNPDF30_
==== POWHEG (OR POSSIBLY HERWIG) WILL USE LHAPDFv6 ====
LHAPDF 6.1.4 loading /cvmfs/
NNPDF30_
LHAPDF 6.1.4 loading /cvmfs/
NNPDF30_
LHAPDF 6.1.4 loading /cvmfs/
NNPDF30_
LHAPDF 6.1.4 loading /cvmfs/
...
NNPDF30_
LHAPDF 6.1.4 loading /cvmfs/
NNPDF30_
LHAPDF 6.1.4 loading /cvmfs/
NNPDF30_
A PDF is used, so alpha_s(MZ) is going to be modified
Old value of alpha_s from param_card: 0.11799999999999999
New value of alpha_s from PDF lhapdf : 0.11227348553523188
using LHAPDF
Enter event file name
Enter 1 to save all cross sections on tape
0 otherwise
Using:
QES_over_ref: 1.0000000000000000
muR_over_ref: 1.0000000000000000
muF1_over_ref: 1.0000000000000000
muF2_over_ref: 1.0000000000000000
Doing scale reweight:
0.50000000000
0.50000000000
Doing PDF reweight:
Central set id: 260400
Min error set id: 260401
Max error set id: 260500
The events appear to be unweighted
Will store the ratios of recomputed weights
over reference weights
======
process combination map (specified per FKS dir):
1 map 1 2
1 inv. map 1 2
2 map 1 2
2 inv. map 1 2
3 map 1 2
3 inv. map 1 2
4 map 1 2
4 inv. map 1 2
5 map 1 2
5 inv. map 1 2
6 map 1 2
6 inv. map 1 2
7 map 1 1 1 1 2 2 2 2
7 inv. map 4 8
8 map 1 2
8 inv. map 1 2
9 map 1 1 2 2
9 inv. map 2 4
10 map 1 2
10 inv. map 1 2
======
Scale values (may change event by event):
muR, muR_reference: 0.957774D+02 0.957774D+02 1.00
muF1, muF1_reference: 0.957774D+02 0.957774D+02 1.00
muF2, muF2_reference: 0.957774D+02 0.957774D+02 1.00
QES, QES_reference: 0.957774D+02 0.957774D+02 1.00
muR_reference [functional form]:
4*Sqrt[m(b)**2 + pT(b)**2], b=b quark mod 3
muF1_reference [functional form]:
4*Sqrt[m(b)**2 + pT(b)**2], b=b quark mod 3
muF2_reference [functional form]:
4*Sqrt[m(b)**2 + pT(b)**2], b=b quark mod 3
QES_reference [functional form]:
4*Sqrt[m(b)**2 + pT(b)**2], b=b quark mod 3
alpha_s= 0.11140954936075177
#######
#######
So I supsect something is going wrong with the scale reweighting.
However, a run with 10000 events worked fine, and the corresponding antitop generation also runs smoothly (with another setscales.f to pick the b quark there, too, which is just the nexternal-2 and nexternal-3 respectively to account for the different order of particles in the folder_structure).
I checked that I picked the b quark every time, and previous versions of madgraph like 2.1.2 etc worked perfectly fine with scale reweighting, but now I get this error, and I don't know what is the problem :)
Thanks,
Benedikt
Question information
- Language:
- English Edit question
- Status:
- Open
- Assignee:
- Rikkert Frederix Edit question
- Last query:
- Last reply:
Can you help with this problem?
Provide an answer of your own, or ask Benedikt Maier for more information if necessary.