1 / 37

Willman-1 preliminary analysis (data set 01-02/2008) Saverio Lombardi, INFN - Padova 2008/03/20

Willman-1 preliminary analysis (data set 01-02/2008) Saverio Lombardi, INFN - Padova 2008/03/20. Data set 01-02/2008.

yeo-fernandez
Download Presentation

Willman-1 preliminary analysis (data set 01-02/2008) Saverio Lombardi, INFN - Padova 2008/03/20

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Willman-1 preliminary analysis(data set 01-02/2008)Saverio Lombardi, INFN - Padova 2008/03/20

  2. Data set 01-02/2008 ON data:seq. 00329294 (2008-01-17)seq. 00334000 (2008-02-05)seq. 00334125 (2008-02-05)seq. 00335115 (2008-02-07)seq. 00335160 (2008-02-07)seq. 00335518 (2008-02-08)seq. 00335627 (2008-02-08)Total time: 5.6 hours, 22<Zd<30, PSF ~ 11mmOFF data:seq. 00334066 (2008-02-05)seq. 00335155 (2008-02-07)seq. 00335574 (2008-02-08)Total time: 1.9 hours, 22<Zd<26, PSF ~ 11mm

  3. Analysis chain (I) Image cleaning: ABS(6-3) + time info (STD cleaning):MJStar.MImgCleanStd.UseTime: yesMJStar.MImgCleanStd.MaxTimeOff: 4.5MJStar.MImgCleanStd.MaxTimeDiff: 1.5MJStar.MImgCleanStd.CleanLevel1: 6MJStar.MImgCleanStd.CleanLevel2: 3MJStar.MImgCleanStd.CleanMethod: Absolute(search for best settings in progress)Applied cuts:1.5-4.*log10(MNewImagePar.fConc)<log10(MHillas.fSize)MHillas.fSize<80MHillas.fSize>1000000MNewImagePar.fLeakage1>0.2MNewImagePar.fNumCorePixels<4MImagePar.fNumIslands>3MHillasSrc.fDist<40MHillasSrc.fDist>340Image cleaning applied for each sequence separately

  4. Analysis chain (II) Osteria RF /h separation parameters (source-DEP): 1) log10(MHillas.fSize)2) pow((MHillas.fMeanX-MSrcPosCam.fX),2)+pow((MHillas.fMeanY-MSrcPosCam.fY),2)3) MHillas.fWidth4) MHillas.fLength5) log10(MHillas.fSize/(MHillas.fWidth*MHillas.fLength))6) MHillasExt.fM3Long*sgn(MHillasSrc.fCosDeltaAlpha)7) MNewImagePar.fConc8) MHillasSrc.fDist9) MHillasTimeFit.fP1Grad*sgn(MHillasSrc.fCosDeltaAlpha)10) MHillasTime.fRMSTime Osteria RF En.-rec. parameters (source-DEP):1) log10(MHillas.fSize)2) MHillas.fWidth3) MHillas.fLength4) log10(MHillas.fSize/(MHillas.fLength*MHillas.fWidth))5) MNewImagePar.fConc6) MNewImagePar.fLeakage17) MHillasTimeFit.fP1Grad*sgn(MHillasSrc.fCosDeltaAlpha)8) MHillasSrc.fDist9) MPointingPos.fZd10) MMcEvt.fEnergy(search for best settings in progress)

  5. Real data Zd distributions and MC selection ON OFF MC selection: =-2.6, MUX-FADC, PSF=10.6mm, zbins: 8-12 (23°-29°) Train sub-sample 20%, Test sub-sample 80%

  6. DTs (I) 00329294 (ON) 00334000 (ON) 00334066 (OFF) 00334125 (ON) 00335115 (ON) 00335155 (OFF)

  7. DTs (II) 00335160 (ON) 00335518 (ON) 00335574 (OFF) 00335627 (ON) mean DT value ~ 16  1

  8. Rates (I) 00329294 (ON) 00334000 (ON) 00334066 (OFF) 00334125 (ON) 00335115 (ON) 00335155 (OFF)

  9. Rates (II) 00335160 (ON) 00335518 (ON) 00335574 (OFF) 00335627 (ON) mean rate value ~ (20030) Hz

  10. Efficiency of the cleaning (E > 80 GeV, hadcut=0.2)

  11. Phi-COG-Pos-plots: Size80-200, hadcut=1 -plots: 1) ON (red), OFF(blue) 2) MC (black) 3) diff.% ON-OFF(norm.) COG-plots: 1) ON 2) OFF 3) diff. ON-OFF(norm.) Source-Position-plots: (red circle = central pixel) 1) ON 2) OFF 3) diff. ON-OFF(norm.)

  12. Phi-COG-Pos-plots: Size80-200, hadcut=0.2 -plots: 1) ON (red), OFF(blue) 2) MC (black) 3) diff.% ON-OFF(norm.) COG-plots: 1) ON 2) OFF 3) diff. ON-OFF(norm.) Source-Position-plots: 1) ON 2) OFF 3) diff. ON-OFF(norm.)

  13. Phi-COG-Pos-plots: Size200-800, hadcut=1 -plots: 1) ON (red), OFF(blue) 2) MC (black) 3) diff.% ON-OFF(norm.) COG-plots: 1) ON 2) OFF 3) diff. ON-OFF(norm.) Source-Position-plots: 1) ON 2) OFF 3) diff. ON-OFF(norm.)

  14. Phi-COG-Pos-plots: Size200-800, hadcut=0.1 -plots: 1) ON (red), OFF(blue) 2) MC (black) 3) diff.% ON-OFF(norm.) COG-plots: 1) ON 2) OFF 3) diff. ON-OFF(norm.) Source-Position-plots: 1) ON 2) OFF 3) diff. ON-OFF(norm.)

  15. Phi-COG-Pos-plots: Size800-10000, hadcut=1 -plots: 1) ON (red), OFF(blue) 2) MC (black) 3) diff.% ON-OFF(norm.) COG-plots: 1) ON 2) OFF 3) diff. ON-OFF(norm.) Source-Position-plots: 1) ON 2) OFF 3) diff. ON-OFF(norm.)

  16. Phi-COG-Pos-plots: Size800-10000, hadcut=0.1 -plots: 1) ON (red), OFF(blue) 2) MC (black) 3) diff.% ON-OFF(norm.) COG-plots: 1) ON 2) OFF 3) diff. ON-OFF(norm.) Source-Position-plots: 1) ON 2) OFF 3) diff. ON-OFF(norm.)

  17. Phi-COG-Pos-plots: En.80-150[GeV], hadcut=0.25 -plots: 1) ON (red), OFF(blue) 2) MC (black) 3) diff.% ON-OFF(norm.) COG-plots: 1) ON 2) OFF 3) diff. ON-OFF(norm.) Source-Position-plots: 1) ON 2) OFF 3) diff. ON-OFF(norm.)

  18. Phi-COG-Pos-plots: En.80-20000[GeV], hadcut=0.2 -plots: 1) ON (red), OFF(blue) 2) MC (black) 3) diff.% ON-OFF(norm.) COG-plots: 1) ON 2) OFF 3) diff. ON-OFF(norm.) Source-Position-plots: 1) ON 2) OFF 3) diff. ON-OFF(norm.)

  19. ON&OFF Hillas Parameters vs Runs

  20. Rec. En., True En., Size, Had. (MC) Hadcut=0.3 Hadcut=0.3 Hadcut=0.3 Hadcut=0.8 Hadcut=0.8 Hadcut=0.3

  21. Analysis threshold (MC) Hadcut=0.3, True En.<200 GeV, Rec. En.<200 GeV

  22. Alpha-plot: 80 < E[GeV] < 150, hadcut=0.25

  23. Alpha-plot: 150 < E[GeV] < 300, hadcut=0.15

  24. Alpha-plot: 300 < E[GeV] < 500, hadcut=0.1

  25. Alpha-plot: 500 < E[GeV] < 20000, hadcut=0.15

  26. Setting for the calculation of the Upper Limits Energy range: 80 GeV < E < 20000 GeV (~infinity)Hadronness cut: hadcut=0.2Cut efficiency (MC):  = 0.6 (for < 9, E > 80 GeV)Supposed power law spectrum index: -1.5

  27. Collection Area: hadcut=0.2

  28. Alpha-plot: 80 < E[GeV] < 20000, hadcut=0.2 Total Entries: ON 76834 , OFF 26537 ==> ratio = 2.89Normalization factor = Ratio(ON/OFF)>25 = 2.91Entries inside signal region (< 9): ON 8989, OFF(norm.) 9223.65 Excess Upper Limit from AlphaPlot: Rolke u.l. (30%) = 135.9 [evts] Live ON time = 19449 sec = 5.4 h

  29. Flux Upper Limit from Willman-1 (5.4 h): F(E > 80 GeV) = 1.16079e-11 [cm-2 s-1]Projected Flux Upper Limit in 30 h: F(E > 80GeV) = 2.3218e-12 [cm-2 s-1] Projected Flux Upper Limit in 50 h: F(E > 80GeV) = 1.39031e-12 [cm-2 s-1]

  30. Same analysis performed on the only available CrabNebula (ON-OFF) sample with similar data taken conditions:Period: 02-03/2007 (MUX data)ON ~ 5.8 h OFF ~ 2.1 h , ratio ~ 2.8Same mean values of DTs and RatePSF 13mm , 8<Zd<30 :MC selected according to these values

  31. Alpha-plot: 80 < E[GeV] < 150, hadcut=0.25 CRAB

  32. Alpha-plot: 80 < E[GeV] < 150, hadcut=0.25 WILLMAN-1

  33. Alpha-plot: 80 < E[GeV] < 20000, hadcut=0.2 CRAB

  34. Alpha-plot: 80 < E[GeV] < 20000, hadcut=0.2 WILLMAN-1

  35. Alpha-plot: 250 < E[GeV] < 20000, hadcut=0.1 CRAB : max significance ---> ~ 15 /sqrt(h), Sensitivity: 1.6% CRAB

  36. Q-factors from MC and OFF (Alpha<12): 2D

  37. Q-factors from MC and OFF (Alpha<12): 3D

More Related