Error during integration

[LPGaffney]

Hi all,

I am having some trouble during the integration step in Gosia. I have tried for 2 days now and cannot see the reason behind the problem. I had it running perfectly ok with 8 experiments, and recently added a further 2. Since then, I get the error:

Code: Select all

At line 2012 of file gosia_20081208.10.f (unit = 14, file = 'dum.14')
Fortran runtime error: Bad integer for item 2 in list input

It seems to relate to a section of the code which loops over pin diodes even though I do not define any in my experiments and have the azimuthal option set in the EXPT section. The integration runs ok for the original 8 experiments but fails when it reaches the new, 9th experiment. I have tried changing the definitions of the experiments in EXPT and OP,INTI but still the error persists.

I was hoping somebody who uses Gosia more regularly would be able to spot my obvious mistake and help me to move on!!

I have linked to the input files I am using, plus the version of gosia which I edited to include the efficiency curves used (the error occurs even with the release version from the wiki page)
http://ns.ph.liv.ac.uk/~lg/Gosia/


Many thanks,
Liam

[hayes]

Hi Liam,

I tried all the wrong things first, but I think I finally found the problem. The integration seems to be overflowing by trying to access point yields for the wrong detectors.

You have the following in OP,YIEL

Code: Select all

OP,YIEL
0
0,-1
  0.01
.
.
.
experiments 1--8
.
.
.
23
  25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47
   112.9,112.0,132.6, 47.6, 72.4, 67.2, 49.8, 73.7, 62.1,109.2,124.9,130.0,107.4,131.1,119.5, 49.7, 69.6, 55.5,126.7,110.8,104.1, 68.3, 44.5
    126.2,150.7,139.5, 47.4, 54.3, 27.7,135.8,133.5,110.8, 41.3, 59.3, 31.8,251.3,254.5,228.7,239.9,232.0,214.4,330.8,311.9,335.0,335.0,327.7

This defines for experiment 9 logical detector 1 as physical type 25. After this, I think you want to reference only logical detector numbers in OP,RAW. (Check me against the manual or trusted inputs.)

Then in OP,RAW for experiment 9 you have

Code: Select all

OP,RAW
.
.
.
experiments 1--8
.
.
.
9
  0.04348,4.47444,-0.398138,-0.228986,0.17364,-0.0452714,0,0
.
.    "            "             "                "             "            "       
.
  0.04348,4.47444,-0.398138,-0.228986,0.17364,-0.0452714,0,0
   1
    23
25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47


I think that you should replace the line "25, 26,..." with logical detector numbers 1, 2, 3..., so it looks like

Code: Select all

.
.
.
  0.04348,4.47444,-0.398138,-0.228986,0.17364,-0.0452714,0,0
   1   
    23 
1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23


When I did this it gave me integrated yields that didn't overflow. (I stripped out all other experiments except #9 and renumbered it to #1 to isolate the problem, but it looks like experiment #10 has the same problem.)

I hope this helps.

Adam

[LPGaffney]

Hi Adam,

Incredible! Thank you very much. This would've taken me a very long time, maybe forever. I was unaware of the specific differences between logical and physical detectors in Gosia. I still have a lot to learn about this code.

Thanks again for a speedy and accurate response.

Regards,
Liam