status of geant4 veto simulation

Status of Geant4 Status of Geant4 VetoVeto SimulationSimulation

E.Leonardi, M.Serra, P.ValenteE.Leonardi, M.Serra, P.ValenteINFN RomaINFN Roma

NA48/3 Veto group meeting – April 20NA48/3 Veto group meeting – April 20thth, 2005, 2005

20/04/200520/04/2005 E.Leonardi, M.Serra, P.Valente - NA48/3 Veto Group MeetingE.Leonardi, M.Serra, P.Valente - NA48/3 Veto Group Meeting 22

GoalsGoals

Geant4 simulation of Geant4 simulation of large angle large angle vetoveto– CKMCalCKMCal: Lead/Scintillating tiles with : Lead/Scintillating tiles with

WLS readoutWLS readout– SpaCalSpaCal: Lead/Scintillating fibers with : Lead/Scintillating fibers with

lightguide readoutlightguide readout Use layout in ProposalUse layout in Proposal


GeometryGeometry

CKMCalCKMCal– Two veto rings, 16 trapezoidsTwo veto rings, 16 trapezoids– Tile planes Tile planes to beam axis to beam axis– Tile planes staggered in Tile planes staggered in – Aluminum support enclosureAluminum support enclosure

SpaCalSpaCal– Fibers Fibers to beam axis to beam axis– Aluminum support enclosureAluminum support enclosure


Beam Axis

27.7

5cm

42.5 cm

30

cm

Aluminum ScintillatorLead

1 mm

1 cm

CKMCalCKMCal

Flat bottom face…


Beam Axis

Aluminum ScintillatorLead

1 cm

1 mm

SpaCalSpaCal

2.5 m

24 c

m

29

.9

cm

C-shape module straightened…


Geometrical Geometrical ParametersParameters Aluminum supportAluminum support

– Covers front, back, and bottom surfaces of brickCovers front, back, and bottom surfaces of brick– Thickness=0.1, 1, 0.1 cmThickness=0.1, 1, 0.1 cm

FibersFibers– KLOE-style hexagonal arrangementKLOE-style hexagonal arrangement =1 mm, =1 mm, ΔΔx=1.35 mm, x=1.35 mm, ΔΔy=1.17 mmy=1.17 mm– ~70000 fibers/brick~70000 fibers/brick 20 X20 X00

TilesTiles– 5 mm Scintillator, 1 mm Lead5 mm Scintillator, 1 mm Lead– 68 tiles/brick68 tiles/brick 16 X16 X00


CKMCal readoutCKMCal readout

4 PMT’s per sector4 PMT’s per sector 64 PMT’s per veto ring64 PMT’s per veto ring


SpaCal readoutSpaCal readout

4x4 or 4x6 PMT’s per 4x4 or 4x6 PMT’s per sideside

64 or 96 PMT’s per veto 64 or 96 PMT’s per veto ringring

Winston lightguideWinston lightguide


Positions, Positions, and E and E

small small 0.1 – 1 GeV0.1 – 1 GeV

small small 0.1 – 1 GeV0.1 – 1 GeV

large large 0.05 – 0.1 GeV0.05 – 0.1 GeV

large large 0.05 – 0.1 GeV0.05 – 0.1 GeV

11

22

33

4455


Digitization: CKMCalDigitization: CKMCal

LLyy==1010 p p.e..e./MeV/MeV– ConservativeConservative light light yield from CKM yield from CKM

measurementsmeasurements– Includes light production, attenuation, Includes light production, attenuation,

collection and conversioncollection and conversion

NNp.e.p.e.==i i ΔΔEEi i xx L Lyy


CKMCal photoelectronsCKMCal photoelectrons

thr.=5 p.e.

PMT 0 PMT 1

PMT 2 PMT 3

100 MeV 100 MeV , position 2, position 2

13131414

1515

1212


CKMCal photoelectronsCKMCal photoelectrons100 MeV 100 MeV , position 2, position 2


CKMCal photoelectronsCKMCal photoelectrons100 MeV 100 MeV , position 2, position 2

1213

PMT 0PMT 0

PMT 1PMT 1

PMT 2PMT 2PMT 3PMT 3


CKMCal photoelectronsCKMCal photoelectrons100 MeV 100 MeV , position 2, position 2Slice 12 Slice 13

PMT 0PMT 0 PMT 0PMT 0PMT 1PMT 1 PMT 1PMT 1

PMT 2PMT 2 PMT 2PMT 2PMT 3PMT 3 PMT 3PMT 3


CKMCal p.e.CKMCal p.e.

Highest photoelectron signal distributionInefficiency vs. threshold

No digi=49 events


Digitization: SpaCalDigitization: SpaCal

L=25 ph/MeVL=25 ph/MeV Light productionLight production

==320 cm320 cmLight attenuation lengthLight attenuation length

C=0.93C=0.93 Light-guide efficiencyLight-guide efficiency

Q=0.2Q=0.2 Photocathode quantum Photocathode quantum efficiencyefficiency

NNp.e.p.e.==i i ΔΔEEi i xx L L xx exp(-d exp(-dii//) ) xx C C xx Q Q

L and L and for Kuraray SCSF-81 from for Kuraray SCSF-81 from A.Antonelli et al. NIM A 370 (1996) 367A.Antonelli et al. NIM A 370 (1996) 367


SpaCal photoelectronsSpaCal photoelectrons

100 MeV 100 MeV , position 1, position 1 100 MeV 100 MeV , position 2, position 2


ResultsResults

veto 0 veto 1

veto 0 veto 1 SpaC

al

SpaC

al

CK

MC

al

CK

MC

al


PositionPosition EnergyEnergy No digiNo digi

in 0in 0No digiNo digi

in 0+1in 0+1No No interactioninteraction

0+10+1

InefficiencyInefficiency

11 100 MeV100 MeV 608608 608608 606606 6.08e-36.08e-3

1 GeV1 GeV 163163 163163 163163 1.63e-31.63e-3

22 100 MeV100 MeV 5353 4949 4444 4.9e-44.9e-4

1 GeV1 GeV 66 55 55 5e-55e-5

33 100 MeV100 MeV 13291329 767767 00 7.67e-37.67e-3

1 GeV1 GeV 1010 11 00 1e-51e-5

44 100 MeV100 MeV 57965796 11841184 33 1.18e-21.18e-2

1 GeV1 GeV 18911891 33 11 3e-53e-5

55 100 MeV100 MeV 9974999749 22882288 4343 2.29e-22.29e-2

1 GeV1 GeV 9918199181 1414 99 1.4e-41.4e-4

11

22

33

4 54 5

small small PositionPosition EnergyEnergy No digiNo digi



0+10+1


11 100 MeV100 MeV 2121 1919 1515 1.9e-41.9e-4

1 GeV1 GeV 11 11 00 1e-51e-5

22 100 MeV100 MeV 1414 1313 1313 1.3e-41.3e-4

1 GeV1 GeV 11 11 11 1e-51e-5

33 100 MeV100 MeV 29832983 28092809 00 2.81e-22.81e-2

1 GeV1 GeV 4949 55 00 5e-55e-5

44 100 MeV100 MeV 1767617676 53775377 22 5.38e-25.38e-2

1 GeV1 GeV 97399739 1414 00 1.4e-41.4e-4

55 100 MeV100 MeV 9082690826 41484148 1313 4.15e-24.15e-2

1 GeV1 GeV 8893888938 2727 33 2.7e-42.7e-4N=10N=1055

SpaC

al

SpaC

al

CK

MC

al

CK

MC

al

20/04/200520/04/2005 E.Leonardi, M.Serra, P.Valente - NA48/3 Veto Group MeetingE.Leonardi, M.Serra, P.Valente - NA48/3 Veto Group Meeting 2121N=10N=1055

11

22

33

4 54 5

large large




0+10+1


11 50 MeV50 MeV 9508095080 9508095080 9429494294 9.51e-19.51e-1

100 MeV100 MeV 9420694206 9420694206 9375193751 9.42e-19.42e-1

22 50 MeV50 MeV 1640516405 1640516405 1270112701 1.64e-11.64e-1

100 MeV100 MeV 1287512875 1287512875 95219521 1.29e-11.29e-1

33 50 MeV50 MeV 11001100 11001100 986986 1.1e-21.1e-2

100 MeV100 MeV 549549 549549 545545 5.49e-35.49e-3

44 50 MeV50 MeV 452452 452452 443443 4.52e-34.52e-3

100 MeV100 MeV 168168 168168 168168 1.68e-31.68e-3

55 50 MeV50 MeV 9965099650 9962999629 9044090440 9.96e-19.96e-1

100 MeV100 MeV 9977299772 9975999759 8976089760 9.98e-19.98e-1




0+10+1


11 50 MeV50 MeV 9873198731 9873198731 8893488934 9.87e-19.87e-1

100 MeV100 MeV 9895798957 9895798957 8777087770 9.9e-19.9e-1

22 50 MeV50 MeV 17621762 17621762 872872 1.76e-21.76e-2

100 MeV100 MeV 587587 587587 397397 5.87e-35.87e-3

33 50 MeV50 MeV 411411 411411 99 4.11e-34.11e-3

100 MeV100 MeV 5454 5454 22 5.4e-45.4e-4

44 50 MeV50 MeV 44 44 33 4e-54e-5

100 MeV100 MeV 00 00 00 00

55 50 MeV50 MeV 9944699446 9944099440 8676686766 9.94e-19.94e-1

100 MeV100 MeV 9960599605 9960099600 8553385533 9.96e-19.96e-1

SpaC

al

SpaC

al

CK

MC

al

CK

MC

al





0+10+1


11 100 MeV100 MeV 608608 608608 606606 6.08e-36.08e-3

1 GeV1 GeV 163163 163163 163163 1.63e-31.63e-3

22 100 MeV100 MeV 5353 4949 4444 4.9e-44.9e-4

1 GeV1 GeV 66 55 55 5e-55e-5

5 GeV5 GeV 88 88 77 8e-58e-5

33 100 MeV100 MeV 13291329 767767 00 7.67e-37.67e-3

1 GeV1 GeV 1010 11 00 1e-51e-5

5 GeV5 GeV 66 00 00 00

44 100 MeV100 MeV 57965796 11841184 33 1.18e-21.18e-2

1 GeV1 GeV 18911891 33 11 3e-53e-5

5 GeV5 GeV 16271627 00 00 00

55 100 MeV100 MeV 9974999749 22882288 4343 2.29e-22.29e-2

1 GeV1 GeV 9918199181 1414 99 1.4e-41.4e-4

5 GeV5 GeV 9761497614 88 66 8e-58e-5

11

22

33

4 54 5

small small

N=10N=1055

CK

MC

al

CK

MC

al

No significative change going from 1 to 5 GeVNo significative change going from 1 to 5 GeV


(Old)(Old) ToDo list ToDo list 1 1

More statistics More statistics from 10from 1044 to 10 to 1055 events events (OK?)(OK?) MMore test configurationsore test configurations

– More points on brick (bottom surface) More points on brick (bottom surface) OKOK– More realistic photon beam More realistic photon beam OK (for SpaCal)OK (for SpaCal)– Energy/Angular distribution? Energy/Angular distribution? Yes... still roughYes... still rough

Improve geometry/set-upImprove geometry/set-up– Use multiple phototubes per brick Use multiple phototubes per brick OKOK– Add separatAdd separationion layerlayer (adjacent bricks) (adjacent bricks) Not Not

yet…yet… Improve Improve light collection/light collection/digitization digitization

realistic parameters, realistic parameters, but not yet full digitizationbut not yet full digitization


(Old)(Old) ToDo list ToDo list 2 2

More complex setupMore complex setup– Add more bricks along the beam axis Add more bricks along the beam axis

2 rings now…2 rings now… (OK? easy to add more)(OK? easy to add more)– Use real brick shape Use real brick shape OK OK (for TileCal)(for TileCal)

Test different layouts to optimize detector Test different layouts to optimize detector efficiencyefficiency– Pb/Sci ratio Pb/Sci ratio Not yet…Not yet… but now readybut now ready– Aluminum thickness Aluminum thickness OK?OK?– Read-out Read-out Not yet…Not yet… but now readybut now ready


(New)(New) ToDo list 3 ToDo list 3

Use FLYO output as Use FLYO output as -beam gun -beam gun (z along tube, (z along tube, , E), E)

Compare with existing data?Compare with existing data?

status of geant4 veto simulation

Documents