Skip to content

Commit 2a3ce8c

Browse files
T-BrittonT-Britton
authored
Merge pull request #49 from JeffersonLab/DEV
Dev
2 parents 5aef677 + aebc5ac commit 2a3ce8c

File tree

14 files changed

+1238
-183
lines changed

14 files changed

+1238
-183
lines changed

‎MCwrapper/Generators/gen_amp/examples/gen_2k.cfg‎

Lines changed: 57 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,57 @@
1+
#
2+
#The above line is used to provide additional options to the generator within the wrapper. These options do NOT include output file names, run numbers, and numbewr of events which are controlled by the wrapper
3+
#####################################
4+
#### THIS IS A CONFIG FILE ####
5+
#####################################
6+
##
7+
## Blank lines or lines beginning with a "#" are ignored.
8+
##
9+
## Double colons (::) are treated like a space.
10+
## This is sometimes useful for grouping (for example,
11+
## grouping strings like "reaction::sum::amplitudeName")
12+
##
13+
## All non-comment lines must begin with one of the following keywords.
14+
##
15+
## (note: <word> means necessary
16+
## (word) means optional)
17+
##
18+
## include <file>
19+
## define <word> (defn1) (defn2) (defn3) ...
20+
## fit <fitname>
21+
## keyword <keyword> <min arguments> <max arguments>
22+
## reaction <reaction> <particle1> <particle2> (particle3) ...
23+
## data <reaction> <class> (arg1) (arg2) (arg3) ...
24+
## genmc <reaction> <class> (arg1) (arg2) (arg3) ...
25+
## accmc <reaction> <class> (arg1) (arg2) (arg3) ...
26+
## normintfile <reaction> <file>
27+
## sum <reaction> <sum> (sum2) (sum3) ...
28+
## amplitude <reaction> <sum> <amp> <class> (arg1) (arg2) ([par]) ...
29+
## initialize <reaction> <sum> <amp> <"events"/"polar"/"cartesian">
30+
## <value1> <value2> ("fixed"/"real")
31+
## scale <reaction> <sum> <amp> <value or [parameter]>
32+
## constrain <reaction1> <sum1> <amp1> <reaction2> <sum2> <amp2> ...
33+
## permute <reaction> <sum> <amp> <index1> <index2> ...
34+
## parameter <par> <value> ("fixed"/"bounded"/"gaussian")
35+
## (lower/central) (upper/error)
36+
## DEPRECATED:
37+
## datafile <reaction> <file> (file2) (file3) ...
38+
## genmcfile <reaction> <file> (file2) (file3) ...
39+
## accmcfile <reaction> <file> (file2) (file3) ...
40+
##
41+
#####################################
42+
43+
#Useful definition
44+
define phi 1.020 0.043
45+
46+
fit twok
47+
48+
reaction K+K- Beam Proton K+ K-
49+
50+
# consider just x polarized amplitudes
51+
sum K+K- xpol
52+
53+
# Currently not using any input parameters for TwoPiAngles in the generator
54+
amplitude K+K-::xpol::res BreitWigner phi 0 2 3
55+
56+
initialize K+K-::xpol::res cartesian 1.0 0.0
57+

‎MCwrapper/Generators/gen_amp/examples/gen_2pi_amp.cfg‎

Lines changed: 209 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,209 @@
1+
#
2+
#The above line is used to provide additional options to the generator within the wrapper. These options do NOT include output file names, run numbers, and numbewr of events which are controlled by the wrapper
3+
#####################################
4+
#### THIS IS A CONFIG FILE ####
5+
#####################################
6+
##
7+
## Blank lines or lines beginning with a "#" are ignored.
8+
##
9+
## Double colons (::) are treated like a space.
10+
## This is sometimes useful for grouping (for example,
11+
## grouping strings like "reaction::sum::amplitudeName")
12+
##
13+
## All non-comment lines must begin with one of the following keywords.
14+
##
15+
## (note: <word> means necessary
16+
## (word) means optional)
17+
##
18+
## include <file>
19+
## define <word> (defn1) (defn2) (defn3) ...
20+
## fit <fitname>
21+
## keyword <keyword> <min arguments> <max arguments>
22+
## reaction <reaction> <particle1> <particle2> (particle3) ...
23+
## data <reaction> <class> (arg1) (arg2) (arg3) ...
24+
## genmc <reaction> <class> (arg1) (arg2) (arg3) ...
25+
## accmc <reaction> <class> (arg1) (arg2) (arg3) ...
26+
## normintfile <reaction> <file>
27+
## sum <reaction> <sum> (sum2) (sum3) ...
28+
## amplitude <reaction> <sum> <amp> <class> (arg1) (arg2) ([par]) ...
29+
## initialize <reaction> <sum> <amp> <"events"/"polar"/"cartesian">
30+
## <value1> <value2> ("fixed"/"real")
31+
## scale <reaction> <sum> <amp> <value or [parameter]>
32+
## constrain <reaction1> <sum1> <amp1> <reaction2> <sum2> <amp2> ...
33+
## permute <reaction> <sum> <amp> <index1> <index2> ...
34+
## parameter <par> <value> ("fixed"/"bounded"/"gaussian")
35+
## (lower/central) (upper/error)
36+
## DEPRECATED:
37+
## datafile <reaction> <file> (file2) (file3) ...
38+
## genmcfile <reaction> <file> (file2) (file3) ...
39+
## accmcfile <reaction> <file> (file2) (file3) ...
40+
##
41+
#####################################
42+
43+
#Take useful definitions from gen_3pi.cfg
44+
define rho 0.775 0.146
45+
# some definitions for adjusting the beam polarization. In the equations beamX: Polarization>0, beamY: Polarization<0.
46+
# Phi is the azimuthal angle of the polarization vector in degrees (in the lab coordinates)
47+
define phipol TEMPPOLANGLE
48+
define polFrac TEMPPOLFRAC
49+
define beamX 0 polFrac
50+
define beamY 90 polFrac
51+
52+
# Uniform angles: flat=1; YLMs: flat=0;
53+
define flat 0
54+
55+
# fit twopi_amp
56+
57+
reaction Pi+Pi- Beam Proton Pi+ Pi-
58+
59+
normintfile Pi+Pi- twopi_ni.txt
60+
61+
# sum for helicity of N'=+1/2. Amplitudes should be duplicated for N'=-1/2.
62+
sum Pi+Pi- helplusN+
63+
sum Pi+Pi- helnegN+
64+
sum Pi+Pi- helplusN-
65+
sum Pi+Pi- helnegN-
66+
67+
# genmc Pi+Pi- ROOTDataReader tree_AmpToolsFormatThrown_flat_amp.root
68+
# accmc Pi+Pi- ROOTDataReader tree_AmpToolsFormatThrown_flat_amp.root
69+
# data Pi+Pi- ROOTDataReader tree_AmpToolsFormatThrown_amp.root
70+
71+
# amplitude Pi+Pi-::helplusN+::rho1 TwoPiAngles_amp (phipol, pol fraction, Jz for rho M=+/-1 and 0, prefix factor, flat)
72+
# phipol is the lab azimuthal angle of the polarization vector.
73+
# Prefix factor 0: 0.5*sqrt(1-polFrac)*(cos(Phi) - i*sin(Phi)), M=M
74+
# Prefix factor 1: 0.5*sqrt(1+polFrac)*(cos(Phi) - i*sin(Phi)), M=M
75+
# Prefix factor 2: 0.5*sqrt(1-polFrac)*(cos(Phi) + i*sin(Phi)), M=M
76+
# Prefix factor 3: -0.5*sqrt(1+polFrac)*(cos(Phi) + i*sin(Phi)), M=M
77+
# Prefix factor 4: 0.5*(-1)^M*sqrt(1-polFrac)*(cos(Phi) - i*sin(Phi)), M=-M
78+
# Prefix factor 5: 0.5*(-1)^M*sqrt(1+polFrac)*(cos(Phi) - i*sin(Phi)), M=-M
79+
# Prefix factor 6: 0.5*(-1)^M*sqrt(1-polFrac)*(cos(Phi) + i*sin(Phi)), M=-M
80+
# Prefix factor 7: -0.5*(-1)^M*sqrt(1+polFrac)*(cos(Phi) + i*sin(Phi)), M=-M
81+
# flat=1 generates distribution uniform in angles. flat=0 use YLMs
82+
83+
amplitude Pi+Pi-::helplusN+::g1VM1 TwoPiAngles_amp phipol polFrac 1 0 flat
84+
amplitude Pi+Pi-::helplusN+::g1VM1 BreitWigner rho 1 2 3
85+
86+
amplitude Pi+Pi-::helplusN+::g1VM0 TwoPiAngles_amp phipol polFrac 0 0 flat
87+
amplitude Pi+Pi-::helplusN+::g1VM0 BreitWigner rho 1 2 3
88+
89+
amplitude Pi+Pi-::helplusN+::g1VM-1 TwoPiAngles_amp phipol polFrac -1 0 flat
90+
amplitude Pi+Pi-::helplusN+::g1VM-1 BreitWigner rho 1 2 3
91+
92+
amplitude Pi+Pi-::helplusN+::g-1VM1 TwoPiAngles_amp phipol polFrac 1 2 flat
93+
amplitude Pi+Pi-::helplusN+::g-1VM1 BreitWigner rho 1 2 3
94+
95+
amplitude Pi+Pi-::helplusN+::g-1VM0 TwoPiAngles_amp phipol polFrac 0 2 flat
96+
amplitude Pi+Pi-::helplusN+::g-1VM0 BreitWigner rho 1 2 3
97+
98+
amplitude Pi+Pi-::helplusN+::g-1VM-1 TwoPiAngles_amp phipol polFrac -1 2 flat
99+
amplitude Pi+Pi-::helplusN+::g-1VM-1 BreitWigner rho 1 2 3
100+
101+
102+
amplitude Pi+Pi-::helnegN+::g1VM1 TwoPiAngles_amp phipol polFrac 1 1 flat
103+
amplitude Pi+Pi-::helnegN+::g1VM1 BreitWigner rho 1 2 3
104+
105+
amplitude Pi+Pi-::helnegN+::g1VM0 TwoPiAngles_amp phipol polFrac 0 1 flat
106+
amplitude Pi+Pi-::helnegN+::g1VM0 BreitWigner rho 1 2 3
107+
108+
amplitude Pi+Pi-::helnegN+::g1VM-1 TwoPiAngles_amp phipol polFrac -1 1 flat
109+
amplitude Pi+Pi-::helnegN+::g1VM-1 BreitWigner rho 1 2 3
110+
111+
amplitude Pi+Pi-::helnegN+::g-1VM1 TwoPiAngles_amp phipol polFrac 1 3 flat
112+
amplitude Pi+Pi-::helnegN+::g-1VM1 BreitWigner rho 1 2 3
113+
114+
amplitude Pi+Pi-::helnegN+::g-1VM0 TwoPiAngles_amp phipol polFrac 0 3 flat
115+
amplitude Pi+Pi-::helnegN+::g-1VM0 BreitWigner rho 1 2 3
116+
117+
amplitude Pi+Pi-::helnegN+::g-1VM-1 TwoPiAngles_amp phipol polFrac -1 3 flat
118+
amplitude Pi+Pi-::helnegN+::g-1VM-1 BreitWigner rho 1 2 3
119+
120+
121+
amplitude Pi+Pi-::helplusN-::g-1VM1 TwoPiAngles_amp phipol polFrac 1 4 flat
122+
amplitude Pi+Pi-::helplusN-::g-1VM1 BreitWigner rho 1 2 3
123+
124+
amplitude Pi+Pi-::helplusN-::g-1VM0 TwoPiAngles_amp phipol polFrac 0 4 flat
125+
amplitude Pi+Pi-::helplusN-::g-1VM0 BreitWigner rho 1 2 3
126+
127+
amplitude Pi+Pi-::helplusN-::g-1VM-1 TwoPiAngles_amp phipol polFrac -1 4 flat
128+
amplitude Pi+Pi-::helplusN-::g-1VM-1 BreitWigner rho 1 2 3
129+
130+
amplitude Pi+Pi-::helplusN-::g1VM1 TwoPiAngles_amp phipol polFrac 1 6 flat
131+
amplitude Pi+Pi-::helplusN-::g1VM1 BreitWigner rho 1 2 3
132+
133+
amplitude Pi+Pi-::helplusN-::g1VM0 TwoPiAngles_amp phipol polFrac 0 6 flat
134+
amplitude Pi+Pi-::helplusN-::g1VM0 BreitWigner rho 1 2 3
135+
136+
amplitude Pi+Pi-::helplusN-::g1VM-1 TwoPiAngles_amp phipol polFrac -1 6 flat
137+
amplitude Pi+Pi-::helplusN-::g1VM-1 BreitWigner rho 1 2 3
138+
139+
140+
amplitude Pi+Pi-::helnegN-::g-1VM1 TwoPiAngles_amp phipol polFrac 1 5 flat
141+
amplitude Pi+Pi-::helnegN-::g-1VM1 BreitWigner rho 1 2 3
142+
143+
amplitude Pi+Pi-::helnegN-::g-1VM0 TwoPiAngles_amp phipol polFrac 0 5 flat
144+
amplitude Pi+Pi-::helnegN-::g-1VM0 BreitWigner rho 1 2 3
145+
146+
amplitude Pi+Pi-::helnegN-::g-1VM-1 TwoPiAngles_amp phipol polFrac -1 5 flat
147+
amplitude Pi+Pi-::helnegN-::g-1VM-1 BreitWigner rho 1 2 3
148+
149+
amplitude Pi+Pi-::helnegN-::g1VM1 TwoPiAngles_amp phipol polFrac 1 7 flat
150+
amplitude Pi+Pi-::helnegN-::g1VM1 BreitWigner rho 1 2 3
151+
152+
amplitude Pi+Pi-::helnegN-::g1VM0 TwoPiAngles_amp phipol polFrac 0 7 flat
153+
amplitude Pi+Pi-::helnegN-::g1VM0 BreitWigner rho 1 2 3
154+
155+
amplitude Pi+Pi-::helnegN-::g1VM-1 TwoPiAngles_amp phipol polFrac -1 7 flat
156+
amplitude Pi+Pi-::helnegN-::g1VM-1 BreitWigner rho 1 2 3
157+
158+
159+
initialize Pi+Pi-::helplusN+::g1VM1 cartesian 500.0 0.0 real
160+
initialize Pi+Pi-::helplusN+::g1VM0 cartesian 0.0 0.0
161+
initialize Pi+Pi-::helplusN+::g1VM-1 cartesian 0.0 0.0
162+
initialize Pi+Pi-::helplusN+::g-1VM1 cartesian 0.0 0.0
163+
initialize Pi+Pi-::helplusN+::g-1VM0 cartesian 0.0 0.0
164+
initialize Pi+Pi-::helplusN+::g-1VM-1 cartesian 500.0 0.0
165+
166+
initialize Pi+Pi-::helnegN+::g1VM1 cartesian 500.0 0.0
167+
initialize Pi+Pi-::helnegN+::g1VM0 cartesian 0.0 0.0
168+
initialize Pi+Pi-::helnegN+::g1VM-1 cartesian 0.0 0.0
169+
initialize Pi+Pi-::helnegN+::g-1VM1 cartesian 0.0 0.0
170+
initialize Pi+Pi-::helnegN+::g-1VM0 cartesian 0.0 0.0
171+
initialize Pi+Pi-::helnegN+::g-1VM-1 cartesian 500.0 0.0
172+
173+
initialize Pi+Pi-::helplusN-::g1VM1 cartesian 500.0 0.0
174+
initialize Pi+Pi-::helplusN-::g1VM0 cartesian 0.0 0.0
175+
initialize Pi+Pi-::helplusN-::g1VM-1 cartesian 0.0 0.0
176+
initialize Pi+Pi-::helplusN-::g-1VM1 cartesian 0.0 0.0
177+
initialize Pi+Pi-::helplusN-::g-1VM0 cartesian 0.0 0.0
178+
initialize Pi+Pi-::helplusN-::g-1VM-1 cartesian 500.0 0.0
179+
180+
initialize Pi+Pi-::helnegN-::g1VM1 cartesian 500.0 0.0
181+
initialize Pi+Pi-::helnegN-::g1VM0 cartesian 0.0 0.0
182+
initialize Pi+Pi-::helnegN-::g1VM-1 cartesian 0.0 0.0
183+
initialize Pi+Pi-::helnegN-::g-1VM1 cartesian 0.0 0.0
184+
initialize Pi+Pi-::helnegN-::g-1VM0 cartesian 0.0 0.0
185+
initialize Pi+Pi-::helnegN-::g-1VM-1 cartesian 500.0 0.0
186+
187+
# constrain Pi+Pi-::helplusN+::g1VM1 Pi+Pi-::helnegN+::g1VM1
188+
# constrain Pi+Pi-::helplusN+::g1VM0 Pi+Pi-::helnegN+::g1VM0
189+
# constrain Pi+Pi-::helplusN+::g1VM-1 Pi+Pi-::helnegN+::g1VM-1
190+
# constrain Pi+Pi-::helplusN+::g-1VM1 Pi+Pi-::helnegN+::g-1VM1
191+
# constrain Pi+Pi-::helplusN+::g-1VM0 Pi+Pi-::helnegN+::g-1VM0
192+
# constrain Pi+Pi-::helplusN+::g-1VM-1 Pi+Pi-::helnegN+::g-1VM-1
193+
194+
# constrain Pi+Pi-::helplusN+::g1VM1 Pi+Pi-::helplusN-::g1VM1
195+
# constrain Pi+Pi-::helplusN+::g1VM0 Pi+Pi-::helplusN-::g1VM0
196+
# constrain Pi+Pi-::helplusN+::g1VM-1 Pi+Pi-::helplusN-::g1VM-1
197+
# constrain Pi+Pi-::helplusN+::g-1VM1 Pi+Pi-::helplusN-::g-1VM1
198+
# constrain Pi+Pi-::helplusN+::g-1VM0 Pi+Pi-::helplusN-::g-1VM0
199+
# constrain Pi+Pi-::helplusN+::g-1VM-1 Pi+Pi-::helplusN-::g-1VM-1
200+
201+
# constrain Pi+Pi-::helplusN+::g1VM1 Pi+Pi-::helnegN-::g1VM1
202+
# constrain Pi+Pi-::helplusN+::g1VM0 Pi+Pi-::helnegN-::g1VM0
203+
# constrain Pi+Pi-::helplusN+::g1VM-1 Pi+Pi-::helnegN-::g1VM-1
204+
# constrain Pi+Pi-::helplusN+::g-1VM1 Pi+Pi-::helnegN-::g-1VM1
205+
# constrain Pi+Pi-::helplusN+::g-1VM0 Pi+Pi-::helnegN-::g-1VM0
206+
# constrain Pi+Pi-::helplusN+::g-1VM-1 Pi+Pi-::helnegN-::g-1VM-1
207+
208+
209+

‎MCwrapper/Generators/gen_amp/examples/gen_2pi_mom.cfg‎

Lines changed: 100 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,100 @@
1+
#
2+
#The above line is used to provide additional options to the generator within the wrapper. These options do NOT include output file names, run numbers, and numbewr of events which are controlled by the wrapper
3+
#####################################
4+
#### THIS IS A CONFIG FILE ####
5+
#####################################
6+
##
7+
## Blank lines or lines beginning with a "#" are ignored.
8+
##
9+
## Double colons (::) are treated like a space.
10+
## This is sometimes useful for grouping (for example,
11+
## grouping strings like "reaction::sum::amplitudeName")
12+
##
13+
## All non-comment lines must begin with one of the following keywords.
14+
##
15+
## (note: <word> means necessary
16+
## (word) means optional)
17+
##
18+
## include <file>
19+
## define <word> (defn1) (defn2) (defn3) ...
20+
## fit <fitname>
21+
## keyword <keyword> <min arguments> <max arguments>
22+
## reaction <reaction> <particle1> <particle2> (particle3) ...
23+
## data <reaction> <class> (arg1) (arg2) (arg3) ...
24+
## genmc <reaction> <class> (arg1) (arg2) (arg3) ...
25+
## accmc <reaction> <class> (arg1) (arg2) (arg3) ...
26+
## normintfile <reaction> <file>
27+
## sum <reaction> <sum> (sum2) (sum3) ...
28+
## amplitude <reaction> <sum> <amp> <class> (arg1) (arg2) ([par]) ...
29+
## initialize <reaction> <sum> <amp> <"events"/"polar"/"cartesian">
30+
## <value1> <value2> ("fixed"/"real")
31+
## scale <reaction> <sum> <amp> <value or [parameter]>
32+
## constrain <reaction1> <sum1> <amp1> <reaction2> <sum2> <amp2> ...
33+
## permute <reaction> <sum> <amp> <index1> <index2> ...
34+
## parameter <par> <value> ("fixed"/"bounded"/"gaussian")
35+
## (lower/central) (upper/error)
36+
## DEPRECATED:
37+
## datafile <reaction> <file> (file2) (file3) ...
38+
## genmcfile <reaction> <file> (file2) (file3) ...
39+
## accmcfile <reaction> <file> (file2) (file3) ...
40+
##
41+
#####################################
42+
43+
#Take useful definitions from gen_2pi_amp.cfg
44+
define sigma 0.500 0.400
45+
define rho 0.775 0.146
46+
define f2 1.275 0.187
47+
48+
# Uniform angles: flat=1; YLMs: flat=0;
49+
define flat 0
50+
51+
# fit twopi_amp
52+
53+
reaction Pi+Pi- Beam Proton Pi+ Pi-
54+
55+
normintfile Pi+Pi- twopi_ni.txt
56+
57+
# sum for helicity of N'=+1/2. Amplitudes should be duplicated for N'=-1/2.
58+
sum Pi+Pi- Positive
59+
sum Pi+Pi- Negative
60+
61+
# genmc Pi+Pi- ROOTDataReader tree_AmpToolsFormatThrown_flat_amp.root
62+
# accmc Pi+Pi- ROOTDataReader tree_AmpToolsFormatThrown_flat_amp.root
63+
# data Pi+Pi- ROOTDataReader tree_AmpToolsFormatThrown_amp.root
64+
65+
# Negative Reflectivity Waves
66+
67+
amplitude Pi+Pi-::Negative::S0- TwoPSHelicity 0 0 -1
68+
amplitude Pi+Pi-::Negative::S0- BreitWigner rho 1 2 3
69+
70+
amplitude Pi+Pi-::Negative::P0- TwoPSHelicity 1 0 -1
71+
amplitude Pi+Pi-::Negative::P0- BreitWigner rho 1 2 3
72+
73+
amplitude Pi+Pi-::Negative::P1- TwoPSHelicity 1 1 -1
74+
amplitude Pi+Pi-::Negative::P1- BreitWigner rho 1 2 3
75+
76+
amplitude Pi+Pi-::Negative::D0- TwoPSHelicity 2 0 -1
77+
amplitude Pi+Pi-::Negative::D0- BreitWigner rho 1 2 3
78+
79+
amplitude Pi+Pi-::Negative::D1- TwoPSHelicity 2 1 -1
80+
amplitude Pi+Pi-::Negative::D1- BreitWigner rho 1 2 3
81+
82+
# Positive Refectivity Waves:
83+
84+
amplitude Pi+Pi-::Positive::P1+ TwoPSHelicity 1 1 1
85+
amplitude Pi+Pi-::Positive::P1+ BreitWigner rho 1 2 3
86+
87+
amplitude Pi+Pi-::Positive::D1+ TwoPSHelicity 2 1 1
88+
amplitude Pi+Pi-::Positive::D1+ BreitWigner rho 1 2 3
89+
90+
# Initialize: One of each set can be real
91+
92+
initialize Pi+Pi-::Negative::S0- cartesian 0.0 0.0
93+
initialize Pi+Pi-::Negative::P0- cartesian 0.0 0.0
94+
initialize Pi+Pi-::Negative::P1- cartesian 500.0 0.0 real
95+
initialize Pi+Pi-::Negative::D0- cartesian 0.0 0.0
96+
initialize Pi+Pi-::Negative::D1- cartesian 0.0 0.0
97+
98+
initialize Pi+Pi-::Positive::P1+ cartesian 500.0 0.0 real
99+
initialize Pi+Pi-::Positive::D1+ cartesian 0.0 0.0
100+

0 commit comments

Comments
 (0)