OPTION, ECHO = FALSE;
OPTION, PSDUMPFREQ  = 10;
OPTION, SPTDUMPFREQ = 10;
OPTION, PSDUMPEACHTURN = FALSE;
OPTION, REPARTFREQ = 1000;
OPTION, CZERO = FALSE;
OPTION, TELL = TRUE;
Option, VERSION=10900;

TITLE, string="Elektrischer Deflektor Test";

//---------------------------------------------------------------------------------------------//                          
//                                      Physical Parameters                                    //                          
//---------------------------------------------------------------------------------------------//                          
REAL Ebeam         = 1*0.001;                      // injection energy of the beam                                    
REAL gamma         = (Ebeam+PMASS)/PMASS;           // relativistic gamma
REAL beta          = sqrt(1-(1/gamma^2));           // relativistic beta
REAL gammabeta     = gamma*beta;                    // self explanitory       
REAL P0            = PMASS*gammabeta;               // Initial momentum of the particle                                                        
VALUE,{gamma,beta,gammabeta,Ebeam,CLIGHT};         // Prints values in output to check                                
                                                   // initial conditions 

//---------------------------------------------------------------------------------------------//                          
//                                      Elements on the tranfer line                                    //                          
//---------------------------------------------------------------------------------------------// 

//Drift 1

D1: DRIFT, ELEMEDGE = 0.0, L = 0.5;

// Elektrostatischer Deflektor 1

EDEF1: SBEND, ELEMEDGE = 0.5, ANGLE = 0.68, E1 = 0. ,E2 = 0. , FMAPFN= "eBend.fmap", DESIGNENERGY = Ebeam;

// Example Bend

Bend: RBend, ANGLE = 30.0 * Pi / 180.0,
         FMAPFN = "1DPROFILE1-DEFAULT",
         ELEMEDGE = 0.5,
         DESIGNENERGY = 10.0,
         L = 0.5,
         GAP = 0.02;

// Drift 2
D2: DRIFT, ELEMEDGE = 1.0 , L = 0.5;

L1: LINE = (D1, Bend ,D2);
//L1: LINE = (D1, EDEF1, D2);

//---------------------------------------------------------------------------------------------//                          
//                                      Characteristics of the beam                                    //                          
//---------------------------------------------------------------------------------------------// 

Beam1: BEAM, PARTICLE = PROTON, CHARGE = 1.0, PC = P0, NPART = 1;

Dist1: DISTRIBUTION, TYPE = FROMFILE, FNAME = "sp.dist", INPUTMOUNITS = EV;

//---------------------------------------------------------------------------------------------//                          
//               Field solver (its off since we only have one particle)                                    //                          
//---------------------------------------------------------------------------------------------//


Fs1:FIELDSOLVER, FSTYPE=NONE, MX=32, MY=32, MT=32,
                 PARFFTX=true, PARFFTY=true, PARFFTT=false,
                 BCFFTX=open, BCFFTY=open, BCFFTT=open,
                 BBOXINCR=1, GREENSF=STANDARD;

//---------------------------------------------------------------------------------------------//                          
//                                      Tracking and Computation                               //                          
//---------------------------------------------------------------------------------------------//  

TRACK, LINE = L1, BEAM = Beam1, MAXSTEPS = 10000, DT = 1.0e-15;
RUN, METHOD = "PARALLEL-T", BEAM = Beam1, FIELDSOLVER = Fs1, 
     DISTRIBUTION = Dist1;
ENDTRACK;
STOP;