How to use the DecayTreeFitter tools

There is a twiki here

Here is an example below for applying a J/psi constraint

 1000Particle* part;
 1010DecayTreeFitter::Fitter fitter(**iterC);
 1020fitter.setMassConstraint(LHCb::ParticleID(443));
 1030fitter.fit();
 1040
 1050// getting some of the results
 1060if (DecayTreeFitter::Fitter::Success == fitter.status()){
 1070     LHCb::DecayTree tree = fitter.getFittedTree();
 1080     const LHCb::Particle* head = tree.head();
 1090     mReFit = M(head); 
 1100     kineChi2 = fitter.chiSquare() ; 
 1110     kinedof = fitter.nDof();
 1120     kstatus = fitter.status();
 1130     emReFit2 = Gaudi::Math::sigmamass(head->momentum(), head->momCovMatrix());
 1140}
 1150else {
 1160    kstatus = fitter.status();
 1170} 
 1180

-- MatthewNeedham - 14-Oct-2010

Edit | Attach | Watch | Print version | History: r1 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r1 - 2010-10-14 - MatthewNeedham
 
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    Main All webs login

This site is powered by the TWiki collaboration platform Powered by PerlCopyright & 2008-2021 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
or Ideas, requests, problems regarding TWiki? use Discourse or Send feedback