#include <TFitParticleSpher.h>
Public Member Functions | |
virtual TLorentzVector * | calc4Vec (const TMatrixD *params) |
virtual TAbsFitParticle * | clone (TString newname="") const |
virtual TMatrixD * | getDerivative () |
virtual void | setIni4Vec (const TLorentzVector *pini) |
TFitParticleSpher () | |
TFitParticleSpher (TLorentzVector *pini, const TMatrixD *theCovMatrix) | |
TFitParticleSpher (const TFitParticleSpher &fitParticle) | |
TFitParticleSpher (const TString &name, const TString &title, TLorentzVector *pini, const TMatrixD *theCovMatrix) | |
virtual TMatrixD * | transform (const TLorentzVector &vec) |
virtual | ~TFitParticleSpher () |
Protected Member Functions | |
void | init (TLorentzVector *pini, const TMatrixD *theCovMatrix) |
Definition at line 10 of file TFitParticleSpher.h.
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.
{ }
TLorentzVector * TFitParticleSpher::calc4Vec | ( | const TMatrixD * | params | ) | [virtual] |
Implements TAbsFitParticle.
Definition at line 92 of file TFitParticleSpher.cc.
References TAbsFitParticle::_nPar, TAbsFitParticle::_pini, 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 | ( | TString | newname = "" | ) | 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, 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.