#*******************************************************************
# 15 March 2010
# File for checking TRacy is running well
#
# input script to set the file and bool flag
#
#  comments:
# 1) Blank line and line starting with "#" (comment) are ignored
# 2) Mandatory parameter: lat_file,
#     lat file name is without the extension of .lat
# 3) All the other parameters are optional, and all the default 
#      values for flags are false
#
#  Test verison for the development of Tracy 3, always the nearest
#        input file verison.
#
#                     09/2011    Jianfeng Zhang 
#
#******************************************************************

#*************************
#  mandatory setting
**************************
# relative directory of the files in the script
 in_dir ./
#in_dir  /home/zhang/codes/TracyIII/lattice/example/

# lattice file must be wihout .lat extension
  lat_file   example_lattice  

# H-corrector name in lattice file
  h_corr        HCM    
# V-corrector name in lattice file
  v_corr        VCM    
# girder start symbol
  gs            GS     
# girder end symbol
  ge            GE     
# BPM name in lattice file
  bpm_name      BPM    
# skew quadrupole
  qt            QT


#*******************************
#  option setting 
#*******************************
# GLOBAL PARAMETER PART
# Activate quadrupole Fringefield
 QuadFringeOnFlag 


# VACCUUM CHAMBER PART
# read vacuum chamber from chamber_file
  ReadChamberFlag    example_chamber.dat    

# set RF voltage, if false, then use the definition value
# of RF cavity in the lattice
# name, flag, RF voltage [V]
  RFvoltageFlag   3000000

  
#####################
# Read the field error, and replacement the corresponding components
# of the element.
#  Generic features
####################

ReadfefileFlag   example_dip.fe





# fit tunes
  FitTuneFlag     q7   q9  18.202 10.317
# fit chromaticities
  FitChromFlag   sx9 sx10  2.0    2.0

# NOT finished
#PrintGirderFlag testgirder.out

######################
# alignment error file, after read the misalignment, the orbit is corrected
#####################
  ReadaefileFlag    example.ae
#####################
#  orbit correction, if ReadaefieldFlag is on, then must set the parameters to do orbit correction
#####################
# file to specify which correctors are used for H orbit correction  
  hcorr_file  example_hcorr_56nom.state
# file to specify which correctors are used for V orbit correction  
  vcorr_file  example_vcorr_56nom.state

 #number of statistics (seeds) --> N flat-files
   n_stat        2    
 #number of times to scale the errors
  n_scale       1     
 #number of iterations for orbit correction, if 0, no orbit correction 
  n_orbit       0     
#number of singular values in H-plane, must be not larger than the number of correctors used for orbit correction
nwh           60    
#number of singular values in V-plane, must be not larger  than the number of correctors used for orbit correction
nwv           60    



#######################################
#  print the twiss paramters to a file
#######################################
  PrintTwissFlag     twisstest.out

########################
# print close orbit to a file
##########################
  PrintCODFlag codtest.out

#print  coordinates tracking around COD at each element
PrintTrackFlag  tracktest.out  0.001 0.000  0.0  0.0  0.0  0.0   50


# tune obtained by tracking
  TuneTracFlag   
# chromaticity obtained by tracking
  ChromTracFlag   

###
### FITTING FACTORY
###
# fit tune
# name, flag, tuned quadrople,tuned quadrople,targetnudx    targetnudz
  FitTuneFlag     q7 q9   18.202 10.317

# fit chromaticity
# name,flag, tuned sextupole, tuned sextupole, targetksix    targetksiz
  FitChromFlag   sx9 sx10  2.0    2.6


# TUNESHIFT WITH AMPLITUDE
# tune shift with amplitude and energy
#    AmplitudeTuneShiftFlag flag Nbx Nby TurnNumber xmax(m) ymax(m) delta
#    Nbx and Nby Number of amplitude values in H et V planes
#    TurnNumber number of turns for tracking  
#    xmax ymax maximum amplitude value (sqrt law)
#    delta energy offset
 AmplitudeTuneShiftFlag nudxtest.out nudztest.out   50 30 516 0.035 0.02 0.0

# EnergyTuneShiftFlag  Nb TurnNumber emax
#    Nb point number
#    TurnNumber number of turn for tracking  
#    emax maximum energy for tracking
 EnergyTuneShiftFlag nudptest.out  31 1026 0.06 


# FMA calculate frequency map 
#    Nbx Nby number of point in h and v planes
#    TurnNumber Number of turns for tracking
#    xmax, ymax maximum amplitude
#    energyoffset 
#    diffusionflag to compute tunediffusion 
 FmapFlag  fmaptest.out 31 21 516 0.025 0.005 0.0 true

#  FMA dp calculate frequency map
#     Nbx:  number of point in h planes
#     Nbe: number of point for the energy
#     TurnNumber: Number of turns for tracking
#     xmax: maximum amplitude in the horizontal plane
#     emax: maximum  energyoffset
#        y: vertical amplitude 
#     diffusionflag to compute tunediffusion 
 FmapdpFlag  fmapdptest.out 101 121 1026 25e-3 0.06 0.3e-3 true

# Error coupling
# add random rotation error to the full quadrupole
#     random seed number, RMS value of the rotation angle 
#    of the quadrupole
  ErrorCouplingFlag   0  0.0007
# coupling flag
#    calcualte the coupling factor in the lattice 
  CouplingFlag  

# calculate momentum acceptance for Touscheck lifetime
#   Dim  istart istop eminp emaxp nstepp eminn emaxn nstepn nturns zinitial
#   Dim: 4D/6D tracking, default value is 6D
  MomentumAccFlag momacctest.out  6D 1 189 0.01 0.05 100 -0.01 -0.05 100 1026 0.0002

# calculate phase space
# parameters are:  
# 4D/6D dimension to track the phase space; 
# x[m],x',y[m],y',delta,c*tau, Number of turns, flag to turn on/off
# damping during the tracking
  PhaseSpaceFlag phasetest.out 6D 1e-6 0.0 1e-6 0.0 0.012 0.0 1000 false




###
### Touschek lifetime
###
# calculate Touschek lifetime
# name, flag
  #TouschekFlag Qb eps[X_] eps[Y_] n_turns sigma_delta sigma_s

# Intra Beam Scattering
  #IBSFlag  
# Track momentum acceptance and then get Touschek lifetime
  #TousTrackFlag