diff --git a/Projects/NebulaPlus/.gitignore b/Projects/NebulaPlus/.gitignore
new file mode 100644
index 0000000000000000000000000000000000000000..0d62deb8e3ce2f9649e855aad5762a3c2ea6b88d
--- /dev/null
+++ b/Projects/NebulaPlus/.gitignore
@@ -0,0 +1,2 @@
+./root
+./raw_root
diff --git a/Projects/NebulaPlus/Analysis.cxx b/Projects/NebulaPlus/Analysis.cxx
new file mode 100644
index 0000000000000000000000000000000000000000..83eb160239b7abe57137db7b858c010d24590548
--- /dev/null
+++ b/Projects/NebulaPlus/Analysis.cxx
@@ -0,0 +1,83 @@
+/*****************************************************************************
+ * Copyright (C) 2009-2016 this file is part of the NPTool Project *
+ * *
+ * For the licensing terms see $NPTOOL/Licence/NPTool_Licence *
+ * For the list of contributors see $NPTOOL/Licence/Contributors *
+ *****************************************************************************/
+
+/*****************************************************************************
+ * Original Author: XAUTHORX contact address: XMAILX *
+ * *
+ * Creation Date : XMONTHX XYEARX *
+ * Last update : *
+ *---------------------------------------------------------------------------*
+ * Decription: *
+ * This class describe Nebula analysis project *
+ * *
+ *---------------------------------------------------------------------------*
+ * Comment: *
+ * *
+ *****************************************************************************/
+
+#include<iostream>
+using namespace std;
+#include"Analysis.h"
+#include"TGraph.h"
+#include"TCanvas.h"
+#include"TPad.h"
+#include"NPAnalysisFactory.h"
+#include"NPDetectorManager.h"
+////////////////////////////////////////////////////////////////////////////////
+Analysis::Analysis(){
+}
+////////////////////////////////////////////////////////////////////////////////
+Analysis::~Analysis(){
+}
+
+////////////////////////////////////////////////////////////////////////////////
+void Analysis::Init(){
+ NebulaPlus= (TNebulaPlusPhysics*) m_DetectorManager->GetDetector("NEBULAPLUS");
+}
+
+////////////////////////////////////////////////////////////////////////////////
+void Analysis::TreatEvent(){
+/*
+ TGraph* g = new TGraph();
+ for(int i=0; i<NebulaPlus->DetectorNumber.size(); i++){
+ //if(i==0)std::cout << "------------------" << std::endl;
+ //std::cout << NebulaPlus->GetPos(i).X() << std::endl;
+ g->SetPoint(i,NebulaPlus->GetPos(i).X(),NebulaPlus->GetPos(i).Y())
+ }
+*/
+// TGraph* g = new TGraph(NebulaPlus->DetectorNumber.size(),&(NebulaPlus->PosX[0]),&(NebulaPlus->PosY[0]));
+ // g->Fit("pol1");
+
+}
+
+
+////////////////////////////////////////////////////////////////////////////////
+void Analysis::End(){
+}
+
+
+////////////////////////////////////////////////////////////////////////////////
+// Construct Method to be pass to the DetectorFactory //
+////////////////////////////////////////////////////////////////////////////////
+NPL::VAnalysis* Analysis::Construct(){
+ return (NPL::VAnalysis*) new Analysis();
+}
+
+////////////////////////////////////////////////////////////////////////////////
+// Registering the construct method to the factory //
+////////////////////////////////////////////////////////////////////////////////
+extern "C"{
+class proxy{
+ public:
+ proxy(){
+ NPL::AnalysisFactory::getInstance()->SetConstructor(Analysis::Construct);
+ }
+};
+
+proxy p;
+}
+
diff --git a/Projects/NebulaPlus/Analysis.h b/Projects/NebulaPlus/Analysis.h
new file mode 100644
index 0000000000000000000000000000000000000000..ca51067ee431814ada236f4e582e205f2cf0456f
--- /dev/null
+++ b/Projects/NebulaPlus/Analysis.h
@@ -0,0 +1,42 @@
+#ifndef Analysis_h
+#define Analysis_h
+/*****************************************************************************
+ * Copyright (C) 2009-2016 this file is part of the NPTool Project *
+ * *
+ * For the licensing terms see $NPTOOL/Licence/NPTool_Licence *
+ * For the list of contributors see $NPTOOL/Licence/Contributors *
+ *****************************************************************************/
+
+/*****************************************************************************
+ * Original Author: XAUTHORX contact address: XMAILX *
+ * *
+ * Creation Date : XMONTHX XYEARX *
+ * Last update : *
+ *---------------------------------------------------------------------------*
+ * Decription: *
+ * This class describe Nebula analysis project *
+ * *
+ *---------------------------------------------------------------------------*
+ * Comment: *
+ * *
+ *****************************************************************************/
+
+#include"NPVAnalysis.h"
+#include"TNebulaPlusPhysics.h"
+class Analysis: public NPL::VAnalysis{
+ public:
+ Analysis();
+ ~Analysis();
+
+ public:
+ void Init();
+ void TreatEvent();
+ void End();
+
+ static NPL::VAnalysis* Construct();
+
+ private:
+ TNebulaPlusPhysics* NebulaPlus;
+
+};
+#endif
diff --git a/Projects/NebulaPlus/CMakeLists.txt b/Projects/NebulaPlus/CMakeLists.txt
new file mode 100644
index 0000000000000000000000000000000000000000..22c74affdfc45019bdda2594f8439c52d4ab97ec
--- /dev/null
+++ b/Projects/NebulaPlus/CMakeLists.txt
@@ -0,0 +1,5 @@
+cmake_minimum_required (VERSION 2.8)
+# Setting the policy to match Cmake version
+cmake_policy(VERSION ${CMAKE_MAJOR_VERSION}.${CMAKE_MINOR_VERSION})
+# include the default NPAnalysis cmake file
+include("../../NPLib/ressources/CMake/NPAnalysis.cmake")
diff --git a/Projects/NebulaPlus/Makefile b/Projects/NebulaPlus/Makefile
new file mode 100644
index 0000000000000000000000000000000000000000..7e183603a93b2331c3b57ac00feb94dede691b8e
--- /dev/null
+++ b/Projects/NebulaPlus/Makefile
@@ -0,0 +1,178 @@
+# CMAKE generated file: DO NOT EDIT!
+# Generated by "Unix Makefiles" Generator, CMake Version 3.16
+
+# Default target executed when no arguments are given to make.
+default_target: all
+
+.PHONY : default_target
+
+# Allow only one "make -f Makefile2" at a time, but pass parallelism.
+.NOTPARALLEL:
+
+
+#=============================================================================
+# Special targets provided by cmake.
+
+# Disable implicit rules so canonical targets will work.
+.SUFFIXES:
+
+
+# Remove some rules from gmake that .SUFFIXES does not remove.
+SUFFIXES =
+
+.SUFFIXES: .hpux_make_needs_suffix_list
+
+
+# Suppress display of executed commands.
+$(VERBOSE).SILENT:
+
+
+# A target that is always out of date.
+cmake_force:
+
+.PHONY : cmake_force
+
+#=============================================================================
+# Set environment variables for the build.
+
+# The shell in which to execute make rules.
+SHELL = /bin/sh
+
+# The CMake executable.
+CMAKE_COMMAND = /usr/bin/cmake
+
+# The command to remove a file.
+RM = /usr/bin/cmake -E remove -f
+
+# Escaping for special characters.
+EQUALS = =
+
+# The top-level source directory on which CMake was run.
+CMAKE_SOURCE_DIR = /home/expand/Experiments/Test_SAMURAI/nptool/Projects/NebulaPlus
+
+# The top-level build directory on which CMake was run.
+CMAKE_BINARY_DIR = /home/expand/Experiments/Test_SAMURAI/nptool/Projects/NebulaPlus
+
+#=============================================================================
+# Targets provided globally by CMake.
+
+# Special rule for the target rebuild_cache
+rebuild_cache:
+ @$(CMAKE_COMMAND) -E cmake_echo_color --switch=$(COLOR) --cyan "Running CMake to regenerate build system..."
+ /usr/bin/cmake -S$(CMAKE_SOURCE_DIR) -B$(CMAKE_BINARY_DIR)
+.PHONY : rebuild_cache
+
+# Special rule for the target rebuild_cache
+rebuild_cache/fast: rebuild_cache
+
+.PHONY : rebuild_cache/fast
+
+# Special rule for the target edit_cache
+edit_cache:
+ @$(CMAKE_COMMAND) -E cmake_echo_color --switch=$(COLOR) --cyan "No interactive CMake dialog available..."
+ /usr/bin/cmake -E echo No\ interactive\ CMake\ dialog\ available.
+.PHONY : edit_cache
+
+# Special rule for the target edit_cache
+edit_cache/fast: edit_cache
+
+.PHONY : edit_cache/fast
+
+# The main all target
+all: cmake_check_build_system
+ $(CMAKE_COMMAND) -E cmake_progress_start /home/expand/Experiments/Test_SAMURAI/nptool/Projects/NebulaPlus/CMakeFiles /home/expand/Experiments/Test_SAMURAI/nptool/Projects/NebulaPlus/CMakeFiles/progress.marks
+ $(MAKE) -f CMakeFiles/Makefile2 all
+ $(CMAKE_COMMAND) -E cmake_progress_start /home/expand/Experiments/Test_SAMURAI/nptool/Projects/NebulaPlus/CMakeFiles 0
+.PHONY : all
+
+# The main clean target
+clean:
+ $(MAKE) -f CMakeFiles/Makefile2 clean
+.PHONY : clean
+
+# The main clean target
+clean/fast: clean
+
+.PHONY : clean/fast
+
+# Prepare targets for installation.
+preinstall: all
+ $(MAKE) -f CMakeFiles/Makefile2 preinstall
+.PHONY : preinstall
+
+# Prepare targets for installation.
+preinstall/fast:
+ $(MAKE) -f CMakeFiles/Makefile2 preinstall
+.PHONY : preinstall/fast
+
+# clear depends
+depend:
+ $(CMAKE_COMMAND) -S$(CMAKE_SOURCE_DIR) -B$(CMAKE_BINARY_DIR) --check-build-system CMakeFiles/Makefile.cmake 1
+.PHONY : depend
+
+#=============================================================================
+# Target rules for targets named NPAnalysis
+
+# Build rule for target.
+NPAnalysis: cmake_check_build_system
+ $(MAKE) -f CMakeFiles/Makefile2 NPAnalysis
+.PHONY : NPAnalysis
+
+# fast build rule for target.
+NPAnalysis/fast:
+ $(MAKE) -f CMakeFiles/NPAnalysis.dir/build.make CMakeFiles/NPAnalysis.dir/build
+.PHONY : NPAnalysis/fast
+
+Analysis.o: Analysis.cxx.o
+
+.PHONY : Analysis.o
+
+# target to build an object file
+Analysis.cxx.o:
+ $(MAKE) -f CMakeFiles/NPAnalysis.dir/build.make CMakeFiles/NPAnalysis.dir/Analysis.cxx.o
+.PHONY : Analysis.cxx.o
+
+Analysis.i: Analysis.cxx.i
+
+.PHONY : Analysis.i
+
+# target to preprocess a source file
+Analysis.cxx.i:
+ $(MAKE) -f CMakeFiles/NPAnalysis.dir/build.make CMakeFiles/NPAnalysis.dir/Analysis.cxx.i
+.PHONY : Analysis.cxx.i
+
+Analysis.s: Analysis.cxx.s
+
+.PHONY : Analysis.s
+
+# target to generate assembly for a file
+Analysis.cxx.s:
+ $(MAKE) -f CMakeFiles/NPAnalysis.dir/build.make CMakeFiles/NPAnalysis.dir/Analysis.cxx.s
+.PHONY : Analysis.cxx.s
+
+# Help Target
+help:
+ @echo "The following are some of the valid targets for this Makefile:"
+ @echo "... all (the default if no target is provided)"
+ @echo "... clean"
+ @echo "... depend"
+ @echo "... rebuild_cache"
+ @echo "... edit_cache"
+ @echo "... NPAnalysis"
+ @echo "... Analysis.o"
+ @echo "... Analysis.i"
+ @echo "... Analysis.s"
+.PHONY : help
+
+
+
+#=============================================================================
+# Special targets to cleanup operation of make.
+
+# Special rule to run CMake to check the build system integrity.
+# No rule that depends on this can have commands that come from listfiles
+# because they might be regenerated.
+cmake_check_build_system:
+ $(CMAKE_COMMAND) -S$(CMAKE_SOURCE_DIR) -B$(CMAKE_BINARY_DIR) --check-build-system CMakeFiles/Makefile.cmake 0
+.PHONY : cmake_check_build_system
+
diff --git a/Projects/NebulaPlus/NebulaPlus.detector b/Projects/NebulaPlus/NebulaPlus.detector
new file mode 100644
index 0000000000000000000000000000000000000000..029a0a119feaa716a710b482d9e4e2f0c29650e5
--- /dev/null
+++ b/Projects/NebulaPlus/NebulaPlus.detector
@@ -0,0 +1,68 @@
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Target
+ THICKNESS= 0.49 cm
+ RADIUS= 28 mm
+ MATERIAL= Vacuum
+ ANGLE= 0 deg
+ X= 3.5 cm
+ Y= 0 mm
+ Z= 0 mm
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+NEBULAPLUS
+ POS= 0 0 0 cm
+ NumberOfModules= 22
+ Veto= 0
+ Frame= 0
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+NEBULAPLUS
+ POS= 0 0 13 cm
+ NumberOfModules= 22
+ Veto= 0
+ Frame= 0
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+NEBULAPLUS
+ POS= 0 0 84.5 cm
+ NumberOfModules= 23
+ Veto= 0
+ Frame= 0
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+NEBULAPLUS
+ POS= 0 0 97.5 cm
+ NumberOfModules= 25
+ Veto= 0
+ Frame= 0
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
+% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%Samurai
+% POS= 0 0 0 mm
+% ANGLE= 30 deg
+% %
+% METHOD= 1
+% % 0=RungeKutta, 1=EliaOmar
+% FIELDMAP=3T.table.bin
+% %fieldmap/3T.table.bin
+% % fieldmap path
+% STEPS_PER_METER= 1000
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%SAMURAIFDC1
+% %XML= db/SAMURAIFDC0_20200109.xml
+% %Offset= -0.00666226 0.102191 -3370.01 mm
+% Offset= -0.00666226 0.102191 -4370.01 mm
+% %InvertX= 1
+% %InvertY= 0
+% %InvertD= 1
+%
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%SAMURAIFDC2
+% %XML= db/SAMURAIFDC2.xml
+% Offset= 1000 200 5000 mm
+% OffAngle= 59.930 deg
+% %-252.416 -0.228477 4122.57 mm
+% %InvertX= 0
+% %InvertY= 1
+% %InvertD= 1
+
+
diff --git a/Projects/NebulaPlus/PhysicsListOption.txt b/Projects/NebulaPlus/PhysicsListOption.txt
new file mode 100644
index 0000000000000000000000000000000000000000..2a32a3b1722a1ec55d040d6e46faf8b513c453a0
--- /dev/null
+++ b/Projects/NebulaPlus/PhysicsListOption.txt
@@ -0,0 +1,11 @@
+EmPhysicsList Option4
+DefaultCutOff 1
+IonBinaryCascadePhysics 0
+NPIonInelasticPhysics 0
+EmExtraPhysics 0
+HadronElasticPhysics 0
+StoppingPhysics 0
+OpticalPhysics 0
+HadronPhysicsINCLXX 0
+HadronPhysicsQGSP_BIC_HP 1
+Decay 0
diff --git a/Projects/NebulaPlus/RunToTreat.txt b/Projects/NebulaPlus/RunToTreat.txt
new file mode 100755
index 0000000000000000000000000000000000000000..18635cc9c11300c5e957e63dfe44159135575752
--- /dev/null
+++ b/Projects/NebulaPlus/RunToTreat.txt
@@ -0,0 +1,4 @@
+TTreeName
+ RawTree
+RootFileName
+ ./raw_root/run_38_0001.root
diff --git a/Projects/NebulaPlus/SimpleNeutrons.reaction b/Projects/NebulaPlus/SimpleNeutrons.reaction
new file mode 100755
index 0000000000000000000000000000000000000000..6f5aea35fb4485aedc285ab8ba581ecc955b4c9a
--- /dev/null
+++ b/Projects/NebulaPlus/SimpleNeutrons.reaction
@@ -0,0 +1,14 @@
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Beam
+ Particle= n
+ Energy= 196 MeV
+ SigmaEnergy= 0.5 MeV
+ SigmaThetaX= 0.01 deg
+ SigmaPhiY= 0.01 deg
+ SigmaX= 0.5 mm
+ SigmaY= 0.5 mm
+ MeanThetaX= 0 deg
+ MeanPhiY= 0 deg
+ MeanX= 3.5 cm
+ MeanY= 0 mm
+ ZEmission= -1000 mm
diff --git a/Projects/NebulaPlus/calibration.txt b/Projects/NebulaPlus/calibration.txt
new file mode 100644
index 0000000000000000000000000000000000000000..3e5bd1d0f20981ae8758f2eba2ebc7bd84007ff1
--- /dev/null
+++ b/Projects/NebulaPlus/calibration.txt
@@ -0,0 +1,4 @@
+CalibrationFilePath
+./calibration/NebulaPlus_Up_Energy.cal
+./calibration/NebulaPlus_Down_Energy.cal
+./calibration/NebulaPlus_Tdiff_Offset.cal
diff --git a/Projects/NebulaPlus/calibration/NebulaPlus_Down_Energy.cal b/Projects/NebulaPlus/calibration/NebulaPlus_Down_Energy.cal
new file mode 100644
index 0000000000000000000000000000000000000000..dc432020a0080222dc32d6ab314e3ffc1e42830d
--- /dev/null
+++ b/Projects/NebulaPlus/calibration/NebulaPlus_Down_Energy.cal
@@ -0,0 +1,91 @@
+NebulaPlus_D1_ENERGY 0 1
+NebulaPlus_D2_ENERGY 0 1
+NebulaPlus_D3_ENERGY 0 1
+NebulaPlus_D4_ENERGY 0 1
+NebulaPlus_D5_ENERGY 0 1
+NebulaPlus_D6_ENERGY 0 1
+NebulaPlus_D7_ENERGY 0 1
+NebulaPlus_D8_ENERGY 0 1
+NebulaPlus_D9_ENERGY 0 1
+NebulaPlus_D10_ENERGY 0 1
+NebulaPlus_D11_ENERGY 0 1
+NebulaPlus_D12_ENERGY 0 1
+NebulaPlus_D13_ENERGY 0 1
+NebulaPlus_D14_ENERGY 0 1
+NebulaPlus_D15_ENERGY 0 1
+NebulaPlus_D16_ENERGY 0 1
+NebulaPlus_D17_ENERGY 0 1
+NebulaPlus_D18_ENERGY 0 1
+NebulaPlus_D19_ENERGY 0 1
+NebulaPlus_D20_ENERGY 0 1
+NebulaPlus_D21_ENERGY 0 1
+NebulaPlus_D22_ENERGY 0 1
+NebulaPlus_D23_ENERGY 0 1
+NebulaPlus_D24_ENERGY 0 1
+NebulaPlus_D25_ENERGY 0 1
+NebulaPlus_D26_ENERGY 0 1
+NebulaPlus_D27_ENERGY 0 1
+NebulaPlus_D28_ENERGY 0 1
+NebulaPlus_D29_ENERGY 0 1
+NebulaPlus_D30_ENERGY 0 1
+NebulaPlus_D31_ENERGY 0 1
+NebulaPlus_D32_ENERGY 0 1
+NebulaPlus_D33_ENERGY 0 1
+NebulaPlus_D34_ENERGY 0 1
+NebulaPlus_D35_ENERGY 0 1
+NebulaPlus_D36_ENERGY 0 1
+NebulaPlus_D37_ENERGY 0 1
+NebulaPlus_D38_ENERGY 0 1
+NebulaPlus_D39_ENERGY 0 1
+NebulaPlus_D40_ENERGY 0 1
+NebulaPlus_D41_ENERGY 0 1
+NebulaPlus_D42_ENERGY 0 1
+NebulaPlus_D43_ENERGY 0 1
+NebulaPlus_D44_ENERGY 0 1
+NebulaPlus_D45_ENERGY 0 1
+NebulaPlus_D46_ENERGY 0 1
+NebulaPlus_D47_ENERGY 0 1
+NebulaPlus_D48_ENERGY 0 1
+NebulaPlus_D49_ENERGY 0 1
+NebulaPlus_D50_ENERGY 0 1
+NebulaPlus_D51_ENERGY 0 1
+NebulaPlus_D52_ENERGY 0 1
+NebulaPlus_D53_ENERGY 0 1
+NebulaPlus_D54_ENERGY 0 1
+NebulaPlus_D55_ENERGY 0 1
+NebulaPlus_D56_ENERGY 0 1
+NebulaPlus_D57_ENERGY 0 1
+NebulaPlus_D58_ENERGY 0 1
+NebulaPlus_D59_ENERGY 0 1
+NebulaPlus_D60_ENERGY 0 1
+NebulaPlus_D61_ENERGY 0 1
+NebulaPlus_D62_ENERGY 0 1
+NebulaPlus_D63_ENERGY 0 1
+NebulaPlus_D64_ENERGY 0 1
+NebulaPlus_D65_ENERGY 0 1
+NebulaPlus_D66_ENERGY 0 1
+NebulaPlus_D67_ENERGY 0 1
+NebulaPlus_D68_ENERGY 0 1
+NebulaPlus_D69_ENERGY 0 1
+NebulaPlus_D70_ENERGY 0 1
+NebulaPlus_D71_ENERGY 0 1
+NebulaPlus_D72_ENERGY 0 1
+NebulaPlus_D73_ENERGY 0 1
+NebulaPlus_D74_ENERGY 0 1
+NebulaPlus_D75_ENERGY 0 1
+NebulaPlus_D76_ENERGY 0 1
+NebulaPlus_D77_ENERGY 0 1
+NebulaPlus_D78_ENERGY 0 1
+NebulaPlus_D79_ENERGY 0 1
+NebulaPlus_D80_ENERGY 0 1
+NebulaPlus_D81_ENERGY 0 1
+NebulaPlus_D82_ENERGY 0 1
+NebulaPlus_D83_ENERGY 0 1
+NebulaPlus_D84_ENERGY 0 1
+NebulaPlus_D85_ENERGY 0 1
+NebulaPlus_D86_ENERGY 0 1
+NebulaPlus_D87_ENERGY 0 1
+NebulaPlus_D88_ENERGY 0 1
+NebulaPlus_D89_ENERGY 0 1
+NebulaPlus_D90_ENERGY 0 1
+
diff --git a/Projects/NebulaPlus/calibration/NebulaPlus_Energy.cal b/Projects/NebulaPlus/calibration/NebulaPlus_Energy.cal
new file mode 100644
index 0000000000000000000000000000000000000000..77372aef6c2ff5914e4c177473c50884ecc713ad
--- /dev/null
+++ b/Projects/NebulaPlus/calibration/NebulaPlus_Energy.cal
@@ -0,0 +1,91 @@
+NebulaPlus_U1_ENERGY 0 1
+NebulaPlus_U2_ENERGY 0 1
+NebulaPlus_U3_ENERGY 0 1
+NebulaPlus_U4_ENERGY 0 1
+NebulaPlus_U5_ENERGY 0 1
+NebulaPlus_U6_ENERGY 0 1
+NebulaPlus_U7_ENERGY 0 1
+NebulaPlus_U8_ENERGY 0 1
+NebulaPlus_U9_ENERGY 0 1
+NebulaPlus_U10_ENERGY 0 1
+NebulaPlus_U11_ENERGY 0 1
+NebulaPlus_U12_ENERGY 0 1
+NebulaPlus_U13_ENERGY 0 1
+NebulaPlus_U14_ENERGY 0 1
+NebulaPlus_U15_ENERGY 0 1
+NebulaPlus_U16_ENERGY 0 1
+NebulaPlus_U17_ENERGY 0 1
+NebulaPlus_U18_ENERGY 0 1
+NebulaPlus_U19_ENERGY 0 1
+NebulaPlus_U20_ENERGY 0 1
+NebulaPlus_U21_ENERGY 0 1
+NebulaPlus_U22_ENERGY 0 1
+NebulaPlus_U23_ENERGY 0 1
+NebulaPlus_U24_ENERGY 0 1
+NebulaPlus_U25_ENERGY 0 1
+NebulaPlus_U26_ENERGY 0 1
+NebulaPlus_U27_ENERGY 0 1
+NebulaPlus_U28_ENERGY 0 1
+NebulaPlus_U29_ENERGY 0 1
+NebulaPlus_U30_ENERGY 0 1
+NebulaPlus_U31_ENERGY 0 1
+NebulaPlus_U32_ENERGY 0 1
+NebulaPlus_U33_ENERGY 0 1
+NebulaPlus_U34_ENERGY 0 1
+NebulaPlus_U35_ENERGY 0 1
+NebulaPlus_U36_ENERGY 0 1
+NebulaPlus_U37_ENERGY 0 1
+NebulaPlus_U38_ENERGY 0 1
+NebulaPlus_U39_ENERGY 0 1
+NebulaPlus_U40_ENERGY 0 1
+NebulaPlus_U41_ENERGY 0 1
+NebulaPlus_U42_ENERGY 0 1
+NebulaPlus_U43_ENERGY 0 1
+NebulaPlus_U44_ENERGY 0 1
+NebulaPlus_U45_ENERGY 0 1
+NebulaPlus_U46_ENERGY 0 1
+NebulaPlus_U47_ENERGY 0 1
+NebulaPlus_U48_ENERGY 0 1
+NebulaPlus_U49_ENERGY 0 1
+NebulaPlus_U50_ENERGY 0 1
+NebulaPlus_U51_ENERGY 0 1
+NebulaPlus_U52_ENERGY 0 1
+NebulaPlus_U53_ENERGY 0 1
+NebulaPlus_U54_ENERGY 0 1
+NebulaPlus_U55_ENERGY 0 1
+NebulaPlus_U56_ENERGY 0 1
+NebulaPlus_U57_ENERGY 0 1
+NebulaPlus_U58_ENERGY 0 1
+NebulaPlus_U59_ENERGY 0 1
+NebulaPlus_U60_ENERGY 0 1
+NebulaPlus_U61_ENERGY 0 1
+NebulaPlus_U62_ENERGY 0 1
+NebulaPlus_U63_ENERGY 0 1
+NebulaPlus_U64_ENERGY 0 1
+NebulaPlus_U65_ENERGY 0 1
+NebulaPlus_U66_ENERGY 0 1
+NebulaPlus_U67_ENERGY 0 1
+NebulaPlus_U68_ENERGY 0 1
+NebulaPlus_U69_ENERGY 0 1
+NebulaPlus_U70_ENERGY 0 1
+NebulaPlus_U71_ENERGY 0 1
+NebulaPlus_U72_ENERGY 0 1
+NebulaPlus_U73_ENERGY 0 1
+NebulaPlus_U74_ENERGY 0 1
+NebulaPlus_U75_ENERGY 0 1
+NebulaPlus_U76_ENERGY 0 1
+NebulaPlus_U77_ENERGY 0 1
+NebulaPlus_U78_ENERGY 0 1
+NebulaPlus_U79_ENERGY 0 1
+NebulaPlus_U80_ENERGY 0 1
+NebulaPlus_U81_ENERGY 0 1
+NebulaPlus_U82_ENERGY 0 1
+NebulaPlus_U83_ENERGY 0 1
+NebulaPlus_U84_ENERGY 0 1
+NebulaPlus_U85_ENERGY 0 1
+NebulaPlus_U86_ENERGY 0 1
+NebulaPlus_U87_ENERGY 0 1
+NebulaPlus_U88_ENERGY 0 1
+NebulaPlus_U89_ENERGY 0 1
+NebulaPlus_U90_ENERGY 0 1
+
diff --git a/Projects/NebulaPlus/calibration/NebulaPlus_Tdiff_Offset.cal b/Projects/NebulaPlus/calibration/NebulaPlus_Tdiff_Offset.cal
new file mode 100644
index 0000000000000000000000000000000000000000..24c21ef2b51b23ca4bf7124bf914ab7fd78adef7
--- /dev/null
+++ b/Projects/NebulaPlus/calibration/NebulaPlus_Tdiff_Offset.cal
@@ -0,0 +1,92 @@
+NebulaPlus_Tdiff_Offset_1 5.17288
+NebulaPlus_Tdiff_Offset_2 0.568133
+NebulaPlus_Tdiff_Offset_3 2.47937
+NebulaPlus_Tdiff_Offset_4 3.43329
+NebulaPlus_Tdiff_Offset_5 2.26801
+NebulaPlus_Tdiff_Offset_6 -3.03789
+NebulaPlus_Tdiff_Offset_7 4.90289
+NebulaPlus_Tdiff_Offset_8 3.9337
+NebulaPlus_Tdiff_Offset_9 -4.71495
+NebulaPlus_Tdiff_Offset_10 1.5473
+NebulaPlus_Tdiff_Offset_11 3.9492
+NebulaPlus_Tdiff_Offset_12 0.99885
+NebulaPlus_Tdiff_Offset_13 -5.47822
+NebulaPlus_Tdiff_Offset_14 5.5472
+NebulaPlus_Tdiff_Offset_15 2.048
+NebulaPlus_Tdiff_Offset_16 -0.440609
+NebulaPlus_Tdiff_Offset_17 -1.87889
+NebulaPlus_Tdiff_Offset_18 1.06144
+NebulaPlus_Tdiff_Offset_19 2.83913
+NebulaPlus_Tdiff_Offset_20 2.55256
+NebulaPlus_Tdiff_Offset_21 1.30455
+NebulaPlus_Tdiff_Offset_22 0.996235
+NebulaPlus_Tdiff_Offset_23 4.05148
+NebulaPlus_Tdiff_Offset_24 -5.01333
+NebulaPlus_Tdiff_Offset_25 -5.04853
+NebulaPlus_Tdiff_Offset_26 -6.92471
+NebulaPlus_Tdiff_Offset_27 -6.88351
+NebulaPlus_Tdiff_Offset_28 0.96321
+NebulaPlus_Tdiff_Offset_29 3.12615
+NebulaPlus_Tdiff_Offset_30 -5.6342
+NebulaPlus_Tdiff_Offset_31 -0.798026
+NebulaPlus_Tdiff_Offset_32 0.287694
+NebulaPlus_Tdiff_Offset_33 1.66875
+NebulaPlus_Tdiff_Offset_34 -6.99295
+NebulaPlus_Tdiff_Offset_35 -3.19359
+NebulaPlus_Tdiff_Offset_36 0.02765
+NebulaPlus_Tdiff_Offset_37 -0.880016
+NebulaPlus_Tdiff_Offset_38 4.7604
+NebulaPlus_Tdiff_Offset_39 -3.06803
+NebulaPlus_Tdiff_Offset_40 0.907779
+NebulaPlus_Tdiff_Offset_41 -4.67336
+NebulaPlus_Tdiff_Offset_42 0.829625
+NebulaPlus_Tdiff_Offset_43 0.28951
+NebulaPlus_Tdiff_Offset_44 3.63011
+NebulaPlus_Tdiff_Offset_45 0.932612
+NebulaPlus_Tdiff_Offset_46 -7.92853
+NebulaPlus_Tdiff_Offset_47 -0.989877
+NebulaPlus_Tdiff_Offset_48 3.52937
+NebulaPlus_Tdiff_Offset_49 0.236056
+NebulaPlus_Tdiff_Offset_50 -2.87001
+NebulaPlus_Tdiff_Offset_51 1.85101
+NebulaPlus_Tdiff_Offset_52 -4.24166
+NebulaPlus_Tdiff_Offset_53 4.159
+NebulaPlus_Tdiff_Offset_54 5.85506
+NebulaPlus_Tdiff_Offset_55 3.56968
+NebulaPlus_Tdiff_Offset_56 -0.89544
+NebulaPlus_Tdiff_Offset_57 -0.826875
+NebulaPlus_Tdiff_Offset_58 4.09274
+NebulaPlus_Tdiff_Offset_59 -1.75094
+NebulaPlus_Tdiff_Offset_60 0.493957
+NebulaPlus_Tdiff_Offset_61 -0.260756
+NebulaPlus_Tdiff_Offset_62 6.46713
+NebulaPlus_Tdiff_Offset_63 -4.96601
+NebulaPlus_Tdiff_Offset_64 -5.46893
+NebulaPlus_Tdiff_Offset_65 11.4451
+NebulaPlus_Tdiff_Offset_66 1.35551
+NebulaPlus_Tdiff_Offset_67 1.1528
+NebulaPlus_Tdiff_Offset_68 0.0792745
+NebulaPlus_Tdiff_Offset_69 -0.0368659
+NebulaPlus_Tdiff_Offset_70 5.84517
+NebulaPlus_Tdiff_Offset_71 -2.02714
+NebulaPlus_Tdiff_Offset_72 3.68631
+NebulaPlus_Tdiff_Offset_73 -5.92629
+NebulaPlus_Tdiff_Offset_74 -3.25264
+NebulaPlus_Tdiff_Offset_75 -4.746
+NebulaPlus_Tdiff_Offset_76 3.08691
+NebulaPlus_Tdiff_Offset_77 -5.22465
+NebulaPlus_Tdiff_Offset_78 1.80784
+NebulaPlus_Tdiff_Offset_79 3.56686
+NebulaPlus_Tdiff_Offset_80 -5.3148
+NebulaPlus_Tdiff_Offset_81 -6.70568
+NebulaPlus_Tdiff_Offset_82 -0.436624
+NebulaPlus_Tdiff_Offset_83 3.85383
+NebulaPlus_Tdiff_Offset_84 4.01968
+NebulaPlus_Tdiff_Offset_85 0.279455
+NebulaPlus_Tdiff_Offset_86 0.456427
+NebulaPlus_Tdiff_Offset_87 3.23678
+NebulaPlus_Tdiff_Offset_88 0.784784
+NebulaPlus_Tdiff_Offset_89 -5.3909
+NebulaPlus_Tdiff_Offset_90 0
+NebulaPlus_Tdiff_Offset_91 3.12457
+NebulaPlus_Tdiff_Offset_92 0
diff --git a/Projects/NebulaPlus/calibration/NebulaPlus_Up_Energy.cal b/Projects/NebulaPlus/calibration/NebulaPlus_Up_Energy.cal
new file mode 100644
index 0000000000000000000000000000000000000000..77372aef6c2ff5914e4c177473c50884ecc713ad
--- /dev/null
+++ b/Projects/NebulaPlus/calibration/NebulaPlus_Up_Energy.cal
@@ -0,0 +1,91 @@
+NebulaPlus_U1_ENERGY 0 1
+NebulaPlus_U2_ENERGY 0 1
+NebulaPlus_U3_ENERGY 0 1
+NebulaPlus_U4_ENERGY 0 1
+NebulaPlus_U5_ENERGY 0 1
+NebulaPlus_U6_ENERGY 0 1
+NebulaPlus_U7_ENERGY 0 1
+NebulaPlus_U8_ENERGY 0 1
+NebulaPlus_U9_ENERGY 0 1
+NebulaPlus_U10_ENERGY 0 1
+NebulaPlus_U11_ENERGY 0 1
+NebulaPlus_U12_ENERGY 0 1
+NebulaPlus_U13_ENERGY 0 1
+NebulaPlus_U14_ENERGY 0 1
+NebulaPlus_U15_ENERGY 0 1
+NebulaPlus_U16_ENERGY 0 1
+NebulaPlus_U17_ENERGY 0 1
+NebulaPlus_U18_ENERGY 0 1
+NebulaPlus_U19_ENERGY 0 1
+NebulaPlus_U20_ENERGY 0 1
+NebulaPlus_U21_ENERGY 0 1
+NebulaPlus_U22_ENERGY 0 1
+NebulaPlus_U23_ENERGY 0 1
+NebulaPlus_U24_ENERGY 0 1
+NebulaPlus_U25_ENERGY 0 1
+NebulaPlus_U26_ENERGY 0 1
+NebulaPlus_U27_ENERGY 0 1
+NebulaPlus_U28_ENERGY 0 1
+NebulaPlus_U29_ENERGY 0 1
+NebulaPlus_U30_ENERGY 0 1
+NebulaPlus_U31_ENERGY 0 1
+NebulaPlus_U32_ENERGY 0 1
+NebulaPlus_U33_ENERGY 0 1
+NebulaPlus_U34_ENERGY 0 1
+NebulaPlus_U35_ENERGY 0 1
+NebulaPlus_U36_ENERGY 0 1
+NebulaPlus_U37_ENERGY 0 1
+NebulaPlus_U38_ENERGY 0 1
+NebulaPlus_U39_ENERGY 0 1
+NebulaPlus_U40_ENERGY 0 1
+NebulaPlus_U41_ENERGY 0 1
+NebulaPlus_U42_ENERGY 0 1
+NebulaPlus_U43_ENERGY 0 1
+NebulaPlus_U44_ENERGY 0 1
+NebulaPlus_U45_ENERGY 0 1
+NebulaPlus_U46_ENERGY 0 1
+NebulaPlus_U47_ENERGY 0 1
+NebulaPlus_U48_ENERGY 0 1
+NebulaPlus_U49_ENERGY 0 1
+NebulaPlus_U50_ENERGY 0 1
+NebulaPlus_U51_ENERGY 0 1
+NebulaPlus_U52_ENERGY 0 1
+NebulaPlus_U53_ENERGY 0 1
+NebulaPlus_U54_ENERGY 0 1
+NebulaPlus_U55_ENERGY 0 1
+NebulaPlus_U56_ENERGY 0 1
+NebulaPlus_U57_ENERGY 0 1
+NebulaPlus_U58_ENERGY 0 1
+NebulaPlus_U59_ENERGY 0 1
+NebulaPlus_U60_ENERGY 0 1
+NebulaPlus_U61_ENERGY 0 1
+NebulaPlus_U62_ENERGY 0 1
+NebulaPlus_U63_ENERGY 0 1
+NebulaPlus_U64_ENERGY 0 1
+NebulaPlus_U65_ENERGY 0 1
+NebulaPlus_U66_ENERGY 0 1
+NebulaPlus_U67_ENERGY 0 1
+NebulaPlus_U68_ENERGY 0 1
+NebulaPlus_U69_ENERGY 0 1
+NebulaPlus_U70_ENERGY 0 1
+NebulaPlus_U71_ENERGY 0 1
+NebulaPlus_U72_ENERGY 0 1
+NebulaPlus_U73_ENERGY 0 1
+NebulaPlus_U74_ENERGY 0 1
+NebulaPlus_U75_ENERGY 0 1
+NebulaPlus_U76_ENERGY 0 1
+NebulaPlus_U77_ENERGY 0 1
+NebulaPlus_U78_ENERGY 0 1
+NebulaPlus_U79_ENERGY 0 1
+NebulaPlus_U80_ENERGY 0 1
+NebulaPlus_U81_ENERGY 0 1
+NebulaPlus_U82_ENERGY 0 1
+NebulaPlus_U83_ENERGY 0 1
+NebulaPlus_U84_ENERGY 0 1
+NebulaPlus_U85_ENERGY 0 1
+NebulaPlus_U86_ENERGY 0 1
+NebulaPlus_U87_ENERGY 0 1
+NebulaPlus_U88_ENERGY 0 1
+NebulaPlus_U89_ENERGY 0 1
+NebulaPlus_U90_ENERGY 0 1
+
diff --git a/Projects/NebulaPlus/cmake_install.cmake b/Projects/NebulaPlus/cmake_install.cmake
new file mode 100644
index 0000000000000000000000000000000000000000..93b19c768b6e7f17709bdd02c9caf0564cf60fb9
--- /dev/null
+++ b/Projects/NebulaPlus/cmake_install.cmake
@@ -0,0 +1,49 @@
+# Install script for directory: /home/expand/Experiments/Test_SAMURAI/nptool/Projects/NebulaPlus
+
+# Set the install prefix
+if(NOT DEFINED CMAKE_INSTALL_PREFIX)
+ set(CMAKE_INSTALL_PREFIX "/usr/local")
+endif()
+string(REGEX REPLACE "/$" "" CMAKE_INSTALL_PREFIX "${CMAKE_INSTALL_PREFIX}")
+
+# Set the install configuration name.
+if(NOT DEFINED CMAKE_INSTALL_CONFIG_NAME)
+ if(BUILD_TYPE)
+ string(REGEX REPLACE "^[^A-Za-z0-9_]+" ""
+ CMAKE_INSTALL_CONFIG_NAME "${BUILD_TYPE}")
+ else()
+ set(CMAKE_INSTALL_CONFIG_NAME "Release")
+ endif()
+ message(STATUS "Install configuration: \"${CMAKE_INSTALL_CONFIG_NAME}\"")
+endif()
+
+# Set the component getting installed.
+if(NOT CMAKE_INSTALL_COMPONENT)
+ if(COMPONENT)
+ message(STATUS "Install component: \"${COMPONENT}\"")
+ set(CMAKE_INSTALL_COMPONENT "${COMPONENT}")
+ else()
+ set(CMAKE_INSTALL_COMPONENT)
+ endif()
+endif()
+
+# Install shared libraries without execute permission?
+if(NOT DEFINED CMAKE_INSTALL_SO_NO_EXE)
+ set(CMAKE_INSTALL_SO_NO_EXE "1")
+endif()
+
+# Is this installation the result of a crosscompile?
+if(NOT DEFINED CMAKE_CROSSCOMPILING)
+ set(CMAKE_CROSSCOMPILING "FALSE")
+endif()
+
+if(CMAKE_INSTALL_COMPONENT)
+ set(CMAKE_INSTALL_MANIFEST "install_manifest_${CMAKE_INSTALL_COMPONENT}.txt")
+else()
+ set(CMAKE_INSTALL_MANIFEST "install_manifest.txt")
+endif()
+
+string(REPLACE ";" "\n" CMAKE_INSTALL_MANIFEST_CONTENT
+ "${CMAKE_INSTALL_MANIFEST_FILES}")
+file(WRITE "/home/expand/Experiments/Test_SAMURAI/nptool/Projects/NebulaPlus/${CMAKE_INSTALL_MANIFEST}"
+ "${CMAKE_INSTALL_MANIFEST_CONTENT}")
diff --git a/Projects/NebulaPlus/energy_loss/.gitignore b/Projects/NebulaPlus/energy_loss/.gitignore
new file mode 100644
index 0000000000000000000000000000000000000000..e69de29bb2d1d6434b8b29ae775ad8c2e48c5391
diff --git a/Projects/NebulaPlus/flat.txt b/Projects/NebulaPlus/flat.txt
new file mode 100644
index 0000000000000000000000000000000000000000..09b15128730aecce2f1126bf7937af827eba3712
--- /dev/null
+++ b/Projects/NebulaPlus/flat.txt
@@ -0,0 +1,3 @@
+0 1
+90 1
+180 1
diff --git a/Projects/NebulaPlus/macro/.gitignore b/Projects/NebulaPlus/macro/.gitignore
new file mode 100644
index 0000000000000000000000000000000000000000..e69de29bb2d1d6434b8b29ae775ad8c2e48c5391
diff --git a/Projects/NebulaPlus/macro/FitTrack.C b/Projects/NebulaPlus/macro/FitTrack.C
new file mode 100644
index 0000000000000000000000000000000000000000..3801c4b3e51b43c338f9c1ad5b10443962eaf630
--- /dev/null
+++ b/Projects/NebulaPlus/macro/FitTrack.C
@@ -0,0 +1,157 @@
+TChain* chain=NULL ;
+
+////////////////////////////////////////////////////////////////////////////////
+void LoadChain(){
+ chain = new TChain("PhysicsTree");
+ //chain->Add("./root/analysis/testmult16.root");
+ chain->Add("./root/analysis/run_38_0001.root");
+}
+double fpol1(double* x, double* par){
+ double result = par[0]+par[1]*x[0];
+ return result;
+}
+
+////////////////////////////////////////////////////////////////////////////////
+void FitTrack(){
+ LoadChain();
+
+ TNebulaPlusPhysics* m_NebulaPlus = new TNebulaPlusPhysics();
+ chain->SetBranchAddress("NebulaPlus", &m_NebulaPlus);
+ Long64_t nentries = chain->GetEntries();
+ TH1F* h = new TH1F("h","h",300,0,300);
+ TH1F* hTheta = new TH1F("hTheta","hTheta",100,0,100);
+ TH1F** hdT = new TH1F*[90];
+ for ( int i=0; i<90; i++) {
+ hdT[i]=new TH1F(Form("hdT%d",i+1),Form("TOF%d - TOF%d:TOF(ns)",i+1,i+2),400,-10,10);
+ }
+
+ TCanvas* c = new TCanvas("c","c",0,0,700,700);
+ c->cd();
+ gStyle->SetOptFit(111111111);
+ h->GetXaxis()->SetRangeUser(-300,300);
+ h->GetYaxis()->SetRangeUser(-300,300);
+ double Theta;
+
+ for(Long64_t i = 0; i < nentries; i++){
+ chain->GetEntry(i);
+ int mult= m_NebulaPlus->DetectorNumber.size();
+
+ /// build flag to select single layer event ///
+ double tmp=-99;
+ bool multilayer=false;
+ for(int j=0;j<mult;j++){
+ if(j==0) tmp = m_NebulaPlus->PosZ[j];
+ if(m_NebulaPlus->PosZ[j]!=tmp) multilayer=true;
+ }
+
+ if(!multilayer){
+ if(mult>6){
+
+ ///////////////////////////////////////////////
+ //// Track construction and Fit in Y vs X /////
+ ///////////////////////////////////////////////
+ TGraph* g = new TGraph(
+ mult,
+ &(m_NebulaPlus->PosX[0]),
+ &(m_NebulaPlus->PosY[0])
+ );
+ g->SetMarkerStyle(3);
+ h->Draw("");
+ g->Draw("Psame");
+
+ TF1* toto = new TF1("f1",fpol1,0,300,2);
+ g->Fit("f1","QR");
+ double offset = toto->GetParameter(0);
+ double slope = toto->GetParameter(1);
+
+ //std::cout << offset << std::endl;;
+ //std::cout << slope << std::endl;;
+
+ TVector2 a(0,offset);
+ TVector2 b(1,offset+slope);
+ TVector2 c = b-a;
+ Theta = fabs(TMath::ATan(c.Y()/c.X())*180./TMath::Pi());
+ std::cout << Theta << std::endl;
+ hTheta->Fill(Theta);
+ gPad->WaitPrimitive();
+
+ ///////////////////////////////////////////////
+ //// Time analysis between adjacent bars //////
+ ///////////////////////////////////////////////
+
+ vector<int> ID(mult);
+ vector<int> IDsort(mult);
+ vector<int> idx(mult);
+ vector<double> T(mult);
+ double dT=0;
+
+ ID = m_NebulaPlus->DetectorNumber;
+
+ //std::cout << "---------" << std::endl;
+ //std::cout << ID[0] << std::endl;
+ //std::cout << m_NebulaPlus->DetectorNumber[0] << std::endl;
+ TMath::Sort(mult,&ID[0],&idx[0],false);
+
+ //std::cout << "---------" << std::endl;
+
+ for(int k=0;k<mult;k++){
+ T[k] = m_NebulaPlus->TOF[idx[k]];
+ IDsort[k] = m_NebulaPlus->DetectorNumber[idx[k]];
+ //if(k=0) refT = T[0]; refID = ID[0];
+ if(k>0 && IDsort[k-1]==IDsort[k]-1 ){
+ // if(T[k]>T[k-1]) dT = T[k]-T[k-1];
+ // if(T[k]<T[k-1]) dT = T[k-1]-T[k];
+
+ dT = T[k]-T[k-1];
+ hdT[IDsort[k-1]]->Fill(dT);
+
+ //std::cout <<IDsort[k]<< "\t" << hdT[IDsort[k-1]]->GetName() << "\t" << dT<< std::endl;
+ }
+ }
+
+
+
+ }
+ }
+
+ }
+
+ c->cd();
+ hTheta->Draw();
+ c->Print("TOF.pdf[");
+ TGraph *res1 = new TGraph(90);
+ TGraph *res2 = new TGraph(90);
+ res1->SetName("res1");
+ res2->SetName("res2");
+ for (int i=0; i<90; i++){
+ TF1* f_gauss = new TF1("f_gauss","gaus(0)+gaus(3)",-10,10);
+ double par[6] = {hdT[i]->GetMaximum(), hdT[i]->GetMean()-.5, 0.5,hdT[i]->GetMaximum(), hdT[i]->GetMean()+0.5 , 0.5 };
+ f_gauss->SetParameters(par);
+ f_gauss->SetNpx(1000);
+ hdT[i]->Fit(f_gauss,"Q");
+ c->Print("TOF.pdf");
+ res1->SetPoint(i,i+1,TMath::Abs(f_gauss->GetParameter(2)));
+ res2->SetPoint(i,i+1,TMath::Abs(f_gauss->GetParameter(5)));
+
+ }
+ res1->SetMarkerStyle(kStar);
+ res1->SetMarkerColor(kRed);
+
+
+ res2->SetMarkerStyle(kStar);
+ res2->SetMarkerColor(kBlue);
+ TMultiGraph* mg = new TMultiGraph();
+ mg->Add(res1,"p");
+ mg->Add(res2,"p");
+ mg->SetTitle(" ; ID ; sigma value");
+ mg->Draw("a");
+
+
+ c->Print("TOF.pdf");
+ c->Print("TOF.pdf]");
+
+
+
+
+
+}
diff --git a/Projects/NebulaPlus/macro/ShowTrack.C b/Projects/NebulaPlus/macro/ShowTrack.C
new file mode 100644
index 0000000000000000000000000000000000000000..d62e352c7c7260d594aa636324647e032f08e4b7
--- /dev/null
+++ b/Projects/NebulaPlus/macro/ShowTrack.C
@@ -0,0 +1,40 @@
+TChain* chain=NULL ;
+
+////////////////////////////////////////////////////////////////////////////////
+void LoadChain(){
+ chain = new TChain("PhysicsTree");
+ chain->Add("./root/analysis/testmult16.root");
+}
+
+////////////////////////////////////////////////////////////////////////////////
+void ShowTrack(){
+ LoadChain();
+ TH1 *h[20];
+ TGraph *g[20];
+ string command;
+ string histname;
+
+ TCanvas* c1 = new TCanvas("YX","YX",0,0,1200,600);
+ TCanvas* c2 = new TCanvas("ZX","ZX",0,0,1200,600);
+ c1->Divide(4,4);
+ c2->Divide(4,4);
+ for(int i=0; i<16; i++){
+ c1->cd(i+1);
+ command ="PosY:PosX>>hYX_"+to_string(i)+"(100,0,300,400,-200,200)";
+ chain->Draw(command.c_str(),"","",1,i);
+ //histname="hYX_"+to_string(i);
+ //h[i]= (TH1F*) gDirectory->FindObjectAny(histname.c_str());
+ //g[i] = new TGraph(h[i]);
+ //g[i]->SetMarkerStyle(3) ;
+ //g[i]->Draw("AP") ;
+
+ c2->cd(i+1);
+ command ="PosZ/10.:PosX>>hZX_"+to_string(i)+"(100,0,300,120,-10,110)";
+ chain->Draw(command.c_str(),"","*",1,i);
+ //histname="hWall2_"+to_string(i);
+ //h[i-1]= (TH1F*) gDirectory->FindObjectAny(histname.c_str());
+ //h[i-1]->SetMarkerColor(kRed) ;
+
+ //gPad->WaitPrimitive();
+ }
+}
diff --git a/Projects/NebulaPlus/macro/TdiffCal.C b/Projects/NebulaPlus/macro/TdiffCal.C
new file mode 100644
index 0000000000000000000000000000000000000000..b495a48e01ffd753c64ae10bce3759d8453802bd
--- /dev/null
+++ b/Projects/NebulaPlus/macro/TdiffCal.C
@@ -0,0 +1,45 @@
+TFile* file=NULL ;
+
+////////////////////////////////////////////////////////////////////////////////
+void LoadFile(){
+ file = new TFile("../root/analysis/hTdiff_ID_run18.root");
+}
+
+////////////////////////////////////////////////////////////////////////////////
+void TdiffCal(bool savehist=true){
+ LoadFile();
+ TH2F* hQu_ID = (TH2F*) gDirectory->FindObjectAny("hTdiff_ID");
+
+ TCanvas* c1 = new TCanvas("Raw","Raw",0,0,1200,600);
+ c1->cd();
+ //c1->SetLogy();
+
+ TH1F* hQu[90];
+ Double_t mean=0;
+ string histname;
+ bool found = false;
+
+ //RANGE FOR THE MUON PEAK SEARCH
+ double range_inf= -30;
+ double range_sup= 30;
+
+ ofstream peaks_file;
+ peaks_file.open("tmp.txt");
+ for(int i=1; i<93; i++){
+
+ histname = "hQu_"+to_string(i);
+ cout<<i<<"histname\t"<<histname<<endl;
+ hQu_ID->ProjectionY(histname.c_str(),i+1,i+1);
+ hQu[i-1] = (TH1F*) gDirectory->FindObjectAny(histname.c_str());
+ hQu[i-1]->Draw("");
+ hQu[i-1]->SetTitle(histname.c_str());
+ hQu[i-1]->GetXaxis()->SetRangeUser(range_inf,range_sup);
+ mean = hQu[i-1]->GetMean(1);
+ //if(!found) cout<<"BarID:"<<i<<"UP, problem in finding peaks"<<endl;
+ //gPad->WaitPrimitive();
+ peaks_file<<"NebulaPlus_Tdiff_Offset_"<<i<<"\t"<<mean<<endl;
+ cout<<i<<"mean:\t"<<mean<<