Updates on 0 rejection and cosmic ray simulation in ProtoDUNE-DP - - PowerPoint PPT Presentation

Updates on π0 rejection and cosmic ray simulation in ProtoDUNE-DP

Andrea Scarpelli APC - Paris Université Paris Diderot (VII) ProtoDUNE sim/reco meeting 14/06/2017

• Give some updates on preliminary results on π0 rejection based on gap detection in

dual-phase geometry

• Updates on some effort for generating cosmic events in the ProtoDUNE-DP geometry

both for the 6x6x6 m3 and the 3x1x1 m3 detectors.

Updates on π0 rejection in ProtoDUNE-DP

Very preliminary

π0 in Beam events

3

In beam events from Pions and (less often from Protons) there is a certain number

• f π0 appearing with a quite clear topology

Andrea Scarpelli

Updates on pi0 rejection and cosmic ray simulation in ProtoDUNE-DP

14/06/17

Raw linecluster Pi0 shower mip track mip track Pi0 shower Spurious em activity Spurious em activity Pi0 vetex Pi0 vetex

Strategy:

4 Andrea Scarpelli

Updates on pi0 rejection and cosmic ray simulation in ProtoDUNE-DP

14/06/17

Small em clusters Large em cluster not related to interesting vertex delta ray pi0 production point (truth) Track (m.i.p particle) Large Em-region, directions related to vertex

Identification of large region of electromagnetic activity in the event with a quite precise directionality towards a point 1.CNN output to label hits in 2D clusters (“linecluster”)

• 2. Size and purity (accumulated probability of all the hits in the cluster)
• 3. Direction with respect to vertex

4.Check if any hits of the cascade candidates cluster is around the considered vertex

Original particle gun: Charged pions @ 2GeV

Cluster selection on size and purity

5 Andrea Scarpelli

Updates on pi0 rejection and cosmic ray simulation in ProtoDUNE-DP

14/06/17

Large em cluster not related to interesting vertex pi0 production point (truth) Track (m.i.p particle) Large Em-region, directions related to vertex

Identification of large region of electromagnetic activity in the event with a quite precise directionality towards a point 1.CNN output to label hits in 2D clusters

• 2. Size and purity (accumulated probability of all the hits in the cluster)
• 3. Direction with respect to vertex

4.Check if any hits of the cascade candidates cluster is around the considered vertex

Cluster selection on direction

6 Andrea Scarpelli

Updates on pi0 rejection and cosmic ray simulation in ProtoDUNE-DP

14/06/17

Large em cluster not related to interesting vertex Track (m.i.p particle) Large Em-region, directions related to vertex

Identification of large region of electromagnetic activity in the event with a quite precise directionality towards a point 1.CNN output to label hits in 2D clusters

• 2. Size and purity (accumulated probability of all the hits in the cluster)
• 3. Direction with respect to vertex

4.Check if any hits of the cascade candidates cluster is around the considered vertex

Direction between cluster Center of mass and true pi0 vertex Direction between cluster start and cluster end

Gap identification

7 Andrea Scarpelli

Updates on pi0 rejection and cosmic ray simulation in ProtoDUNE-DP

14/06/17

Track (m.i.p particle)

Cascade candidates

Identification of large region of electromagnetic activity in the event with a quite precise directionality towards a point 1.CNN output to label hits in 2D clusters

• 2. Size and purity (accumulated probability of all the hits in the cluster)
• 3. Direction with respect to vertex

4.Check if any hits of the cascade candidates cluster is around the considered vertex

Gap detection around the vertex on cascade candidates only

Parameters choice

8 Andrea Scarpelli

Updates on pi0 rejection and cosmic ray simulation in ProtoDUNE-DP

14/06/17

Name Size Purity Direction (cosθ) Vertex distance (Wires)

Value 15 > 0.5 > 50∘ 2

• Parameters optimization on sigle pi0

events

• Dot product between the Cluster

direction (hit end - hit start) and Cluster Center of mass - Vertex direction —— size 10 —— size 15 —— size 20

Performances and comments

9 Andrea Scarpelli

Updates on pi0 rejection and cosmic ray simulation in ProtoDUNE-DP

14/06/17

Efficiency Overall ~0.74* Event-wise ~0.69*

* Almost no differences between the two views

• Keeping only events with at least one

pi0 decaying (2775 events over 5000 primary pions)

• Truth knowledge of the gap
• No cosmic overlapping or beam halo
• Mono-energetic particles (no hadronic

showers)

• A. Failure rate is mostly related to the size of the cluster (low energy showers)
• Identification of em-region to cluster smaller cluster with similar

directionality

• Use of different clustering algorithms (?)
• B. Selection on direction si quite high: better estimator for the cluster direction must be

envisaged (linear fit between the hits position in a cluster).

• C. Moving to 3D objects could be interesting and solve potential ambiguities

Cosmic ray generation in ProtoDUNE-DP

Remind on DP geometry status in LArSoft

11

Dual Phase geometry has drift direction parallel to gravity.

• Efforts are made to rotate the full simulation in LArSoft
• In the meantime, the rotation of the cosmic simulation flux can be

made to match the correct detector geometry:

• Naive approach, simply switched the x and y in larsim/larsim/

EventGenerator/CORSIKAGen_module.cc (coherent changes with this choice was made in all the file)

Andrea Scarpelli

Updates on pi0 rejection and cosmic ray simulation in ProtoDUNE-DP

14/06/17 3x1x1 Geometry (K. Fusshoeller) 6x6x6 Geometry (C. Alt)

Issue with geometry

12 Andrea Scarpelli

Updates on pi0 rejection and cosmic ray simulation in ProtoDUNE-DP

14/06/17

To achieve a reasonable result I had to make some changes in the geometry files for both protodune 6x6x6 and 3x1x1.

• Definition of the Worldbox from function geo::GeometryCore::WorldBox

was in conflict with the cryostat coordinates definition, in particular, there was an overlapping along z

• This caused weird asymmetries in the generated flux, and sometimes

some warnings not compatible with the full simulation

x, y z WorldBox Cryostat

geo::GeometryCore:: WorldBox geo::GeometryCore:: CryostatBoundaries World length as in the .gdml file

13 Andrea Scarpelli

Updates on pi0 rejection and cosmic ray simulation in ProtoDUNE-DP

14/06/17

• For the moment I made a choice of parameters in analogy with single

phase and put it into the .gdml file to make things working

• Some comments here?

<box name="World" lunit="cm" x=“7350.24" y=“7215.04" z=“7415.04"/> <box name="World" lunit="cm" x=“3483.54" y=“3447.94" z=“3723.64”/> <box name="World" lunit="cm" x="1350.24+2*3000" y="1215.04+2*3000" z="1415.04+2*3000"/> <box name="World" lunit="cm" x="483.54+2*0" y="447.94+2*0" z="723.64+2*0”/>

6x6x6 old 6x6x6 my correction 3x1x1 old 3x1x1 my correction Maybe too large!!

Issue with geometry

Event display - 3x1x1

14 Andrea Scarpelli

Updates on pi0 rejection and cosmic ray simulation in ProtoDUNE-DP

14/06/17

Raw data Orto3D and MCTruth particles

Event display - 6x6x6

15 Andrea Scarpelli

Updates on pi0 rejection and cosmic ray simulation in ProtoDUNE-DP

14/06/17

Raw data Orto3D with Truth particles track

Angular distribution - 3x1x1

16 Andrea Scarpelli

Updates on pi0 rejection and cosmic ray simulation in ProtoDUNE-DP

14/06/17

Rotated geometry: cosθ = -Px/|P| Φ = tg-1(Pz/Py) Non rotated geometry: cosθ = -Py/|P| Φ = tg-1(Pz/Px) Direction of the primary cr particle Polar angle distribution

— Rotated — Non Rotated

Just a first experiment: not matching with the real trigger on the 3x1x1, still some parameters to be tuned

Angular distribution - 6x6x6

17 Andrea Scarpelli

Updates on pi0 rejection and cosmic ray simulation in ProtoDUNE-DP

14/06/17

Rotated geometry: cosθ = -Px/|P| Φ = tg-1(Pz/Py) Non rotated geometry: cosθ = -Py/|P| Φ = tg-1(Pz/Px) Direction of the primary cr particle Polar angle distribution

The new rotated coordinates seem in agreement with the expected angular distribution:

• Match with the real trigger with the beam spill
• push module and .fcl in the repository and prepare a sample of cosmics to store on eos before the

end of June

— Rotated — Non Rotated

18 Andrea Scarpelli

Pion Rejection

✓

Discrete efficiency is obtained from a selection based on clusters size and direction

• Parameters optimization needed
• Reco vertex from mip tracks and sample with more particles
• Move to 3D objects

Cosmic rays sim

✓

Rotated the cosmic generation

✓

Applicable for 6x6x6

• Trigger on 3x1x1
• Push module and generation of a small sample
• Goal: have this in place before the LArSoft Hackatron days at the end
• f June

Updates on pi0 rejection and cosmic ray simulation in ProtoDUNE-DP

14/06/17

Conclusion and outlooks

Thank you!