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