#************************************************************ #* MadGraph5_aMC@NLO * #* * #* * * * #* * * * * * #* * * * * 5 * * * * * #* * * * * * #* * * * #* * #* * #* VERSION 2.3.2.2 2015-09-06 * #* * #* The MadGraph5_aMC@NLO Development Team - Find us at * #* https://server06.fynu.ucl.ac.be/projects/madgraph * #* and * #* http://amcatnlo.cern.ch * #* * #************************************************************ #* * #* Command File for aMCatNLO * #* * #* run as ./bin/aMCatNLO.py filename * #* * #************************************************************ launch auto Traceback (most recent call last): File "/afs/cern.ch/user/s/sgreydan/MadGraph5/madgraph/interface/extended_cmd.py", line 889, in onecmd return self.onecmd_orig(line, **opt) File "/afs/cern.ch/user/s/sgreydan/MadGraph5/madgraph/interface/extended_cmd.py", line 882, in onecmd_orig return func(arg, **opt) File "/afs/cern.ch/user/s/sgreydan/MadGraph5/madgraph/interface/amcatnlo_run_interface.py", line 1218, in do_launch self.run_mcatnlo(evt_file) File "/afs/cern.ch/user/s/sgreydan/MadGraph5/madgraph/interface/amcatnlo_run_interface.py", line 2629, in run_mcatnlo files.mv(os.path.join(rundir, 'events.%s.gz' % ext), hep_file) File "/afs/cern.ch/user/s/sgreydan/MadGraph5/madgraph/iolibs/files.py", line 172, in mv shutil.move(path1, path2) File "/usr/lib64/python2.6/shutil.py", line 260, in move copy2(src, real_dst) File "/usr/lib64/python2.6/shutil.py", line 95, in copy2 copyfile(src, dst) File "/usr/lib64/python2.6/shutil.py", line 50, in copyfile with open(src, 'rb') as fsrc: IOError: [Errno 2] No such file or directory: '/afs/cern.ch/user/s/sgreydan/mg5_runs/MG5_cteq6m_NLO/MCatNLO/RUN_HERWIG6_1/events.hep.gz' Value of current Options: text_editor : None 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 : /usr/libexec/CERNLIB/2006a/src/mclibs/pythia hwpp_path : None golem : None pythia-pgs_path : None td_path : None delphes_path : None thepeg_path : None cluster_type : condor exrootanalysis_path : None OLP : MadLoop applgrid : applgrid-config eps_viewer : None fastjet : None run_mode : 2 web_browser : None automatic_html_opening : False cluster_temp_path : None cluster_size : 100 cluster_queue : None syscalc_path : None madanalysis_path : None lhapdf : lhapdf-config stdout_level : 20 f2py_compiler : None hepmc_path : None amcfast : amcfast-config cluster_retry_wait : 300 output_dependencies : external fortran_compiler : None nb_core : 8 loop_color_flows : False auto_update : 7 cluster_nb_retry : 1 mg5_path : /afs/cern.ch/user/s/sgreydan/MadGraph5 timeout : 60 gauge : unitary complex_mass_scheme : False cpp_compiler : None max_npoint_for_channel : 0 #************************************************************ #* MadGraph5_aMC@NLO * #* * #* * * * #* * * * * * #* * * * * 5 * * * * * #* * * * * * #* * * * #* * #* * #* VERSION 2.3.2.2 2015-09-06 * #* * #* 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_color_flows False set gauge unitary set complex_mass_scheme False set max_npoint_for_channel 0 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_sm generate p p > h mu+ mu- [QCD] output MG5_cteq6m_NLO ###################################################################### ## PARAM_CARD AUTOMATICALY GENERATED BY MG5 #### ###################################################################### ################################### ## INFORMATION FOR LOOP ################################### BLOCK LOOP # 1 9.118800e+01 # mu_r ################################### ## INFORMATION FOR MASS ################################### BLOCK MASS # 5 4.700000e+00 # mb 6 1.725000e+02 # mt 15 1.777000e+00 # mta 23 9.118760e+01 # mz 25 1.250900e+02 # mh 1 0.000000e+00 # d : 0.0 2 0.000000e+00 # u : 0.0 3 0.000000e+00 # s : 0.0 4 0.000000e+00 # c : 0.0 11 0.000000e+00 # e- : 0.0 12 0.000000e+00 # ve : 0.0 13 0.000000e+00 # mu- : 0.0 14 0.000000e+00 # vm : 0.0 16 0.000000e+00 # vt : 0.0 21 0.000000e+00 # g : 0.0 22 0.000000e+00 # a : 0.0 24 8.041900e+01 # w+ : cmath.sqrt(mz__exp__2/2. + cmath.sqrt(mz__exp__4/4. - (aew*cmath.pi*mz__exp__2)/(gf*sqrt__2))) ################################### ## INFORMATION FOR SMINPUTS ################################### BLOCK SMINPUTS # 1 1.325070e+02 # aewm1 2 1.166370e-05 # gf 3 1.180000e-01 # as ################################### ## INFORMATION FOR YUKAWA ################################### BLOCK YUKAWA # 5 4.700000e+00 # ymb 6 1.730000e+02 # ymt 15 1.777000e+00 # ymtau ################################### ## INFORMATION FOR DECAY ################################### DECAY 6 1.491500e+00 # wt DECAY 23 2.495200e+00 # wz DECAY 24 2.088720e+00 # ww DECAY 25 4.210000e-03 # wh DECAY 1 0.000000e+00 # d : 0.0 DECAY 2 0.000000e+00 # u : 0.0 DECAY 3 0.000000e+00 # s : 0.0 DECAY 4 0.000000e+00 # c : 0.0 DECAY 5 0.000000e+00 # b : 0.0 DECAY 11 0.000000e+00 # e- : 0.0 DECAY 12 0.000000e+00 # ve : 0.0 DECAY 13 0.000000e+00 # mu- : 0.0 DECAY 14 0.000000e+00 # vm : 0.0 DECAY 15 0.000000e+00 # ta- : 0.0 DECAY 16 0.000000e+00 # vt : 0.0 DECAY 21 0.000000e+00 # g : 0.0 DECAY 22 0.000000e+00 # a : 0.0 #*********************************************************************** # MadGraph5_aMC@NLO * # * # run_card.dat aMC@NLO * # * # This file is used to set the parameters of the run. * # * # Some notation/conventions: * # * # Lines starting with a hash (#) are info or comments * # * # mind the format: value = variable ! comment * #*********************************************************************** # #******************* # Running parameters #******************* # #*********************************************************************** # Tag name for the run (one word) * #*********************************************************************** tag_1 = run_tag ! name of the run #*********************************************************************** # Number of LHE events (and their normalization) and the required * # (relative) accuracy on the Xsec. * # These values are ignored for fixed order runs * #*********************************************************************** 100000 = nevents ! Number of unweighted events requested -1.0 = req_acc ! Required accuracy (-1=auto determined from nevents) -1 = nevt_job ! Max number of events per job in event generation. ! (-1= no split). #*********************************************************************** # Normalize the weights of LHE events such that they sum or average to * # the total cross section * #*********************************************************************** average = event_norm ! average or sum #*********************************************************************** # Number of points per itegration channel (ignored for aMC@NLO runs) * #*********************************************************************** 0.01 = req_acc_FO ! Required accuracy (-1=ignored, and use the ! number of points and iter. below) # These numbers are ignored except if req_acc_FO is equal to -1 5000 = npoints_FO_grid ! number of points to setup grids 4 = niters_FO_grid ! number of iter. to setup grids 10000 = npoints_FO ! number of points to compute Xsec 6 = niters_FO ! number of iter. to compute Xsec #*********************************************************************** # Random number seed * #*********************************************************************** 0 = iseed ! rnd seed (0=assigned automatically=default)) #*********************************************************************** # Collider type and energy * #*********************************************************************** 1 = lpp1 ! beam 1 type (0 = no PDF) 1 = lpp2 ! beam 2 type (0 = no PDF) 7000 = ebeam1 ! beam 1 energy in GeV 7000 = ebeam2 ! beam 2 energy in GeV #*********************************************************************** # PDF choice: this automatically fixes also alpha_s(MZ) and its evol. * #*********************************************************************** cteq6_m = pdlabel ! PDF set 244600 = lhaid ! if pdlabel=lhapdf, this is the lhapdf number #*********************************************************************** # Include the NLO Monte Carlo subtr. terms for the following parton * # shower (HERWIG6 | HERWIGPP | PYTHIA6Q | PYTHIA6PT | PYTHIA8) * # WARNING: PYTHIA6PT works only for processes without FSR!!!! * #*********************************************************************** HERWIG6 = parton_shower #*********************************************************************** # Renormalization and factorization scales * # (Default functional form for the non-fixed scales is the sum of * # the transverse masses of all final state particles and partons. This * # can be changed in SubProcesses/set_scales.f) * #*********************************************************************** False = fixed_ren_scale ! if .true. use fixed ren scale False = fixed_fac_scale ! if .true. use fixed fac scale 91.118 = muR_ref_fixed ! fixed ren reference scale 91.118 = muF1_ref_fixed ! fixed fact reference scale for pdf1 91.118 = muF2_ref_fixed ! fixed fact reference scale for pdf2 -1 = dynamical_scale_choice ! Choose one of the preselected dynamical choices #*********************************************************************** # Renormalization and factorization scales (advanced and NLO options) * #*********************************************************************** False = fixed_QES_scale ! if .true. use fixed Ellis-Sexton scale 91.118 = QES_ref_fixed ! fixed Ellis-Sexton reference scale 1.0 = muR_over_ref ! ratio of current muR over reference muR 1.0 = muF1_over_ref ! ratio of current muF1 over reference muF1 1.0 = muF2_over_ref ! ratio of current muF2 over reference muF2 1.0 = QES_over_ref ! ratio of current QES over reference QES #*********************************************************************** # Reweight flags to get scale dependence and PDF uncertainty * # For scale dependence: factor rw_scale_up/down around central scale * # For PDF uncertainty: use LHAPDF with supported set * #*********************************************************************** True = reweight_scale ! reweight to get scale dependence 0.5 = rw_Rscale_down ! lower bound for ren scale variations 2.0 = rw_Rscale_up ! upper bound for ren scale variations 0.5 = rw_Fscale_down ! lower bound for fact scale variations 2.0 = rw_Fscale_up ! upper bound for fact scale variations False = reweight_PDF ! reweight to get PDF uncertainty 244601 = PDF_set_min ! First of the error PDF sets 244700 = PDF_set_max ! Last of the error PDF sets #*********************************************************************** # ickkw parameter: * # 0: No merging * # 3: FxFx Merging - WARNING! Applies merging only at the hard-event * # level. After showering an MLM-type merging should be applied as * # well. See http://amcatnlo.cern.ch/FxFx_merging.htm for details. * # 4: UNLOPS merging (with pythia8 only). No interface from within * # MG5_aMC available, but available in Pythia8. * # -1: NNLL+NLO jet-veto computation. See arxiv:1412.8408 [hep-ph]. * #*********************************************************************** 0 = ickkw #*********************************************************************** # #*********************************************************************** # BW cutoff (M+/-bwcutoff*Gamma) * #*********************************************************************** 15.0 = bwcutoff #*********************************************************************** # Cuts on the jets * # Jet clustering is performed by FastJet. # When matching to a parton shower, these generation cuts should be * # considerably softer than the analysis cuts. * # (more specific cuts can be specified in SubProcesses/cuts.f) * #*********************************************************************** 1.0 = jetalgo ! FastJet jet algorithm (1=kT, 0=C/A, -1=anti-kT) 0.7 = jetradius ! The radius parameter for the jet algorithm 10.0 = ptj ! Min jet transverse momentum -1.0 = etaj ! Max jet abs(pseudo-rap) (a value .lt.0 means no cut) #*********************************************************************** # Cuts on the charged leptons (e+, e-, mu+, mu-, tau+ and tau-) * # (more specific gen cuts can be specified in SubProcesses/cuts.f) * #*********************************************************************** 20.0 = ptl ! Min lepton transverse momentum -1.0 = etal ! Max lepton abs(pseudo-rap) (a value .lt.0 means no cut) 0.0 = drll ! Min distance between opposite sign lepton pairs 0.0 = drll_sf ! Min distance between opp. sign same-flavor lepton pairs 0.0 = mll ! Min inv. mass of all opposite sign lepton pairs 30.0 = mll_sf ! Min inv. mass of all opp. sign same-flavor lepton pairs #*********************************************************************** # Photon-isolation cuts, according to hep-ph/9801442 * # When ptgmin=0, all the other parameters are ignored * #*********************************************************************** 20.0 = ptgmin ! Min photon transverse momentum -1.0 = etagamma ! Max photon abs(pseudo-rap) 0.1 = R0gamma ! Radius of isolation code 1.0 = xn ! n parameter of eq.(3.4) in hep-ph/9801442 1.0 = epsgamma ! epsilon_gamma parameter of eq.(3.4) in hep-ph/9801442 True = isoEM ! isolate photons from EM energy (photons and leptons) #*********************************************************************** # Maximal PDG code for quark to be considered a jet when applying cuts.* # At least all massless quarks of the model should be included here. * #*********************************************************************** 4 = maxjetflavor #*********************************************************************** # For aMCfast+APPLGRID use in PDF fitting (http://amcfast.hepforge.org)* #*********************************************************************** 0 = iappl ! aMCfast switch (0=OFF, 1=prepare APPLgrids, 2=fill grids) #***********************************************************************