status of analysis software
Skip this Video
Download Presentation
Status of Analysis Software

Loading in 2 Seconds...

play fullscreen
1 / 57

Status of Analysis Software - PowerPoint PPT Presentation

  • Uploaded on

Status of Analysis Software. Kanglin He [email protected] Outline. Physics Tools Project Compose Tools Particle ID Analysis Event Data Model (BParticle Project) Dst Event Data Model Global Link and Event Assembly Analysis Event Data Model Physics Event Data Model Event Tag Project

I am the owner, or an agent authorized to act on behalf of the owner, of the copyrighted work described.
Download Presentation

PowerPoint Slideshow about 'Status of Analysis Software' - baby

An Image/Link below is provided (as is) to download presentation

Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.

- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -
Presentation Transcript
  • Physics Tools Project
    • Compose Tools
    • Particle ID
  • Analysis Event Data Model (BParticle Project)
    • Dst Event Data Model
    • Global Link and Event Assembly
    • Analysis Event Data Model
    • Physics Event Data Model
  • Event Tag Project
  • Physics Performance Check
  • Next Plan
compose tools
Compose Tools
  • Telesis (TKinematicFit)
    • Part of BESII telesis package, by SDU
  • Klams(BKlamsVTX)
    • BESII klams0 program, by B.J. Liu
  • VertexFit:
    • new stuff, by K.L He
    • VertexFit
    • SecondVertexFit
    • KinematicFit


TKinematicFit* kmfit = TKinematicFit::instance();

DstMdcTrack *pipTrk=xxx;

DstMdcTrack *pimTrk = xxx;

DstEmcTrack *g1Trk = xxx;

DstEmcTrack *g2Trk = xxx;


kmfit->AddTrack(1, pipTrk, mpi, 0);

kmfit->AddTrack(2, pimTrk, mpi, 0);

kmfit->AddTrack(3, g1Trk, 0.0, 0);

kmfit->AddTrack(4, g2Trk, 0.0, 0);

kmfit->AddResonance(1, 0.135, 3, 4);

If(kmfit->fit()) {

HepLorentzVector ppip = kmfit->pfit(1);

HepLorentzVector ppim = kmfit->pfit(2);

HepLorentzVector ppi0 = kmfit->pfit(3) + kmfit->pfit(4);

// user code, fill histogram/ntuple


// see Analysis/Physics/RhopiAlg/RhopiAlg-00-00-02 for detail



Liu Beijiang’s talk



for detail

vertexfit and kinematicfit
VertexFit and KinematicFit
  • Formulas(adopt from Paul Avery’s Lessons)
  • WTrackParameter
  • VertexParameter
  • KinematicConstraints
  • VertexConstraints
  • BuildVirtualParticle
  • SecondVertexFit
wtrackparameter and vertexparameter
WTrackParameter and VertexParameter

HepLorentzVector (px, py, pz, E)

HepPoint3D (x, y, z)

HepSymMatrix Ew

HepPoint3D (vx, vy, vz)

HepSymMatrix Ex

addtrack utility
AddTrack Utility

// event data model dependent

void AddTrack(const int number, const double mass, const DstMdcTrack *trk);

void AddTrack(const int number, const double mass, const DstMdcKalTrack *trk);

void AddTrack(const int number, const double mass, const DstEmcTrack *trk);

// event data model independent

void AddTrack(const int number, const WTrackParameter wtrk);

kinematic constraints
Kinematic Constraints
  • Invariant mass constraints
    • AddResonance(number, mass, n1, n2, ...)
  • Total energy constraints
    • AddTotalEnergy(number, energy, n1, n2, …)
  • Total Momentum constraints
    • AddTotalMomentum(number, ptot, n1, n2, …)
  • Total 3-momentum constraints
    • AddThreeMomentum(number, Hep3Vector p3, n1, n2, …)
  • Total 4-momentum constraints
    • AddFourMomentum(number, HepLorentzVector p4, n1, n2, …)
    • AddFourMomentum(number, etot, n1, n2, …)
  • Equal mass constraints
    • Not well designed
vertex constraints
Vertex Constraints
  • Vertex constraint to a fixed position
    • Beam Fit
  • Vertex constraint to an unknown position
    • SecondVertex
    • Find Beam position
  • General Vertex fit
build virtual particle
Build Virtual Particle
  • Virtual Particle from KinematicFit
    • Pi0, eta, …
  • Virtual Particle from VertexFit
    • Ks, Λ, …
  • Can be treated as a normal particle in further kinematic/Vertex fit
  • BuildVirtualParticle(number) in KinematicFit / VertexFit
how to use compose tools
How to use compose tools

TKinematicFit *kmfit=TKinematicFit::instance();

VertexFit *vtxfit=VertexFit::instance();

SecondVertexFit *vtxfit = SecondVertexFit::instance();

KinematicFit *kmfit=KinematicFit::instance();


// Track and constraints

if(kmfit->Fit()) {//extract kinematic fit information;}


// Track and constraints

if(vtxfit->Fit()) {// extract (Second)VertexFit information;}

particle id
Particle ID
  • DedxPID: dE/dx info
    • Ready now, correction was made
  • Tof1PID: inner barrel TOF
    • Ready now, correction was made
  • Tof2PID: outer barrel TOF
    • Ready now, correction was made
  • TofEPID: Endcap TOF
    • Not readdy now, due to reconstruction
  • TofQPID: Q of Barrel TOF
    • Not ready now, need detail study
  • EmcPID: energy deposit in EMC
    • Not ready now, Wang Liangliang is working at
  • ParticleID
    • Combine the above info
  • MuonID
    • Develop by PKU
  • ParticleIDBase
    • For developer only
particleidbase class i common interface
ParticleIDBase class (I)common interface


virtual bool IsPidInfoValid() = 0;

virtual double chi(int n) const = 0;

virtual double prob(int n) const = 0;

virtual int ndof() const = 0;

DstTrkList* PidTrk() const {return m_trk;}


// PID calculation routines(include PID correction) virtual int particleIDCalculation() = 0;

sub pidsys class design
Sub-PidSys class design

class Tof1PID : public ParticleIDBase {

public: // constructor & desconstructor

Tof1PID(DstTrkList *dstTrk, double chi=4);


bool IsPidInfoValid() {return (m_ndof > 0); }

double chi(int n) const {return m_chi[n];}

double prob(int n) const {return m_prob[n];}

int ndof() const {return m_ndof;}


int particleIDCalculation();

double offsetTof1(int n, int cntr, double ptrk, double ztof, double ph);

double sigmaTof1(int n, int cntr, double ptrk, double ztof, double ph);




how to use particleid
How to use ParticleID

// Be careful to avoid memory leakage

ParticleID *pid = 0;

for (int i = 0; i < ncharg; i++) {

DstTrkList *dstTrk = *(dstTrkListCol->begin() + i);

if(pid) delete pid;

pid = new ParticleID(dstTrk);

pid->usePidSys(pid->useDedx() | pid->useTof1() | pid->useTof2() )|;

// pid->usePidSys(pid->useDedx());

// pid->usePidSys(pid->useTof1());

// pid->usePidSys(pid->useTof2());

pid->identify(pid->onlyPion() | pid->onlyKaon());

// pid->identify(pid->onlyPion());

// pid->identify(pid->onlyKaon());


if(!(pid->IsPidInfoValid())) continue;

// user's selections


if(pid) delete pid;

analysis event data model
Analysis Event Data Model
  • Dst Event Data Model
    • From Reconstruction
  • Dst Track List
    • Global Link and Dst Event Assembly
  • Analysis Event Data Model
    • BParticle Project
  • Physics Event Data Model
    • Physics Group
tracklist sub detector level
TrackList(sub-detector level)
  • DstMdcTrack
  • DstMdcKalTrack
  • DstDedx
  • DstTofTrack
  • DstEmcTrack
  • DstMucTrack
  • DstExtTrack
  • Event/DstEvent subdirctory
tracklist track level
TrackList (track level)

int trackID() const {return m_trackID;}

bool IsMdcTrkValid() {return (m_originMdcTrack != 0);}

bool IsMdcKalTrkValid() {return (m_originMdcKalTrack != 0);}

bool IsDedxValid() {return (m_originDedx != 0);}

bool IsTofTrkValid() {return (m_originTofTrack != 0);}

bool IsEmcTrkValid() {return (m_originEmcTrack != 0);}

bool IsMucTrkValid() {return (m_originMucTrack != 0);}

bool IsExtTrkValid() {return (m_originExtTrack != 0);}

DstMdcTrack* MdcTrk() {return m_originMdcTrack;}

DstMdcKalTrack* MdcKalTrk() {return m_originMdcKalTrack;}

DstDedx* DedxTrk() {return m_originDedx;}

DstTofTrack* TofTrk() {return m_originTofTrack;}

DstEmcTrack* EmcTrk() {return m_originEmcTrack;}

DstMucTrack* MucTrk() {return m_originMucTrack;}

DstExtTrack* ExtTrk() {return m_originExtTrack;}

how to access dst event data
How to access Dst Event Data

In analysis source code:

#include “Event/DstTrkList.h”

SmartDataPtr dstCol ( eventSvc(), EventModel::Dst::DstTrackListCol)

// To get Energy deposit in Emc

DstTrkListCol::iterater itTrk = dstCol->begin()

For(; itTrk < dstCol->end(); itTrk++) {

if((*itTrk)->IsEmcTrkValid()) {

DstEmcTrack *emcTrk = (*itTrk)->EmcTrk();

double energy = emcTrk->energy();

// piece of analysis code, fill histogram/ntuple



In analysis job option file:

#include “$DSTEVENTASSEMBLYALGROOT/jobOptions_DstEventAssembly.txt”

bparticle project
BParticle 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
  • BParticle
    • Event Analysis Event Data Model
  • BParticleID
    • Particle identification data
  • BVertex
    • (2nd)Vertex Reconstruction data
  • BStableParticle
    • Photon, Electron, Muon, Pion, Kaon, Proton
  • BComposeParticle
    • Pi0, Eta, K*, phi, Charm Meson, new Particle,……
analysis event data model bparticle
Analysis Event Data ModelBParticle

// track List ID

int m_trackID;


int m_particleID;

// reference point(0, 0, 0)

HepPoint3D m_refpoint;

// WTrackparameters(charge, px, py, pz, e, x, y, z)

WTrackParameter m_wtrk;

SmartRef m_dstTrk; // Dst TrackList

SmartRef m_pid; // ParticleID info

SmartRef m_vertex; // Vertex info

SmartRef m_bcompart; // ComposeParticle

SmartRef m_mother; // mother

SmartRefVector m_daughters; //decay daughters


int m_trackID; // TrackList ID

int m_type; // 1: electron 2: muon 3: pi/K/p

int m_ndof;

double m_chiDedx[5];

double m_chiTof1[5];

double m_chiTof2[5];

double m_chiTofE[5];

double m_chiTofQ[5];

double m_chiEmc[5];

double m_prob[5];

double m_chisq[5];


int m_vertexID;

VertexType m_vertexType;

VertexParameter m_vpar;

double m_chisq;

int m_ndof;

WTrackParameter m_wtrk;

double m_lxyz;

double m_lxyz_err;

SmartRefVector m_outgo;


Pi0, Eta

Ks, Lambda

D, Ds, …

K*, phi, omega, …

Any “new” Particle

bparticle frame work
BParticle Frame work
  • Tracklist duplication
  • VeeVtxReconstruction
    • KShortReconstruction, by B.J. Liu
    • LambdaReconstruction, by K.L. He
  • GoodTrackSelection
    • GoodPhotonSelection, by M.Yang & M.S. Chen
    • GoodElectronSelection
    • GoodMuonSelection
    • GoodHadronSelection, by K.L. He
  • Pi0,Eta Reconstruction by M.S. Chen and M.Yang
  • CharmMesonList, J.Y. Zhang, Y.Z. Sun & S.S.Sun
  • Analysis Cuts are controlled through job option file
bparticle in analysis
BParticle in Analysis

#include “BParticle/BParticle.h”

#include “BParticle/BParticleID.h”

#include “BParticle/BVertex.h”

#include “BParticle/StableParticle.h”

#include “BParticle/BComposedParticle.h”

// get pion Lists

SmartDataPtr pionCol(eventSvc(), EventModel::Analysis::BPionCol);

// get Kaon Lists

SmartDataPtr kaonCol(eventSvc(), EventModel::Analysis::BKaonCol);

// get Photon Lists

SmartDataPtr photonCol(eventSvc(), EventModel::Analysis::BPhotonCol);

//get KShort Lists

SmartDataPtr ksCol(eventSvc(),EventModel::Analysis::BKShortCol);

//get pi0 Lists

SmartDataPtr pi0Col(eventSvc(), EventModel::Analysis::BPi0Col);

//DO not forget to add following line in your job option file

#include “$PARTICLESELECTIONALGROOT/share/jobOPtions_ParticleSelection.txt”

physics event data model
Physics Event Data Model
  • Designed and developed by Physics Group
    • D0, D+, Ds tags
    • Psi’pi pi J/psi
physics performance
Physics Performance
  • Neutrals
    • Spatial resolution of photons (Yang Ming)
    • Photon detection efficiencies (Yang Ming)
    • Pi0 mass resolution (K.L. He, H.B. Li)
  • Particle ID performance
    • TOF (Wang Liangliang
    • dE/dx(Wang Liangliang
    • Muon (Wang Zheng)
  • Charged Tracks
    • Ks, Lambda mass resolution (Liu Beijiang)
    • Phi, K*, Omega in Inclusive Production (K.L. He, H.B. Li)
    • Tracking efficiencies (not ready)
  • Physics Analysis
    • D (J.Y Zhang, Y.Z Sun)
    • J/psi -> gamma X, X->pi0 pi0 (Chen Mingshui)
    • Psip -> pi pi (eta) J/psi (Li Gang)



From Yang M.

J/psi->rho pi

M_gg, raw

J/psi->rho pi

M_gg, 4C fit

From Li H.B


Yang M

Chen M S




Yang M

Chen M S




Yang M

Chen M S


From Chen M S

J/psi->gamma X, X(1270)->pi0 + pi0

Pi0 Reconstruction Efficiencies

Mean_chi, dE/dx

Sigma_chi, dE/dx

From Wang LL etal

M_gg, ~5.1MeV

From Li H B



From Liu B J






From Li H B


From Li H B


D0K pi

From Sun YZ, Zhang J Y

D0K pi pi0

D+K pi pi

From Sun YZ, Zhang J Y

D0 K pi pi pi

next plan
Next Plan
  • Finish EmcPID
  • More Manpower in ParticleID
  • Tracking and PID efficiencies check
  • Applying Kalman filter in KinematicFit and VertexFit
  • Start 1st Vertex Reconstruction
  • Continue the BParticle project
  • Design and develop the EVENT TAG project
  • Some Physics Analysis Tools (BWFIT, ROOFit,… etc)
  • Applying more generator in future