Calorimetry (high precision @ LSD) FroST16 workshop FNAL, Chicago - - PowerPoint PPT Presentation

calorimetry
SMART_READER_LITE
LIVE PREVIEW

Calorimetry (high precision @ LSD) FroST16 workshop FNAL, Chicago - - PowerPoint PPT Presentation

Calorimetry (high precision @ LSD) FroST16 workshop FNAL, Chicago USA March 2016 Marco GRASSI (IHEP , China) Anatael Cabrera CNRS / IN2P3 APC Laboratory (Paris, FR) LNCA Underground Laboratory (Chooz, FR) todays LSD detectors 2


slide-1
SLIDE 1

Anatael Cabrera

CNRS / IN2P3 APC Laboratory (Paris, FR) LNCA Underground Laboratory (Chooz, FR)

Calorimetry

(high precision @ LSD)

FroST16 workshop FNAL, Chicago USA — March 2016

Marco GRASSI (IHEP , China)

slide-2
SLIDE 2 Anatael Cabrera (CNRS-IN2P3 & APC)

today’s LSD detectors…

2

today neutrino detectors (historical ordering)…

  • Liquid Scintillator [LSD] (ex. Double Chooz, KamLAND,etc)
  • Water Cherenkov (ex.Super-KamiokaNDE, etc)
  • Liquid-Argon TPC (ex. ICARUS, etc)

Liquid Scintillator detector (LSD) features… ✓signal ✓energy (→excellent) ✗ background: no event-by-event ID

  • radioactivity purity
↳ long & very expensive process
  • shielding & underground
↳ even up to ~50% cost (→ even more, if new underground lab needed)

✗ no doping or little (few % or ‰) [→limited physics programme] ✓cost (driven by PMT) energy measurement irreducible BG characterisation and/or subtraction

  • the name of game is more light & resolution (not just light)
  • low energies range (<MeV) and/or better calorimetry

⇒ this talk: a 2 topics (first time presented) on high resolution calorimetry

LSD: transparent & “large” PMTs largest LSDs: KamLAND & SNO+

simple transparent* homogenous* large size* composition (# protons)

~1k ton LSD @ ~2km underground

slide-3
SLIDE 3 Anatael Cabrera (CNRS-IN2P3 & APC)

LSD calorimetry in action…

3

Borexino (solar ν’s) KamLAND-ZEN phase-1 (2β0ν signal) Double Chooz (reactor ν’s) [DC-IV: θ13 & 8Li+9He measured at once] LSD signal⊕BGs rate+shape simultaneous measurement

  • thanks to energy resolution & linearity (i.e. calorimetry)
  • irreducible BGs→ only way to handle (if lucky)
slide-4
SLIDE 4 Anatael Cabrera (CNRS-IN2P3 & APC)

LSD’s PID in action (example PSD)…

4

LSD’s ability PSD ability (very limited→ little pattern recognition) depends light level

  • even if calorimetric is poor (due to dominant non-stochastic terms)
  • better if high precision calorimetry is possible→ likely better PID
slide-5
SLIDE 5 Anatael Cabrera (CNRS-IN2P3 & APC)

electronics electronics γ γ

pe λatte≤25m λscat≥10m

a necessary but not sufficient condition…light!

5

an LSD logic

Energy→ γ(scint) → γ(ws) → pe’s → (charge,time)

Liquid Scintillator wavelength-shifter (→fluors) PMT electronics

some fraction light loss→ >1% resolution PMT detection efficiency (~30%)→ ≳2% resolution charge digitisation (bias, non-linearities, etc)→ ~3% resolution ~10,000γ/MeV→ ≲1% resolution

how about response systematics?

better PMTs (more PEs) double calorimetry (better systematics)

slide-6
SLIDE 6 Anatael Cabrera (CNRS-IN2P3 & APC)

better PMTs…

6

(more & higher quality light)

slide-7
SLIDE 7

11" ETEL PMT

ETEL Development of 11 inch PMT started in 2013 for LBNE

Barros et al (including Svoboda) @ arXiv:1512.06916v2

7

15 prototypes produced and tested at Penn, UC Davis and Drexel.

slide-8
SLIDE 8

trigger ETEL PMT ETEL-11-PMT SPE Charge Spectrum

8
  • envelop made with Schott

8250 glass tube instead of EU glass

  • discolored spot due to the

manufacturing process – corrected in later version

  • characterisation of initial

batch corrects for this

  • 12 dynode structure: 107 @

~1300V (nominal)

  • TTS: 1.8ns (RMS→ ?)

ETEL-11-PMT SPE Time Spectrum

slide-9
SLIDE 9

charge: 11-inch similar 12-inch timing (TTS): 11-inch worse 12-inch ⇒ lower operating voltage may play into this (~1300V as

  • pposed to ~2000V)
9
slide-10
SLIDE 10

Relative Efficiency

Quantum x Collection Efficiency per cm2 comparable to Hamamatsu 12-inch and 10-inch HQE PMTs

10
slide-11
SLIDE 11 Anatael Cabrera (CNRS-IN2P3 & APC)

(within JUNO) new 20” MCP based PMT…

11
  • 20” PMT (but also 8”) developed in China
  • increase in effective detection efficiency

(comparable to Hamamatsu20”) → increase of light level is highest priority → large peak-to-valley (SPE efficiency)

  • TTS ~10ns not best (location of MCP)

→ JUNO hybrid system with other PMTs @JUNO decision…

  • 15,000 MCP-PMT 20”
  • 5,000 Hamamatsu PMT 20”
slide-12
SLIDE 12 Anatael Cabrera (CNRS-IN2P3 & APC)

control of systematics…

12

(i.e. non-stochastic effects)

slide-13
SLIDE 13 Anatael Cabrera (CNRS-IN2P3 & APC)

JUNO location…

13

simplistic schedule: data-taking by 2020

slide-14
SLIDE 14 Anatael Cabrera (CNRS-IN2P3 & APC) 14

~18,000 PMTs (20” diameter)→ Large-PMT system (LPMT) ~36,000 PMTs (3” diameter)→ Small-PMT system (SPMT)

slide-15
SLIDE 15 Anatael Cabrera (CNRS-IN2P3 & APC) 15

largest photo-cathode density ever built ⇒ highest precision calorimetry ever built largest light level ever detected ~1200PE/MeV ⇒ stochastic resolution <3% @ 1MeV control of non-stochastic resolution extremely demanding→ ≲1% (driven by SPMT)

slide-16
SLIDE 16 Anatael Cabrera (CNRS-IN2P3 & APC) >15 laboratories so far… Brasil
  • FABC (Sao Paulo)
  • PUC (Rio de Janeiro)
Belgium
  • UBL (Brussels)
Chile
  • PUC (Santiago)
China
  • IHEP (Beijing)
  • SYSU (Guangzhou)
France
  • APC (Paris)(coordination)
  • CPPM (Marseille)
  • LLR (Paris)
  • OMEGA (Paris)
  • SUBATECH (Nantes)
Italy
  • Padova-INFN (Padova)
Taiwan
  • National Taiwan University NTU (Taipei)
  • National Chiao Tung University NCTU (Hsinchu)
  • National United University NUU (Miaoli)
A few more institutions joining…
  • ur (very international) team…
16
slide-17
SLIDE 17 Anatael Cabrera (CNRS-IN2P3 & APC)

motivation…

17

— why the SPMT? —

slide-18
SLIDE 18 Anatael Cabrera (CNRS-IN2P3 & APC)

Visible Energy (MeV) 1 2 3 4 5 6 7 8 9 Energy Resolution 0.00 0.02 0.04 0.06 0.08 0.10 0.12 0.14

Cs 137 Ge 68 Cf) 252 H ( Co 60 ) ν Cf, 252 Gd ( C (GC, spall. n) data MC volume source (data) volume source (MC) DC-III (Gd-n) Preliminary

DC with 1200PE/MeV

non-stochastic terms (i.e. b & c): very sensitive to high energy level arm (understood?)

DC: ~200PE/MeV RMS=0.35%

Elapsed Days 100 200 300 400 500 600 700 Variation 0.97 0.98 0.99 1.00 1.01 1.02 1.03 1.04 1.05 1.06 n-H captures BiPo212 n-Gd captures DC-III (Gd-n) Preliminary after stability calibration

(BiPo poor stats)

control of response uniformity ±1%

DC as prototype for JUNO…

18

control of response stability

slide-19
SLIDE 19 Anatael Cabrera (CNRS-IN2P3 & APC)

no perfect world…

19
  • if perfect light measurement: σ(E)2non-stoch→0 (i.e. LS⊕PMT⊕electronics no dispersive effects)
  • if perfect calibration: σ(E)2non-stoch→0 (i.e. perfect correction of dispersive effects)

(unfortunately) none is true!!

JUNO* JUNO* JUNO* JUNO* Visible Energy (MeV) JUNO* [1.2kPE/MeV only stochastic] JUNO* [non-stochastic: a la DC] JUNO* [non-stochastic: half DC] JUNO* [non-stochastic: “negligible”] σ(E)2 = σ(E)2stoch + σ(E)2non-stoch ⟹ empiric formulation:

(1200PE/MeV) (??%)

~1.2k PEs σ(E)stoch < 3% the impact of σ(E)non-stoch dominates!!

slide-20
SLIDE 20 Anatael Cabrera (CNRS-IN2P3 & APC)

the double calorimetry…

20

σ(E)2 = σ(E)2stoch + σ(E)2non-stoch (1200PE @ 1MeV) if σ(E)2≤3.0%⇒ σ(E)2stoch=2.89% & + σ(E)2non-stoch=0.82% (remaining) now consider (1200±50)PEs @ 1MeV (same condition as before)⇒

  • +50PEs implies σ(E)2stoch=2.83% & + σ(E)2non-stoch=1.00% (remaining)
  • -50PEs implies σ(E)2stoch=2.95% & + σ(E)2non-stoch=0.55% (remaining)

small difference in light level (>1150PE/MeV)⇒ major impact to σ(E)2non-stoch: most challenging!!

~2x @DC: σ(E)2non-stoch≳2%

≥1300PE/MeV (→σnon-stoch≥1.0%)

“double-calorimetry”

articulate 2 energy estimators (different behaviours) Energy(photon-counting) i.e. digital (PS) Energy(charge integration) i.e. digital (QI)

⇒ E(response,x,y,z)DC = E(PS)⊕E(QI)

[via NN, correction, etc] control/reduction σ(E)2non-stoch & redundancy [if ±Δm2→ convince JUNO can]

ρ position (mm)

response uniformity

Response (normalised @ ⊙)

Response(QI) Response(PS)

slide-21
SLIDE 21 Anatael Cabrera (CNRS-IN2P3 & APC)

the JUNO challenge…

21

HIGHEST precision calorimetry (≤3% @ 1MeV) ⊕ LARGEST dynamic range in calorimetry (channel-wise) [⇒ uniformity⊕linearity⊕stability] (λ⦿≈0.28)

mean illumination per channel (PE/PMT) if λ≲0.5⇒ ~photon-counting regime

KamLAND

1880PMTs ~250PE/MeV

(λ⦿≈0.35) (λ⦿≈1.0) (λ⦿≈0.13)

DB

190PMTs ~180PE/MeV

DC

390PMTs ~180PE/MeV

Bx

2212PMTs ~500PE/MeV

JUNO

17000PMTs ~1200PE/MeV

~2x ~3x

λ⦿ = mean illumination per channel @ center

@1MeV

~4.5m buffer ≤4x NT GC (λ⦿≈0.07) ~100x LPMT ≤4x SPMT

slide-22
SLIDE 22 Anatael Cabrera (CNRS-IN2P3 & APC)

Elapsed Days 100 200 300 400 500 600 α ∆

  • 0.05
  • 0.04
  • 0.03
  • 0.02
  • 0.01

0.00 0.01 0.02 0.03 0.04 0.05

Double Chooz Preliminary

PS vs QI in action…

22

DC-III (data)

“digital” response stability @ 2.2MeV (zero tracking⊕other effect) (invisible to charge integration estimator alone)

Energy(PC) & Energy(QI) are highly complementary!!

response stability

slide-23
SLIDE 23 Anatael Cabrera (CNRS-IN2P3 & APC)

Photon-Counting vs Charge-Integration…

23

Readout Window Readout Window

PE discrimination threshold (a fraction of a PE) RECO-INFO @PC…
  • 1hit=1PE→no reconstruction (extreme: no gain needed!)
  • time-stamp (HW→ TDC-like)
  • charge info (HW→ high dynamics ADC)
RECO-INFO @QI…
  • 1hit≠1PE→reconstruction is a must! [QI-reco]
  • time-stamp the digitised readout window ~[0.3,1.0]μs
  • time & charge upon SW-reco (on/off-line)
time-stamp (self-trigger) (pulse @ threshold→ discriminator) time-stamp (global trigger) (readout-window sampled, example @ start)

@LPMT @SPMT

slide-24
SLIDE 24 Anatael Cabrera (CNRS-IN2P3 & APC)

the SPMT & LPMT calorimetry regimes…

24

~50% statistics

PE Maximum @1MeV

LPMT has dramatic variation across volume (→ systematics and/or biasses) (wildest variation in region with large fraction of statistics) (opposite) SPMT has FLAT response across volume (by construction) (SPMT ideal input for Trigger) ≲2x ≲5x ~%

~25% statistics ≲3% statistics

slide-25
SLIDE 25 Anatael Cabrera (CNRS-IN2P3 & APC)

(illustration) response/channel vs position…

25

LPMT only

1PE [2,5]PE >5PE

PMT fraction Charge fraction

1PE [2,5]PE >5PE

Large PMTs can detect up to 100pe for an IBD event in the last shell (20% of events)

small bias in few LPMTs⇒ large impact to over calorimetry!

@center (≤4m) @edge (≥16m) IBD

slide-26
SLIDE 26 Anatael Cabrera (CNRS-IN2P3 & APC) Energy [MeV] 2.4 2.6 2.8 Co60 Hits 20" Photomultiplier True Hits Reco Hits Reco / True 0.99 0.995 1 Radius [m] 5 10 15 J15v1r2 pre1 10 20 30 40 50 60 70 80 True charge [p.e.] 1 2 3 4 5 6 7 8 9 10 Reconstructed / true charge 0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 1.3 1.4 1.5 Time [ns] 150 200 250 300 350 400 450 500 550 600 650 Charge [p.e.] 0.1 0.2 0.3 0.4 Integrated signal

energy reconstruction bias estimation (1)…

26

non-linearity (channel-wise) non-uniformity (position-wise) [QI regime variations] worsens resolution (full detector)

realistic pulse reco (QI)

non-linearity (QI) calibration mimicking 20%→5% (no gain bias)

slide-27
SLIDE 27 Anatael Cabrera (CNRS-IN2P3 & APC) Radius (m) 200 400 600 800 1000 1200 1400 1600 1800 Energy (MeV) 1 2 3 4 5 6 7 8 p2_rot_SPmt Entries 10000 Mean x 776.7 Mean y 3.547 RMS x 451.1 RMS y 2.027 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1 p2_rot_SPmt Entries 10000 Mean x 776.7 Mean y 3.547 RMS x 451.1 RMS y 2.027

linearity⊕uniformity crosstalk handling…

27

SPMT only

Radius (cm) 200 400 600 800 1000 1200 1400 1600 1800 Energy (MeV) 1 2 3 4 5 6 7 p2_rot_LPmt Entries 10000 Mean x 1164 Mean y 2.995 RMS x 302.2 RMS y 1.713 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1 p2_rot_LPmt Entries 10000 Mean x 1164 Mean y 2.995 RMS x 302.2 RMS y 1.713

LPMT only

if linearity⊕uniformity⇒ LPMT 3D-maps a must!

3 D m a p s ? n

  • n

e e d ( ≈ fl a t )

SPMT: uniformity map & linearity⇒ (independent) 3D-map validation (simpler, complementary & robust→ unique, if SPMT)

slide-28
SLIDE 28 Anatael Cabrera (CNRS-IN2P3 & APC)

(illustration) LPMT 3D calibration maps…

28

Radius (cm)

200 400 600 800 1000 1200 1400 1600 1800

Energy (MeV)

1 2 3 4 5 6 7 8 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1

LPMT 3D map (easy to say), but which source?

slide-29
SLIDE 29 Anatael Cabrera (CNRS-IN2P3 & APC)

response summary…

29

LPMT: uniformity • linearity • stability ≠ 0

(i.e. not orthogonal bias/systematics)

SPMT: uniformity • linearity • stability ≈ 0

(i.e. effective orthogonal bias/systematics)

vs

(far more knowledge when combining)

slide-30
SLIDE 30 Anatael Cabrera (CNRS-IN2P3 & APC)

JUNO upgrade…

30

JUNO (before) JUNO (now) double calorimetric single-calorimetric

slide-31
SLIDE 31 Anatael Cabrera (CNRS-IN2P3 & APC)

SPMT system: much more…

31
slide-32
SLIDE 32 Anatael Cabrera (CNRS-IN2P3 & APC)

SPMT: excellent μ-physics…

32
slide-33
SLIDE 33 Anatael Cabrera (CNRS-IN2P3 & APC)

improving multi-μ identification…?

33

μ: ≤300PE per SPMT (no saturation whatsoever) LPMT (no saturation) LPMT (saturation at 4000PE) SPMT

evidently so…

saturation model very complex (not uniform, no flat, etc)

slide-34
SLIDE 34 Anatael Cabrera (CNRS-IN2P3 & APC) 34

…less is more! (→SPMT) when dazzling… (i.e. saturation) when dealing with μ’s…

slide-35
SLIDE 35 Anatael Cabrera (CNRS-IN2P3 & APC)

SPMT as an “aider” to the LPMT…

35
  • A. high precision calorimetry response systematics IBD physics

(highest priority: aide ≤3% @ 1MeV resolution)

  • B. improve inner-detector μ-reconstruction resolution

(highest priority: aide 12B/9Li/8He tagging/vetoing)

  • C. high rate SN pile-up (if very near)

(medium priority: minimise bias in absolute rate & energy spectrum)

  • D. vital complementarity: time resolution, dynamic range & trigger

(articulate additional complementary to LPMT system: better/simpler)

how about neutrino physics?

slide-36
SLIDE 36 Anatael Cabrera (CNRS-IN2P3 & APC)

high precision (θ12,δm2) also with SPMT?

36

Energy Visible (MeV) Δm2 (i.e. period) sin2(2θ13) sin2(2θ12) “atmospheric”

  • scillations

“solar”

  • scillations

δm2 JUNO several δm2 (<1% precision)… (only 2 fully independent) (δm2)SPMT independent (digital calorimetry) (δm2)LPMT independent (integration calorimetry) (δm2)LPMT⊕SPMT independent (double calorimetry) use (δm2)SPMT to validate linearity (or bias) of (δm2)LPMT & (δm2)LPMT⊕SPMT (use solar disappearance to cross-calibrate calorimetry for Mass Ordering precision & accuracy)

slide-37
SLIDE 37 Anatael Cabrera (CNRS-IN2P3 & APC)

what to remember?

  • LSD calorimetry so far…
  • most valuable asset of LSD technology (→ vital even for BG reduction)
  • (so far) experiments ≲500PE/MeV (higher to reach lower energies: Bx→>100keV)
  • most experiments dominated by photon-statistics: non-stochastic ~few% @ 1MeV ~irrelevant
  • energy resolution is not a vital criteria for physics so far [most experiment don’t even ever show!!]
  • higher precision calorimetry…
  • <5% energy resolution has hardly been tried as criterial for physics (to my knowledge) [ex. JUNO]
  • no show-stoppers, but novelties are needed to handle non-stochastic terms: double-calorimetry
  • ~3% [current technology: DC, DYB, KamLAND, Bx, etc] likely impossible due to non-stochastic
  • beyond Order(3%) resolution is likely impossible(?) with PMT⊕LS(~10,000γ/MeV)→ stochastic term
  • leave alone the non-stochastic terms (already @ few ‰)…
  • beyond high precision calorimetry…
  • more is always better, as this implies more information [evident!]
  • impossible better calorimetry (non-stochastic terms)→ typical case of segmentation [dead-material]
  • more light use for PID, such as PSD, etc [strongly recommend for high BG environment]
37