-
Notifications
You must be signed in to change notification settings - Fork 105
NC Coherent Gamma production #123
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
base: master
Are you sure you want to change the base?
Changes from 250 commits
2ddcc37
5964a1a
8f9cd91
3fea82b
c2e43a1
a527f99
bd3deab
754ffc5
da168cd
97cfa51
883227a
60d52dd
8f65584
3fe30da
2c551f6
e415836
c9bcfc5
7f1727e
2e5f420
adba501
0a2edc2
eace982
1fa9cd4
c109ddd
7abbf8a
5ca1e3c
14b3160
3522032
6d50479
0538ae1
d40db41
ee77214
bc682e7
138c06c
75b0134
3f93d2d
d8039aa
4331cd6
77d049f
c909875
1bba81f
b19a235
82f0757
ee75378
7393982
f839237
3955c41
4da49f0
6a6ebbd
9f8045d
a7982f8
3fae956
c78a37b
1ecef55
da45574
3225d9f
17e0576
81f902c
a2aa190
0da3090
7a3a5b5
315b0a8
ff32829
6e17e8f
e4c9611
b40bc19
bb22727
b69ff06
243efc2
1aebecf
7fbdb2d
a880768
a266b12
2d15098
cd47ca5
403b497
613f0a2
5863e0d
50773e7
cacd0d8
01d947b
42979f8
6327d67
58fac75
2f38b02
b4f1c46
f5e5e4a
7319256
9c29176
cf2a1b6
f8fb08d
e274fcf
dbf4369
91af15c
c6f8f17
63d5744
95ca193
e06ea5e
1458730
5f7a913
9a9075d
7db130f
0c9a899
2016da8
7e7a257
0966820
2a23f43
9db7d9a
20b4b8f
77f7206
740d39f
c399e32
1ce6514
b3af309
af6db3a
c036db0
47019c7
f725ddb
bb7663e
0fe6c03
7642950
bf824ac
5a676ac
ac879af
d7145ac
6dca4a8
1c05712
2653e0f
b8b377a
1642369
4c4745e
8d1e9e5
11fb087
8fa48d9
2564bf0
16027cb
830bdaf
7b94bc1
2215155
d4d9002
f4bd90e
43e2ef4
03a3cc3
df80f76
a4df6a0
6c8719f
acb5533
72bc284
6f867b7
9a0deaf
bd6ebc6
89042b8
bf5084b
f4d8543
c4e2f0c
4f40399
0d6ce86
567a2e6
3e29516
350e569
a486b1c
951cd89
521280d
aa4c877
2880820
5faefb2
bfb09d1
49fb076
75475b2
a9fee75
5ec597d
c00fbfa
2c853dc
0e2f14d
81b672c
d5b51ec
3a8db67
bb35ac3
6f619b0
0deb930
a183172
15f8c2b
43794c6
28552b2
a598f48
61e667b
3954d65
5511631
710c7eb
dafef00
4884789
bd247a7
8fe8067
bf76e85
19305f7
f53a903
b602bda
50bee4f
70a26a6
cca6659
945e386
bf01390
d54a738
1d4a808
7b0f533
f2ec7af
96056d9
46779bc
eefed53
b88d53a
d6dd41a
1a8b621
23259c5
224d2e6
f8debf0
8f5a192
f62a7cf
379115d
e315783
46d0b32
8f99b6a
fb9da98
57bbe2c
6b9ca47
34cd86b
ae54a93
a6d733a
471fd3e
5add956
45db318
9e0fa9e
9f59c10
0e710fc
9891881
faea194
b568fd0
f3a1982
bfcd3d8
7e28865
5a873bb
aefffa1
59a0432
995bf4b
063568a
e65bac2
f5a085c
a1749e7
a657aaf
4c1518d
a2525e9
File filter
Filter by extension
Conversations
Jump to
Diff view
Diff view
There are no files selected for viewing
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,43 @@ | ||
| <?xml version="1.0" encoding="ISO-8859-1"?> | ||
|
|
||
| <!-- | ||
|
|
||
| Configuration sets for the AlvarezRusoSalaCOHGammaPXSec | ||
|
|
||
| Configurable Parameters: | ||
| ......................................................................................................................... | ||
| Name Type Optional Comment Default | ||
| ......................................................................................................................... | ||
| COH-HadronicCurrents vec-alg No Vector of hadronic current sub algos - | ||
| XSec-Integrator alg No Cross section integrator sub algo - | ||
| COH-FormFactor alg yes DeVries form facor sub algo The one of the integrator | ||
| AsymptoticMaxGammaEnergy double No Configuration parameter for maximum W for taken from subalgo | ||
| the cross section. | ||
| --> | ||
|
|
||
|
|
||
|
|
||
|
|
||
| <alg_conf> | ||
|
|
||
| <param_set name="Default"> | ||
|
|
||
| <param type="vec-alg" name="COH-HadronicCurrents" delim=";"> | ||
| genie::COHDeltaCurrent/Default | ||
| </param> | ||
|
|
||
| <param type="alg" name="XSec-Integrator"> genie::COHXSecAR/Gamma </param> | ||
|
|
||
| </param_set> | ||
|
|
||
| <param_set name="Consistent"> | ||
|
|
||
| <param type="vec-alg" name="COH-HadronicCurrents" delim=";"> | ||
| genie::COHDeltaCurrent/FermiMomentumTable | ||
| </param> | ||
|
|
||
| </param_set> | ||
|
|
||
|
|
||
|
|
||
| </alg_conf> |
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,33 @@ | ||
| <?xml version="1.0" encoding="ISO-8859-1"?> | ||
|
|
||
| <!-- | ||
|
|
||
| Configuration sets for the COHDeltaCurrent | ||
|
|
||
| Configurable Parameters: | ||
| .......................................................................................................... | ||
| Name Type Optional Comment Default | ||
| .......................................................................................................... | ||
| COH-Delta-FormFactor alg No Delta Transistion form facort sub algo - | ||
| COH-Delta-Corrections alg No Delta Corrections in medium sub algo - | ||
| --> | ||
|
|
||
|
|
||
|
|
||
| <alg_conf> | ||
|
|
||
| <param_set name="Default"> | ||
|
|
||
| <param type="alg" name="COH-Delta-FormFactor"> genie::DeltaTransitionFormFactor/Default </param> | ||
| <param type="alg" name="COH-Delta-Corrections"> genie::DeltaInMediumCorrections/Default </param> | ||
|
|
||
| </param_set> | ||
|
|
||
| <param_set name="FermiMomentumTable"> | ||
|
|
||
| <param type="alg" name="COH-Delta-Corrections"> genie::DeltaInMediumCorrections/FermiMomentumTable </param> | ||
|
|
||
| </param_set> | ||
|
|
||
|
|
||
| </alg_conf> |
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,32 @@ | ||
| <?xml version="1.0" encoding="ISO-8859-1"?> | ||
|
|
||
|
|
||
|
|
||
| <!-- | ||
|
|
||
| Configuration sets for the COHGammaIntegrationLimits | ||
|
|
||
| Configurable Parameters: | ||
| ..................................................................................................................... | ||
| Name Type Optional Comment Default | ||
| ..................................................................................................................... | ||
| COH-FormFactor alg No form facor sub algo - | ||
| AsymptoticMaxGammaEnergy double No Energy cut off for gamma energy [GeV] | ||
| in the limit for infinit mass nuclei | ||
| MaxGammaTheta double yes Theta cut off for gamma [deg] 180-epsilon | ||
|
|
||
| --> | ||
|
|
||
|
|
||
|
|
||
| <alg_conf> | ||
|
|
||
| <param_set name="Default"> | ||
| <param type="double" name="AsymptoticMaxGammaEnergy"> 1.4 </param> | ||
| <param type="double" name="MaxGammaTheta"> 150. </param> | ||
|
|
||
| <param type="alg" name="COH-FormFactor"> genie::COHProtonFormFactorInterpolation/Default </param> | ||
|
|
||
| </param_set> | ||
|
|
||
| </alg_conf> |
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,45 @@ | ||
| <?xml version="1.0" encoding="ISO-8859-1"?> | ||
|
|
||
| <!-- | ||
| Configuration for the COHKinematicsGenerator EventRecordVisitorI | ||
| Algorithm Configurable Parameters: | ||
| ....................................................................................................................... | ||
| Name Type Opt Comment Default | ||
| ....................................................................................................................... | ||
| UniformOverPhaseSpace bool Yes kinematics uniformly over allowd phase space false | ||
| wgt = (phase_space_volume)*(diff_xsec)/(xsec) | ||
| MaxXSec-SafetyFactor double Yes multiplies max xsec in rejection method 1.6 | ||
| MaxXSec-DiffTolerance double Yes max allowed 200*(xsec-xsecmax)/(xsec+xsecmax) 999999.00 (disable) | ||
| if xsec>xsecmax | ||
| Cache-MinEnergy double Yes minimum energy for which max xsec is cached 0.00 | ||
| MinimScanPoint vec-int yes number of points in each phase space all 1 | ||
| coordinate to be used for a scan | ||
| for Maximum XSec | ||
| IntegrationLimits algo No Specificies integration limit algorithm | ||
| It is required for COH gamma configurations | ||
| tPhaseSpace bool yes Generate the events in the t phase space false | ||
| instead of theta_l | ||
| MaxIterations int yes Maximum number of iterations to try controls::kRjMaxIterations | ||
| acceptance / rejection method | ||
| --> | ||
|
|
||
| <alg_conf> | ||
|
|
||
| <param_set name="Default"> | ||
|
|
||
| <param type="double" name="MaxXSec-SafetyFactor"> 1.05 </param> | ||
|
|
||
| <param type="int" name="MaxIterations"> 10000 </param> | ||
|
|
||
| <param type="vec-int" name="MinimScanPoint" delim=";" > 10 ; 200 ; 10 ; 1 </param> | ||
|
|
||
| <param type="alg" name ="IntegrationLimits" > genie::COHGammaIntegrationLimits/Default </param> | ||
|
|
||
| </param_set> | ||
|
|
||
|
|
||
| </alg_conf> | ||
|
|
||
This file was deleted.
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,38 @@ | ||
| <?xml version="1.0" encoding="ISO-8859-1"?> | ||
|
|
||
| <!-- | ||
| Configuration for the COHPionKinematicsGenerator EventRecordVisitorI | ||
|
|
||
| Algorithm Configurable Parameters: | ||
| ....................................................................................................................... | ||
| Name Type Opt Comment Default | ||
| ....................................................................................................................... | ||
| UniformOverPhaseSpace bool Yes kinematics uniformly over allowd phase space false | ||
| wgt = (phase_space_volume)*(diff_xsec)/(xsec) | ||
| MaxXSec-SafetyFactor double Yes multiplies max xsec in rejection method 1.6 | ||
| MaxXSec-DiffTolerance double Yes max allowed 200*(xsec-xsecmax)/(xsec+xsecmax) 999999.00 (disable) | ||
| if xsec>xsecmax | ||
| Cache-MinEnergy double Yes minimum energy for which max xsec is cached 0.00 | ||
|
|
||
| COH-Ro double No Nuclear size scale CommonParam[Coherent] | ||
| COH-Q2-min double No Minimum considered Q^2 for Berger-Sehgal CommonParam[Coherent] | ||
| coherent reactions when estimating the max | ||
| cross section. Units in GeV^2. | ||
| COH-Q2-max double No Maximum considered Q^2 for Berger-Sehgal CommonParam[Coherent] | ||
| coherent reactions when estimating the max | ||
| cross section. Units in GeV^2. | ||
| COH-t-max double No Maximum considered t for Berger-Sehgal CommonParam[Coherent] | ||
| coherent reactions when estimating | ||
| the max cross section. Units in GeV^2. | ||
| --> | ||
|
|
||
| <alg_conf> | ||
|
|
||
| <param_set name="Default"> | ||
|
|
||
| <param type="string" name="CommonParam"> Coherent </param> | ||
|
|
||
| </param_set> | ||
|
|
||
| </alg_conf> | ||
|
|
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,27 @@ | ||
| <?xml version="1.0" encoding="ISO-8859-1"?> | ||
|
|
||
| <!-- | ||
|
|
||
| Configuration sets for the COHProtonFormFactorInterpolation | ||
|
|
||
| Configurable Parameters: | ||
| ......................................................................................................................... | ||
| Name Type Optional Comment Default | ||
| ......................................................................................................................... | ||
| BaseCOHFormFactor alg No Base form factor algorithm | ||
| AllowExtrapolation bool yes Allow map extrapolation at higher pdg than map maximum false | ||
| --> | ||
|
|
||
|
|
||
|
|
||
| <alg_conf> | ||
|
|
||
| <param_set name="Default"> | ||
| <param type="alg" name="BaseCOHFormFactor" > | ||
| genie::DeVriesFormFactorMap/Default | ||
| </param> | ||
|
|
||
| <param type="bool" name="AllowExtrapolation"> false </param> | ||
| </param_set> | ||
|
|
||
| </alg_conf> |
| Original file line number | Diff line number | Diff line change |
|---|---|---|
|
|
@@ -12,6 +12,17 @@ split-integral bool yes If true, breaks 4d integral into 3d+ | |
| gsl-integration-type string yes Algorithm to use for multidimensional integral vegas | ||
| gsl-relative-tolerance double yes Desired numerical accuracy for each integral 0.01 | ||
| gsl-max-eval int yes Max limit of evaluations for multidimensional integral 20000 | ||
|
|
||
| IsCOHPion bool yes This integrator can be used for both COH Pion and Gamma false | ||
|
There was a problem hiding this comment. Choose a reason for hiding this commentThe reason will be displayed to describe this comment to others. Learn more. Is this description still right in light of IsCOHGamma? (also case for IsCOHGamma vs ISCOHGamma) There was a problem hiding this comment. Choose a reason for hiding this commentThe reason will be displayed to describe this comment to others. Learn more. Yes. I didn't like that a change in the name of the algorithms can case a change in the behaviour so I introduced dedicated flags for every type of interaction. This algo can then integrate both the COH Pion cross section from Luis and the COH Gamma as well. If we don't like them we can make a dedicated integrator class for the COH Gamma as we did for the generator. Shall I go that way? |
||
| ISCOHGamma bool yes false | ||
| OmegaPhaseSpace bool Yes Only used for the COH Gamma case false | ||
| It perform the integration in the dOmega_g dOmega_l | ||
| tPhaseSpace bool Yes Only used for the COH Gamma case false | ||
| It perform the integration in the t phase space | ||
| instead of theta_l | ||
| IntegrationLimits algo Maybe Specificies integration limit algorithm | ||
| It is required for COH gamma configurations | ||
|
|
||
| .................................................................................................... | ||
|
|
||
| --> | ||
|
|
@@ -22,18 +33,58 @@ gsl-max-eval int yes Max limit of evaluations for multidi | |
| <param type="string" name ="gsl-integration-type"> miser </param> | ||
| <param type="int" name ="gsl-max-eval" > 40000 </param> | ||
| <param type="double" name ="gsl-relative-tolerance" > 0.01 </param> | ||
| <param type="bool" name ="split-integral" > true </param> | ||
| <param type="bool" name ="split-integral" > true </param> | ||
mroda88 marked this conversation as resolved.
Outdated
Show resolved
Hide resolved
|
||
|
|
||
| </param_set> | ||
|
|
||
| <param_set name="Fast"> | ||
|
|
||
| <param_set name="Pion"> | ||
| <param type="bool" name ="IsCOHPion" > true </param> | ||
| </param_set> | ||
|
|
||
|
|
||
| <param_set name="Fast_Pion"> | ||
| <!-- <param type="alg" name="Integrator"> genie::Simpson2D/Default-Logarithmic-FixedNBins-101 </param> --> | ||
|
|
||
| <param type="string" name ="gsl-integration-type"> miser </param> | ||
| <param type="int" name ="gsl-max-eval" > 10000 </param> | ||
| <param type="double" name ="gsl-relative-tolerance" > 0.05 </param> | ||
|
|
||
| <param type="bool" name ="IsCOHPion" > true </param> | ||
|
|
||
| </param_set> | ||
|
|
||
| <param_set name="Gamma"> | ||
|
|
||
| <param type="bool" name ="IsCOHGamma" > true </param> | ||
|
|
||
| <param type="bool" name ="OmegaPhaseSpace" > false </param> | ||
|
|
||
| <param type="bool" name ="split-integral" > false </param> | ||
| <param type="string" name ="gsl-integration-type"> vegas </param> | ||
|
|
||
| <param type="int" name ="gsl-max-eval" > 400000 </param> | ||
| <param type="double" name ="gsl-relative-tolerance" > 0.005 </param> | ||
|
|
||
| <param type="alg" name ="IntegrationLimits" > genie::COHGammaIntegrationLimits/Default </param> | ||
|
|
||
| </param_set> | ||
|
|
||
|
|
||
| <param_set name="Gamma_OmegaPhaseSpace"> | ||
| <param type="bool" name ="IsCOHGamma" > true </param> | ||
| <param type="bool" name ="OmegaPhaseSpace" > true </param> | ||
| <param type="bool" name ="split-integral" > false </param> | ||
| <param type="string" name ="gsl-integration-type"> vegas </param> | ||
|
|
||
| <param type="int" name ="gsl-max-eval" > 400000 </param> | ||
| <param type="double" name ="gsl-relative-tolerance" > 0.005 </param> | ||
|
|
||
|
|
||
| </param_set> | ||
|
|
||
|
|
||
|
|
||
|
|
||
| </alg_conf> | ||
|
|
||
Uh oh!
There was an error while loading. Please reload this page.