s-channel subtraction in single-top: problem with MadSpin

[Alberto Orso Maria Iorio]

Hi all,

I don't know if it's a bug or a feature, but I report it here nevertheless:

Following the example on the aMC@NLO page I am generating events for single-top in the s-channel, t-channel 4FS and t-channel 5FS, respectively:

generate p p > t b~ [QCD]

generate p p > t b~ j $$ w+ w- [QCD]

generate p p > t j $$ w+ w- [QCD] (using loop_sm-no_b_mass and adding b's to the p and j definition)

I manage to decay the s-channel with MadSpin, but when I try to use the $$ w+ w- option in any process I get an error like this (this in particular is for t-channel 4FS):

=== Error on screen:
Command "add process p p > t j $$ w+ w- j " interrupted in sub-command:
"add process p p > t j $$ w+ w- j" with error:
InvalidCmd : Duplicate process Process: u b > t d WEIGHTED=4 $$ w+ w- g u c d s b u~ c~ d~ s~ b~ found. Please check your processes.
INFO: Done 3.939
INFO: generating the full square matrix element (with decay)
INFO: generate p p > t j $$ w+ w- , (t > b w+ , w+ > all all );add process p p > t j $$ w+ w- j, (t > b w+ , w+ > all all ) ;
Command "add process p p > t j $$ w+ w- j, (t > b w+ , w+ > all all ) " interrupted in sub-command:
"add process p p > t j $$ w+ w- j, (t > b w+ , w+ > all all )" with error:
InvalidCmd : Duplicate process Process: u b > t d WEIGHTED=4 $$ w+ w- g u c d s b u~ c~ d~ s~ b~
 Process: c b > t s WEIGHTED=4 $$ w+ w- g u c d s b u~ c~ d~ s~ b~
 Process: b u > t d WEIGHTED=4 $$ w+ w- g u c d s b u~ c~ d~ s~ b~
 Process: b c > t s WEIGHTED=4 $$ w+ w- g u c d s b u~ c~ d~ s~ b~
 Process: b d~ > t u~ WEIGHTED=4 $$ w+ w- g u c d s b u~ c~ d~ s~ b~
 Process: b s~ > t c~ WEIGHTED=4 $$ w+ w- g u c d s b u~ c~ d~ s~ b~
 Process: d~ b > t u~ WEIGHTED=4 $$ w+ w- g u c d s b u~ c~ d~ s~ b~
 Process: s~ b > t c~ WEIGHTED=4 $$ w+ w- g u c d s b u~ c~ d~ s~ b~
 Decays:
   Process: t > b w+ WEIGHTED=2
   Decays:
     Process: w+ > u d~ WEIGHTED=2
     Process: w+ > c s~ WEIGHTED=2
     Process: w+ > b~ t WEIGHTED=2
     Process: w+ > a w+ WEIGHTED=2
            Process: w+ > ve e+ WEIGHTED=2
     Process: w+ > vm mu+ WEIGHTED=2
     Process: w+ > vt ta+ WEIGHTED=2
     Process: w+ > z w+ WEIGHTED=2
     Process: w+ > w+ h WEIGHTED=2 found. Please check your processes.
INFO: Done 4.897
INFO: generate matrix element for decay only (1 - > N).
INFO: output standalone_ms /afs/cern.ch/work/o/oiorio/CMSSW_5_3_8_patch1/src/Madgraph/MadGraph5_v2_0_0_beta3/t_5FS/decay_me
INFO: Done 3.129
INFO: Compiling code
INFO: detect independant decays
INFO: Done in 0.108659982681s
INFO:
INFO: Estimating the maximum weight
INFO: *****************************
INFO: Probing the first 75 events
INFO: with 400 phase space points
INFO:
INFO: Event 1/75 : 0.57s
            decay_config 1 [t_wpb_wp_udx] : 8.47933449263e-06
           decay_config 2 [t_wpb_wp_vttap] : 2.88049828644e-06
Command "launch auto " interrupted with error:
KeyError : ((4, 21), (-5, 3, 6))
Please report this bug on https://bugs.launchpad.net/madgraph5
More information is found in '/afs/cern.ch/work/o/oiorio/CMSSW_5_3_8_patch1/src/Madgraph/MadGraph5_v2_0_0_beta3/t_5FS/run_01_tag_1_debug.log'.
Please attach this file to your report.
===

At the end of the mail I attach the content of the debug.log

I also tried to generate t+s channel, i.e.

generate p p > t j [QCD]

With those MadSpin works fine.

 Cheers and many thanks in advance,
  Orso

content of t_5FS/run_01_tag_1_debug.log

#************************************************************
#* MadGraph/aMC@NLO 5 *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.0.0.beta3 2013-02-14 *
#* *
#* The MadGraph Development Team - Please visit us at *
#* https://server06.fynu.ucl.ac.be/projects/madgraph *
#* *
#************************************************************
#* *
#* Command File for aMCatNLO *
#* *
#* run as ./bin/madevent.py filename *
#* *
#************************************************************
launch auto
Traceback (most recent call last):
  File "/afs/cern.ch/work/o/oiorio/CMSSW_5_3_8_patch1/src/Madgraph/MadGraph5_v2_0_0_beta3/madgraph/interface/extended_cmd.py", line 825, in one
cmd
    return self.onecmd_orig(line, **opt)
  File "/afs/cern.ch/work/o/oiorio/CMSSW_5_3_8_patch1/src/Madgraph/MadGraph5_v2_0_0_beta3/madgraph/interface/extended_cmd.py", line 818, in one
cmd_orig
    return func(arg, **opt)
  File "/afs/cern.ch/work/o/oiorio/CMSSW_5_3_8_patch1/src/Madgraph/MadGraph5_v2_0_0_beta3/madgraph/interface/amcatnlo_run_interface.py", line 1
062, in do_launch
    self.exec_cmd('decay_events -from_cards', postcmd=False)
  File "/afs/cern.ch/work/o/oiorio/CMSSW_5_3_8_patch1/src/Madgraph/MadGraph5_v2_0_0_beta3/madgraph/interface/extended_cmd.py", line 865, in exe
c_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/afs/cern.ch/work/o/oiorio/CMSSW_5_3_8_patch1/src/Madgraph/MadGraph5_v2_0_0_beta3/madgraph/interface/extended_cmd.py", line 818, in one
cmd_orig
    return func(arg, **opt)
  File "/afs/cern.ch/work/o/oiorio/CMSSW_5_3_8_patch1/src/Madgraph/MadGraph5_v2_0_0_beta3/madgraph/interface/common_run_interface.py", line 132
1, in do_decay_events
    madspin_cmd.import_command_file(path)
  File "/afs/cern.ch/work/o/oiorio/CMSSW_5_3_8_patch1/src/Madgraph/MadGraph5_v2_0_0_beta3/madgraph/interface/extended_cmd.py", line 975, in imp
ort_command_file
    self.exec_cmd(line, precmd=True)
  File "/afs/cern.ch/work/o/oiorio/CMSSW_5_3_8_patch1/src/Madgraph/MadGraph5_v2_0_0_beta3/madgraph/interface/extended_cmd.py", line 865, in exe
c_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/afs/cern.ch/work/o/oiorio/CMSSW_5_3_8_patch1/src/Madgraph/MadGraph5_v2_0_0_beta3/madgraph/interface/extended_cmd.py", line 818, in one
cmd_orig
    return func(arg, **opt)
  File "/afs/cern.ch/work/o/oiorio/CMSSW_5_3_8_patch1/src/Madgraph/MadGraph5_v2_0_0_beta3/MadSpin/interface_madspin.py", line 378, in do_launch
    generate_all.run()
  File "/afs/cern.ch/work/o/oiorio/CMSSW_5_3_8_patch1/src/Madgraph/MadGraph5_v2_0_0_beta3/MadSpin/decay.py", line 2478, in run
    self.get_max_weight_from_event(decay_mapping)
  File "/afs/cern.ch/work/o/oiorio/CMSSW_5_3_8_patch1/src/Madgraph/MadGraph5_v2_0_0_beta3/MadSpin/decay.py", line 3234, in get_max_weight_from_
event
    production_tag, event_map = self.load_event()
  File "/afs/cern.ch/work/o/oiorio/CMSSW_5_3_8_patch1/src/Madgraph/MadGraph5_v2_0_0_beta3/MadSpin/decay.py", line 3391, in load_event
    P_order = self.all_ME[production_tag]['tag2order'][production_tag]
KeyError: ((4, 21), (-5, 3, 6))
Value of current Options:
              web_browser : None
              text_editor : None
        cluster_temp_path : None
       group_subprocesses : Auto
            cluster_queue : None
         madanalysis_path : None
                   lhapdf : lhapdf-config
ignore_six_quark_processes : False
                 mg5_path : /afs/cern.ch/work/o/oiorio/CMSSW_5_3_8_patch1/src/Madgraph/MadGraph5_v2_0_0_beta3
    loop_optimized_output : True
    cluster_status_update : (600, 30)
             pythia8_path : None
   automatic_html_opening : True
             stdout_level : 20
          pythia-pgs_path : None
                  td_path : None
             delphes_path : None
              auto_update : 7
             cluster_type : condor
         fortran_compiler : None
   MCatNLO-utilities_path : None
      exrootanalysis_path : None
               eps_viewer : None
                  timeout : 60
                  fastjet : /afs/cern.ch/work/o/oiorio/CMSSW_5_3_8_patch1/src/Madgraph/fastjet-3.0.3/fastjet-config
                    gauge : unitary
      complex_mass_scheme : False
                  nb_core : 16
                 run_mode : 2

[Olivier Mattelaer (olivier-mattelaer)]

Thanks for reporting,
I will look at it, and should have a fix very quickly.

Thanks,

Olivier

[Olivier Mattelaer (olivier-mattelaer)]

Hi Alberto,

I've attach a patch for this bug. As you pointed correctly MS was not supporting correctly the presence of $$ (or /) options.
This patch fixes the problem.

In order to apply the patch you need to do:
patch =p1 < bug1176874.patch

Cheers,

Olivier