diff --git a/700xxx/700398/mc.Sh_2211_mumugamma.py b/700xxx/700398/mc.Sh_2211_mumugamma.py
new file mode 100644
index 0000000000000000000000000000000000000000..41e127c39dc9fa8236443476e8584325e6bce81e
--- /dev/null
+++ b/700xxx/700398/mc.Sh_2211_mumugamma.py
@@ -0,0 +1,67 @@
+include("Sherpa_i/Base_Fragment.py")
+include("Sherpa_i/NNPDF30NNLO.py")
+include("Sherpa_i/EW_scheme_sinthetaW_mZ.py")
+
+evgenConfig.description = "Sherpa mumugamma + 0,1j@NLO + 2,3,4j@LO with 7<pT_y"
+evgenConfig.keywords = ["SM", "muon", "neutrino", "photon", "NLO" ]
+evgenConfig.contact  = [ "atlas-generators-sherpa@cern.ch", "matthew.gignac@cern.ch" ]
+evgenConfig.nEventsPerJob = 1000
+
+genSeq.Sherpa_i.RunCard="""
+(run){
+
+  %scales, tags for scale variations
+  FSF:=1.; RSF:=1.; QSF:=1.;
+  SCALES STRICT_METS{FSF*MU_F2}{RSF*MU_R2}{QSF*MU_Q2};
+
+  # Shower improvements
+  NLO_SUBTRACTION_SCHEME=2;
+
+  # me generator settings
+  ME_SIGNAL_GENERATOR Comix Amegic LOOPGEN;
+  LOOPGEN:=OpenLoops
+
+  # tags for process setup
+  NJET:=4; LJET:=3,4; QCUT:=20;
+
+  # EW corrections setup
+  ASSOCIATED_CONTRIBUTIONS_VARIATIONS=EW EW|LO1;
+
+  # speed and neg weight fraction improvements
+  PP_RS_SCALE VAR{sqr(sqrt(H_T2)-PPerp(p[2])-PPerp(p[3])+MPerp(p[2]+p[3]))/4};
+  METS_BBAR_MODE=5
+  NLO_CSS_PSMODE=1
+
+}(run)
+
+(processes){
+  Process 93 93 -> 13 -13 22 93{NJET}
+  Enhance_Function VAR{pow(max(PPerp(p[2]+p[3]),PPerp(p[4]))/20.0,2)}
+  Associated_Contributions EW|LO1 {LJET};
+  Order (*,3); CKKW sqr(QCUT/E_CMS)
+  NLO_QCD_Mode MC@NLO {LJET}
+  ME_Generator Amegic {LJET}
+  RS_ME_Generator Comix {LJET}
+  Loop_Generator LOOPGEN {LJET}
+  PSI_ItMin 20000 {4}
+  Integration_Error 0.99 {4}
+  PSI_ItMin 50000 {5,6,7}
+  Integration_Error 0.99 {5,6,7}
+  End process
+}(processes)
+
+(selector){
+  PTNLO  22  7  E_CMS
+  IsolationCut  22  0.1  2  0.10
+  DeltaRNLO  22  90  0.1 1000.0
+  Mass  90  90  10.0  E_CMS
+}(selector)
+"""
+
+genSeq.Sherpa_i.NCores = 128
+genSeq.Sherpa_i.OpenLoopsLibs = [ "pplla","ppllaj","pplla_ew","ppllaj_ew"]
+
+#factor alpha(0)/alpha(mZ) for the photon, using Sherpa default values
+#KFACTOR VAR{128.802/137.03599976};
+genSeq.Sherpa_i.CrossSectionScaleFactor = 0.939913601
+
diff --git a/700xxx/700398/mc_13TeV.Sh_2211_mumugamma.GRID.tar.gz b/700xxx/700398/mc_13TeV.Sh_2211_mumugamma.GRID.tar.gz
new file mode 120000
index 0000000000000000000000000000000000000000..30abe9f66d13834e8308b0848df4f4b7321700ef
--- /dev/null
+++ b/700xxx/700398/mc_13TeV.Sh_2211_mumugamma.GRID.tar.gz
@@ -0,0 +1 @@
+/eos/user/m/mgignac/mc/2211/Zmumuy_final_sin2ThetaW_newSeed/mc_13TeV.Sh_2211_mumugamma.GRID.tar.gz
\ No newline at end of file
diff --git a/700xxx/700399/mc.Sh_2211_eegamma.py b/700xxx/700399/mc.Sh_2211_eegamma.py
new file mode 100644
index 0000000000000000000000000000000000000000..b62e01363e28b41b6e0cf1565a8d7034b400b537
--- /dev/null
+++ b/700xxx/700399/mc.Sh_2211_eegamma.py
@@ -0,0 +1,66 @@
+include("Sherpa_i/Base_Fragment.py")
+include("Sherpa_i/NNPDF30NNLO.py")
+include("Sherpa_i/EW_scheme_sinthetaW_mZ.py")
+
+evgenConfig.description = "Sherpa mumugamma + 0,1j@NLO + 2,3,4j@LO with 7<pT_y"
+evgenConfig.keywords = ["SM", "electron", "neutrino", "photon", "NLO" ]
+evgenConfig.contact  = [ "atlas-generators-sherpa@cern.ch", "matthew.gignac@cern.ch" ]
+evgenConfig.nEventsPerJob = 1000
+
+genSeq.Sherpa_i.RunCard="""
+(run){
+
+  %scales, tags for scale variations
+  FSF:=1.; RSF:=1.; QSF:=1.;
+  SCALES STRICT_METS{FSF*MU_F2}{RSF*MU_R2}{QSF*MU_Q2};
+
+  # Shower improvements
+  NLO_SUBTRACTION_SCHEME=2;
+
+  # me generator settings
+  ME_SIGNAL_GENERATOR Comix Amegic LOOPGEN;
+  LOOPGEN:=OpenLoops
+
+  # tags for process setup
+  NJET:=4; LJET:=3,4; QCUT:=20;
+
+  # EW corrections setup
+  ASSOCIATED_CONTRIBUTIONS_VARIATIONS=EW EW|LO1; 
+
+  # speed and neg weight fraction improvements
+  PP_RS_SCALE VAR{sqr(sqrt(H_T2)-PPerp(p[2])-PPerp(p[3])+MPerp(p[2]+p[3]))/4};
+  METS_BBAR_MODE=5
+  NLO_CSS_PSMODE=1
+
+}(run)
+
+(processes){
+  Process 93 93 -> 11 -11 22 93{NJET}
+  Enhance_Function VAR{pow(max(PPerp(p[2]+p[3]),PPerp(p[4]))/20.0,2)}
+  Associated_Contributions EW|LO1 {LJET};
+  Order (*,3); CKKW sqr(QCUT/E_CMS)
+  NLO_QCD_Mode MC@NLO {LJET}
+  ME_Generator Amegic {LJET}
+  RS_ME_Generator Comix {LJET}
+  Loop_Generator LOOPGEN {LJET}
+  PSI_ItMin 20000 {4}
+  Integration_Error 0.99 {4}
+  PSI_ItMin 50000 {5,6,7}
+  End process
+}(processes)
+
+(selector){
+  PTNLO  22  7  E_CMS
+  IsolationCut  22  0.1  2  0.10
+  DeltaRNLO  22  90  0.1 1000.0
+  Mass  90  90  10.0  E_CMS
+}(selector)
+"""
+
+genSeq.Sherpa_i.NCores = 128
+genSeq.Sherpa_i.OpenLoopsLibs = [ "pplla","ppllaj","pplla_ew","ppllaj_ew"]
+
+#factor alpha(0)/alpha(mZ) for the photon, using Sherpa default values
+#KFACTOR VAR{128.802/137.03599976};
+genSeq.Sherpa_i.CrossSectionScaleFactor = 0.939913601
+
diff --git a/700xxx/700399/mc_13TeV.Sh_2211_eegamma.GRID.tar.gz b/700xxx/700399/mc_13TeV.Sh_2211_eegamma.GRID.tar.gz
new file mode 120000
index 0000000000000000000000000000000000000000..ed88a64f09bdcd949d21caacd83ecef22aba3b28
--- /dev/null
+++ b/700xxx/700399/mc_13TeV.Sh_2211_eegamma.GRID.tar.gz
@@ -0,0 +1 @@
+/eos/user/m/mgignac/mc/2211/Vy_final_sin2ThetaW/mc_13TeV.Sh_2211_eegamma.GRID.tar.gz
\ No newline at end of file
diff --git a/700xxx/700400/mc.Sh_2211_tautaugamma.py b/700xxx/700400/mc.Sh_2211_tautaugamma.py
new file mode 100644
index 0000000000000000000000000000000000000000..096d31832ff5e9247b9318776e3ac22abaebab9c
--- /dev/null
+++ b/700xxx/700400/mc.Sh_2211_tautaugamma.py
@@ -0,0 +1,69 @@
+include("Sherpa_i/Base_Fragment.py")
+include("Sherpa_i/NNPDF30NNLO.py")
+include("Sherpa_i/EW_scheme_sinthetaW_mZ.py")
+
+evgenConfig.description = "Sherpa tautaugamma + 0,1j@NLO + 2,3,4j@LO with 7<pT_y"
+evgenConfig.keywords = ["SM", "tau", "neutrino", "photon", "NLO" ]
+evgenConfig.contact  = [ "atlas-generators-sherpa@cern.ch", "matthew.gignac@cern.ch" ]
+evgenConfig.nEventsPerJob = 500
+
+genSeq.Sherpa_i.RunCard="""
+(run){
+
+  %scales, tags for scale variations
+  FSF:=1.; RSF:=1.; QSF:=1.;
+  SCALES STRICT_METS{FSF*MU_F2}{RSF*MU_R2}{QSF*MU_Q2};
+
+  # Shower improvements
+  NLO_SUBTRACTION_SCHEME=2;
+
+  # me generator settings
+  ME_SIGNAL_GENERATOR Comix Amegic LOOPGEN;
+  LOOPGEN:=OpenLoops
+
+  # tags for process setup
+  NJET:=4; LJET:=3,4; QCUT:=20;
+
+  # EW corrections setup
+  ASSOCIATED_CONTRIBUTIONS_VARIATIONS=EW EW|LO1; 
+
+  SOFT_SPIN_CORRELATIONS=1;
+
+  # speed and neg weight fraction improvements
+  PP_RS_SCALE VAR{sqr(sqrt(H_T2)-PPerp(p[2])-PPerp(p[3])+MPerp(p[2]+p[3]))/4};
+  METS_BBAR_MODE=5
+  NLO_CSS_PSMODE=1
+
+}(run)
+
+(processes){
+  Process 93 93 -> 15 -15 22 93{NJET}
+  Enhance_Function VAR{pow(max(PPerp(p[2]+p[3]),PPerp(p[4]))/20.0,2)}
+  Associated_Contributions EW|LO1 {LJET};
+  Order (*,3); CKKW sqr(QCUT/E_CMS)
+  NLO_QCD_Mode MC@NLO {LJET}
+  ME_Generator Amegic {LJET}
+  RS_ME_Generator Comix {LJET}
+  Loop_Generator LOOPGEN {LJET}
+  PSI_ItMin 20000 {4}
+  Integration_Error 0.99 {4}
+  PSI_ItMin 50000 {5,6,7}
+  Integration_Error 0.99 {5,6,7}
+  End process
+}(processes)
+
+(selector){
+  PTNLO  22  7  E_CMS
+  IsolationCut  22  0.1  2  0.10
+  DeltaRNLO  22  90  0.1 1000.0
+  Mass  90  90  10.0  E_CMS
+}(selector)
+"""
+
+genSeq.Sherpa_i.NCores = 128
+genSeq.Sherpa_i.OpenLoopsLibs = [ "pplla","ppllaj","pplla_ew","ppllaj_ew"]
+
+#factor alpha(0)/alpha(mZ) for the photon, using Sherpa default values
+#KFACTOR VAR{128.802/137.03599976};
+genSeq.Sherpa_i.CrossSectionScaleFactor = 0.939913601
+
diff --git a/700xxx/700400/mc_13TeV.Sh_2211_tautaugamma.GRID.tar.gz b/700xxx/700400/mc_13TeV.Sh_2211_tautaugamma.GRID.tar.gz
new file mode 120000
index 0000000000000000000000000000000000000000..30a63551e00891172bda2b49d2e29dba321eaf9b
--- /dev/null
+++ b/700xxx/700400/mc_13TeV.Sh_2211_tautaugamma.GRID.tar.gz
@@ -0,0 +1 @@
+/eos/user/m/mgignac/mc/2211/Vy_final_sin2ThetaW/mc_13TeV.Sh_2211_tautaugamma.GRID.tar.gz
\ No newline at end of file
diff --git a/700xxx/700401/mc.Sh_2211_nunugamma.py b/700xxx/700401/mc.Sh_2211_nunugamma.py
new file mode 100644
index 0000000000000000000000000000000000000000..5d60f292218831fec747be3b13f406941b96e1da
--- /dev/null
+++ b/700xxx/700401/mc.Sh_2211_nunugamma.py
@@ -0,0 +1,62 @@
+include("Sherpa_i/Base_Fragment.py")
+include("Sherpa_i/NNPDF30NNLO.py")
+include("Sherpa_i/EW_scheme_sinthetaW_mZ.py")
+
+evgenConfig.description = "Sherpa mumugamma + 0,1j@NLO + 2,3,4j@LO with 7<pT_y"
+evgenConfig.keywords = ["SM", "neutrino", "photon", "NLO" ]
+evgenConfig.contact  = [ "atlas-generators-sherpa@cern.ch", "matthew.gignac@cern.ch" ]
+evgenConfig.nEventsPerJob = 1000
+
+genSeq.Sherpa_i.RunCard="""
+(run){
+
+  %scales, tags for scale variations
+  FSF:=1.; RSF:=1.; QSF:=1.;
+  SCALES STRICT_METS{FSF*MU_F2}{RSF*MU_R2}{QSF*MU_Q2};
+
+  # Shower improvements
+  NLO_SUBTRACTION_SCHEME=2;
+
+  # me generator settings
+  ME_SIGNAL_GENERATOR Comix Amegic LOOPGEN;
+  LOOPGEN:=OpenLoops
+
+  # tags for process setup
+  NJET:=4; LJET:=3,4; QCUT:=20;
+
+  # EW corrections setup
+  ASSOCIATED_CONTRIBUTIONS_VARIATIONS=EW EW|LO1; 
+  PP_RS_SCALE VAR{sqr(sqrt(H_T2)-PPerp(p[2])-PPerp(p[3])+MPerp(p[2]+p[3]))/4};
+  METS_BBAR_MODE=5
+  NLO_CSS_PSMODE=1
+
+}(run)
+
+(processes){
+  Process 93 93 -> 91 91 22 93{NJET}
+  Enhance_Function VAR{pow(max(PPerp(p[2]+p[3]),PPerp(p[4]))/20.0,2)}
+  Associated_Contributions EW|LO1 {LJET};
+  Order (*,3); CKKW sqr(QCUT/E_CMS)
+  NLO_QCD_Mode MC@NLO {LJET}
+  ME_Generator Amegic {LJET}
+  RS_ME_Generator Comix {LJET}
+  Loop_Generator LOOPGEN {LJET}
+  PSI_ItMin 20000 {4}
+  Integration_Error 0.99 {4}
+  PSI_ItMin 50000 {5,6,7}
+  Integration_Error 0.99 {5,6,7}
+  End process
+}(processes)
+
+(selector){
+  PTNLO  22  7  E_CMS
+  IsolationCut  22  0.1  2  0.10
+}(selector)
+"""
+
+genSeq.Sherpa_i.NCores = 128
+genSeq.Sherpa_i.OpenLoopsLibs = [ "pplla","ppllaj","pplla_ew","ppllaj_ew"]
+
+#factor alpha(0)/alpha(mZ) for the photon, using Sherpa default values
+#KFACTOR VAR{128.802/137.03599976};
+genSeq.Sherpa_i.CrossSectionScaleFactor = 0.939913601
diff --git a/700xxx/700401/mc_13TeV.Sh_2211_nunugamma.GRID.tar.gz b/700xxx/700401/mc_13TeV.Sh_2211_nunugamma.GRID.tar.gz
new file mode 120000
index 0000000000000000000000000000000000000000..b0fd3d46d725ba760770de6654fa644c0f01ad38
--- /dev/null
+++ b/700xxx/700401/mc_13TeV.Sh_2211_nunugamma.GRID.tar.gz
@@ -0,0 +1 @@
+/eos/user/m/mgignac/mc/2211/Vy_final_sin2ThetaW/mc_13TeV.Sh_2211_nunugamma.GRID.tar.gz
\ No newline at end of file
diff --git a/700xxx/700402/mc.Sh_2211_munugamma.py b/700xxx/700402/mc.Sh_2211_munugamma.py
new file mode 100644
index 0000000000000000000000000000000000000000..9329d4b0522fa6e8c42048af1d0b91e7fc3ef18a
--- /dev/null
+++ b/700xxx/700402/mc.Sh_2211_munugamma.py
@@ -0,0 +1,80 @@
+
+include("Sherpa_i/Base_Fragment.py")
+include("Sherpa_i/NNPDF30NNLO.py")
+
+evgenConfig.description = "Sherpa munugamma + 0,1j@NLO + 2,3,4j@LO with 7<pT_y"
+evgenConfig.keywords = ["SM", "muon", "neutrino", "photon", "NLO" ]
+evgenConfig.contact  = [ "atlas-generators-sherpa@cern.ch", "matthew.gignac@cern.ch" ]
+evgenConfig.nEventsPerJob = 1000
+
+genSeq.Sherpa_i.RunCard="""
+(run){
+
+  %scales, tags for scale variations
+  FSF:=1.; RSF:=1.; QSF:=1.;
+  SCALES STRICT_METS{FSF*MU_F2}{RSF*MU_R2}{QSF*MU_Q2};
+
+  # Shower improvements
+  NLO_SUBTRACTION_SCHEME=2;
+
+  # me generator settings
+  ME_SIGNAL_GENERATOR Comix Amegic LOOPGEN;
+  LOOPGEN:=OpenLoops
+
+  # tags for process setup
+  NJET:=4; LJET:=3,4; QCUT:=20;
+
+  # EW corrections setup
+  OL_PARAMETERS=ew_scheme 2 ew_renorm_scheme 1
+  ASSOCIATED_CONTRIBUTIONS_VARIATIONS=EW EW|LO1; 
+  EW_SCHEME=3
+  GF=1.166397e-5
+  METS_BBAR_MODE=5
+
+  # speed and neg weight fraction improvements
+  PP_RS_SCALE VAR{sqr(sqrt(H_T2)-PPerp(p[2])-PPerp(p[3])+MPerp(p[2]+p[3]))/4};
+  NLO_CSS_PSMODE=1
+
+}(run)
+
+(processes){
+  Process 93 93 -> 22 13 -14 93{NJET}
+  Enhance_Function VAR{pow(max(PPerp(p[2]+p[3]),PPerp(p[4]))/20.0,2)}
+  Associated_Contributions EW|LO1 {LJET};
+  Order (*,3); CKKW sqr(QCUT/E_CMS)
+  NLO_QCD_Mode MC@NLO {LJET}
+  ME_Generator Amegic {LJET}
+  RS_ME_Generator Comix {LJET}
+  Loop_Generator LOOPGEN {LJET}
+  PSI_ItMin 20000 {4}
+  Integration_Error 0.99 {4}
+  PSI_ItMin 50000 {5,6,7}
+  Integration_Error 0.99 {5,6,7}
+  End process
+
+  Process 93 93 -> 22 -13 14 93{NJET}
+  Enhance_Function VAR{pow(max(PPerp(p[2]+p[3]),PPerp(p[4]))/20.0,2)}
+  Associated_Contributions EW|LO1 {LJET};
+  Order (*,3); CKKW sqr(QCUT/E_CMS)
+  NLO_QCD_Mode MC@NLO {LJET}
+  ME_Generator Amegic {LJET}
+  RS_ME_Generator Comix {LJET}
+  Loop_Generator LOOPGEN {LJET}
+  PSI_ItMin 20000 {4}
+  Integration_Error 0.99 {4}
+  PSI_ItMin 50000 {5,6,7}
+  Integration_Error 0.99 {5,6,7}
+  End process
+}(processes)
+
+(selector){
+  PTNLO  22  7  E_CMS
+  IsolationCut  22  0.1  2  0.10
+  DeltaRNLO  22  90  0.1 1000.0
+  Mass  90  91  1.0  E_CMS
+}(selector)
+"""
+
+genSeq.Sherpa_i.NCores = 128
+genSeq.Sherpa_i.Parameters += [ "OL_PARAMETERS=ew_scheme=2 ew_renorm_scheme=1 write_parameters=1" ]
+genSeq.Sherpa_i.Parameters += [ "EW_SCHEME=3", "GF=1.166397e-5" ]
diff --git a/700xxx/700402/mc_13TeV.Sh_2211_munugamma.GRID.tar.gz b/700xxx/700402/mc_13TeV.Sh_2211_munugamma.GRID.tar.gz
new file mode 120000
index 0000000000000000000000000000000000000000..6b45725e6c66e4f51e1c263725bc8d5a23291a33
--- /dev/null
+++ b/700xxx/700402/mc_13TeV.Sh_2211_munugamma.GRID.tar.gz
@@ -0,0 +1 @@
+/eos/user/m/mgignac/mc/2211/Vy_final/mc_13TeV.Sh_2211_munugamma.GRID.tar.gz
\ No newline at end of file
diff --git a/700xxx/700403/mc.Sh_2211_enugamma.py b/700xxx/700403/mc.Sh_2211_enugamma.py
new file mode 100644
index 0000000000000000000000000000000000000000..d30d39a3145305e68fa55eae88db6d63e4be210d
--- /dev/null
+++ b/700xxx/700403/mc.Sh_2211_enugamma.py
@@ -0,0 +1,80 @@
+
+include("Sherpa_i/Base_Fragment.py")
+include("Sherpa_i/NNPDF30NNLO.py")
+
+evgenConfig.description = "Sherpa enugamma + 0,1j@NLO + 2,3,4j@LO with 7<pT_y"
+evgenConfig.keywords = ["SM", "electron", "neutrino", "photon", "NLO" ]
+evgenConfig.contact  = [ "atlas-generators-sherpa@cern.ch", "matthew.gignac@cern.ch" ]
+evgenConfig.nEventsPerJob = 1000
+
+genSeq.Sherpa_i.RunCard="""
+(run){
+
+  %scales, tags for scale variations
+  FSF:=1.; RSF:=1.; QSF:=1.;
+  SCALES STRICT_METS{FSF*MU_F2}{RSF*MU_R2}{QSF*MU_Q2};
+
+  # Shower improvements
+  NLO_SUBTRACTION_SCHEME=2;
+
+  # me generator settings
+  ME_SIGNAL_GENERATOR Comix Amegic LOOPGEN;
+  LOOPGEN:=OpenLoops
+
+  # tags for process setup
+  NJET:=4; LJET:=3,4; QCUT:=20;
+
+  # EW corrections setup
+  OL_PARAMETERS=ew_scheme 2 ew_renorm_scheme 1
+  ASSOCIATED_CONTRIBUTIONS_VARIATIONS=EW EW|LO1; 
+  EW_SCHEME=3
+  GF=1.166397e-5
+  METS_BBAR_MODE=5
+
+  # speed and neg weight fraction improvements
+  PP_RS_SCALE VAR{sqr(sqrt(H_T2)-PPerp(p[2])-PPerp(p[3])+MPerp(p[2]+p[3]))/4};
+  NLO_CSS_PSMODE=1
+
+}(run)
+
+(processes){
+  Process 93 93 -> 22 11 -12 93{NJET}
+  Enhance_Function VAR{pow(max(PPerp(p[2]+p[3]),PPerp(p[4]))/20.0,2)} 
+  Associated_Contributions EW|LO1 {LJET};
+  Order (*,3); CKKW sqr(QCUT/E_CMS)
+  NLO_QCD_Mode MC@NLO {LJET}
+  ME_Generator Amegic {LJET}
+  RS_ME_Generator Comix {LJET}
+  Loop_Generator LOOPGEN {LJET}
+  PSI_ItMin 20000 {4}
+  Integration_Error 0.99 {4}
+  PSI_ItMin 50000 {5,6,7}
+  Integration_Error 0.99 {5,6,7}
+  End process
+
+  Process 93 93 -> 22 -11 12 93{NJET}
+  Enhance_Function VAR{pow(max(PPerp(p[2]+p[3]),PPerp(p[4]))/20.0,2)}
+  Associated_Contributions EW|LO1 {LJET};
+  Order (*,3); CKKW sqr(QCUT/E_CMS)
+  NLO_QCD_Mode MC@NLO {LJET}
+  ME_Generator Amegic {LJET}
+  RS_ME_Generator Comix {LJET}
+  Loop_Generator LOOPGEN {LJET}
+  PSI_ItMin 20000 {4}
+  Integration_Error 0.99 {4}
+  PSI_ItMin 50000 {5,6,7}
+  Integration_Error 0.99 {5,6,7}
+  End process
+}(processes)
+
+(selector){
+  PTNLO  22  7  E_CMS
+  IsolationCut  22  0.1  2  0.10
+  DeltaRNLO  22  90  0.1 1000.0
+  Mass  90  91  1.0  E_CMS
+}(selector)
+"""
+
+genSeq.Sherpa_i.NCores = 128
+genSeq.Sherpa_i.Parameters += [ "OL_PARAMETERS=ew_scheme=2 ew_renorm_scheme=1 write_parameters=1" ]
+genSeq.Sherpa_i.Parameters += [ "EW_SCHEME=3", "GF=1.166397e-5" ]
diff --git a/700xxx/700403/mc_13TeV.Sh_2211_enugamma.GRID.tar.gz b/700xxx/700403/mc_13TeV.Sh_2211_enugamma.GRID.tar.gz
new file mode 120000
index 0000000000000000000000000000000000000000..ab9a9ffa1e1bcdb94829cdee9f888a1914af3c5d
--- /dev/null
+++ b/700xxx/700403/mc_13TeV.Sh_2211_enugamma.GRID.tar.gz
@@ -0,0 +1 @@
+/eos/user/m/mgignac/mc/2211/Vy_final/mc_13TeV.Sh_2211_enugamma.GRID.tar.gz
\ No newline at end of file
diff --git a/700xxx/700404/mc.Sh_2211_taunugamma.py b/700xxx/700404/mc.Sh_2211_taunugamma.py
new file mode 100644
index 0000000000000000000000000000000000000000..1116321b83860b667f4aee6d930dbcf1994573fe
--- /dev/null
+++ b/700xxx/700404/mc.Sh_2211_taunugamma.py
@@ -0,0 +1,82 @@
+
+include("Sherpa_i/Base_Fragment.py")
+include("Sherpa_i/NNPDF30NNLO.py")
+
+evgenConfig.description = "Sherpa taunugamma + 0,1j@NLO + 2,3,4j@LO with 7<pT_y"
+evgenConfig.keywords = ["SM", "tau", "neutrino", "photon", "NLO" ]
+evgenConfig.contact  = [ "atlas-generators-sherpa@cern.ch", "matthew.gignac@cern.ch" ]
+evgenConfig.nEventsPerJob = 1000
+
+genSeq.Sherpa_i.RunCard="""
+(run){
+
+  %scales, tags for scale variations
+  FSF:=1.; RSF:=1.; QSF:=1.;
+  SCALES STRICT_METS{FSF*MU_F2}{RSF*MU_R2}{QSF*MU_Q2};
+
+  # Shower improvements
+  NLO_SUBTRACTION_SCHEME=2;
+
+  # me generator settings
+  ME_SIGNAL_GENERATOR Comix Amegic LOOPGEN;
+  LOOPGEN:=OpenLoops
+
+  # tags for process setup
+  NJET:=4; LJET:=3,4; QCUT:=20;
+
+  # EW corrections setup
+  OL_PARAMETERS=ew_scheme 2 ew_renorm_scheme 1
+  ASSOCIATED_CONTRIBUTIONS_VARIATIONS=EW EW|LO1; 
+  EW_SCHEME=3
+  GF=1.166397e-5
+  METS_BBAR_MODE=5
+
+  SOFT_SPIN_CORRELATIONS=1;
+
+  # speed and neg weight fraction improvements
+  PP_RS_SCALE VAR{sqr(sqrt(H_T2)-PPerp(p[2])-PPerp(p[3])+MPerp(p[2]+p[3]))/4};
+  NLO_CSS_PSMODE=1
+
+}(run)
+
+(processes){
+  Process 93 93 -> 22 15 -16 93{NJET}
+  Enhance_Function VAR{pow(max(PPerp(p[2]+p[3]),PPerp(p[4]))/20.0,2)}
+  Associated_Contributions EW|LO1 {LJET};
+  Order (*,3); CKKW sqr(QCUT/E_CMS)
+  NLO_QCD_Mode MC@NLO {LJET}
+  ME_Generator Amegic {LJET}
+  RS_ME_Generator Comix {LJET}
+  Loop_Generator LOOPGEN {LJET}
+  PSI_ItMin 20000 {4}
+  Integration_Error 0.99 {4}
+  PSI_ItMin 50000 {5,6,7}
+  Integration_Error 0.99 {5,6,7}
+  End process
+
+  Process 93 93 -> 22 -15 16 93{NJET}
+  Enhance_Function VAR{pow(max(PPerp(p[2]+p[3]),PPerp(p[4]))/20.0,2)}
+  Associated_Contributions EW|LO1 {LJET};
+  Order (*,3); CKKW sqr(QCUT/E_CMS)
+  NLO_QCD_Mode MC@NLO {LJET}
+  ME_Generator Amegic {LJET}
+  RS_ME_Generator Comix {LJET}
+  Loop_Generator LOOPGEN {LJET}
+  PSI_ItMin 20000 {4}
+  Integration_Error 0.99 {4}
+  PSI_ItMin 50000 {5,6,7}
+  Integration_Error 0.99 {5,6,7}
+  End process
+}(processes)
+
+(selector){
+  PTNLO  22  7  E_CMS
+  IsolationCut  22  0.1  2  0.10
+  DeltaRNLO  22  90  0.1 1000.0
+  Mass  90  91  1.0  E_CMS
+}(selector)
+"""
+
+genSeq.Sherpa_i.NCores = 128
+genSeq.Sherpa_i.Parameters += [ "OL_PARAMETERS=ew_scheme=2 ew_renorm_scheme=1 write_parameters=1" ]
+genSeq.Sherpa_i.Parameters += [ "EW_SCHEME=3", "GF=1.166397e-5" ]
diff --git a/700xxx/700404/mc_13TeV.Sh_2211_taunugamma.GRID.tar.gz b/700xxx/700404/mc_13TeV.Sh_2211_taunugamma.GRID.tar.gz
new file mode 120000
index 0000000000000000000000000000000000000000..8917e1adfcec3161a018c1b0dc9608faca77cb93
--- /dev/null
+++ b/700xxx/700404/mc_13TeV.Sh_2211_taunugamma.GRID.tar.gz
@@ -0,0 +1 @@
+/eos/user/m/mgignac/mc/2211/Vy_final/mc_13TeV.Sh_2211_taunugamma.GRID.tar.gz
\ No newline at end of file