SBND (UK)* software and physics report Dom Brailsford on behalf of - - PowerPoint PPT Presentation

sbnd uk software and physics report
SMART_READER_LITE
LIVE PREVIEW

SBND (UK)* software and physics report Dom Brailsford on behalf of - - PowerPoint PPT Presentation

SBND (UK)* software and physics report Dom Brailsford on behalf of SBND UK DUNE UK meeting at Manchester 06/03/19 Covers topics not already presented at this meeting* The Short Baseline Near Detector (SBND) Serves as near detector of the


slide-1
SLIDE 1

SBND (UK)* software and physics report

Dom Brailsford on behalf of SBND UK DUNE UK meeting at Manchester 06/03/19

Covers topics not already presented at this meeting*

slide-2
SLIDE 2

The Short Baseline Near Detector (SBND)

  • LArTPC
  • 112 t fiducial volume
  • 110 m from BNB target
  • Serves as near detector of the

short baseline neutrino program

  • Sterile neutrino search

2

slide-3
SLIDE 3

The Short Baseline Near Detector (SBND)

Motivation

νμ CC 0π 3p signal event in the SBND MC sample μ- p p p

Neutrino vertex Time, TDC Wire number νμ CC inclusive event rate breakdown, not stacked νμ CC exclusive event rate breakdown, stacked Proton multiplicity in the true νμ CC final state

SBND will observe huge statistics at bubble-chamber resolutions (3 mm). Ensuring the software is capable of detecting particles in such detail is crucial for

  • Neutrino energy reconstruction

○ Low momentum particles ○ High multiplicities

  • Neutrino interaction studies

○ 2p-2h cross-section measurement

  • Detecting and analysing rare final states

SBND

GENIE v02.12.10, Default+MEC

  • SBND will observe millions of

neutrino interactions

  • Measurements of many exclusive

cross section channels

  • BSM physics
  • R. Jones
  • R. Jones

3

slide-4
SLIDE 4

SBND physics organisation

  • Geared towards SBND’s own physics program i.e. what can

SBND measure on its own

  • Cross sections
  • BSM physics
  • Detector measurements
  • Organised via the Physics and Analysis Tools (PAT) group
  • 2/4 convenors are UK-based
  • Andrzej Szelc
  • Costas Andreopoulos
  • Heavy UK involvement

4

slide-5
SLIDE 5

10 sub-groups

SBN physics

  • rganisation
  • Shift to more formal SBN

group structure

  • 10 groups
  • 1 SBND and 1 ICARUS

convener per group

  • Aim: Unify sim./reco./

analysis for the multi- detector physics program

5

slide-6
SLIDE 6

10 sub-groups

SBN physics

  • rganisation

UK leadership Strong UK contribution

6

slide-7
SLIDE 7

MC production

  • 3 production campaigns runs between 2018-2019
  • We took a production processing hiatus in mid-2018 to switch from

project.py-based productions to POMS-based productions

  • Manual submission -> automatic submission
  • Manual job recovery -> automatic recovery
  • Complete handling/bookkeeping of Monte Carlo file metadata
  • Automatic transfer of files to tape for permanent storage
  • Processing on the open science grid as well as FNAL
  • SBND is currently running its largest ever production campaign
  • ~1,000,000 events for physics analysis
  • D. Brailsford

7

slide-8
SLIDE 8

SBND-POMS production workflow

POMS CPU running fife-wrap SAM FTS box Tape

Initial jobs submitted File metadata declared Output files copied Files copied using SAM bookkeeping

Jobs auto-submitted using SAM datasets from tape

8

Jobs report finished

  • D. Brailsford
slide-9
SLIDE 9

BNB timing structure in LArSoft

Time structure of our current neutrino simulation The actual time structure of the BNB

  • Modelling the BNB timing structure in

LArSoft would allow us to study:

  • Detector timing resolution
  • BSM physics
  • First implementation feature complete
  • Timing structure correctly propagates

through the entire simulation change

  • Now finalising before making available to

the community

  • A. Ezeribe

9

slide-10
SLIDE 10

Cosmic Ray Tagger (CRT) system reconstruction

Cosmic ray tagger system reconstruction

The SBND TPC is almost completely covered by a number of cosmic ray taggers formed of plastic scintillator. They can be used to reconstruct CRT hits (the position and time

  • f particle-CRT intersections) and CRT tracks (the trajectory of

through going particles). CRT reconstruction CRT-TPC matching TPC reconstructed tracks can be matched to CRT hits by projecting their ends on to the CRT taggers. They can be matched to CRT tracks by comparing angles and start/end positions.

  • T. Brooks

10

slide-11
SLIDE 11

CRT cosmic background removal

CRT cosmic background removal

Fiducial volume:

  • Remove tracks with start and end points

within 10 cm of TPC walls. TPC topology cosmic ID:

  • Tracks reconstructed outside of TPC.
  • Tracks which start outside of fiducial volume

and stop inside it.

  • Tracks which match CRT tracks.

t0 tagging cosmic ID:

  • t0 from stitching tracks across the CPA.
  • t0 from tracks which match CRT hits.
  • t0 from matching tracks with APA crossing

points. Background removal Sample of 5,000 TPC contained neutrino events with a corsika cosmic overlay. Only using CRT matching and basic TPC information (no light or Pandora reconstruction). Able to go from a 1:14 neutrino muon to cosmic muon ratio to 2:1. Preliminary results

  • T. Brooks

11

slide-12
SLIDE 12

Shower reconstruction validation

  • SBND is currently investigating

various combinatorics of reconstruction algorithms

  • Blurred cluster (M. Wallbank)
  • EMShower (M. Wallbank)
  • Pandora (Team pandora)
  • SBND has developed a shower

validation module

  • Easily compare different

reconstruction algorithms

  • Easily benchmark reconstruction

performance using metrics emshowerNew: Pandora + EMShower pandoraShower: Full pandora reco. emshowerBLUR: Blur. cluster + EMShower

  • D. Barker

12

slide-13
SLIDE 13

Example electron reconstruction

emshowerNew: Pandora + EMShower pandoraShower: Full pandora reco. emshowerBLUR: Blur. cluster + EMShower

Shower reconstruction validation

  • D. Barker

13

slide-14
SLIDE 14
  • Aspects of pandora are

tune-able by end users

  • Track-shower separation
  • Vertex identification
  • Currently running a phase
  • f exploratory track-shower

separation tunings

  • Focussed on shower

segmentation for electron particle gun

  • We are now ramping up a

metric-based tuning approach with a more realistic topology

Before After

  • E. Tyley

Reconstruction tuning

14

slide-15
SLIDE 15

Recombination

  • Aim is to understand how

recombination affects calorimetric reconstruction in showers

  • Truth-based study using particle gun

electrons and muons

  • A flat recombination factor for electrons

is appropriate, which differs to muons

  • Study extended to hit-based calorimetry

using the same samples

  • Sums reco. energy for every hit
  • Reconstructed energy as a function of

true electron momentum shows a modest dependence

  • E. Tyley

15

slide-16
SLIDE 16

PID via track follow-down

  • Calorimetry-based PID struggles to

distinguish similar mass particles such as muons and pions

  • Other options are needed to separate

muons and pions

  • ~70% of muons are captured before

decay - simple straight line track topology

  • Decaying muons can be tagged by the

Michel electron

  • Pions are more inclined to scatter in

the TPC

  • Goal is to develop a pion/muon separator

to improve the CC νμ1π selection

Muon vs pion ID in LArTPC detectors

Need to define a PID method for muons and pions based on a track follow down procedure

  • Calorimetry PID methods are not good to

distinguish particles with similar mass

  • However, muons and pions leave different

signatures in the TPC ○ 73% of the times the muon is absorbed - straight line track ○ Michel electrons can be used to tag muon tracks ○ Pions may interact with the environment by absorption, elastic or inelastic scattering

  • The final goal is to apply this to 𝜉𝜈CC1π+ events to

improve the analysis

1 J.Tena Vidal - University of Liverpool

  • J. Tena-Vidal

16

slide-17
SLIDE 17
  • Pandora provides a parent/child hierarchy
  • Method initially developed by

MicroBooNE for Michel electron searches

  • Based on calculating Pearson coefficient

for hits within a search window

PID via track follow-down

  • J. Tena-Vidal

17

slide-18
SLIDE 18

CC νμ selection

π π p

Fiducial volume of the TPC

p n ν

μ

μ n

  • Use calorimetric information to

distinguish protons from π & μ

  • Along with geometrical

information for MIPs

  • Initial priority on

μ tagging

  • Small amount of

focus on protons Selection

  • ptimised for

μ & proton purity

Total BNB-only events with a single contained, reconstructed neutrino vertex 63,830

True vertex also contained 96.3% Maximum 1 escaping track 99.9% Exactly 1 escaping track 5.5%

Of these, only the true muon escapes 95.9%

2

5.3% of events are basically free, guaranteed muons!

Topological selection

↓ True / Reco → CC Inc. CC 0π CC 0π

39,100 32,650

CC 1π

8,386 3,218

CC Other

658 70

NC

2,967 2,130

Efficiency

92.0% 76.9%

Purity

94.2% 85.8%

Single interaction νμ BNB-only

Purity: Signaltopology / Total selectedtopology Efficiency: Signaltopology / Total truetopology

  • R. Jones

18

slide-19
SLIDE 19

π π p

Fiducial volume of the TPC

p n ν

μ

μ n

  • Use calorimetric information to

distinguish protons from π & μ

  • Along with geometrical

information for MIPs

  • Initial priority on

μ tagging

  • Small amount of

focus on protons Selection

  • ptimised for

μ & proton purity

Total BNB-only events with a single contained, reconstructed neutrino vertex 63,830

True vertex also contained 96.3% Maximum 1 escaping track 99.9% Exactly 1 escaping track 5.5%

Of these, only the true muon escapes 95.9%

3

5.3% of events are basically free, guaranteed muons!

Topological selection

↓ True / Reco → CC Inc. CC 0π CC 0π

1,831 1,471

CC 1π

404 174

CC Other

30 6

NC

190 142

Efficiency

92.1% (+0.1%) 74.1% (-2.8%)

Purity

92.3% (-1.9%) 82.0% (-3.8%)

νμ BNB+cosmic overlay - cheated neutrino ID

Low statistics makes it difficult to judge how the inclusion of cosmics affects the performance of the selection. A 5-10% efficiency loss is expected: Table 4: arXiv:1708.03135

CC νμ selection

  • R. Jones

19

π π p

Fiducial volume of the TPC

p n ν

μ

μ n

  • Use calorimetric information to

distinguish protons from π & μ

  • Along with geometrical

information for MIPs

  • Initial priority on

μ tagging

  • Small amount of

focus on protons Selection

  • ptimised for

μ & proton purity

Total BNB-only events with a single contained, reconstructed neutrino vertex 63,830

True vertex also contained 96.3% Maximum 1 escaping track 99.9% Exactly 1 escaping track 5.5%

Of these, only the true muon escapes 95.9%

3

5.3% of events are basically free, guaranteed muons!

Topological selection

↓ True / Reco → CC Inc. CC 0π CC 0π

1,831 1,471

CC 1π

404 174

CC Other

30 6

NC

190 142

Efficiency

92.1% (+0.1%) 74.1% (-2.8%)

Purity

92.3% (-1.9%) 82.0% (-3.8%)

νμ BNB+cosmic overlay - cheated neutrino ID

Low statistics makes it difficult to judge how the inclusion of cosmics affects the performance of the selection. A 5-10% efficiency loss is expected: Table 4: arXiv:1708.03135

slide-20
SLIDE 20
  • t
  • and
  • DM Production Channels:

DM Scattering Channels: BNB Target SBND 110m χ χϯ p This analysis Focuses on the electron scattering channel.

18 40𝐵𝑠

χ e

  • Fixed target neutrino experiments

could be a viable means of detector light dark matter (DM) in the sub-GeV mass range

  • Benchmark model of vector portal

light dark matter used

  • DM is produced at the BNB target

and propagates to SBND. BdNMC used as Monte Carlo DM event generator

  • Background events are due to BNB

neutrinos and cosmic rays. GENIE generates neutrino background and CORSIKA generates cosmic rays

Light dark matter search

  • E. Sandford and C. Hasnip

20

slide-21
SLIDE 21

Selection Cuts:

  • Apply selection cuts similar to searches for intrinsic

electron neutrino - scattering.

  • KE*𝜄 cut applied at each DM mass.

Prel relim imin inary ry SB SBND Se Sensi sitivity

  • SBND sensitivity calculated as a function of

parameter Y, a function of coupling and dark portal particle masses.

  • This is a work in progress. Running the BNB in

beam dump mode could help increase sensitivity. Timing cut:

  • Expect DM to take longer to

travel to SBND than neutrinos.

Red = Neutrino Blue = DM 𝑛χ= 0.3 GeV

Preliminary

  • E. Sandford and C. Hasnip

Light dark matter search

21

slide-22
SLIDE 22

Summary

  • SBND is a LArTPC that will serve as the near detector in the short

baseline neutrino oscillation program as well as having its own rich physics program

  • A more formal SBN-based working group structure has formed over

the past year

  • SBND-UK plays a very key role in leadership and active contribution

to that working group structure

  • SBND’s software and physics program is undergoing very active

development

  • Monte Carlo production
  • Simulation
  • Reconstruction
  • Physics analysis

22

slide-23
SLIDE 23

Backups

slide-24
SLIDE 24

Single Plane Semantic Segmentation

  • Pixel-level tagging of cosmic/neutrinos
  • Convolutional Neural Network
  • U-Resnet architecture (developed by MicroBooNE)
  • 95% success rate in pixel tagging of νe CC + cosmic MC
  • See talk by S. Dennis yesterday
  • Single Plane Semantic Segmentation using U-Resnet for Cosmic Tagging

○ nueCC Dataset Training (Single CC interation with cosmics) ■ 95% Efficiency in tagging event pixels out of cosmics and background

Work Summary

Meeting Jaggar Henzerling 1 Segmentation Image (Energy Deposition) Prediction 100 cm 100 cm

CC interaction CC interaction

  • J. Henzerling

24