Stopping Target Optimization for COMET Phase II Giorgi - - PowerPoint PPT Presentation

▶
stopping target optimization for comet phase ii
SMART_READER_LITE
LIVE PREVIEW

Stopping Target Optimization for COMET Phase II Giorgi - - PowerPoint PPT Presentation

Jun 28, 2023 175 likes •451 views

Stopping Target Optimization for COMET Phase II Giorgi Kumsiashvili M2 Year-end Presentation 2019 1 COMET Phase II Beamline Electron Spectrometer Production Muon Target Stopping Target Detector Region Torus 1 Torus 2

slide-1
SLIDE 1

Stopping Target Optimization 
 for COMET Phase II

Giorgi Kumsiashvili M2

1

Year-end Presentation 2019

slide-2
SLIDE 2

COMET Phase II Beamline

Electron Spectrometer Detector 
 Region Muon
 Stopping
 Target Torus 2 Torus 1 Production 
 Target

2

slide-3
SLIDE 3

Original Target Design

htemp Entries 73341 Mean 101.8 Std Dev 11.76 20 40 60 80 100 sqrt(pxa[0]**2+pya[0]**2+pza[0]**2) 5000 10000 15000 20000 25000 30000 35000 40000 45000 htemp Entries 73341 Mean 101.8 Std Dev 11.76

sqrt(pxa[0]**2+pya[0]**2+pza[0]**2) {hit0}

Target Radius 100mm Beam Blocker - Vacuum DIO Blocker - Vacuum Default Gradient Mag. Field

3

slide-4
SLIDE 4

Signal Electrons

Signal Electrons that are
 later accepted by the detector

4

slide-5
SLIDE 5

Distribution in Straw Tracker

5

slide-6
SLIDE 6

Test Designs

1st Disk - 50mm, 9th - 100mm, 17th - 50mm 1st Disk - 150mm,17th - 50mm 1st Disk - 50mm, 17th - 150mm 1st Disk - 100mm, 9th - 50mm, 17th - 100mm

6

slide-7
SLIDE 7

Distributions In Straw Tracker

800 − 600 − 400 − 200 − 200 400 600 800 ya 10200 10400 10600 10800 11000 11200 11400 xa 200 400 600 800 1000 1200 1400 1600 1800 2000 2200

xa:ya {hit0&&hit1&&hit2&&hit3&&hit4&&xa>10000}

50mm-100mm-50mm

800 − 600 − 400 − 200 − 200 400 600 800 ya 10200 10400 10600 10800 11000 11200 11400 11600 xa 200 400 600 800 1000 1200 1400 1600 1800 2000

xa:ya {hit0&&hit1&&hit2&&hit3&&hit4&&xa>10000}

50mm-150mm

800 − 600 − 400 − 200 − 200 400 600 800 ya 10200 10400 10600 10800 11000 11200 11400 xa 200 400 600 800 1000 1200 1400 1600 1800 2000

xa:ya {hit0&&hit1&&hit2&&hit3&&hit4&&xa>10000}

100mm-50mm-100mm

800 − 600 − 400 − 200 − 200 400 600 800 ya 10200 10400 10600 10800 11000 11200 11400 xa 200 400 600 800 1000 1200 1400 1600 1800

xa:ya {hit0&&hit1&&hit2&&hit3&&hit4&&xa>10000}

150mm-50mm

7

slide-8
SLIDE 8

Momentum

htemp Entries 74078 Mean 102.3 Std Dev 10.88 20 40 60 80 100 sqrt(pxa[0]**2+pya[0]**2+pza[0]**2) 5000 10000 15000 20000 25000 30000 35000 40000 45000 htemp Entries 74078 Mean 102.3 Std Dev 10.88

sqrt(pxa[0]**2+pya[0]**2+pza[0]**2) {hit0}

htemp Entries 72286 Mean 102 Std Dev 11.22 20 40 60 80 100 sqrt(pxa[0]**2+pya[0]**2+pza[0]**2) 5000 10000 15000 20000 25000 30000 35000 40000 45000 htemp Entries 72286 Mean 102 Std Dev 11.22

sqrt(pxa[0]**2+pya[0]**2+pza[0]**2) {hit0}

50mm-100mm-50mm 50mm-150mm

htemp Entries 74059 Mean 102.2 Std Dev 11.22 20 40 60 80 100 sqrt(pxa[0]**2+pya[0]**2+pza[0]**2) 5000 10000 15000 20000 25000 30000 35000 40000 45000 htemp Entries 74059 Mean 102.2 Std Dev 11.22

sqrt(pxa[0]**2+pya[0]**2+pza[0]**2) {hit0}

100mm-50mm-100mm

htemp Entries 72908 Mean 101.9 Std Dev 11.91 20 40 60 80 100 sqrt(pxa[0]**2+pya[0]**2+pza[0]**2) 5000 10000 15000 20000 25000 30000 35000 40000 htemp Entries 72908 Mean 101.9 Std Dev 11.91

sqrt(pxa[0]**2+pya[0]**2+pza[0]**2) {hit0}

150mm-50mm 74.07% 72.28% 74.05% 72.9%

8

slide-9
SLIDE 9

Initial Distribution

100’000 Electrons

9

slide-10
SLIDE 10

Mu- Momentum in Disks

hMomentumAll_StoppingTargetPhaseIITargetDisk_0_13

Entries 6222 Mean 51.39 Std Dev 14.04

20 40 60 80 100 120 140 160 180 Momentum (MeV/c) 20 40 60 80 100 120 140

hMomentumAll_StoppingTargetPhaseIITargetDisk_0_13

Entries 6222 Mean 51.39 Std Dev 14.04

Momentum at StoppingTargetPhaseIITargetDisk_0 for mu-

In this example: Number of disks - 17 Disk radius - 100mm Disk 1 Disk 17

hMomentumAll_StoppingTargetPhaseIITargetDisk_16_13

Entries 4012 Mean 48.72 Std Dev 15.27 20 40 60 80 100 120 140 160 180 200 Momentum (MeV/c) 20 40 60 80 100 120 140

hMomentumAll_StoppingTargetPhaseIITargetDisk_16_13

Entries 4012 Mean 48.72 Std Dev 15.27 Momentum at StoppingTargetPhaseIITargetDisk_16 for mu-

Around 64.5%

slide-11
SLIDE 11

Mu- Momentum in Disks

hMomentumAll_StoppingTargetPhaseIITargetDisk_0_13

Entries 6374 Mean 51.32 Std Dev 14.18 20 40 60 80 100 120 140 160 180 Momentum (MeV/c) 20 40 60 80 100 120 140

hMomentumAll_StoppingTargetPhaseIITargetDisk_0_13

Entries 6374 Mean 51.32 Std Dev 14.18

Momentum at StoppingTargetPhaseIITargetDisk_0 for mu-

In this example: Number of disks - 34 Disk radius - 100mm Disk 1 Disk 34

hMomentumAll_StoppingTargetPhaseIITargetDisk_33_13

Entries 2374 Mean 49.22 Std Dev 15.59 20 40 60 80 100 120 140 160 180 200 Momentum (MeV/c) 20 40 60 80 100 120 140

hMomentumAll_StoppingTargetPhaseIITargetDisk_33_13

Entries 2374 Mean 49.22 Std Dev 15.59 Momentum at StoppingTargetPhaseIITargetDisk_33 for mu-

Around 37%

slide-12
SLIDE 12

Summary and future work

  • Run larger scale simulations for usage in signal acceptance tests.
  • See if using 34 target disks hurt signal acceptance or not. (Smaller disks?.)

12

  • The original muon stopping target design is not optimal.
  • Test target designs yielded marginal improvements. There is a need to obtain


realistic muon stopping distribution in target disks. (in progress)

  • Significant number of Muons go through the target disks and downstream.

What’s next.. Thank you for your attention!

slide-13
SLIDE 13

13

slide-14
SLIDE 14

BACKUP

14

slide-15
SLIDE 15

Initial Distribution (That hit all 5 stations)

72’908 Electrons

15

slide-16
SLIDE 16

COMET Phase II Background Estimations

From COMET Phase II note

18

slide-17
SLIDE 17

Decay-in-orbit

  • DIO Spectrum

From Czernicki’s Calculations Logarithmic Scale Linear Scale

19

slide-18
SLIDE 18

Analysis Workflow

Simulation Using ICEDUST Analysis of Simulation output 
 using Analysis Code Drawing Histograms and calculation 
 using CERN ROOT Raw Output .root file Analysed .root file (TTree)

20

slide-19
SLIDE 19

Signal and DIO Simulation

  • 500000 Muons uniformly distributed in 17 disks
  • Beam Blocker removed
  • DIO Blocker removed

htemp

Entries 495730 Mean 1435 RMS 244.8

lastposz 1000 1200 1400 1600 1800 5000 10000 15000 20000 25000 30000

htemp

Entries 495730 Mean 1435 RMS 244.8

The number of stopped muons in each disks

2164 2 2 27360 1 2 5 29111 2 1 1 1 29197 1 3 29152 2 2 3 29665 3 2 1 29249 1 2 3 1 29168 1 4 2 2 29161 4 3 1 29300 3 3 2 2 29096 3 2 1 3 29371 4 1 4 29357 2 2 2 2 29346 2 1 1 3 29427 2 2 1 29369 1 3 1 28832 3 3 2 27278 2 1 1 2 1 1

htemp

Entries 495730 Mean 1435 RMS 244.8

htemp Entries 500000 Mean 1437 RMS 245.6

iniposz 1000 1100 1200 1300 1400 1500 1600 1700 1800 1900 5000 10000 15000 20000 25000 30000

htemp Entries 500000 Mean 1437 RMS 245.6

Initial electron with 105MeV

996 28342 29573 29563 29616 29466 29242 29444 29207 29601 29594 29269 29068 29289 29426 29357 29455 29492

htemp Entries 500000 Mean 1437 RMS 245.6

21

slide-20
SLIDE 20

DIO Electrons

22

slide-21
SLIDE 21

DIO Distribution in Spectrometer

Plots by Wu Chen

23

slide-22
SLIDE 22

DIO Spectrum Analysis

dio2017 Entries 2200 Mean 99.88 Std Dev 0.6887 100 101 102 103 104 105 Momentum [MeV/c]

24 −

10

23 −

10

22 −

10

21 −

10

20 −

10

19 −

10

18 −

10

17 −

10

16 −

10

15 −

10

14 −

10

13 −

10 dio2017 Entries 2200 Mean 99.88 Std Dev 0.6887

DIO 2017

DIO Spectrum with logarithmic scale DIO Spectrum recorded in the Straw Tracker From Czarnecki’s Calculations

h1

Entries 260000 Mean 101.1 Std Dev 0.571

100.5 101 101.5 102 102.5 103 103.5 104 104.5 105

14 −

10

13 −

10

12 −

10

11 −

10

10 −

10 h1

Entries 260000 Mean 101.1 Std Dev 0.571

DIO Hist 25

slide-23
SLIDE 23

Smeared DIO

DIO Spectrum smeared with a Gaussian 
 with sigma of 150 KeV

h1

Entries 260000 Mean 101.1 Std Dev 0.571

100.5 101 101.5 102 102.5 103 103.5 104 104.5 105

14 −

10

13 −

10

12 −

10

11 −

10

10 −

10 h1

Entries 260000 Mean 101.1 Std Dev 0.571

DIO Hist

Big thanks to Chen Wu
 for helping with the code! MeV/c

26

slide-24
SLIDE 24

Acceptance and Improvement factors

28

slide-25
SLIDE 25

Sensitivity

29

slide-26
SLIDE 26

Stopping Rates

30

slide-27
SLIDE 27

DIO Momentum Distribution on each straw station

31