# ***************************************************************************
#			    NSL RETINA MODEL 
#
#	Retina model description:
#	------------------------ 
#	Quantitative modeling of Responses of Anuran retina: Stimulus Shape and 
#	Size Dependency. In this simulations we test the model's ability to
#	reproduce quantitatively tabulated data on the dependency on stimulus
#	shape and size (Ewert 1976).
#	Input to the model is Light on the receptors (40X40 to simulate the 
#	receptieve field of a single ganglion cell of each type).
#	The model also simulates "simple" Horizontal and Bipolar cells.
# 	The output of the model the temporal firing rate of a  Ganglion cell of 
#   each type (R2, R3, R4). So the graphs given by NSL show the firing rate 
#	of a ganglion cell vs. the time.
#    
#	Notes:
#	-----
# 	No Trailing edge effect since Ewert computed his data with only
#	the response to the leading edge.
# 	No Leaky Integrators for the Ganglion Cells.
#
# 	Goal: Match Ewert's quantitative data on the Toad's retinal ganglion 
# 	cells.
#  	Status: goal achieved, paper submitted to Vision Research.
#
# 	File : "retina.nsl" 
#	Contents: Parameter values and Graphics set up. (see "retina.mod")
# 	Last update: Feb 1st, 2000
#
# ***************************************************************************

nsl set system.model RetinaModel

nsl set system.simDelta 0.06579  	;#  **** Simulation Time Step = 65.79 msec.
nsl set system.simEndTime 7.0	;#  Total simulation time = 7 sec

;# ______________________________ PARAMETERS _________________________________

nsl set RetinaModel.retina.amacrineOff.axh_tm  0.3  ;# AMACRINES time constant
 													;# 0.05 for velocity experiments.
nsl set RetinaModel.retina.amacrineOn.axd_tm 0.3  

nsl set RetinaModel.retina.horizontal.hor_tm  0.1   ;# HORIZONTALS  time constant

nsl set RetinaModel.retina.horizontal.H_level 0     ;# Uniform horizontal cell potential
      												;# 0 if the background is bright, 1 if dark.
nsl set RetinaModel.retina.r3.p  0.5	;#0.3		;# Effect of trailing edge on R3

nsl set RetinaModel.retina.r2.trailing 0			;# Effect of trailing edge on R2.

# Constants of integration  ** NOT USED **
nsl set RetinaModel.retina.r2.tm 0.05
nsl set RetinaModel.retina.r3.tm 0.05
nsl set RetinaModel.retina.r4.tm 0.1

# ************************** MASKS ***********************

nsl set RetinaModel.retina.r2.erf_dia 4.0
nsl set RetinaModel.retina.r2.irf_dia 19.5
nsl set RetinaModel.retina.r2.erf_sig 2.4
nsl set RetinaModel.retina.r2.irf_sig 4.0
nsl set RetinaModel.retina.r2.erf_wgt 1.0
nsl set RetinaModel.retina.r2.irf_wgt 2.3		;# Teeters was 2.1	         .	

nsl set RetinaModel.retina.r3.erf_dia 8.0
nsl set RetinaModel.retina.r3.irf_dia 19.5
nsl set RetinaModel.retina.r3.erf_sig 2.0		;# Teeters was 5		
nsl set RetinaModel.retina.r3.irf_sig 10.0
nsl set RetinaModel.retina.r3.erf_wgt 1.15		;# Teeters was 1
nsl set RetinaModel.retina.r3.irf_wgt 2.38		;# Teeters was 2                     			    

nsl set RetinaModel.retina.r4.erf_dia 15.5                       
nsl set RetinaModel.retina.r4.erf_sig 3.5
nsl set RetinaModel.retina.r4.erf_wgt 1.0

nsl set RetinaModel.retina.r2.thresh 0.001		;# Thresholds for the Ganglion Cells.
nsl set RetinaModel.retina.r3.thresh 0.001
nsl set RetinaModel.retina.r4.thresh 0.001

nsl set RetinaModel.retina.r2.k 43.8 			;# Scaling Factors for the ganglion cells.				
nsl set RetinaModel.retina.r3.k 44.0
nsl set RetinaModel.retina.r4.k 37.5

nsl set RetinaModel.retina.r2.pbh_erf 0.3
nsl set RetinaModel.retina.r2.pbh_irf 0

nsl load retina_graph

# nsl set input_layer RetinaModel.retina.r
# distance between the center of two squares (dx,dy) should be (0.5,0.5). 
# We set them to (2,2) for visual purposes. All stimuli values should multiplied by 4 accordingly
nsl set RetinaModel.retina.visin.in.dx 2     ;# 
nsl set RetinaModel.retina.visin.in.dy 2     ;# 
nsl set RetinaModel.retina.visin.in.xz 0	 ;#       ******* initial Position
nsl set RetinaModel.retina.visin.in.yz 20	 ;#18  

# The model is used for different size and shape stimuli. In particular we 
# tried to reproduce Ewert's data on presentation of Squares, Worms and 
# Antiworms from 2 to 32 visual degrees.
#
# Load one of the following:
#
# 1. nsl load retina_worm.nsl
# 2. nsl load retina_antiworm.nsl
# 3. nsl load retina_square.nsl