Update Upd te o on the n the dual dual phase phase light ght - - PowerPoint PPT Presentation

▶
update upd te o on the n the dual dual phase phase light
SMART_READER_LITE
LIVE PREVIEW

Update Upd te o on the n the dual dual phase phase light ght - - PowerPoint PPT Presentation

Oct 31, 2023 263 likes •394 views

Update Upd te o on the n the dual dual phase phase light ght simulati si ation i on in l n larsoft arsoft J. Soto - 1 / 12 - 12 - | 09 09/10 10/18 Some activities ongoing Towards the TDR: Working on generating a sample

slide-1
SLIDE 1

| 09 09/10 10/18

  • J. Soto
  • 1 / 12

12 -

Upd Update te o

  • n the

n the dual dual phase phase light ght si simulati ation i

  • n in l

n larsoft arsoft

slide-2
SLIDE 2

| 09 09/10 10/18 18

  • J. Soto
  • 2 / 12

12 -

Some activities ongoing

  • Towards the TDR:

– Working on generating a sample of events to work

with.

  • Improving the simulation:

– S1 photon library with timing. – S2 photon library.

slide-3
SLIDE 3

| 09 09/10 10/18 18

  • J. Soto
  • 3 / 12

12 -

Goal of event samples

Including with charge and light:

  • 1. Radiological Bkg
  • 2. Mono energetic SN nu (energy scan)
  • 3. Mono energetic e- (energy scan)
  • 4. 2 +1
  • 5. 3 + 1
  • 6. NDK (channel p->nubarK+)
  • 7. 6 + 1

✔ ✔

slide-4
SLIDE 4

| 09 09/10 10/18 18

  • J. Soto
  • 4 / 12

12 -

Radiological Background

  • Thanks to Jason Stock!
  • We adapted the parameters needed to the dual phase

geometry.

  • Background is generated randomly in a window of 2

times the drift time.

  • It includes all argon inside the cryostat (liquid and gas).
  • Charge and light is included.
slide-5
SLIDE 5

| 09 09/10 10/18 18

  • J. Soto
  • 5 / 12

12 -

Main parameters for the samples

  • Main parameters:

– Charge:

  • Purity / tau = 3ms (0.9 less electrons in 12m drift)
  • Gain = 20
  • Birks Recomb (for 0,5kV/cm) ~0.7 for a MIP
  • We = 23.6eV (4e7 electrons/GeV)
  • dE/dx ~ 2.12 MeV/cm (tipically for a MIP)
  • No noise included!!!.

– Light

  • Light yield of 24k photons/MeV.
  • For EM saturations:

fNumScintPhotons = VisibleEnergyDeposition * scintYield

– http://www.apc.univ-paris7.fr/~franco/g4doxy/html/classG4EmSaturation.htm

l#4c7405e2bb881dc85dc9c958aec0ad67

  • Gain of 1.e7, 1ms readout window, 250MHz sampling, 50kHz DC,

0.12QE including TPB.

slide-6
SLIDE 6

| 09 09/10 10/18 18

  • J. Soto
  • 6 / 12

12 -

Benchmark

1 radiological evt Memory Peak (Mb) CPU time (s) Output file size (Mb) gen 1300 5 G4 (no S1prop) 6300 1500 detsim 6700 1700 reco 7500 150 440Mb ana 1800 1 ~200Mb 1 SN evt Memory Peak (Mb) CPU time (s) Output file size gen 1000 0.7 G4 (no S1prop) 1200 6 detsim 2800 1200 reco 1100 2 0.5M ana 1900 3 30k

slide-7
SLIDE 7

| 09 09/10 10/18 18

  • J. Soto
  • 7 / 12

12 -

S1 extended photon library

  • Using cryostat boundaries (we want to generate light all over the

volume), also outside the TPC active volume.

  • 15x15x62 voxels (1m3 voxels).
  • 1M photons per voxel.
  • Including timing (extended PL).
  • Photon library already in pardata.
slide-8
SLIDE 8

| 09 09/10 10/18 18

  • J. Soto
  • 8 / 12

12 -

S1 extended photon library

  • Timing is exported using the SimPhotonCounter

module.

  • Then a root macro is used to do the Landau fit

to every pair voxel PMT.

  • I only do the fit if I get more than 50 photons.
  • Thus we only get a TF1 for 102249 pairs voxel-

PMT (~10%).

  • Then, to use the library with the timing we move

from 1.5kMb of RAM to 6.5kMb! High but still affordable.

  • Issues:

– No timing interpolation between voxels, and they are

too big.

– Timing range is hard coded (fixed in the next

release).

– If we are interpolating visibility between voxels, and

timing is active, code crash… (working in a solution).

  • Still some things to improve… So timing won’t

be included for now.

slide-9
SLIDE 9

| 09 09/10 10/18 18

  • J. Soto
  • 9 / 12

12 -

S2 photon library

  • Visibility in the top of the detector is very very low.
  • Usual code doesn’t work, since we need much more statistics of photons, to have a

reasonable value.

  • To solve this I run several libraries in parallel, and then merged them.
  • This was relatively fast: We need a lot of memory to simulate 3M photons, but most of them

got killed very fast (as they are generated very close to the lems, i.e. solid copper).

  • In this way, I generated a 18M photons/voxel photon library.
  • 1x12x60m voxels of 1m2.
  • Se still have too low statistics to do the timing fit!

All voxels! All voxels!

slide-10
SLIDE 10

| 09 09/10 10/18 18

  • J. Soto
  • 10 / 12

12 -

S2 photon library

  • A possible/fast solution would be to parametrize the

timing for all voxels and all PMTs and provide a general TF1 for all the S2:

  • 31k photons out of 18Mph/voxel * 720voxels:

– Global S2 visibility of 2e-6 (2 photons/Mph) in all PMTs. Landau fit does not work...

slide-11
SLIDE 11

| 09 09/10 10/18 18

  • J. Soto
  • 11 / 12

12 -

Next steps

  • Test the hybrid library:

– To do this we need a high statistics and high voxel

density photon library to work with:

  • Generating now a 1Mph/25cm3 voxel PL.
  • But it is taking too much time in the grid… jobs don’t get

to run… too long and heavy… Exploring other alternatives now.

slide-12
SLIDE 12

| 09 09/10 10/18 18

  • J. Soto
  • 12 / 12

12 -

Backup - PL Benchmark

1evt, 6GeV electron # photons # voxels Voxel size Mb CPU Time (s) Lib 300k (only tpc active) 300k 16x16x81 0.75m3 1400 18 Lib 1M (all cryostat) 1M 15x15x62 1m3 1400 16 Lib 1M with timing (all cryostat) 1M 15x15x62 1m3 6400 34