ProtoDUNE Dual Phase light data analysis ! e t a d p u Ana - - PowerPoint PPT Presentation

▶
protodune dual phase light data analysis
SMART_READER_LITE
LIVE PREVIEW

ProtoDUNE Dual Phase light data analysis ! e t a d p u Ana - - PowerPoint PPT Presentation

Oct 31, 2022 374 likes •544 views

ProtoDUNE Dual Phase light data analysis ! e t a d p u Ana Gallego Ros DPPD mee/ng, 26-11-2019 Contents 2 S1 signal analysis Background analysis Update on the PMT calibration status next meeting! S1 signal analysis 1)

slide-1
SLIDE 1

ProtoDUNE Dual Phase light data analysis

Ana Gallego Ros DPPD mee/ng, 26-11-2019

u p d a t e !

slide-2
SLIDE 2

✦ S1 signal analysis ✦ Background analysis ✦ Update on the PMT calibration status… next meeting!

Contents

2

slide-3
SLIDE 3

S1 signal analysis

1) Rate vs amplitude 2) Charge vs amplitude

slide-4
SLIDE 4

4

S1 signal: 1) rate vs amplitude

✦ Goal: obtain the rate of muons (S1 signals) crossing the detector with no fields ✦ Next steps:

  • S1 rate vs different field conditions
  • Run with TPB response ~ PEN response
  • **Long run (6-7/11/2019): S1 rate stability**

✦ Run #1625 (18/10/2019)

  • Random trigger
  • Fields off
  • 30k events (250 subruns)
  • Time window: 1 ms
  • Sampling: 16 ns
  • PMTs at G=1e7 (20190912)
slide-5
SLIDE 5

5

S1 signal: 1) rate vs amplitude

example

5 ADC 7-10 ADC 10 ADC

SPE amplitude (G = 1.5e7)

✦ S1 identification algorithm:

  • pedestal = 20 first samples with

ped_RMS<3 ADC

  • amplitude = pedestal - ADC_value
  • S1 signal if amplitude > 20 ADC
  • If S1 is found, the 16 μs window

after is ignored!

Pedestal RMS (ADC) 0.5 1 1.5 2 2.5 3 Events 20 40 60 80 100 120 140 160 FA0104 FA0130 FA0132 FA0139

Peak counting:

Event 18 --> S1 signal! Time: 328 us; Amplitude: 47 ADC; Charge: 1.74e+08 (e) Event 18 --> S1 signal! Time: 463 us; Amplitude: 145 ADC; Charge: 4.64e+08 (e) Event 18 --> S1 signal! Time: 845 us; Amplitude: 292 ADC; Charge: 1.05e+09 (e) Event 18 --> S1 signal! Time: 942 us; Amplitude: 27 ADC; Charge: 7.66e+07 (e)

✔

reminder!

slide-6
SLIDE 6

2000 4000 6000 8000 10000 12000 14000 16000 18000 20000

S1 rate (Hz)

6

S1 signal: 1) rate vs amplitude

ZOOM

✦ S1 rate = #S1 signals / effective time ✦ effective time = (time window x #events) - (#S1 signals x 16 μs) ✦ Results (average ± STD DEV) for S1 signals > 20 ADC:

  • TPB PMTs: S1 rate = 18 ± 3 kHz
  • PEN PMTs: S1 rate = 4.9 ± 0.9 kHz

(last meeting: only 4 central PMTs!)

Amplitude (ADC) 500 1000 1500 2000 2500 3000 3500 4000 S1 rate (Hz) 1 10

2

10

3

10

4

10

5

10 FA0111 FA0150 FA0113 FA0133 FA0129 FA0116 FA0105 FA0112 FA0157 FA0139 FA0104 FA0107 FA0153 FA0151 FA0137 FA0134 Amplitude (ADC) 500 1000 1500 2000 2500 3000 3500 4000 S1 rate (Hz) 1 10

2

10

3

10

4

10

5

10

6

10 FA0146 FA0121 FA0115 FA0136 FC0005 FC0004 FA0155 FA0149 FA0122 FA0114 FA0148 FA0135 FA0110 FA0124 FA0132 FA0130 Amplitude (ADC) 500 1000 1500 2000 2500 3000 3500 4000 S1 rate (Hz) 1 10

2

10

3

10

4

10 FA0111 FA0150 FA0113 FA0133 FA0129 FA0116 FA0105 FA0112 FA0157 FA0139 FA0104 FA0107 FA0153 FA0151 FA0137 FA0134 FA0146 FA0121 FA0115 FA0136 FC0005 FC0004 FA0155 FA0149 FA0122 FA0114 FA0148 FA0135 FA0110 FA0124 FA0132 FA0130

> ! PMTs 1-16 PMTs 17-32 32 PMTs

new! Run #1625 (2340/30k ev)

FA0114 FA0114

slide-7
SLIDE 7

7

S1 signal: 2) charge vs amplitude

✦ Goals:

  • Obtain average charge deposited by muons (S1 signals) vs amplitude
  • Detect possible ADC and/or PMT saturations

✦ Next step: compare results from different integration windows ✦ Run #1625 (18/10/2019)

  • Random trigger
  • Fields off
  • 30k events (250 subruns)
  • Time window: 1 ms
  • Sampling: 16 ns
  • PMTs at G=1e7 (20190912)
slide-8
SLIDE 8

S1 signal: 2) charge vs amplitude

8

1 2 3 4 5 6 3 10 × FA0133 Charge (e) 5 10 15 20 25 30 9 10 × Amplitude (ADC) 0.5 1 1.5 2 2.5 3 3.5 4 3 10 × FA0133 1 2 3 4 5 3 10 × FA0150 Charge (e) 5 10 15 20 25 30 9 10 × Amplitude (ADC) 0.5 1 1.5 2 2.5 3 3.5 4 3 10 × FA0150 1 2 3 4 5 3 10 × FA0116 Charge (e) 5 10 15 20 25 30 9 10 × Amplitude (ADC) 0.5 1 1.5 2 2.5 3 3.5 4 3 10 × FA0116 1 2 3 4 5 6 3 10 × FA0129 Charge (e) 5 10 15 20 25 30 9 10 × Amplitude (ADC) 0.5 1 1.5 2 2.5 3 3.5 4 3 10 × FA0129 1 2 3 4 5 3 10 × FA0112 Charge (e) 5 10 15 20 25 30 9 10 × Amplitude (ADC) 0.5 1 1.5 2 2.5 3 3.5 4 3 10 × FA0112 1 2 3 4 5 6 3 10 × FA0105 Charge (e) 5 10 15 20 25 30 9 10 × Amplitude (ADC) 0.5 1 1.5 2 2.5 3 3.5 4 3 10 × FA0105 1 2 3 4 5 3 10 × FA0111 Charge (e) 5 10 15 20 25 30 9 10 × Amplitude (ADC) 0.5 1 1.5 2 2.5 3 3.5 4 3 10 × FA0111 1 2 3 4 5 6 3 10 × FA0113 Charge (e) 5 10 15 20 25 30 9 10 × Amplitude (ADC) 0.5 1 1.5 2 2.5 3 3.5 4 3 10 × FA0113 1 2 3 4 5 6 3 10 × FA0134 Charge (e) 5 10 15 20 25 30 9 10 × Amplitude (ADC) 0.5 1 1.5 2 2.5 3 3.5 4 3 10 × FA0134 1 2 3 4 5 6 7 8 3 10 × FA0137 Charge (e) 5 10 15 20 25 30 9 10 × Amplitude (ADC) 0.5 1 1.5 2 2.5 3 3.5 4 3 10 × FA0137 1 2 3 4 5 6 7 8 3 10 × FA0139 Charge (e) 5 10 15 20 25 30 9 10 × Amplitude (ADC) 0.5 1 1.5 2 2.5 3 3.5 4 3 10 × FA0139 1 2 3 4 5 6 3 10 × FA0151 Charge (e) 5 10 15 20 25 30 9 10 × Amplitude (ADC) 0.5 1 1.5 2 2.5 3 3.5 4 3 10 × FA0151 1 2 3 4 5 6 3 10 × FA0153 Charge (e) 5 10 15 20 25 30 9 10 × Amplitude (ADC) 0.5 1 1.5 2 2.5 3 3.5 4 3 10 × FA0153 1 2 3 4 5 6 7 3 10 × FA0157 Charge (e) 5 10 15 20 25 30 9 10 × Amplitude (ADC) 0.5 1 1.5 2 2.5 3 3.5 4 3 10 × FA0157 1 2 3 4 5 3 10 × FC0004 Charge (e) 5 10 15 20 25 30 9 10 × Amplitude (ADC) 0.5 1 1.5 2 2.5 3 3.5 4 3 10 × FC0004 0.5 1 1.5 2 2.5 3 10 × FC0005 Charge (e) 5 10 15 20 25 30 9 10 × Amplitude (ADC) 0.5 1 1.5 2 2.5 3 3.5 4 3 10 × FC0005 1 2 3 4 5 6 3 10 × FA0115 Charge (e) 5 10 15 20 25 30 9 10 × Amplitude (ADC) 0.5 1 1.5 2 2.5 3 3.5 4 3 10 × FA0115 1 2 3 4 5 3 10 × FA0121 Charge (e) 5 10 15 20 25 30 9 10 × Amplitude (ADC) 0.5 1 1.5 2 2.5 3 3.5 4 3 10 × FA0121 1 2 3 4 5 3 10 × FA0136 Charge (e) 5 10 15 20 25 30 9 10 × Amplitude (ADC) 0.5 1 1.5 2 2.5 3 3.5 4 3 10 × FA0136 1 2 3 4 5 3 10 × FA0146 Charge (e) 5 10 15 20 25 30 9 10 × Amplitude (ADC) 0.5 1 1.5 2 2.5 3 3.5 4 3 10 × FA0146 1 2 3 4 5 3 10 × FA0149 Charge (e) 5 10 15 20 25 30 9 10 × Amplitude (ADC) 0.5 1 1.5 2 2.5 3 3.5 4 3 10 × FA0149 1 2 3 4 5 6 3 10 × FA0155 Charge (e) 5 10 15 20 25 30 9 10 × Amplitude (ADC) 0.5 1 1.5 2 2.5 3 3.5 4 3 10 × FA0155 1 2 3 4 5 6 7 8 3 10 × FA0107 Charge (e) 5 10 15 20 25 30 9 10 × Amplitude (ADC) 0.5 1 1.5 2 2.5 3 3.5 4 3 10 × FA0107 1 2 3 4 5 6 7 8 3 10 × FA0110 Charge (e) 5 10 15 20 25 30 9 10 × Amplitude (ADC) 0.5 1 1.5 2 2.5 3 3.5 4 3 10 × FA0110 1 2 3 4 5 6 7 3 10 × FA0114 Charge (e) 5 10 15 20 25 30 9 10 × Amplitude (ADC) 0.5 1 1.5 2 2.5 3 3.5 4 3 10 × FA0114 1 2 3 4 5 6 3 10 × FA0122 Charge (e) 5 10 15 20 25 30 9 10 × Amplitude (ADC) 0.5 1 1.5 2 2.5 3 3.5 4 3 10 × FA0122 1 2 3 4 5 6 7 8 3 10 × FA0124 Charge (e) 5 10 15 20 25 30 9 10 × Amplitude (ADC) 0.5 1 1.5 2 2.5 3 3.5 4 3 10 × FA0124 1 2 3 4 5 6 3 10 × FA0130 Charge (e) 5 10 15 20 25 30 9 10 × Amplitude (ADC) 0.5 1 1.5 2 2.5 3 3.5 4 3 10 × FA0130 1 2 3 4 5 6 7 8 3 10 × FA0132 Charge (e) 5 10 15 20 25 30 9 10 × Amplitude (ADC) 0.5 1 1.5 2 2.5 3 3.5 4 3 10 × FA0132 1 2 3 4 5 3 10 × FA0135 Charge (e) 5 10 15 20 25 30 9 10 × Amplitude (ADC) 0.5 1 1.5 2 2.5 3 3.5 4 3 10 × FA0135 1 2 3 4 5 6 3 10 × FA0148 Charge (e) 5 10 15 20 25 30 9 10 × Amplitude (ADC) 0.5 1 1.5 2 2.5 3 3.5 4 3 10 × FA0148
  • sci
  • sci
  • sci
  • sci

TPB TPB TPB TPB TPB TPB

1 2 3 4 5 6 7 3 10 × FA0114 Charge (e) 5 10 15 20 25 30 9 10 × Amplitude (ADC) 0.5 1 1.5 2 2.5 3 3.5 4 3 10 × FA0114

TPB

1 2 3 4 5 3 10 × FA0149 Charge (e) 5 10 15 20 25 30 9 10 × Amplitude (ADC) 0.5 1 1.5 2 2.5 3 3.5 4 3 10 × FA0149

PEN

Integration window: 96 ns (-2bins,+4bins) = (-32ns,+64ns)

1 2 3 4 5 6 7 3 10 × FA0104 Charge (e) 5 10 15 20 25 30 9 10 × Amplitude (ADC) 0.5 1 1.5 2 2.5 3 3.5 4 3 10 × FA0104

* ”free riders” (José)

* * * *

Run #1625: 2340/30k events!

1 2 3 4 5 3 10 × FC0004 Charge (e) 0.02 0.04 0.06 0.08 0.1 0.12 12 10 × Amplitude (ADC) 0.5 1 1.5 2 2.5 3 3.5 4 3 10 × FC0004 1 2 3 4 5 3 10 × FC0004 Charge (e) 5 10 15 20 25 30 9 10 × Amplitude (ADC) 0.5 1 1.5 2 2.5 3 3.5 4 3 10 × FC0004

Integration window: 4 μs (-30bins,+220bins) = (-480ns,+3520ns)

✦ 3x1x1 comet effect seen mostly in PMT with NB?

3x1x1

slide-9
SLIDE 9

S1 signal: 3) more statistics?

9

1 2 3 4 5 6 7 8 3 10 ×

FA0139

Charge (e) 5 10 15 20 25 30 9 10 × Amplitude (ADC) 0.5 1 1.5 2 2.5 3 3.5 4 3 10 ×

FA0139

1 2 3 4 5 6 3 10 ×

FA0130

Charge (e) 5 10 15 20 25 30 9 10 × Amplitude (ADC) 0.5 1 1.5 2 2.5 3 3.5 4 3 10 ×

FA0130

1 2 3 4 5 6 7 8 3 10 ×

FA0132

Charge (e) 5 10 15 20 25 30 9 10 × Amplitude (ADC) 0.5 1 1.5 2 2.5 3 3.5 4 3 10 ×

FA0132

Run #1625

10 20 30 40 50 60 70 80 90 3 10 ×

FA0104

Charge (e) 5 10 15 20 25 30 9 10 × Amplitude (ADC) 0.5 1 1.5 2 2.5 3 3.5 4 3 10 ×

FA0104

10 20 30 40 50 60 70 80 3 10 ×

FA0130

Charge (e) 5 10 15 20 25 30 9 10 × Amplitude (ADC) 0.5 1 1.5 2 2.5 3 3.5 4 3 10 ×

FA0130

20 40 60 80 100 3 10 ×

FA0132

Charge (e) 5 10 15 20 25 30 9 10 × Amplitude (ADC) 0.5 1 1.5 2 2.5 3 3.5 4 3 10 ×

FA0132

20 40 60 80 100 3 10 ×

FA0139

Charge (e) 5 10 15 20 25 30 9 10 × Amplitude (ADC) 0.5 1 1.5 2 2.5 3 3.5 4 3 10 ×

FA0139

1 2 3 4 5 6 7 3 10 ×

FA0104

Charge (e) 5 10 15 20 25 30 9 10 × Amplitude (ADC) 0.5 1 1.5 2 2.5 3 3.5 4 3 10 ×

FA0104

2340 events 30k events

TPB TPB TPB TPB

✦ S1 rate (kHz) for S1 signals > 20 ADC:

2340 events 30k events FA0132 21.0 21.6 FA0139 20.3 21.0 FA0130 5.6 5.7 FA0104 6.0 6.0

✔

slide-10
SLIDE 10

Background analysis

1) SPE amplitude vs PMT gain 2) SPE rate

v e r y p r e l i m i n a r !

slide-11
SLIDE 11

11

Background: 1) SPE amplitude vs PMT gain

✦ Goal: determine the SPE amplitude for a certain PMT gain ✦ Runs used so far:

  • #1006-1013 (20190912) —> PMT calibration (G vs HV)
  • 20k events
  • Fields off
  • Time window: 1.2 us
  • Sampling: 16 ns

✦ Other runs that can be used:

  • #862-897 (20190903)
  • #1284-1291 (20191003)

✦ More PMTs!

slide-12
SLIDE 12

PMT Gain 0.05 0.1 0.15 0.2 0.25 0.3

9

10 × 50 100 150 200 250 300 350 400

SPE amplitude (ADC) - FA0150

2e-08) x G ± SPE [ADC] = (1.23e-06

12

Background: 1) SPE amplitude vs PMT gain

✦ PMT gain obtained from HV & calibration results (A, B) ✦ Linear fit: SPEamp = factor x G (blue) ✦ Result: SPEamp(1e7) = 7 ± 2 ADC (11 / 32 PMTs)

SN SPE amp. @ 1e7 Error (ADC) (ADC) FA0111 7 0.07 FA0113 7 0.09 FA0147 7 0.2 FA0150 12 0.2 FA0156 7 0.06 FA0133 7 0.1 FA0105 7 0.1 FA0112 7 0.3 FA0116 6 0.1 FA0129 7 0.2 FA0157 7 0.06

Runs #1006-1013

unzoom!

example

✦ Maximum amplitude per event (8 HVs)

slide-13
SLIDE 13

s) µ Time ( 210 215 220 225 230 ADC 3992 3994 3996 3998 4000 4002 4004

Run #1625_0001, FA0104, ADC ch #23 - WF #3

13

Background: 2) SPE rate

example!

✦ Goal: obtain the background rate in the detector (PMT dark current + background) ✦ “SPE signal” identification algorithm:

  • pedestal = 20 first samples with ped_RMS<3

ADC

  • amplitude = pedestal - ADC_value
  • “SPE signal” if amplitude < 20 ADC
  • If S1 is found, the 16 μs window after is

ignored!

Similar to S1 peak counting!

slide-14
SLIDE 14

SPE amplitude (ADC) 2 4 6 8 10 12 14 16 18 20 Events

3

10

4

10

5

10

6

10

7

10 FA0111 FA0150 FA0113 FA0133 FA0129 FA0116 FA0105 FA0112 SPE amplitude (ADC) 2 4 6 8 10 12 14 16 18 20 Events

3

10

4

10

5

10

6

10

7

10 FA0157 FA0139 FA0104 FA0107 FA0153 FA0151 FA0137 FA0134 SPE amplitude (ADC) 2 4 6 8 10 12 14 16 18 20 Events

3

10

4

10

5

10

6

10

7

10 FA0122 FA0114 FA0148 FA0135 FA0110 FA0124 FA0132 FA0130

14

Background: 2) SPE rate

PMTs 1-8 PMTs 9-16

SPE ~ 7 ADC SPE ~ 7 ADC Electronic noise ~ 1 ADC Electronic noise ~ 1 ADC

SPE amplitude (ADC) 2 4 6 8 10 12 14 16 18 20 Events

3

10

4

10

5

10

6

10

7

10 FA0146 FA0121 FA0115 FA0136 FC0005 FC0004 FA0155 FA0149

PMTs 17-24

SPE ~ 7 ADC Electronic noise ~ 1 ADC

✦ Baseline PMT dark current: “…at CT, DC rate increases up to 1.7 ± 0.3 kHz, and some PMTs reach up to 2.5 kHz.” JINST 13 (2018) T10006: arXiv: 1806.04571v3 
 ✦ Background sources (SPE level) in ProtoDUNE-DP: Ar-39, Ar-42…

PMTs 25-32

SPE ~ 7 ADC Electronic noise ~ 1 ADC

Run #1625 (2340/30k ev)

slide-15
SLIDE 15

Amplitude (ADC) 2 4 6 8 10 12 14 16 18 20 SPE rate (Hz)

3

10

4

10

5

10

6

10

7

10 FA0146 FA0121 FA0115 FA0136 FC0005 FC0004 FA0155 FA0149 FA0122 FA0114 FA0148 FA0135 FA0110 FA0124 FA0132 FA0130

15

Background: 2) SPE rate

Amplitude (ADC) 2 4 6 8 10 12 14 16 18 20 SPE rate (Hz)

3

10

4

10

5

10

6

10

7

10 FA0111 FA0150 FA0113 FA0133 FA0129 FA0116 FA0105 FA0112 FA0157 FA0139 FA0104 FA0107 FA0153 FA0151 FA0137 FA0134

1-16 PMTs

SN Threshold (ADC) Rate (Hz) Error rate (Hz) FA0112 6 73059 192 FA0112 7 47793 156 FA0112 8 29980 123

—> SPE rate @ 1e7 (7 ADC) ≈ rate(th@6ADC) - rate (th@8ADC) ≈ 85.1 ± 0.3 kHz

example!

17-32 PMTs > ! > !

Run #1625 (2340/30k ev)