Calorimeter respons Helga Holmestad 11. April 2013 Helga Holmestad - - PowerPoint PPT Presentation

▶

Jun 22, 2023 272 likes •346 views

Calorimeter respons Helga Holmestad 11. April 2013 Helga Holmestad DHCal 11. April 2013 1 / 22 Particle selection Distinguish between muons/electrons/pions How: Cherenkov counters (electrons) Density (Number of hits/ Active

SLIDE 1

Calorimeter respons

Helga Holmestad

11. April 2013

Helga Holmestad DHCal

11. April 2013

1 / 22

SLIDE 2

Particle selection

Distinguish between muons/electrons/pions
How:
Cherenkov counters (electrons)
Density (Number of hits/ Active Layers)
Total RMS/Active Layers
Make histogram of these variables. Then two clearly separable

regions should be seen. These plots decide where to put the cuts.

Helga Holmestad DHCal

11. April 2013

2 / 22

SLIDE 3

7 GeV beam

TotalNumberOfHits/activeLayers 2 4 6 8 10 12 2000 4000 6000 8000 10000 12000 14000 16000 18000 20000 density density Entries 165526 Mean 3.474 RMS 1.754 density TotalRMS/activeLayers 10 20 30 40 50 60 70 1000 2000 3000 4000 5000 6000 7000 8000 9000

TotalRMS/ActiveLayers

Entries 165526 Mean 11.97 RMS 10.32

TotalRMS/ActiveLayers

270 GeV beam

TotalNumberOfHits/activeLayers 10 20 30 40 50 60 70 80 90 500 1000 1500 2000 2500 3000 3500 4000 density density Entries 30099 Mean 31.93 RMS 20.73 density TotalRMS/activeLayers 10 20 30 40 50 60 70 80 90 100 200 300 400 500 600 700 800

TotalRMS/ActiveLayers

Entries 30099 Mean 44.34 RMS 26.57

TotalRMS/ActiveLayers

Helga Holmestad DHCal

11. April 2013

3 / 22

SLIDE 4

The cuts

The cuts used are showed in the table below.
Because the tailcatcher hits are not used I tried to take away

events which could leak energy. (Only done for runs with higher energy than 50 GeV).

This was done by excluding the events which had less than 25

hits in the 10 first layers, this should exclude the events with late shower start. Energy density cut centerRMS cut 7 GeV 2 1.8 10 GeV 2 1.8 50 GeV 3 3 270 GeV 12 23 300 GeV 15 25

Helga Holmestad DHCal

11. April 2013

4 / 22

SLIDE 5

Fitting a Gaussian

The number of hits from the pion events is assumed to be

Gaussian distributed, when no saturation occurs.

I have fitted a Gaussian function in the interval ± RMS from

the arithmetic mean.

This is done to be independent of outliers
Outliers whit large energy looks like box events not found by

the Box Event finder.

Outliers with low energy?

7 GeV

Number of hits 20 40 60 80 100 120 140 160 1000 2000 3000 4000 5000 6000 7000 8000

Hit distribution

Hit distribution Entries 165526 Mean 67.1 RMS 19.75 Before cuts After cuts

Hit distribution

270 GeV

Number of hits 500 1000 1500 2000 2500 3000 500 1000 1500 2000 2500 3000 3500

Hit distribution

Hit distribution Entries 30099 Mean 1162 RMS 762.3 Before cuts After cuts

Hit distribution

Helga Holmestad DHCal

11. April 2013

5 / 22

SLIDE 6

For each energy I take the mean of the fitted Gaussian.
Need to do more...

Helga Holmestad DHCal

11. April 2013

6 / 22

SLIDE 7

A first look at the MC data

10 GeV

Number of hits 20 40 60 80 100 120 140 160 180 50 100 150 200 250 300 350

Hit distribution

Hit distribution Entries 9966 Mean 100.8 RMS 24.06 Before cuts After cuts

Hit distribution

50 GeV

Number of hits 100 200 300 400 500 600 700 800 20 40 60 80 100 120

Hit distribution

Hit distribution Entries 2948 Mean 395.5 RMS 89.23 Before cuts After cuts

Hit distribution

Energy efficiency muons going trough the cuts 5 GeV 97.81 % 50 GeV 97.76 % 300 GeV 0.275%

Helga Holmestad DHCal

11. April 2013

7 / 22