#### If true fSIDIS=+fDY, (wrong)
#### if false fSIDIS=-fDY (correct)
useWrongSign = False
#%%
#######################################
#Initialize artemide
#######################################
import harpy
path_to_constants = MAINPATH + "FittingPrograms/Sivers20/Constants-files/"
#harpy.initialize(path_to_constants+"const-Sivers20_lo")
# harpy.initialize(path_to_constants+"const-Sivers20_nnlo_piK")
if (useOrder == "nnlo"):
harpy.initialize(path_to_constants + "const-Sivers20_nnlo")
#### All=0 Case
harpy.setNPparameters_TMDR([2., 0.0398333])
harpy.setNPparameters_uTMDPDF([
0.185239, 6.22706, 580.946, 2.44166, -2.53161, 0., 0.0014, 0.442, 4.14
])
harpy.setNPparameters_uTMDFF([0.279443, 0.460015, 0.435955, 0.551302])
# # #### All=0 case piK
# harpy.setNPparameters_TMDR([2., 0.0394095])
# harpy.setNPparameters_uTMDPDF([0.180718, 4.38119, 426.208, 2.22347, -0.0646396, 0., 0.17, 0.48, 2.15])
# harpy.setNPparameters_uTMDFF([0.293548, 0.462093, 0.442867, 0.590596, 0.427915, 0.462578, 0.304421,1.18113])
elif (useOrder == "n3lo"):
harpy.initialize(path_to_constants + "const-Sivers20_n3lo")
#### All=0 Case n3lo
# importing libraries
#######################################
import numpy
import time
import DataProcessor.harpyInterface
import DataProcessor.DataMultiSet
import DataProcessor.ArtemideReplicaSet
#%%
#######################################
#Initialize artemide
#######################################
import harpy
harpy.initialize(PathToConstantsFile)
#%%
rSet=DataProcessor.ArtemideReplicaSet.ReadRepFile(\
"/home/vla18041/LinkData2/WorkingFiles/TMD/Fit_Notes/TMD<->PDF/REPS/SV21-"+PDFinUse+"-nnlo-"+CASE+".rep")
rSet.SetReplica()
#%%
#######################################
# read the list of files and return the list of DataSets
#######################################
def loadThisData(listOfNames):
import DataProcessor.DataSet
dataCollection = []
for name in listOfNames:
sys.path.append("/home/vla18041/LinkData2/arTeMiDe_Repository/DataProcessor")
import DataProcessor.harpyInterface
import DataProcessor.DataMultiSet
import DataProcessor.DataSet
import DataProcessor.ArtemideReplicaSet
MAINPATH="/home/vla18041/LinkData2/arTeMiDe_Repository/DataProcessor/OtherPrograms/EIC_impact_Study/"
#%%
rSet=DataProcessor.ArtemideReplicaSet.ReadRepFile("/home/vla18041/LinkData2/arTeMiDe_Repository/artemide/Models/SV19/Replicas/DY+SIDIS/SV19_nnlo_all=0.rep")
#%%
#######################################
#Initialize artemide
#######################################
import harpy
harpy.initialize(MAINPATH+"/Constants-files/"+"const-NNPDF31+DSS_nnlo_reduce")
rSet.SetReplica(0)
#%%
def cutFuncFUU(p):
from math import log10, floor
def round_2dg(x):
return round(x, 1-int(floor(log10(abs(x)))))
#p["<x>"]=numpy.round(numpy.mean(p["x"]),5)
#p["<z>"]=numpy.round(numpy.mean(p["z"]),3)
#p["<pT>"]=numpy.round(numpy.mean(p["pT"]),2)
p["<x>"]=round_2dg(numpy.mean(p["x"]))
p["<z>"]=round_2dg(numpy.mean(p["z"]))
p["<pT>"]=round_2dg(numpy.mean(p["pT"]))
import sys
import numpy
sys.path.append("/home/vla18041/LinkData2/arTeMiDe_Repository/DataProcessor")
import DataProcessor.harpyInterface
import DataProcessor.ArtemideReplicaSet
import DataProcessor.DataMultiSet
MAINPATH = "/home/vla18041/LinkData2/arTeMiDe_Repository/DataProcessor/"
#%%
#######################################
#Initialize artemide with a replica -2
#######################################
import harpy
path_to_constants = MAINPATH + "FittingPrograms/SV19/Constants-files/"
harpy.initialize(path_to_constants + "DY_nnlo/const-NNPDF31_NNLO")
harpy.setNPparameters([
2.0340, 0.0299, 0.2512, 7.7572, 334.6108, 2.4543, -4.8203, 0.1000, 0.0000
])
#%%
### read the list of files and return the list of DataSets
import DataProcessor.DataSet
data = DataProcessor.DataSet.LoadCSV(
"/home/vla18041/LinkData2/arTeMiDe_Repository/DataProcessor/OtherPrograms/Predictions/ATLAS8/A8-predict.csv"
)
############## I drop all points after ~0.25, last two bins are not trastful
import sys
import time
import numpy
sys.path.append("/home/vla18041/LinkData2/arTeMiDe_Repository/DataProcessor")
import DataProcessor.harpyInterface
import DataProcessor.ArtemideReplicaSet
MAINPATH="/home/vla18041/LinkData2/arTeMiDe_Repository/DataProcessor/"
#%%
#######################################
#Initialize artemide with a replica -2
#######################################
import harpy
path_to_constants=MAINPATH+"FittingPrograms/SV19/Constants-files/"
harpy.initialize(path_to_constants+"DY_n3lo/const-NNPDF31_n3lo")
#harpy.initialize(path_to_constants+"DY_nnlo/const-HERA20_NNLO")
#harpy.initialize(path_to_constants+"DY_n3lo/const-HERA20_n3lo")
#harpy.initialize(path_to_constants+"DY_nnlo/const-MMHT14_NNLO")
#harpy.initialize(path_to_constants+"DY_nnlo/const-CT14_NNLO")
#harpy.initialize(path_to_constants+"DY_nnlo/const-PDF4LHC_NNLO")
harpy.setNPparameters([2.0340, 0.0299, 0.2512, 7.7572,334.6108, 2.4543,-4.8203, 0.1000, 0.0000])
#harpy.setNPparameters_TMDR(-2)
#harpy.setNPparameters_uTMDPDF(-2)
#%%
rSet=DataProcessor.ArtemideReplicaSet.ReadRepFile("/home/vla18041/LinkData2/arTeMiDe_Repository/artemide/Models/SV19/Replicas/"+
"DY_n3lo/DY_NNPDF31_n3lo.rep")
# "Sivers20_model9case1(noDY-n3lo).rep")
sys.path.append("/home/vla18041/LinkData2/arTeMiDe_Repository/DataProcessor")
import DataProcessor.harpyInterface
import DataProcessor.DataMultiSet
import DataProcessor.ArtemideReplicaSet
MAINPATH = "/home/vla18041/LinkData2/arTeMiDe_Repository/DataProcessor/"
#%%
#######################################
#Initialize artemide
#######################################
import harpy
path_to_constants = MAINPATH + "FittingPrograms/SV19/Constants-files/" + "DY+SIDIS_nnlo/const-DY+SIDIS_NNPDF31+DSS_nnlo"
harpy.initialize(path_to_constants)
harpy.setNPparameters_TMDR([1.93, 0.0434])
harpy.setNPparameters_uTMDPDF(
[0.253434, 9.04351, 346.999, 2.47992, -5.69988, 0.1, 0.])
harpy.setNPparameters_uTMDFF([0.264, 0.479, 0.459, 0.539])
#%%
##################Cut function
def cutFunc(p):
if p["type"] == "DY":
delta = p["<qT>"] / p["<Q>"]
return delta < 0.25, p
with open(logFile, 'a') as file:
file.write(PCname + ' : ' + time.ctime() + ' : [+' + str(hours) +
':' + str(minutes) + ':' + str(seconds) + ' ]\n')
file.write(' --> ' + text + '\n')
file.write('\n')
savedTime = time.time()
#%%
#######################################
#Initialize artemide
#######################################
import harpy
path_to_constants = MAINPATH + "FittingPrograms/SV19/Constants-files/"
harpy.initialize(
path_to_constants +
"DY+SIDIS_nnlo_piK_all=0/const-DY+SIDIS_NNPDF31+DSS_nnlo_piK_all=0")
harpy.setNPparameters_TMDR([1.93, 0.0434])
harpy.setNPparameters_uTMDPDF(
[0.253434, 9.04351, 346.999, 2.47992, -5.69988, 0.1, 0.])
harpy.setNPparameters_uTMDFF(
[0.264, 0.479, 0.459, 0.539, 0.264, 0.479, 0.459, 0.539])
#%%
### read the list of files and return the list of DataSets
def loadThisDataDY(listOfNames):
import DataProcessor.DataSet
path_to_data = "/home/vla18041/LinkData2/arTeMiDe_Repository/DataProcessor/DataLib/unpolDY/"
import DataProcessor.ArtemideReplicaSet
#prefix="ff_K_n3lo_all=0"
prefix = "n3lo_all=0"
MAINPATH = "/home/vla18041/LinkData2/arTeMiDe_Repository/DataProcessor/OtherPrograms/TMDGrids/"
SAVEPATH = "/home/vla18041/LinkData2/arTeMiDe_Repository/DataProcessor/OtherPrograms/TMDGrids/Grids/SV19_" + prefix + "/SV19_" + prefix
#%%
#######################################
# Initializing artemide
#######################################
import harpy
path_to_constants = MAINPATH + "Constants/"
#harpy.initialize(path_to_constants+"const-TMDPDF-NNPDF")
harpy.initialize(path_to_constants + "const-TMDPDF-NNPDF-n3lo")
#harpy.initialize(path_to_constants+"const-TMDFF-DSS-n3lo")
#rSet=DataProcessor.ArtemideReplicaSet.ReadRepFile("/home/vla18041/LinkData2/arTeMiDe_Repository/artemide/Models/SV19/Replicas/DY+SIDIS/SV19_nnlo_all=0.rep")
#rSet=DataProcessor.ArtemideReplicaSet.ReadRepFile("/home/vla18041/LinkData2/arTeMiDe_Repository/artemide/Models/SV19/Replicas/DY+SIDIS/SV19_nnlo.rep")
#rSet=DataProcessor.ArtemideReplicaSet.ReadRepFile("/home/vla18041/LinkData2/arTeMiDe_Repository/artemide/Models/SV19/Replicas/DY+SIDIS/SV19_n3lo.rep")
rSet = DataProcessor.ArtemideReplicaSet.ReadRepFile(
"/home/vla18041/LinkData2/arTeMiDe_Repository/artemide/Models/SV19/Replicas/DY+SIDIS/SV19_n3lo_all=0.rep"
)
rSet.SetReplica(num=0, part="uTMDPDF")
#%%
#######################################
# Save Grid specification
#######################################
with open(SAVEPATH + '.info', 'w') as outfile:
outfile.write("SetDesc: Unpolarized TMDFF (kaon+) SV19_" + prefix + "\n")
#######################################
# Initializing artemide
#######################################
import harpy
path_to_constants=MAINPATH+"Constants/"
###
# prefix="n3lo_all=0"
# cases=["_","_ff_pi_","_ff_K_"]
# harpy.initialize(path_to_constants+"const-SV19-n3lo")
# rSet=DataProcessor.ArtemideReplicaSet.ReadRepFile("/home/vla18041/LinkData2/arTeMiDe_Repository/artemide/Models/SV19/Replicas/DY+SIDIS/SV19_n3lo_all=0_hesse.rep")
###
prefix="n3lo"
cases=["_","_ff_pi_","_ff_K_"]
harpy.initialize(path_to_constants+"const-SV19-n3lo")
rSet=DataProcessor.ArtemideReplicaSet.ReadRepFile("/home/vla18041/LinkData2/arTeMiDe_Repository/artemide/Models/SV19/Replicas/DY+SIDIS/SV19_n3lo_hesse.rep")
###
# prefix="nnlo_all=0"
# cases=["_"]
# harpy.initialize(path_to_constants+"const-SV19-nnlo")
# rSet=DataProcessor.ArtemideReplicaSet.ReadRepFile("/home/vla18041/LinkData2/arTeMiDe_Repository/artemide/Models/SV19/Replicas/DY+SIDIS/SV19_nnlo_all=0_hesse.rep")
##
# prefix="nnlo"
# cases=["_"]
# harpy.initialize(path_to_constants+"const-SV19-nnlo")
# rSet=DataProcessor.ArtemideReplicaSet.ReadRepFile("/home/vla18041/LinkData2/arTeMiDe_Repository/artemide/Models/SV19/Replicas/DY+SIDIS/SV19_nnlo_hesse.rep")
rSet.SetReplica(0)
import DataProcessor.ArtemideReplicaSet
prefix="pion"
#prefix="pion_all=0"
MAINPATH="/home/vla18041/LinkData2/arTeMiDe_Repository/DataProcessor/OtherPrograms/TMDGrids/"
SAVEPATH="/home/vla18041/LinkData2/arTeMiDe_Repository/DataProcessor/OtherPrograms/TMDGrids/Grids/SV19_"+prefix+"/SV19_"+prefix
#%%
#######################################
# Initializing artemide
#######################################
import harpy
path_to_constants=MAINPATH+"Constants/"
#harpy.initialize(path_to_constants+"const-TMDPDF-NNPDF")
harpy.initialize(path_to_constants+"const-Vpion19_nnlo")
#rSet=DataProcessor.ArtemideReplicaSet.ReadRepFile("/home/vla18041/LinkData2/arTeMiDe_Repository/artemide/Models/SV19/Replicas/DY+SIDIS/SV19_n3lo_all=0.rep")
rSet=DataProcessor.ArtemideReplicaSet.ReadRepFile("/home/vla18041/LinkData2/arTeMiDe_Repository/artemide/Models/SV19/Replicas/DY+SIDIS/SV19_n3lo.rep")
rSet.SetReplica(num=0,part="TMDR")
harpy.setNPparameters_uTMDPDF([0.17975, 3.9081, 453.883, 2.07199, 1.60774, 0, 0.0014, 0.442, 4.14])
#%%
#######################################
# Save Grid specification
#######################################
with open(SAVEPATH+'.info', 'w') as outfile:
outfile.write("SetDesc: Unpolarized TMDPDF (pion+) SV19_"+prefix+". The sets n=1,2,...,6 are"
+" boundaries of Hessian uncertanty band coresponding to 3 NP parameters of the model. \n")
outfile.write("Authors: A.Vladimirov"+"\n"),
import sys
import time
import numpy
sys.path.append("/home/vla18041/LinkData2/arTeMiDe_Repository/DataProcessor")
import DataProcessor.harpyInterface
import DataProcessor.DataMultiSet
MAINPATH = "/home/vla18041/LinkData2/arTeMiDe_Repository/DataProcessor/"
#%%
#######################################
#Initialize artemide with a replica -2
#######################################
import harpy
path_to_constants = MAINPATH + "FittingPrograms/PDF-TMD/Constants-files/SV19/"
harpy.initialize(path_to_constants + "SV19-NNPDF31_NNLO")
#harpy.initialize(path_to_constants+"SV19-CT18_NNLO")
#harpy.initialize(path_to_constants+"SV19-MSHT20_NNLO")
harpy.setNPparameters([
2.0340, 0.0299, 0.2512, 7.7572, 334.6108, 2.4543, -4.8203, 0.1000, 0.0000
])
#harpy.setNPparameters_TMDR(-2)
#harpy.setNPparameters_uTMDPDF(-2)
#%%
### read the list of files and return the list of DataSets
def loadThisData(listOfNames):
import DataProcessor.DataSet
