diff --git a/Armory/anasenMS.cpp b/Armory/anasenMS.cpp
index 93b5fe4..96a344b 100644
--- a/Armory/anasenMS.cpp
+++ b/Armory/anasenMS.cpp
@@ -96,15 +96,15 @@ int main(int argc, char **argv){
   if( argc >= 2 ) numEvent = atoi(argv[1]);
   TransferReaction transfer;
 
-  transfer.SetA(18, 10, 0);            // 18Ne projectile
+  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_Ne-18.dat"); 
   transfer.Seta(4, 2);                 // 4He target
   transfer.Setb(1, 1);                 // outgoing proton from the primary transfer
-  transfer.SetB(21, 11);               // 21Na* heavy product
-  const double beamA = 18.0;                   // mass number of 18Ne beam
+  transfer.SetB(30, 14);               // 21Na* heavy product
+  const double beamA = 27.0;                   // mass number of 27Al beam
 
-  bool enableSequentialDecay = true; // turning to false to disable sequential decay for now, can be set to true to enable
+  bool enableSequentialDecay = false; // turning to false to disable sequential decay for now, can be set to true to enable
   const int decayDaughterA = 20;
   const int decayDaughterZ = 10;
   const int decayEjectA = 1;
@@ -112,7 +112,7 @@ int main(int argc, char **argv){
 
   // Excited state lists (projectile and heavy-product excitation states)
   std::vector<float> ExAList = {0};              // 18Ne projectile excitations in MeV
-  std::vector<float> ExList = {2.797, 2.829, 3.544, 3.862, 4.294};             // 21Na* excitation in MeV
+  std::vector<float> ExList = {0};             // 21Na* excitation in MeV
 
   // define vertex position uniform distribution ranges (mm)
   double vertexXRange[2] = {   -5, 5}; // mm - 5, 5