rtang/ANASEN_analysis
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

refining energy loss with Lise

Browse Source
This commit is contained in:
James Szalkie 2026-06-17 12:03:07 -04:00
parent bcaa540a95
commit 4db59641cd
1 changed files with 4 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
Armory/anasenMS.cpp
View File

@ -97,12 +97,12 @@ int main(int argc, char **argv){
if( argc >= 2 ) numEvent = atoi(argv[1]);
TransferReaction transfer;
//To set beam energy loss, use energy loss app, and create table with target isotope, set Initial beam energy as max energy
transfer.SetA(27, 13, 0); // 18Ne projectile
//transfer.SetIncidentEnergyAngle(0, 0, 0); // KEA in MeV/u, theta and phi in rad
TGraph* elossBeam = LoadELoss("../ELoss/E_vs_x_Al-27.dat");
TGraph* elossBeam = LoadELoss("../ELoss/HeLoss/E_vs_x_Al-27.dat");
transfer.Seta(4, 2); // 4He target
transfer.Setb(2, 1); // outgoing proton from the primary transfer
transfer.SetB(29, 14); // 21Na* heavy product
transfer.Setb(1, 1); // outgoing proton from the primary transfer
transfer.SetB(30, 14); // 21Na* heavy product
const double beamA = 27; // mass number of 27Al beam
bool enableSequentialDecay = false; // turning to false to disable sequential decay for now, can be set to true to enable