* Updatating Nana project

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  ;

Projects/Nana/Eu.mac

+/run/beamOn 38400000

Projects/Nana/NIM.cxx

+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();
+  }

Projects/Nana/PhysicsListOption.txt

+EmPhysicsList Option4 
+DefaultCutOff 0.1 
+IonBinaryCascadePhysics 0
+EmExtraPhysics 0 
+HadronElasticPhysics 0
+StoppingPhysics 0
+OpticalPhysics 0
+HadronPhysicsQGSP_BIC_HP 0
+Decay 1