Coverage for src / lstautorta / Analyse_DL3.py: 0%

1import warnings 

2 

3# import DL2_Analysis_Wobble 

4from gammapy.data import DataStore 

5from SourceAnalyse import event 

6 

7warnings.filterwarnings("ignore") 

8warnings.simplefilter("ignore") 

9 

10 

11def CreatePlots(id_=0): 

12 directory = "/home/sami.caroff/DL3_Fake_RTA/mergedDL3/" 

13 

14 datastore = DataStore.from_dir(directory) 

15 

16 Analysis = event(directory, [9952]) 

17 # Analysis = event(directory,[9917,9918,9919,9920]) 

18 

19 # Analysis.create_on_region(330.68,42.28) #BL Lac 

20 Analysis.create_on_region(166.108, 38.21) # Mrk421 

21 # Analysis.create_on_region() 

22 Analysis.create_exclusion_mask() 

23 Analysis.reduction_chain() 

24 Analysis.data_maker() 

25 Analysis.data_run() 

26 Analysis.signal_info() 

27 Analysis.fit_spectrum() 

28 

29 Analysis.plot_timee(id_=id_) 

30 excess_reflected, excess_reflected_uncertainty, c, d, sig_reflected, sig_livetime = Analysis.signal_info() 

31 Analysis.flux_plot(plot=False, id_=id_) 

32 

33 Analysis.calculate_acceptance_model(plot=False) 

34 Analysis.add_acceptance_map() 

35 Analysis.map_geometry(plot=False) 

36 ring_excess, ring_sqrt_ts, ring_bg, uncertainty_excess_ring = Analysis.ring_data_estimation() 

37 Analysis.excess_significance_plot(plot=True, id_=id_) 

38 # Analysis.off_distribution_plot(plot=False) 

39 return ( 

40 excess_reflected, 

41 excess_reflected_uncertainty, 

42 sig_reflected, 

43 sig_livetime, 

44 ring_excess, 

45 uncertainty_excess_ring, 

46 ring_sqrt_ts, 

47 )