{ Linac900Ext4Xc2 = Linac900Ext6tc2 with 4 blocks and only accepting more accurate v4.45+ results May 14th, 2014 } !MINVER 4.45 {Parameter Value} t_max 3.5us n_decays 5 seed 0 seedpitch 1 {Cylinder-blockage Zsize radius Zhere Z rotZX AbsorptionLength} BackFace 0.10 0.40 #0.10/2# -0.15 0 0 {Annulus-blockage Zsize OuterRadius InnerRadius Z AbsorptionLength} BackAnnulus 0.20 0.40 0.30 -0.1 0 BackEnd: BackFace,BackAnnulus; {RodSource File X Y Z RodAngle} TantalumRod pi_plus_10GeV.csv 0 0 [0,0.45] [0,0.5] {Synchronous-particle t_start Energy Z} SyncMuon 5ns 180MeV #TantalumRod.Z# #Blocks=4; {Parameter Value} DecayCells [1,50]; #// 17-44 possible in original, PhaseRotC optimum was 23 #DecayCells=int(DecayCells); #for (n=1;n<=9 && n<=DecayCells;n++) #{ {Solenoid Length Radius Field} # if (n==1) #S1 [0.2,0.45] 0.1 [0,20];#; # else {S[n]# [0.2,0.9] [0.1,0.4] #if (n<=4) #[-5,5];#; else #[-4,4];#;} {Drift Length Bore} #D[n]# [0.5,1] 0.5 #} #for (bn=1;bn<=Blocks;bn++) #{ {Solenoid Length Radius Field} #SB[bn]# [0.2,0.9] [0.1,0.4] [-4,4] {Drift Length Bore} #DB[bn]# [0.5,1] 0.5 #} DecayChannel: #for (n=1;n<=9 && n<=DecayCells;n++) {S[n]#,#D[n] if (n=bn+1) bn++; # SB[bn] if (n=bn+1) bn++; LD#bn#a,RF#bn#,LD#bn#b,LS#bn if (n