fermi-lat / Likelihood

BSD 3-Clause "New" or "Revised" License
2 stars 1 forks source link

ERROR: Caught St13runtime_error at the top level: RoiCuts::sortCuts: #111

Open nyctophile-aayu opened 2 years ago

nyctophile-aayu commented 2 years ago

Hi, I am getting the similar issue, even when I work without python in Unbinned LIkelihood Analysis. Please check the error below:

gtexpmap The exposure maps generated by this tool are meant to be used for unbinned likelihood analysis only. Do not use them for binned analyses. Event data file[] 1633+382_region_filtered_gti.fits Spacecraft data file[] spacecraft.fits Exposure hypercube file[] 1633+382_ltcube.fits output file name[] 1633+382_unbin_expmap.fits Response functions[CALDB] Radius of the source region (in degrees)[30] 25 Number of longitude points (2:1000) [120] 120 Number of latitude points (2:1000) [120] 120 Number of energies (2:100) [20] 30 Caught St13runtime_error at the top level: RoiCuts::sortCuts: There should be exactly one energy range cut, one acceptance cone cut, and at least one time range and/or GTI cut. The event file contains the following DSS selections:

DSTYP1: BIT_MASK(EVENT_CLASS,128,P8R3) DSUNI1: DIMENSIONLESS DSVAL1: 1:1

DSTYP2: POS(RA,DEC) DSUNI2: deg DSVAL2: CIRCLE(248.814554,38.134583,15)

DSTYP3: POS(RA,DEC) DSUNI3: deg DSVAL3: CIRCLE(248.815,38.1346,15)

DSTYP4: BIT_MASK(EVENT_TYPE,3,P8R3) DSUNI4: DIMENSIONLESS DSVAL4: 1:1

DSTYP5: ENERGY DSUNI5: MeV DSVAL5: 100:300000

DSTYP6: ZENITH_ANGLE DSUNI6: deg DSVAL6: 0:90

DSTYP7: TIME DSUNI7: s DSVAL7: TABLE DSREF7: :GTI

GTIs: (suppressed)

For both of the analysis, I am using these values for source: 1633+382

Equatorial coordinates (degrees) | (248.815,38.1346) Time range (MET) | (660799526,676351526) Time range (Gregorian) | (2021-12-10 03:25:21,2022-06-08 03:25:21) Energy range (MeV) | (100,300000) Search radius (degrees) | 15

I've tried as suggested in a closed post "https://github.com/fermi-lat/Likelihood/issues/107#issuecomment-1150241157". But, the error persists. These errors are on repeated mode, even if I change the values. Please suggest what and where I am going wrong. Plus, why am I getting runtime error so frequently?

Thanks in advance -Ayush

nmirabal commented 2 years ago

You have two ROI selections in the header of the 1633+382_region_filtered_gti.fits file. Try deleting one of the following two entries from the header(see below) or redo the analysis from the original event files keeping the same RA,Dec throughout your analysis. DSTYP2: POS(RA,DEC) DSUNI2: deg DSVAL2: CIRCLE(248.814554,38.134583,15)

DSTYP3: POS(RA,DEC) DSUNI3: deg DSVAL3: CIRCLE(248.815,38.1346,15)