TRIGGER PRIMITIVE RATES IN SIMULATION Karl Warburton Iowa State - - PowerPoint PPT Presentation

trigger primitive rates in simulation
SMART_READER_LITE
LIVE PREVIEW

TRIGGER PRIMITIVE RATES IN SIMULATION Karl Warburton Iowa State - - PowerPoint PPT Presentation

TRIGGER PRIMITIVE RATES IN SIMULATION Karl Warburton Iowa State University 20 th June 2019 1 MOTIVATION FOR LOOKING INTO THIS At a meeting on 15 th October, Phil showed some Trigger Primitive rates from ProtoDUNE data. Slides here, and some


slide-1
SLIDE 1

TRIGGER PRIMITIVE RATES IN SIMULATION

Karl Warburton Iowa State University 20th June 2019

1

slide-2
SLIDE 2

MOTIVATION FOR LOOKING INTO THIS

  • At a meeting on 15th October, Phil showed some Trigger Primitive rates from ProtoDUNE data.
  • Slides here, and some screenshots from this presentation are shown below.
  • It was noted that it would be interesting to see how this compares with simulation.
  • Unfortunately, the most similar samples are in the 1x2x6 geometry, not ProtoDUNE, so the plots that I’m about to show have to

be taken with a slight pinch of salt.

2

slide-3
SLIDE 3

WHAT I’M SHOWING TODAY / TOOLS USED

  • Apologies for a couple of very terse slides…
  • Using the “Tools for

Trigger Primitive Studies” framework, which is described in detail here.

  • Assuming that this is the most up-to-date version of code. I

am using dunetpc released on 12th April (v08_15_00).

  • Used the default noise levels, and the default trigger

primitive algorithm finders.

  • I did add a number of different trigger primitive threshold

modules though (seen right).

  • I don’t know what ADC values the different thresholds

correspond to though.

  • Did anyone come up with a conversion matrix for them?

3

  • In case anyone wants to look at any of the work

presented, all of my code resides in the various subdirectories contained in

/dune/app/users/karlwarb/LArDevelop/WorkArea_DAQ/TrigPrim_ProtoDUNE/Simulation-1x2x6/

slide-4
SLIDE 4

WHAT I’M SHOWING TODAY / TOOLS USED

  • Used outputs of samples the DAQ group requested in

MCC11, here.

  • Specifically the MARLEY sample (containing radio bkg) and

the 5 times Radio Bkg only sample.

  • I’m only showing plots from the baseline radiological rate

today, though I can quickly spin up the other sample.

  • The FlatTrees that are written out by theTrigPrim Ana

module are very large, so this represents only 10% of the full stats, but the flat tree is still 80 GB!

  • Could very well reduce the amount of info stored so that I

can process the full sample.

4

  • The TrigPrimAna module identifies the generator

responsible for each Trigger Primitive.

  • Be that noise, a radiological mode, or MARLEY, etc.
  • Can therefore “just” run the module, and then after

getting the total number of TrigPrim’s per event, calculate an expected rate for continuous running.

  • Not showing any MARLEY numbers here, but NB:
  • Sample has 3 interactions per event, so need to scale.
  • Obviously not going to be having SN interactions in every

drift window, so it’s more an indication of what we can expect if a SN does occur.

slide-5
SLIDE 5

NUMBER OF ARGON 39 TRIGGER PRIMITIVES FOR DIFF THRESHOLDS

  • Forgot a legend, but;
  • The lowest threshold (800) is black, mean of around 70.
  • The highest threshold (3200) is reddy-brown, mean < 10.
  • Argon39 is by far the highest number of Trigger Prims, and

we can’t really tell much from this plot, so won’t be showing any more of these plots.

  • All you can tell is that the number of trig prims is reduced

quite sharply by an increasing threshold.

  • T
  • p Plot isTPC 0, Bottom Plot isTPC 1.

5

20 40 60 80 100 120 140 Number of Hits 20 40 60 80 100

3

10 ´ Number of events

Hits_Per_TPC_1_Ar39 Hits_Per_TPC_1_Ar39

20 40 60 80 100 120 Number of Hits 20 40 60 80 100

3

10 ´ Number of events

Hits_Per_TPC_0_Ar39 Hits_Per_TPC_0_Ar39

slide-6
SLIDE 6

NUMBER OF ARGON 39 TRIGGER PRIMITIVES FOR DIFF THRESHOLDS

  • Forgot a legend, but;
  • The lowest threshold (800) is black, mean of around 1700.
  • The highest threshold (3200) is reddy-brown, mean ~30.
  • Again, can’t really tell much other than the obvious,

however the separation is quite cool.

  • Note how large the variation in NTrigPrims is for the lower

thresholds (the RHS of the plot).

  • Hints as to later, doubling the threshold from 800 (black) to

1600 (green) reduces the rate by a factor of ~4.

6

200 400 600 800 1000 1200 1400 1600 1800 2000 Number of Hits 10000 20000 30000 40000 50000 60000 70000 80000 Number of events

Hits_On_All_TPCs_Ar39 Hits_On_All_TPCs_Ar39

2

10

3

10 Number of Hits 1 10

2

10

3

10

4

10 Number of events

Hits_On_All_TPCs_Ar39

slide-7
SLIDE 7

CALCULATING THE TRIGGER PRIMITIVE RATE

  • No error bars on this plot, as taking the RMS of

the distribution seems to give too large errors (at least aesthetically).

  • Note that this is the whole detector, so divide by

24 to get the per TPC rates.

  • Ar39 dominates at all thresholds, noise is

removed quite quickly, and there were never any Trigger Primitives for the Polonium generator.

  • 1 event = 4492 ticks (500 ns) = 2.246e-3 s.
  • T
  • tal of 1e5 evts à 224.6 s in this study.

7

800 1000 1200 1400 1600 1800 2000 2400 2800 3200

Threshold

1

  • 10

1 10

2

10

3

10

4

10

5

10

6

10 Rate (Hz)

TriggerPrimRate_AllTPC

Ar39 Unknown APA CPA Neut Kryp Radon Ar42

TriggerPrimRate_AllTPC

slide-8
SLIDE 8

TRIGGER PRIMITIVE RATES PER TPC

8

800 1000 1200 1400 1600 1800 2000 2400 2800 3200

Threshold

1

  • 10

1 10

2

10

3

10

4

10 Rate (Hz)

TriggerPrimRate_TPC0

Ar39 Unknown APA CPA Neut Kryp Radon Ar42

TriggerPrimRate_TPC0

800 1000 1200 1400 1600 1800 2000 2400 2800 3200

Threshold

1

  • 10

1 10

2

10

3

10

4

10 Rate (Hz)

TriggerPrimRate_TPC1

Ar39 Unknown APA CPA Neut Kryp Radon Ar42

TriggerPrimRate_TPC1

800 1000 1200 1400 1600 1800 2000 2400 2800 3200

Threshold

1

  • 10

1 10

2

10

3

10

4

10 Rate (Hz)

TriggerPrimRate_TPC2

Ar39 Unknown APA CPA Neut Kryp Radon Ar42

TriggerPrimRate_TPC2

800 1000 1200 1400 1600 1800 2000 2400 2800 3200

Threshold

1

  • 10

1 10

2

10

3

10

4

10 Rate (Hz)

TriggerPrimRate_TPC3

Ar39 Unknown APA CPA Neut Kryp Radon Ar42

TriggerPrimRate_TPC3

  • Busy slide, but interest.

to have all together.

  • Most gen’s show

same pattern per TPC as for full TPC, and the rates are very stable.

  • However, the CPA

gen only makesTrig Prims in oddTPCs…

slide-9
SLIDE 9

SUMMARY / NEXT STEPS

  • Open questions:
  • Need to have a think about the CPA only making TrigPrims

in odd TPCs. Maybe the radiological group expects this?

  • Fix the error bars on the rate plots, will be interesting to

see what the event to event variation is around the mean.

  • How do the rates scale when radiological background is 5

times larger?

  • Niggles:
  • If we want to see more plots like Slides 5/6, I’ll have to stack

the hists or do a slightly better binning.

  • Do we expect Polonium to not produce any TrigPrims?
  • Conversion of Trigger Threshold into ADCs.

9

  • Not shown:
  • The average number of Marley Trig Prims per event (3

interactions) goes from 35 (Thr: 800) to 21 (Thr: 3200).

  • Raw numbers of Trig Prims per TPC for anything other than

Argon39, as above, I don’t think we can tell much though.

  • Large scale steps:
  • Josh expressed a need to reproduce Phil’s plots for “higher-

quality runs” as these were made for runs just after the detector turned off. Did that ever get done?

  • Don’t have any cosmics in here. Did we request any cosmic

samples? What generator do we use?

slide-10
SLIDE 10

BACKUP SLIDES

10