Updates of DISC Software Mustafa Schmidt, Simon Bodenschatz, Michael - - PowerPoint PPT Presentation

▶

Jul 28, 2023 192 likes •392 views

Updates of DISC Software Mustafa Schmidt, Simon Bodenschatz, Michael D uren, Erik Etzelm uller, Klaus F ohl, Avetik Hayrapetyan, Oliver Merle, Julian Rieke September 13, 2016 0/18 Mustafa Schmidt Updates of DISC Software 0 / 18

SLIDE 1

0/18

Updates of DISC Software

Mustafa Schmidt, Simon Bodenschatz, Michael D¨ uren, Erik Etzelm¨ uller, Klaus F¨

hl, Avetik Hayrapetyan, Oliver Merle,

Julian Rieke September 13, 2016

Mustafa Schmidt Updates of DISC Software 0 / 18

SLIDE 2

1/18

Full Simulation Chain

PndDisc() ${PANDAROOT}/disc/ ${PANDAROOT}/pnddata/DiscData/ Simulation PndDiscTaskDigitization() Digitization PndDiscTaskReconstruction() Reconstruction PndDiscTaskPID() PID PndDiscParticleMCPoint()

Information about charged particles

PndDiscSensorMCPoint()

Information about sensor hits

PndDiscDigitizedHit()

Sensor information (pixel number, TDC time etc.)

PndDiscReconResult()

Theoretical hit pattern

PndDiscCherenkovResult()

Reconstructed Cherenkov Angle

PndDiscPID()

Log(Likelihood) values for π, K, p

PhotoSensor()

Sensor Description Mustafa Schmidt Updates of DISC Software 1 / 18

SLIDE 3

2/18

Simulation Output

PndDiscSensorMCPoint() Photon Hit Position: fX, fY, fZ Photon Entering Position: photon entering pos (fX, fY, fZ) Photon Hit Momentum: fPx, fPy, fPz Photon Entering Momentum: photon entering mom (fX, fY, fZ) Reflection Angle: internal reflecting angle PndDiscParticleMCPoint() Particle Incoming Position: fX, fY, fZ Particle Outgoing Position: pos out (fX, fY, fZ) Particle Incoming Momentum: fPx, fPy, fPz Particle Outgoing Momentum: mom out (fX, fY, fZ) PDG Code: pdgCode Particle Charge: charge Particle Mass: mass

Mustafa Schmidt Updates of DISC Software 2 / 18

SLIDE 4

3/18

Digitization, Reconstruction, PID Output

PndDiscTaskDigitization() Detector ID: detector id Readout ID: readout id Sensor ID: sensor id Pixel Number: pixel number Pixel Position: pixel position TDC Time: tdc time PndDiscReconResult() Sensor ID: sensor Pixe ID: pixel Time: time Hypothesis: hypothesis PndDiscCherenkovResult() Cherenkov angle: cherenkov angle PndDiscPID() Likelihood π: loglikepion Likelihood K: loglikekaon Likelihood p: loglikeproton

Mustafa Schmidt Updates of DISC Software 3 / 18

SLIDE 5

4/18

Adding Event Display

Mustafa Schmidt Updates of DISC Software 4 / 18

SLIDE 6

5/18

Fixed Problems

Problem in ROOT::Math::Interpolation::kLINEAR → PDE inside digitization could not be used until now Library libMathMore added to the dependencies in CMakeLists.txt Other Problems: DEPENDENCIES was written wrong in CMakeLists.txt

[nm] λ

300 400 500 600 700

#Efficiency [%]

5 10 15 20 25 MCP Efficiency Mustafa Schmidt Updates of DISC Software 5 / 18

SLIDE 7

6/18

Changes in Track Reconstruction

Possibility to use Monte-Carlo tracks or reconstructed tracks with GenFit: PndDiscTaskReconstruction :: UseTrueTracks(kTRUE/kFALSE) Getting position and momentum information with: TClonesArray *tclarr_track (TClonesArray*) io_manager->GetObject("SttMvdGemGenTrack"); PndTrack *track = (PndTrack*)tclarr_track->At(i); FairTrackParP par = track->GetParamLast(); Calculating position on radiator disk assuming straight line:

xDisc

yDisc

= zTrack − zDisc

pz ·

Mustafa Schmidt

Updates of DISC Software 6 / 18

SLIDE 8

7/18

Resolution Studies for GenFit

Spatial and momentum resolution of Protons with p = 4 GeV/c and θ = 15◦

Resolution [cm] 0.15 − 0.1 − 0.05 − 0.05 0.1 0.15 0.2 Entries 10 20 30 40 50 60 70 80 90 Directions x y Spatial Resolution of PANDA Tracking System at Disc DIRC Resolution [mrad] 0.01 − 0.008 − 0.006 − 0.004 − 0.002 − 0.002 0.004 0.006 0.008 0.01 Entries 20 40 60 80 100

Polar Angle Resolution

Resolution [GeV/c] 1 − 0.8 − 0.6 − 0.4 − 0.2 − 0.2 0.4 0.6 0.8 1 Entries 20 40 60 80 100 120

Momentum Resolution of PANDA Tracking System

Mustafa Schmidt Updates of DISC Software 7 / 18

SLIDE 9

8/18

Adding Calibration Script

Script added for calibration of Focusing Elements using Monte-Carlo hits with random ϕ angles:

Pixel # 20 40 60 80 100 [rad] ϕ 0.35 0.4 0.45 0.5 0.55 0.6 0.65 0.7 Data Points Linear Fit

MC FEL Calibration

Values are used for the hit pattern calculation

Mustafa Schmidt Updates of DISC Software 8 / 18

SLIDE 10

9/18

Cherenkov Angle Reconstruction

Input Parameters xp, yp, θp, φp, tp Hit Pattern nFEL, z, t Mass Hypotheses mπ, mK , mp Cherenkov Angle θc Time Cut θ′

Photon Propagation sph, tph Removed Unphysical Angles θ′′

Hypothesis Cut θ′′′

; |θ′′

c − θc,theo| < θcthresh

Average Cherenkov Angle ¯ θc =

1 N

θ′′′

Theoretical Cherenkov Angle θc,theo Mustafa Schmidt Updates of DISC Software 9 / 18

SLIDE 11

10/18

Reconstructed Chernekov Angle

Particle momentum p = 3 GeV/c, polar angle θ = 15◦

[rad]

θ 0.7 0.72 0.74 0.76 0.78 0.8 0.82 0.84 0.86 0.88 0.9 Entries 10 20 30 40 50 Primary Particle Pion Kaon Proton

Cherenkov Angle Distribution

Mustafa Schmidt Updates of DISC Software 10 / 18

SLIDE 12

11/18

Hitpattern Matching

Input Parameters xp, yp, θp, φp, tp Mass Hypotheses mπ, mK , mp Hit Pattern nFEL, z, t Theoretical Hitpattern nFEL,theo, ztheo Photon Propagation sph, tph Spatial Cut |z − ztheo| < zthresh Time Cut |t − ttheo| < tthresh Likelihood Value Mustafa Schmidt Updates of DISC Software 11 / 18

SLIDE 13

12/18

Reconstruction Results

Hit pattern matching for particle momentum p = 3 GeV and polar angle θ = 15◦:

Sensor ID 5 10 15 20 25 30 Pixel Number 10 20 30 40 50 60 70 80 90 100 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1

Simulated Hitpattern

Sensor ID 5 10 15 20 25 30 35 Pixel Number 10 20 30 40 50 60 70 80 90 100

Predicted Hitpattern

Pion Kaon Proton

Used bandpass filter 360–465 nm including PDE

Mustafa Schmidt Updates of DISC Software 12 / 18

SLIDE 14

13/18

Particle Identification

Momentum p = 4 GeV/c, polar angle θ = 15◦

))-log(L(K)) π log(L( 60 − 40 − 20 − 20 40 60 Events 10 20 30 40 50 60 70 80 90 Separation Power: 2.64618 Primary Particles Pions Kaons Gaussian Fit

Likelihood Distribution

Mustafa Schmidt Updates of DISC Software 13 / 18

SLIDE 15

14/18

Conclusion & Outlook

PID in principle working in PandaRoot (results not sufficient until now) Cleaning of unnecesary code Better implementation of code into PandaRoot framework Improvement of reconstruction algorithm Testing PID with specific benchmark channel

Mustafa Schmidt Updates of DISC Software 14 / 18

SLIDE 16

15/18

Thank you very much for your attention!

Mustafa Schmidt Updates of DISC Software 15 / 18

SLIDE 17

16/18

Backup Slides

Mustafa Schmidt Updates of DISC Software 16 / 18

SLIDE 18

17/18

Theoretical Model

Angle Definitions: particle φrel αFEL Radiator Disk FEL Particle θc ϕ tan ϕ′ = tan ϕ cos αFEL ϕ′

Mustafa Schmidt Updates of DISC Software 17 / 18

SLIDE 19

18/18

Theoretical Model

Calculation of the Cherenkov angle: θc = arccos(sin θp cos φrel cos ϕ + cos θp sin ϕ) (1) θp: θ angle of particle φrel: angular difference between φ angle of particle and photon ϕ: Angle between total reflected photon and radiator disk surface Calculation of ϕ if θc is known: cos ϕ = A cos θc B ±

cos2 θp − cos2 θc

B +

A cos θc

with A = sin θp cos φrel and B = A2 + cos2 θp

Mustafa Schmidt Updates of DISC Software 18 / 18