{
SolenoidsTo15cm - September 26th, 2003
The solenoid decay channel, optimised in high generality.
Identical to SolenoidsOnly except last solenoid fixed in radius.
}

{RodSource File X Y Z RodAngle}
TantalumRod 'pi_plus_2G2eV' 0 0 [0,0.45] [0,0.5]

#z=0; ok=1;
#for (n=1;ok;n++)
#{
	{Solenoid Length Radius Field}
#	if (n==1) #S1 [0.2,0.45] 0.1 [0,20];#;
#	else {S[n]# [0.2,0.6] [0.1,0.4] #if (n<=4) #[-5,5];#; else #[-4,4];#;}
#	z+=S[n].Length;
	{Drift Length}
	#D[n]# [0.5,1]
#	z+=D[n].Length;
#	if (z>=30) ok=0;
#}
#MaxSol=n-1;
{Solenoid Length Radius Field}
Sfinal [0.2,0.6] 0.15 [-4,4]

{Finish-aperture radius} END 0.5

DecayChannel: #for (n=1;n<=MaxSol;n++) {S[n]#,#D[n]#,#}#Sfinal;

TantalumRod,DecayChannel,END;