12 TKinFitter fitter(
"CandMassFit",
"CandMassFit");
15 vector<TMatrixD>
errors(daus, TMatrix(3,3));
16 vector<TVector3> momenta(daus);
17 vector<TFitParticleMCCart *>
particles(daus, 0);
19 "MassConstraint", 0, 0 ,
mass_);
20 for (
size_t i = 0;
i < daus; ++
i ) {
24 TVector3 & mom = momenta[
i];
25 mom = TVector3( p4.px(), p4.py(), p4.pz() );
38 fitter.setMaxNbIter( 30 );
39 fitter.setMaxDeltaS( 1
e-2 );
40 fitter.setMaxF( 1
e-1 );
41 fitter.setVerbosity( 0 );
44 TLorentzVector sum( 0, 0, 0, 0 );
45 for(
size_t i = 0;
i < daus; ++
i ) {
54 return FitQuality( fitter.getS(), fitter.getNDF());
virtual void setP4(const LorentzVector &p4)=0
set 4-momentum
virtual const Candidate * daughter(size_type i) const =0
return daughter at a given position, i = 0, ... numberOfDaughters() - 1 (read only mode) ...
virtual const LorentzVector & p4() const =0
four-momentum Lorentz vector
virtual double mass() const =0
mass
const TLorentzVector * getCurr4Vec()
T get() const
get a component
virtual size_type numberOfDaughters() const =0
number of daughters
math::PtEtaPhiELorentzVectorF LorentzVector