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

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

Full Simulation Chain $ { PANDAROOT } /disc/ $ { PANDAROOT } /pnddata/DiscData/ PndDiscParticleMCPoint() Simulation Information about charged particles PndDisc() PndDiscSensorMCPoint() Information about sensor hits Digitization PhotoSensor() PndDiscDigitizedHit() PndDiscTaskDigitization() Sensor Description Sensor information (pixel number, TDC time etc.) PndDiscReconResult() Theoretical hit pattern Reconstruction PndDiscTaskReconstruction() PndDiscCherenkovResult() Reconstructed Cherenkov Angle PID PndDiscPID() PndDiscTaskPID() Log(Likelihood) values for π , K, p 1/18 Mustafa Schmidt Updates of DISC Software 1 / 18

Simulation Output PndDiscParticleMCPoint() PndDiscSensorMCPoint() Particle Incoming Position: Photon Hit Position: fX, fY, fZ fX, fY, fZ Particle Outgoing Position: Photon Entering Position: pos out (fX, fY, fZ) photon entering pos (fX, fY, fZ) Particle Incoming Momentum: Photon Hit Momentum: fPx, fPy, fPz fPx, fPy, fPz Particle Outgoing Photon Entering Momentum: Momentum: mom out (fX, fY, fZ) photon entering mom (fX, fY, fZ) PDG Code: pdgCode Reflection Angle: Particle Charge: charge internal reflecting angle Particle Mass: mass 2/18 Mustafa Schmidt Updates of DISC Software 2 / 18

Digitization, Reconstruction, PID Output PndDiscReconResult() PndDiscTaskDigitization() Sensor ID: sensor Detector ID: detector id Pixe ID: pixel Readout ID: readout id Time: time Sensor ID: Hypothesis: hypothesis sensor id Pixel Number: pixel number PndDiscCherenkovResult() Pixel Position: pixel position Cherenkov angle: cherenkov angle TDC Time: tdc time PndDiscPID() Likelihood π : loglikepion Likelihood K : loglikekaon Likelihood p : loglikeproton 3/18 Mustafa Schmidt Updates of DISC Software 3 / 18

Adding Event Display 4/18 Mustafa Schmidt Updates of DISC Software 4 / 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 MCP Efficiency 25 #Efficiency [%] 20 15 10 5 0 300 400 500 600 700 λ [nm] 5/18 Mustafa Schmidt Updates of DISC Software 5 / 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: � � � � x Disc = z Track − z Disc p x · y Disc p y p z 6/18 Mustafa Schmidt Updates of DISC Software 6 / 18

Resolution Studies for GenFit Spatial and momentum resolution of Protons with p = 4 GeV/c and θ = 15 ◦ Polar Angle Resolution Spatial Resolution of PANDA Tracking System at Disc DIRC Entries Entries 90 100 Directions 80 x 70 80 y 60 60 50 40 40 30 20 20 10 0 0 − 0.15 − 0.1 − 0.05 0 0.05 0.1 0.15 0.2 − 0.01 − 0.008 − 0.006 − 0.004 − 0.002 0 0.002 0.004 0.006 0.008 0.01 Resolution [cm] Resolution [mrad] Momentum Resolution of PANDA Tracking System Entries 120 100 80 60 40 20 0 − 1 − 0.8 − 0.6 − 0.4 − 0.2 0 0.2 0.4 0.6 0.8 1 Resolution [GeV/c] 7/18 Mustafa Schmidt Updates of DISC Software 7 / 18

Adding Calibration Script Script added for calibration of Focusing Elements using Monte-Carlo hits with random ϕ angles: MC FEL Calibration [rad] ϕ 0.7 Data Points 0.65 Linear Fit 0.6 0.55 0.5 0.45 0.4 0.35 0 20 40 60 80 100 Pixel # Values are used for the hit pattern calculation 8/18 Mustafa Schmidt Updates of DISC Software 8 / 18

Cherenkov Angle Reconstruction Input Parameters Hit Pattern Mass Hypotheses x p , y p , θ p , φ p , t p n FEL , z , t m π , m K , m p Photon Propagation Theoretical Cherenkov Angle Cherenkov Angle s ph , t ph θ c , theo θ c Time Cut θ ′ c Removed Unphysical Angles θ ′′ c Hypothesis Cut θ ′′′ ; | θ ′′ c − θ c , theo | < θ cthresh c Average Cherenkov Angle � ¯ 1 θ ′′′ θ c = c N 9/18 Mustafa Schmidt Updates of DISC Software 9 / 18

Reconstructed Chernekov Angle Particle momentum p = 3 GeV/c, polar angle θ = 15 ◦ Cherenkov Angle Distribution Entries 50 Primary Particle Pion Kaon 40 Proton 30 20 10 0 0.7 0.72 0.74 0.76 0.78 0.8 0.82 0.84 0.86 0.88 0.9 θ [rad] c 10/18 Mustafa Schmidt Updates of DISC Software 10 / 18

Hitpattern Matching Input Parameters Mass Hypotheses Hit Pattern x p , y p , θ p , φ p , t p m π , m K , m p n FEL , z , t Photon Propagation Theoretical Hitpattern s ph , t ph n FEL , theo , z theo Spatial Cut | z − z theo | < z thresh Time Cut | t − t theo | < t thresh Likelihood Value 11/18 Mustafa Schmidt Updates of DISC Software 11 / 18

Reconstruction Results Hit pattern matching for particle momentum p = 3 GeV and polar angle θ = 15 ◦ : Simulated Hitpattern Predicted Hitpattern 100 1 100 Pixel Number Pixel Number 90 0.9 90 Pion 80 0.8 80 Kaon Proton 70 0.7 70 60 0.6 60 50 0.5 50 40 0.4 40 30 0.3 30 20 0.2 20 10 0.1 10 0 0 0 0 5 10 15 20 25 30 5 10 15 20 25 30 35 Sensor ID Sensor ID Used bandpass filter 360–465 nm including PDE 12/18 Mustafa Schmidt Updates of DISC Software 12 / 18

Particle Identification Momentum p = 4 GeV/c, polar angle θ = 15 ◦ Likelihood Distribution Events 90 80 Primary Particles 70 Pions Separation Power: 2.64618 Kaons 60 Gaussian Fit 50 40 30 20 10 0 − − − 60 40 20 0 20 40 60 π log(L( ))-log(L(K)) 13/18 Mustafa Schmidt Updates of DISC Software 13 / 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 14/18 Mustafa Schmidt Updates of DISC Software 14 / 18

Thank you very much for your attention! 15/18 Mustafa Schmidt Updates of DISC Software 15 / 18

Backup Slides 16/18 Mustafa Schmidt Updates of DISC Software 16 / 18

Theoretical Model Angle Definitions : FEL θ c ϕ Particle Radiator Disk tan ϕ tan ϕ ′ = α FEL particle cos α FEL φ rel ϕ ′ 17/18 Mustafa Schmidt Updates of DISC Software 17 / 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 2 θ p − cos 2 θ c cos ϕ = A cos θ c � A cos θ c � + (2) ± B B B with A = sin θ p cos φ rel and B = A 2 + cos 2 θ p 18/18 Mustafa Schmidt Updates of DISC Software 18 / 18