Dear Carla.
I will cross-check the coherent weights and get back to you.
Regarding the problematic events when reweighing the formation zone parameter.
The reweighing doesn’t work for particles that we can not propagate through the nucleus (because the weight
calculation depends on a mean free path calculation). In this particular case it chokes on K_L0.
Thanks for sending the event. There is nothing wrong you did. This is a limitation of our intranuclear cascade.
There is no easy fix short of adding methods to handle these particles in our cascade code.
Yes, you can override the default userPhysicsOptions.xml. In fact, you can override any file in the config folder.
If I recall correctly, $GENIE/config is the _last_ directory where GENIE is looking at for config files.
If you create a different directory with config files and set the GXMLPATH env var to point to it,
then GENIE will look at if first. You can put there all the files you want to override. I think this should work.
We are looking to provide additional functionality, as we want to be able to move easily between
different tunes. This new functionality will come with v3.0.0.
Yes, try with Emin at 1 MeV. It should work. Let us know if something goes wrong. It should be relatively
easy to fix (from the point of view of avoiding crashes). A bigger problem is the validity of the model at this range.
There the lower threshold is probably closer to 100-200 MeV.
cheers
Costas
--
Dr. Constantinos Andreopoulos
Reader (Associate Professor) in Experimental Particle Physics
University of Liverpool & STFC Rutherford Appleton Lab
+44-(0)7540-847333 (Mobile)
+44-(0)1235-445091 (Office/RAL)
+44-(0)1517-943201 (Office/Liverpool)
http://costas.andreopoulos.eu
On 18 Feb 2016, at 15:17, carla distefano <[log in to unmask]<mailto:[log in to unmask]>> wrote:
Dear all,
I’m trying to use genie2.10.2 to evaluate the systematics.
I generated muon neutrinos and anti-neutrinos (CC+NC, energy range 1-1000 GeV) and used the ReWeight package to compute the systematic weights referring to single input physics parameters P (setting x_P=-1 and +1). Then I used the systematic weights to find the modified neutrino energy spectra and the ratio of them to the nominal one.
For MaCOHpi and R0COHpi parameters, x_P=1 and x_P=-1 gives the same result. I find a similar behavior also for AGKYpT1pi. This does not happen for the other parameters (e.g MaCCQE). I attached the plots.
Besides when I try to run the ReWeight package for the FormZone parameter, I have a lot of problematic events for which systematic weights are not computed. An example of these problematic events are in the attached file: FormZone.txt
Am I wrong in something?
I have also some questions about technical issues:
- is it possible to override the file UserPhysicsOptions.xml? I’d like to change some settings (for instance the particle decay flags) without modifying the file $GENIE/config/UserPhysicsOptions.xml
- the GENIE validity range (in $GENIE/config/UserPhysicsOptions.xml) is set with Emin=0.010. Is it possible to use GENIE with Emin=0.001?
thanks,
Carla Distefano
<R0COHpi-ratio.gif>
<MaCOHpi-ratio.gif>
<AGKYpT1pi-ratio.gif>
<MaCCQE-ratio.gif>
<FormZone.txt>
|