#include <SecondVertexFit.h>

Inheritance diagram for SecondVertexFit:

Public Member Functions
	~SecondVertexFit ()

void	init ()

void	setIterNumber (const int niter=10)

void	setChisqCut (const double chicut=500, const double chiter=1.0e-2)

void	setVpar (const VertexParameter vpar)

void	setPrimaryVertex (const VertexParameter vpar)

bool	Fit ()

double	chisq () const

double	decayLength () const

double	decayLengthError () const

double	ctau () const

double	ctauError () const

HepLorentzVector	p4par () const

WTrackParameter	wpar () const

VertexParameter	vpar () const

HepPoint3D	crossPoint () const

Hep3Vector	crossVector () const

Public Member Functions inherited from TrackPool
	TrackPool ()

	~TrackPool ()

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

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

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

void	AddTrack (const int number, const WTrackParameter wtrk)

void	AddMissTrack (const int number, const double mass)

void	AddMissTrack (const int number, const double mass, const HepLorentzVector p4)

void	AddMissTrack (const int number, const double mass, const RecEmcShower *trk)

void	AddMissTrack (const int number, const RecEmcShower *trk)

void	AddMissTrack (const int number, const HepLorentzVector p4)

void	AddTrackVertex (const int number, const double mass, const RecEmcShower *trk)

std::vector< int >	AddList (int n1)

std::vector< int >	AddList (int n1, int n2)

std::vector< int >	AddList (int n1, int n2, int n3)

std::vector< int >	AddList (int n1, int n2, int n3, int n4)

std::vector< int >	AddList (int n1, int n2, int n3, int n4, int n5)

std::vector< int >	AddList (int n1, int n2, int n3, int n4, int n5, int n6)

std::vector< int >	AddList (int n1, int n2, int n3, int n4, int n5, int n6, int n7)

std::vector< int >	AddList (int n1, int n2, int n3, int n4, int n5, int n6, int n7, int n8)

std::vector< int >	AddList (int n1, int n2, int n3, int n4, int n5, int n6, int n7, int n8, int n9)

std::vector< int >	AddList (int n1, int n2, int n3, int n4, int n5, int n6, int n7, int n8, int n9, int n10)

std::vector< int >	AddList (int n1, int n2, int n3, int n4, int n5, int n6, int n7, int n8, int n9, int n10, int n11)

std::vector< int >	AddList (int n1, int n2, int n3, int n4, int n5, int n6, int n7, int n8, int n9, int n10, int n11, int n12)

std::vector< WTrackParameter >	wTrackOrigin () const

std::vector< WTrackParameter >	wTrackInfit () const

std::vector< int >	wTrackList () const

WTrackParameter	wTrackOrigin (int n) const

WTrackParameter	wTrackInfit (int n) const

int	wTrackList (int n) const

int	numberWTrack () const

std::vector< GammaShape >	GammaShapeValue () const

std::vector< int >	GammaShapeList () const

GammaShape	GammaShapeValue (int n) const

int	GammaShapeList (int n) const

int	numberGammaShape () const

void	setWTrackOrigin (const int n, const WTrackParameter wtrk)

void	setWTrackInfit (const int n, const WTrackParameter wtrk)

void	setWTrackOrigin (const WTrackParameter wtrk)

void	setWTrackInfit (const WTrackParameter wtrk)

void	setWTrackList (const int n)

void	clearWTrackOrigin ()

void	clearWTrackInfit ()

void	clearWTrackList ()

void	clearone ()

void	cleartwo ()

int	numberone () const

int	numbertwo () const

vector< int >	mappositionA () const

vector< int >	mappositionB () const

vector< int >	mapkinematic () const

void	clearMapkinematic ()

void	clearMappositionA ()

void	clearMappositionB ()

void	setMapkinematic (const int n)

void	setMappositionA (const int n)

void	setMappositionB (const int n)

void	setGammaShape (const int n, const GammaShape gammashape)

void	setGammaShape (const GammaShape gammashape)

void	setGammaShapeList (const int n)

void	clearGammaShape ()

void	clearGammaShapeList ()

void	setBeamPosition (const HepPoint3D BeamPosition)

void	setVBeamPosition (const HepSymMatrix VBeamPosition)

HepPoint3D	getBeamPosition () const

HepSymMatrix	getVBeamPosition () const

Static Public Member Functions
static SecondVertexFit *	instance ()

Detailed Description

Definition at line 18 of file SecondVertexFit.h.

Constructor & Destructor Documentation

◆ ~SecondVertexFit()

SecondVertexFit::~SecondVertexFit ( )

Definition at line 64 of file SecondVertexFit.cxx.

{
//  if(m_pointer) delete m_pointer;
}

Member Function Documentation

◆ chisq()

double SecondVertexFit::chisq ( ) const

inline

Definition at line 41 of file SecondVertexFit.h.

41{return m_chisq;}

Referenced by DQAKsInfo::calculate(), inclks::execute(), incllambda::execute(), DQAKsKpi::execute(), Fit(), K0kk::MTotal(), K0kpi::MTotal(), K0pi0::MTotal(), K0pipi::MTotal(), and K0pipipi0::MTotal().

◆ crossPoint()

HepPoint3D SecondVertexFit::crossPoint ( ) const

inline

Definition at line 75 of file SecondVertexFit.h.

{
    return HepPoint3D(m_crxyz[0], m_crxyz[1], m_crxyz[2]);
}

◆ crossVector()

Hep3Vector SecondVertexFit::crossVector ( ) const

inline

Definition at line 80 of file SecondVertexFit.h.

{
    return Hep3Vector(m_crxyz[0] - m_vpar_primary.Vx()[0], 
                      m_crxyz[1] - m_vpar_primary.Vx()[1],
                      m_crxyz[2] - m_vpar_primary.Vx()[2]);
}

◆ ctau()

double SecondVertexFit::ctau ( ) const

inline

Definition at line 44 of file SecondVertexFit.h.

44{return m_ctau;}

Referenced by inclks::execute(), incllambda::execute(), and DQAKsKpi::execute().

◆ ctauError()

double SecondVertexFit::ctauError ( ) const

inline

Definition at line 45 of file SecondVertexFit.h.

45{return m_ctau_error;}

◆ decayLength()

double SecondVertexFit::decayLength ( ) const

inline

Definition at line 42 of file SecondVertexFit.h.

42{return m_lxyz;}

Referenced by DQAKsInfo::calculate(), inclks::execute(), incllambda::execute(), DQAKsKpi::execute(), K0kk::MTotal(), K0kpi::MTotal(), K0pi0::MTotal(), K0pipi::MTotal(), and K0pipipi0::MTotal().

◆ decayLengthError()

double SecondVertexFit::decayLengthError ( ) const

inline

Definition at line 43 of file SecondVertexFit.h.

43{return m_lxyz_error;}

Referenced by DQAKsInfo::calculate(), inclks::execute(), incllambda::execute(), and DQAKsKpi::execute().

◆ Fit()

bool SecondVertexFit::Fit ( )

Definition at line 86 of file SecondVertexFit.cxx.

{
//  const double alpha = -0.00299792458;
//  const double field = 1.0;
//  double field_p = VertexFitBField::instance()->getBFieldZ(m_vpar_primary.vx());
//  double field_s = VertexFitBField::instance()->getBFieldZ(m_vpar_secondary.vx());
//  double field = (field_p + field_s) * 0.5 ;
//  std::cout << "field in SeconVertexFit= " << field << std::endl;
 
    bool okfit = false;
 
    HepVector aOrigin(10, 0);
    HepVector aInfit(10, 0);
    HepSymMatrix VaOrigin(10, 0);
    HepSymMatrix VaInfit(10, 0);
    aOrigin.sub(1, wTrackOrigin(0).w());
    aOrigin.sub(8, m_vpar_primary.Vx());
    VaOrigin.sub(1, wTrackOrigin(0).Ew());
    VaOrigin.sub(8, m_vpar_primary.Evx());
    HepVector ctOrigin(1, 0);
    HepVector ctInfit(1, 0);
    HepSymMatrix Vct(1, 0);
    aInfit  = aOrigin;
    ctInfit = ctOrigin;
 
    std::vector<double> chisq;
    chisq.clear();
    double chi2 = 999;
    for(int it = 0; it < m_niter; it++)
    {
        HepMatrix D(3, 10, 0);
        HepLorentzVector p4par = HepLorentzVector(aInfit[0], aInfit[1], aInfit[2], aInfit[3]);
        HepMatrix E(3,1,0);
        HepVector d(3, 0);
        if (wTrackOrigin(0).charge() == 0) 
        {
            D[0][4] = -1.0;
            D[0][7] = 1.0;
            D[1][5] = -1.0;
            D[1][8] = 1.0;
            D[2][6] = -1.0;
            D[2][9] = 1.0;
 
            E[0][0] = p4par.px()/p4par.m();
            E[1][0] = p4par.py()/p4par.m();
            E[2][0] = p4par.pz()/p4par.m();
 
            d[0] = aInfit[7]-aInfit[4]+ctInfit[0]*p4par.px()/p4par.m();
            d[1] = aInfit[8]-aInfit[5]+ctInfit[0]*p4par.py()/p4par.m();
            d[2] = aInfit[9]-aInfit[6]+ctInfit[0]*p4par.pz()/p4par.m();
        }
        else 
        {
            double afield = VertexFitBField::instance()->getCBz(m_vpar_primary.Vx(), m_vpar_secondary.Vx());
            double a = afield * wTrackOrigin(0).charge();
            D[0][4] = -1.0;
            D[0][7] = 1.0;
            D[1][5] = -1.0;
            D[1][8] = 1.0;
            D[2][6] = -1.0;
            D[2][9] = 1.0;
 
            E[0][0] = p4par.px() / p4par.m() * cos(a * ctInfit[0] / p4par.m()) + p4par.py() / p4par.m() * sin(a * ctInfit[0] / p4par.m());
            E[1][0] = p4par.py() / p4par.m() * cos(a * ctInfit[0] / p4par.m()) - p4par.px() / p4par.m() * sin(a * ctInfit[0] / p4par.m());
            E[2][0] = p4par.pz() / p4par.m();
 
            d[0] = aInfit[7] - aInfit[4]+p4par.px()/a * sin(a * ctInfit[0]/p4par.m()) + p4par.py()/a*(1-cos(a*ctInfit[0]/p4par.m()));
            d[1] = aInfit[8] - aInfit[5]+p4par.py()/a * sin(a * ctInfit[0]/p4par.m()) - p4par.px()/a*(1-cos(a*ctInfit[0]/p4par.m()));
            d[2] = aInfit[9] - aInfit[6]+ctInfit[0]*p4par.pz()/p4par.m();
        }
        
        HepSymMatrix VD(3, 0);
        HepVector dela0(10, 0);
        HepVector lambda0(3, 0);
        HepVector delct(1, 0);
        HepVector lambda(3, 0);
        int ifail;
 
        VD      = (VaOrigin.similarity(D)).inverse(ifail);
        dela0   = aOrigin - aInfit;
        lambda0 = VD*(D*dela0 + d);
        Vct     = (VD.similarity(E.T())).inverse(ifail);
        delct   = -(Vct * E.T()) * lambda0;
        ctInfit = ctInfit + delct;
        lambda  = lambda0 + (VD * E) * delct;
        aInfit  = aOrigin - (VaOrigin * D.T()) * lambda;
        chi2    = dot(lambda, D*dela0 + d);
        VaInfit = VaOrigin - (VD.similarity(D.T())).similarity(VaOrigin);
        VaInfit = VaInfit + (((Vct.similarity(E)).similarity(VD)).similarity(D.T())).similarity(VaOrigin);
 
        chisq.push_back(chi2);
 
        if(it > 0) 
        {
            double delchi = chisq[it] - chisq[it-1];
            if (fabs(delchi) < m_chiter) 
            {
//              std::cout <<"Fit(): pass all " << it <<" , " << delchi << std::endl;
                break;
            }
        }
    }
    if (chi2 < 0 || chi2 > m_chicut)
        return okfit;
 
    HepLorentzVector p4par = HepLorentzVector(aInfit[0], aInfit[1], aInfit[2], aInfit[3]);
    m_ctau       = ctInfit[0];
    m_ctau_error = sqrt(Vct[0][0]);
    m_lxyz       = ctInfit[0] * p4par.rho() / p4par.m();
    m_lxyz_error = sqrt(Vct[0][0]) * p4par.rho() / p4par.m();
    m_chisq      = chi2;
    m_p4par      = p4par;
    for(int i = 0; i < 3; i++)
        m_crxyz[i] = aInfit[4+i];
    HepVector w(7, 0);
    HepSymMatrix Ew(7, 0);
    for(int i = 0; i < 7; i++) 
    {
        w[i] = aInfit[i];
        for(int j = 0; j < 7; j++) 
        {
            Ew[i][j] = VaInfit[i][j];
        }
    }
    m_wtrk.setW(w);
    m_wtrk.setEw(Ew);
    m_wtrk.setCharge(wTrackOrigin(0).charge());
    okfit = true;
    return okfit;
}

Referenced by DQAKsInfo::calculate(), inclks::execute(), incllambda::execute(), DQAKsKpi::execute(), JsiLL::execute(), K0kk::MTotal(), K0kpi::MTotal(), K0pi0::MTotal(), K0pipi::MTotal(), and K0pipipi0::MTotal().

◆ init()

void SecondVertexFit::init ( )

Definition at line 69 of file SecondVertexFit.cxx.

{
    clearWTrackOrigin();
    clearWTrackInfit();
    clearWTrackList();
    m_vpar_secondary = VertexParameter();
    m_lxyz       = 0;
    m_lxyz_error = 0;
    m_p4par      = HepLorentzVector(0, 0, 0, 0);
    m_crxyz      = HepVector(3, 0);
    m_chisq      = 9999;
    m_wtrk       = WTrackParameter();
    m_niter      = 10;
    m_chicut     = 500;
    m_chiter     = 1.0e-2;
}

Referenced by DQAKsInfo::calculate(), inclks::execute(), incllambda::execute(), DQAKsKpi::execute(), JsiLL::execute(), DQAKsInfo::init(), K0kk::MTotal(), K0kpi::MTotal(), K0pi0::MTotal(), K0pipi::MTotal(), and K0pipipi0::MTotal().

◆ instance()

SecondVertexFit * SecondVertexFit::instance ( )

static

Definition at line 45 of file SecondVertexFit.cxx.

{
    if (m_pointer) 
        return m_pointer;
    m_pointer = new SecondVertexFit();
    return m_pointer;
}

Referenced by DQADtag::DQADtag(), DQAKsInfo::DQAKsInfo(), inclks::execute(), incllambda::execute(), DQAKsKpi::execute(), JsiLL::execute(), K0kk::MTotal(), K0kpi::MTotal(), K0pi0::MTotal(), K0pipi::MTotal(), and K0pipipi0::MTotal().

◆ p4par()

HepLorentzVector SecondVertexFit::p4par ( ) const

inline

Definition at line 47 of file SecondVertexFit.h.

47{return m_p4par;}

Referenced by DQAKsInfo::calculate(), inclks::execute(), incllambda::execute(), Fit(), K0kk::MTotal(), K0kpi::MTotal(), K0pi0::MTotal(), K0pipi::MTotal(), and K0pipipi0::MTotal().

◆ setChisqCut()

void SecondVertexFit::setChisqCut	(	const double	chicut = `500`,
		const double	chiter = `1.0e-2`
	)

inline

Definition at line 31 of file SecondVertexFit.h.

31{m_chicut = chicut;m_chiter=chiter;}

Referenced by DQAKsInfo::calculate().

◆ setIterNumber()

void SecondVertexFit::setIterNumber ( const int niter = 10 )

inline

Definition at line 30 of file SecondVertexFit.h.

30{m_niter = niter;}

◆ setPrimaryVertex()

void SecondVertexFit::setPrimaryVertex ( const VertexParameter vpar )

inline

Definition at line 35 of file SecondVertexFit.h.

35{m_vpar_primary.setVx(vpar.Vx());m_vpar_primary.setEvx(vpar.Evx());}

SecondVertexFit::vpar

VertexParameter vpar() const

Definition: SecondVertexFit.h:49

VertexParameter::setEvx

void setEvx(const HepSymMatrix &eVx)

Definition: VertexParameter.h:29

VertexParameter::setVx

void setVx(const HepPoint3D &vx)

Definition: VertexParameter.h:54

Referenced by DQAKsInfo::calculate(), DQAKsKpi::execute(), K0kk::MTotal(), K0kpi::MTotal(), K0pi0::MTotal(), K0pipi::MTotal(), and K0pipipi0::MTotal().

◆ setVpar()

void SecondVertexFit::setVpar ( const VertexParameter vpar )

inline

Definition at line 33 of file SecondVertexFit.h.

33{m_vpar_secondary = vpar;}

Referenced by DQAKsInfo::calculate(), inclks::execute(), incllambda::execute(), DQAKsKpi::execute(), JsiLL::execute(), K0kk::MTotal(), K0kpi::MTotal(), K0pi0::MTotal(), K0pipi::MTotal(), and K0pipipi0::MTotal().

◆ vpar()

VertexParameter SecondVertexFit::vpar ( ) const

inline

Definition at line 49 of file SecondVertexFit.h.

49{return m_vpar_secondary;}

Referenced by setPrimaryVertex(), and setVpar().

◆ wpar()

WTrackParameter SecondVertexFit::wpar ( ) const

inline

Definition at line 48 of file SecondVertexFit.h.

48{return m_wtrk;}

Referenced by DQAKsKpi::execute(), JsiLL::execute(), and K0pipi::MTotal().

The documentation for this class was generated from the following files:

Public Member Functions

Static Public Member Functions

Detailed Description

Constructor & Destructor Documentation

◆ ~SecondVertexFit()

Member Function Documentation

◆ chisq()

◆ crossPoint()

◆ crossVector()

◆ ctau()

◆ ctauError()

◆ decayLength()

◆ decayLengthError()

◆ Fit()

◆ init()

◆ instance()

◆ p4par()

◆ setChisqCut()

◆ setIterNumber()

◆ setPrimaryVertex()

◆ setVpar()

◆ vpar()

◆ wpar()