#
#####################################
####	THIS IS A CONFIG FILE	 ####
#####################################
##
##  Blank lines or lines beginning with a "#" are ignored.
##
##  Double colons (::) are treated like a space.
##     This is sometimes useful for grouping (for example,
##     grouping strings like "reaction::sum::amplitudeName")
##
##  All non-comment lines must begin with one of the following keywords.
##
##  (note:  <word> means necessary 
##	    (word) means optional)
##
##  include	  <file>
##  define	  <word> (defn1) (defn2) (defn3) ...
##  fit 	  <fitname>
##  keyword	  <keyword> <min arguments> <max arguments>
##  reaction	  <reaction> <particle1> <particle2> (particle3) ...
##  data	  <reaction> <class> (arg1) (arg2) (arg3) ...
##  genmc	  <reaction> <class> (arg1) (arg2) (arg3) ...
##  accmc	  <reaction> <class> (arg1) (arg2) (arg3) ...
##  normintfile   <reaction> <file>
##  sum 	  <reaction> <sum> (sum2) (sum3) ...
##  amplitude	  <reaction> <sum> <amp> <class> (arg1) (arg2) ([par]) ... 
##  initialize    <reaction> <sum> <amp> <"events"/"polar"/"cartesian">
##		    <value1> <value2> ("fixed"/"real")
##  scale	  <reaction> <sum> <amp> <value or [parameter]>
##  constrain	  <reaction1> <sum1> <amp1> <reaction2> <sum2> <amp2> ...
##  permute	  <reaction> <sum> <amp> <index1> <index2> ...
##  parameter	  <par> <value> ("fixed"/"bounded"/"gaussian") 
##		    (lower/central) (upper/error)
##    DEPRECATED:
##  datafile	  <reaction> <file> (file2) (file3) ...
##  genmcfile	  <reaction> <file> (file2) (file3) ...
##  accmcfile	  <reaction> <file> (file2) (file3) ...
##
#####################################

# resonance parameters
parameter azeromass 0.980 fixed
parameter azerowidth 0.075 fixed
parameter atwomass 1.318 fixed
parameter atwowidth 0.111 fixed
parameter afourmass 1.995 fixed
parameter afourwidth 0.257 fixed
parameter pionemass 1.354 fixed
parameter pionewidth 0.330 fixed

fit etapi0

reaction EtaPi0 Beam Proton Eta Pi0

genmc EtaPi0 ROOTDataReader flat.root
accmc EtaPi0 ROOTDataReader flat.root
data  EtaPi0 ROOTDataReader sample1.root

normintfile EtaPi0 etapi0_ni.txt

# sum for helicity of N'=+1/2. Amplitudes should be duplicated for N'=-1/2.
sum EtaPi0 Negative
sum EtaPi0 Positive

# negative reflectivity
amplitude EtaPi0::Negative::S0- TwoPSAngles 0 0 -1
amplitude EtaPi0::Negative::S0- BreitWigner [azeromass] [azerowidth] 0 2 3

# positive reflectivity
amplitude EtaPi0::Positive::P1+ TwoPSAngles 1 1 1
amplitude EtaPi0::Positive::P1+ BreitWigner [pionemass] [pionewidth] 1 2 3

amplitude EtaPi0::Positive::D1+ TwoPSAngles 2 1 1
amplitude EtaPi0::Positive::D1+ BreitWigner [atwomass] [atwowidth] 2 2 3

amplitude EtaPi0::Positive::G1+ TwoPSAngles 4 1 1
amplitude EtaPi0::Positive::G1+ BreitWigner [afourmass] [afourwidth] 4 2 3

# initialize
initialize EtaPi0::Negative::S0- cartesian 10.0 0.0 real

initialize EtaPi0::Positive::P1+ cartesian 10.0 10.0
initialize EtaPi0::Positive::D1+ cartesian 10.0 0.0 real
initialize EtaPi0::Positive::G1+ cartesian 10.0 10.0