def repchi_2(x): startT=time.time() nn=[0.184739, 6.22437, 588.193, 2.44327, -2.51106, 0., x[0],x[1],x[2]] harpy.setNPparameters_uTMDPDF(nn) print('np set =',["{:8.3f}".format(i) for i in x], end =" ") ccDY2,cc3=DataProcessor.harpyInterface.ComputeChi2(repDataDY) cc=ccDY2/totalNnew endT=time.time() print(':->',cc,' t=',endT-startT) return ccDY2
####################################### #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 harpy.setNPparameters_TMDR([2., 0.0442327]) harpy.setNPparameters_uTMDPDF( [0.17975, 3.9081, 453.883, 2.07199, 1.60774, 0, 0.0014, 0.442, 4.14]) harpy.setNPparameters_uTMDFF([0.26994, 0.456091, 0.423312, 0.615092])
def SetUnReplica(n): unSet.SetReplica(n, part="TMDR") unSet.SetReplica(n, part="uTMDFF") rr = unSet.GetReplica(n, part="uTMDPDF") harpy.setNPparameters_uTMDPDF(rr[:-1] + [0.0014, 0.442, 4.14])
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 unSet = DataProcessor.ArtemideReplicaSet.ReadRepFile( "/home/vla18041/LinkData2/arTeMiDe_Repository/artemide/Models/SV19/Replicas/DY+SIDIS/SV19_nnlo_all=0.rep" ) unSet.SetReplica(0, part="TMDR") unSet.SetReplica(0, part="uTMDFF") rr = unSet.GetReplica(0, part="uTMDPDF") harpy.setNPparameters_uTMDPDF(rr[:-1] + [0.0014, 0.442, 4.14]) # # #### 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 unSet = DataProcessor.ArtemideReplicaSet.ReadRepFile( "/home/vla18041/LinkData2/arTeMiDe_Repository/artemide/Models/SV19/Replicas/DY+SIDIS/SV19_n3lo_all=0.rep" ) unSet.SetReplica(0, part="TMDR") unSet.SetReplica(0, part="uTMDFF") rr = unSet.GetReplica(0, part="uTMDPDF")
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 #%% ### Loading the DY data set theData = DataProcessor.DataSet.LoadCSV( "/home/vla18041/LinkData2/arTeMiDe_Repository/DataProcessor/OtherPrograms/LHC_impact_Study/LHC13TeV-Empty.csv"
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"), outfile.write("Reference: arXiv:1907.10356"+"\n"), outfile.write("SetIndex: 710000"+"\n"), outfile.write("TMDScheme: Pavia TMDs"+"\n"), outfile.write("TMDType: pdf"+"\n"), outfile.write("CollDist: JAM18PionPDFnlo"+"\n"),
import DataProcessor.harpyInterface import DataProcessor.DataMultiSet MAINPATH="/home/vla18041/LinkData2/arTeMiDe_Repository/DataProcessor/FittingPrograms/Vpion19/" #%% ####################################### #Initialize artemide with a replica -2 ####################################### import harpy path_to_constants=MAINPATH+"Constants-files/" #harpy.initialize(path_to_constants+"DY_nnlo/const-Vpion19_nnlo") harpy.initialize(path_to_constants+"const-Vpion19_nnlo_all=0") #### All=0 Case harpy.setNPparameters_TMDR([2., 0.0398333]) harpy.setNPparameters_uTMDPDF([0.184739, 6.22437, 588.193, 2.44327, -2.51106, 0., 0.17, 0.48, 2.15]) #%% ### read the list of files and return the list of DataSets def loadThisData(listOfNames): import DataProcessor.DataSet path_to_data="/home/vla18041/LinkData2/arTeMiDe_Repository/DataProcessor/DataLib/piDY/" dataCollection=[] for name in listOfNames: loadedData=DataProcessor.DataSet.LoadCSV(path_to_data+name+".csv") dataCollection.append(loadedData)
def SetReplica(self, num=0, part='full'): """ Set the replica according to the list Parameters ---------- num : -1 = initla replica 0 = mean replica 1.... = replica from list The default is 0. part : string, optional Specification which part of the replica to be set. The default is "full". Possible values: 'full', 'TMDR', 'uTMDPDF', 'uTMDFF', 'lpTMDPDF', 'SiversTMDPDF', etc Returns ------- None. """ import harpy if (num == 0): r = self.meanReplica elif (num == -1): r = self.initialReplica else: r = self.replicaList[num - 1] ## resolvinf the part-condition if (part == "full"): doTMDR = True douTMDPDF = True douTMDFF = True dolpTMDPDF = True doSiversTMDPDF = True else: doTMDR = False douTMDPDF = False douTMDFF = False dolpTMDPDF = False doSiversTMDPDF = False if (part == 'TMDR'): doTMDR = True if (part == 'uTMDPDF'): douTMDPDF = True if (part == 'uTMDFF'): douTMDFF = True if (part == 'lpTMDPDF'): dolpTMDPDF = True if (part == 'SiversTMDPDF'): doSiversTMDPDF = True if (doTMDR and self._TMDRend >= self._TMDRstart + 1 > 0): harpy.setNPparameters_TMDR(r[self._TMDRstart:self._TMDRend]) if (douTMDPDF and self._uTMDPDFend >= self._uTMDPDFstart + 1 > 0): if (self._c_uTMDPDFend >= self._c_uTMDPDFstart + 1 > 0): harpy.setPDFreplica(r[self._c_uTMDPDFstart]) harpy.setNPparameters_uTMDPDF( r[self._uTMDPDFstart:self._uTMDPDFend]) if (douTMDFF and self._uTMDFFend >= self._uTMDFFstart + 1 > 0): if (self._c_uTMDFFend >= self._c_uTMDFFstart + 1 > 0): print("Modification of FF replica is not implimented") harpy.setNPparameters_uTMDFF(r[self._uTMDFFstart:self._uTMDFFend]) if (dolpTMDPDF and self._lpTMDPDFend >= self._lpTMDPDFstart + 1 > 0): if (self._c_lpTMDPDFend >= self._c_lpTMDPDFstart + 1 > 0): print("Modification of lpPDF replica is not implimented") harpy.setNPparameters_lpTMDPDF( r[self._lpTMDPDFstart:self._lpTMDPDFend]) if (doSiversTMDPDF and self._SiversTMDPDFend >= self._SiversTMDPDFstart + 1 > 0): if (self._SiversTMDPDFend >= self._SiversTMDPDFstart + 1 > 0): print("Modification of Sivers replica is not implimented") harpy.setNPparameters_SiversTMDPDF( r[self._SiversTMDPDFstart:self._SiversTMDPDFend])