Fix.

Author: npierre

HadronCount.jl

@@ -0,0 +1,65 @@
+using DelimitedFiles
+
+HadronN = readdlm("data/multiplicities_kaon_norich.txt")
+HadronM = readdlm("data/Kaon_RC_VM_20190327_noRICHunfold.txt")
+
+PionN = readdlm("data/multiplicities_kaon_norich.txt")
+PionM = readdlm("data/Kaon_RC_VM_20190327_noRICHunfold.txt")
+
+KaonN = readdlm("data/multiplicities_kaon_norich.txt")
+KaonM = readdlm("data/Kaon_RC_VM_20190327_noRICHunfold.txt")
+
+ProtonN = readdlm("data/multiplicities_kaon_norich.txt")
+ProtonM = readdlm("data/Kaon_RC_VM_20190327_noRICHunfold.txt")
+
+CountH = zeros(540,9)
+CountP = zeros(540,9)
+CountK = zeros(540,9)
+CountPr = zeros(540,9)
+
+for i in 1:540
+    CountH[i,1] = HadronN[i,1]
+    CountH[i,2] = HadronN[i,2]
+    CountH[i,3] = HadronN[i,3]
+    CountH[i,4] = HadronN[i,11]-HadronM[i,10]
+    CountH[i,5] = HadronN[i,12]-HadronM[i,11]
+    CountH[i,6] = HadronN[i,21]-HadronM[i,12]
+    CountH[i,7] = HadronN[i,22]-HadronM[i,13]
+    CountH[i,8] = HadronN[i,9]>0 ? (HadronN[i,8]-HadronM[i,4])/HadronN[i,9] : 0
+    CountH[i,9] = HadronN[i,18]>0 ? (HadronN[i,18]-HadronM[i,6])/HadronN[i,19] : 0
+
+    CountP[i,1] = PionN[i,1]
+    CountP[i,2] = PionN[i,2]
+    CountP[i,3] = PionN[i,3]
+    CountP[i,4] = PionN[i,11]-PionM[i,10]
+    CountP[i,5] = PionN[i,12]-PionM[i,11]
+    CountP[i,6] = PionN[i,21]-PionM[i,12]
+    CountP[i,7] = PionN[i,22]-PionM[i,13]
+    CountP[i,8] = PionN[i,9]>0 ? (PionN[i,8]-PionM[i,4])/PionN[i,9] : 0
+    CountP[i,9] = PionN[i,18]>0 ? (PionN[i,18]-PionM[i,6])/PionN[i,19] : 0
+
+    CountK[i,1] = KaonN[i,1]
+    CountK[i,2] = KaonN[i,2]
+    CountK[i,3] = KaonN[i,3]
+    CountK[i,4] = KaonN[i,11]-KaonM[i,10]
+    CountK[i,5] = KaonN[i,12]-KaonM[i,11]
+    CountK[i,6] = KaonN[i,21]-KaonM[i,12]
+    CountK[i,7] = KaonN[i,22]-KaonM[i,13]
+    CountK[i,8] = KaonN[i,9]>0 ? (KaonN[i,8]-KaonM[i,4])/KaonN[i,9] : 0
+    CountK[i,9] = KaonN[i,18]>0 ? (KaonN[i,18]-KaonM[i,6])/KaonN[i,19] : 0
+
+    CountPr[i,1] = ProtonN[i,1]
+    CountPr[i,2] = ProtonN[i,2]
+    CountPr[i,3] = ProtonN[i,3]
+    CountPr[i,4] = ProtonN[i,11]-ProtonM[i,10]
+    CountPr[i,5] = ProtonN[i,12]-ProtonM[i,11]
+    CountPr[i,6] = ProtonN[i,21]-ProtonM[i,12]
+    CountPr[i,7] = ProtonN[i,22]-ProtonM[i,13]
+    CountPr[i,8] = ProtonN[i,9]>0 ? (ProtonN[i,8]-ProtonM[i,4])/ProtonN[i,9] : 0
+    CountPr[i,9] = ProtonN[i,18]>0 ? (ProtonN[i,18]-ProtonM[i,6])/ProtonN[i,19] : 0
+end
+
+writedlm("plots/CountHadron.txt", CountH)
+writedlm("plots/CountPion.txt", CountP)
+writedlm("plots/CountKaon.txt", CountK)
+writedlm("plots/CountProton.txt", CountPr)

KaonCount.jl

@@ -0,0 +1,14 @@
+using DelimitedFiles
+
+HadronN = readdlm("test.txt")
+HadronM = readdlm("out1.txt")
+list = Int64[]
+
+for i in 1:4280
+    if (abs(HadronN[i,1]-HadronM[i,1])> 0.00001 || abs(HadronN[i,2]-HadronM[i,2]) > 0.01 || abs(HadronN[i,3]-HadronM[i,3]) > 0.0001
+        || abs(HadronN[i,4]-HadronM[i,4]) > 0.01 || abs(HadronN[i,5]-HadronM[i,5]) > 0.0001)
+        push!(list,i)
+    end
+end
+
+writedlm("plots/QuickSort.txt", list)