SiPixelLAHarvestMCS::FitFPixMuH Class Reference

Inheritance diagram for SiPixelLAHarvestMCS::FitFPixMuH:

Public Member Functions
void	add (const FPixMuH &fmh)
int	fit (double initVal)
	FitFPixMuH ()
double	getMuH ()
double	getMuHErr ()
unsigned int	size ()
	~FitFPixMuH () override=default

Private Member Functions
double	calcChi2 (double par_0)

Private Attributes
std::vector< FPixMuH >	cmvar
TVirtualFitter *	minuit
double	muH
double	muHErr

Friends
void	fcn_func (Int_t &npar, Double_t *deriv, Double_t &f, Double_t *par, Int_t flag)

Detailed Description

Definition at line 77 of file SiPixelLorentzAnglePCLHarvesterMCS.cc.

Constructor & Destructor Documentation

FitFPixMuH()

SiPixelLAHarvestMCS::FitFPixMuH::FitFPixMuH ( )

inline

Definition at line 79 of file SiPixelLorentzAnglePCLHarvesterMCS.cc.

: muH(0), muHErr(0) {}

~FitFPixMuH()

SiPixelLAHarvestMCS::FitFPixMuH::~FitFPixMuH ( )

overridedefault

Member Function Documentation

add()

void SiPixelLAHarvestMCS::FitFPixMuH::add ( const FPixMuH &  fmh )

inline

Definition at line 86 of file SiPixelLorentzAnglePCLHarvesterMCS.cc.

References cmvar.

Referenced by SiPixelLorentzAnglePCLHarvesterMCS::dqmEndJob(), counter.Counter::register(), SequenceTypes._TaskBase::remove(), and SequenceTypes._TaskBase::replace().

{ cmvar.push_back(fmh); }

calcChi2()

double SiPixelLAHarvestMCS::FitFPixMuH::calcChi2 ( double  par_0 )

inlineprivate

Definition at line 116 of file SiPixelLorentzAnglePCLHarvesterMCS.cc.

References nano_mu_local_reco_cff::chi2, cmvar, mps_fire::i, dqmiodumpmetadata::n, and findQualityFiles::v.

                                  {
      double tanlorpertesla = par_0;
      double tlpt2 = tanlorpertesla * tanlorpertesla;
      double v[3], xshift, yshift;

      int n = cmvar.size();

      double chi2 = 0.0;
      for (int i = 0; i < n; i++) {
        v[0] = -(tanlorpertesla * cmvar[i].bfield[1] + tlpt2 * cmvar[i].bfield[0] * cmvar[i].bfield[2]);
        v[1] = tanlorpertesla * cmvar[i].bfield[0] - tlpt2 * cmvar[i].bfield[1] * cmvar[i].bfield[2];
        v[2] = -(1. + tlpt2 * cmvar[i].bfield[2] * cmvar[i].bfield[2]);

        xshift = v[0] / v[2];
        yshift = v[1] / v[2];

        chi2 += (xshift - cmvar[i].shiftx) * (xshift - cmvar[i].shiftx) / cmvar[i].shiftx_err / cmvar[i].shiftx_err +
                (yshift - cmvar[i].shifty) * (yshift - cmvar[i].shifty) / cmvar[i].shifty_err / cmvar[i].shifty_err;
      }
      return chi2;
    }

fit()

int SiPixelLAHarvestMCS::FitFPixMuH::fit ( double  initVal )

inline

Definition at line 87 of file SiPixelLorentzAnglePCLHarvesterMCS.cc.

References fcn_func, HLTSiStripMonitoring_cff::Fitter, minuit, muH, muHErr, mps_update::status, and up.

Referenced by SiPixelLorentzAnglePCLHarvesterMCS::dqmEndJob(), and trackingPlots.Iteration::modules().

                            {
      minuit = TVirtualFitter::Fitter(this, 1);
      minuit->SetFCN(fcn_func);
      double tanlorpertesla = initVal;
      minuit->SetParameter(0, "muH", tanlorpertesla, 0.1, 0., 0.);

      double arglist[100];
      arglist[0] = 3.;
      minuit->ExecuteCommand("SET PRINT", arglist, 0);
      double up = 1.0;
      minuit->SetErrorDef(up);
      arglist[0] = 100000.;
      int status = minuit->ExecuteCommand("MIGRAD", arglist, 0);
      muH = minuit->GetParameter(0);
      muHErr = minuit->GetParError(0);

      return status;
    }

getMuH()

double SiPixelLAHarvestMCS::FitFPixMuH::getMuH ( )

inline

Definition at line 106 of file SiPixelLorentzAnglePCLHarvesterMCS.cc.

References muH.

Referenced by SiPixelLorentzAnglePCLHarvesterMCS::dqmEndJob().

{ return muH; }

getMuHErr()

double SiPixelLAHarvestMCS::FitFPixMuH::getMuHErr ( )

inline

Definition at line 107 of file SiPixelLorentzAnglePCLHarvesterMCS.cc.

References muHErr.

Referenced by SiPixelLorentzAnglePCLHarvesterMCS::dqmEndJob().

{ return muHErr; }

size()

unsigned int SiPixelLAHarvestMCS::FitFPixMuH::size ( void  )

inline

Definition at line 108 of file SiPixelLorentzAnglePCLHarvesterMCS.cc.

References cmvar.

Referenced by ntupleDataFormat._Collection::__iter__(), ntupleDataFormat._Collection::__len__(), and SiPixelLorentzAnglePCLHarvesterMCS::dqmEndJob().

{ return cmvar.size(); }

Friends And Related Function Documentation

fcn_func

void fcn_func ( Int_t &  npar, Double_t *  deriv, Double_t &  f, Double_t *  par, Int_t  flag  )

friend

Definition at line 139 of file SiPixelLorentzAnglePCLHarvesterMCS.cc.

Referenced by fit().

                                                                                      {
    f = ((dynamic_cast<FitFPixMuH*>((TVirtualFitter::GetFitter())->GetObjectFit()))->calcChi2(par[0]));
  }

Member Data Documentation

cmvar

std::vector<FPixMuH> SiPixelLAHarvestMCS::FitFPixMuH::cmvar

private

Definition at line 112 of file SiPixelLorentzAnglePCLHarvesterMCS.cc.

Referenced by add(), calcChi2(), and size().

minuit

TVirtualFitter* SiPixelLAHarvestMCS::FitFPixMuH::minuit

private

Definition at line 111 of file SiPixelLorentzAnglePCLHarvesterMCS.cc.

Referenced by fit().

muH

double SiPixelLAHarvestMCS::FitFPixMuH::muH

private

Definition at line 113 of file SiPixelLorentzAnglePCLHarvesterMCS.cc.

Referenced by fit(), and getMuH().

muHErr

double SiPixelLAHarvestMCS::FitFPixMuH::muHErr

private

Definition at line 114 of file SiPixelLorentzAnglePCLHarvesterMCS.cc.

Referenced by fit(), and getMuHErr().