Adding reaction file

Inputs/DetectorConfiguration/hira.detector

@@ -7,10 +7,10 @@ GeneralTarget
 % Radius in mm
 % Temperature in K, Pressure in bar
 Target
-	THICKNESS= 25
+	THICKNESS= 6.26
   ANGLE= 0
 	RADIUS=	10
-	MATERIAL= CH2
+	MATERIAL= CD2
 	X= 0
 	Y= 0
 	Z= 0

Inputs/EventGenerator/34Ar_pd.reaction

@@ -21,7 +21,7 @@ TwoBodyReaction
 	Heavy= 33Ar
 	ExcitationEnergyHeavy= 0.0
 	ExcitationEnergyLight= 0.0
-  CrossSectionPath= CS_2.txt CS
+  CrossSectionPath= flat.txt CS
 	ShootLight= 1
 	ShootHeavy= 0
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

Inputs/EventGenerator/46Ar_pd.reaction

@@ -21,7 +21,7 @@ TwoBodyReaction
 	Heavy= 45Ar
 	ExcitationEnergyHeavy= 0.542
 	ExcitationEnergyLight= 0.0
-  CrossSectionPath= CS_2.txt CS
+  CrossSectionPath= flat.txt CS
 	ShootLight= 1
 	ShootHeavy= 0
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

NPAnalysis/Hira/Analysis.cxx

@@ -34,8 +34,8 @@ int main(int argc, char** argv)
     //	Instantiate some Reaction
     NPL::Reaction*  TransfertReaction = new Reaction								;
     //TransfertReaction	->	ReadConfigurationFile("34Ar_pd.reaction")	;
-    TransfertReaction	->	ReadConfigurationFile("46Ar_pd.reaction")	;
-    //TransfertReaction	->	ReadConfigurationFile("11Be_d3He.reaction")	;
+    //TransfertReaction	->	ReadConfigurationFile("46Ar_pd.reaction")	;
+    TransfertReaction	->	ReadConfigurationFile("11Be_d3He.reaction")	;
     
     //Get Detector pointer :
     THiraPhysics* Hira 	  = (THiraPhysics*) 			myDetector -> GetDetector("HIRAArray")		;

NPAnalysis/Hira/Analysis.h

@@ -124,8 +124,8 @@ using namespace NPL ;
 namespace ENERGYLOSS
 	{
 	
-            EnergyLoss EL_deuteron_CH2 =  EnergyLoss("deuteron_CH2.G4table","G4Table",100 );
-         //EnergyLoss EL_deuteron_CH2 =  EnergyLoss("He3_CH2.G4table","G4Table",100 );
+            //EnergyLoss EL_deuteron_CH2 =  EnergyLoss("deuteron_CH2.G4table","G4Table",100 );
+         EnergyLoss EL_deuteron_CH2 =  EnergyLoss("He3_CD2.G4table","G4Table",100 );
 		//	Declare your Energy loss here	:
 	/*		EnergyLoss DeuerontTarget = EnergyLoss 	(	"CD2.txt" 	,
 														100 		,

NPAnalysis/Hira/RunToTreat.txt

 TTreeName 
 	SimulatedTree
 RootFileName 
-	/Users/pierremorfouace/Physics/NPTool/nptool/Outputs/Simulation/46Ar_pd_gs.root
- 	/Users/pierremorfouace/Physics/NPTool/nptool/Outputs/Simulation/46Ar_pd_1st.root
+	%/Users/pierremorfouace/Physics/NPTool/nptool/Outputs/Simulation/46Ar_pd_gs.root
+ 	%/Users/pierremorfouace/Physics/NPTool/nptool/Outputs/Simulation/46Ar_pd_1st.root
 
-	%/Users/pierremorfouace/Physics/NPTool/nptool/Outputs/Simulation/34Ar_pd_gs.root
+	/Users/pierremorfouace/Physics/NPTool/nptool/Outputs/Simulation/11Be_d3He.root
 
 	
 	

NPAnalysis/Hira/Show.C

@@ -41,14 +41,14 @@ void Show()
     
     
     mainC1->cd(1);
-    chain->Draw("E_ThickSi:E_CsI>>h1(1000,0,200,1000,0,25)","","colz");
+    chain->Draw("E_ThickSi:E_CsI>>h1(1000,0,600,1000,0,30)","","colz");
     TH2F* h1 = (TH2F*)gDirectory->FindObjectAny("h1");
     h1->GetXaxis()->SetTitle("E_{CsI} (MeV)");
     h1->GetYaxis()->SetTitle("E_{Si} (MeV)");
     h1->SetTitle("");
 
     mainC1->cd(2);
-    chain->Draw("E_ThinSi:E_ThickSi>>h12(1000,0,25,1000,0,5)","","colz");
+    chain->Draw("E_ThinSi:E_ThickSi>>h12(1000,0,30,1000,0,10)","","colz");
     TH2F* h12 = (TH2F*)gDirectory->FindObjectAny("h12");
     h12->GetXaxis()->SetTitle("E_{Si} (MeV)");
     h12->GetYaxis()->SetTitle("#Delta E (MeV)");
@@ -56,16 +56,16 @@ void Show()
 
     
     mainC2->cd();
-    chain->Draw("ELab:ThetaLab>>h2(1000,0,50,1000,0,200)","","colz");
+    chain->Draw("ELab:ThetaLab>>h2(1000,0,50,1000,0,50)","","colz");
     h2->SetMinimum(1);
     TH2F* h2 = (TH2F*)gDirectory->FindObjectAny("h2");
     h2->GetXaxis()->SetTitle("#theta_{lab} (deg)");
     h2->GetYaxis()->SetTitle("E (MeV)");
     h2->SetTitle("");
     
-    //NPL::Reaction *r = new NPL::Reaction("11Be(d,3He)10Li@770");
+    NPL::Reaction *r = new NPL::Reaction("11Be(d,3He)10Li@770");
     //NPL::Reaction *r = new NPL::Reaction("34Ar(p,d)33Ar@2380");
-    NPL::Reaction *r = new NPL::Reaction("46Ar(p,d)45Ar@3220");
+    //NPL::Reaction *r = new NPL::Reaction("46Ar(p,d)45Ar@3220");
     kin = r->GetKinematicLine3();
     
     kin->SetLineColor(2);
@@ -122,8 +122,8 @@ void InitChain()
     chain = new TChain("AnalysedTree","");
     
     //chain->Add("/Users/pierremorfouace/Physics/NPTool/nptool/Outputs/Analysis/34Ar_pd.root");
-    chain->Add("/Users/pierremorfouace/Physics/NPTool/nptool/Outputs/Analysis/46Ar_pd.root");
-    //chain->Add("/Users/pierremorfouace/Physics/NPTool/nptool/Outputs/Analysis/11Be_d3He_gs.root");
+    //chain->Add("/Users/pierremorfouace/Physics/NPTool/nptool/Outputs/Analysis/46Ar_pd.root");
+    chain->Add("/Users/pierremorfouace/Physics/NPTool/nptool/Outputs/Analysis/11Be_d3He_1mg.root");
     
     return;
 }

NPAnalysis/newHira/CMakeLists.txt

-cmake_minimum_required (VERSION 2.8) 
-#Finding NPTool
-set(NPTOOL "$ENV{NPTOOL}")
-set(NPLIB "${NPTOOL}/NPLib")
-set(NPTOOL_INCLUDE_DIR "${NPLIB}/include")
-set(NPTOOL_LIB_DIR "${NPLIB}/lib")
-  
-include("${NPLIB}/FindROOT.cmake")  
-
-project (NPAnalysis)
-set(CMAKE_BUILD_TYPE Release)  
-# Add root to the link and include directories
-include_directories( ${ROOT_INCLUDE_DIR})
-link_directories( ${ROOT_LIBRARY_DIR})
-include_directories( ${NPTOOL_INCLUDE_DIR})
-link_directories( ${NPTOOL_LIB_DIR})
-
-# Get the compilator flag from root to assure consistancy
-EXEC_PROGRAM(${ROOT_CONFIG_EXECUTABLE}
-     ARGS "--cflags"
-     OUTPUT_VARIABLE root_cflags )
-
-set(CMAKE_CXX_FLAGS "${CMAKE_CXX_FLAGS} ${root_cflags}")
-
-# If the compiler is Clang, silence the unrecognised flags
-if(${CMAKE_CXX_COMPILER_ID} MATCHES ".*Clang.*")
-  set(CMAKE_CXX_FLAGS "${CMAKE_CXX_FLAGS} -Qunused-arguments -undefined dynamic_lookup")
-endif()
-
-add_library(NPAnalysis SHARED Analysis.cxx)
-target_link_libraries(NPAnalysis ${ROOT_LIBRARIES} -L${NPLIB}/lib -lNPCore -lNPPhysics)