diff --git a/NPSimulation/SSSD/SSSD.hh b/NPSimulation/SSSD/SSSD.hh
index 857b9b96094b8e55fa23f3159c591f4c982b2a09..df4e9ddfcab003fd53cb739f67bddb9db5eccf08 100644
--- a/NPSimulation/SSSD/SSSD.hh
+++ b/NPSimulation/SSSD/SSSD.hh
@@ -48,7 +48,7 @@ namespace SSSD_LOCAL{
    const G4double EnergyThreshold = 100*keV;
    // Geometry
    const G4double DetectorSize      = 68*mm           ;
-   const G4double SiliconThickness  = 15*micrometer   ;
+   const G4double SiliconThickness  = 20*micrometer   ;
    const G4double FrameThickness    = 3*mm            ;
    const G4double SiliconSize       = 50*mm           ;
    const G4double AluThickness      = 0.4*micrometer  ;
diff --git a/Projects/Nana/Eu.mac b/Projects/Nana/Eu.mac
new file mode 100644
index 0000000000000000000000000000000000000000..14f30e12ee17434838ee3e2fede71fe258105cbc
--- /dev/null
+++ b/Projects/Nana/Eu.mac
@@ -0,0 +1 @@
+/run/beamOn 38400000
diff --git a/Projects/Nana/NIM.cxx b/Projects/Nana/NIM.cxx
new file mode 100644
index 0000000000000000000000000000000000000000..319247c6a76201f2ca4870ecd5052d90337a9745
--- /dev/null
+++ b/Projects/Nana/NIM.cxx
@@ -0,0 +1,51 @@
+void NIM(){
+  TCanvas* c = new TCanvas("c","c",900,900);
+
+
+//  TFile* fexp = new TFile("Nana_exp.root","READ");
+//  TH1* hexp = (TH1*) fexp->FindObjectAny("Detector3");
+  TFile* fexp = new TFile("Eu_4hr.root","READ");
+  TTree* treex = (TTree*) fexp->FindObjectAny("Tree"); 
+  treex->Draw("fNANA_LaBr3_EnergyLong>>Detector3x(4096,0,4096)","fNANA_LaBr3_DetectorNbr@.size()==1 && fNANA_LaBr3_DetectorNbr==3","");
+  TH1* hexp = (TH1*) fexp->FindObjectAny("Detector3x");
+  TH1* hexp_o = (TH1*)hexp->Clone("hexp_o");
+
+  TFile* fbgd = new TFile("Bg_4hr_nu.root","READ");
+  TTree* treeb = (TTree*) fbgd->FindObjectAny("Tree"); 
+  treeb->Draw("fNANA_LaBr3_EnergyLong>>Detector3(4096,0,4096)","fNANA_LaBr3_DetectorNbr@.size()==1 && fNANA_LaBr3_DetectorNbr==3","same");
+  TH1* hbgd = (TH1*) fbgd->FindObjectAny("Detector3");
+
+  double min2= 1700;
+  double max2 = 3500;
+  double normalisation_bgd = hexp->Integral(hexp->FindBin(min2),hexp->FindBin(max2))/hbgd->Integral(hbgd->FindBin(min2),hbgd->FindBin(max2));
+  cout << "Normalisation = " << normalisation_bgd << endl;
+  hbgd->Scale(normalisation_bgd);
+  hexp->Add(hbgd,-1);
+
+  TFile* fsim = new TFile("../../Outputs/Analysis/Eu_5min.root","READ");
+  TTree* tree = (TTree*) fsim->FindObjectAny("PhysicsTree"); 
+ 
+  tree->Draw("LaBr_E*1000.>>sum2(4096,0,4096)","Nana.DetectorNumber@.size()==1 && Nana.DetectorNumber==3","E1 same");
+
+  TH1* hsim = (TH1*) fsim->FindObjectAny("sum2");
+  double min1= 400;
+  double max1 = 1400;
+  double normalisation = hexp->Integral(hexp->FindBin(min1),hexp->FindBin(max1))/hsim->Integral(hsim->FindBin(min1),hsim->FindBin(max1));
+  cout << "Normalisation = " << normalisation << endl;
+  hsim->Scale(normalisation);
+ 
+  hsim->SetLineColor(kOrange-3);
+  hsim->SetMarkerSize(0.7);
+  hexp->Draw();
+  hexp->GetXaxis()->SetTitle("LaBr3 Energy (keV)");
+  hexp->GetYaxis()->SetTitle("Counts per 4 keV");
+  hsim->Draw("same");
+  hbgd->SetFillColor(kRed);
+  hbgd->SetLineColor(kRed);
+  //hbgd->Draw("same");
+  hexp_o->SetFillColor(kGreen);
+  hexp_o->SetLineColor(kGreen);
+  //hexp_o->Draw("same"); 
+  
+  gPad->SetLogy();
+  }
diff --git a/Projects/Nana/PhysicsListOption.txt b/Projects/Nana/PhysicsListOption.txt
new file mode 100644
index 0000000000000000000000000000000000000000..58ffd6c8a1328ed252f511d65d607e949170f24d
--- /dev/null
+++ b/Projects/Nana/PhysicsListOption.txt
@@ -0,0 +1,9 @@
+EmPhysicsList Option4 
+DefaultCutOff 0.1 
+IonBinaryCascadePhysics 0
+EmExtraPhysics 0 
+HadronElasticPhysics 0
+StoppingPhysics 0
+OpticalPhysics 0
+HadronPhysicsQGSP_BIC_HP 0
+Decay 1