#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)

void	setMagCorrFactor (const double factor=1.000)

bool	Fit ()

double	chisq () const

double	decayLength () const

double	decayLengthError () const

double	ctau () const

double	ctauError () const

double	factor () 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 20 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 45 of file SecondVertexFit.h.

45{return m_chisq;}

Referenced by DQAKsInfo::calculate(), inclks::execute(), incllambda::execute(), DQAKsKpi::execute(), DQAKsKpiDEDX::execute(), Signal::execute(), Fit(), K0kk::MTotal(), K0kpi::MTotal(), K0pi0::MTotal(), K0pipi::MTotal(), K0pipipi0::MTotal(), LocalKsSelector::operator()(), LocalLambdaSelector::operator()(), utility::SecondaryVFit(), and utility::SecondaryVFit_Lambda().

◆ crossPoint()

HepPoint3D SecondVertexFit::crossPoint ( ) const

inline

Definition at line 81 of file SecondVertexFit.h.

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

◆ crossVector()

Hep3Vector SecondVertexFit::crossVector ( ) const

inline

Definition at line 86 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 48 of file SecondVertexFit.h.

48{return m_ctau;}

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

◆ ctauError()

double SecondVertexFit::ctauError ( ) const

inline

Definition at line 49 of file SecondVertexFit.h.

49{return m_ctau_error;}

◆ decayLength()

double SecondVertexFit::decayLength ( ) const

inline

Definition at line 46 of file SecondVertexFit.h.

46{return m_lxyz;}

Referenced by DQAKsInfo::calculate(), inclks::execute(), incllambda::execute(), DQAKsKpi::execute(), DQAKsKpiDEDX::execute(), Signal::execute(), K0kk::MTotal(), K0kpi::MTotal(), K0pi0::MTotal(), K0pipi::MTotal(), K0pipipi0::MTotal(), LocalKsSelector::operator()(), LocalLambdaSelector::operator()(), utility::SecondaryVFit(), and utility::SecondaryVFit_Lambda().

◆ decayLengthError()

double SecondVertexFit::decayLengthError ( ) const

inline

Definition at line 47 of file SecondVertexFit.h.

47{return m_lxyz_error;}

Referenced by DQAKsInfo::calculate(), inclks::execute(), incllambda::execute(), DQAKsKpi::execute(), DQAKsKpiDEDX::execute(), Signal::execute(), LocalKsSelector::operator()(), LocalLambdaSelector::operator()(), utility::SecondaryVFit(), and utility::SecondaryVFit_Lambda().

◆ factor()

double SecondVertexFit::factor ( ) const

inline

Definition at line 50 of file SecondVertexFit.h.

50{return m_factor;}

Referenced by setMagCorrFactor().

◆ Fit()

bool SecondVertexFit::Fit ( )

Definition at line 87 of file SecondVertexFit.cxx.

{
    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 afield = m_factor * 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) 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(), DQAKsKpiDEDX::execute(), JsiLL::execute(), Signal::execute(), K0kk::MTotal(), K0kpi::MTotal(), K0pi0::MTotal(), K0pipi::MTotal(), K0pipipi0::MTotal(), LocalKsSelector::operator()(), LocalLambdaSelector::operator()(), utility::SecondaryVFit(), and utility::SecondaryVFit_Lambda().

◆ 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;
    m_factor     = 1.000;
}

Referenced by DQAKsInfo::calculate(), inclks::execute(), incllambda::execute(), DQAKsKpi::execute(), DQAKsKpiDEDX::execute(), JsiLL::execute(), Signal::execute(), DQAKsInfo::init(), K0kk::MTotal(), K0kpi::MTotal(), K0pi0::MTotal(), K0pipi::MTotal(), K0pipipi0::MTotal(), LocalKsSelector::operator()(), LocalLambdaSelector::operator()(), utility::SecondaryVFit(), and utility::SecondaryVFit_Lambda().

◆ 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(), DQAKsKpiDEDX::execute(), JsiLL::execute(), Signal::execute(), K0kk::MTotal(), K0kpi::MTotal(), K0pi0::MTotal(), K0pipi::MTotal(), K0pipipi0::MTotal(), LocalKsSelector::operator()(), LocalLambdaSelector::operator()(), utility::SecondaryVFit(), and utility::SecondaryVFit_Lambda().

◆ p4par()

HepLorentzVector SecondVertexFit::p4par ( ) const

inline

Definition at line 52 of file SecondVertexFit.h.

52{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 33 of file SecondVertexFit.h.

33{m_chicut = chicut;m_chiter=chiter;}

Referenced by DQAKsInfo::calculate().

◆ setIterNumber()

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

inline

Definition at line 32 of file SecondVertexFit.h.

32{m_niter = niter;}

◆ setMagCorrFactor()

void SecondVertexFit::setMagCorrFactor ( const double factor = 1.000 )

inline

Definition at line 39 of file SecondVertexFit.h.

39{m_factor = factor;}

SecondVertexFit::factor

double factor() const

Definition: SecondVertexFit.h:50

◆ setPrimaryVertex()

void SecondVertexFit::setPrimaryVertex ( const VertexParameter vpar )

inline

Definition at line 37 of file SecondVertexFit.h.

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

SecondVertexFit::vpar

VertexParameter vpar() const

Definition: SecondVertexFit.h:54

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(), DQAKsKpiDEDX::execute(), Signal::execute(), K0kk::MTotal(), K0kpi::MTotal(), K0pi0::MTotal(), K0pipi::MTotal(), K0pipipi0::MTotal(), LocalKsSelector::operator()(), LocalLambdaSelector::operator()(), utility::SecondaryVFit(), and utility::SecondaryVFit_Lambda().

◆ setVpar()

void SecondVertexFit::setVpar ( const VertexParameter vpar )

inline

Definition at line 35 of file SecondVertexFit.h.

35{m_vpar_secondary = vpar;}

Referenced by DQAKsInfo::calculate(), inclks::execute(), incllambda::execute(), DQAKsKpi::execute(), DQAKsKpiDEDX::execute(), JsiLL::execute(), Signal::execute(), K0kk::MTotal(), K0kpi::MTotal(), K0pi0::MTotal(), K0pipi::MTotal(), K0pipipi0::MTotal(), LocalKsSelector::operator()(), LocalLambdaSelector::operator()(), utility::SecondaryVFit(), and utility::SecondaryVFit_Lambda().

◆ vpar()

VertexParameter SecondVertexFit::vpar ( ) const

inline

Definition at line 54 of file SecondVertexFit.h.

54{return m_vpar_secondary;}

Referenced by setPrimaryVertex(), and setVpar().

◆ wpar()

WTrackParameter SecondVertexFit::wpar ( ) const

inline

Definition at line 53 of file SecondVertexFit.h.

53{return m_wtrk;}

Referenced by DQAKsKpi::execute(), DQAKsKpiDEDX::execute(), JsiLL::execute(), Signal::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()

◆ factor()

◆ Fit()

◆ init()

◆ instance()

◆ p4par()

◆ setChisqCut()

◆ setIterNumber()

◆ setMagCorrFactor()

◆ setPrimaryVertex()

◆ setVpar()

◆ vpar()

◆ wpar()