Highlights from ARGO-YBJ G. Di Sciascio INFN Sez. Roma Tor Vergata - - PowerPoint PPT Presentation

highlights from argo ybj
SMART_READER_LITE
LIVE PREVIEW

Highlights from ARGO-YBJ G. Di Sciascio INFN Sez. Roma Tor Vergata - - PowerPoint PPT Presentation

Highlights from ARGO-YBJ G. Di Sciascio INFN Sez. Roma Tor Vergata On behalf of the ARGO-YBJ Collaboration Vulcano Workshop 2010 May 25, 2010 The ARGO-YBJ experiment An unconventional EAS-array exploiting the Longitude 90 31


slide-1
SLIDE 1

Highlights from ARGO-YBJ

  • G. Di Sciascio

INFN – Sez. Roma “Tor Vergata” On behalf of the ARGO-YBJ Collaboration Vulcano Workshop 2010 – May 25, 2010

slide-2
SLIDE 2

Vulcano Workshop 2010

  • G. Di Sciascio

2

Longitude 90° 31’ 50” East Latitude 30° 06’ 38” North 90 Km North from Lhasa (Tibet) An unconventional EAS-array exploiting the full coverage approach at very high altitude, with the aim of studying

The ARGO-YBJ experiment

Tibet ASγ ARGO

The Yangbajing Cosmic Ray Laboratory

 VHE γ-Ray Astronomy  Gamma Ray Burst Physics  Cosmic Ray Physics

4300 m above the sea level

slide-3
SLIDE 3

The ARGO-YBJ experiment

Vulcano Workshop 2010

  • G. Di Sciascio

3

INFN and Dpt. di Fisica Università, Lecce INFN and Dpt. di Fisica Universita’, Napoli INFN and Dpt. di Fisica Universita’, Pavia INFN and Dpt di Fisica Università “Roma Tre”, Roma INFN and Dpt. di Fisica Univesità “Tor Vergata”, Roma INAF/IFSI and INFN, Torino INAF/IASF, Palermo and INFN, Catania IHEP, Beijing Shandong University, Jinan South West Jiaotong University, Chengdu Tibet University, Lhasa Yunnan University, Kunming ZhengZhou University, ZhengZhou Hong Kong University, Hong Kong

International Collaboration:

 Chinese Academy of Science (CAS)  Istituto Nazionale di Fisica Nucleare (INFN)

slide-4
SLIDE 4

Experimental Hall & Detector Layout

Vulcano Workshop 2010

  • G. Di Sciascio

4

Single layer of Resistive Plate Chambers (RPCs) with a full coverage (92% active surface) of a large area (5600 m2) + sampling guard ring (6700 m2 in total)

time resolution ~1-2 ns (pad) space resolution = strip

10 Pads (56 x 62 cm2) for each RPC 8 Strips (6.5 x 62 cm2) for each Pad 1 CLUSTER = 12 RPCs

78 m 111 m 99 m 74 m (5.7 7.6 m2)

Gas Mixture: Ar/ Iso/TFE = 15/10/75 HV = 7200 V Central Carpet:

130 Clusters 1560 RPCs 124800 Strips

slide-5
SLIDE 5

Vulcano Workshop 2010

  • G. Di Sciascio

5

Operational Modes

Object:

  • flaring phenomena (high energy tail of GRBs, solar flares)
  • detector and environment monitor

Recording the counting rates (Nhit ≥1, ≥2, ≥3, ≥4) for each cluster at fixed time intervals (every 0.5 s) lowers the energy threshold down to ≈ 1 GeV. No information on the arrival direction and spatial distribution of the detected particles.

  • Scaler Mode:

Detection of Extensive Air Showers (direction, size, core …) Coincidence of different detector units (pads) within 420 ns Trigger : ≥ 20 fired pads on the central carpet (rate ~3.6 kHz)

Object:

  • Cosmic Ray physics (above ~1 TeV)
  • VHE γ-astronomy (above ~300 GeV)
  • Shower Mode:

INDIPENDENT DAQ

slide-6
SLIDE 6

Number of fired Strips

Shower mode:

Space pixel: 7 62 cm2 (single strip) Time pixel: 56 62 cm2 (8 ORed strips = 1 Pad) Time resolution: ≈ 1 ns The size of pixels, the time resolution and the full coverage allow reconstruction with unprecedent details

slide-7
SLIDE 7

Vulcano Workshop 2010

  • G. Di Sciascio

7

  • Stable data taking since Nov. 2007 with full detector
  • Average duty cycle ~ 90%
  • Trigger rate ~3.6 kHz @ 20 pad threshold
  • Dead time 4%
  • 220 GB/day transferred to IHEP/CNAF data centers
slide-8
SLIDE 8

Vulcano Workshop 2010

  • G. Di Sciascio

8

Detector performance

Moon Shadow & Angular Resolution

slide-9
SLIDE 9

Vulcano Workshop 2010

  • G. Di Sciascio

9

The Moon Shadow

Geomagnetic Field: positively charged particles are deflected towards the West. Ion spectrometer The observation of the Moon shadow can provide a direct check

  • f the relation between size and primary energy

Cosmic rays are hampered by the Moon Deficit of cosmic rays in the direction of the Moon

) ( 6 . 1 TeV E Z

Moon diameter ~0.5 deg

 Size of the deficit  Position of the deficit Angular Resolution Pointing Error Energy calibration  West displacement

slide-10
SLIDE 10

Vulcano Workshop 2010

  • G. Di Sciascio

10

All data: 2006 → 2009

N > 100

55 s.d. 9 standard deviations / month PSF of the detector 3200 hours on-source

θ < 50°

The deficit surface is the convolution of the PSF

  • f the detector and the widespread Moon disc.
slide-11
SLIDE 11

Vulcano Workshop 2010

  • G. Di Sciascio

11

Moon Shadow analysis

Measured angular resolution Measured EW displacement

slide-12
SLIDE 12

Sun shadow

Vulcano Workshop 2010

  • G. Di Sciascio

12

Displacement of the Sun shadow correlates with the SMMF

The displacement of the Sun shadow is a good measurement of the IMF, especially in a this particular quiet phase between 23th and 24th cycles. Sun at maximum → shadow is washed out Sun at minimum → good shadow & SMF symmetric between NS EW shift due to GMF NS shift due to IMF

slide-13
SLIDE 13

Vulcano Workshop 2010

  • G. Di Sciascio

13

Cosmic Ray Physics

  • Antiproton/proton ratio measurement
  • Light-component spectrum of primary CRs
slide-14
SLIDE 14

p/p ratio at TeV energies

Vulcano Workshop 2010

  • G. Di Sciascio

14

Using data on Moon shadow, limits on antiparticle flux can be derived. Protons are deflected towards West, antiprotons are deflected towards East

→ 2 symmetric shadows expected.

If the displacement is large and the angular resolution small enough we can distinguish between the 2 shadows. If no event deficit on the antimatter side is observed an upper limit on antiproton content can be calculated.

slide-15
SLIDE 15

Upper limit on p/p by ARGO-YBJ

Vulcano Workshop 2010

  • G. Di Sciascio

15

Preliminary

slide-16
SLIDE 16

Light-component spectrum of CRs

Vulcano Workshop 2010

  • G. Di Sciascio

16

Measurement of the light-component (p+He) spectrum of primary CRs in the energy region (5 – 250) TeV via a Bayesian unfolding procedure CNO < 2%

ARGO data agree with CREAM results

Evidence that the proton spectrum is flatter than in the lower energy region CREAM p+He EAS-TOP + MACRO Horandel p+He CREAM p CREAM He

ARGO preliminary

slide-17
SLIDE 17

Vulcano Workshop 2010

  • G. Di Sciascio

17

Anisotropy & γ-ray astronomy

  • Large & Medium Scale Anisotropy
  • Search for point like sources
  • AGN Follow Up – Flaring
slide-18
SLIDE 18

Vulcano Workshop 2010

  • G. Di Sciascio

18

Cosmic rays excess 0.06% 0.1% N PAD > 40

584 days: 2007 Dec. – 2009 Nov.

Proton median energy 2 TeV

Smoothing radius = 5

Intermediate scale anisotropy

9 ·1010 events

r.a.=0 r.a.=360

slide-19
SLIDE 19

Vulcano Workshop 2010

  • G. Di Sciascio

19

Proton median energy 2 TeV

ARGO-YBJ

He Heliotail tail Ge Gemi minga Ga Galacti tic Pl Plane

~6 ·10-4 ~4 ·10-4 Proton median energy 10 TeV

MILAGRO

Multiple explanations were proposed:

Salvati & Sacco, A&A 485 (2008) 527 Drury & Aharonian, Astrop. Phys. 29 (2008) 420.

  • K. Munakata ,AIP Conf Proc Vol 932, page 283

Salvati, A&A 513 (2010) A28

slide-20
SLIDE 20

Large scale anisotropy

Tail-in Loss-cone Cygnus region

ARGO-YBJ DATA: 2008 and 2009

Vulcano Workshop 2010 20

  • G. Di Sciascio
slide-21
SLIDE 21

Vulcano Workshop 2010

  • G. Di Sciascio

21

Tibet AS

  • M. Amenomori et.al. Science, 2006
slide-22
SLIDE 22

Fit function:

Agree with diffusion model: larger amplitude for higher energy.

1D anisotropy for different energies

  • G. Guillian et. al. 2007 PRD

0.7 TeV 1.5 TeV 3.9 TeV

Vulcano Workshop 2010 22

  • G. Di Sciascio
slide-23
SLIDE 23

All sky survey result

Vulcano Workshop 2010

  • G. Di Sciascio

23

  • 3 sources with significance >5σ in 825 days (Jul. 06 – Oct. 09)
  • Crab 14.5 σ, Mrk421 12 σ, MGRO1908+06 5.4 σ

Crab Mrk421

Mean = -9.3 2.1 10-3 Sigma = 1.008 0.002

slide-24
SLIDE 24

Vulcano Workshop 2010

  • G. Di Sciascio

24

dN/dE = (3.73 0.80) ·10-11 · E –2.67 0.25 ev cm –2 s –1 TeV –1

N PAD Events /day Emed (TeV) 40 – 100 128 24 0.85 100 – 300 17.9 6.3 1.8 > 300 9.2 2.3 5.2

Crab Nebula

~14.5 s.d. in ~800 days

~50 % Crab/year

NO selection NO γ/h discrimination Absolute measurement

slide-25
SLIDE 25

Mrk421

Vulcano Workshop 2010

  • G. Di Sciascio

25

Total significance: 12 s.d. in ~800 days

Different flares observed by ARGO-YBJ in TeV range

Full DAQ TEST data

Mrk421 is characterized by a strong flaring activity both in X-rays and in TeV γ–rays. Swift (15-50 keV)

slide-26
SLIDE 26

Vulcano Workshop 2010

  • G. Di Sciascio

26

Mrk421 ASM / RXTE

Mrk 421: July 2006 flare

ARGO started recording data with the full central carpet during the X-ray flare of Mrk421 in July 2006  Commissioning Phase Evidence for TeV emission in coincidence with a X-ray flare

Npad > 60 Eγ

50 ~1.6 TeV, Eγ mode ~500 GeV

days 187-245 (110 hours) R.a. (deg)

≈6σ

δ (deg)

Flux 3-4 Crab

slide-27
SLIDE 27

Vulcano Workshop 2010

  • G. Di Sciascio

27

NPAD > 100

ARGO ASM/RXTE days from 1-1-08

10 days average

Mrk421: 2008 emission

Correlation coefficient = 0.64

Gamma Rays vs X-rays

slide-28
SLIDE 28

Vulcano Workshop 2010

  • G. Di Sciascio

28

Power law spectrum + EBL absorption :

dN/dE = (7.5 1.7) · 10-11 E –2.51 0.29 e- (E) ev cm –2 s –1 TeV –1

EBL from Primack et al. AIP conf. Proc. 745, 23, 2005

Integral flux (E > 1 TeV) 4.9 2.0 10-11 ev cm –2 s –1

2 Crab

Mrk421 energy spectrum

2008 days 41 – 180: when the source was in active state

slide-29
SLIDE 29

Vulcano Workshop 2010

  • G. Di Sciascio

29

GASP- WEBT SWIFT RXTE SWIFT SuperAGILE AGILE MAGIC VERITAS

Mrk421 - June 2008 flare

A set of simultaneous measurements covering 12 decades from optical to TeV energies was performed during the flaring activity in first half of June 2008 data from:

 GASP-WEBT (R-band; May 24 to June 23)  SWIFT (UVOT & XRT; June 12-13)  AGILE (E > 100 MeV; June 9-15)  MAGIC and VERITAS (E> 400 GeV; May 27 - June 8)

in Donnarumma et al., ApJ 691 L13 (2009) complemented by publicy-available data from RossiXTE/ASM (2-12 keV) and Swift/BAT (15-50 keV).

No VHE data after June 8

the moonlight hampered the Cherenkov telescopes measurements ! 2 flaring episodes were reported: June 3-8 and June 9-15

slide-30
SLIDE 30

Vulcano Workshop 2010

  • G. Di Sciascio

30

Mrk421 - June 2008 flare

Expected from theoretical SED

~18 h on source

NPAD > 100 3 days average June 5-7, 3.0 June 11-13, 4.2 1 day average

slide-31
SLIDE 31

Mrk421 - June 2008 flare

Vulcano Workshop 2010

  • G. Di Sciascio

31

dN/dE = (3.2 1.0) · 10-11 (E/2.5) –2.1 0.7 e- (E) ev cm –2 s –1 TeV –1

EBL from Raue, M. & Mazin, D., 2008, Int. J. Mod. Phys. D17, 1515

The spectrum slope is consistent with that measured by Whipple in 2000/2001

  • bserving a similar flare, suggesting that

the relation spectral index – flux is a long term property of the source. The integral flux above 1 TeV during June 11-13 is ~6 Crab units:

  • ne of the most powerful flare ever
  • bserved by Mrk421.
  • G. Aielli et al. – ApJL 714 (2010) L208
slide-32
SLIDE 32
  • Feb. 2010: strong Mrk421 flare

Vulcano Workshop 2010

  • G. Di Sciascio

32

 VERITAS reported similar

  • bservation in Atel #2443.

 The flux exceeded 3 Crab units for the duration of the

  • bservations. The peak flux

(16 Feb.) is >10 Crab units.  The flux measured by ARGO is consistent with VERITAS evaluation.  For the first time an EAS- array observed a TeV flare at 5σ on a daily basis. 16-18 Feb. 16 Feb. 17 Feb. 18 Feb.  ARGO observed a strong TeV γ-ray flare from Mrk421 during 16-18 Feb. 2010 at 6σ.

slide-33
SLIDE 33

MGRO1908+06

Vulcano Workshop 2010

  • G. Di Sciascio

33

  • Pulsar Wind Nebula discovered by Milagro,

confirmed by HESS and VERITAS.

  • HESS: intrinsic extension is 0.34 deg, spectrum

index is -2.1 up to 20TeV without cutoff.

  • Milagro: extension is <2.6 deg, spectrum cut-off

at about 14 TeV and its flux is higher than HESS result.

  • ARGO preliminary: extension ~0.6-1 deg. Our

data seems to support the Milagro spectrum.

slide-34
SLIDE 34

Cygnus region

Vulcano Workshop 2010

  • G. Di Sciascio

34

VERITAS

Milagro Fermi LAT

ARGO-YBJ

Tibet ASγ R=1.5 deg

No detection at 5σ, but with 2 years data ARGO is observing some signals from 2 TeV sources.

slide-35
SLIDE 35

Conclusions

Vulcano Workshop 2010

  • G. Di Sciascio

35

ARGO-YBJ detector taking data since November 2007 with excellent performance First results on γ-ray astronomy

Crab Nebula spectrum in agreement with other measurements Continuous monitor of Markarian 421 All sky survey Limits on 1-100 GeV fluence from GRBs

Cosmic Rays & Particle Physics

Large & Medium scale anisotropies Solar physics p-p cross section Limit on antiproton flux Light component spectrum measurement Multicore EAS with large transverse momentum

slide-36
SLIDE 36

Vulcano Workshop 2010

  • G. Di Sciascio

36

slide-37
SLIDE 37

IMF measurement with ARGO-YBJ

Vulcano Workshop 2010

  • G. Di Sciascio

37

slide-38
SLIDE 38
  • Long duration GRBs (>2s): 73
  • Short duration GRBs (≤2s): 10

Scaler mode

  • Number of GRBs analysed: 83
  • With known redshift: 14

Vulcano Workshop 2010

  • G. Di Sciascio

38

GRB in the ARGO FOV since Dec 2004 to Sep 2009

No evidence of coincident signal over the GRB T90 duration In stacked analysis no evidence for any integral effect

slide-39
SLIDE 39

GRBs with known redshift

Upper Limits in the 1-100 GeV Energy Range

Fluence upper limits (99% c.l.) obtained with differential spectral indexes ranging from the value measured by satellites to 2.5 .

Vulcano Workshop 2010

  • G. Di Sciascio

39

slide-40
SLIDE 40

Proton-Air cross section

Idea: measurement of shower rate at fixed energy as a function of the zenith angle

1 sec

) ( ) (

  • h

e I I

Where is not the p interaction lenght mainly because of inelasticity, shower fluctuations and detector resolution. = k

int , where k is determined by MC and

depends on:

  • hadronic interactions
  • detector features and location (atm. depth)
  • actual set of experimental observables
  • analysis cuts
  • energy threshold etc….

int p-Air (mb)=2.4 104 / int(g/cm2)

h0 q Data Full MC

Vulcano Workshop 2010

  • G. Di Sciascio

40

slide-41
SLIDE 41
  • Phys. Review D80 (2009) 092004.

Extending the energy range with the analog readout Vulcano Workshop 2010

  • G. Di Sciascio

41

slide-42
SLIDE 42

Vulcano Workshop 2010

  • G. Di Sciascio

42

slide-43
SLIDE 43

Vulcano Workshop 2010

  • G. Di Sciascio

43

Moon shadow: East-West projection

The higher the energy the lower the West displacement induced by the Earth magnetic field.