SecondVertexFit Class Reference

#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
	~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
	~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
	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
	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 ()
static SecondVertexFit *	instance ()
static SecondVertexFit *	instance ()

Detailed Description

Definition at line 20 of file Analysis/VertexFit/include/VertexFit/SecondVertexFit.h.

Constructor & Destructor Documentation

~SecondVertexFit() [1/3]

SecondVertexFit::~SecondVertexFit

(

)

Definition at line 61 of file SecondVertexFit.cxx.

                                  {
  // if(m_pointer) delete m_pointer;
}

~SecondVertexFit() [2/3]

SecondVertexFit::~SecondVertexFit

(

)

~SecondVertexFit() [3/3]

SecondVertexFit::~SecondVertexFit

(

)

Member Function Documentation

chisq() [1/3]

double SecondVertexFit::chisq ( ) const

inline

Definition at line 49 of file Analysis/VertexFit/include/VertexFit/SecondVertexFit.h.

{ return m_chisq; }

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

chisq() [2/3]

double SecondVertexFit::chisq ( ) const

inline

Definition at line 49 of file InstallArea/x86_64-el9-gcc13-dbg/include/VertexFit/SecondVertexFit.h.

{ return m_chisq; }

chisq() [3/3]

double SecondVertexFit::chisq ( ) const

inline

Definition at line 49 of file InstallArea/x86_64-el9-gcc13-opt/include/VertexFit/SecondVertexFit.h.

{ return m_chisq; }

crossPoint() [1/3]

HepPoint3D SecondVertexFit::crossPoint ( ) const

inline

Definition at line 84 of file Analysis/VertexFit/include/VertexFit/SecondVertexFit.h.

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

crossPoint() [2/3]

HepPoint3D SecondVertexFit::crossPoint ( ) const

inline

crossPoint() [3/3]

HepPoint3D SecondVertexFit::crossPoint ( ) const

inline

crossVector() [1/3]

Hep3Vector SecondVertexFit::crossVector ( ) const

inline

Definition at line 88 of file Analysis/VertexFit/include/VertexFit/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] );
}

crossVector() [2/3]

Hep3Vector SecondVertexFit::crossVector ( ) const

inline

crossVector() [3/3]

Hep3Vector SecondVertexFit::crossVector ( ) const

inline

ctau() [1/3]

double SecondVertexFit::ctau ( ) const

inline

Definition at line 52 of file Analysis/VertexFit/include/VertexFit/SecondVertexFit.h.

{ return m_ctau; }

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

ctau() [2/3]

double SecondVertexFit::ctau ( ) const

inline

Definition at line 52 of file InstallArea/x86_64-el9-gcc13-dbg/include/VertexFit/SecondVertexFit.h.

{ return m_ctau; }

ctau() [3/3]

double SecondVertexFit::ctau ( ) const

inline

Definition at line 52 of file InstallArea/x86_64-el9-gcc13-opt/include/VertexFit/SecondVertexFit.h.

{ return m_ctau; }

ctauError() [1/3]

double SecondVertexFit::ctauError ( ) const

inline

Definition at line 53 of file Analysis/VertexFit/include/VertexFit/SecondVertexFit.h.

{ return m_ctau_error; }

ctauError() [2/3]

double SecondVertexFit::ctauError ( ) const

inline

Definition at line 53 of file InstallArea/x86_64-el9-gcc13-dbg/include/VertexFit/SecondVertexFit.h.

{ return m_ctau_error; }

ctauError() [3/3]

double SecondVertexFit::ctauError ( ) const

inline

Definition at line 53 of file InstallArea/x86_64-el9-gcc13-opt/include/VertexFit/SecondVertexFit.h.

{ return m_ctau_error; }

decayLength() [1/3]

double SecondVertexFit::decayLength ( ) const

inline

Definition at line 50 of file Analysis/VertexFit/include/VertexFit/SecondVertexFit.h.

{ return m_lxyz; }

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

decayLength() [2/3]

double SecondVertexFit::decayLength ( ) const

inline

Definition at line 50 of file InstallArea/x86_64-el9-gcc13-dbg/include/VertexFit/SecondVertexFit.h.

{ return m_lxyz; }

decayLength() [3/3]

double SecondVertexFit::decayLength ( ) const

inline

Definition at line 50 of file InstallArea/x86_64-el9-gcc13-opt/include/VertexFit/SecondVertexFit.h.

{ return m_lxyz; }

decayLengthError() [1/3]

double SecondVertexFit::decayLengthError ( ) const

inline

Definition at line 51 of file Analysis/VertexFit/include/VertexFit/SecondVertexFit.h.

{ return m_lxyz_error; }

Referenced by DQAKsKpi::execute(), DQAKsKpiDEDX::execute(), inclks::execute(), incllambda::execute(), KsKpi::execute(), LocalKsSelector::operator()(), LocalLambdaSelector::operator()(), utility::SecondaryVFit(), utility::SecondaryVFit_Lambda(), utility::SecondaryVFit_Lambdaref(), and utility::SecondaryVFitref().

decayLengthError() [2/3]

double SecondVertexFit::decayLengthError ( ) const

inline

Definition at line 51 of file InstallArea/x86_64-el9-gcc13-dbg/include/VertexFit/SecondVertexFit.h.

{ return m_lxyz_error; }

decayLengthError() [3/3]

double SecondVertexFit::decayLengthError ( ) const

inline

Definition at line 51 of file InstallArea/x86_64-el9-gcc13-opt/include/VertexFit/SecondVertexFit.h.

{ return m_lxyz_error; }

factor() [1/3]

double SecondVertexFit::factor ( ) const

inline

Definition at line 54 of file Analysis/VertexFit/include/VertexFit/SecondVertexFit.h.

{ return m_factor; }

Referenced by setMagCorrFactor().

factor() [2/3]

double SecondVertexFit::factor ( ) const

inline

Definition at line 54 of file InstallArea/x86_64-el9-gcc13-dbg/include/VertexFit/SecondVertexFit.h.

{ return m_factor; }

factor() [3/3]

double SecondVertexFit::factor ( ) const

inline

Definition at line 54 of file InstallArea/x86_64-el9-gcc13-opt/include/VertexFit/SecondVertexFit.h.

{ return m_factor; }

Fit() [1/3]

bool SecondVertexFit::Fit

(

)

Definition at line 82 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 DQAKsKpi::execute(), DQAKsKpiDEDX::execute(), inclks::execute(), incllambda::execute(), JsiLL::execute(), KsKpi::execute(), K0kk::MTotal(), K0kpi::MTotal(), K0pi0::MTotal(), K0pipi::MTotal(), K0pipipi0::MTotal(), LocalKsSelector::operator()(), LocalLambdaSelector::operator()(), utility::SecondaryVFit(), utility::SecondaryVFit_Lambda(), utility::SecondaryVFit_Lambdaref(), and utility::SecondaryVFitref().

Fit() [2/3]

bool SecondVertexFit::Fit

(

)

Fit() [3/3]

bool SecondVertexFit::Fit

(

)

init() [1/3]

void SecondVertexFit::init

(

)

Definition at line 65 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 DQAKsKpi::execute(), DQAKsKpiDEDX::execute(), inclks::execute(), incllambda::execute(), JsiLL::execute(), KsKpi::execute(), K0kk::MTotal(), K0kpi::MTotal(), K0pi0::MTotal(), K0pipi::MTotal(), K0pipipi0::MTotal(), LocalKsSelector::operator()(), LocalLambdaSelector::operator()(), utility::SecondaryVFit(), utility::SecondaryVFit_Lambda(), utility::SecondaryVFit_Lambdaref(), and utility::SecondaryVFitref().

init() [2/3]

void SecondVertexFit::init

(

)

init() [3/3]

void SecondVertexFit::init

(

)

instance() [1/3]

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(), DQAKsInfo::DQAKsInfo(), DQAKsKpi::execute(), DQAKsKpiDEDX::execute(), inclks::execute(), incllambda::execute(), JsiLL::execute(), KsKpi::execute(), K0kk::MTotal(), K0kpi::MTotal(), K0pi0::MTotal(), K0pipi::MTotal(), K0pipipi0::MTotal(), LocalKsSelector::operator()(), LocalLambdaSelector::operator()(), utility::SecondaryVFit(), utility::SecondaryVFit_Lambda(), utility::SecondaryVFit_Lambdaref(), and utility::SecondaryVFitref().

instance() [2/3]

SecondVertexFit * SecondVertexFit::instance ( )

static

instance() [3/3]

SecondVertexFit * SecondVertexFit::instance ( )

static

p4par() [1/3]

HepLorentzVector SecondVertexFit::p4par ( ) const

inline

Definition at line 56 of file Analysis/VertexFit/include/VertexFit/SecondVertexFit.h.

{ return m_p4par; }

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

p4par() [2/3]

HepLorentzVector SecondVertexFit::p4par ( ) const

inline

Definition at line 56 of file InstallArea/x86_64-el9-gcc13-dbg/include/VertexFit/SecondVertexFit.h.

{ return m_p4par; }

p4par() [3/3]

HepLorentzVector SecondVertexFit::p4par ( ) const

inline

Definition at line 56 of file InstallArea/x86_64-el9-gcc13-opt/include/VertexFit/SecondVertexFit.h.

{ return m_p4par; }

setChisqCut() [1/3]

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

inline

Definition at line 31 of file Analysis/VertexFit/include/VertexFit/SecondVertexFit.h.

                                                                          {
    m_chicut = chicut;
    m_chiter = chiter;
  }

setChisqCut() [2/3]

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

inline

Definition at line 31 of file InstallArea/x86_64-el9-gcc13-dbg/include/VertexFit/SecondVertexFit.h.

                                                                          {
    m_chicut = chicut;
    m_chiter = chiter;
  }

setChisqCut() [3/3]

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

inline

Definition at line 31 of file InstallArea/x86_64-el9-gcc13-opt/include/VertexFit/SecondVertexFit.h.

                                                                          {
    m_chicut = chicut;
    m_chiter = chiter;
  }

setIterNumber() [1/3]

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

inline

Definition at line 30 of file Analysis/VertexFit/include/VertexFit/SecondVertexFit.h.

{ m_niter = niter; }

setIterNumber() [2/3]

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

inline

Definition at line 30 of file InstallArea/x86_64-el9-gcc13-dbg/include/VertexFit/SecondVertexFit.h.

{ m_niter = niter; }

setIterNumber() [3/3]

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

inline

Definition at line 30 of file InstallArea/x86_64-el9-gcc13-opt/include/VertexFit/SecondVertexFit.h.

{ m_niter = niter; }

setMagCorrFactor() [1/3]

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

inline

Definition at line 43 of file Analysis/VertexFit/include/VertexFit/SecondVertexFit.h.

{ m_factor = factor; }

setMagCorrFactor() [2/3]

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

inline

Definition at line 43 of file InstallArea/x86_64-el9-gcc13-dbg/include/VertexFit/SecondVertexFit.h.

{ m_factor = factor; }

setMagCorrFactor() [3/3]

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

inline

Definition at line 43 of file InstallArea/x86_64-el9-gcc13-opt/include/VertexFit/SecondVertexFit.h.

{ m_factor = factor; }

setPrimaryVertex() [1/3]

void SecondVertexFit::setPrimaryVertex ( const VertexParameter vpar )

inline

Definition at line 38 of file Analysis/VertexFit/include/VertexFit/SecondVertexFit.h.

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

Referenced by DQAKsKpi::execute(), DQAKsKpiDEDX::execute(), KsKpi::execute(), K0kk::MTotal(), K0kpi::MTotal(), K0pi0::MTotal(), K0pipi::MTotal(), K0pipipi0::MTotal(), LocalKsSelector::operator()(), LocalLambdaSelector::operator()(), utility::SecondaryVFit(), utility::SecondaryVFit_Lambda(), utility::SecondaryVFit_Lambdaref(), and utility::SecondaryVFitref().

setPrimaryVertex() [2/3]

void SecondVertexFit::setPrimaryVertex ( const VertexParameter vpar )

inline

Definition at line 38 of file InstallArea/x86_64-el9-gcc13-dbg/include/VertexFit/SecondVertexFit.h.

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

setPrimaryVertex() [3/3]

void SecondVertexFit::setPrimaryVertex ( const VertexParameter vpar )

inline

Definition at line 38 of file InstallArea/x86_64-el9-gcc13-opt/include/VertexFit/SecondVertexFit.h.

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

setVpar() [1/3]

void SecondVertexFit::setVpar ( const VertexParameter vpar )

inline

Definition at line 36 of file Analysis/VertexFit/include/VertexFit/SecondVertexFit.h.

{ m_vpar_secondary = vpar; }

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

setVpar() [2/3]

void SecondVertexFit::setVpar ( const VertexParameter vpar )

inline

Definition at line 36 of file InstallArea/x86_64-el9-gcc13-dbg/include/VertexFit/SecondVertexFit.h.

{ m_vpar_secondary = vpar; }

setVpar() [3/3]

void SecondVertexFit::setVpar ( const VertexParameter vpar )

inline

Definition at line 36 of file InstallArea/x86_64-el9-gcc13-opt/include/VertexFit/SecondVertexFit.h.

{ m_vpar_secondary = vpar; }

vpar() [1/3]

VertexParameter SecondVertexFit::vpar ( ) const

inline

Definition at line 58 of file Analysis/VertexFit/include/VertexFit/SecondVertexFit.h.

{ return m_vpar_secondary; }

Referenced by setPrimaryVertex(), and setVpar().

vpar() [2/3]

VertexParameter SecondVertexFit::vpar ( ) const

inline

Definition at line 58 of file InstallArea/x86_64-el9-gcc13-dbg/include/VertexFit/SecondVertexFit.h.

{ return m_vpar_secondary; }

vpar() [3/3]

VertexParameter SecondVertexFit::vpar ( ) const

inline

Definition at line 58 of file InstallArea/x86_64-el9-gcc13-opt/include/VertexFit/SecondVertexFit.h.

{ return m_vpar_secondary; }

wpar() [1/3]

WTrackParameter SecondVertexFit::wpar ( ) const

inline

Definition at line 57 of file Analysis/VertexFit/include/VertexFit/SecondVertexFit.h.

{ return m_wtrk; }

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

wpar() [2/3]

WTrackParameter SecondVertexFit::wpar ( ) const

inline

Definition at line 57 of file InstallArea/x86_64-el9-gcc13-dbg/include/VertexFit/SecondVertexFit.h.

{ return m_wtrk; }

wpar() [3/3]

WTrackParameter SecondVertexFit::wpar ( ) const

inline

Definition at line 57 of file InstallArea/x86_64-el9-gcc13-opt/include/VertexFit/SecondVertexFit.h.

{ return m_wtrk; }

---

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

• Analysis/VertexFit/include/VertexFit/SecondVertexFit.h
• InstallArea/x86_64-el9-gcc13-dbg/include/VertexFit/SecondVertexFit.h
• InstallArea/x86_64-el9-gcc13-opt/include/VertexFit/SecondVertexFit.h
• SecondVertexFit.cxx