Updating sofia simulation project

NPSimulation/Detectors/Sofia/SofTofW.cc

@@ -147,7 +147,7 @@ G4AssemblyVolume* SofTofW::BuildTOFDetector(){
     Tv.setZ(0);
     for(unsigned int i=0; i<28; i++){
       int k = -14+i;
-      Tv.setY(k*(SofTofW_NS::tof_plastic_width+0.5));
+      Tv.setY(k*(SofTofW_NS::tof_plastic_width+0.1));
       m_TofWall->AddPlacedVolume(m_PlasticTof, Tv, Rv);
     }
 

Projects/Sofia/macro/ShowResultSimu.C

@@ -8,6 +8,12 @@ void LoadRootFile(){
   chain->Add("../../../Outputs/Simulation/sofia_simu_2.root");
   chain->Add("../../../Outputs/Simulation/sofia_simu_3.root");
   chain->Add("../../../Outputs/Simulation/sofia_simu_4.root");
+  chain->Add("../../../Outputs/Simulation/sofia_simu_5.root");
+  chain->Add("../../../Outputs/Simulation/sofia_simu_6.root");
+  chain->Add("../../../Outputs/Simulation/sofia_simu_7.root");
+  chain->Add("../../../Outputs/Simulation/sofia_simu_8.root");
+  chain->Add("../../../Outputs/Simulation/sofia_simu_9.root");
+  chain->Add("../../../Outputs/Simulation/sofia_simu_10.root");
 }
 
 //////////////////////////////////////////////////////////////////
@@ -15,13 +21,13 @@ void ShowResultSimu()
 {
   LoadRootFile();
 
-  chain->Draw("fDetected_Position_Y:fDetected_Position_X>>hpos(500,-1700,-700,500,-350,350)","fTOF_Energy@.size()==2","colz");
+  chain->Draw("fDetected_Position_Y:fDetected_Position_X>>hpos(500,-1900,-900,500,-350,350)","fTOF_Energy@.size()>0","colz");
   TH2F* hpos = (TH2F*)gDirectory->FindObjectAny("hpos");
 
   chain->Draw("fFC_Fragment_Z>>h1(35,30,65)");
   TH1F* h1 = (TH1F*)gDirectory->FindObjectAny("h1");
 
-  chain->Draw("fDetected_Z>>h2(35,30,65)","fTOF_Energy@.size()==2");
+  chain->Draw("fDetected_Z>>h2(35,30,65)","fTOF_Energy@.size()>0");
   TH1F* h2 = (TH1F*)gDirectory->FindObjectAny("h2");
 
   h1->Sumw2();

Projects/Sofia/sofia.detector

@@ -12,17 +12,17 @@ SofTofW
  R= 4 m
  THETA= -20 deg
  PHI= 0 deg
- Build_GLAD= 0
+ Build_GLAD= 1
  GLAD_TiltAngle= 14 deg
  GLAD_DistanceFromTarget= 3.05 m
  Build_MagneticField= 1 
  GLAD_MagField= 1.70 T 
- Build_VacuumPipe= 1
+ Build_VacuumPipe= 0
  VacuumPipeX= 0 cm
  VacuumPipeY= 0 cm
  VacuumPipeZ= 0.9 m
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-SofTwim
+%SofTwim
  R= 1.85 m
  THETA= 0
  PHI= 0