Distribution-Gauss-1.in 1.65 KB
Newer Older
gsell's avatar
gsell committed
1 2 3
OPTION, ECHO=FALSE; //ECHO Controls printing of an echo of input lines on the standard error file.
OPTION, INFO=FALSE; //INFO If this option is turned off, OPAL suppresses all information messages.
OPTION, STATDUMPFREQ=1;
4
Option, VERSION=10900;
gsell's avatar
gsell committed
5 6 7

TITLE, STRING="Gauss Distribution Protons";

8 9
REAL Edes   = 0.500;    // 0.5 GeV
REAL QB     = 0.77e-12;   // charge in Cb
10
REAL BF     = 1300.0e6;   // frequency in Mhz
11
REAL BC     = QB*BF;    // bunch current
gsell's avatar
gsell committed
12

13 14 15 16
REAL gamma=(Edes+PMASS)/PMASS;
REAL beta=sqrt(1-(1/gamma^2));
REAL gambet=gamma*beta;
REAL P0 = gamma*beta*PMASS;
gsell's avatar
gsell committed
17 18 19 20 21 22


Q1: QUADRUPOLE, L=1.0, ELEMEDGE=0.0, K1=0.0;

lattice: LINE = (Q1);

23
Fs1:FIELDSOLVER, FSTYPE = FFT, MX = 4, MY = 4, MT = 4,
gsell's avatar
gsell committed
24 25 26 27 28
PARFFTX = TRUE, PARFFTY = TRUE, PARFFTT = TRUE,
BCFFTX  = OPEN, BCFFTY  = OPEN, BCFFTT  = OPEN,
BBOXINCR = 1, GREENSF = INTEGRATED;


29
//Dist2:DISTRIBUTION, DISTRIBUTION = GAUSS,
30 31
//SIGMAX = 0.00001, SIGMAPX=0.00001,
//SIGMAY = 0.00001, SIGMAPY=0.00001,
32 33
//SIGMAT = 1e-8 , PT=Edes*1e9;

34 35 36
REAL sigpx = sqrt((0.00001*1.0e-9/PMASS + 1.0)^2 - 1.0);
REAL sigpy = sqrt((0.00001*1.0e-9/PMASS + 1.0)^2 - 1.0);
REAL pzave = sqrt((Edes/PMASS + 1.0)^2 - 1.0);
37

38
Dist2:DISTRIBUTION, TYPE = GAUSS,
39
		    SIGMAX = 0.00001,
40
		    SIGMAPX = 0.00001,
41
		    SIGMAY = 0.00001,
42
		    SIGMAPY = 0.00001,
43
		    SIGMAT = 1e-8 ,
adelmann's avatar
update  
adelmann committed
44
		    CORRX = -0.0,
45
		    CORRY = 0.5,
46
		    INPUTMOUNITS = EVOVERC;
47

gsell's avatar
gsell committed
48 49

// pc:   momentum per particle in GeV/c
50
beam1: BEAM, PARTICLE = PROTON, PC = P0, NPART = 1000000, BFREQ = BF * 1e-6, BCURRENT = BC;
gsell's avatar
gsell committed
51 52 53

SELECT, LINE=lattice;

54
TRACK, LINE=lattice, BEAM=beam1, MAXSTEPS=1, DT=1e-11, ZSTOP=0.25;
gsell's avatar
gsell committed
55 56
 RUN, METHOD = "PARALLEL-T", BEAM = beam1, FIELDSOLVER = Fs1, DISTRIBUTION = Dist2;
ENDTRACK;
57
QUIT;