Status of Analysis Software - PowerPoint PPT Presentation

Status of analysis software
1 / 57

  • Uploaded on
  • Presentation posted in: General

Status of Analysis Software. Kanglin He 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

Status of Analysis Software

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

Status of analysis software

Status of Analysis Software

Kanglin He



  • 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

Formulas for kinematic fitting

Formulas for Kinematic Fitting

Formulas for kinematic fitting continue

Formulas for Kinematic Fitting (continue)

Formulas for vertex fitting

Formulas for Vertex Fitting

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

Swimming after vertexfit

Swimming after VertexFit

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<DstTrkListCol> 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<DstTrkList> m_dstTrk; // Dst TrackList

SmartRef<BParticleID> m_pid; // ParticleID info

SmartRef<BVertex> m_vertex; // Vertex info

SmartRef<BComposeParticle> m_bcompart; // ComposeParticle

SmartRef<BParticle> m_mother; // mother

SmartRefVector<BParticle> 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<BParticle> 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<BPionCol> pionCol(eventSvc(), EventModel::Analysis::BPionCol);

// get Kaon Lists

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

// get Photon Lists

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

//get KShort Lists

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

//get pi0 Lists

SmartDataPtr<BPi0Col> 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)

Status of analysis software


From Yang. M

Status of analysis software


From Yang. M

Status of analysis software




From Yang M.

Status of analysis software

Angle between gamma gamma


From Yang M

Status of analysis software

J/psi->rho pi

M_gg, raw

J/psi->rho pi

M_gg, 4C fit

From Li H.B

Status of analysis software


Yang M

Chen M S


Status of analysis software



Yang M

Chen M S

Status of analysis software




Yang M

Chen M S

Status of analysis software


From Chen M S

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

Pi0 Reconstruction Efficiencies

Status of analysis software

Mean_chi, dE/dx

Sigma_chi, dE/dx

From Wang LL etal

Status of analysis software

TOF 1, From Wang LL etal

Status of analysis software

TOF 2, From WangLL etal

Status of analysis software

M_gg, ~5.1MeV

From Li H B

Status of analysis software



From Liu B J



Status of analysis software




From Li H B

Status of analysis software


From Li H B


Status of analysis software

D0K pi

From Sun YZ, Zhang J Y

D0K pi pi0

Status of analysis software

D+K pi pi

From Sun YZ, Zhang J Y

D0 K pi pi pi

Status of analysis software

Psi’pi+ pi- J/Psi, J/Psi  mu+ mu-, raw

From Li G

Status of analysis software

Psi’pi+ pi- J/Psi, J/Psi  mu+ mu-, 4C Fit

From Li. G

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

  • Login