CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Groups Pages
List of all members | Public Types | Public Member Functions | Protected Member Functions
MuonResidualsBfieldAngleFitter Class Reference

#include <MuonResidualsBfieldAngleFitter.h>

Inheritance diagram for MuonResidualsBfieldAngleFitter:
MuonResidualsFitter

Public Types

enum  {
  kAngle = 0, kBfrompt, kBfrompz, kdEdx,
  kSigma, kGamma, kNPar
}
 
enum  { kResidual = 0, kQoverPt, kQoverPz, 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
 
 MuonResidualsBfieldAngleFitter (int residualsModel, int minHitsPerRegion, 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:26
Revision:
1.5
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 MuonResidualsBfieldAngleFitter.h.

Member Enumeration Documentation

anonymous enum
Enumerator
kAngle 
kBfrompt 
kBfrompz 
kdEdx 
kSigma 
kGamma 
kNPar 

Definition at line 14 of file MuonResidualsBfieldAngleFitter.h.

14 { kAngle = 0, kBfrompt, kBfrompz, kdEdx, kSigma, kGamma, kNPar };
anonymous enum
Enumerator
kResidual 
kQoverPt 
kQoverPz 
kNData 

Definition at line 16 of file MuonResidualsBfieldAngleFitter.h.

16 { kResidual = 0, kQoverPt, kQoverPz, kNData };

Constructor & Destructor Documentation

MuonResidualsBfieldAngleFitter::MuonResidualsBfieldAngleFitter ( int  residualsModel,
int  minHitsPerRegion,
int  useResiduals,
bool  weightAlignment = true 
)
inline

Definition at line 18 of file MuonResidualsBfieldAngleFitter.h.

19  : MuonResidualsFitter(residualsModel, minHitsPerRegion, useResiduals, weightAlignment) {}
align_cfg.weightAlignment
tuple weightAlignment
Definition: align_cfg.py:30
MuonResidualsFitter::MuonResidualsFitter
MuonResidualsFitter(int residualsModel, int minHits, int useResiduals, bool weightAlignment=true)
Definition: MuonResidualsFitter.cc:151
align_cfg.useResiduals
list useResiduals
Definition: align_cfg.py:36

Member Function Documentation

bool MuonResidualsBfieldAngleFitter::fit ( Alignable *  ali)
overridevirtual

Implements MuonResidualsFitter.

Definition at line 42 of file MuonResidualsBfieldAngleFitter.cc.

References MuonResidualsFitter::dofit(), MuonResidualsFitter::initialize_table(), kAngle, kBfrompt, kBfrompz, kdEdx, kGamma, MuonResidualsFitter::kGaussPowerTails, MuonResidualsFitter::kPureGaussian, kResidual, kSigma, MuonResidualsFitter::m_minHits, SiStripPI::mean, MuonResidualsBfieldAngleFitter_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().

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 qoverpt = (*resiter)[kQoverPt];
52 
53  if (fabs(residual) < 0.1) { // truncate at 100 mrad
54  sum_x += residual;
55  sum_xx += residual * residual;
56  N++;
57  }
58  }
59 
60  if (N < m_minHits)
61  return false;
62 
63  // truncated mean and stdev to seed the fit
64  double mean = sum_x / double(N);
65  double stdev = sqrt(sum_xx / double(N) - pow(sum_x / double(N), 2));
66 
67  // refine the standard deviation calculation
68  sum_x = 0.;
69  sum_xx = 0.;
70  N = 0;
71  for (std::vector<double *>::const_iterator resiter = residuals_begin(); resiter != residuals_end(); ++resiter) {
72  const double residual = (*resiter)[kResidual];
73  if (mean - 1.5 * stdev < residual && residual < mean + 1.5 * stdev) {
74  sum_x += residual;
75  sum_xx += residual * residual;
76  N++;
77  }
78  }
79  mean = sum_x / double(N);
80  stdev = sqrt(sum_xx / double(N) - pow(sum_x / double(N), 2));
81 
82  sum_x = 0.;
83  sum_xx = 0.;
84  N = 0;
85  for (std::vector<double *>::const_iterator resiter = residuals_begin(); resiter != residuals_end(); ++resiter) {
86  const double residual = (*resiter)[kResidual];
87  if (mean - 1.5 * stdev < residual && residual < mean + 1.5 * stdev) {
88  sum_x += residual;
89  sum_xx += residual * residual;
90  N++;
91  }
92  }
93  mean = sum_x / double(N);
94  stdev = sqrt(sum_xx / double(N) - pow(sum_x / double(N), 2));
95 
96  std::vector<int> parNum;
97  std::vector<std::string> parName;
98  std::vector<double> start;
99  std::vector<double> step;
100  std::vector<double> low;
101  std::vector<double> high;
102 
103  parNum.push_back(kAngle);
104  parName.push_back(std::string("angle"));
105  start.push_back(mean);
106  step.push_back(0.1);
107  low.push_back(0.);
108  high.push_back(0.);
109  parNum.push_back(kBfrompt);
110  parName.push_back(std::string("bfrompt"));
111  start.push_back(0.);
112  step.push_back(0.1 * stdev / 0.05);
113  low.push_back(0.);
114  high.push_back(0.);
115  parNum.push_back(kBfrompz);
116  parName.push_back(std::string("bfrompz"));
117  start.push_back(0.);
118  step.push_back(0.1 * stdev / 0.05);
119  low.push_back(0.);
120  high.push_back(0.);
121  parNum.push_back(kdEdx);
122  parName.push_back(std::string("dEdx"));
123  start.push_back(0.);
124  step.push_back(0.1 * stdev / 0.05);
125  low.push_back(0.);
126  high.push_back(0.);
127  parNum.push_back(kSigma);
128  parName.push_back(std::string("sigma"));
129  start.push_back(stdev);
130  step.push_back(0.1 * stdev);
131  low.push_back(0.);
132  high.push_back(0.);
133  if (residualsModel() != kPureGaussian && residualsModel() != kGaussPowerTails) {
134  parNum.push_back(kGamma);
135  parName.push_back(std::string("gamma"));
136  start.push_back(stdev);
137  step.push_back(0.1 * stdev);
138  low.push_back(0.);
139  high.push_back(0.);
140  }
141 
142  return dofit(&MuonResidualsBfieldAngleFitter_FCN, parNum, parName, start, step, low, high);
143 }
MuonResidualsFitter::residuals_begin
std::vector< double * >::const_iterator residuals_begin() const
Definition: MuonResidualsFitter.h:190
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
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
command_line.start
start
Definition: command_line.py:165
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:19
N
#define N
Definition: blowfish.cc:9
MuonResidualsFitter::residuals_end
std::vector< double * >::const_iterator residuals_end() const
Definition: MuonResidualsFitter.h:191
step
step
Definition: StallMonitor.cc:94
MuonResidualsBfieldAngleFitter_FCN
void MuonResidualsBfieldAngleFitter_FCN(int &npar, double *gin, double &fval, double *par, int iflag)
Definition: MuonResidualsBfieldAngleFitter.cc:7
funct::pow
Power< A, B >::type pow(const A &a, const B &b)
Definition: Power.h:29
void MuonResidualsBfieldAngleFitter::inform ( TMinuit *  tMinuit)
overrideprotectedvirtual

Implements MuonResidualsFitter.

Definition at line 5 of file MuonResidualsBfieldAngleFitter.cc.

References MuonResidualsBfieldAngleFitter_TMinuit.

5 { MuonResidualsBfieldAngleFitter_TMinuit = tMinuit; }
MuonResidualsBfieldAngleFitter_TMinuit
static TMinuit * MuonResidualsBfieldAngleFitter_TMinuit
Definition: MuonResidualsBfieldAngleFitter.cc:3
int MuonResidualsBfieldAngleFitter::ndata ( )
inlineoverridevirtual

Implements MuonResidualsFitter.

Definition at line 33 of file MuonResidualsBfieldAngleFitter.h.

References kNData.

33 { return kNData; }
int MuonResidualsBfieldAngleFitter::npar ( )
inlineoverridevirtual

Implements MuonResidualsFitter.

Definition at line 23 of file MuonResidualsBfieldAngleFitter.h.

References cms::cuda::assert(), MuonResidualsFitter::kGaussPowerTails, kNPar, MuonResidualsFitter::kPowerLawTails, MuonResidualsFitter::kPureGaussian, MuonResidualsFitter::kROOTVoigt, and MuonResidualsFitter::residualsModel().

23  {
24  if (residualsModel() == kPureGaussian || residualsModel() == kGaussPowerTails)
25  return kNPar - 1;
26  else if (residualsModel() == kPowerLawTails)
27  return kNPar;
28  else if (residualsModel() == kROOTVoigt)
29  return kNPar;
30  else
31  assert(false);
32  }
cms::cuda::assert
assert(be >=bs)
double MuonResidualsBfieldAngleFitter::plot ( std::string  name,
TFileDirectory *  dir,
Alignable *  ali 
)
overridevirtual

Implements MuonResidualsFitter.

Definition at line 145 of file MuonResidualsBfieldAngleFitter.cc.

References kAngle, kBfrompt, kBfrompz, kdEdx, kGamma, MuonResidualsFitter::kGaussPowerTails, MuonResidualsFitter::kPowerLawTails, MuonResidualsFitter::kPureGaussian, kQoverPt, kQoverPz, kResidual, MuonResidualsFitter::kROOTVoigt, kSigma, TFileDirectory::make(), MuonResidualsFitter_GaussPowerTails_TF1(), MuonResidualsFitter_powerLawTails_TF1(), MuonResidualsFitter_pureGaussian_TF1(), MuonResidualsFitter_ROOTVoigt_TF1(), MuonResidualsFitter::numResiduals(), MuonResidualsFitter::residuals_begin(), MuonResidualsFitter::residuals_end(), MuonResidualsFitter::residualsModel(), AlCaHLTBitMon_QueryRunRegistry::string, and MuonResidualsFitter::value().

145  {
146  std::stringstream raw_name, narrowed_name, qoverpt_name, qoverpz_name, psquared_name;
147  raw_name << name << "_raw";
148  narrowed_name << name << "_narrowed";
149  qoverpt_name << name << "_qoverpt";
150  qoverpz_name << name << "_qoverpz";
151  psquared_name << name << "_psquared";
152 
153  TH1F *raw_hist =
154  dir->make<TH1F>(raw_name.str().c_str(), (raw_name.str() + std::string(" (mrad)")).c_str(), 100, -100., 100.);
155  TH1F *narrowed_hist = dir->make<TH1F>(
156  narrowed_name.str().c_str(), (narrowed_name.str() + std::string(" (mrad)")).c_str(), 100, -100., 100.);
157  TProfile *qoverpt_hist = dir->make<TProfile>(
158  qoverpt_name.str().c_str(), (qoverpt_name.str() + std::string(" (mrad)")).c_str(), 100, -0.05, 0.05);
159  TProfile *qoverpz_hist = dir->make<TProfile>(
160  qoverpz_name.str().c_str(), (qoverpz_name.str() + std::string(" (mrad)")).c_str(), 100, -0.05, 0.05);
161  TProfile *psquared_hist = dir->make<TProfile>(
162  psquared_name.str().c_str(), (psquared_name.str() + std::string(" (mrad)")).c_str(), 100, -0.05, 0.05);
163 
164  narrowed_name << "fit";
165  qoverpt_name << "fit";
166  qoverpz_name << "fit";
167  psquared_name << "fit";
168 
169  double scale_factor = double(numResiduals()) * (100. - -100.) / 100; // (max - min)/nbins
170 
171  TF1 *narrowed_fit = nullptr;
172  if (residualsModel() == kPureGaussian) {
173  narrowed_fit = new TF1(narrowed_name.str().c_str(), MuonResidualsFitter_pureGaussian_TF1, -100., 100., 3);
174  narrowed_fit->SetParameters(scale_factor, value(kAngle) * 1000., value(kSigma) * 1000.);
175  narrowed_fit->Write();
176  } else if (residualsModel() == kPowerLawTails) {
177  narrowed_fit = new TF1(narrowed_name.str().c_str(), MuonResidualsFitter_powerLawTails_TF1, -100., 100., 4);
178  narrowed_fit->SetParameters(scale_factor, value(kAngle) * 1000., value(kSigma) * 1000., value(kGamma) * 1000.);
179  narrowed_fit->Write();
180  } else if (residualsModel() == kROOTVoigt) {
181  narrowed_fit = new TF1(narrowed_name.str().c_str(), MuonResidualsFitter_ROOTVoigt_TF1, -100., 100., 4);
182  narrowed_fit->SetParameters(scale_factor, value(kAngle) * 1000., value(kSigma) * 1000., value(kGamma) * 1000.);
183  narrowed_fit->Write();
184  } else if (residualsModel() == kGaussPowerTails) {
185  narrowed_fit = new TF1(narrowed_name.str().c_str(), MuonResidualsFitter_GaussPowerTails_TF1, -100., 100., 3);
186  narrowed_fit->SetParameters(scale_factor, value(kAngle) * 1000., value(kSigma) * 1000.);
187  narrowed_fit->Write();
188  }
189 
190  TF1 *qoverpt_fit = new TF1(qoverpt_name.str().c_str(), "[0]+x*[1]", -0.05, 0.05);
191  qoverpt_fit->SetParameters(value(kAngle) * 1000., value(kBfrompt) * 1000.);
192  qoverpt_fit->Write();
193 
194  TF1 *qoverpz_fit = new TF1(qoverpz_name.str().c_str(), "[0]+x*[1]", -0.05, 0.05);
195  qoverpz_fit->SetParameters(value(kAngle) * 1000., value(kBfrompz) * 1000.);
196  qoverpz_fit->Write();
197 
198  TF1 *psquared_fit = new TF1(psquared_name.str().c_str(), "[0]+[1]*x**2", -0.05, 0.05);
199  psquared_fit->SetParameters(value(kAngle) * 1000., value(kdEdx) * 1000.);
200  psquared_fit->Write();
201 
202  for (std::vector<double *>::const_iterator resiter = residuals_begin(); resiter != residuals_end(); ++resiter) {
203  const double raw_residual = (*resiter)[kResidual];
204  const double qoverpt = (*resiter)[kQoverPt];
205  const double qoverpz = (*resiter)[kQoverPz];
206  const double psquared = (1. / qoverpt / qoverpt + 1. / qoverpz / qoverpz) * (qoverpt > 0. ? 1. : -1.);
207 
208  double qoverpt_correction = value(kBfrompt) * qoverpt;
209  double qoverpz_correction = value(kBfrompz) * qoverpz;
210  double dEdx_correction = value(kdEdx) * psquared;
211  double corrected_residual = raw_residual - qoverpt_correction - qoverpz_correction - dEdx_correction;
212 
213  raw_hist->Fill(raw_residual * 1000.);
214  narrowed_hist->Fill(corrected_residual * 1000.);
215 
216  qoverpt_hist->Fill(qoverpt, (raw_residual - qoverpz_correction - dEdx_correction) * 1000.);
217  qoverpz_hist->Fill(qoverpz, (raw_residual - qoverpt_correction - dEdx_correction) * 1000.);
218  psquared_hist->Fill(psquared, (raw_residual - qoverpt_correction - qoverpz_correction) * 1000.);
219  }
220 
221  return 0.;
222 }
MuonResidualsFitter::residuals_begin
std::vector< double * >::const_iterator residuals_begin() const
Definition: MuonResidualsFitter.h:190
MuonResidualsFitter_powerLawTails_TF1
Double_t MuonResidualsFitter_powerLawTails_TF1(Double_t *xvec, Double_t *par)
Definition: MuonResidualsFitter.cc:116
MuonResidualsFitter::value
double value(int parNum)
Definition: MuonResidualsFitter.h:145
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
MuonResidualsFitter_GaussPowerTails_TF1
Double_t MuonResidualsFitter_GaussPowerTails_TF1(Double_t *xvec, Double_t *par)
Definition: MuonResidualsFitter.cc:142
TFileDirectory::make
T * make(const Args &...args) const
make new ROOT object
Definition: TFileDirectory.h:53
MuonResidualsFitter_pureGaussian_TF1
Double_t MuonResidualsFitter_pureGaussian_TF1(Double_t *xvec, Double_t *par)
Definition: MuonResidualsFitter.cc:38
MuonResidualsFitter::residuals_end
std::vector< double * >::const_iterator residuals_end() const
Definition: MuonResidualsFitter.h:191
MuonResidualsFitter_ROOTVoigt_TF1
Double_t MuonResidualsFitter_ROOTVoigt_TF1(Double_t *xvec, Double_t *par)
Definition: MuonResidualsFitter.cc:124
double MuonResidualsBfieldAngleFitter::sumofweights ( )
inlineoverridevirtual

Implements MuonResidualsFitter.

Definition at line 36 of file MuonResidualsBfieldAngleFitter.h.

References MuonResidualsFitter::numResiduals().

36 { return numResiduals(); }
int MuonResidualsBfieldAngleFitter::type ( ) const
inlineoverridevirtual

Implements MuonResidualsFitter.

Definition at line 21 of file MuonResidualsBfieldAngleFitter.h.

References MuonResidualsFitter::kAngleBfieldFitter.

21 { return MuonResidualsFitter::kAngleBfieldFitter; }
Generated for CMSSW Reference Manual by doxygen 1.8.5