search

warr / gosia

Gosia Coulomb excitation code
1 stars 1 forks source link

Fitting experiments with no gamma-ray yields copies the next experiment #2

Open jackwoodside opened 1 year ago

jackwoodside commented 1 year ago

I performed an experiment with eight particle detectors, however one of the particle detectors broke so no yields were obtained. I already had the gosia input file set up for eight Coulomb excitation experiments, and wanted to see the trend still in the end results.

When entering the measured yields in file 3, for the broken particle detector I set ND to 0 and did not enter any experimental yields. After running OP,MINI, the output file shows that for the ND = 0 detector it instead just uses the experimental yields for the next detector, causing the fit to change significantly. If the ND = 0 detector is listed last, then the program hangs when running OP,MINI.

lpgaff commented 1 year ago

I suppose Gosia isn't designed to handle the case where an experiment is declared but the number of data are zero. It should probably produce a warning or error, rather than continue with something weird.

Anyway, if you don't want to remove this experiment from the file, I would suggest to add a dummy yield in file 3 and set the weight to zero instead.

jackwoodside commented 1 year ago

I agree that an error or even a warning would be good.

For the mean time I've just made a separate input file without the experiment, but thanks for the tip about setting the weighting, that would be simpler if this happens again in future.