CMS 3D CMS Logo

List of all members | Public Types | Public Member Functions | Protected Member Functions
MuonResidualsPositionFitter Class Reference

#include <MuonResidualsPositionFitter.h>

Inheritance diagram for MuonResidualsPositionFitter:
MuonResidualsFitter

Public Types

enum  {
  kPosition = 0, kZpos, kPhiz, kScattering,
  kSigma, kGamma, kNPar
}
 
enum  {
  kResidual = 0, kAngleError, kTrackAngle, kTrackPosition,
  kNData
}
 
- Public Types inherited from MuonResidualsFitter
enum  {
  kPureGaussian, kPowerLawTails, kROOTVoigt, kGaussPowerTails,
  kPureGaussian2D
}
 
enum  {
  k1DOF, k5DOF, k6DOF, k6DOFrphi,
  kPositionFitter, kAngleFitter, kAngleBfieldFitter
}
 
enum  {
  k1111, k1110, k1100, k1010,
  k0010, k1000, k0100
}
 

Public Member Functions

bool fit (Alignable *ali) override
 
 MuonResidualsPositionFitter (int residualsModel, int minHits, int useResiduals, bool weightAlignment=true)
 
int ndata () override
 
int npar () override
 
double plot (std::string name, TFileDirectory *dir, Alignable *ali) override
 
double sumofweights () override
 
int type () const override
 
- Public Member Functions inherited from MuonResidualsFitter
void computeHistogramRangeAndBinning (int which, int &nbins, double &a, double &b)
 
virtual void correctBField ()=0
 
virtual void correctBField (int idx_momentum, int idx_q)
 
TMatrixDSym correlationMatrix ()
 
double covarianceElement (int parNum1, int parNum2)
 
TMatrixDSym covarianceMatrix ()
 
void eraseNotSelectedResiduals ()
 
double errorerror (int parNum)
 
void fiducialCuts (double xMin=-80.0, double xMax=80.0, double yMin=-80.0, double yMax=80.0, bool fidcut1=false)
 
void fill (double *residual)
 
void fix (int parNum, bool dofix=true)
 
bool fixed (int parNum)
 
void histogramChi2GaussianFit (int which, double &fit_mean, double &fit_sigma)
 
double loglikelihood ()
 
 MuonResidualsFitter (int residualsModel, int minHits, int useResiduals, bool weightAlignment=true)
 
int nfixed ()
 
long numResiduals () const
 
long numsegments ()
 
int parNum2parIdx (int parNum)
 
void plotsimple (std::string name, TFileDirectory *dir, int which, double multiplier)
 
void plotweighted (std::string name, TFileDirectory *dir, int which, int whichredchi2, double multiplier)
 
void read (FILE *file, int which=0)
 
std::vector< double * >::const_iterator residuals_begin () const
 
std::vector< double * >::const_iterator residuals_end () const
 
int residualsModel () const
 
std::vector< bool > & selectedResidualsFlags ()
 
void selectPeakResiduals (double nsigma, int nvar, int *vars)
 
void selectPeakResiduals_simple (double nsigma, int nvar, int *vars)
 
void setInitialValue (int parNum, double value)
 
void setPrintLevel (int printLevel)
 
void setStrategy (int strategy)
 
int useRes (int pattern=-1)
 
double value (int parNum)
 
void write (FILE *file, int which=0)
 
virtual ~MuonResidualsFitter ()
 

Protected Member Functions

void inform (TMinuit *tMinuit) override
 
- Protected Member Functions inherited from MuonResidualsFitter
bool dofit (void(*fcn)(int &, double *, double &, double *, int), std::vector< int > &parNum, std::vector< std::string > &parName, std::vector< double > &start, std::vector< double > &step, std::vector< double > &low, std::vector< double > &high)
 
void initialize_table ()
 

Additional Inherited Members

- Protected Attributes inherited from MuonResidualsFitter
double m_center [20]
 
TMatrixDSym m_cov
 
std::vector< double > m_error
 
std::vector< bool > m_fixed
 
double m_loglikelihood
 
int m_minHits
 
std::map< int, double > m_parNum2InitValue
 
std::map< int, int > m_parNum2parIdx
 
int m_printLevel
 
double m_radii [20]
 
std::vector< double * > m_residuals
 
std::vector< bool > m_residuals_ok
 
int m_residualsModel
 
int m_strategy
 
int m_useResiduals
 
std::vector< double > m_value
 
bool m_weightAlignment
 

Detailed Description

Date
2010/03/12 22:23:34
Revision
1.8
Author
J. Pivarski - Texas A&M University pivar.nosp@m.ski@.nosp@m.physi.nosp@m.cs.t.nosp@m.amu.e.nosp@m.du

Definition at line 12 of file MuonResidualsPositionFitter.h.

Member Enumeration Documentation

◆ anonymous enum

anonymous enum
Enumerator
kPosition 
kZpos 
kPhiz 
kScattering 
kSigma 
kGamma 
kNPar 

Definition at line 14 of file MuonResidualsPositionFitter.h.

◆ anonymous enum

anonymous enum
Enumerator
kResidual 
kAngleError 
kTrackAngle 
kTrackPosition 
kNData 

Definition at line 16 of file MuonResidualsPositionFitter.h.

Constructor & Destructor Documentation

◆ MuonResidualsPositionFitter()

MuonResidualsPositionFitter::MuonResidualsPositionFitter ( int  residualsModel,
int  minHits,
int  useResiduals,
bool  weightAlignment = true 
)
inline

Member Function Documentation

◆ fit()

bool MuonResidualsPositionFitter::fit ( Alignable ali)
overridevirtual

Implements MuonResidualsFitter.

Definition at line 42 of file MuonResidualsPositionFitter.cc.

42  {
43  initialize_table(); // if not already initialized
44 
45  double sum_x = 0.;
46  double sum_xx = 0.;
47  int N = 0;
48 
49  for (std::vector<double *>::const_iterator resiter = residuals_begin(); resiter != residuals_end(); ++resiter) {
50  const double residual = (*resiter)[kResidual];
51  // const double angleerror = (*resiter)[kAngleError];
52  // const double trackangle = (*resiter)[kTrackAngle];
53  // const double trackposition = (*resiter)[kTrackPosition];
54 
55  if (fabs(residual) < 10.) { // truncate at 100 mm
56  sum_x += residual;
57  sum_xx += residual * residual;
58  N++;
59  }
60  }
61 
62  if (N < m_minHits)
63  return false;
64 
65  // truncated mean and stdev to seed the fit
66  double mean = sum_x / double(N);
67  double stdev = sqrt(sum_xx / double(N) - pow(sum_x / double(N), 2));
68 
69  // refine the standard deviation calculation
70  sum_x = 0.;
71  sum_xx = 0.;
72  N = 0;
73  for (std::vector<double *>::const_iterator resiter = residuals_begin(); resiter != residuals_end(); ++resiter) {
74  const double residual = (*resiter)[kResidual];
75  if (mean - 1.5 * stdev < residual && residual < mean + 1.5 * stdev) {
76  sum_x += residual;
77  sum_xx += residual * residual;
78  N++;
79  }
80  }
81  mean = sum_x / double(N);
82  stdev = sqrt(sum_xx / double(N) - pow(sum_x / double(N), 2));
83 
84  sum_x = 0.;
85  sum_xx = 0.;
86  N = 0;
87  for (std::vector<double *>::const_iterator resiter = residuals_begin(); resiter != residuals_end(); ++resiter) {
88  const double residual = (*resiter)[kResidual];
89  if (mean - 1.5 * stdev < residual && residual < mean + 1.5 * stdev) {
90  sum_x += residual;
91  sum_xx += residual * residual;
92  N++;
93  }
94  }
95  mean = sum_x / double(N);
96  stdev = sqrt(sum_xx / double(N) - pow(sum_x / double(N), 2));
97 
98  std::vector<int> parNum;
99  std::vector<std::string> parName;
100  std::vector<double> start;
101  std::vector<double> step;
102  std::vector<double> low;
103  std::vector<double> high;
104 
105  parNum.push_back(kPosition);
106  parName.push_back(std::string("position"));
107  start.push_back(mean);
108  step.push_back(0.1);
109  low.push_back(0.);
110  high.push_back(0.);
111  parNum.push_back(kZpos);
112  parName.push_back(std::string("zpos"));
113  start.push_back(0.);
114  step.push_back(0.1);
115  low.push_back(0.);
116  high.push_back(0.);
117  parNum.push_back(kPhiz);
118  parName.push_back(std::string("phiz"));
119  start.push_back(0.);
120  step.push_back(0.1);
121  low.push_back(0.);
122  high.push_back(0.);
123  parNum.push_back(kScattering);
124  parName.push_back(std::string("scattering"));
125  start.push_back(0.);
126  step.push_back(0.1 * 1000.);
127  low.push_back(0.);
128  high.push_back(0.);
129  parNum.push_back(kSigma);
130  parName.push_back(std::string("sigma"));
131  start.push_back(stdev);
132  step.push_back(0.1 * stdev);
133  low.push_back(0.);
134  high.push_back(0.);
136  parNum.push_back(kGamma);
137  parName.push_back(std::string("gamma"));
138  start.push_back(stdev);
139  step.push_back(0.1 * stdev);
140  low.push_back(0.);
141  high.push_back(0.);
142  }
143 
144  return dofit(&MuonResidualsPositionFitter_FCN, parNum, parName, start, step, low, high);
145 }

References MuonResidualsFitter::dofit(), LaserClient_cfi::high, MuonResidualsFitter::initialize_table(), kGamma, MuonResidualsFitter::kGaussPowerTails, kPhiz, kPosition, MuonResidualsFitter::kPureGaussian, kResidual, kScattering, kSigma, kZpos, LaserClient_cfi::low, MuonResidualsFitter::m_minHits, SiStripPI::mean, MuonResidualsPositionFitter_FCN(), N, funct::pow(), MuonResidualsFitter::residuals_begin(), MuonResidualsFitter::residuals_end(), MuonResidualsFitter::residualsModel(), mathSSE::sqrt(), command_line::start, plotscripts::stdev(), AlCaHLTBitMon_QueryRunRegistry::string, makeMuonMisalignmentScenario::sum_x, and makeMuonMisalignmentScenario::sum_xx.

Referenced by trackingPlots.Iteration::modules().

◆ inform()

void MuonResidualsPositionFitter::inform ( TMinuit *  tMinuit)
overrideprotectedvirtual

◆ ndata()

int MuonResidualsPositionFitter::ndata ( )
inlineoverridevirtual

Implements MuonResidualsFitter.

Definition at line 33 of file MuonResidualsPositionFitter.h.

33 { return kNData; }

References kNData.

◆ npar()

int MuonResidualsPositionFitter::npar ( )
inlineoverridevirtual

◆ plot()

double MuonResidualsPositionFitter::plot ( std::string  name,
TFileDirectory dir,
Alignable ali 
)
overridevirtual

Implements MuonResidualsFitter.

Definition at line 147 of file MuonResidualsPositionFitter.cc.

147  {
148  std::stringstream raw_name, narrowed_name, angleerror_name, trackangle_name, trackposition_name;
149  raw_name << name << "_raw";
150  narrowed_name << name << "_narrowed";
151  angleerror_name << name << "_angleerror";
152  trackangle_name << name << "_trackangle";
153  trackposition_name << name << "_trackposition";
154 
155  TH1F *raw_hist =
156  dir->make<TH1F>(raw_name.str().c_str(), (raw_name.str() + std::string(" (mm)")).c_str(), 100, -100., 100.);
157  TH1F *narrowed_hist = dir->make<TH1F>(
158  narrowed_name.str().c_str(), (narrowed_name.str() + std::string(" (mm)")).c_str(), 100, -100., 100.);
159  TProfile *angleerror_hist = dir->make<TProfile>(
160  angleerror_name.str().c_str(), (angleerror_name.str() + std::string(" (mm)")).c_str(), 100, -30., 30.);
161  TProfile *trackangle_hist = dir->make<TProfile>(
162  trackangle_name.str().c_str(), (trackangle_name.str() + std::string(" (mm)")).c_str(), 100, -0.5, 0.5);
163  TProfile *trackposition_hist = dir->make<TProfile>(
164  trackposition_name.str().c_str(), (trackposition_name.str() + std::string(" (mm)")).c_str(), 100, -300., 300.);
165 
166  angleerror_hist->SetAxisRange(-100., 100., "Y");
167  trackangle_hist->SetAxisRange(-10., 10., "Y");
168  trackposition_hist->SetAxisRange(-10., 10., "Y");
169 
170  narrowed_name << "fit";
171  angleerror_name << "fit";
172  trackangle_name << "fit";
173  trackposition_name << "fit";
174 
175  double scale_factor = double(numResiduals()) * (100. - -100.) / 100; // (max - min)/nbins
176 
177  TF1 *narrowed_fit = nullptr;
178  if (residualsModel() == kPureGaussian) {
179  narrowed_fit = new TF1(narrowed_name.str().c_str(), MuonResidualsFitter_pureGaussian_TF1, -100., 100., 3);
180  narrowed_fit->SetParameters(scale_factor, value(kPosition) * 10., value(kSigma) * 10.);
181  narrowed_fit->Write();
182  } else if (residualsModel() == kPowerLawTails) {
183  narrowed_fit = new TF1(narrowed_name.str().c_str(), MuonResidualsFitter_powerLawTails_TF1, -100., 100., 4);
184  narrowed_fit->SetParameters(scale_factor, value(kPosition) * 10., value(kSigma) * 10., value(kGamma) * 10.);
185  narrowed_fit->Write();
186  } else if (residualsModel() == kROOTVoigt) {
187  narrowed_fit = new TF1(narrowed_name.str().c_str(), MuonResidualsFitter_ROOTVoigt_TF1, -100., 100., 4);
188  narrowed_fit->SetParameters(scale_factor, value(kPosition) * 10., value(kSigma) * 10., value(kGamma) * 10.);
189  narrowed_fit->Write();
190  } else if (residualsModel() == kGaussPowerTails) {
191  narrowed_fit = new TF1(narrowed_name.str().c_str(), MuonResidualsFitter_GaussPowerTails_TF1, -100., 100., 3);
192  narrowed_fit->SetParameters(scale_factor, value(kPosition) * 10., value(kSigma) * 10.);
193  narrowed_fit->Write();
194  }
195 
196  TF1 *angleerror_fit = new TF1(angleerror_name.str().c_str(), "[0]+x*[1]", -30., 30.);
197  angleerror_fit->SetParameters(value(kPosition) * 10., value(kScattering) * 10. / 1000.);
198  angleerror_fit->Write();
199 
200  TF1 *trackangle_fit = new TF1(trackangle_name.str().c_str(), "[0]+x*[1]", -0.5, 0.5);
201  trackangle_fit->SetParameters(value(kPosition) * 10., value(kZpos) * 10.);
202  trackangle_fit->Write();
203 
204  TF1 *trackposition_fit = new TF1(trackposition_name.str().c_str(), "[0]+x*[1]", -300., 300.);
205  trackposition_fit->SetParameters(value(kPosition) * 10., value(kPhiz) * 10.);
206  trackposition_fit->Write();
207 
208  for (std::vector<double *>::const_iterator resiter = residuals_begin(); resiter != residuals_end(); ++resiter) {
209  const double raw_residual = (*resiter)[kResidual];
210  const double angleerror = (*resiter)[kAngleError];
211  const double trackangle = (*resiter)[kTrackAngle];
212  const double trackposition = (*resiter)[kTrackPosition];
213 
214  double angleerror_correction = value(kScattering) * angleerror;
215  double trackangle_correction = value(kZpos) * trackangle;
216  double trackposition_correction = value(kPhiz) * trackposition;
217 
218  double corrected_residual = raw_residual - angleerror_correction - trackangle_correction - trackposition_correction;
219 
220  raw_hist->Fill(raw_residual * 10.);
221  narrowed_hist->Fill(corrected_residual * 10.);
222 
223  angleerror_hist->Fill(angleerror * 1000., (raw_residual - trackangle_correction - trackposition_correction) * 10.);
224  trackangle_hist->Fill(trackangle, (raw_residual - angleerror_correction - trackposition_correction) * 10.);
225  trackposition_hist->Fill(trackposition, (raw_residual - angleerror_correction - trackangle_correction) * 10.);
226  }
227 
228  return 0.;
229 }

References DeadROC_duringRun::dir, kAngleError, kGamma, MuonResidualsFitter::kGaussPowerTails, kPhiz, kPosition, MuonResidualsFitter::kPowerLawTails, MuonResidualsFitter::kPureGaussian, kResidual, MuonResidualsFitter::kROOTVoigt, kScattering, kSigma, kTrackAngle, kTrackPosition, kZpos, MuonResidualsFitter_GaussPowerTails_TF1(), MuonResidualsFitter_powerLawTails_TF1(), MuonResidualsFitter_pureGaussian_TF1(), MuonResidualsFitter_ROOTVoigt_TF1(), Skims_PA_cff::name, MuonResidualsFitter::numResiduals(), MuonResidualsFitter::residuals_begin(), MuonResidualsFitter::residuals_end(), MuonResidualsFitter::residualsModel(), AlCaHLTBitMon_QueryRunRegistry::string, and MuonResidualsFitter::value().

◆ sumofweights()

double MuonResidualsPositionFitter::sumofweights ( )
inlineoverridevirtual

Implements MuonResidualsFitter.

Definition at line 36 of file MuonResidualsPositionFitter.h.

36 { return numResiduals(); }

References MuonResidualsFitter::numResiduals().

◆ type()

int MuonResidualsPositionFitter::type ( ) const
inlineoverridevirtual
MuonResidualsFitter_pureGaussian_TF1
Double_t MuonResidualsFitter_pureGaussian_TF1(Double_t *xvec, Double_t *par)
Definition: MuonResidualsFitter.cc:38
start
Definition: start.py:1
MuonResidualsPositionFitter::kGamma
Definition: MuonResidualsPositionFitter.h:14
MuonResidualsPositionFitter::kScattering
Definition: MuonResidualsPositionFitter.h:14
SiStripPI::mean
Definition: SiStripPayloadInspectorHelper.h:169
MuonResidualsFitter::MuonResidualsFitter
MuonResidualsFitter(int residualsModel, int minHits, int useResiduals, bool weightAlignment=true)
Definition: MuonResidualsFitter.cc:151
MuonResidualsPositionFitter::kNData
Definition: MuonResidualsPositionFitter.h:16
MuonResidualsPositionFitter::kAngleError
Definition: MuonResidualsPositionFitter.h:16
MuonResidualsFitter::value
double value(int parNum)
Definition: MuonResidualsFitter.h:145
align_cfg.useResiduals
useResiduals
Definition: align_cfg.py:36
step
step
Definition: StallMonitor.cc:94
MuonResidualsPositionFitter::kTrackPosition
Definition: MuonResidualsPositionFitter.h:16
MuonResidualsPositionFitter::kPhiz
Definition: MuonResidualsPositionFitter.h:14
MuonResidualsFitter::dofit
bool dofit(void(*fcn)(int &, double *, double &, double *, int), std::vector< int > &parNum, std::vector< std::string > &parName, std::vector< double > &start, std::vector< double > &step, std::vector< double > &low, std::vector< double > &high)
Definition: MuonResidualsFitter.cc:272
MuonResidualsFitter::numResiduals
long numResiduals() const
Definition: MuonResidualsFitter.h:125
cms::cuda::assert
assert(be >=bs)
MuonResidualsPositionFitter::kResidual
Definition: MuonResidualsPositionFitter.h:16
MuonResidualsFitter::kGaussPowerTails
Definition: MuonResidualsFitter.h:86
MuonResidualsPositionFitter_TMinuit
static TMinuit * MuonResidualsPositionFitter_TMinuit
Definition: MuonResidualsPositionFitter.cc:3
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:19
MuonResidualsFitter::kPowerLawTails
Definition: MuonResidualsFitter.h:86
MuonResidualsPositionFitter::kZpos
Definition: MuonResidualsPositionFitter.h:14
MuonResidualsFitter_powerLawTails_TF1
Double_t MuonResidualsFitter_powerLawTails_TF1(Double_t *xvec, Double_t *par)
Definition: MuonResidualsFitter.cc:116
N
#define N
Definition: blowfish.cc:9
MuonResidualsPositionFitter::kPosition
Definition: MuonResidualsPositionFitter.h:14
MuonResidualsFitter_ROOTVoigt_TF1
Double_t MuonResidualsFitter_ROOTVoigt_TF1(Double_t *xvec, Double_t *par)
Definition: MuonResidualsFitter.cc:124
MuonResidualsFitter::kPureGaussian
Definition: MuonResidualsFitter.h:86
MuonResidualsFitter::residuals_end
std::vector< double * >::const_iterator residuals_end() const
Definition: MuonResidualsFitter.h:191
MuonResidualsFitter::kPositionFitter
Definition: MuonResidualsFitter.h:88
MuonResidualsFitter::m_minHits
int m_minHits
Definition: MuonResidualsFitter.h:225
MuonResidualsPositionFitter_FCN
void MuonResidualsPositionFitter_FCN(int &npar, double *gin, double &fval, double *par, int iflag)
Definition: MuonResidualsPositionFitter.cc:7
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
MuonResidualsPositionFitter::kSigma
Definition: MuonResidualsPositionFitter.h:14
MuonResidualsFitter::kROOTVoigt
Definition: MuonResidualsFitter.h:86
plotscripts.stdev
def stdev(xlist)
Definition: plotscripts.py:69
LaserClient_cfi.high
high
Definition: LaserClient_cfi.py:50
align_cfg.weightAlignment
tuple weightAlignment
Definition: align_cfg.py:30
MuonResidualsPositionFitter::kNPar
Definition: MuonResidualsPositionFitter.h:14
Skims_PA_cff.name
name
Definition: Skims_PA_cff.py:17
MuonResidualsFitter::initialize_table
void initialize_table()
Definition: MuonResidualsFitter.cc:205
MuonResidualsFitter_GaussPowerTails_TF1
Double_t MuonResidualsFitter_GaussPowerTails_TF1(Double_t *xvec, Double_t *par)
Definition: MuonResidualsFitter.cc:142
funct::pow
Power< A, B >::type pow(const A &a, const B &b)
Definition: Power.h:29
MuonResidualsPositionFitter::kTrackAngle
Definition: MuonResidualsPositionFitter.h:16
command_line.start
start
Definition: command_line.py:167
makeMuonMisalignmentScenario.sum_x
int sum_x
More diagnostics.
Definition: makeMuonMisalignmentScenario.py:176
MuonResidualsFitter::residualsModel
int residualsModel() const
Definition: MuonResidualsFitter.h:124
makeMuonMisalignmentScenario.sum_xx
int sum_xx
Definition: makeMuonMisalignmentScenario.py:183
MuonResidualsFitter::residuals_begin
std::vector< double * >::const_iterator residuals_begin() const
Definition: MuonResidualsFitter.h:190
LaserClient_cfi.low
low
Definition: LaserClient_cfi.py:52
reco_skim_cfg_mod.minHits
minHits
Definition: reco_skim_cfg_mod.py:177
DeadROC_duringRun.dir
dir
Definition: DeadROC_duringRun.py:23