Add a function to the IEventList interface to extract the event weights

[Martyn Gigg]

and add it to the Python interface.

Similarly to how the getTofs method works.

[Nick Draper]

re #7208 Add getWeight and getWeightErrors to eventlist

This is also accessible via Python

el = ws.getEventList(0)
weightList = el.getWeights()

print len(weightList) # should contain the same number of weights as events in the event list
print weightList[0] #first value

Changeset: ae82a054862a4b4be586130b1f423db4f162b95d

[Nick Draper]

Events without weights will always return event weights of 1 and errors of 1. Weighted events will return the weight and weight error values.

[Nick Draper]

To test:

1. check the unit tests pass, particularly the eventlist unit tests and python unit tests.
2. run the following script

   ws=Load("CNCS_7860_event.nxs")
   el=ws.getEventList(37184)
   wl=el.getWeights()
   wel=el.getWeightErrors()
   print "count =",el.getNumberEvents(), "=", len(wl), "=", len(wel)
   print "weight sum =",numpy.sum(wl)
   print "error sum =",numpy.sum(wel)
   cws=CompressEvents(ws)
   cel=cws.getEventList(37184)
   cwl=cel.getWeights()
   cwel=cel.getWeightErrors()
   print "Compressed count =",cel.getNumberEvents(), "=", len(cwl), "=", len(cwel)
   print "Compressed weight sum =",numpy.sum(cwl)
   print "Compressed error sum =",numpy.sum(cwel)
   

3. you should get these answers

   ws=Load("CNCS_7860_event.nxs")
   el=ws.getEventList(37184)
   wl=el.getWeights()
   wel=el.getWeightErrors()
   print "count =",el.getNumberEvents(), "=", len(wl), "=", len(wel)
   print "weight sum =",numpy.sum(wl)
   print "error sum =",numpy.sum(wel)
   cws=CompressEvents(ws)
   cel=cws.getEventList(37184)
   cwl=cel.getWeights()
   cwel=cel.getWeightErrors()
   print "Compressed count =",cel.getNumberEvents(), "=", len(cwl), "=", len(cwel)
   print "Compressed weight sum =",numpy.sum(cwl)
   print "Compressed error sum =",numpy.sum(cwel)
   

The errors differ due the the processing in CompressEvents that cannot be simply mimicked.

[Peter Peterson]

Two reasons to bounce so far:

• getTofs() no longer works in python
• The returns should use the numpy helper to convert everything into numpy arrays

[Nick Draper]

re #7208 wrapped arrays in numpy and fixed getTofs

GetTofs fails due a a line deleted by mistake (I'm suprised it compiled). Wrapped getweight, getweighterrors and getTofs in numpy arrays Added py unit test for getTofs

Changeset: 878792ae547fc2b8b044bdcbbd08f15f92fa8552

[Nick Draper]

To test:

1. check the unit tests pass, particularly the eventlist unit tests and python unit tests.
2. run the following script

   ws=Load("CNCS_7860_event.nxs")
   el=ws.getEventList(37184)
   wl=el.getWeights()
   wel=el.getWeightErrors()
   tofs=el.getTofs()
   print "count =",el.getNumberEvents(), "=", len(wl), "=", len(wel)
   print "weight sum =",numpy.sum(wl)
   print "error sum =",numpy.sum(wel)
   print "tofs sum =",numpy.sum(tofs)
   cws=CompressEvents(ws)
   cel=cws.getEventList(37184)
   cwl=cel.getWeights()
   cwel=cel.getWeightErrors()
   ctofs=cel.getTofs()
   print "Compressed count =",cel.getNumberEvents(), "=", len(cwl), "=", len(cwel)
   print "Compressed weight sum =",numpy.sum(cwl)
   print "Compressed error sum =",numpy.sum(cwel)
   print "Compressed tofs sum =",numpy.sum(ctofs)
   

3. should give the answers

   count = 175 = 175 = 175
   weight sum = 175.0
   error sum = 175.0
   tofs sum = 9194645.59766
   Compressed count = 155 = 155 = 155
   Compressed weight sum = 175.0
   Compressed error sum = 163.18789493
   Compressed tofs sum = 8143774.69531
   

The errors differ due the the processing in CompressEvents that cannot be simply mimicked.

[Peter Peterson]

All good now.

[Stuart Campbell]

This ticket has been transferred to github issue ​8054