Kalorimetry calorimeters

J. Žáček Experimentální metody jaderné a subjaderné fyziky

1. Proces energetických ztrát je statistický DE/E ~ 1/√ E process of energy losses is statistical ies2. Rozměry kalorimetrů ~ ln (E0) dimensions of calorimeters 3. Není třeba magnetického pole no magnetic field is needed4. Lze je segmentovat, tj. lze měřit i směr pohybu částice they can be segmented i.e. the direction of particle motion can be measured 5. Lze je použít na spouštění trigrů they are used in triggers6. Radiační poškození Radiation demnage

Measurement of energies of particles at higher energies, when a cascade process (i.e. a shower) is initiated

Kalorimetr obsahuje: -pasivní materiál (většinou) -aktivní materiál (vždy) v němž se detekuje deponovaná energie tzv. „viditelná energie“

sprška

A. Homogeneous calorimeters

B. Sampling calorimeters B.1 electromagnetic - electromagnetic showers B.2 hadronic - hadronic showers B.3. calorimeter response B.4. calibration B.5. improvement of the resolution of hadron calorimeters

C. Calorimeter with liquid Ar D. SPACAL calorimeter with scintillation fibres

used for detection of particles which interact electromagnetically i. e. electrons, positrons, photons, (muons ?)

A. Homogeneous calorimeters

Homogenní kalorimetr-olovnaté sklo

Homogeneous calorimeter - lead glass

OPAL experiment, CERN, collider LEP electrons (~50 GeV) vs positrons (~50 GeV)

B. Sampling kalorimetry sampling calorimeters

Elektromagnetické : detekce fotonů, elektronů Jejich princip je založen na šíření elektromagnetických spršek

Electromagnetic : detection of photons and electrons method : electromagnetic showers

Hadronové: detekce hadronů (piony, protony, neutrony ..) hadronové spršky

Hadron calorimeters : hadron detection (pions, protons, neutrons..) hadronic showers

Absorbers – shower developmentActive (detection) layer – no shower development, ionization only

B.1 elektromagnetické “sampling” kalorimetry - elektromagnetické electromagnetic calorimeters electromagnetic showers spršky

Jednoduchý model : a simplified model of shower developement in a medium

1) Každý e+, e- s E > Ec urazí 1 radiační délku a vyzáří 1/2 své energie

each e+, e- > transverses 1 radiation length end emits ½ of its energy

2) Každý foton s E > 2 urazí 1 radiační délku s a vytvoří symetrický pár e+, e-

each photon with E > 2 transverses 1 radiation length and creates symmetric pair

3) Pro E<Ec žádné brzdné záření, pouze ztráty ionizací for E<Ec no bremstrahlung, ionization losses only

4) Pro E > Ec zanedbáme ionizační ztráty for E > Ec ionization losses are neglected

1) Počet částic v hloubce t N(t)=2 t number of particles at depth t

2) E(t) / částici = E0 / 2 t , E/particle

3) Hloubka v níž částice mají energii E´ depth where partciles have energy E´ t(E´) = ln(E0 / E´) / ln2

4) Maximální hloubka pro E(t)=Ec

maximum depth for tmax = ln(E0 / Ec) / ln2 5) Celková délka drah nabitých částic total length of charged particles E0

T ≈ --- X0

Emission angles are neglectedAll tracks are parallel

T = 2X0 / 3 * Σ 2j + S0 * 2 / 3 * Nmax

S0 je průměrná dráha nabitých částic s energií menší než je kritická energie is average track length of charged particles with E<Ec

T = 2X0 / 3 * ( 2 tmax+1 - 1) + S0 * 2 / 3 * Nmax ≈ 4/3 X0 * E0 / Ec + S0 * 2 / 3 * E0 / Ec

T ~ E0 / Ec * X0

Pokud detekujeme elektrony od nějaké energie Ed jeIf electrons are detected from the energy Ed

T = X0 * F(k) * E0/Ec

k je parametr, definovaný např. jakok is a parameter defined e.g. as

k=2.29 * Ed / Ec

F(k) je experimentálně nalezená formule is determined experimentally

F(k) = ek ( 1 + k * ln ( k / 1.526))

Pro k → 0 je F(k) → 1For is

Podélný profil spršek longitudinal shower profile

Parametrizace deponované energie v závislosti na hloubce t (X0)získaná ze simulací a z měření

dE ba = E0 ta e -b t

dt G(a)

a, b parametry, G - gama funkce parameters, gamma function

tmax = (a -1)/b

Deposited dE/dt energy at the depth t

DeponovanáEnergie(lib. jednotky)

Deposited energy(arbtr. units)

deposited energy/ total energy

Příčný profil elektromagnetických spršek

Energie v příčném směruEnergy in the transverse directionE(r) ≈ C * exp(-4r/RM)C konstanta, a constant

Přesněji dvě komponenty In fact 2 componentsE(r) = A * exp(-br) + C *exp(-dr), A,b,C,d konstanty constants

Transverse (lateral) profile

photon) = 9/7 )

n vrstev detekčního prostředí „s“, n detection (active) layers “s”

Deponovaná energie, ( dE/dx jsou ionizační ztráty v detekčním prostředíDeposited energy , ( dE/dx ion. loss of in the detection material

Es = ( dE/dx * ds ( N1 + N2 + N3 + …… Nn )

Ni počet nabitých částic ve vrstvě „i“ number of charged particles in the layer “ï”Es = (dE/dx * ds * N

N je celkový počet nabitých částic prošlých detekčním prostředím total number of charged particles traversing detection layers

absorbátor

Active layer:Detekční médium svelkým X0

Detection material with high X0

absorber

Assumption (dE/dx constant

Rozlišení resolution

Rozvoj spršky probíhá tak dlouho, až energie elektronů a pozitronů v absorbátoru poklesne na kritickou energii, tj pro Fe je cca 23 MeV, pro Pb je cca 8 MeV. (kritická energie Ec ≈ 660/ (Z + 1.24) MeV )

V detekčním prostředí elektrony a pozitrony pouze ionizují. Vzhledem k jejich energii jsou jejich ionizační ztráty na minimum a málo závisí na energii. Proto v Es je dE/dx konstantní.

Showers end if positron and electron energies < critical energy

In the active layer electron and positron ionize only. Their ionization lossesare at minimum because their energies are sufficiently high.

( dE/dx ≈ constant

Ionization losses electrons – no brehmsstrahlung

Absorbátor: Celkový počet částic vytvořených vAbsorber: absorbátoru a které projdou do detekčního prostředítotal number of particles produced in the absorber and which enter detection mediumN = T/da = E0 / Ec * X0 / da

Energetické rozlišení v deponované energii je určeno fluktuacemi v N, energy resolution of deposited energy correspond to the fluctuation of N tj. DN=√N

Resolution: DEs / Es = DN / N = 1 / √N = 1/√E0 * √da * Ec / X0

Výše uvedený vztah platí za předpokladu, že všechny dráhy částice ve spršce jsou paralelní se směrem dopadající částice a že detekujemečástice všech energiích. Jestliže úhel emise sekundárních částice je q a detekujeme částice od nějaké prahové energie dostaneme

DEs / Es = 1 / √ E0 * √ da * Ec / [X0 * F(k) * < cos q >]

Modification of the resolution if energies are detected from a threshould energy and the emission angle of secondary particlesare taken into account

Celkové energetické rozlišení total energy resolution

ΔE E ( )2

√ E0

( ) ( )2 2+ A2 ln(E0)( )+

2 A3 * s0.5

( )+ + A4

Fluktuace ve spršce

Shower fluctuation

Fluktuace vdetekčnímSystému fluctuationof the detectionsystem

Ztráty na podélný rozměrLosses due to longitudinal leakage

Celkový elektronický šum, S- je šum vjednotkách energieTotal electronic noise, S- noise in the unit of energy

Nehomogenita non-homogenuity

homogenní kalorimetry

sampling kalorimetry

B.2. Hadron sampling calorimeters - hadronic showers hadronové „sampling“ kalorimetry – hadronové spršky

Electromagnetic shower e

Hadronic shower

Příspěvek k dE/dx pro protony s energií 5 GeV v železecontribution to dE/dx in iron for protons at 5 GeV

„viditelná“ energie visible energy „neviditelná“ energie invisible energy

Nabité piony, protony 40% vazbová energie 18 % charged pions, protons binding energy

Neutrální piony π0 - 2γ 17% neutrony +ostatní 17 % neutral pions neutrons + others

Jaderné fragmenty 8 %nuclear frgments

Podélný a příčný profil hadronových spršek longitudinal and lateral profile of hadron showers

charakteristická veličina je interakční délka λI = 0.35 A1/3 g cm-2

characteristic quantity interaction lenght nebo absorbční délka nezahrnuje elastický rozpty absorption length Počet částic ve spršce <n> = A0.1 ln(E2

tot )Number of particles in a showerVětšinou piony cca 90% . Mostly pions ~90%

Podélný profil longitudinal profile Energetické ztráty ve vzdálenosti l od počátku spršky na délce dl energy loss at the distance l from the start of shower in the lengh dl

dE(l) = E0 { (1 – c0 ) H(x) dx + c0 F(y) dy }

x= ahad * l / λI , y= ael * l / X0 , c0 frakce neutrálních pionů parameters fraction of neutral pions

(platí pokud se neutr. piony produkují pouze v primárním vrcholu)(valid if neutral pions are produced in the primary vertex)

Transverse profile r ⊥ direction of motionconstants

Hadron energy EValid for

Resolution of hadron calorimeters

Sampling term constant term

Vliv různých procesů na rozlišení hadronových kalorimetrůInfluence of various processes on the resolution of hadron cal.

Influence of various processes on the resolution of hadron cal.

Process characteristic feature influence on resolution

Hadron multiplicity of secondary Interactions particles fluctuation

Nuclear evaporation energy ~10 % loss of binding de-excitation binding energy 10 % energy neutrons 40 % bad detection protons 40 % of slow protons and slow neutrons decays energy loss of neutrinos and muons

Elektromagnetická část40 vrstev (Pb + scint.)Electromagnetic section

Hadronová část55 vrstev (Fe + scint.)Hadronic section

Sampling kalorimetr

Profil spršky pro piony o energii 270 GeV v kalorimetru o 90 vrstvácholova o tlouštce ¾ inche

Layer numberČíslo vrstvy

, součet přes mnoho pionů

Profil vzhledem k počátku spršky

Longitudinalshower profile forpions at 270 GeVin a calorimeter of 90 layers of Pb, thickness ¾ inch, sum over many pions

profile wrt. to the shower start

energie/vrstvu

Profil jednoho pionu o energii 270 GeV

energie/vrstvu

číslo vrstvy

One pion

Profil dalšíchdvou pionů

Other two pions

elektromagnetická sprška hadronová sprškaElectromagnetic shower hadronic shower longitudinal profile

Délka v podélném směru , kde je maximální deponovaná energie

X0 tmax ≈ [ 0.2 lnE (GeV)+ 0.7] ∙ λI

Délka v podélném směru , kde je deponováno 95 % energie

t95% = tmax + ( 0.08Z + 9.6) X0 t95% ≈ a lnE + b

Příčný profil lateral profile

95% energie je ve válci o poloměru energy in a cylinder of radius r=2 RM = 21 MeV/ Ec ∙ X0

Způsoben mnohonásobným rozptylem due to multiple scattering

95% energie je ve válci o poloměru r ~ λI

Způsoben fixní hodnotou příčné hybnostisekundárních částic ≈ 0.4 GeV Fixed transverse momentum of secondaryparticles ~ 0.4 GeV

Železo, X0 =1.76 cm, λI = 16.8 cm, Ec = 22.5 MeV, částice s E=100 GeV: a=9.4, b=39,

tmax = 21 cm, t95% ≈ 42 cmr = 3.2 cm

tmax ≈ 27 cm, t95% ≈ 80 cm r ~ 16.8 cm

Hadronové kalorimetry jsouHadron calorimeters are

- nelineární non-linear

- mají špatné rozlišení , kolem 0.9 / √E , (E v GeV)- they have bad resolution. about- - mají jinou odezvu pro hadrony než elektrony při stejné primární energii- they have different response for hadrons and electrons at the same primary energy

B.3 Odezva kalorimetrů Calorimeter response

Calorimeter response homogeneous em. calorimeters

t thickness of detection (active) material, of absorber (passive)

mips: minimum ioniz. particle

Calorimeter response sampling em. calorimeters

Miony jsou většinou na minimu ionizace

𝒇 𝑹(𝒆)/ 𝒇 𝑹(µ)

Calorimeter response hadron sampling calorimeters

Usually non compensating calorimeters,Hadronové kalorimetry jsou obvykle nekompenzační i.e. the response is different for electrons an hadrons tj. Odezva je různá pro elektrony a hadrony

Calorimeter signal : T pion signal:

Particle energy 21 GeVCalorimeter: Fe + scint.

Normalized to the pion signal

e/h ratio: Kompenzační kalorimetr

e electromag. energy, h hadron energy f fractions

mip: energy of a minimum ionizing particle usually muon, (rel, p, n corresponding energies

rel energie deponovaná relativistickými částicemi p energie dep. pomalými protony n neutrony inv neviditelná energie

f frakce

(tj. většinou pomalé protony)

( velmi silný signál, obvykle dosahuje maximální hodnoty, tj nezávisí na energie protonů – saturace)

Měření energií při velké multiplicitě sekundárních částic

Je třeba odlišit částice vstupující do kalorimetru

Vysoké primární energie hodně částic, měří se energie shluků částic tzv. jetů

Nutná příčná segmentace detekční části signál z tzv. buněk (cel)

buňka

metoda: klastry

Klustr, sečte se energie ze všech buněk stanoví se „centrum klustru“, tj. jejich x,y

Izolovaná buňka – ignoruje se

- Spojí se odpovídající klustry v jednotlivých vrstvách, sečte se energie- tj. dostaneme energii částice či jetu.

tj. klustry které mají centrum s téměř stejnými hodnotami souřadnic

Determination of energy from measured signals i.e. E=A + B * S S měřený signál např. napěťový puls, nebo náboj, measured chargeS- obvykle digitalizovaný analogový signál, (s pomocí amplitud-digital-convertor, ADC), což je číslo v jednotkách ADC, tj. udává např. celou plochu analogového signálu A, B, kalibrační parametry, které je třeba určit calibration parameters, have to be determined from experiments E je známá energie částice, použije se několik různých energií known particle energy, several energies are used kalibrace radioaktivními zdroji v jaderné fyzice calibration with radioactive sources in nuclear physics kalibrace svazky elektronů a hadronů, with electron and hadron beams kalibrace miony with muons

a) Kalibrace svazky elektronů je jednodušší , z ní se stanový tzv.elektromagnetická škála tj. A a B. The calibration with electrons electromagnetic scale., i.e. parameters A and B. S těmito parametry se určí odezva na hadrony a tím se určí sampling poměrA response to hadrons is determined with these parameters e/h ratiob) Testování odezvy detekčního prostředí v různých jeho částech, např. v rozích scintilátorů, uprostřed atd. Testing of the response in various part of active cells e.g. edges of scintillators etc.

B.4 Calibration

Kalibrace s miony

Muon calibrationwith muon beams at accelerators, muon energy is known

Energetickéztráty mionů

Muon energy cca 1 GeV to cca 100 GeV ≈ nearly at ionization minimum These muons are suitable for the calibration and testing of homogeneity ofcalorimeter cells. The energy losses are described by Landau formula.

Měření těchto ztrát v kalorimetru v jedné buňce (absorbátor Fe- scintilátor)Energy loss measurement in a cell of a calorimeter (Fe + scint.)The measured signal A in ADC (amplitude to digital convertor) units.For each muon which penetrates the cell ⟹ one ADC value

Pozadíbackground

nejpravd. hodnota Amip

The most probable value

1. Cely mají různé hodnoty Amip. Tyto hodnoty se překalibrují na nějakou střední hodnotu, tj každá hodnota se násobí konstantou C tak, aby Amip ∙ C = < Amip > cells have various values of Amip. These values are recalibrated to some mean value by a constant C, Amip ∙ C = < Amip >

uniform cell response

2. Energetická kalibrace: Miony se obvykle plně neabsorbují. Proto se energetické ztráty přesně spočítají, tj v jedné cele je ΔE = Cen ∙ < Amip >, kde Cen je energetická kalibrační konstanta energy calibration: muons are not absorbed ⟹ their energy losses are exactly calculated, for a cell ΔE = Cen ∙ < Amip >, where Cen is the calibration constant

Odezva kalorimetru na elektrony, piony a miony o energii8 GeV . Zkalibrováno na elektromagnetickou škálu

Response of a calorimeter to electrons, pions and muons at primary energy 8 GeVCalibrated on em. scale.

B.5 Zlepšení rozlišení hadronových kalorimetrů Improvement of the resolution of hadron calorimeters

- Hardware

- Software, tzv. metoda vážení, použitelná pro segmentované kalorimetry

Metoda: kompenzace odezvy kalorimetru na jednotlivé složky signálu (em, protony, neutrony..)

Kompenzaci lze získat: redukcí elektromagnetické odezvy zvýšením hadronové odezvy

e/h ~ 1

závisí: na energii, materiálu, vlastnostech „samplingu“

Způsob realizace kompenzace:

Hardware compensation reduction of the e-response and increase of the h-response

Resolutions hadronic ΔE/E ≈ 0.3/

Účinné průřezy interakce neutronů s uranem a vodíkem

i.e. weighting method, suitable for segmented calorimeters

E = Σ ci Ei, ci is the weighting constants in the cell I with the energy sum over all cells for c , ci = 1, electromagnetic scale other constants are determined from the simulation to get the best resolution

could be achieved hadronic resolution ΔE/E ≈ 0.5/

Software compensation

iron sampling calorimeter(absorber Fe)

software compensation

Příklady sampling kalorimetrů

Spršky ve vzduchu, tlak 1 atmShowers in air, pressure 1 atm

X0 = 304 m λI = 745 m, tmax = 1200 m r=745 m

C. Kalorimetr s kapalným argonem calorimeter with liquid Ar

Kapalný Ar : liquid Ar hustota 1g/cm3 density nezachycuje elektrony electronegative pohyblivost elektronů 5 ∙ 105 cm/s při napětí 1 kV/mm electron mobility at the voltage ionizační potenciál 26.5 eV ionization potential (dE/dx)min = 2.11 MeV/cm nízká teplota 86o K low temperature

Jaký je indukovaný náboj od ionizačních elektronů ?What is the induced charge of ionization electrons?

Absorbátor absorber

Elektrodaelectrod

Δq = Q Δx/d

Induced charge

Induced current in the external circuit :

Calorimeter cell

Primární částice primary particle

ionizační elektronyionization electrons

V čase t=0 je celkový náboj ionizačních elektronů Q0

At t = 0 total electron charge Q0

Počet elektronů se mění,neboť se pohybují ke kladné elektrodě, kde jsou neutralizovány v je rychlost elektronů v is electron velocity

počet

Induced current at time t:

is initial number of electrons, is number of electrons at time t

A decrease of electrons dn during dt passing an unit area

is electron density between electrodes

= 3.8 ∙10-15 C

Induced current at time t,

Total induced charge

example

total collection time of electrons d/v

Příklad strukturyhadronového kalorimetru z kapalného Ar

the structure ofa hadronic calorimeter

D. SPACAL kalorimetr ze scintilačních vláken scintillating fibre calorimeter

Olovo a scintilační vlákna, Pb + scintillating fibres, (vlákna průměr 0.5 mm, délka 30 cm, fibres: diameter 0.5 mm, length 30 cm Pb listy o tlouštce 0.8 mm, 40x40x200 mm) Pb sheets: thickness 0.8.mm, 50x40x200 mm)

• dobře měřený příčný profil lateral profile well measured• neměřený podélný profil longitudinal profile not measured• dobrá identifikace elektronů (p/e ~10-4 ) good electron identification • nekompenzační, poměr signálu e/h=1.3 , non-compenzating , e/h=1.3• dobré rozlišení pro elektrony ~ 7 % good resolution for electrons ~7%

2 fotonásobiče

Příčný řez kalorimetrem SPACALTransverse cut of the calorimeter SPACAL

Průměr 1.5 m, diameter 1.5. m

Urychlovací trubiceAccelerator pipe

Submodul 8x4 cm

Spacal montáž do aparatury experimentu H1Spacal installation into the detection system of the experiment H1

Kalorimetry calorimeters

Documents

Transcript of Kalorimetry calorimeters

Exploitation of Jet Properties for Energy Scale Corrections for the CMS Calorimeters · 2011. 3. 2. · Exploitation of Jet Properties for Energy Scale Corrections for the CMS Calorimeters

Precision Crystal Calorimeterszhu/talks/ryz_080402_crystals.pdf · 2008. 4. 5. · Fermilab Colloquium, April 2, 2008 Precision Crystal Calorimeters in High Energy Physics: Past,

Crystal Calorimeters at Linear Colliderhep.caltech.edu › ~zhu › talks › ryz_lc_0201.pdf · Crystal Calorimeters at Linear Collider Ren-yuan Zhu California Institute of Technology-DQ

Year end presentation 20171228 - Osaka Universityosksn2.hep.sci.osaka-u.ac.jp/2017Nenmatsu/slides/...2017/12/28 · # of CsI Calorimeters: 2240(Small) + 476(Large) 5 cm 5 cm 2.5 cm

Performance of the LHCb calorimeters · Y~7m X~8.5m Z~2.7m PS/SPD LHCb Calorimeters 20131125 V CPAN Days - Xvc 2 10#250&mrad 10#300&mrad VELO Vertexing MAGNET Calorimeters PID: e,γ,

Metallic Magnetic Calorimeters for spectrometry applications

CALIBRATION BOARDS FOR THE LAr CALORIMETERS ATLAS N. Dumont-Dayot, M. Moynot, P. Perrodo, G. Perrot, I. Wingerter-Seez Laboratoire d’Annecy-Le-Vieux de.

Spring, 2009Phys 521A1 Hadronic calorimeters Recall that λ I > X 0 for dense materials hadronic calorimeters must be longer than EM calorimeters This,

Calorimeters for Precision Timing Measurements in High ...bornheim/html/CALOR2014_Bornheim.pdf · Calorimeters for Precision Timing ... High energy physics collider experiments ...

Precision Crystal Calorimeters in High Energy Physics Past ... · April 4, 2006 International Symposium on Detector Development, SLAC, USA 1 Precision Crystal Calorimeters in High

Crystal Calorimeters in the Next Decadezhu/talks/ryz_080526_calor.pdfInternational Conference on Calorimetry in Particle Physics, Pavia, Italy Crystal Calorimeters in the Next Decade

BNDMarcelRaas - TU Dresden · 2009. 9. 15. · 0.01 0.015 0.02 0.025 0.03 0.035 0.04 0.045 0.05 Muon Identification in the ATLAS calorimeters 9/15/09 11 . Outlook & conclusions Muon

IDEA Meeting Cuoricino status, current issues...IDEA Meeting Cuoricino status, current issues Luca Gironi Pi Paris – 19 - 20 N b 200720 November 2007. Cuoricino TeO 2 thermal calorimeters

1 Calorimeters of the Very Forward Region Iftach Sadeh Tel Aviv University DESY Collaboration High precision design March 5 th 2008.

Light collection in “shashlik” calorimeters

1 국제 Working Group 활동 보고 (Detector) Brief Report on LCWS05 at SLAC ILC Detector concept Variety of Trackers, Calorimeters, Muons Beam tests and Organization.

Design and R&D of very forward calorimeters for detectors at future e

Role of the Andreev current in bolometers and calorimeters …ltd16.grenoble.cnrs.fr/IMG/UserFiles/Images/ORAUX/12_LTD16_faivre.pdf · Role of the Andreev current in bolometers and

Performance of the LHCb calorimeters during the period 2010-2012

Calorimeters for the linear collider detector