diff --git a/Armory/anasenMS.cpp b/Armory/anasenMS.cpp index 2c09133..2be207b 100644 --- a/Armory/anasenMS.cpp +++ b/Armory/anasenMS.cpp @@ -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