1

update of reconstruction

Dorota Stefan CERN/NCBJ Robert Sulej NCBJ

protoDUNE SP measurement group 15 February 2016

D.Stefan R.Sulej

Outline

1. Chain of reconstruction for protoDUNE.

2. Angular track/shower reconstruction as an input for beam window discussion.

2 D.Stefan R.Sulej

3

Current works on reconstruction chain

Output from PMA (spatial reconstruction) in the same convention as Pandora:

• PFParticle hierarchy: 1 parent, many daughters

set in fcl: standard_pmalgtrackmaker.MakePFPs: true

• Each particle starts with 3D vertex object, even if isolated from other parts of event

set: standard_pmalgtrackmaker.RunVertexing: true (default in 35t and FD) if you need to see only results of non-trivial vertexing work, set in fcl: standard_pmalgtrackmaker.SaveOnlyBranchingVtx: true

• Minor difference: vertex in the interaction point, not on the first 3D point of a daughter, so e.g.

• single vertex object shared by many secondary tracks • vertex of delta on the muon trajectory even if first visible hits of electron are detached

Metadata of art assns used to save index and dx of each hit associated to tracks: ready to be used by calorimetry module and dE/dx-based PID module (implementing now calculation of dE/dx sequence to single PID value in LArSoft’s algorithms what we presented for e/gamma).

All MC-based tests code removed from the PMA module (uBoone asked for this).

indicated if a fcl option is default, otherwise you need to set if

ongoing work marked with blue

4

Current works on reconstruction chain

Particle hierarchy and full reconstruction chain options:

• Usually used up to now:

• Hits ClusterCrawler set of independent 2D clusters

• PMA: aggregate clusters, fit 3D tracks apply corrections to 2D find vertices, reoptimize track-vertex structure find primary track/vertex (new, developing)

standard_pmalgtrackmaker.CluMatchingAlg: 1 (default) standard_pmalgtrackmaker.FlipToBeam: true (naive, set upstream vtx as a root) standard_pmalgtrackmaker.AutoFlip_dQdx: true (can override naive flips) …todo neutrino vtx selection

• Use input from any pattern recognition that can save results in PFParticles:

• Hits e.g. Pandora, Cluster3D (never tried), eventually Wire-Cell fit 3D tracks to clusters associated with each PFP proceed (or not) with vertex finding and primary selection

standard_pmalgtrackmaker.HitModuleLabel: hits_used_by_pat_recog

standard_pmalgtrackmaker.ClusterModuleLabel: clu_and_pfp_label

standard_pmalgtrackmaker.CluMatchingAlg: 3

Pandora is tagging shower-like objects with electron PDG code, this is good, no need to fit trajectory there (or todo: implement fitting cascade beginning only). standard_pmalgtrackmaker.TrackingSkipPdg: [11] (now default) standard_pmalgtrackmaker.TrackingOnlyPdg: [0] (set pdg if you want fit only selected particles for any reason, and if they are tagged by pattern recognition of course)

• If 3D fit based corrections to pattern recognition output are desired, please use:

standard_pmalgtrackmaker.CluMatchingAlg: 2

5

simulation: kaon(480MeV/c) muon electron, first event in the file, three reco chains …and none is completely correct

ClusterCrawler + full PMA Pandora + PMA 3D fit only (asked to skip electrons)

Pandora 2D & 3D fits

2D

3D

6

1. We can setup fcl for testing various reco chains and produce some reco for people interested in using it.

• OK?

• Use simulations from Elizabeth or make new one with more recent ideas for beam orientation?

• Some measurement goals can be excersised now, some still need work on tools. We are (slowly) progressing on all topics, in contact with various groups (sorry for sometimes delayed replies).

2. Time to advertise ideas from the proposal.

• Some were only indicated briefly, while we had good (we hope) procedures for tests, calibrations, etc – no need to reinvent, of course welcome to improve.

• Work at a distance is OK, but it is good to start together. We are at CERN, maybe not much working space now, but we need reasons to ask for it ;) …so visitors are welcome, this would accelerate all the measurement works!

7

Beam window studies

P. Sala: talk on Neutrino Platform 11-02-2016

Question: what is the best angle to choose from the reconstruction point of view?

D.Stefan R.Sulej

3D reco 3D reco

2D raw plane Z

plane V

plane U

XZ

YZ

f

2GeV/c proton

reco vertex

p and p+ 2GeV/c simulated to test effect of beam orientation • beam angle range: ±15° in XZ, ±10° in YZ • entry point smearing: ±10cm in XY, ±0.5cm in Z

8

Track reconstruction – reminder

D.Stefan R.Sulej

9 D. Stefan, R. Sulej, 22/10/2015

protons 2GeV/c

RMS: 2.8°

Initial direction of the primary beam particle (protons), and also secondary particles some effect for beam XZ angle < ~3-4°; angular resolution correlated to track length (#hits)

track length [cm]

beam XZ angle [deg]

df [deg]

df [deg]

df [deg]

#hits (all planes)

initial direction: secondaries initial direction: primary

initial direction: primary initial direction: primary

[deg]

10

Shower direction reconstruction

Simulation:

1. Various angular bins: 90o – 0o wrt wire plane.

2. Energy range: 200MeV – 2 GeV.

Reconstruction:

1. Hits.

2. 3D segment, two approaches:

• use all hits from a shower, fix one endpoint of segment at mc primary vertex and reconstruct direction using hits from all views,

• the same but use only hits from initial part of a shower (5cm).

coll ind2 ind1 D.Stefan R.Sulej

11

Shower direction reconstruction (all hits in a shower)

90o 40o 10o RMS: 6o RMS: 4o RMS: 6o

8o RMS: 20o 5o RMS: 21o 3o RMS: 18o

0o RMS: 20o

Direction reconstruction:

- Segment attached to MC primary vtx and optimized to hits from all views.

- Precision depends on shower angle w.r.t. wire plane.

vertical YZ 6o

D.Stefan R.Sulej cosine (angle between mc truth and reco dir)

cosine cosine cosine

cosine cosine cosine

12

Shower direction reconstruction (all and part of hits)

Horizontal: 12o

vtx-reco: RMS: 8o

Horizontal: 8o

vtx-reco: RMS: 10o

use hits of initial part only (5cm from each view):

Horizontal: 3o

vtx-reco: RMS: 13o

use hits of full shower:

Horizontal: 12o

vtx-reco: RMS: 6o

• Vtx reco means that we know exactly the starting point of a shower: direction of segment is reconstructed being attached to this vertex.

• Reco means that segment is reconstructed without being attached to any vertex.

vertical YZ 11o

Horizontal: 8o

vtx-reco: RMS: 12o

Horizontal: 3o

vtx-reco: RMS: 17o

D.Stefan R.Sulej

cosine cosine cosine

cosine cosine cosine

13

Conclusion

• Various reconstruction chains can be now processed to start work on measurements which are in the proposal.

• Beam window: for the reconstruction preferably angle is 12o.

D.Stefan R.Sulej