ANASEN_analysis/results/compareE.C
Sudarsan Balakrishnan 3b99dc05d1 Bugfixes/Changelog:
- Changed Kinematics.h to use TLorentzVector - doesn't affect Ex calculations at all
	- 27Al, 17F analyses fleshed out to identify vertexz regions sensitive to a, d, p ejectiles
	- Better handling of dE using detector hyperboloids, dE4 introduced. Marginal improvement over dE3. This is a function inside PC_Stepladder_Correction.h which needs to be renamed.
	- Anode-wire #23 seems to have disproportionally large # of counts compared to the other wires.
	- Move towards wire-wise dE/E gates. Not yet fully implemented.
	- This push sees some Ex plots that appear to be from (a,p) on 27Al(a,p). 17F(a,p) is a lot less clear-cut.
	- Updated beam starting energies to improved guesses, for p+a, 17F and 27Al data.
2026-07-07 14:40:18 -04:00

33 lines
1.1 KiB
C

{
TGraph kinematics("a(p,p)a_kinematics_6.88MeV_p.txt","%*lf %*lf %lf %lf");
kinematics.Scale(4.0);
kinematics.SetLineColor(kBlue);
kinematics.SetLineWidth(2.);
//TFile fin("4He(pp)_candidate_kinematic_curve.root");
//TH2F *h2 = (TH2F*)(fin.Get("Ef_thetaf_AnodeQQQR_TC1_PC1_pidlow0_1"));
// TFile fin("out_zminus5_catima.root");
TFile fin("out.root");
//TH2F *h2 = (TH2F*)(fin.Get("Ef_thetaf_AnodeQQQR_TC1_PC1_pidlow0_1"));
// TFile fin("run15_sx3_qqq_alphas_kinematic.root");
TH2F *h2 = (TH2F*)(fin.Get("a(p,p)/ap_Ef_vs_theta_qqq_a1c2"));
h2->Draw("box");
h2->RebinY(8);
kinematics.Draw("L SAME");
gPad->Modified();
gPad->Update();
while(gPad->WaitPrimitive());
TGraph kinematicsp("a(p,p)a_kinematics_7MeV_p.txt","%lf %lf %*lf %*lf");
kinematicsp.SetLineWidth(2.);
TH2F *h2p = (TH2F*)(fin.Get("a(p,p)/ap_Ef_vs_theta_sx3_a1c2"));
h2p->RebinY(8);
h2p->SetLineColor(kBlack);
//h2p->SetLineColor(kRed);
h2p->Draw("box same");
kinematicsp.Draw("L SAME");
gPad->SaveAs("c1.png");
}