Merge pull request #3545 from zhiyuanlcern/BSM-Htautau-Signal

BSM Htautau signal gridpacks production code

bin/Powheg/production/Run3/13p6TeV/gg_H_2HDM_13p6TeV/CardSpace/createCard.sh

@@ -0,0 +1,24 @@
+#!/bin/bash
+
+for mass in 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
+        for hfactmode 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  1: mass   2: channel: incl, t, b or tb # Combinations are (contrib-scale): incl-t/b/tb, t-t, b-b, tb-tb, t-tb, b-tb (!!!now t only)  3. width depending on mass  4.hfactmode
+            ./setup.sh $mass t $width $hfactmode; 
+        done;
+    done;
+done;

bin/Powheg/production/Run3/13p6TeV/gg_H_2HDM_13p6TeV/CardSpace/scales-higgs-mass-scan.dat

@@ -0,0 +1,97 @@
+60	33	8	4
+65	35	9	5
+70	36	10	5
+75	37	11	6
+80	38	12	6
+85	39	13	6
+90	40	13	6
+95	41  	14	7
+100	42	14	7
+105	43	15	8
+110	44	16	8
+115	46	17	8
+120	47	17	8
+125	48	18	9
+130	50	18	9
+135	51	19	10
+140	52	19	10
+150	55	21	11
+160	58	22	11
+170	61	23	12
+180	64	24	13
+190	67	26	14
+200	71	27	14
+210	75	28	15
+220	80	29	16
+230	92	30	17
+240	103	31	18
+250	108	32	18
+260	111	34	19
+270	112	35	20
+280	112	36	21
+290	112	37	22
+300	111	38	23
+310	108	39	23
+320	105	40	24
+330	101	41	25
+340	95	42	26
+350	87	43	26
+360	82	44	25
+370	82	46	25
+380	81	47	24
+390	81	48	24
+400	81	49	23
+410	82	50	23
+420	83	51	23
+430	85	52	22
+440	86	53	21
+450	87	53	21
+460	89	54	20
+470	91	55	19
+480	92	56	19
+490	94	57	18
+500	96	58	17
+510	97	59	16
+520	99	60	15
+530	101	61	14
+540	102	62	12
+550	104	63	11
+560	106	64	10
+570	107	65	8
+580	109	66	5
+590	111	67	2
+600	113	68	6
+610	115	69	8
+620	117	70	10
+630	119	71	12
+640	121	72	14
+650	127	73	16
+660	129	74	17
+670	131	75	19
+680	133	76	21
+690	135	77	22
+700	137	78	24
+710	139	79	25
+720	141	79	27
+730	143	80	28
+740	146	81	30
+750	148	82	32
+760	150	83	33
+770	152	84	35
+780	154	85	44
+790	156	86	45
+800	158	87	46
+900	180	95	65
+1000	202	103	279
+1100	223	113	277    
+1200	243	122	275
+1400	284	138	282
+1500	305	147	287
+1600	324	154	292
+1800	362	171	305
+2000	398	189	318
+2300	451	213	340
+2600	502	237	358
+2900	547	260	381
+3200	589	282	404
+3500	640	306	425

bin/Powheg/production/Run3/13p6TeV/gg_H_2HDM_13p6TeV/CardSpace/scripts/scales-higgs-mass-scan.dat

@@ -0,0 +1 @@
+../scales-higgs-mass-scan.dat

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

@@ -0,0 +1,71 @@
+#!/bin/bash
+
+export VO_CMS_SW_DIR=/cvmfs/cms.cern.ch
+source $VO_CMS_SW_DIR/cmsset_default.sh
+
+mass=${1}
+contrib=${2} # t-only
+width=${3}
+hdampmode=${4}
+hfactmode=${5}
+
+# set hfact depending on mass and contrib
+#e.g, for h125 t = 48, b = 18, t+b = 9
+hfact=9
+
+inline=$(cat scripts/scales-higgs-mass-scan.dat | grep -P "^${mass}\t")
+if [ "$inline" == "" ]; then
+  echo "Error: There are no hfact damping scales for your chosen mass point, exiting."
+  exit
+fi
+stringarray=($inline)
+if [ "${contrib}" == "t" ]; then
+  hfact=${stringarray[1]}
+elif [ "${contrib}" == "b" ]; then
+  hfact=${stringarray[2]}
+elif [ "${contrib}" == "tb" ]; then
+  hfact=${stringarray[3]}
+else
+  echo "Error: Your chosen contribution must can only be t, b, tb, exiting"
+  exit
+fi
+
+if [ "${hfactmode}" == 1 ]; then hfact=`awk "BEGIN {print 2*${hfact}}"`; fi
+if [ "${hfactmode}" == 2 ]; then hfact=`awk "BEGIN {print 0.5*${hfact}}"`; fi
+
+echo "hfact set to: " $hfact
+
+workdir=$PWD/../
+template=$workdir/gg_H_2HDM_template/
+massdir=$workdir/m${mass}_${contrib}
+
+if [[ "${hfactmode}" == 1 ]]; then massdir=$workdir/m${mass}_${contrib}_hfactUp; fi
+if [[ "${hfactmode}" == 2 ]]; then massdir=$workdir/m${mass}_${contrib}_hfactDown; fi
+
+echo "creating directories"
+mkdir -p $massdir
+echo "Copying files:"
+cp -v $template/pwg-rwl.dat $massdir
+cp -v $template/ggH_mssm.input-* $massdir
+cp -v $template/JHUGen.input $massdir
+cp -rv $template/lib $massdir
+
+sed -i "s/XHMASSX/${mass}/g" $massdir/ggH_mssm.input-*
+sed -i "s/XHWIDTHX/${width}/g" $massdir/ggH_mssm.input-*
+sed -i "s/XHFACTX/${hfact}/g" $massdir/ggH_mssm.input-*
+sed -i "s/XHDAMPX/${hdamp}/g" $massdir/ggH_mssm.input-*
+
+# tanb = 15 for h and A
+# tanb = 50 for H
+# alpha always pi/4
+
+tanb=15
+echo tanb set to ${tanb}
+
+sed -i "s/XTANBX/${tanb}/g" $massdir/ggH_mssm.input-*
+sed -i "s/XALPHAX/0.785398163397448/g" $massdir/ggH_mssm.input-*
+
+sed -i "s/XNOBOTX/1/g" $massdir/ggH_mssm.input-*
+sed -i "s/XNOTOPX/0/g" $massdir/ggH_mssm.input-*
+
+mv $massdir/ggH_mssm.input-base $massdir/ggH_mssm.input

bin/Powheg/production/Run3/13p6TeV/gg_H_2HDM_13p6TeV/gg_H_2HDM_template/ggH_mssm.input-base

@@ -0,0 +1,103 @@
+numevts NEVENTS   ! number of events to be generated
+ih1   1           ! hadron 1 (1 for protons, -1 for antiprotons)
+ih2   1           ! hadron 2 (1 for protons, -1 for antiprotons)
+ebeam1 6800     ! energy of beam 1
+ebeam2 6800    ! energy of beam 2
+
+! To be set only if using LHA pdfs
+lhans1 325300       ! pdf set for hadron 1 (LHA numbering)
+lhans2 325300       ! pdf set for hadron 2 (LHA numbering)
+
+! To be set only if using different pdf sets for the two incoming hadrons
+! QCDLambda5  0.25 ! for not equal pdf sets 
+
+! Parameters to allow or not the use of stored data
+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  100000  ! number of calls for initializing the integration grid
+itmx1    1     ! number of iterations for initializing the integration grid
+ncall2  100000  ! number of calls for computing the integral and finding upper bound
+itmx2     5    ! number of iterations for computing the integral and finding upper bound
+nubound 200000  ! number of bbarra calls to setup norm of upper bounding function
+foldcsi   2    ! number of folds on csi integration
+foldy     5    ! number of folds on  y  integration
+foldphi   2    ! number of folds on phi integration
+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
+
+! OPTIONAL PARAMETERS
+
+# renscfact  1   ! (default 1d0) ren scale factor: muren  = muref * renscfact
+# facscfact  1   ! (default 1d0) fac scale factor: mufact = muref * facscfact
+# ptsupp     0d0   ! (default 0d0)  mass param for Born suppression factor (generation cut) If < 0 suppfact = 1
+# bornonly   1      ! (default 0) if 1 do Born only
+# smartsig   0      ! (default 1) remember equal amplitudes (0 do not remember)
+# withsubtr  0      ! (default 1) subtract real counterterms (0 do not subtract)
+# withdamp    1      ! (default 0, do not use) use Born-zero damping factor
+# ptsqmin    0.8   ! (default 0.8 GeV) minimum pt for generation of radiation
+# charmthr   1.5    ! (default 1.5 GeV) charm treshold for gluon splitting 
+# bottomthr  5.0    ! (default 5.0 GeV) bottom treshold for gluon splitting
+testplots  1      ! (default 0, do not) do NLO and PWHG distributions
+# testsuda  1       ! (default 0, do not test) test Sudakov form factor
+# radregion 1       ! (default all regions) only generate radiation in the selected singular region  
+# charmthrpdf  1.5  ! (default 1.5 GeV) pdf charm treshold
+# bottomthrpdf 5.0  ! (default 5.0 GeV) pdf bottom treshold
+withnegweights 1 ! use negative weights
+
+runningscale 1    ! 0 = scales equal to the Higgs pole mass; 1  = scales equal to the Higgs virtuality; 
+                  ! 2 = scales equal to the Higgs pole mass for Born-like configuration and to the transverse mass for real emission contribution
+# When using 2 uncomment the following option
+# btlscalereal 1
+
+iseed    SEED      ! initialize random number sequence 
+maxseeds 4999
+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
+
+# iupperisr 1 ! (default 1) choice of ISR upper bounding functional form 
+# iupperfsr 2 ! (default 2) choice of FSR upper bounding functional form 
+
+! GGF_H production:
+! **** Mandatory parameters for ALL models ****
+zerowidth 1          ! Control if the Higgs boson is to be produced on-shell or not: 1 = On-Shell; 0 = Off-shell with Breit-Wigner
+bwshape 1	     ! BW shape; 1 = running width; 0 = fixed width
+ew 0                 ! ew = 0 disable EW corrections - ew = 1 enable EW corrections
+fastew 1	     ! mass sampling of ew corrections
+gfermi 0.116637D-04  ! GF
+hdecaymode -1   ! PDG code for first decay product of the higgs
+                ! -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     
+masswindow 10d0      !(default 10d0) number of widths around hmass in the BW for an off-shell Higgs boson
+scheme 0	     ! Mass renormalization scheme: 0 = OS (Default); 1 = MSBAR; 2 = DRBAR
+topmass 172.5       ! top quark mass
+bottommass 4.75d0    ! bottom quark mass - if defined it enables the bottom quark
+!charmmass 1.5d0     ! char quark mass - if defined it enables the charm quark
+! Optional
+hdecaywidth 0        ! If equals to 1 read total decay width from HDECAY sm.br2 file
+! 2HDM parameters
+2HDMtype 2
+alpha XALPHAXd0      ! 0.785398163397448d0
+
+higgstype 1            ! 1:h, 2:H, 3:A
+hmass XHMASSXd0
+hwidth XHWIDTHXd0
+tanb XTANBXd0
+hfact XHFACTX
+nobot XNOBOTX
+notop XNOTOPX
+
+storeinfo_rwgt 1   ! store info for reweighting events
+compute_rwgt 0     ! store info for reweighting events
+pdfreweight 1

bin/Powheg/production/Run3/13p6TeV/gg_H_2HDM_13p6TeV/gg_H_2HDM_template/lib/libchaplin.la

@@ -0,0 +1,41 @@
+# libchaplin.la - a libtool library file
+# Generated by ltmain.sh (GNU libtool) 2.2.6b Debian-2.2.6b-2ubuntu3
+#
+# Please DO NOT delete this file!
+# It is necessary for linking the library.
+
+# The name that we can dlopen(3).
+dlname='libchaplin.so.0'
+
+# Names of this library.
+library_names='libchaplin.so.0.0.0 libchaplin.so.0 libchaplin.so'
+
+# The name of the static archive.
+old_library='libchaplin.a'
+
+# Linker flags that can not go in dependency_libs.
+inherited_linker_flags=''
+
+# Libraries that this one depends upon.
+dependency_libs=''
+
+# Names of additional weak libraries provided by this library
+weak_library_names=''
+
+# Version information for libchaplin.
+current=0
+age=0
+revision=0
+
+# Is this an already installed library?
+installed=yes
+
+# Should we warn about portability when linking against -modules?
+shouldnotlink=no
+
+# Files to dlopen/dlpreopen
+dlopen=''
+dlpreopen=''
+
+# Directory that this library needs to be installed in:
+libdir='/vols/cms/dw515/ggh_mc/ggh-mssm/m125/CMSSW_10_2_3/src/powheg_new/CMSSW_9_3_0/src/genproductions/bin/Powheg/my_ggH/POWHEG-BOX/gg_H_2HDM/chaplin-1.2/../lib'

bin/Powheg/production/Run3/13p6TeV/gg_H_2HDM_13p6TeV/gg_H_2HDM_template/pwg-rwl.dat

@@ -0,0 +1,49 @@
+<initrwgt>
+<weightgroup name='scale_variation' combine='envelope' >
+<weight id='1001'> lhapdf=90900 renscfact=1d0 facscfact=1d0 </weight>
+<weight id='1002'> lhapdf=90900 renscfact=1d0 facscfact=2d0 </weight>
+<weight id='1003'> lhapdf=90900 renscfact=1d0 facscfact=0.5d0 </weight>
+<weight id='1004'> lhapdf=90900 renscfact=2d0 facscfact=1d0 </weight>
+<weight id='1005'> lhapdf=90900 renscfact=2d0 facscfact=2d0 </weight>
+<weight id='1006'> lhapdf=90900 renscfact=2d0 facscfact=0.5d0 </weight>
+<weight id='1007'> lhapdf=90900 renscfact=0.5d0 facscfact=1d0 </weight>
+<weight id='1008'> lhapdf=90900 renscfact=0.5d0 facscfact=2d0 </weight>
+<weight id='1009'> lhapdf=90900 renscfact=0.5d0 facscfact=0.5d0 </weight>
+<weight id='1010'> lhapdf=90900 renscfact=0.25d0 facscfact=0.25d0 </weight>
+<weight id='1011'> lhapdf=90900 renscfact=0.5d0 facscfact=0.25d0 </weight>
+<weight id='1012'> lhapdf=90900 renscfact=0.25d0 facscfact=0.5d0 </weight>
+</weightgroup>
+<weightgroup name='PDF_variation1' combine='hessian' >
+<weight id='10000'> lhapdf=90900 </weight>
+<weight id='10001'> lhapdf=90901 </weight>
+<weight id='10002'> lhapdf=90902 </weight>
+<weight id='10003'> lhapdf=90903 </weight>
+<weight id='10004'> lhapdf=90904 </weight>
+<weight id='10005'> lhapdf=90905 </weight>
+<weight id='10006'> lhapdf=90906 </weight>
+<weight id='10007'> lhapdf=90907 </weight>
+<weight id='10008'> lhapdf=90908 </weight>
+<weight id='10009'> lhapdf=90909 </weight>
+<weight id='10010'> lhapdf=90910 </weight>
+<weight id='10011'> lhapdf=90911 </weight>
+<weight id='10012'> lhapdf=90912 </weight>
+<weight id='10013'> lhapdf=90913 </weight>
+<weight id='10014'> lhapdf=90914 </weight>
+<weight id='10015'> lhapdf=90915 </weight>
+<weight id='10016'> lhapdf=90916 </weight>
+<weight id='10017'> lhapdf=90917 </weight>
+<weight id='10018'> lhapdf=90918 </weight>
+<weight id='10019'> lhapdf=90919 </weight>
+<weight id='10020'> lhapdf=90920 </weight>
+<weight id='10021'> lhapdf=90921 </weight>
+<weight id='10022'> lhapdf=90922 </weight>
+<weight id='10023'> lhapdf=90923 </weight>
+<weight id='10024'> lhapdf=90924 </weight>
+<weight id='10025'> lhapdf=90925 </weight>
+<weight id='10026'> lhapdf=90926 </weight>
+<weight id='10027'> lhapdf=90927 </weight>
+<weight id='10028'> lhapdf=90928 </weight>
+<weight id='10029'> lhapdf=90929 </weight>
+<weight id='10030'> lhapdf=90930 </weight>
+</weightgroup>
+</initrwgt>