diff --git a/NPAnalysis/Example1/Analysis.cxx b/NPAnalysis/Example1/Analysis.cxx index 94bb8bba9c7f50bb2d79f32f250ddbf68c3185a3..1e3b646b0435d7e39cee7d04ea1cffb449b8311f 100755 --- a/NPAnalysis/Example1/Analysis.cxx +++ b/NPAnalysis/Example1/Analysis.cxx @@ -109,7 +109,8 @@ int main(int argc, char** argv) int SiNumber = SSSD->DetectorNumber[countSSSD]; /************************************************/ - if(TelescopeNumber==SiNumber ){ + // Matching between Thin Si and MUST2, and Forward Telescope Only + if(TelescopeNumber==SiNumber && TelescopeNumber<5){ DetectorNumber = TelescopeNumber ; /************************************************/ // Part 1 : Impact Angle diff --git a/NPAnalysis/Example1/RunToTreat.txt b/NPAnalysis/Example1/RunToTreat.txt index 8ffc0450433a69fe4515234f458c3bdac6e3535f..82908222e590c8760e20a5e47347834912b6c19c 100755 --- a/NPAnalysis/Example1/RunToTreat.txt +++ b/NPAnalysis/Example1/RunToTreat.txt @@ -1,15 +1,5 @@ TTreeName SimulatedTree RootFileName - ../../Outputs/Simulation/He10_10.root -% ../../Outputs/Simulation/He10_1.root -% ../../Outputs/Simulation/He10_2.root -% ../../Outputs/Simulation/He10_3.root -% ../../Outputs/Simulation/He10_4.root -% ../../Outputs/Simulation/He10_5.root -% ../../Outputs/Simulation/He10_6.root -% ../../Outputs/Simulation/He10_7.root -% ../../Outputs/Simulation/He10_8.root -% ../../Outputs/Simulation/He10_9.root -% ../../Outputs/Simulation/He10_10.root + ../../Outputs/Simulation/Example1.root diff --git a/NPAnalysis/Example1/ShowResult.C b/NPAnalysis/Example1/ShowResult.C new file mode 100644 index 0000000000000000000000000000000000000000..6a137a24d4ab6e888decd97517ea35c683f70734 --- /dev/null +++ b/NPAnalysis/Example1/ShowResult.C @@ -0,0 +1,65 @@ +#include"NPReaction.h" + +TCutG* ETOF=NULL; +TCutG* EDE=NULL; +TChain* chain=NULL ; +TCanvas* c1 = NULL; + +//////////////////////////////////////////////////////////////////////////////// +void LoadCuts(){ +TFile* File_ETOF = new TFile("cuts/ETOF.root","READ"); +ETOF = (TCutG*) File_ETOF->FindObjectAny("ETOF"); + +TFile* File_EDE = new TFile("cuts/EDE.root","READ"); +EDE= (TCutG*) File_EDE->FindObjectAny("EDE"); +} + +//////////////////////////////////////////////////////////////////////////////// +void LoadChain(){ +chain = new TChain("ResultTree"); +chain->Add("../../Outputs/Analysis/Example1.root"); +} + +//////////////////////////////////////////////////////////////////////////////// +void ShowResult(){ +LoadChain(); +LoadCuts(); + +c1 = new TCanvas("Example1","Example1",0,0,600,600); +c1->Divide(2,2); + +// Light Particle ID // +// E-DE +c1->cd(1); +chain->Draw("SSSD.Energy:MUST2.Si_E>>hIDE(1000,0,35,1000,0,5)","MUST2.CsI_E<0 && MUST2.TelescopeNumber<5","colz"); +EDE->Draw("same"); + +// E-TOF +c1->cd(2); +chain->Draw("-MUST2.Si_T:SSSD.Energy+MUST2.Si_E>>hIDT(1000,0,35,1000,-15,0)","MUST2.CsI_E<0 && MUST2.TelescopeNumber<5","colz"); +ETOF->Draw("same"); + +// Kinematical Line // +c1->cd(3); +chain->Draw("ELab:ThetaLab>>hKine(1000,0,90,400,0,40)","MUST2.CsI_E<0 && MUST2.TelescopeNumber<5 && EDE && ETOF","colz"); + +NPL::Reaction r("11Li(d,3He)10He@553"); +TGraph* Kine = r.GetKinematicLine3(); +Kine->SetLineWidth(2); +Kine->SetLineColor(kOrange-3); +Kine->Draw("c"); + +// Excitation Energy // +c1->cd(4); +int bin=100; +double Emin = -10; +double Emax = 10; + +chain->Draw(Form("Ex>>hEx(%d,%f,%f)",bin,Emin,Emax),"MUST2.CsI_E<0 && MUST2.TelescopeNumber<5 && EDE && ETOF"); +TH1F* hEx = (TH1F*) gDirectory->FindObjectAny("hEx"); +hEx->GetYaxis()->SetTitle(Form("counts / %d keV",(int) (1000*(Emax-Emin)/bin))); +hEx->GetXaxis()->SetTitle("E_{10He}"); +hEx->SetFillStyle(1001); +hEx->SetLineColor(kAzure+7); +hEx->SetFillColor(kAzure+7); +} diff --git a/NPAnalysis/Example1/cuts/EDE.root b/NPAnalysis/Example1/cuts/EDE.root new file mode 100644 index 0000000000000000000000000000000000000000..1d18123c0647fc8d0dce0661e0b1cabf411b40a0 Binary files /dev/null and b/NPAnalysis/Example1/cuts/EDE.root differ diff --git a/NPAnalysis/Example1/cuts/ETOF.root b/NPAnalysis/Example1/cuts/ETOF.root new file mode 100644 index 0000000000000000000000000000000000000000..462bead6ba3db3a5e7e68426a892a94f0aa23d75 Binary files /dev/null and b/NPAnalysis/Example1/cuts/ETOF.root differ