Status of Analysis SoftwareStatus of Analysis Software

Kanglin HeKanglin He

For Analysis Software GroupFor Analysis Software Group

OutlineOutline Physics Tools Project

Particle ID Composition Tools

Analysis Event Data Model Global Link and Event Assembly Analysis Event Data Model BParticle Project

Event/Physics Tags Project Event Classification Physics Tags Event Collection, Filter

Physics Performance Check Summary

Status of Particle IDStatus of Particle ID DedxPID

Ready now TofPID

Barrel inner TOF PID Ready now

Barrel outer TOF PID Ready now

MuonID Ready now, by PKU

EmcPID Not ready

Endcap TOF PID Not ready

TOFQ PID Not ready

See Dr. Ji X.B ‘s talk for performance check

Composition ToolsComposition Tools Telesis

BESII Kinematic Fit package C++ implemented by SDU

BKlamsVTX BESII Ks, Λ reconstruction package C++ implemented by Dr. B.J. Liu

VertexFit VertexFit SecondVertexFit KinematicFit

Status of TelesisStatus of Telesis

Track presentation (φ,κ,λ) for charged track (φ,λ, sqrt(E)) for neutral track

Constraints Resonance Momentum-energy conservation Miss track

Worked smoothly in BOSS

Status of BKlamsVTXStatus of BKlamsVTX

Get the point (s) of intersection of two helices in transverse plane

Check if the 2nd vertex candidate is located in reasonable region

Worked smoothly in BOSS

See B. J. Liu’s talk for detail

Status of Vertex FitStatus of Vertex Fit Formulas (adopt from Paul Avery’s Lessons) Include: VertexFit, KinematicFit and 3-D

SecondVertexFit Support Composed Particle

Reconstructed π0, η Reconstructed Ks, Λ

Particle and Vertex presentation

W

zyx

V

z)y,x,E,,p,p,(pW

charge

Vx

z)y,(x,

Kinematic ConstraintsKinematic Constraints

Invariant mass constraints Total energy constraints Total Momentum constraints Total 3-momentum constraints Total 4-momentum constraints Equal mass constraints

Vertex ConstraintsVertex Constraints Vertex constraint to a fixed position Vertex constraint to an unknown position General Vertex fit Swimming Track parameters and Error matrix to vertex

position

2

1

22

hp

)hx)(h(pΔxpasin

a

hphΔx0

)hx(Δx2

ahΔx)(p0

SecondVertexFitSecondVertexFit Find decay vertex position of Ks/Λ through

VertexFit Performing a 2-C fit by applying the alignment

constraints Calculating the decay length and its error

0cτm

pzz

0cτm

pyy

0cτm

pxx

zdp

ydp

xdp

Xp: production position (beam position)Xd: decay position

Summary of Analysis ToolsSummary of Analysis Tools

Particle ID package Hadron identification by combined TOF and

dE/dx information Muon identification by MuonID Electron identification is not completed

All Composition Tools packaged are worked smoothly and stably in BOSS

Analysis Event Data ModelAnalysis Event Data Model

Dst Event Data Model From Reconstruction DstTrackList

Global Link and Dst Event Assembly

Analysis Event Data Model BParticle Project

Physics Event Data Model Physics Group

Analysis Event Data FlowAnalysis Event Data Flow

MdcTrackCol

MdcKalTrackCol

DedxCol

TofTrackCol

EmcTrackCol

MucTrackCol

ExtTrackCol

Track #1MdcTrack

MdcKalTrackDedx

TofTrackEmcTrackMucTrackExtTrackTrack #2

MdcTrackMdcKalTrack

DedxTofTrackEmcTrackMucTrackExtTrackTrack #...MdcTrack

MdcKalTrackDedx

TofTrackEmcTrackMucTrackExtTrackTrack #NMdcTrack

MdcKalTrackDedx

TofTrackEmcTrackMucTrackExtTrack

BParticle

BTrack

BParticleID

BVertex

BPhotonCol

BElectronCol

BMuonCol

BPionCol

BKaonCol

BProtonCol

BPi0Col

...BKShortCol

BLambaCol

Bhabha Event

Dimu Event

GammaGammaEvent

...

Hadron Event

Tau PairEvent

D0 Tags

D+ Tags

DsTags

...

J/Psi Tags

Other PhysicsTags

Dst EventData Model

Analysis EventData Model

Physics/EventTag

Global LinkGlobal Link MDC reconstruction

dE/dx, Kalman fit track Match to TOF and Muon

counter Track Applying Track Extrapolatio

information, done in Reconstruction

Match to EMC Position difference of extTrk

and emc Track, treat as Neutral tracks if not matched

See Wang L.L ‘s talk for detail

Dst Track List CollectionDst Track List Collection trackID() MdcTrk() DstMdcTrack MdcKalTrk() DstMdcKalTrack Dedx() DstDedx TofTrk() DstTofTrack EmcTrk() DstEmcTrack MucTrk() DstMucTrack ExtTrk() DstExtTrack

Access sub-detector info through a set of pointers

B(es)Particle ProjectB(es)Particle Project

Write Particle (Charged tracks, neutrals, Vertice and composed particle) information into TDS. Analysis Algorithm can access these info.

May develop to a common standard platform, make Analysis job easily and friendly

Allow Physics groups to design Physics Data model easily

ContentsContents BParticle Collection

Event Analysis Event Data Model BParticleID Collection

Particle identification data BTrack Collection

Charged and Neutral Track infomation data BVertex Collection

(2nd)Vertex Reconstruction data BStableParticle Collection

Photon, Electron, Muon, Pion, Kaon, Proton BComposedParticle Collection

π0,η,Ks,Λ etc

B(es)Particle ProjectB(es)Particle Project

BParticleCol

BVertexBTrackColBParticleID

BComposedParticle

BKShortCol

BLambdaCol

BPi0Col

BEtaCol

BPhotonCol

BElectronCol

BMuonCol

BPionCol

BKaonCol

BProtonCol

Status of BTrackStatus of BTrack Photon selection algorithm

Energy Threshold Angle between photon and Charged Track Angle between photons

Charged Track Selection algorithm Impact parameters (Rxy and Z0),Pxy ,Polar angle Particle ID information

Selection criteria controlled by job option file Good Tracks are put into Photon, electron, Muon, Pion,

Kaon, Proton Collection

See X.B Ji and M.S Chen’s talk for detail

Status of BVertexStatus of BVertex Ks and Λ decay vertex are implemented

BKlamsVTX SecondVertexFit

Gamma conversion vertex will come out soon 1st Vertex reconstruction will be done in Physics

analysis level Well reconstructed Ks and Λ are put into BParticle

Collection (selection criteria controlled by job Option file)

See B.J Liu’s talk for detail

Status of Pi0/Eta ReconstructionStatus of Pi0/Eta Reconstruction

π0/ηγγ Raw mass region Chisq of mass constraint fit Decay angle cut

Reconstructed π0/η are put into Pi0/Eta Collection

Selection criteria controlled by job Option file

See M.S Chen’s talk for detail

Summary of Analysis Summary of Analysis Event Data ModelEvent Data Model

Global Link and Dst Event Assembly Integrated into BOSS beta Release First analysis J/Ψρπ come out in the middle of October 2005

After BOSS beta release, a prototype of BESIII analysis event data model was designed BParticle Collection BTrack/BParticleID Collection BVertex Collection BStableParticle/BComposedParticle Collection

Analysis Event Data Model was implemented by BParticle Project Charged/Neutral Track Selection Algorithm Vertex Reconstruction/Selection Algorithm π0/η Reconstruction/Selection Algorithm

BParticle Project will be developed to Standard BESIII Analysis Framework

Status of Event/Physics Tags Status of Event/Physics Tags ProjectProject

Event Classification (will be started) Bhabha/Dimu/Hadron/etc. X.H Mo is working on

Physics Tags (started) Physics Tag Event Data Model

Designed by Physics Group Charm Meson Tags (D0, D+, Ds)

J.Y Zhang, Y.Z Sun, S.S Sun are working on J/Ψ tags (Ψ’ππJ/Ψ)

X.H Mo, Li Gang are working on Other Physics Tags

Event Collection/Filter (not started) Event Collection Filter: sub-sets of data containing selected tags

Physics Performance CheckPhysics Performance Check Charged Track

Hadron Performance (X.B Ji’s talk) MuonID Performance (PKU)

Photon(π0/η ) performance M.S Chen’s talk

Decay Vertex Reconstruction B.J Liu’s talk

J/Ψ,Ψ’, D/Ds Analysis X.H Mo, X.B Ji, J.Y Zhang ‘s talks

Energy threshold at 40MeV

GeV

Eff

icie

ncie

s(10

0%)

Photon detection Efficiencies as a function of E

20MeV30MeV

40MeV50MeV

60MeV70MeV

From Li H.B

ρπJ/ψ

M

M

ρπJ/ψ

No 4C fit

After 4C fit

γγπ0

γγπ0

From Li H B

M

γγπ0

50 K

Inc

lusi

ve Ψ

” M

C d

ata

~3MeV

~1.2MeV

ππKS

pπ

50 K

Inc

lusi

ve Ψ

” M

C d

ata

KKφ

0πππω

50 K

Inc

lusi

ve Ψ

” M

C d

ata

50K

Inc

lusi

ve Ψ

” M

C d

ata

KπK*0

πKK S*

From Y.Z Sun

πKD0 πππKD0

π

00 ππKD ππKD S

0

From J. Y Zhang

0S ππKD

πKD S

πππKD S

ππKD

J/ψππψ

recM

From Li Gang

50K inclusive Ψ’ MC data

Inclusive photon Energy

0cχ

1cχ2

cχ

γJ/ψχ1,2c

E

μμJ/ψ

SummarySummary Analysis Software in Good but not perfect status,

physics analysis can be started Need More Manpower in ParticleID package A prototype of BESIII Analysis Event Data Model has

been designed and developed Event/Physics Tags project will be developed in 2006 Many(but not enough) Physics Performance checks

have been done

Thank You

Welcome to join us