Bug #1361

GCTAPointing not initialized

Added by Lu Chia-Chun over 9 years ago. Updated over 9 years ago.

Status:ClosedStart date:11/10/2014
Priority:HighDue date:
Assigned To:Knödlseder Jürgen% Done:

100%

Category:-
Target version:00-09-00
Duration:

Description

After recent update, I find GCTAPointing not initialized when an observation is loaded.

For example,

obs = GObservations('obs.xml’)
pnt = obs0.pointing()
pntdir = pnt.dir()
pntdir.dec_deg()

the last line gives 0

It also shows in ctlike log message that
2014-11-10T12:53:27: === GCTAPointing ===
2014-11-10T12:53:27: Pointing direction ........: (RA,Dec)=(not initialised)

Is this a bug?


Recurrence

No recurrence.

History

#1 Updated by Mayer Michael over 9 years ago

  • Description updated (diff)

The problem you describe is in GCTAObservation::read(const GXmlElement& xml). Due to a recent change, line 694 (where the events are loaded) is commented out. A quick fix for you might be to uncomment and recompile.
I am not sure but I think this line has been removed due to a problem with OpenMP and loading events simultaneously.

We might want to change the current way that instead of loading the events, we could only read the fits header to initialise observational attributes like pointing, livetime, etc. on reading. Instead of load, we could call therefore call read_attributes() to initialise all necessary information.

#2 Updated by Lu Chia-Chun over 9 years ago

  • Status changed from New to Feedback

Hi Michael,

thank you. Uncommenting out that line works and ctlike doesn’t crash due to the old problem of openMP.?So I don’t know why this line was commented out. Let’s wait for Juergen’s reply.

And I agree with your suggestion!

#3 Updated by Knödlseder Jürgen over 9 years ago

The reason for the change was a reduction of memory needs. Now the events are only loaded when needed. I agree that we need to add a functionality for loading the pointing information from the header. I did already add some code to load the energy and ROI, but I forgot about the pointing ...

#4 Updated by Knödlseder Jürgen over 9 years ago

  • Status changed from Feedback to Closed
  • Assigned To set to Knödlseder Jürgen
  • Target version set to 00-09-00
  • % Done changed from 0 to 100

I corrected the code and merged it into the devel branch.

Also available in: Atom PDF