willman 1 preliminary analysis data set 01 04 2008 saverio lombardi infn padova 2008 04 14 n.
Skip this Video
Loading SlideShow in 5 Seconds..
Willman-1 preliminary analysis (data set 01-04/2008) Saverio Lombardi, INFN - Padova 2008/04/14 PowerPoint Presentation
Download Presentation
Willman-1 preliminary analysis (data set 01-04/2008) Saverio Lombardi, INFN - Padova 2008/04/14

Loading in 2 Seconds...

  share
play fullscreen
1 / 15
Download Presentation

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

bunme
70 Views
Download Presentation

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

- - - - - - - - - - - - - - - - - - - - - - - - - - - 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]