// example of tune calculation with closed orbit finder for MK 90 cyclotron.

OPTION, VERSION = 20000;

OPTION, ECHO= FALSE;
OPTION, PSDUMPFREQ= 24500000;
OPTION, SPTDUMPFREQ= 50;
OPTION, PSDUMPEACHTURN= FALSE;
OPTION, PSDUMPFRAME= GLOBAL;

// For Tune Calculation
OPTION, CLOTUNEONLY =true;

TITLE,STRING= "MK90 CYCL";

REAL Kr86MASS=86*0.93149432;//GeV
REAL Edes = 1.0192/1000; // GeV
REAL r0   = 38.3;         //mm
REAL pr0  = 0.003645; //m_0*c
REAL f1   = 37;
 
// calculate initial momentum in units beta gamma
REAL gamma  = (Edes+Kr86MASS)/Kr86MASS;
REAL beta   = sqrt(1-(1/gamma^2));
REAL gambet = gamma*beta;
REAL P0     = gamma*beta*Kr86MASS;

// Cyclotron definition
// FMLOWE and FMHIGHE needed for range of closed orbit finder
Ring: CYCLOTRON, TYPE=CARBONCYCL, CYHARMON=4, PHIINIT=108.8, PRINIT=pr0, RINIT=r0, SYMMETRY=4, RFFREQ=f1, FMAPFN="MainMagnetField.dat", FMLOWE =1.0192, FMHIGHE = 9.74*86;

// Beamline
L1:   LINE = (Ring);

//For Tune Calculation with OPAL
DistTO: DISTRIBUTION,  TYPE=GAUSSMATCHED, LINE=L1,
        NSTEPS=1440, DENERGY=0.001, MAXSTEPSCO = 100, NSECTORS=4, SECTOR=FALSE;
// Field solver
Fs1:FIELDSOLVER, FSTYPE=NONE, MX=64, MY=64, MT=64, PARFFTX=true, PARFFTY=true, PARFFTT=false,
                 BCFFTX=open, BCFFTY=open, BCFFTT=open;

Beam1: BEAM, PARTICLE=Kr86, PC=P0, NPART=1, BCURRENT=1.0E-6, CHARGE=12.0, BFREQ= f1;

SELECT, LINE=L1;

TRACK, LINE=L1, BEAM= Beam1, MAXSTEPS= 720, STEPSPERTURN= 720;

RUN, METHOD= "CYCLOTRON-T", BEAM= Beam1, FIELDSOLVER= Fs1, DISTRIBUTION= DistTO;

ENDTRACK;

STOP;