test_redmine.py
| 1 |
from ctools import * # this cools is modified by katagiri to obtain loglikelihood value |
|---|---|
| 2 |
import sys |
| 3 |
|
| 4 |
model_name = "model_onlydiffuse_test.xml"
|
| 5 |
#events_name = "events_onlydiffuse_5hr.fits"
|
| 6 |
events_name = "events_onlydiffuse.fits"
|
| 7 |
caldb = "./"
|
| 8 |
irf = "IFAE_2A_Aar_merged_50h.dat"
|
| 9 |
cntmap_name = "ccube_onlydiffuse.fits"
|
| 10 |
|
| 11 |
# Perform maximum likelihood fitting
|
| 12 |
name="onlydiffuse"
|
| 13 |
#logfile="ctlike_"+str(name)+".log"
|
| 14 |
result_name="fit_"+str(name)+".xml" |
| 15 |
|
| 16 |
like = ctlike() |
| 17 |
#like.logFileOpen(logfile)
|
| 18 |
like["infile"].filename(events_name)
|
| 19 |
#like["infile"].filename(cntmap_name)
|
| 20 |
like["srcmdl"].filename(model_name)
|
| 21 |
like["outmdl"].filename(result_name)
|
| 22 |
like["caldb"].string(caldb)
|
| 23 |
like["irf"].string(irf)
|
| 24 |
|
| 25 |
|
| 26 |
"""
|
| 27 |
like.obs().models(model_name)
|
| 28 |
like.obs().models()['DMline']['Normalization'].free()
|
| 29 |
like.get_parameters()
|
| 30 |
like.execute()
|
| 31 |
print "loglike",like.get_loglike()
|
| 32 |
print "Normalization",like.obs().models()['DMline']['Normalization'].value()
|
| 33 |
like.obs().models()['DMline']['Normalization'].value(2.)
|
| 34 |
like.obs().models()['DMline']['Normalization'].fix()
|
| 35 |
like.get_parameters()
|
| 36 |
like.execute()
|
| 37 |
print "loglike",like.get_loglike()
|
| 38 |
print "Normalization",like.obs().models()['DMline']['Normalization'].value()
|
| 39 |
"""
|