Alignment/MuonAlignmentAlgorithms/interface/MuonResidualsPositionFitter.h - master - githubmirror/cmssw

Alignment/MuonAlignmentAlgorithms/interface/MuonResidualsPositionFitter.h

43 строки2 KB

Cms Build

Clang-Format

16 май 2019, 00:01

16 май 2019, 00:0192500ad

#ifndef Alignment_MuonAlignmentAlgorithms_MuonResidualsPositionFitter_H
#define Alignment_MuonAlignmentAlgorithms_MuonResidualsPositionFitter_H

/** \class MuonResidualsPositionFitter
 *  $Date: 2010/03/12 22:23:34 $
 *  $Revision: 1.8 $
 *  \author J. Pivarski - Texas A&M University <pivarski@physics.tamu.edu>
 */

#include "Alignment/MuonAlignmentAlgorithms/interface/MuonResidualsFitter.h"

class MuonResidualsPositionFitter : public MuonResidualsFitter {
public:
  enum { kPosition = 0, kZpos, kPhiz, kScattering, kSigma, kGamma, kNPar };

  enum { kResidual = 0, kAngleError, kTrackAngle, kTrackPosition, kNData };

  MuonResidualsPositionFitter(int residualsModel, int minHits, int useResiduals, bool weightAlignment = true)
      : MuonResidualsFitter(residualsModel, minHits, useResiduals, weightAlignment) {}

  int type() const override { return MuonResidualsFitter::kPositionFitter; }

  int npar() override {
    if (residualsModel() == kPureGaussian || residualsModel() == kGaussPowerTails)
      return kNPar - 1;
    else if (residualsModel() == kPowerLawTails)
      return kNPar;
    else if (residualsModel() == kROOTVoigt)
      return kNPar;
    else
      assert(false);
  }
  int ndata() override { return kNData; }

  bool fit(Alignable *ali) override;
  double sumofweights() override { return numResiduals(); }
  double plot(std::string name, TFileDirectory *dir, Alignable *ali) override;

protected:
  void inform(TMinuit *tMinuit) override;
};

#endif  // Alignment_MuonAlignmentAlgorithms_MuonResidualsPositionFitter_H