Merge branch '21.6_superchic_datacard2' into '21.6'

21.6 superchic datacard2

See merge request !62990

Generators/Superchic_i/python/DatRepository.py

+#########################################################################
+# File containing default Superchic dat from the source tree which are not deployed to Genser
+######################################################################
+
+import os
+
+def param_card_dat():
+
+    return """
+    ######################################################################
+    ## 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 MASS
+    ###################################
+    Block mass 
+        5 4.700000e+00 # MB 
+        6 1.730000e+02 # MT 
+    15 1.777000e+00 # MTA 
+    23 9.118800e+01 # MZ 
+    25 1.250000e+02 # MH 
+    ## 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.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.166390e-05 # Gf 
+        3 1.180000e-01 # aS (Note that Parameter not used if you use a PDF set) 
+
+    ###################################
+    ## 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.441404e+00 # WZ 
+    DECAY  24 2.047600e+00 # WW 
+    DECAY  25 6.382339e-03 # WH 
+    ## 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.
+    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 
+    """
+
+
+def writeTextToFile(fname, text):
+    """
+    Generic text dumper to file
+    """
+
+    try:
+        with  open(fname, 'w') as f:
+            for line in text:
+                f.write(line)
+    except:
+        raise Exception("File {} is not writable".format(fname))
+
+
+def writeDatRepository(directory):
+
+    if not os.path.exists(directory):
+        os.makedirs(directory)
+
+    writeTextToFile(os.path.join(directory, 'param_card.dat'), param_card_dat())
+

Generators/Superchic_i/python/SuperChicUtils.py

@@ -2,6 +2,8 @@ import subprocess, os, shlex, re
 
 from AthenaCommon import Logging
 
+from DatRepository import writeDatRepository
+
 ## Get handle to Athena logging
 logger = Logging.logging.getLogger("Superchic_i")
 
@@ -371,6 +373,7 @@ def SuperChicRun(Init, genSeq):
     genSeq.SuperChicConfig = Init
     
     writeInputDAT(Init)
+    writeDatRepository('Cards')
     SuperChicInitialize(Init)
     SuperChicExecute(Init)
     

Generators/Superchic_i/share/jobOptions.SuperChicPy8_yyMuMu_DD_surv.py

+evgenConfig.description = "SuperChic4 MC gamma + gamma pp collisions at 13000 GeV to 2 muons with dissociation on single side"
+evgenConfig.keywords = ["2photon","2muon","dissociation"]
+evgenConfig.contact = ["william.patrick.mccormack.iii@cern.ch", "oldrich.kepka@cern.ch"]
+
+from Superchic_i.SuperChicUtils import SuperChicConfig, SuperChicRun
+
+#class with the superchic initialization parameters.  Please see SuperChicUtils for a complete list of tunable parameters.
+scConfig = SuperChicConfig(runArgs)
+
+scConfig.isurv = 4                  # Model of soft survival (from 1 -> 4, corresponding to arXiv:1306.2149)
+scConfig.PDFname = 'MMHT2014lo68cl' # PDF set name
+scConfig.PDFmember = 0              # PDF member
+scConfig.proc = 57                  # Process number (59 = gg->gg, 56 = gg->ee, 68 = gg->a->gg ); Please consult Superchic Manual https://superchic.hepforge.org/
+scConfig.beam = 'prot'              # Beam type ('prot', 'ion')
+scConfig.sfaci = True             # Include soft survival effects
+scConfig.diff = 'dd'                # interaction: elastic ('el'), single ('sd','sda','sdb') and double ('dd') dissociation.
+scConfig.genunw  = True
+scConfig.ymin  = -5.0              # Minimum object rapidity Y_X
+scConfig.ymax  = 5.0               # Maximum object rapidity Y_X
+scConfig.mmin  = 20                 # Minimum object mass M_X
+scConfig.mmax  = 2000               # Maximum object mass M_X
+scConfig.gencuts  = True           # Generate cuts below
+scConfig.ptxmax  = 100         #cut on proton pt
+scConfig.ptamin  = 10.0             # Minimum pT of outgoing object a 
+scConfig.ptbmin  = 10.             # Minimum pT of outgoing object b 
+scConfig.etaamin  = -2.5           # Minimum eta of outgoing object a
+scConfig.etaamax   = 2.5           # Maximum eta of outgoing object a
+scConfig.etabmin  = -2.5           # Minimum eta of outgoing object b
+scConfig.etabmax   = 2.5           # Maximum eta of outgoing object b
+scConfig.acoabmax  = 100
+
+
+SuperChicRun(scConfig, genSeq)
+include('Superchic_i/Pythia8_DD_Common.py')