Newer
Older
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
#include "NPReaction.h"
#include <Rtypes.h>
#include <TCanvas.h>
#include <TString.h>
void PlotTheoreticalKinematics()
{
double beam10Be {28. * 10.013534};
double beam12Be {30. * 12.026922};
//simply compare reaction kinematics
NPL::Reaction r10Be(TString::Format("10Be(d,3He)9Li@%.f", beam10Be).Data());
NPL::Reaction r12Be (TString::Format("12Be(d,3He)11Li@%.f", beam12Be).Data());
//get graphs
auto* g9Li {r10Be.GetKinematicLine3()};
auto* g11Li {r12Be.GetKinematicLine3()};
//also Elab vs thetaCM
auto* gcm9Li {r10Be.GetELabVersusThetaCM()};
auto* gcm11Li {r12Be.GetELabVersusThetaCM()};
//theta3 vs theta4 in LAB
auto* gtheta9Li {r10Be.GetTheta3VsTheta4()};
auto* gtheta11Li {r12Be.GetTheta3VsTheta4()};
//plotting
auto* c1 {new TCanvas("c1", "Kinematic lines comparaison")};
c1->DivideSquare(4);
c1->cd(1);
g11Li->SetTitle(";#theta_{3}^{LAB} [degree];T_{3} [MeV]");
g9Li->SetLineWidth(2); g9Li->SetLineColor(kRed);
g11Li->SetLineWidth(2); g11Li->SetLineColor(kBlue);
g11Li->Draw("apl");
g9Li->Draw("pl same");
c1->cd(2);
gcm11Li->SetTitle(";#theta_{3}^{CM} [degree];T_{3} [MeV]");
gcm9Li->SetLineWidth(2); gcm9Li->SetLineColor(kRed);
gcm11Li->SetLineWidth(2); gcm11Li->SetLineColor(kBlue);
gcm11Li->Draw("apl");
gcm9Li->Draw("pl same");
c1->cd(3);
gtheta9Li->SetTitle(";#theta_{3}^{LAB} [degree];#theta_{4}^{LAB} [degree]");
gtheta9Li->SetLineWidth(2); gtheta9Li->SetLineColor(kRed);
gtheta11Li->SetLineWidth(2); gtheta11Li->SetLineColor(kBlue);
gtheta9Li->Draw("apl");
gtheta11Li->Draw("pl same");
}