Enable durham kt cut

Asked by Alec Drobac on 2020-07-06

I am receiving errors when I try to enable a durham kt cut. My setup is the following:
 - Version 2.7.3. (I have also tried with 2.7.2 and the v3 beta)
 - install boost, lahpdf6, pythia8 (in that order)

Then I launch MG with ./bin/mg5_aMC and do the following:

generate p p > z, z > e+ e-
output zeeDurham
launch zeeDurham
2 (to open run_card.dat)
 - Here I add the following lines:
#*********************************************************************
# Apply pt/E/eta/dr/mij/kt_durham cuts on decay products or not
# (note that etmiss/ptll/ptheavy/ht/sorted cuts always apply)
#*********************************************************************
   False = cut_decays ! Cut decay products
   ## WHAT I ADDED##
   0 = xqcut
   30d0 = ktdurham ! minimum durham kt jet measure between partons
   0.4 = dparameter ! D parameter for hadronic kT
   T = dokt ! T: enable durham kt cut
   ##############
3 (to open pythia8_card.dat)
  - Set Merging:Process = pp>(z>e+e-)

Doing this gives me the following error:
-------------------------------------------
Error detected in "generate_events run_01"
write debug file /afs/cern.ch/user/a/adrobac/MadGraph/MG5_aMC_v3_0_2/zeeDurham/run_01_tag_1_debug.log
If you need help with this issue please contact us on https://answers.launchpad.net/mg5amcnlo
MadGraph5Error : A compilation Error occurs when trying to compile /afs/cern.ch/user/a/adrobac/MadGraph/MG5_aMC_v3_0_2/zeeDurham/Source.
        The compilation fails with the following output message:
            gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o StringCast.o StringCast.f
            gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o ranmar.o ranmar.f
            gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o alfas_functions.o alfas_functions.f
            gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o transpole.o transpole.f
            gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o invarients.o invarients.f
            gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o hfill.o hfill.f
            gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o pawgraphs.o pawgraphs.f
            gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o ran1.o ran1.f
            gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o rw_events.o rw_events.f
            gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o rw_routines.o rw_routines.f
            gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o kin_functions.o kin_functions.f
            gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o open_file.o open_file.f
            gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o basecode.o basecode.f
            gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o setrun.o setrun.f
            gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o run_printout.o run_printout.f
            gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o dgauss.o dgauss.f
            gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o readgrid.o readgrid.f
            run_card.inc:345.10:
                Included at setrun.f:78:

                  DOKT = ' T '
                      1
            Error: Symbol 'dokt' at (1) has no IMPLICIT type
            gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o getissud.o getissud.f
            make: *** [setrun.o] Error 1
            make: *** Waiting for unfinished jobs....

        Please try to fix this compilations issue and retry.
        Help might be found at https://answers.launchpad.net/mg5amcnlo.
        If you think that this is a bug, you can report this at https://bugs.launchpad.net/mg5amcnlo
--------------------------------------------------------
When I comment out the T = dokt line, the generation runs fine, so it does seem to be a problem with just that line. I also tried replacing
T = dokt ! T: enable durham kt cut
with
T = do_kt_durham ! T: enable durham kt cut
and doing so also causes the generation to fail, but with a different error:
---------------------------------------------------------
Error detected in "generate_events run_01"
write debug file /afs/cern.ch/user/a/adrobac/MadGraph/MG5_aMC_v3_0_2/zeeDurham/run_01_tag_1_debug.log
If you need help with this issue please contact us on https://answers.launchpad.net/mg5amcnlo
MadGraph5Error : A compilation Error occurs when trying to compile /afs/cern.ch/user/a/adrobac/MadGraph/MG5_aMC_v3_0_2/zeeDurham/Source.
        The compilation fails with the following output message:
            gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o DiscreteSampler.o DiscreteSampler.f
            gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o setrun.o setrun.f
            cd DHELAS; make
            cd PDF; make
            cd MODEL; make
            cd CERNLIB; make
            cd BIAS/dummy; make
            make[1]: Entering directory `/afs/cern.ch/user/a/adrobac/MadGraph/MG5_aMC_v3_0_2/zeeDurham/Source/DHELAS'
            make[1]: warning: jobserver unavailable: using -j1. Add `+' to parent make rule.
            make[1]: Entering directory `/afs/cern.ch/user/a/adrobac/MadGraph/MG5_aMC_v3_0_2/zeeDurham/Source/CERNLIB'
            make[1]: warning: jobserver unavailable: using -j1. Add `+' to parent make rule.
            make[1]: Entering directory `/afs/cern.ch/user/a/adrobac/MadGraph/MG5_aMC_v3_0_2/zeeDurham/Source/PDF'
            make[1]: warning: jobserver unavailable: using -j1. Add `+' to parent make rule.
            gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o abend.o abend.f
            gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o aloha_functions.o aloha_functions.f
            gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o Ctq6Pdf.o Ctq6Pdf.f
            make[1]: Entering directory `/afs/cern.ch/user/a/adrobac/MadGraph/MG5_aMC_v3_0_2/zeeDurham/Source/BIAS/dummy'
            make[1]: warning: jobserver unavailable: using -j1. Add `+' to parent make rule.
            gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o dummy.o dummy.f
            make[1]: Entering directory `/afs/cern.ch/user/a/adrobac/MadGraph/MG5_aMC_v3_0_2/zeeDurham/Source/MODEL'
            make[1]: warning: jobserver unavailable: using -j1. Add `+' to parent make rule.
            gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o couplings.o couplings.f
            run_card.inc:345.21:
                Included at setrun.f:78:

                  DO_KT_DURHAM = ' T '
                                 1
            Error: Can't convert CHARACTER(1) to LOGICAL(4) at (1)
            make: *** [setrun.o] Error 1
            make: *** Waiting for unfinished jobs....
            gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o dlsqp2.o dlsqp2.f
            ar cru ../../../lib/libbias.a dummy.o
            gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o lenocc.o lenocc.f
            gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o lha_read.o lha_read.f
            gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o mtlprt.o mtlprt.f
            ranlib ../../../lib/libbias.a
            gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o mtlset.o mtlset.f
            gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o pdfwrap.o pdfwrap.f
            make[1]: Leaving directory `/afs/cern.ch/user/a/adrobac/MadGraph/MG5_aMC_v3_0_2/zeeDurham/Source/BIAS/dummy'
            gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o opendata.o opendata.f
            gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o printout.o printout.f
            gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o radmul.o radmul.f
            gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o pdf.o pdf.f
            gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o rw_para.o rw_para.f
            gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o PhotonFlux.o PhotonFlux.f
            ar cru ../../lib/libcernlib.a abend.o dlsqp2.o lenocc.o mtlprt.o mtlset.o radmul.o
            gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o model_functions.o model_functions.f
            gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o pdg2pdf.o pdg2pdf.f
            ranlib ../../lib/libcernlib.a
            gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o couplings1.o couplings1.f
            gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o couplings2.o couplings2.f
            make[1]: Leaving directory `/afs/cern.ch/user/a/adrobac/MadGraph/MG5_aMC_v3_0_2/zeeDurham/Source/CERNLIB'
            gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o NNPDFDriver.o NNPDFDriver.f
            ar cru ../../lib/libmodel.a couplings.o lha_read.o printout.o rw_para.o model_functions.o couplings1.o couplings2.o
            ranlib ../../lib/libmodel.a
            gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o FFV4_3.o FFV4_3.f
            make[1]: Leaving directory `/afs/cern.ch/user/a/adrobac/MadGraph/MG5_aMC_v3_0_2/zeeDurham/Source/MODEL'
            ar cru ../../lib/libpdf.a Ctq6Pdf.o pdfwrap.o opendata.o pdf.o PhotonFlux.o pdg2pdf.o NNPDFDriver.o
            gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o FFV2_3.o FFV2_3.f
            ranlib ../../lib/libpdf.a
            gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o FFV2_0.o FFV2_0.f
            gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o FFV5_0.o FFV5_0.f
            make[1]: Leaving directory `/afs/cern.ch/user/a/adrobac/MadGraph/MG5_aMC_v3_0_2/zeeDurham/Source/PDF'
            gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o FFV3_0.o FFV3_0.f
            ar cru ../../lib/libdhelas.a aloha_functions.o FFV4_3.o FFV2_3.o FFV2_0.o FFV5_0.o FFV3_0.o
            ranlib ../../lib/libdhelas.a
            make[1]: Leaving directory `/afs/cern.ch/user/a/adrobac/MadGraph/MG5_aMC_v3_0_2/zeeDurham/Source/DHELAS'

        Please try to fix this compilations issue and retry.
        Help might be found at https://answers.launchpad.net/mg5amcnlo.
        If you think that this is a bug, you can report this at https://bugs.launchpad.net/mg5amcnlo
--------------------------------------------------------
Looking at examples of others enabling durham kt cuts, I am not seeing any differences in the syntax I'm using, so I am wondering if there are some other options I need to also be setting. Apologies if I am missing something basic, but any help would be appreciated.

Thanks,
 - Alec

Question information

Language:
English Edit question
Status:
Solved
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Solved by:
Olivier Mattelaer
Solved:
2020-07-14
Last query:
2020-07-14
Last reply:
2020-07-10

Sorry for the late reply,

You just need to set the value of
ktdurham
to a value different of -1

They are no "dokt" parameter define in the run_card.

I would suggest in this case to generate your process like
generate p p > z, z > e+ e-
add process p p > z j , z > e+ e-
This should change the content of the default run_card to include section for MLM (activated by default) and for CKKW.
You can use that run_card for the "generate p p > z, z > e+ e-" case.

Cheers,

Olivier

Alec Drobac (adrobac) said : #2

Thanks Olivier Mattelaer, that solved my question.