Interference gen config

Configurations/HWWSemiLepHighMass/Gen2016_v7/config_medmass.py

@@ -0,0 +1,42 @@
+# Configuration file to produce initial root files
+
+treeName = 'Events'
+
+tag = 'hm_gen'
+
+# used by mkShape to define output directory for root files
+outputDir = 'rootFile_vs'
+
+# file with TTree aliases
+aliasesFile = 'aliases_medmass.py'
+
+# file with list of variables
+variablesFile = 'variables_medmass.py'
+# variablesFile = 'variablesPlot.py'
+# variablesFile = 'variablesLimits.py'
+
+# file with list of cuts
+cutsFile = 'cuts_medmass.py'
+
+# file with list of samples
+samplesFile = 'samples_medmass.py'
+
+# file with list of samples
+plotFile = 'plot_medmass.py'
+
+# luminosity to normalize to (in 1/fb)
+lumi = 35.9
+
+# used by mkPlot to define output directory for plots
+# different from "outputDir" to do things more tidy
+outputDirPlots = 'plot_last'
+
+
+# used by mkDatacards to define output directory for datacards
+outputDirDatacard = 'Datacards_last'
+
+# structure file for datacard
+structureFile = 'structure_medmass.py'
+
+# nuisances file for mkDatacards and for mkShape
+#nuisancesFile = 'nuisances_medmass.py'

Configurations/HWWSemiLepHighMass/Gen2016_v7/cuts_medmass.py

@@ -0,0 +1,39 @@
+# cuts
+cuts = {}
+# this should be checked in postprocessing, just to be sure
+supercut = '\
+    LeptonGen_pt[0]>25 \
+&&    (Alt$(LeptonGen_pt[1],0)<=10 \
+&&  abs(LeptonGen_eta[0])<2.5) \
+'
+#&& GenEvtFlag>0 \
+'''
+cuts['ggH'] = {
+    'expr': 'mjjGen_max < 150.',
+    'categories': {
+      'boosted': 'GenEvtFlag==1',
+      'resolved': 'GenEvtFlag==2'  
+    }
+}
+
+cuts['qqH'] = {
+    'expr': 'mjjGen_max > 150.',
+    'categories': {
+      'boosted': 'GenEvtFlag==1',
+      'resolved': 'GenEvtFlag==2'
+    }
+
+}
+
+
+cuts['all'] = {
+    'expr': '1.',
+    'categories': {
+      'boosted': 'GenEvtFlag==1',
+      'resolved': 'GenEvtFlag==2'
+    }
+
+}
+'''
+
+cuts['supercut'] = {'expr': '1.'}

Configurations/HWWSemiLepHighMass/Gen2016_v7/massesAndModels.py

@@ -0,0 +1,23 @@
+massggh = ['115', '120', '124', '125', '126', '130', '135', '140', '145', '150', '155', '160', '165', '170', '175', '180', '190', '200', '210', '230', '250', '270', '300', '350', '400', '450', '500', '550', '600', '650', '700', '750', '800', '900', '1000', '1500', '2000', '2500', '3000', '4000', '5000']
+
+massvbf = ['115', '120', '124', '125', '126', '130', '135', '140', '145', '150', '155', '160', '165', '170', '175', '180', '190', '200', '210', '230', '250', '270', '300', '350', '400', '450', '500', '550', '600', '650', '700', '750', '800', '900', '1000', '1500', '2000', '2500', '3000', '4000', '5000']
+
+massggh = ['500', '2000']  #'200', '300', '400', 
+massvbf = ['500', '2000'] #'200', '300', '400', 
+#massvbf = ['500'] 
+#massggh = ['500'] 
+
+models = ['RelW10.0', 'RelW100.0']
+'''
+model_name = '_'+model.replace('.','')
+model_I = model+'_I'
+model_I_H = model+'_I_Honly'
+model_I_B = model+'_I_Bonly'
+model_B = model+'_B'
+model_H = model+'_H'
+'''
+
+#plotmasses = ['500', '800', '1000', '1500' ,'2000', '3000']
+plotmasses = ['500']
+plotmodels = ['RelW1000']
+#plotmasses = ['500', '1000', '3000']