Rhosigchi assumtions

[vetlewi]

Rhosigchi will do one or more assumptions about the resolution of the detectors that are used. This should be looked at...

[fzeiser]

See eg:

oslo-method-software/prog/rhosigchi.f

Lines 69 to 71 in 5aafac3

	dE_p = 150. ! Particle resolution
	dE_g1MeV = 60. ! NaI resolution at 1 MeV
	dE_g8MeV = 300. ! NaI resolution at 8 MeV

further used here:

oslo-method-software/prog/rhosigchi.f

Lines 141 to 153 in 5aafac3

	C Now calculating max energy for NaI to cover resolution. Max is 800 keV
	a0 = a0 - INT((800./a1) + 0.5)*a1
	iu0 = INT((ABS(a0/a1) + 0.5))
	DO j=0, 511
	Ex = a0 + j*a1
	dE_g = dE_g1MeV + ((dE_g8MeV-dE_g1MeV)/(8000.-1000.))*(Ex-1000.)
	egmax(j) = SQRT(dE_g*dE_g + dE_p*dE_p)
	IF(egmax(j).GT.600.)egmax(j)=600.
	IF(egmax(j).LT.200.)egmax(j)=200.
	igmax(j) = j + INT((egmax(j)/a1)+0.5)
	IF(igmax(j).GT.511)igmax(j)=511
	egmax(j) = Ex + egmax(j)
	ENDDO

which enters at many places

[cecgresonant]

These parameters will in general vary for each experiment and for different gamma-ray detectors. Probably these parameters should be set by the user and not hardcoded? How is this solved in OMpy?

[fzeiser]

We're still discussing how we want to solve it, but probably we want to
i) force the user to provide these values, or, maybe even more likely
ii) issue a warning if some default values are chosen.

[vetlewi]

We're still discussing how we want to solve it, but probably we want to
i) force the user to provide these values, or, maybe even more likely
ii) issue a warning if some default values are chosen.

I would prefer option ii, as most of the use will be experiments from OCL and in those cases the resolution is more or less the same for each experiment. Just make it very clear to the user that it is using 'default' values.