MadGraph5_aMC@NLO

loop MSSM with Herwig++ problem

Asked by Mihoko_Nojiri

I am trying to generate p p > t1 t1~ [QCD ]with loop_MSSM. Herwig6 works fine but I cannot
make it work with Herwig++, though I have set path to the HepMC and so on. It seems
it is HepMC problem. If you can suggest something that would be very very helpful.
The system is Linux. gcc is version 4.9.1 (Red Hat 4.9.1-10)

I think only change I did from default run_card is HERWIG6-> HERWIGPP.

launch auto
Traceback (most recent call last):
  File "/disk1/nojiri/MG5_aMC_v2_3_3/madgraph/interface/extended_cmd.py", line 908, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/disk1/nojiri/MG5_aMC_v2_3_3/madgraph/interface/extended_cmd.py", line 897, in onecmd_orig
    return func(arg, **opt)
  File "/disk1/nojiri/MG5_aMC_v2_3_3/madgraph/interface/amcatnlo_run_interface.py", line 1225, in do_launch
    self.run_mcatnlo(evt_file)
  File "/disk1/nojiri/MG5_aMC_v2_3_3/madgraph/interface/amcatnlo_run_interface.py", line 3011, in run_mcatnlo
    files.mv(os.path.join(rundir, 'events.%s.gz' % ext), hep_file)
  File "/disk1/nojiri/MG5_aMC_v2_3_3/madgraph/iolibs/files.py", line 172, in mv
    shutil.move(path1, path2)
  File "/opt/rh/python27/root/usr/lib64/python2.7/shutil.py", line 301, in move
:
    copy2(src, real_dst)
  File "/opt/rh/python27/root/usr/lib64/python2.7/shutil.py", line 130, in copy2
    copyfile(src, dst)
  File "/opt/rh/python27/root/usr/lib64/python2.7/shutil.py", line 82, in copyfile
    with open(src, 'rb') as fsrc:
IOError: [Errno 2] No such file or directory: '/disk1/nojiri/MG5_aMC_v2_3_3/testnlo/MCatNLO/RUN_HERWIGPP_2/events.hepmc.gz'
Value of current Options:
              text_editor : None
      notification_center : True
                    pjfry : None
       cluster_local_path : /cvmfs/cp3.uclouvain.be/madgraph/
       group_subprocesses : Auto
ignore_six_quark_processes : False
    loop_optimized_output : True
    cluster_status_update : (600, 30)
             pythia8_path : /disk1/nojiri/pythia8210
                hwpp_path : /usr/local
                    golem : None
          pythia-pgs_path : /disk1/nojiri/MG5_aMC_v2_3_3/pythia-pgs
                  td_path : /disk1/nojiri/MG5_aMC_v2_3_3/td
             delphes_path : /disk1/nojiri/MG5_aMC_v2_3_3/Delphes
              thepeg_path : /usr/local/lib
             cluster_type : condor
      exrootanalysis_path : /disk1/nojiri/MG5_aMC_v2_3_3/ExRootAnalysis
                      OLP : MadLoop
                 applgrid : applgrid-config
               eps_viewer : None
                  fastjet : fastjet-config
                 run_mode : 2
              web_browser : None
   automatic_html_opening : False
        cluster_temp_path : None
             cluster_size : 100
            cluster_queue : None
             syscalc_path : /disk1/nojiri/MG5_aMC_v2_3_3/SysCalc
         madanalysis_path : /disk1/nojiri/MG5_aMC_v2_3_3/MadAnalysis
                   lhapdf : lhapdf-config
             stdout_level : 20
            f2py_compiler : None
               hepmc_path : /usr/local/lib
                  amcfast : amcfast-config
       cluster_retry_wait : 300
      output_dependencies : external
         fortran_compiler : None
                  nb_core : 40
         loop_color_flows : False
              auto_update : 7
         cluster_nb_retry : 1
                 mg5_path : /disk1/nojiri/MG5_aMC_v2_3_3
                  timeout : 60
                    gauge : unitary
      complex_mass_scheme : False
             cpp_compiler : None
   max_npoint_for_channel : 0
#************************************************************
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
:
         VERSION 2.3.3 2015-10-25 *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https://server06.fynu.ucl.ac.be/projects/madgraph *
#* *
#************************************************************
#* *
#* Command File for MadGraph5_aMC@NLO *
#* *
#* run as ./bin/mg5_aMC filename *
#* *
#************************************************************
set group_subprocesses Auto
set ignore_six_quark_processes False
set loop_optimized_output True
set loop_color_flows False
set gauge unitary
set complex_mass_scheme False
set max_npoint_for_channel 0
import model sm
define p = g u c d s u~ c~ d~ s~
define j = g u c d s u~ c~ d~ s~
define l+ = e+ mu+
define l- = e- mu-
define vl = ve vm vt
define vl~ = ve~ vm~ vt~
import model loop_MSSM
generate p p > t1~ t1 [QCD]
output testnlo

######################################################################
## PARAM_CARD AUTOMATICALY GENERATED BY MG5 FOLLOWING UFO MODEL ####
######################################################################
## ##
## Width set on Auto will be computed following the information ##
## present in the decay.py files of the model. ##
## See arXiv:1402.1178 for more details. ##
## ##
######################################################################

###################################
## INFORMATION FOR FRALPHA
###################################
Block fralpha
    1 -1.138252e-01 # alp

###################################
## INFORMATION FOR HMIX
###################################
Block hmix
    1 3.576810e+02 # RMUH
    2 9.748624e+00 # tb

###################################
## INFORMATION FOR LOOP
###################################
Block loop
    1 9.118800e+01 # MU_R

###################################
## INFORMATION FOR MASS
###################################
Block mass
    5 4.880000e+00 # MB
    6 1.750000e+02 # MT
   15 1.777000e+00 # Mta
   23 9.118760e+01 # MZ
   24 7.982901e+01 # MW
   25 1.108991e+02 # MH01
   35 3.999601e+02 # MH02
   36 3.995839e+02 # MA0
   37 4.078790e+02 # MH
  1000001 5.684411e+02 # set of param :1*Msd1, 1*Msd2
  1000002 5.611190e+02 # set of param :1*Msu1, 1*Msu2
  1000005 5.130652e+02 # Msd3
  1000006 3.996685e+02 # Msu3
:
  1000011 2.029157e+02 # set of param :1*Msl1, 1*Msl2
  1000012 1.852583e+02 # set of param :1*Msn1, 1*Msn2
  1000015 1.344909e+02 # Msl3
  1000016 1.847085e+02 # Msn3
  1000021 6.077137e+02 # Mgo
  1000022 9.668807e+01 # Mneu1
  1000023 1.810882e+02 # Mneu2
  1000024 1.816965e+02 # Mch1
  1000025 -3.637560e+02 # Mneu3
  1000035 3.817294e+02 # Mneu4
  1000037 3.799393e+02 # Mch2
  2000001 5.452285e+02 # set of param :1*Msd4, 1*Msd5
  2000002 5.492593e+02 # set of param :1*Msu4, 1*Msu5
  2000005 5.437267e+02 # Msd6
  2000006 5.857858e+02 # Msu6
  2000011 1.441028e+02 # set of param :1*Msl4, 1*Msl5
  2000015 2.068678e+02 # Msl6
## Dependent parameters, given by model restrictions.
## Those values should be edited following the
## analytical expression. MG5 ignores those values
## but they are important for interfacing the output of MG5
## to external program such as Pythia.
  1 0.000000 # d : 0.0
  2 0.000000 # u : 0.0
  3 0.000000 # s : 0.0
  4 0.000000 # c : 0.0
  11 0.000000 # e- : 0.0
  12 0.000000 # ve : 0.0
  13 0.000000 # mu- : 0.0
  14 0.000000 # vm : 0.0
  16 0.000000 # vt : 0.0
  21 0.000000 # g : 0.0
  22 0.000000 # a : 0.0
  1000014 185.258300 # svm : Msn1
  1000013 202.915700 # mul- : Msl1
  2000013 144.102800 # mur- : Msl4
  1000004 561.119000 # cl : Msu1
  2000004 549.259300 # cr : Msu4
  1000003 568.441100 # sl : Msd1
  2000003 545.228500 # sr : Msd4
  9000003 91.187600 # ghz : MZ
  9000004 79.829010 # ghwp : MW
  9000005 79.829010 # ghwm : MW
...

Question information

Language:
English Edit question
Status:
Answered
For:
MadGraph5_aMC@NLO Edit question
Assignee:
Paolo Torrielli Edit question
Last query:
Last reply:
Revision history for this message
Paolo Torrielli (paolo-torrielli) said : 		#1

Hi Mihoko,
in RUN_HERWIGPP_2 there should be log files that help understanding
the problem. For example You should find mcatnlo_run.log and/or
Herwigpp.log (or similar names). Could you tell me what these file say?
Thanks.
Cheers.
Paolo

Revision history for this message
Mihoko_Nojiri (nojiri) said : 		#2

Thank you Paolo, the structure of directories are very different
 and it was not easy for me to find log files.
Now I opened mcatnlo_run.log.

It says

Parameters for the mixing of B_sbar0 and B_s0
x = 25.9622 y = 0.0740952
Integrated mixing probability = 0.499263
Parameters for the mixing of Bbar0 and B0
x = 0.777938 y = 0
Integrated mixing probability = 0.18851
Parameters for the mixing of B0 and Bbar0
x = 0.777938 y = 0
Integrated mixing probability = 0.18851
Parameters for the mixing of B_s0 and B_sbar0
x = 25.9622 y = 0.0740952
Integrated mixing probability = 0.499263
Error: LesHouchesReader 'LHEReader' found unknown particle ID 1000006 in Les Houches common block structure.
You need to define the new particle in an input file.
Herwig++: EventGenerator not available.
Check if 'MCATNLO_HERWIGPP.run' is a valid run file.

...wow
upper output lines are standard Herwig++ output and
then it says it does not know about 1000006.
 I need herwig++ setup for MSSM... Maybe I can fix it
if I can change default ***.in file for Herwig++. Do you have suggestion?

Mihoko

Revision history for this message
Paolo Torrielli (paolo-torrielli) said : 		#3

Hi Mihoko,

thanks for the feedback, I’ll ask one of the Herwig++ authors and
get back to you.

Cheers.
Paolo

On 25 Jan 2016, at 13:02, Mihoko_Nojiri <email address hidden> wrote:

> Question #283312 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/283312
>
> Status: Needs information => Open
>
> Mihoko_Nojiri gave more information on the question:
> Thank you Paolo, the structure of directories are very different
> and it was not easy for me to find log files.
> Now I opened mcatnlo_run.log.
>
> It says
>
> Parameters for the mixing of B_sbar0 and B_s0
> x = 25.9622 y = 0.0740952
> Integrated mixing probability = 0.499263
> Parameters for the mixing of Bbar0 and B0
> x = 0.777938 y = 0
> Integrated mixing probability = 0.18851
> Parameters for the mixing of B0 and Bbar0
> x = 0.777938 y = 0
> Integrated mixing probability = 0.18851
> Parameters for the mixing of B_s0 and B_sbar0
> x = 25.9622 y = 0.0740952
> Integrated mixing probability = 0.499263
> Error: LesHouchesReader 'LHEReader' found unknown particle ID 1000006 in Les Houches common block structure.
> You need to define the new particle in an input file.
> Herwig++: EventGenerator not available.
> Check if 'MCATNLO_HERWIGPP.run' is a valid run file.
>
> ...wow
> upper output lines are standard Herwig++ output and
> then it says it does not know about 1000006.
> I need herwig++ setup for MSSM... Maybe I can fix it
> if I can change default ***.in file for Herwig++. Do you have suggestion?
>
>
> Mihoko
>
> --
> You received this question notification because you are subscribed to
> the question.

Revision history for this message
Paolo Torrielli (paolo-torrielli) said : 		#4

Dear Mihoko,

could you try and see if a normal SM process works fine?
Namely generate a trivial p p > W+ [QCD] or similar, and
see if the HW++ reader complains in that case as well.

If it does, the problem may be related to something else.
If instead it goes through, then could you attach the Herwig++
input file and a very small MSSM event sample (like 100 events
or so)?

Thanks.
Cheers.
Paolo

On 25 Jan 2016, at 13:02, Mihoko_Nojiri <email address hidden> wrote:

> Question #283312 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/283312
>
> Status: Needs information => Open
>
> Mihoko_Nojiri gave more information on the question:
> Thank you Paolo, the structure of directories are very different
> and it was not easy for me to find log files.
> Now I opened mcatnlo_run.log.
>
> It says
>
> Parameters for the mixing of B_sbar0 and B_s0
> x = 25.9622 y = 0.0740952
> Integrated mixing probability = 0.499263
> Parameters for the mixing of Bbar0 and B0
> x = 0.777938 y = 0
> Integrated mixing probability = 0.18851
> Parameters for the mixing of B0 and Bbar0
> x = 0.777938 y = 0
> Integrated mixing probability = 0.18851
> Parameters for the mixing of B_s0 and B_sbar0
> x = 25.9622 y = 0.0740952
> Integrated mixing probability = 0.499263
> Error: LesHouchesReader 'LHEReader' found unknown particle ID 1000006 in Les Houches common block structure.
> You need to define the new particle in an input file.
> Herwig++: EventGenerator not available.
> Check if 'MCATNLO_HERWIGPP.run' is a valid run file.
>
> ...wow
> upper output lines are standard Herwig++ output and
> then it says it does not know about 1000006.
> I need herwig++ setup for MSSM... Maybe I can fix it
> if I can change default ***.in file for Herwig++. Do you have suggestion?
>
>
> Mihoko
>
> --
> You received this question notification because you are subscribed to
> the question.

Can you help with this problem?

Provide an answer of your own, or ask Mihoko_Nojiri for more information if necessary.

To post a message you must log in.