1 / 15

Willman-1 preliminary analysis (data set 01-04/2008) Saverio Lombardi, INFN - Padova 2008/04/14

Willman-1 preliminary analysis (data set 01-04/2008) Saverio Lombardi, INFN - Padova 2008/04/14. Data set 01-04/2008.

bunme
Download Presentation

Willman-1 preliminary analysis (data set 01-04/2008) Saverio Lombardi, INFN - Padova 2008/04/14

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-04/2008)Saverio Lombardi, INFN - Padova 2008/04/14

  2. Data set 01-04/2008 ON data:Total time: 16.9 hours, 22< Zd <30, PSF ~ 11mmMean rate before cuts ~ 200 HzMean rate after cuts ~ 140 Hz (no hadronness cut)Mean DT = 16OFF data:Total time: 9.34 hours, 22< Zd <26, PSF ~ 11mmMean rate before cuts ~ 200 HzMean rate after cuts ~ 140 Hz (no hadronness cut)Mean DT = 16

  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: AbsoluteApplied 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>340

  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 Est.Energy 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

  5. 80 < E [GeV] <20000, 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.)

  6. 80 < E [GeV] <20000, 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.)

  7. ON&OFF Hillas Parameters vs Runs More stringent data selection still to be applied!

  8. Analysis threshold (from MC): ~ 80 GeV Hadcut=0.3, True En.<200 GeV, Rec. En.<200 GeV

  9. Analysis threshold (from CRAB-test): ~ 100 GeV Hadcut=0.2 , min-E [GeV] = x , max-E [GeV] = 20000 Normalization region: 30 <  < 80, Signal region:  < 9 x = 70 ==> (LIMA) = 19.81, Exc./Sqrt(Bkg) = 43.89, norm. factor = 2.83 x = 80 ==> (LIMA) = 19.22 , Exc./Sqrt(Bkg) = 42.44, norm. factor = 2.88 x = 90 ==> (LIMA) = 20.48 , Exc./Sqrt(Bkg) = 45.89, norm. factor = 2.83 x = 100 ==> (LIMA) = 21.26 , Exc./Sqrt(Bkg) = 48.23, norm. factor = 2.82 x = 110 ==> (LIMA) = 20.93, Exc./Sqrt(Bkg) = 47.97, norm. factor = 2.84 Hadcut=0.2 , min-E [GeV] = x , max-E [GeV] = 120 Normalization region: 30 <  < 80, Signal region:  < 18 x = 70 ==> (LIMA) = 5.18, Exc./Sqrt(Bkg) = 10.61, norm. factor = 2.96 x = 80 ==> (LIMA) = 5.36, Exc./Sqrt(Bkg) = 10.92, norm. factor = 2.85 x = 90 ==> (LIMA) = 5.21, Exc./Sqrt(Bkg) = 10.66, norm. factor = 2.85 x = 100 ==> (LIMA) = 5.74, Exc./Sqrt(Bkg) = 11.86, norm. factor = 2.82 x = 110 ==> (LIMA) = 3.35, Exc./Sqrt(Bkg) = 6.96, norm. factor = 2.94

  10. CRAB-test: 90 < E [GeV] < 110, hadcut=0.2 CRAB-test: 100 < E [GeV] < 110, hadcut=0.2

  11. Upper Limits calculation setting Energy range: 100 GeV < E < 20000 GeV Hadronness cut: hadcut=0.2Cut efficiency (MC):  = 0.65 (for < 9, E > 100 GeV)Supposed power law spectrum index: -1.5

  12. Willman-1: 100 < E [GeV] < 20000, hadcut=0.2 Normalization factor = Ratio(ON/OFF)30 <  < 80 = 1.86Entries inside signal region (< 9): ON 21378 , OFF(normalized) 21688.9  231.8 (  sqrt(21688.9) )Excess Upper Limit from AlphaPlot: Rolke u.l. (30%) = 356.8 [evts]Live ON time = 60877.6 seconds = 16.91 hours (TimeON/TimeOFF = 1.81)

  13. CRAB-test: 100 < E [GeV] < 20000, hadcut=0.2 Normalization factor = Ratio(ON/OFF)30<  <80 = 2.82Entries inside signal region (< 9): ON 9704 , OFF(normalized) 5975.0  147.7(LIMA)/sqrt(h) = 8.92 , Sensitivity (50h) ~ 3.5 (% CRAB) Live ON time = 20464.3 seconds = 5.68 hours (TimeON/TimeOFF = 2.58)

  14. Collection Area: hadcut=0.2 Integral(100 < E [GeV] < 20000) = 168.754 km2

  15. Flux Upper Limit from Willman-1 (ON: 16.91 hours, OFF: 9.34 hours) F(E > 100 GeV) = 7.39271e-12 [cm-2 s-1]

More Related