mssm Run3 bbH card script

bin/Powheg/production/Run3/13p6TeV/bb_H_2HDM_13p6TeV/CardSpace/createCards.sh

@@ -0,0 +1,24 @@
+#!/bin/bash
+
+for mass in 50 60 65 70 75 80 85 90 95 100 105 110 115 120 125 130 135 140 160 180 200 250 300 350 400 450 500 600 700 800 900 1000 1100 1200 1400 1600 1800 2000 2300 2600 2900 3200 3500;
+do
+
+    for hdampmode in "" 1 2; 
+    do 
+
+        if (( $(echo $mass'<'145 | bc -l) )); 
+        then 
+            width=0.0041
+        elif (( $(echo $mass'<'605 | bc -l) )); then 
+            width=0.1
+        elif (( $(echo $mass'<'2005 | bc -l) )); then
+            width=1.0
+        else
+            width=2.0
+        fi
+
+        ./setup.sh $mass $width $hdampmode; 
+
+    done; 
+done;
+

bin/Powheg/production/Run3/13p6TeV/bb_H_2HDM_13p6TeV/CardSpace/setup.sh

@@ -0,0 +1,24 @@
+#!/bin/bash
+
+mass=$1
+width=$2
+hdamp=`awk "BEGIN {print ($mass+8.36)/4}"`
+if [ "$3" == 1 ]; then hdamp=`awk "BEGIN {print 2*($mass+8.36)/4}"`; fi
+if [ "$3" == 2 ]; then hdamp=`awk "BEGIN {print 0.5*($mass+8.36)/4}"`; fi
+
+workdir=$PWD/../
+template=$workdir/bbH_mssm_template/
+if [ "$3" == "" ]; then massdir=$workdir/m${mass}; mkdir -p $massdir; fi
+if [ "$3" == 1 ]; then massdir=$workdir/m${mass}_up; mkdir -p $massdir; fi
+if [ "$3" == 2 ]; then massdir=$workdir/m${mass}_down; mkdir -p $massdir; fi
+
+echo "Copying files:"
+# cp -v $template/pwg-rwl.dat $massdir
+cp -v $template/bbH.input-* $massdir
+# cp -v $template/JHUGen.input $massdir
+
+sed -i "s/XHMASSX/${mass}/g" $massdir/bbH.input-*
+sed -i "s/XHWIDTHX/${width}/g" $massdir/bbH.input-*
+sed -i "s/XHDAMPX/${hdamp}/g" $massdir/bbH.input-*
+
+mv $massdir/bbH.input-base $massdir/bbH.input

bin/Powheg/production/Run3/13p6TeV/bb_H_2HDM_13p6TeV/bbH_mssm_template/bbH.input-base

@@ -0,0 +1,86 @@
+numevts NEVENTS
+
+! Production parameters for ttH production 
+hmass XHMASSXd0      ! mass of the HIggs boson [GeV]
+hwidth XHWIDTHXd0  ! width of the Higgs boson [GeV]
+hdamp XHDAMPXd0     ! h parameter for hdamp function: 1/4 (mH+2mb)
+
+hdecaymode -1   ! -1 no decay                
+                ! 0 all decay channels open
+                ! 1-6 d dbar, u ubar,..., t tbar 
+                ! 7-9 e+ e-, mu+ mu-, tau+ tau-
+                ! 10  W+W-
+                ! 11  ZZ
+                ! 12  gamma gamma     
+
+#pythiatune 320	! PYTHIA tune
+
+ih1   1           ! hadron 1 (1 for protons, -1 for antiprotons)
+ih2   1           ! hadron 2 (1 for protons, -1 for antiprotons)
+ebeam1 6800d0     ! energy of beam 1
+ebeam2 6800d0     ! energy of beam 2
+
+! To be set only if using 'native' pdfs
+!ndns1 131         ! pdf set for hadron 1 (mlm numbering)
+!ndns2 131         ! pdf set for hadron 2 (mlm numbering)
+! To be set only if using LHA pdfs
+! 10050 cteq6m
+! 10042 cteq6l with LO alphas (corr. to vbfnlo LO set)
+! 10041 cteq6l with NLO alphas
+! 21000 MSTW2008 LO
+! 21100 MSTW2008 NLO
+! 23300 MSTW2008 NLO, NF=4
+! 11082 CT10 Nf4
+! 92000 PDF4LHC15 NLO, NF=4 
+
+! lhans1  325500      ! pdf set for hadron 1 (LHA numbering) origin 260000
+! lhans2  325500      ! pdf set for hadron 2 (LHA numbering) origin 260000
+
+lhans1  92000      ! pdf set for hadron 1 (LHA numbering) origin 260000
+lhans2  92000      ! pdf set for hadron 2 (LHA numbering) origin 260000
+
+bornktmin 1d0    ! Generation cut; minimum kt in underlying born
+#bornsuppfact 1d0 ! (default 0d0) switch for Born suppression factor. If < 0 suppfact = 1
+delta_mbbmin 0d0 ! generation cut on mbb (no gencut used, if set to 0)
+
+use-old-grid    1 ! if 1 use old grid if file pwggrids.dat is present (<> 1 regenerate)
+use-old-ubound  1 ! if 1 use norm of upper bounding function stored in pwgubound.dat, if present; <> 1 regenerate
+
+ncall1  25000  ! number of calls for initializing the integration grid
+itmx1    2     ! number of iterations for initializing the integration grid
+ncall2  5000  ! number of calls for computing the integral and finding upper bound
+itmx2     5     ! number of iterations for computing the integral and finding upper bound
+foldcsi   2    ! number of folds on csi integration
+foldy     5    ! number of folds on  y  integration
+foldphi   2    ! number of folds on phi integration
+nubound 10000  ! number of bbarra calls to setup norm of upper bounding function
+icsimax  1     ! <= 100, number of csi subdivision when computing the upper bounds
+iymax    1     ! <= 100, number of y subdivision when computing the upper bounds
+xupbound 2d0   ! increase upper bound for radiation generation
+
+runningscales 3	   ! default 0 (no running scales); 1,2,3: use running scales
+btlscalereal 1d0   ! use real kinematics for scales, if appropriate
+renscfact 1d0   ! (default 1d0) ren scale factor: muren  = muref * renscfact 
+facscfact 1d0   ! (default 1d0) fac scale factor: mufact = muref * facscfact 
+testplots  1       ! (default 0, do not) do NLO and PWHG distributions
+
+#bornonly   1      ! (default 0) if 1 do Born only
+#fakevirt   1      ! (default 0) if 1 use Born for virtuals
+
+iseed    SEED
+maxseeds 9999
+manyseeds  0       ! Used to perform multiple runs with different random
+                   ! seeds in the same directory.
+                   ! If set to 1, the program asks for an integer j;
+                   ! The file pwgseeds.dat at line j is read, and the
+                   ! integer at line j is used to initialize the random
+                   ! sequence for the generation of the event.
+                   ! The event file is called pwgevents-'j'.lhe
+
+storeinfo_rwgt 1   ! store info for reweighting events 
+compute_rwgt 0     ! store info for reweighting events 
+pdfreweight 1
+
+toploop 1	 ! 1 include top loop effects; 0 disregard top loop effects
+msbar  1	 ! 1 MSbar renormalization; 0 On-Shell renormalization
+withnegweights 1 ! default 0