#799 closed task (fixed)
Revisit share/examples/DrellYanMatching
| Reported by: | Juergen Reuter | Owned by: | Bijan Chokoufe Nejad |
|---|---|---|---|
| Priority: | P0 | Milestone: | v2.5.0 |
| Component: | shower | Version: | 2.3.0 |
| Severity: | major | Keywords: | |
| Cc: |
Description
It looks like as they are broken and no longer working, or at least they are not working for LHAPDF6?
Change History (8)
comment:1 Changed 8 years ago by
comment:2 Changed 8 years ago by
This is the outcome, so even with one event it is not working, it's just stuck.
n_events = 2
?rebuild_events = true
$sample = "DrellYanMatchingP"
checkpoint = 2500
| Starting simulation for processes 'qqee' etc.
| Simulate: using integration grids from file 'qqee_m1.vg'
| Simulate: activating parton shower
LHAPDF 6.1.6 loading /usr/local/share/LHAPDF/CT10/CT10_0000.dat
CT10 PDF set, member #0, version 4; LHAPDF ID = 10800
| Shower: interfacing LHAPDF set #1
MSTU(12) changed from 0 to 12345
| Shower: Using PYTHIA6 shower
| Simulate: applying MLM matching
| Simulate: activating hadronization
| Hadronization: Using PYTHIA6 interface for hadronization and decays
| Simulate: using integration grids from file 'qqeeg_m1.vg'
| Simulate: activating parton shower
LHAPDF 6.1.6 loading /usr/local/share/LHAPDF/CT10/CT10_0000.dat
CT10 PDF set, member #0, version 4; LHAPDF ID = 10800
| Shower: interfacing LHAPDF set #1
| Shower: Using PYTHIA6 shower
| Simulate: applying MLM matching
| Simulate: activating hadronization
| Hadronization: Using PYTHIA6 interface for hadronization and decays
| Simulate: using integration grids from file 'qgeeq_m1.vg'
| Simulate: activating parton shower
LHAPDF 6.1.6 loading /usr/local/share/LHAPDF/CT10/CT10_0000.dat
CT10 PDF set, member #0, version 4; LHAPDF ID = 10800
| Shower: interfacing LHAPDF set #1
| Shower: Using PYTHIA6 shower
| Simulate: applying MLM matching
| Simulate: activating hadronization
| Hadronization: Using PYTHIA6 interface for hadronization and decays
| Simulate: using integration grids from file 'gqeeq_m1.vg'
| Simulate: activating parton shower
LHAPDF 6.1.6 loading /usr/local/share/LHAPDF/CT10/CT10_0000.dat
CT10 PDF set, member #0, version 4; LHAPDF ID = 10800
| Shower: interfacing LHAPDF set #1
| Shower: Using PYTHIA6 shower
| Simulate: applying MLM matching
| Simulate: activating hadronization
| Hadronization: Using PYTHIA6 interface for hadronization and decays
| RNG: Initializing TAO random-number generator
| RNG: Setting seed for random-number generator to 14712
| Simulation: requested number of events = 2
| corr. to luminosity [fb-1] = 8.4663E-06
| Events: writing to LHEF file 'DrellYanMatchingP.lhe'
| Events: writing to raw file 'DrellYanMatchingP.evx'
| Events: generating 2 unweighted, unpolarized events ...
| Events: event normalization mode '1'
1****************** PYINIT: initialization of PYTHIA routines *****************
==== LHAPDF6 USING PYTHIA-TYPE LHAGLUE INTERFACE ====
LHAPDF 6.1.6 loading /usr/local/share/LHAPDF/cteq6l1/cteq6l1_0000.dat
cteq6l1 PDF set, member #0, version 4; LHAPDF ID = 10042
==============================================================================
I I
I PYTHIA will be initialized for p+ on pbar- user configuration I
I with 980.000 GeV on 980.000 GeV beam energies I
I I
I corresponding to 1960.000 GeV center-of-mass energy I
I I
==============================================================================
Warning: requested subprocess 96 has vanishing cross-section.
Process switched off!
******** PYMAXI: summary of differential cross-section maximum search ********
==========================================================
I I I
I ISUB Subprocess name I Maximum value I
I I I
==========================================================
I I I
I 4 User process 1 I 1.9064D-07 I
I 6 User process 2 I 3.1516D-08 I
I 7 User process 3 I 7.3852D-09 I
I 9 User process 4 I 6.6894D-09 I
I 96 Semihard QCD 2 -> 2 I 0.0000D+00 I
I I I
==========================================================
****** PYMULT: initialization of multiple interactions for MSTP(82) = 4 ******
pT0 = 1.84 GeV gives sigma(parton-parton) = NaN mb: accepted
****** PYMIGN: initialization of multiple interactions for MSTP(82) = 4 ******
pT0 = 1.84 GeV gives sigma(parton-parton) = NaN mb: accepted
********************** PYINIT: initialization completed **********************
comment:3 Changed 8 years ago by
| | % complete | events generated | events remaining | time remaining
| |===================================================================|
| 0.0 0 3 ???
1****************** PYINIT: initialization of PYTHIA routines *****************
==== LHAPDF6 USING PYTHIA-TYPE LHAGLUE INTERFACE ====
LHAPDF 6.1.6 loading /nfs/theoc/data2/bcho/sl6/install/share/LHAPDF/cteq6l1/cteq6l1_0000.dat
cteq6l1 PDF set, member #0, version 4; LHAPDF ID = 10042
==============================================================================
I I
I PYTHIA will be initialized for p+ on pbar- user configuration I
I with 980.000 GeV on 980.000 GeV beam energies I
I I
I corresponding to 1960.000 GeV center-of-mass energy I
I I
==============================================================================
******** PYMAXI: summary of differential cross-section maximum search ********
==========================================================
I I I
I ISUB Subprocess name I Maximum value I
I I I
==========================================================
I I I
I 4 User process 1 I 1.9181D-07 I
I 6 User process 2 I 3.2429D-08 I
I 7 User process 3 I 6.7862D-09 I
I 9 User process 4 I 7.3435D-09 I
I 96 Semihard QCD 2 -> 2 I 5.2875D+03 I
I I I
==========================================================
****** PYMULT: initialization of multiple interactions for MSTP(82) = 4 ******
pT0 = 1.84 GeV gives sigma(parton-parton) = 3.18D+02 mb: accepted
****** PYMIGN: initialization of multiple interactions for MSTP(82) = 4 ******
pT0 = 1.84 GeV gives sigma(parton-parton) = 1.56D+02 mb: accepted
********************** PYINIT: initialization completed **********************
| 100.0 3 0 0m:00s
| |===================================================================|
|
| ... event sample complete.
| Events: actual unweighting efficiency = 0.49 %
| Events: closing LHEF file 'DrellYanMatchingP.lhe'
| Events: closing raw file 'DrellYanMatchingP.evx'
gqeeq:
7.3434552E+03 +- 6.04E+02 fb
qgeeq:
6.7861768E+03 +- 5.56E+02 fb
qqeeg:
3.2428999E+04 +- 3.52E+03 fb
qqee:
1.9180816E+05 +- 1.06E+03 fb
| WHIZARD run finished.
|=============================================================================|
Thanks for using LHAPDF 6.1.6. Please make sure to cite the paper:
Eur.Phys.J. C75 (2015) 3, 132 (http://arxiv.org/abs/1412.7420)
comment:5 Changed 7 years ago by
| Milestone: | v2.4.0 → v2.4.1 |
|---|
comment:6 Changed 7 years ago by
| Priority: | P3 → P0 |
|---|---|
| Severity: | normal → major |
This is obviously easy to check. BACN believes that all is well here, JRR has to confirm.
comment:7 Changed 7 years ago by
| Resolution: | → fixed |
|---|---|
| Status: | new → closed |
JRR convinced himself that this still runs through. Good to close. It is awfully slow, tho. Reduced the number of events in the example file.
Note: See
TracTickets for help on using
tickets.
Please show the errors. I was just able to run
DrellYanMatchingP.sinandDrellYanMatchingW.sinout-of-the-box (n_eventsis quite high in the files. I set it to 3 for testing). LHAPDF version 6.1.6