diff --git a/Projects/s455/macro/yields/MapAsymmetry.C b/Projects/s455/macro/yields/MapAsymmetry.C
new file mode 100644
index 0000000000000000000000000000000000000000..cafb495b6101b9e0696de331bf2aad6554d05944
--- /dev/null
+++ b/Projects/s455/macro/yields/MapAsymmetry.C
@@ -0,0 +1,452 @@
+double GetAsymmetry(string, int Zb);
+double GetRMS(string, int Zb);
+double GetAsymmetryMoller(string, int Zb);
+double GetAsymmetrySmooth(string, int Zb);
+
+
+TH2I* h1;
+TH2I* h2;
+TH2I* h3;
+TGraph* gStableNuclei;
+TH2I* hStableNuclei;
+
+vector<int> vZ;
+vector<int> vA;
+vector<int> vN;
+vector<string> yield_file;
+
+int NumberOfZ=0;
+int NumberOfN=0;
+
+////////////////////////////////////////////
+void MapAsymmetry()
+{
+ gROOT->SetStyle("pierre_style");
+
+ gStyle->SetPalette(kInvertedDarkBodyRadiator);
+
+ NPL::Particle* npparticle = new NPL::Particle(1,1);
+ gStableNuclei = npparticle->GetStableNuclei();
+ gStableNuclei->SetMarkerStyle(21);
+ gStableNuclei->SetFillColor(1);
+
+ hStableNuclei = new TH2I("hstable","hstable",150,0,150,100,0,100);
+ int Npoint = gStableNuclei->GetN();
+ double gN, gZ;
+ for(int i=0; i<Npoint; i++){
+ gStableNuclei->GetPoint(i,gN,gZ);
+ hStableNuclei->Fill((int)gN, (int)gZ);
+ }
+ hStableNuclei->SetMarkerStyle(1);
+ hStableNuclei->SetFillColor(1);
+
+ h2 = new TH2I("h2","h2",150,0,150,95,0,95);
+ h2->GetXaxis()->SetTitle("Neutron number N");
+ h2->GetYaxis()->SetTitle("Proton number Z");
+ h2->GetXaxis()->CenterTitle();
+ h2->GetYaxis()->CenterTitle();
+ h2->SetMarkerStyle(1);
+ h2->SetMarkerColor(0);
+
+ h3 = new TH2I("h3","h3",150,0,150,100,0,100);
+ h3->GetXaxis()->SetTitle("Neutrons, N");
+ h3->GetYaxis()->SetTitle("Protons, Z");
+ h3->GetZaxis()->SetTitle("Asymmetry");
+ h3->GetXaxis()->CenterTitle();
+ h3->GetYaxis()->CenterTitle();
+ h3->GetZaxis()->CenterTitle();
+
+ h1 = new TH2I("h1","h1",150,0,150,100,0,100);
+ h1->GetXaxis()->SetTitle("N");
+ h1->GetYaxis()->SetTitle("Z");
+ h1->GetXaxis()->CenterTitle();
+ h1->GetYaxis()->CenterTitle();
+
+
+ string filename = "list_all.dat";
+ //string filename = "list_s455.dat";
+ ifstream ifile;
+ ifile.open(filename.c_str());
+
+ ofstream ofile;
+ string output_filename = "asymmetry_file.dat";
+ ofile.open(output_filename.c_str());
+
+ int Z=0;
+ int Zprev=0;
+ int A;
+ int N;
+ NPL::Particle* iso;
+ while(ifile>>Z>>A){
+ N = A-Z;
+ //ifile >> Z >> A;
+ iso = new NPL::Particle(Z,A);
+ string isoname = iso->GetName();
+ if(Z==92)
+ isoname.resize(4);
+ else
+ isoname.resize(5);
+ string yield_filename = "dat/yield_" + isoname + ".dat";
+ yield_file.push_back(yield_filename);
+ //cout << yield_filename << endl;
+ vZ.push_back(Z);
+ vA.push_back(A);
+ vN.push_back(A-Z);
+
+ double asy = GetAsymmetry(yield_filename, Z)*100;
+ double asy_moller = GetAsymmetryMoller(yield_filename, Z)*100;
+ double asy_smooth = GetAsymmetrySmooth(yield_filename, Z);
+ //double asy = GetRMS(yield_filename, Z)/Z*100;
+
+ ofile << Z << " " << A-Z << " " << 102 - asy_smooth*100 << endl;
+
+ h1->Fill(N,Z,asy);
+ //h2->Fill(N,Z,asy_moller);
+ h2->Fill(N,Z,102-asy_smooth*100);
+ //h3->Fill(N,Z,110-GetAsymmetryMoller(yield_filename, Z)*100);
+ h3->Fill(N,Z,101-GetAsymmetrySmooth(yield_filename, Z)*100);
+ }
+
+ TExec* ex1 = new TExec("ex1","gStyle->SetPalette(kBird)");
+ TExec* ex2 = new TExec("ex2","gStyle->SetPalette(kDarkBodyRadiator)");
+
+
+ h2->GetXaxis()->SetRangeUser(90,150);
+ h2->GetYaxis()->SetRangeUser(70,95);
+ TCanvas* c1 = new TCanvas("c1","c1",1000,800);
+ c1->cd(1);
+ //hStableNuclei->Draw("a");
+ h2->Draw("colz");
+ ex2->Draw();
+ hStableNuclei->Draw("colzsame");
+ ex1->Draw();
+ //h2->Draw("surf2");
+ h2->Draw("colzsame");
+ //hStableNuclei->Draw("colsame");
+
+ TGraph* gX[200];
+ TGraph* gY[200];
+ for(int i=0; i<200; i++){
+ gX[i] = new TGraph();
+ gX[i]->SetPoint(0,i,0);
+ gX[i]->SetPoint(1,i,150);
+
+ gY[i] = new TGraph();
+ gY[i]->SetPoint(0,0,i);
+ gY[i]->SetPoint(1,150,i);
+
+ gX[i]->SetLineStyle(1);
+ gX[i]->SetLineWidth(1);
+ gY[i]->SetLineStyle(1);
+ gY[i]->SetLineWidth(1);
+
+ gX[i]->Draw("lsame");
+ gY[i]->Draw("lsame");
+ }
+
+ int magic[6] = {4,20,28,50,82,126};
+ TGraph* gXmagic_down[6];
+ TGraph* gXmagic_up[6];
+ TGraph* gYmagic_down[6];
+ TGraph* gYmagic_up[6];
+ for(int i=0; i<6; i++){
+ gXmagic_down[i] = new TGraph();
+ gXmagic_down[i]->SetPoint(0,magic[i],0);
+ gXmagic_down[i]->SetPoint(1,magic[i],150);
+
+ gXmagic_up[i] = new TGraph();
+ gXmagic_up[i]->SetPoint(0,magic[i]+1,0);
+ gXmagic_up[i]->SetPoint(1,magic[i]+1,150);
+
+ gYmagic_down[i] = new TGraph();
+ gYmagic_down[i]->SetPoint(0,0,magic[i]);
+ gYmagic_down[i]->SetPoint(1,150,magic[i]);
+
+ gYmagic_up[i] = new TGraph();
+ gYmagic_up[i]->SetPoint(0,0,magic[i]+1);
+ gYmagic_up[i]->SetPoint(1,150,magic[i]+1);
+
+ gXmagic_down[i]->SetLineStyle(1);
+ gXmagic_down[i]->SetLineWidth(3);
+ gXmagic_down[i]->SetLineColor(2);
+ gYmagic_down[i]->SetLineStyle(1);
+ gYmagic_down[i]->SetLineWidth(3);
+ gYmagic_down[i]->SetLineColor(2);
+
+ gXmagic_up[i]->SetLineStyle(1);
+ gXmagic_up[i]->SetLineWidth(3);
+ gXmagic_up[i]->SetLineColor(2);
+ gYmagic_up[i]->SetLineStyle(1);
+ gYmagic_up[i]->SetLineWidth(3);
+ gYmagic_up[i]->SetLineColor(2);
+
+ gXmagic_down[i]->Draw("lsame");
+ gYmagic_down[i]->Draw("lsame");
+
+ gXmagic_up[i]->Draw("lsame");
+ gYmagic_up[i]->Draw("lsame");
+ }
+
+
+ h3->GetXaxis()->SetTitle("Neutrons, N");
+ h3->GetYaxis()->SetTitle("Protons, Z");
+ h3->GetZaxis()->SetTitle("Asymmetry");
+ h3->GetXaxis()->CenterTitle();
+ h3->GetYaxis()->CenterTitle();
+ h3->GetZaxis()->CenterTitle();
+ h3->GetXaxis()->SetTitleOffset(1.9);
+ h3->GetYaxis()->SetTitleOffset(1.9);
+
+ hStableNuclei->GetXaxis()->SetRangeUser(90,150);
+ hStableNuclei->GetYaxis()->SetRangeUser(70,100);
+
+ h3->SetFillColor(1);
+ h3->SetLineColor(1);
+ hStableNuclei->SetLineColor(0);
+ //h3->Add(hStableNuclei);
+
+ h3->GetXaxis()->SetRangeUser(90,150);
+ h3->GetYaxis()->SetRangeUser(70,100);
+
+ // X axis //
+ h3->GetXaxis()->SetAxisColor(0);
+ h3->GetXaxis()->SetTickSize(0);
+ h3->GetXaxis()->SetLabelOffset(0.03);
+ h3->GetXaxis()->SetLabelSize(0);
+ h3->GetXaxis()->SetTitle("");
+
+ // Y axis //
+ h3->GetYaxis()->SetAxisColor(0);
+ h3->GetYaxis()->SetTickSize(0);
+ h3->GetYaxis()->SetLabelOffset(0.03);
+ h3->GetYaxis()->SetLabelSize(0);
+ h3->GetYaxis()->SetTitle("");
+
+ // Z axis //
+ h3->GetZaxis()->SetAxisColor(0);
+ h3->GetZaxis()->SetTickSize(0);
+ h3->GetZaxis()->SetLabelColor(0);
+ h3->GetZaxis()->SetTitle("");
+
+ TCanvas* c2 = new TCanvas("c2","c2",1600,1000);
+ c2->cd();
+
+ h3->Draw("0lego2fbbb");
+ double x[2] = {85,130};
+ double y[2] = {80.5,80.5};
+ double z[2] = {0,0};
+
+ TPolyLine3D* line = new TPolyLine3D(2,x,y,z);
+ line->SetLineStyle(2);
+ line->Draw();
+
+ x[0] = 110;
+ x[1] = 145;
+ y[0] = 90.5;
+ y[1] = 90.5;
+ line = new TPolyLine3D(2,x,y,z);
+ line->SetLineStyle(2);
+ line->Draw();
+
+ x[0] = 100.5;
+ x[1] = 100.5;
+ y[0] = 69;
+ y[1] = 87;
+ line = new TPolyLine3D(2,x,y,z);
+ line->SetLineStyle(2);
+ line->Draw();
+
+ x[0] = 120.5;
+ x[1] = 120.5;
+ y[0] = 69;
+ y[1] = 92;
+ line = new TPolyLine3D(2,x,y,z);
+ line->SetLineStyle(2);
+ line->Draw();
+
+ x[0] = 140.5;
+ x[1] = 140.5;
+ y[0] = 69;
+ y[1] = 95;
+ line = new TPolyLine3D(2,x,y,z);
+ line->SetLineStyle(2);
+ line->Draw();
+
+ h3->Draw("0lego2fbbbsame");
+
+ TExec* ex3 = new TExec("ex3","gStyle->SetPalette(kAvocado)");
+ ex3->Draw();
+
+ hStableNuclei->Draw("0lego2fbbbsame");
+
+ ex1->Draw();
+}
+
+////////////////////////////////////////////
+double GetAsymmetry(string filename, int Zb)
+{
+ double asymmetry;
+ ifstream ifile;
+
+ cout << "*** opening " << filename << endl;
+ ifile.open(filename.c_str());
+
+ vector<int> v_Z;
+ vector<double> v_yield;
+ vector<double> v_yield_err;
+
+ int Z;
+ double yield;
+ double yield_err;
+ while(!ifile.eof()){
+ ifile >> Z >> yield >> yield_err;
+ //cout << Z << " " << yield << " " << yield_err << endl;
+ if(yield>0 && yield_err>0){
+ v_Z.push_back(Z);
+ v_yield.push_back(yield);
+ v_yield_err.push_back(yield_err);
+ }
+ }
+ ifile.close();
+
+ unsigned int size = v_Z.size();
+ double norm=0;
+ double denum=0;
+ for(unsigned int i=0; i<size; i++){
+ norm += pow(v_yield[i],2)*pow(v_Z[i]-(double)Zb/2,2);
+ denum += pow(v_yield[i],2);
+ }
+
+ asymmetry = sqrt(norm/denum);
+
+ return asymmetry;
+}
+
+////////////////////////////////////////////
+double GetAsymmetryMoller(string filename, int Zb)
+{
+ double asymmetry;
+ ifstream ifile;
+
+ ifile.open(filename.c_str());
+
+ vector<int> v_Z;
+ vector<double> v_yield;
+ vector<double> v_yield_err;
+ double yield_max=0;
+ double yield_sym=0;
+ int Zmax;
+
+ int Z;
+ double yield;
+ double yield_err;
+ while(!ifile.eof()){
+ ifile >> Z >> yield >> yield_err;
+ //cout << Z << " " << yield << " " << yield_err << endl;
+ if(yield>0 && yield_err>0){
+ v_Z.push_back(Z);
+ v_yield.push_back(yield);
+ v_yield_err.push_back(yield_err);
+ if(yield>yield_max){
+ yield_max = yield;
+ Zmax = Z;
+ }
+ if(Z==int(Zb/2.)){
+ yield_sym = yield;
+ }
+ }
+ }
+ ifile.close();
+
+ asymmetry = yield_sym/yield_max;
+
+ return asymmetry;
+}
+
+////////////////////////////////////////////
+double GetAsymmetrySmooth(string filename, int Zb)
+{
+ double asymmetry;
+ ifstream ifile;
+
+ ifile.open(filename.c_str());
+
+ vector<int> v_Z;
+ vector<double> v_yield;
+ vector<double> v_yield_err;
+ double yield_max=0;
+ double yield_sym=0;
+ int Zmax;
+
+ int Z;
+ double yield;
+ double yield_err;
+ int i=0;
+ int i_max=0;
+ int i_sym=0;
+ while(!ifile.eof()){
+ ifile >> Z >> yield >> yield_err;
+ //cout << Z << " " << yield << " " << yield_err << endl;
+ if(yield>0 && yield_err>0){
+ v_Z.push_back(Z);
+ v_yield.push_back(yield);
+ v_yield_err.push_back(yield_err);
+ if(yield>yield_max){
+ yield_max = yield;
+ Zmax = Z;
+ i_max = i;
+ }
+ if(Z==round(Zb/2)){
+ yield_sym = yield;
+ i_sym = i;
+ }
+ i++;
+ }
+ }
+ ifile.close();
+
+ yield_sym = v_yield[i_sym] + 0.5*(v_yield[i_sym-1] + v_yield[i_sym+1]);
+ yield_max = v_yield[i_max] + 0.5*(v_yield[i_max-1] + v_yield[i_max+1]);
+
+
+ asymmetry = yield_sym/yield_max;
+ if(asymmetry>1)
+ asymmetry=1.;
+
+ cout << Zb << " " << asymmetry << endl;
+ return asymmetry;
+}
+
+
+
+////////////////////////////////////////////
+double GetRMS(string filename, int Zb)
+{
+ double rms;
+ ifstream ifile;
+
+ ifile.open(filename.c_str());
+
+ TGraphErrors* gerr = new TGraphErrors();
+
+ int Z;
+ double yield;
+ double yield_err;
+ int i=0;
+ while(!ifile.eof()){
+ ifile >> Z >> yield >> yield_err;
+ if(yield>0 && yield_err>0){
+ gerr->SetPoint(i,Z,yield);
+ gerr->SetPointError(i,0,yield_err);
+ i++;
+ }
+ }
+ ifile.close();
+
+ rms = gerr->GetRMS();
+
+
+ return rms;
+}
+
+
diff --git a/Projects/s455/macro/yields/asymmetry_file.dat b/Projects/s455/macro/yields/asymmetry_file.dat
new file mode 100644
index 0000000000000000000000000000000000000000..395fafcc87363582d5d0e8d92c58fd0e576cb04c
--- /dev/null
+++ b/Projects/s455/macro/yields/asymmetry_file.dat
@@ -0,0 +1,117 @@
+77 95 27.5901
+77 96 2
+77 97 2
+77 98 8.05876
+78 96 39.4121
+78 97 12.7227
+78 98 14.6431
+78 99 29.3682
+78 100 5.41142
+79 97 33.1098
+79 98 35.7432
+79 99 32.7443
+79 100 23.3638
+79 101 23.0197
+79 102 2
+79 103 9.1768
+80 99 76.6616
+80 100 56.2088
+80 101 50.0469
+80 102 44.4312
+80 103 33.6179
+80 104 28.5165
+80 105 12.1036
+80 106 31.8884
+81 102 63.6252
+81 103 37.1221
+81 104 30.8692
+81 105 24.9899
+81 106 12.3406
+81 107 6.35808
+82 104 25.634
+82 105 15.954
+82 106 17.4312
+82 107 10.1358
+82 108 2
+82 109 2
+82 110 2
+82 111 2
+83 106 2
+83 107 2
+83 108 2
+83 109 2
+83 110 2
+83 111 2
+83 112 2
+83 113 2
+83 114 2
+84 108 7.16565
+84 109 4.29144
+84 110 2
+84 111 2
+84 112 2
+84 113 2
+84 114 2
+84 115 4.63402
+84 116 2
+84 117 2
+85 111 22.2306
+85 112 2
+85 113 2
+85 114 2
+85 115 2
+85 116 2
+85 117 2.68491
+85 118 2
+86 114 2
+86 115 2
+86 116 2
+86 117 2
+86 118 2
+86 119 2.01963
+86 120 7.96152
+87 116 14.4873
+87 117 2
+87 118 2
+87 119 2
+87 120 2
+87 121 2
+87 122 2
+87 123 2
+87 124 2
+88 121 2
+88 122 2
+88 123 2
+88 124 2
+88 125 2
+88 126 2
+89 123 2
+89 124 2
+89 125 2
+89 126 2
+89 127 2
+89 128 2
+89 129 3.81067
+90 126 2
+90 127 2
+90 129 2
+90 130 2
+89 130 2
+89 131 3.30119
+89 134 2
+89 138 2
+90 131 2
+90 132 2
+90 133 2
+90 135 2
+90 136 7.7225
+90 139 72.2939
+90 140 81.2252
+91 137 57.9722
+91 138 71.6338
+91 140 82.1061
+91 141 93.1166
+92 142 90.428
+92 143 88.9814
+93 144 88.8437
+93 145 87.6857
diff --git a/Projects/s455/macro/yields/list_Hg.dat b/Projects/s455/macro/yields/list_Hg.dat
new file mode 100644
index 0000000000000000000000000000000000000000..9315d5b33680ac3e17d799d169ceb9cbf28a3c4b
--- /dev/null
+++ b/Projects/s455/macro/yields/list_Hg.dat
@@ -0,0 +1,2 @@
+80 180
+80 182
diff --git a/Projects/s455/macro/yields/list_all.dat b/Projects/s455/macro/yields/list_all.dat
new file mode 100644
index 0000000000000000000000000000000000000000..9aa355444ce2b18ca1077e67c1ecde333e36d4d7
--- /dev/null
+++ b/Projects/s455/macro/yields/list_all.dat
@@ -0,0 +1,118 @@
+77 172
+77 173
+77 174
+77 175
+78 174
+78 175
+78 176
+78 177
+78 178
+79 176
+79 177
+79 178
+79 179
+79 180
+79 181
+79 182
+80 179
+80 180
+80 181
+80 182
+80 183
+80 184
+80 185
+80 186
+81 183
+81 184
+81 185
+81 186
+81 187
+81 188
+82 186
+82 187
+82 188
+82 189
+82 190
+82 191
+82 192
+82 193
+83 189
+83 190
+83 191
+83 192
+83 193
+83 194
+83 195
+83 196
+83 197
+84 192
+84 193
+84 194
+84 195
+84 196
+84 197
+84 198
+84 199
+84 200
+84 201
+85 196
+85 197
+85 198
+85 199
+85 200
+85 201
+85 202
+85 203
+86 200
+86 201
+86 202
+86 203
+86 204
+86 205
+86 206
+87 203
+87 204
+87 205
+87 206
+87 207
+87 208
+87 209
+87 210
+87 211
+88 209
+88 210
+88 211
+88 212
+88 213
+88 214
+89 212
+89 213
+89 214
+89 215
+89 216
+89 217
+89 218
+90 216
+90 217
+90 219
+90 220
+
+89 219
+89 220
+89 223
+89 227
+90 221
+90 222
+90 223
+90 225
+90 226
+90 229
+90 230
+91 228
+91 229
+91 231
+91 232
+92 234
+92 235
+93 237
+93 238
diff --git a/Projects/s455/macro/yields/list_asy.dat b/Projects/s455/macro/yields/list_asy.dat
new file mode 100644
index 0000000000000000000000000000000000000000..9ee1ec8c747d81d98efc8b87088d0a09d68723fc
--- /dev/null
+++ b/Projects/s455/macro/yields/list_asy.dat
@@ -0,0 +1,58 @@
+77 175
+78 174
+78 175
+78 176
+78 177
+78 178
+79 176
+79 177
+79 178
+79 179
+79 180
+79 182
+80 179
+80 180
+80 181
+80 182
+80 183
+80 184
+80 185
+80 186
+81 183
+81 184
+81 185
+81 186
+81 187
+81 188
+82 186
+82 187
+82 188
+82 189
+82 190
+82 191
+82 192
+82 193
+90 216
+90 217
+90 219
+90 220
+
+89 219
+89 220
+89 223
+89 227
+90 221
+90 222
+90 223
+90 225
+90 226
+90 229
+90 230
+91 228
+91 229
+91 231
+91 232
+92 234
+92 235
+93 237
+93 238
diff --git a/Projects/s455/macro/yields/list_divers.dat b/Projects/s455/macro/yields/list_divers.dat
new file mode 100644
index 0000000000000000000000000000000000000000..1b009089ce7fcb311558fbdc30ae4ddbec199d32
--- /dev/null
+++ b/Projects/s455/macro/yields/list_divers.dat
@@ -0,0 +1,5 @@
+82 189
+81 185
+80 182
+79 177
+78 175
diff --git a/Projects/s455/macro/yields/list_nature.dat b/Projects/s455/macro/yields/list_nature.dat
new file mode 100644
index 0000000000000000000000000000000000000000..3d9ec8c01c4b54c8eda419521612ded9e622afca
--- /dev/null
+++ b/Projects/s455/macro/yields/list_nature.dat
@@ -0,0 +1,33 @@
+83 190
+83 191
+83 192
+83 193
+83 194
+83 195
+83 196
+82 186
+82 187
+82 188
+82 189
+82 190
+82 191
+82 192
+81 183
+81 184
+81 185
+81 186
+81 187
+81 188
+80 180
+80 181
+80 182
+80 183
+80 184
+80 185
+80 186
+79 177
+79 178
+79 179
+79 180
+79 181
+79 182
diff --git a/Projects/s455/macro/yields/list_s455.dat b/Projects/s455/macro/yields/list_s455.dat
new file mode 100644
index 0000000000000000000000000000000000000000..5af7d4cdafedb47ddf888ab64c6089c993253daf
--- /dev/null
+++ b/Projects/s455/macro/yields/list_s455.dat
@@ -0,0 +1,100 @@
+77 172
+77 173
+77 174
+77 175
+78 174
+78 175
+78 176
+78 177
+78 178
+79 176
+79 177
+79 178
+79 179
+79 180
+79 181
+79 182
+80 179
+80 180
+80 181
+80 182
+80 183
+80 184
+80 185
+80 186
+81 183
+81 184
+81 185
+81 186
+81 187
+81 188
+82 186
+82 187
+82 188
+82 189
+82 190
+82 191
+82 192
+82 193
+82 194
+83 189
+83 190
+83 191
+83 192
+83 193
+83 194
+83 195
+83 196
+83 197
+84 192
+84 193
+84 194
+84 195
+84 196
+84 197
+84 198
+84 199
+84 200
+84 201
+85 196
+85 197
+85 198
+85 199
+85 200
+85 201
+85 202
+85 203
+86 200
+86 201
+86 202
+86 203
+86 204
+86 205
+86 206
+87 203
+87 204
+87 205
+87 206
+87 207
+87 208
+87 209
+87 210
+87 211
+88 209
+88 210
+88 211
+88 212
+88 213
+88 214
+89 212
+89 213
+89 214
+89 215
+89 216
+89 217
+89 218
+90 216
+90 217
+90 219
+90 220
+
diff --git a/Projects/s455/macro/yields/list_s455_deformation.dat b/Projects/s455/macro/yields/list_s455_deformation.dat
new file mode 100644
index 0000000000000000000000000000000000000000..af0bd572886abe6f30bd0e788aac76fe88d012a6
--- /dev/null
+++ b/Projects/s455/macro/yields/list_s455_deformation.dat
@@ -0,0 +1,101 @@
+# deformation taken from https://www-phynu.cea.fr/science_en_ligne/carte_potentiels_microscopiques/carte_potentiel_nucleaire.htm
+# Z A beta
+77 172 0.28
+77 173 0.29
+77 174 0.30
+77 175 0.31
+78 174 0.3
+78 175 0.3
+78 176 0.3
+78 177 0.35
+78 178 0.35
+79 176 0.30
+79 177 0.30
+79 178 0.35
+79 179 0.35
+79 180 0.35
+79 181 0.35
+79 182 0.35
+80 180 0.30
+80 181 0.35
+80 182 0.35
+80 183 0.35
+80 184 0.30
+80 185 0.30
+80 186 -0.15
+81 183 0.30
+81 184 0.30
+81 185 -0.15
+81 186 0.30
+81 187 -0.15
+81 188 -0.15
+82 186 0
+82 187 0
+82 188 0
+82 189 0
+82 190 0
+82 191 0
+82 192 0
+82 193 0
+82 194 0
+83 189 0
+83 190 -0.2
+83 191 -0.2
+83 192 -0.2
+83 193 -0.1
+83 194 -0.1
+83 195 -0.1
+83 196 -0.1
+83 197 -0.1
+84 192 -0.2
+84 193 -0.2
+84 194 -0.2
+84 195 -0.2
+84 196 0.1
+84 197 0.1
+84 198 0.1
+84 199 0.1
+84 200 0.05
+84 201 0.01
+85 196 0.1
+85 197 0.1
+85 198 0.1
+85 199 0.1
+85 200 0.1
+85 201 0.1
+85 202 0.1
+85 203 0.1
+86 200 0.1
+86 201 0.1
+86 202 0.1
+86 203 0.1
+86 204 0.05
+86 205 0.05
+86 206 0.05
+87 203 0.15
+87 204 0.15
+87 205 0.1
+87 206 0.1
+87 207 0.05
+87 208 0.05
+87 209 0
+87 210 0
+87 211 0
+88 209 0.05
+88 210 0
+88 211 0
+88 212 0
+88 213 0
+88 214 0
+89 212 0
+89 213 0
+89 214 0
+89 215 0
+89 216 0
+89 217 0
+89 218 0
+90 216 0
+90 217 0
+90 219 0
+90 220 0.1
+
diff --git a/Projects/s455/macro/yields/list_sofia.dat b/Projects/s455/macro/yields/list_sofia.dat
new file mode 100644
index 0000000000000000000000000000000000000000..5af7d4cdafedb47ddf888ab64c6089c993253daf
--- /dev/null
+++ b/Projects/s455/macro/yields/list_sofia.dat
@@ -0,0 +1,100 @@
+77 172
+77 173
+77 174
+77 175
+78 174
+78 175
+78 176
+78 177
+78 178
+79 176
+79 177
+79 178
+79 179
+79 180
+79 181
+79 182
+80 179
+80 180
+80 181
+80 182
+80 183
+80 184
+80 185
+80 186
+81 183
+81 184
+81 185
+81 186
+81 187
+81 188
+82 186
+82 187
+82 188
+82 189
+82 190
+82 191
+82 192
+82 193
+82 194
+83 189
+83 190
+83 191
+83 192
+83 193
+83 194
+83 195
+83 196
+83 197
+84 192
+84 193
+84 194
+84 195
+84 196
+84 197
+84 198
+84 199
+84 200
+84 201
+85 196
+85 197
+85 198
+85 199
+85 200
+85 201
+85 202
+85 203
+86 200
+86 201
+86 202
+86 203
+86 204
+86 205
+86 206
+87 203
+87 204
+87 205
+87 206
+87 207
+87 208
+87 209
+87 210
+87 211
+88 209
+88 210
+88 211
+88 212
+88 213
+88 214
+89 212
+89 213
+89 214
+89 215
+89 216
+89 217
+89 218
+90 216
+90 217
+90 219
+90 220
+