TFitParticleSpher Class Reference

#include <TFitParticleSpher.h>

Inheritance diagram for TFitParticleSpher:

Public Member Functions
virtual TLorentzVector *	calc4Vec (const TMatrixD *params)
virtual TAbsFitParticle *	clone (const TString &newname=TString("")) const
virtual TMatrixD *	getDerivative ()
virtual void	setIni4Vec (const TLorentzVector *pini)
	TFitParticleSpher ()
	TFitParticleSpher (const TFitParticleSpher &fitParticle)
	TFitParticleSpher (TLorentzVector *pini, const TMatrixD *theCovMatrix)
	TFitParticleSpher (const TString &name, const TString &title, TLorentzVector *pini, const TMatrixD *theCovMatrix)
virtual TMatrixD *	transform (const TLorentzVector &vec)
virtual	~TFitParticleSpher ()
Public Member Functions inherited from TAbsFitParticle
virtual void	applycorr (TMatrixD *corrMatrix)
virtual const TMatrixD *	getCovMatrix () const
virtual const TMatrixD *	getCovMatrixDeltaY ()
virtual const TMatrixD *	getCovMatrixFit () const
const TLorentzVector *	getCurr4Vec ()
TString	getInfoString ()
const TLorentzVector *	getIni4Vec ()
Int_t	getNPar () const
const TMatrixD *	getParCurr ()
const TMatrixD *	getParIni ()
virtual const TMatrixD *	getPull ()
const TVector3 *	getu1 ()
const TVector3 *	getu2 ()
const TVector3 *	getu3 ()
virtual void	print ()
virtual void	reset ()
virtual void	setCovMatrix (const TMatrixD *theCovMatrix)
virtual void	setCovMatrixFit (const TMatrixD *theCovMatrixFit)
void	setParIni (const TMatrixD *parini)
	TAbsFitParticle ()
	TAbsFitParticle (const TString &name, const TString &title)
virtual	~TAbsFitParticle ()

Protected Member Functions
void	init (TLorentzVector *pini, const TMatrixD *theCovMatrix)
Protected Member Functions inherited from TAbsFitParticle
void	calcCovMatrixDeltaY ()

Additional Inherited Members
Protected Attributes inherited from TAbsFitParticle
TMatrixD	_covMatrix
TMatrixD	_covMatrixDeltaY
TMatrixD	_covMatrixFit
TMatrixD	_iniparameters
Int_t	_nPar
TMatrixD	_parameters
TLorentzVector	_pcurr
TLorentzVector	_pini
TMatrixD	_pull
TVector3	_u1
TVector3	_u2
TVector3	_u3

Detailed Description

Definition at line 10 of file TFitParticleSpher.h.

Constructor & Destructor Documentation

TFitParticleSpher::TFitParticleSpher

(

)

Definition at line 27 of file TFitParticleSpher.cc.

References init().

Referenced by clone().

  :TAbsFitParticle()  
{
  init(0, 0);
}

TFitParticleSpher::TFitParticleSpher

(

const TFitParticleSpher &

fitParticle

)

Definition at line 33 of file TFitParticleSpher.cc.

References TAbsFitParticle::_covMatrix, TAbsFitParticle::_iniparameters, TAbsFitParticle::_nPar, TAbsFitParticle::_parameters, TAbsFitParticle::_pcurr, TAbsFitParticle::_pini, TAbsFitParticle::_u1, TAbsFitParticle::_u2, and TAbsFitParticle::_u3.

  :TAbsFitParticle( fitParticle.GetName(), fitParticle.GetTitle() )
{

  _nPar = fitParticle._nPar;
  _u1 = fitParticle._u1;
  _u2 = fitParticle._u2;
  _u3 = fitParticle._u3;
  _covMatrix.ResizeTo(  fitParticle._covMatrix );
  _covMatrix = fitParticle._covMatrix;
  _iniparameters.ResizeTo( fitParticle._iniparameters );
  _iniparameters = fitParticle._iniparameters;
  _parameters.ResizeTo( fitParticle._parameters );
  _parameters = fitParticle._parameters;
  _pini = fitParticle._pini;
  _pcurr = fitParticle._pcurr;

}

TFitParticleSpher::TFitParticleSpher	(	TLorentzVector *	pini,
		const TMatrixD *	theCovMatrix
	)

Definition at line 52 of file TFitParticleSpher.cc.

References init().

  :TAbsFitParticle()  
{
  init(pini, theCovMatrix);
}

TFitParticleSpher::TFitParticleSpher	(	const TString &	name,
		const TString &	title,
		TLorentzVector *	pini,
		const TMatrixD *	theCovMatrix
	)

Definition at line 58 of file TFitParticleSpher.cc.

References init().

  :TAbsFitParticle(name, title)  
{
  init(pini, theCovMatrix);
}

TFitParticleSpher::~TFitParticleSpher ( )

virtual

Definition at line 77 of file TFitParticleSpher.cc.

                                      {

}

Member Function Documentation

TLorentzVector * TFitParticleSpher::calc4Vec ( const TMatrixD *  params )

virtual

Implements TAbsFitParticle.

Definition at line 92 of file TFitParticleSpher.cc.

References TAbsFitParticle::_nPar, TAbsFitParticle::_pini, ztail::d, phi, alignCSCRings::r, theta(), X, and Gflash::Z.

                                                                    {
  // Calculates a 4vector corresponding to the given
  // parameter values

  if (params == 0) {
    return 0;
  }

  if ( params->GetNcols() != 1 || params->GetNrows() !=_nPar ) {
    edm::LogError ("WrongMatrixSize")
      << GetName() << "::calc4Vec - Parameter matrix has wrong size.";
    return 0;
  }

  Double_t r = (*params)(0,0);
  Double_t theta = (*params)(1,0);
  Double_t phi = (*params)(2,0);
  Double_t d = (*params)(3,0);

  Double_t X = r*TMath::Cos(phi)*TMath::Sin(theta);
  Double_t Y = r*TMath::Sin(phi)*TMath::Sin(theta);
  Double_t Z = r*TMath::Cos(theta);
  Double_t E = TMath::Sqrt( X*X + Y*Y + Z*Z + d*d*_pini.M2() );
        
  TLorentzVector* vec = new TLorentzVector( X, Y, Z, E );
  return vec;

}

TAbsFitParticle * TFitParticleSpher::clone ( const TString &  newname = TString("") ) const

virtual

Implements TAbsFitParticle.

Definition at line 65 of file TFitParticleSpher.cc.

References TFitParticleSpher().

                                                                        {
  // Returns a copy of itself
  
  TAbsFitParticle* myclone = new TFitParticleSpher( *this );
  if ( newname.Length() > 0 ) myclone->SetName(newname);
  return myclone;

}

TMatrixD * TFitParticleSpher::getDerivative ( )

virtual

Implements TAbsFitParticle.

Definition at line 171 of file TFitParticleSpher.cc.

References TAbsFitParticle::_parameters, TAbsFitParticle::_pcurr, TAbsFitParticle::_pini, ztail::d, phi, alignCSCRings::r, and theta().

                                           {
  // returns derivative dP/dy with P=(p,E) and y=(r, theta, phi, d) 
  // the free parameters of the fit. The columns of the matrix contain 
  // (dP/dr, dP/dtheta, ...).

  TMatrixD* DerivativeMatrix = new TMatrixD(4,4);
  (*DerivativeMatrix) *= 0.;

  Double_t r = _parameters(0,0);
  Double_t theta = _parameters(1,0);
  Double_t phi = _parameters(2,0);
  Double_t d = _parameters(3,0);

  //1st column: dP/dr
  (*DerivativeMatrix)(0,0) = TMath::Cos(phi)*TMath::Sin(theta);
  (*DerivativeMatrix)(1,0) = TMath::Sin(phi)*TMath::Sin(theta);
  (*DerivativeMatrix)(2,0) = TMath::Cos(theta);
  (*DerivativeMatrix)(3,0) = 0.;

  //2nd column: dP/dtheta
  (*DerivativeMatrix)(0,1) = r*TMath::Cos(phi)*TMath::Cos(theta);
  (*DerivativeMatrix)(1,1) = r*TMath::Sin(phi)*TMath::Cos(theta);
  (*DerivativeMatrix)(2,1) = -1.*r*TMath::Sin(theta);
  (*DerivativeMatrix)(3,1) = 0.;

   //3rd column: dP/dphi
  (*DerivativeMatrix)(0,2) = -1.*r*TMath::Sin(phi)*TMath::Sin(theta);
  (*DerivativeMatrix)(1,2) = r*TMath::Cos(phi)*TMath::Sin(theta);;
  (*DerivativeMatrix)(2,2) = 0.;
  (*DerivativeMatrix)(3,2) = 0.;

   //4th column: dP/dm
  (*DerivativeMatrix)(0,3) = 0.;
  (*DerivativeMatrix)(1,3) = 0.;
  (*DerivativeMatrix)(2,3) = 0.;
  (*DerivativeMatrix)(3,3) = _pini.M()*_pini.M()*d/_pcurr.E();

  return DerivativeMatrix;

}

void TFitParticleSpher::init ( TLorentzVector *  pini, const TMatrixD *  theCovMatrix  )

protected

Definition at line 84 of file TFitParticleSpher.cc.

References TAbsFitParticle::_nPar, TAbsFitParticle::setCovMatrix(), and setIni4Vec().

Referenced by TFitParticleSpher().

                                                                                {

  _nPar = 4;
  setIni4Vec(pini);
  setCovMatrix(theCovMatrix);

}

void TFitParticleSpher::setIni4Vec ( const TLorentzVector *  pini )

virtual

Implements TAbsFitParticle.

Definition at line 121 of file TFitParticleSpher.cc.

References TAbsFitParticle::_iniparameters, TAbsFitParticle::_nPar, TAbsFitParticle::_parameters, TAbsFitParticle::_pcurr, TAbsFitParticle::_pini, TAbsFitParticle::_u1, TAbsFitParticle::_u2, TAbsFitParticle::_u3, phi, alignCSCRings::r, and theta().

Referenced by init().

                                                             {
  // Set the initial 4vector. Will also set the 
  // inital parameter values 

  if (pini == 0) {

    _u1.SetXYZ(0., 0., 0.);
    _u3.SetXYZ(0., 0., 0.);
    _u2.SetXYZ(0., 0., 0.);
    _pini.SetXYZT(0., 0., 0., 0.);
    _pcurr = _pini;

    _iniparameters.ResizeTo(_nPar,1);
    _iniparameters(0,0) = 0.;
    _iniparameters(1,0) = 0.;
    _iniparameters(2,0) = 0.;
    _iniparameters(3,0) = 1.;
    
    _parameters.ResizeTo(_nPar,1);
    _parameters(0,0) = 0.;
    _parameters(1,0) = 0.;
    _parameters(2,0) = 0.;   
    _parameters(3,0) = 1.;   
    
  } else {
    
    Double_t r = pini->P();
    Double_t theta = pini->Theta();
    Double_t phi = pini->Phi();
    
    _pini = (*pini);
    _pcurr = _pini;
    
    _iniparameters.ResizeTo(_nPar,1);
    _iniparameters(0,0) = r;
    _iniparameters(1,0) = theta;
    _iniparameters(2,0) = phi;
    _iniparameters(3,0) = 1.;
    
    _parameters.ResizeTo(_nPar,1);
    _parameters = _iniparameters;

    _u1.SetXYZ( TMath::Cos(phi)*TMath::Sin(theta), TMath::Sin(phi)*TMath::Sin(theta), TMath::Cos(theta) );
    _u2.SetXYZ( TMath::Cos(phi)*TMath::Cos(theta), TMath::Sin(phi)*TMath::Cos(theta), -1.*TMath::Sin(theta) );
    _u3.SetXYZ( -1.*TMath::Sin(phi), TMath::Cos(phi), 0. );

  }

}

TMatrixD * TFitParticleSpher::transform ( const TLorentzVector &  vec )

virtual

Implements TAbsFitParticle.

Definition at line 212 of file TFitParticleSpher.cc.

References TAbsFitParticle::_nPar, and TAbsFitParticle::_pini.

                                                                {
  // Returns the parameters corresponding to the given 
  // 4vector

  // retrieve parameters
  TMatrixD* tparams = new TMatrixD( _nPar, 1 );
  (*tparams)(0,0) = vec.P();
  (*tparams)(1,0) = vec.Theta();
  (*tparams)(2,0) = vec.Phi();
  (*tparams)(3,0) = vec.M()/_pini.M();

  return tparams;

}