From 36b230d23a048498a86527c40f5d864641eacb4f Mon Sep 17 00:00:00 2001 From: James Szalkie Date: Mon, 13 Apr 2026 16:04:00 -0400 Subject: [PATCH] energy loss table plot --- Armory/anasenMS.cpp | 7 +++++++ 1 file changed, 7 insertions(+) diff --git a/Armory/anasenMS.cpp b/Armory/anasenMS.cpp index 190de60..c6176f9 100644 --- a/Armory/anasenMS.cpp +++ b/Armory/anasenMS.cpp @@ -11,6 +11,7 @@ #include "TApplication.h" // ROOT app loop #include "ClassTransfer.h" // Reaction kinematics and MC event generation #include "ClassAnasen.h" // ANASEN detector model classes (SX3, PW, etc.) +#include "TCanvas.h" // ROOT canvas for drawing //======== Generate light particle based on reaction // calculate real and reconstructed tracks and Q-value uncertainty @@ -35,6 +36,12 @@ int main(int argc, char **argv){ TGraph* elossLight = LoadELoss("../ELoss/Eloss_alpha"); // for light particle (alpha) TGraph* elossHeavy = LoadELoss("../ELoss/Eloss_p"); // for heavy particle (proton) double density = 0.0000876; // example for aluminum target, adjust as needed + auto c1 = new TCanvas("c1", "Graph Example", 800, 600); + auto g = elossLight; + g->SetTitle("Energy Loss Table;Kinetic Energy (MeV);dE/dx (MeV/(mg/cm^{2}))"); + g->Draw("ALP"); + g->SetLineColor(kRed); + c1->Print("eloss_light.png"); // Reaction setup: projectile + target configuration, energy, and product IDs TransferReaction transfer;