Study of Accidental Activity at the Front Barrel of the KOTO - - PowerPoint PPT Presentation

▶
study of accidental activity at the front barrel of the
SMART_READER_LITE
LIVE PREVIEW

Study of Accidental Activity at the Front Barrel of the KOTO - - PowerPoint PPT Presentation

Study of Accidental Activity at the Front Barrel of the KOTO Detector Ryota Shiraishi Yamanaka Group Kuno-Yamanaka Group Year-End Presentation 2019 Ryota Shiraishi 2019.12.23 Kuno-Yamanaka Group Year-End Presentation 2019 Introduction


slide-1
SLIDE 1

Ryota Shiraishi 2019.12.23 Kuno-Yamanaka Group Year-End Presentation 2019

Study of Accidental Activity at the Front Barrel of the KOTO Detector

Ryota Shiraishi Yamanaka Group Kuno-Yamanaka Group Year-End Presentation 2019

slide-2
SLIDE 2

Ryota Shiraishi 2019.12.23 Kuno-Yamanaka Group Year-End Presentation 2019

Introduction

2

The KOTO experiment

Purpose : To observe the decay . Signal : 2 photons + “nothing”

KL → π0ν¯ ν

➥ detected at the CsI calorimeter ➥ other veto detectors make sure of no extra hits

γ γ π0 ν¯ ν

KL

CsI

slide-3
SLIDE 3

Ryota Shiraishi 2019.12.23 Kuno-Yamanaka Group Year-End Presentation 2019

Accidental Signal Loss

3

Accidental hits on veto detectors coincident with the decay could cause signal loss. Major sources…

  • Other KL decay
  • Neutron from the J-PARC primary beam line

π0

decay signal @CsI

π0

Accidental hit @veto detector

Discarded (Acceptance loss) Accidental hits need to be reduced.

Veto Window

slide-4
SLIDE 4

Ryota Shiraishi 2019.12.23 Kuno-Yamanaka Group Year-End Presentation 2019

Purpose of This Study

4

  • To understand accidental activities at the Front

Barrel of the KOTO detector.

  • To calculate accidental hit rates by using data

taken in 2019.

  • To check consistency of accidental hit rates

between physics-triggered data and TMON- triggerd data.

Physics trigger…trigger to collect data TMON trigger…trigger to reproduce accidental hits

KL → π0ν¯ ν

slide-5
SLIDE 5

Ryota Shiraishi 2019.12.23 Kuno-Yamanaka Group Year-End Presentation 2019

Front Barrel

5

View from the downstream side

  • Sandwich of lead & plastic

scintillators

  • 2.75m long
  • Sampled by125MHz FADC
  • 16 modules
  • 32 readout channels

(inner/outer layers are read separately)

Front Barrel

slide-6
SLIDE 6

Ryota Shiraishi 2019.12.23 Kuno-Yamanaka Group Year-End Presentation 2019

TMON Trigger

6

TMON trigger is…

A random trigger generated from the signals of the Target

  • Monitor. The rate is proportional to the beam intensity.

We use this trigger to reproduce accidental activities and

  • verlay the waveforms on generated waveforms in simulation.

Proton Beam

Au Target 50° 16° KOTO Beam Line Plastic Scintillators

TMON Trigger

slide-7
SLIDE 7

Ryota Shiraishi 2019.12.23 Kuno-Yamanaka Group Year-End Presentation 2019

Energy Distribution

7

TMON data Energy is distributed up to ~600 MeV. Now, set the energy threshold to 2 MeV. Consider only events with FBAREne > 2 MeV. ➤

FBARModID 5 10 15 20 25 30 FBAREne [MeV] 100 200 300 400 500 600 1 10

2

10

3

10

4

10

FBAREne:FBARModID {ExtTrigType==2}

FBARModID FBAREne [MeV]

h Entries 30082 Mean 2.337 RMS 4.122

5 10 15 20 25 30

2

10

3

10

4

10

h Entries 30082 Mean 2.337 RMS 4.122

FBAREne {FBARModID==0 && ExtTrigType==2 && FBAREne<50}

# of events

FBAREne {FBARModID==0 && ExtTrigType==2}

600

FBAREne:FBARModID {ExtTrigType==2}

104 103 102 2

FBAREne [MeV]

30

Higher counts in inner channels Lower counts in outer channels

➞

slide-8
SLIDE 8

10 10 20 30 40 50 60 70

Rate [Hz]

FBAR Rate ch0 (E>2.000000MeV)

100 200 300 400 500 600

3

10 ×

FBAR Rate ch0 (E>2.000000MeV)

FBARPTime [clock]

10 − 10 20 30 40 50 60 70

ratio (phys / tmon)

0.5 1 1.5

− : Physics − : TMON

FBARPTime [clock] Rate [Hz] Ratio (phys./tmon) FBAR Rate ch0 (E > 2 MeV)

600

×103

Ryota Shiraishi 2019.12.23 Kuno-Yamanaka Group Year-End Presentation 2019

Comparison between Phys. & TMON data

8

Rate[Hz] per bin(= 1clock=8ns)

Rate = #events / (∆T×#triggered) (∆T = 1clock = 8ns)

Earlier timing region

  • > consistent

Later timing region

  • > subtle discrepancy exists
slide-9
SLIDE 9

Ryota Shiraishi 2019.12.23 Kuno-Yamanaka Group Year-End Presentation 2019

Rates in 10 - 20 clock timing

9

Ratio = Phys./TMON = 1 ~ 1.1

Energy Threshold : 2MeV Rate = #events/(∆T×#triggered) ∆T = ∆T1

10 10 20 30 40 50 60 70

Rate [Hz]

FBAR Rate ch0 (E>2.000000MeV)

100 200 300 400 500 600

3

10 ×

FBAR Rate ch0 (E>2.000000MeV)

FBARPTime [clock]

10 − 10 20 30 40 50 60 70

ratio (phys / tmon)

0.5 1 1.5

∆T1 FBAR Rate ch0 (E>2 MeV)

FBARPTime [clock] Rate [Hz] Ratio (phys./tmon)

− : Physics − : TMON

➞ Good agreement

5 10 15 20 25 30

Rate [Hz]

50 100 150 200 250 300

3

10 ×

(FBARPTime>10.000000 && FBARPTime<20.000000) FBAR Rate E>2.000000MeV

FBARModID

5 10 15 20 25 30

ratio (phys / tmon)

0.9 1 1.1 1.2

E > 2MeV, 10 < FBARPTime < 20

FBARModID Rate [Hz] Ratio (phys./tmon) Outer Channels Inner Channels 300 ×103

⟶ Primary beam line side

slide-10
SLIDE 10

Ryota Shiraishi 2019.12.23 Kuno-Yamanaka Group Year-End Presentation 2019

How we will reduce accidental hits

10

Accidental activities by neutrons coming from the J-PARC primary beam line

KOTO detector

To reduce the neutron flux, we installed a 33cm-thick iron wall.

Iron Wall

KOTO Detector

neutron

Primary Beam Line Iron Wall

slide-11
SLIDE 11

Ryota Shiraishi 2019.12.23 Kuno-Yamanaka Group Year-End Presentation 2019

Summary / To do

11

  • Confirmed consistency of accidental counting

rates in physics/TMON data.

  • To reduce neutrons from the primary beam

line, we installed an iron wall and will check the reduction effect.

slide-12
SLIDE 12

Ryota Shiraishi 2019.12.23 Kuno-Yamanaka Group Year-End Presentation 2019

Backup

slide-13
SLIDE 13

Ryota Shiraishi 2019.12.23 Kuno-Yamanaka Group Year-End Presentation 2019

Detector

13

View from the downstream side

slide-14
SLIDE 14

Ryota Shiraishi 2019.12.23 Kuno-Yamanaka Group Year-End Presentation 2019

Waveforms at the Front Barrel

14

Waveform Examples

Physics-triggerd event Energy : 13.8 [MeV] Time : 17.7 [clock=8ns] TMON-triggerd event Energy : 3.1 [MeV] Time : 47.3 [clock=8ns]

Consider an energy threshold and timing distributions.

Iteration$ 10 20 30 40 50 60 70 FBAR.Data 500 600 700 800 900

FBAR.Data:Iteration$ {FBAR.ModID==0 && Entry$==256}

FBAR.Data:Iteration$ {FBAR.ModID==0 && Entry$==256} Time [clock=8ns] ADC counts 20

Iteration$ 10 20 30 40 50 60 70 FBAR.Data 485 490 495 500 505 510 515 520

FBAR.Data:Iteration$ {FBAR.ModID==0 && Entry$==706}

FBAR.Data:Iteration$ {FBAR.ModID==0 && Entry$==706} Time [clock=8ns] ADC counts 40

slide-15
SLIDE 15

Ryota Shiraishi 2019.12.23 Kuno-Yamanaka Group Year-End Presentation 2019

Time Distribution

15

Take the moving average -> make waveforms smoother Calculate a parabolic curve using three samples around the peak. The parabola time closer to the nominal time (~31 clocks) for the Front Barrel is selected. Tend to have a structure like a broad hill.

➤

htemp Entries 8044 Mean 28.49 RMS 16.34 FBARPTime 10 20 30 40 50 60 50 100 150 200 250 300 350 htemp Entries 8044 Mean 28.49 RMS 16.34

FBARPTime {FBARModID==0 && FBAREne>2 && FBARPTime>-10 && ExtTrigType==2}

FBARPTime [clock]

# of events

TMON

FABRPTime {FBARModID==0 && FBAREne>2 && FBARPTime>-10 && (ScaledTrigBit&0x1)==0x1} 30

htemp Entries 78920 Mean 28.95 RMS 15.4 FBARPTime 10 20 30 40 50 60 500 1000 1500 2000 2500 3000 3500 htemp Entries 78920 Mean 28.95 RMS 15.4

FBARPTime {FBARModID==0 && FBAREne>2 && FBARPTime>-10 && (ScaledTrigBit&0x1)==0x1}

FBARPTime [clock]

# of events

Physics

FABRPTime {FBARModID==0 && FBAREne>2 && FBARPTime>-10 && ExtTrigType==2} 30

Concentrated around the nominal time Dominated by accidental hits Parabola time Parabola curve

slide-16
SLIDE 16

Ryota Shiraishi 2019.12.23 Kuno-Yamanaka Group Year-End Presentation 2019

Moving Average

16

Calculation Code

Take the average of five consecutive samples. Make waveforms smoother and mitigate local fluctuations.

/sw/koto/e14ana/release/v4.01.10/AnalysisLibrary/UserProjects/ E14ProdLibrary/E14ProdDstConv/src/E14CrateData125MHz.cc

10 20 30 40 50 60 460 480 500 520 540 560 580 600

wfm2

Entries 64 Mean x 31.5 Mean y 503.8 Std Dev x 18.47 Std Dev y 16.24

wfm2

Entries 64 Mean x 31.5 Mean y 503.8 Std Dev x 18.47 Std Dev y 16.24

FBARWfm[0]:Iteration$ {Entry$==153}

10 20 30 40 50 60 460 480 500 520 540 560 580 600

wfm

Entries 60 Mean x 31.5 Mean y 503 Std Dev x 17.32 Std Dev y 15.12

wfm

Entries 60 Mean x 31.5 Mean y 503 Std Dev x 17.32 Std Dev y 15.12

wfm

Averaged Wfm Time [clock] Time [clock] ADC ADC

slide-17
SLIDE 17

Ryota Shiraishi 2019.12.23 Kuno-Yamanaka Group Year-End Presentation 2019

Parabola Interpolation Method

17

ptime

Moving Average Calculation of ptime

Peak Search

Nominal Time Condition

slide-18
SLIDE 18

Ryota Shiraishi 2019.12.23 Kuno-Yamanaka Group Year-End Presentation 2019

Rates in 10 - 20 clock timing

18

Ratio = Phys./TMON = 1 ~ 1.1

Energy Threshold FBAREne > 2MeV

10 10 20 30 40 50 60 70

Rate [Hz]

FBAR Rate ch0 (E>2.000000MeV)

100 200 300 400 500 600

3

10 ×

FBAR Rate ch0 (E>2.000000MeV)

FBARPTime [clock]

10 − 10 20 30 40 50 60 70

ratio (phys / tmon)

0.5 1 1.5

∆T1

FBAR Rate ch0 (E>2 MeV)

FBARPTime [clock] Rate [Hz] Ratio (phys./tmon)

− : Physics − : TMON

Rate = #events/(∆T×#triggered) ∆T = ∆T1

➞ Good agreement

5 10 15 20 25 30

Rate [Hz]

50 100 150 200 250 300

3

10 ×

(FBARPTime>10.000000 && FBARPTime<20.000000) FBAR Rate E>2.000000MeV

FBARModID

5 10 15 20 25 30

ratio (phys / tmon)

0.9 1 1.1 1.2

E > 2MeV, 10 < FBARPTime < 20 E > 2MeV, 10 < FBARPTime < 20 FBARModID Rate [Hz] Ratio (phys./tmon) Outer Channels Inner Channels 300 ×103

slide-19
SLIDE 19

Ryota Shiraishi 2019.12.23 Kuno-Yamanaka Group Year-End Presentation 2019

Rates in 35 - 50 clock timing

19

Energy Threshold FBAREne > 2MeV Rate = #events/(∆T×#triggered) ∆T = ∆T2

10 10 20 30 40 50 60 70

Rate [Hz]

FBAR Rate ch0 (E>2.000000MeV)

100 200 300 400 500 600

3

10 ×

FBAR Rate ch0 (E>2.000000MeV)

FBARPTime [clock]

10 − 10 20 30 40 50 60 70

ratio (phys / tmon)

0.5 1 1.5

∆T2

FBAR Rate ch0 (E>2 MeV)

FBARPTime [clock] Rate [Hz] Ratio (phys./tmon)

Ratio = Phys./TMON = 1.1 ~ 1.2

− : Physics − : TMON

5 10 15 20 25 30

Rate [Hz]

50 100 150 200 250 300

3

10 ×

(FBARPTime>35.000000 && FBARPTime<50.000000) FBAR Rate E>2.000000MeV (FBARPTime>35.000000 && FBARPTime<50.000000) FBAR Rate E>2.000000MeV

FBARModID

5 10 15 20 25 30

ratio (phys / tmon)

0.9 1 1.1 1.2

FBARModID Rate [Hz] Ratio (phys./tmon) E > 2MeV, 35 < FBARPTime < 50 Outer Channels Inner Channels 300 ×103

➞ A little bit worse than earlier timing