CMS 3D CMS Logo
List of all members | Classes | Public Types | Public Member Functions | Static Public Member Functions | Private Attributes
LA_Filler_Fitter Class Reference

#include <LA_Filler_Fitter.h>

Classes

struct  EnsembleSummary
 
struct  Result
 

Public Types

enum  Method {
  WIDTH = 1 << 0, FIRST_METHOD = 1 << 0, PROB1 = 1 << 1, AVGV2 = 1 << 2,
  AVGV3 = 1 << 3, RMSV2 = 1 << 4, RMSV3 = 1 << 5, LAST_METHOD = 1 << 5
}
 

Public Member Functions

poly< std::string > allAndOne (const unsigned width) const
 
void fill (TTree *, Book &) const
 
void fill_one_cluster (Book &, const poly< std::string > &, const unsigned, const float, const float, const float, const float) const
 
poly< std::string > granularity (const SiStripDetId, const float, const Long64_t, const float, const unsigned) const
 
 LA_Filler_Fitter (int methods, bool layer, bool module, float localybin, unsigned stripbin, unsigned max, const TrackerTopology *tTopo)
 
 LA_Filler_Fitter (int methods, int M, int N, double low, double up, unsigned max, const TrackerTopology *tTopo)
 
std::string layerLabel (const SiStripDetId) const
 
void summarize_ensembles (Book &) const
 
poly< std::string > varWidth (const unsigned width) const
 

Static Public Member Functions

static std::map< std::string, std::vector< Result > > ensemble_results (const Book &, const Method)
 
static std::map< std::string, std::vector< EnsembleSummary > > ensemble_summary (const Book &)
 
static unsigned find_rebin (const TH1 *const)
 
static void fit (Book &book)
 
static void fit_width_profile (Book &)
 
static unsigned layer_index (bool TIB, bool stereo, unsigned layer)
 
static std::map< std::string, Result > layer_results (const Book &, const Method)
 
static void make_and_fit_symmchi2 (Book &)
 
static std::string method (Method m, bool fit=true)
 
static std::map< uint32_t, Result > module_results (const Book &, const Method)
 
static std::string moduleLabel (const SiStripDetId)
 
static std::pair< std::pair< float, float >, std::pair< float, float > > offset_slope (const std::vector< EnsembleSummary > &)
 
static float pull (const std::vector< EnsembleSummary > &)
 
static Result result (Method, const std::string name, const Book &)
 
static TH1 * rms_profile (const std::string, const TProfile *const)
 
static std::string subdetLabel (const SiStripDetId)
 
static TH1 * subset_probability (const std::string name, const TH1 *const, const TH1 *const)
 

Private Attributes

const bool byLayer_
 
const bool byModule_
 
const int ensembleBins_
 
const double ensembleLow_
 
const int ensembleSize_
 
const double ensembleUp_
 
const float localYbin_
 
const Long64_t maxEvents_
 
const int methods_
 
const unsigned stripsPerBin_
 
const TrackerTopology * tTopo_
 

Detailed Description

Definition at line 14 of file LA_Filler_Fitter.h.

Member Enumeration Documentation

◆ Method

enum LA_Filler_Fitter::Method
Enumerator
WIDTH 
FIRST_METHOD 
PROB1 
AVGV2 
AVGV3 
RMSV2 
RMSV3 
LAST_METHOD 

Definition at line 50 of file LA_Filler_Fitter.h.

50  {
51  WIDTH = 1 << 0,
52  FIRST_METHOD = 1 << 0,
53  PROB1 = 1 << 1,
54  AVGV2 = 1 << 2,
55  AVGV3 = 1 << 3,
56  RMSV2 = 1 << 4,
57  RMSV3 = 1 << 5,
58  LAST_METHOD = 1 << 5
59  };

Constructor & Destructor Documentation

◆ LA_Filler_Fitter() [1/2]

LA_Filler_Fitter::LA_Filler_Fitter ( int  methods,
int  M,
int  N,
double  low,
double  up,
unsigned  max,
const TrackerTopology *  tTopo 
)
inline

Definition at line 17 of file LA_Filler_Fitter.h.

18  : ensembleSize_(M),
19  ensembleBins_(N),
20  ensembleLow_(low),
21  ensembleUp_(up),
22  byLayer_(true),
23  byModule_(false),
24  localYbin_(0),
25  stripsPerBin_(0),
26  maxEvents_(max),
27  methods_(methods),
28  tTopo_(tTopo) {}

◆ LA_Filler_Fitter() [2/2]

LA_Filler_Fitter::LA_Filler_Fitter ( int  methods,
bool  layer,
bool  module,
float  localybin,
unsigned  stripbin,
unsigned  max,
const TrackerTopology *  tTopo 
)
inline

Definition at line 31 of file LA_Filler_Fitter.h.

38  : ensembleSize_(0),
39  ensembleBins_(0),
40  ensembleLow_(0),
41  ensembleUp_(0),
42  byLayer_(layer),
43  byModule_(module),
44  localYbin_(localybin),
45  stripsPerBin_(stripbin),
46  maxEvents_(max),
47  methods_(methods),
48  tTopo_(tTopo) {}

Member Function Documentation

◆ allAndOne()

poly< std::string > LA_Filler_Fitter::allAndOne ( const unsigned  width) const

Definition at line 78 of file LA_Filler.cc.

78  {
79  poly<std::string> a1("_all");
80  if (width == 1)
81  a1 *= "_w1";
82  return a1;
83 }

References ApeEstimator_cff::width.

Referenced by fill_one_cluster().

◆ ensemble_results()

std::map< std::string, std::vector< LA_Filler_Fitter::Result > > LA_Filler_Fitter::ensemble_results ( const Book &  book,
const Method  m 
)
static

Definition at line 80 of file LA_Results.cc.

81  {
82  std::map<std::string, std::vector<Result> > results;
83  for (Book::const_iterator it = book.begin(".*_sample.*" + method(m)); it != book.end(); ++it) {
84  const std::string name = boost::regex_replace(it->first, boost::regex("sample\\d*_"), "");
85  results[name].push_back(result(m, it->first, book));
86  }
87  return results;
88 }

References Book::begin(), Book::end(), visualization-live-secondInstance_cfg::m, method(), Skims_PA_cff::name, result(), bookConverter::results, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by summarize_ensembles().

◆ ensemble_summary()

std::map< std::string, std::vector< LA_Filler_Fitter::EnsembleSummary > > LA_Filler_Fitter::ensemble_summary ( const Book &  book)
static

Definition at line 116 of file LA_Results.cc.

117  {
118  std::map<std::string, std::vector<EnsembleSummary> > summary;
119  for (Book::const_iterator it = book.begin(".*_ensembleReco"); it != book.end(); ++it) {
120  const std::string base = boost::erase_all_copy(it->first, "_ensembleReco");
121 
122  const TH1* const reco = it->second;
123  const TH1* const measure = book[base + "_measure"];
124  const TH1* const merr = book[base + "_merr"];
125  const TH1* const pull = book[base + "_pull"];
126 
127  EnsembleSummary s;
128  s.samples = reco->GetEntries();
129  s.truth = reco->GetMean();
130  s.meanMeasured = std::make_pair<float, float>(measure->GetFunction("gaus")->GetParameter(1),
131  measure->GetFunction("gaus")->GetParError(1));
132  s.sigmaMeasured = std::make_pair<float, float>(measure->GetFunction("gaus")->GetParameter(2),
133  measure->GetFunction("gaus")->GetParError(2));
134  s.meanUncertainty = std::make_pair<float, float>(merr->GetFunction("gaus")->GetParameter(1),
135  merr->GetFunction("gaus")->GetParError(1));
136  s.pull = std::make_pair<float, float>(pull->GetFunction("gaus")->GetParameter(2),
137  pull->GetFunction("gaus")->GetParError(2));
138 
139  const std::string name = boost::regex_replace(base, boost::regex("ensembleBin\\d*_"), "");
140  summary[name].push_back(s);
141  }
142  return summary;
143 }

References newFWLiteAna::base, Book::begin(), Book::end(), Skims_PA_cff::name, pull(), alignCSCRings::s, AlCaHLTBitMon_QueryRunRegistry::string, and edmLumisInFiles::summary.

Referenced by sistrip::EnsembleCalibrationLA::write_ensembles_text().

◆ fill()

void LA_Filler_Fitter::fill ( TTree *  tree,
Book &  book 
) const

Definition at line 8 of file LA_Filler.cc.

8  {
9  TTREE_FOREACH_ENTRY(tree) {
10  TFE_MAX(maxEvents_);
11  TFE_PRINTSTATUS;
12  std::vector<unsigned> PLEAF(tsostrackmulti, tree);
13  std::vector<unsigned> PLEAF(clusterdetid, tree);
14  std::vector<unsigned> PLEAF(clusterwidth, tree);
15  std::vector<float> PLEAF(clustervariance, tree);
16  std::vector<float> PLEAF(tsosdriftx, tree);
17  std::vector<float> PLEAF(tsosdriftz, tree);
18  std::vector<float> PLEAF(tsoslocaltheta, tree);
19  std::vector<float> PLEAF(tsoslocalphi, tree);
20  std::vector<float> PLEAF(tsosglobalZofunitlocalY, tree);
21 
22  const unsigned N(clusterdetid.size());
23  std::vector<float> BdotY(N, 0);
24  if (!ensembleBins_) {
25  std::vector<float> PLEAF(tsosBdotY, tree);
26  swap(BdotY, tsosBdotY);
27  }
28  std::vector<float> localy(N, 0);
29  if (localYbin_) {
30  std::vector<float> PLEAF(tsoslocaly, tree);
31  swap(localy, tsoslocaly);
32  }
33  std::vector<unsigned> seedstrip(N, 0);
34  if (stripsPerBin_) {
35  std::vector<unsigned> PLEAF(clusterseedstrip, tree);
36  swap(seedstrip, clusterseedstrip);
37  }
38 
39  for (unsigned i = 0; i < N; i++) {
40  const SiStripDetId detid(clusterdetid[i]);
41  if (tsostrackmulti[i] != 1 ||
42  (detid.subDetector() != SiStripDetId::TIB && detid.subDetector() != SiStripDetId::TOB))
43  continue;
44 
45  const int sign = tsosglobalZofunitlocalY[i] < 0 ? -1 : 1;
46  const float tthetaL = sign * tsosdriftx[i] / tsosdriftz[i];
47  const float tthetaT = sign * tan(tsoslocaltheta[i]) * cos(tsoslocalphi[i]);
48 
49  fill_one_cluster(book,
50  granularity(detid, tthetaL, TFE_index, localy[i], seedstrip[i] % 128),
51  clusterwidth[i],
52  clustervariance[i],
53  tthetaL,
54  tthetaT,
55  fabs(BdotY[i]));
56  }
57  }
58 }

References funct::cos(), ensembleBins_, fill_one_cluster(), granularity(), mps_fire::i, localYbin_, maxEvents_, N, Validation_hcalonly_cfi::sign, stripsPerBin_, SiStripDetId::subDetector(), edm::swap(), funct::tan(), SiStripDetId::TIB, and SiStripDetId::TOB.

Referenced by sistrip::EnsembleCalibrationLA::endJob(), and sistrip::MeasureLA::MeasureLA().

◆ fill_one_cluster()

void LA_Filler_Fitter::fill_one_cluster ( Book &  book,
const poly< std::string > &  gran,
const unsigned  width,
const float  variance,
const float  tthetaL,
const float  tthetaT,
const float  BdotY 
) const

Definition at line 60 of file LA_Filler.cc.

66  {
67  book.fill(tthetaL, gran + "_reconstruction", 360, -1.0, 1.0);
68  book.fill(tthetaT - tthetaL, gran + allAndOne(width), 360, -1.0, 1.0);
69  book.fill(tthetaT - tthetaL, variance, gran + varWidth(width), 360, -1.0, 1.0);
70  if (methods_ & WIDTH)
71  book.fill(tthetaT, width, gran + method(WIDTH), 81, -0.6, 0.6);
72  if (!ensembleBins_) {
73  book.fill(BdotY, gran + "_field", 101, 1, 5);
74  book.fill(width, gran + "_width", 10, 0, 10);
75  }
76 }

References allAndOne(), ensembleBins_, Book::fill(), method(), methods_, varWidth(), ApeEstimator_cff::width, and WIDTH.

Referenced by fill().

◆ find_rebin()

unsigned LA_Filler_Fitter::find_rebin ( const TH1 * const  hist)
static

Definition at line 104 of file LA_Fitter.cc.

104  {
105  const double mean = hist->GetMean();
106  const double rms = hist->GetRMS();
107  const int begin = std::min(1, hist->GetXaxis()->FindFixBin(mean - rms));
108  const int end = std::max(hist->GetNbinsX(), hist->GetXaxis()->FindFixBin(mean + rms)) + 1;
109  unsigned current_hole(0), max_hole(0);
110  for (int i = begin; i < end; i++) {
111  if (!hist->GetBinError(i))
112  current_hole++;
113  else if (current_hole) {
114  max_hole = std::max(current_hole, max_hole);
115  current_hole = 0;
116  }
117  }
118  return max_hole + 1;
119 }

References begin, end, compare::hist, mps_fire::i, SiStripPI::max, SiStripPI::mean, min(), and SiStripPI::rms.

Referenced by make_and_fit_symmchi2().

◆ fit()

static void LA_Filler_Fitter::fit ( Book &  book)
inlinestatic

Definition at line 110 of file LA_Filler_Fitter.h.

110  {
111  make_and_fit_symmchi2(book);
112  fit_width_profile(book);
113  }

References fit_width_profile(), and make_and_fit_symmchi2().

Referenced by sistrip::EnsembleCalibrationLA::endJob(), sistrip::MeasureLA::MeasureLA(), and trackingPlots.Iteration::modules().

◆ fit_width_profile()

void LA_Filler_Fitter::fit_width_profile ( Book &  book)
static

Definition at line 8 of file LA_Fitter.cc.

8  {
9  for (Book::iterator it = book.begin(".*" + method(WIDTH)); it != book.end(); ++it) {
10  it->second->SetTitle("Mean Cluster Width;tan#theta_{t}");
11  TH1* const p = it->second;
12  if (p->GetEntries() < 400) {
13  delete p;
14  book[it->first] = nullptr;
15  continue;
16  }
17  p->SetTitle(";tan#theta_{t};");
18  const float min = p->GetMinimum();
19  const float max = p->GetMaximum();
20  float xofmin = p->GetBinCenter(p->GetMinimumBin());
21  if (xofmin > 0.0 || xofmin < -0.15)
22  xofmin = -0.05;
23  const float xofmax = p->GetBinCenter(p->GetMaximumBin());
24 
25  TF1* const fit = new TF1("LA_profile_fit", "[2]*(TMath::Abs(x-[0]))+[1]", -1, 1);
26  fit->SetParLimits(0, -0.15, 0.01);
27  fit->SetParLimits(1, 0.6 * min, 1.25 * min);
28  fit->SetParLimits(2, 0.1, 10);
29  fit->SetParameters(xofmin, min, (max - min) / fabs(xofmax - xofmin));
30 
31  int badfit = p->Fit(fit, "IEQ", "", -.5, .3);
32  if (badfit)
33  badfit = p->Fit(fit, "IEQ", "", -.46, .26);
34  if (badfit) {
35  book.erase(it);
36  }
37  }
38 }

References Book::begin(), Book::end(), Book::erase(), SiStripPI::max, method(), min(), AlCaHLTBitMon_ParallelJobs::p, tools::TF1, and WIDTH.

Referenced by fit().

◆ granularity()

poly< std::string > LA_Filler_Fitter::granularity ( const SiStripDetId  detid,
const float  tthetaL,
const Long64_t  TFE_index,
const float  localy,
const unsigned  apvstrip 
) const

Definition at line 95 of file LA_Filler.cc.

99  {
100  poly<std::string> gran;
101  gran += subdetLabel(detid);
102  if (byLayer_)
103  gran *= layerLabel(detid);
104  if (byModule_)
105  gran *= moduleLabel(detid);
106  if (localYbin_)
107  gran += (localy < 0 ? "_yM" : "_yP") + std::to_string(abs((int)(localy / localYbin_ + (localy < 0 ? -1 : 0))));
108  if (stripsPerBin_)
109  gran += "_strip" +
110  std::to_string(
111  (unsigned)((0.5 + ((apvstrip / 64) ? (127 - apvstrip) : apvstrip) / stripsPerBin_) * stripsPerBin_));
112  if (ensembleBins_) {
113  gran +=
114  "_ensembleBin" + std::to_string((int)(ensembleBins_ * (tthetaL - ensembleLow_) / (ensembleUp_ - ensembleLow_)));
115  gran += "";
116  if (ensembleSize_)
117  gran *= "_sample" + std::to_string(TFE_index % ensembleSize_);
118  }
119  return gran;
120 }

References funct::abs(), byLayer_, byModule_, ensembleBins_, ensembleLow_, ensembleSize_, ensembleUp_, layerLabel(), localYbin_, moduleLabel(), stripsPerBin_, and subdetLabel().

Referenced by fill().

◆ layer_index()

static unsigned LA_Filler_Fitter::layer_index ( bool  TIB,
bool  stereo,
unsigned  layer 
)
inlinestatic

Definition at line 105 of file LA_Filler_Fitter.h.

105  {
106  return layer + (TIB ? 0 : 6) + (stereo ? 1 : 0) + ((layer > 2) ? 1 : (layer == 1) ? -1 : 0);
107  }

References SectorBuilder_Tib_cff::TIB.

Referenced by sistrip::MeasureLA::calibration_key(), and sistrip::EnsembleCalibrationLA::write_ensembles_text().

◆ layer_results()

std::map< std::string, LA_Filler_Fitter::Result > LA_Filler_Fitter::layer_results ( const Book &  book,
const Method  m 
)
static

Definition at line 71 of file LA_Results.cc.

71  {
72  std::map<std::string, Result> results;
73  for (Book::const_iterator it = book.begin(".*layer\\d.*" + method(m)); it != book.end(); ++it) {
74  const std::string name = boost::erase_all_copy(it->first, method(m));
75  results[name] = result(m, it->first, book);
76  }
77  return results;
78 }

References Book::begin(), Book::end(), visualization-live-secondInstance_cfg::m, method(), Skims_PA_cff::name, result(), bookConverter::results, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by sistrip::MeasureLA::process_reports().

◆ layerLabel()

std::string LA_Filler_Fitter::layerLabel ( const SiStripDetId  detid) const

Definition at line 128 of file LA_Filler.cc.

128  {
129  const bool isTIB = detid.subdetId() == StripSubdetector::TIB;
130  unsigned layer = isTIB ? tTopo_->tibLayer(detid) : tTopo_->tobLayer(detid);
131  bool stereo = isTIB ? tTopo_->tibStereo(detid) : tTopo_->tobStereo(detid);
132 
133  return subdetLabel(detid) + "_layer" + std::to_string(layer) + (stereo ? "s" : "a");
134 }

References DetId::subdetId(), subdetLabel(), StripSubdetector::TIB, TrackerTopology::tibLayer(), TrackerTopology::tibStereo(), TrackerTopology::tobLayer(), TrackerTopology::tobStereo(), and tTopo_.

Referenced by granularity(), and sistrip::MeasureLA::summarize_module_muH_byLayer().

◆ make_and_fit_symmchi2()

void LA_Filler_Fitter::make_and_fit_symmchi2 ( Book &  book)
static

Definition at line 40 of file LA_Fitter.cc.

40  {
41  for (Book::iterator it = book.begin(".*_all"); it != book.end(); ++it) {
42  const std::string base = boost::erase_all_copy(it->first, "_all");
43 
44  std::vector<Book::iterator> rebin_hists;
45  const Book::iterator& all = it;
46  rebin_hists.push_back(all);
47  Book::iterator w1 = book.find(base + "_w1");
48  rebin_hists.push_back(w1);
49  Book::iterator var_w2 = book.find(base + method(AVGV2, false));
50  rebin_hists.push_back(var_w2);
51  Book::iterator var_w3 = book.find(base + method(AVGV3, false));
52  rebin_hists.push_back(var_w3);
53 
54  const unsigned rebin = std::max(var_w2 == book.end() ? 0 : find_rebin(var_w2->second),
55  var_w3 == book.end() ? 0 : find_rebin(var_w3->second));
56  for (const auto& it : rebin_hists)
57  if (it != book.end())
58  it->second->Rebin(rebin > 1 ? rebin < 7 ? rebin : 6 : 1);
59 
60  TH1* const prob_w1 =
61  w1 == book.end() ? nullptr : subset_probability(base + method(PROB1, false), w1->second, all->second);
62  TH1* const rmsv_w2 =
63  var_w2 == book.end() ? nullptr : rms_profile(base + method(RMSV2, false), (TProfile* const)var_w2->second);
64  TH1* const rmsv_w3 =
65  var_w3 == book.end() ? nullptr : rms_profile(base + method(RMSV3, false), (TProfile* const)var_w3->second);
66 
67  std::vector<TH1*> fit_hists;
68  if (prob_w1) {
69  book.book(base + method(PROB1, false), prob_w1);
70  fit_hists.push_back(prob_w1);
71  prob_w1->SetTitle("Width==1 Probability;tan#theta_{t}-(dx/dz)_{reco}");
72  }
73  if (var_w2 != book.end()) {
74  book.book(base + method(RMSV2, false), rmsv_w2);
75  fit_hists.push_back(var_w2->second);
76  var_w2->second->SetTitle("Width==2 Mean Variance;tan#theta_{t}-(dx/dz)_{reco}");
77  fit_hists.push_back(rmsv_w2);
78  rmsv_w2->SetTitle("Width==2 RMS Variance;tan#theta_{t}-(dx/dz)_{reco}");
79  }
80  if (var_w3 != book.end()) {
81  book.book(base + method(RMSV3, false), rmsv_w3);
82  fit_hists.push_back(var_w3->second);
83  var_w3->second->SetTitle("Width==3 Mean Variance;tan#theta_{t}-(dx/dz)_{reco}");
84  fit_hists.push_back(rmsv_w3);
85  rmsv_w3->SetTitle("Width==3 RMS Variance;tan#theta_{t}-(dx/dz)_{reco}");
86  }
87 
88  if (fit_hists.empty())
89  continue;
90  const unsigned bins = fit_hists[0]->GetNbinsX();
91  const unsigned guess = fit_hists[0]->FindBin(0);
92  const std::pair<unsigned, unsigned> range(guess - bins / 30, guess + bins / 30);
93 
94  for (auto const& hist : fit_hists) {
95  TH1* const chi2 = SymmetryFit::symmetryChi2(hist, range);
96  if (chi2) {
97  book.book(chi2->GetName(), chi2);
98  chi2->SetTitle("Symmetry #chi^{2};tan#theta_{t}-(dx/dz)_{reco}");
99  }
100  }
101  }
102 }

References python.cmstools::all(), AVGV2, AVGV3, newFWLiteAna::base, Book::begin(), trigObjTnPSource_cfi::bins, Book::book(), hltPixelTracks_cff::chi2, Book::end(), Book::find(), find_rebin(), mps_check::guess, compare::hist, SiStripPI::max, method(), PROB1, FastTimerService_cff::range, rms_profile(), RMSV2, RMSV3, AlCaHLTBitMon_QueryRunRegistry::string, subset_probability(), and SymmetryFit::symmetryChi2().

Referenced by fit().

◆ method()

static std::string LA_Filler_Fitter::method ( Method  m,
bool  fit = true 
)
inlinestatic

Definition at line 61 of file LA_Filler_Fitter.h.

61  {
62  switch (m) {
63  case WIDTH:
64  return "_width_prof";
65  case PROB1:
66  return fit ? SymmetryFit::name(method(m, false)) : "_prob_w1";
67  case AVGV2:
68  return fit ? SymmetryFit::name(method(m, false)) : "_avg_var_w2";
69  case AVGV3:
70  return fit ? SymmetryFit::name(method(m, false)) : "_avg_var_w3";
71  case RMSV2:
72  return fit ? SymmetryFit::name(method(m, false)) : "_rms_var_w2";
73  case RMSV3:
74  return fit ? SymmetryFit::name(method(m, false)) : "_rms_var_w3";
75  default:
76  return "_UNKNOWN";
77  }
78  }

References AVGV2, AVGV3, visualization-live-secondInstance_cfg::m, SymmetryFit::name(), PROB1, RMSV2, RMSV3, and WIDTH.

Referenced by ensemble_results(), fill_one_cluster(), fit_width_profile(), layer_results(), make_and_fit_symmchi2(), module_results(), result(), sistrip::MeasureLA::summarize_module_muH_byLayer(), varWidth(), sistrip::MeasureLA::write_report_plots(), and sistrip::MeasureLA::write_report_text_ms().

◆ module_results()

std::map< uint32_t, LA_Filler_Fitter::Result > LA_Filler_Fitter::module_results ( const Book &  book,
const Method  m 
)
static

Definition at line 61 of file LA_Results.cc.

61  {
62  std::map<uint32_t, Result> results;
63  for (Book::const_iterator it = book.begin(".*_module\\d*" + method(m)); it != book.end(); ++it) {
64  const uint32_t detid = boost::lexical_cast<uint32_t>(
65  boost::regex_replace(it->first, boost::regex(".*_module(\\d*)_.*"), std::string("\\1")));
66  results[detid] = result(m, it->first, book);
67  }
68  return results;
69 }

References Book::begin(), Book::end(), visualization-live-secondInstance_cfg::m, method(), result(), bookConverter::results, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by sistrip::MeasureLA::process_reports(), and sistrip::MeasureLA::summarize_module_muH_byLayer().

◆ moduleLabel()

std::string LA_Filler_Fitter::moduleLabel ( const SiStripDetId  detid)
static

Definition at line 125 of file LA_Filler.cc.

125  {
126  return subdetLabel(detid) + "_module" + std::to_string(detid());
127 }

References subdetLabel().

Referenced by granularity().

◆ offset_slope()

std::pair< std::pair< float, float >, std::pair< float, float > > LA_Filler_Fitter::offset_slope ( const std::vector< EnsembleSummary > &  ensembles)
static

Definition at line 145 of file LA_Results.cc.

146  {
147  try {
148  std::vector<float> x, y, xerr, yerr;
149  for (auto const& ensemble : ensembles) {
150  x.push_back(ensemble.truth);
151  xerr.push_back(0);
152  y.push_back(ensemble.meanMeasured.first);
153  yerr.push_back(ensemble.meanMeasured.second);
154  }
155  TGraphErrors graph(x.size(), &(x[0]), &(y[0]), &(xerr[0]), &(yerr[0]));
156  //Need our own copy for thread safety
157  TF1 pol1("mypol1", "pol1");
158  graph.Fit(&pol1);
159  const TF1* const fit = graph.GetFunction("pol1");
160 
161  return std::make_pair(std::make_pair(fit->GetParameter(0), fit->GetParError(0)),
162  std::make_pair(fit->GetParameter(1), fit->GetParError(1)));
163  } catch (edm::Exception const& e) {
164  std::cerr << "Fitting Line Failed " << std::endl << e << std::endl;
165  return std::make_pair(std::make_pair(0, 0), std::make_pair(0, 0));
166  }
167 }

References beam_dqm_sourceclient-live_cfg::cerr, MillePedeFileConverter_cfg::e, tools::TF1, x, and y.

Referenced by sistrip::EnsembleCalibrationLA::write_ensembles_text().

◆ pull()

float LA_Filler_Fitter::pull ( const std::vector< EnsembleSummary > &  ensembles)
static

Definition at line 169 of file LA_Results.cc.

169  {
170  float p(0), w(0);
171  for (auto const& ensemble : ensembles) {
172  const float unc2 = pow(ensemble.pull.second, 2);
173  p += ensemble.pull.first / unc2;
174  w += 1 / unc2;
175  }
176  return p / w;
177 }

References AlCaHLTBitMon_ParallelJobs::p, funct::pow(), and w.

Referenced by ensemble_summary(), and sistrip::EnsembleCalibrationLA::write_ensembles_text().

◆ result()

LA_Filler_Fitter::Result LA_Filler_Fitter::result ( Method  m,
const std::string  name,
const Book &  book 
)
static

Definition at line 12 of file LA_Results.cc.

12  {
13  Result p;
14  const std::string base = boost::erase_all_copy(name, method(m));
15 
16  const TH1* const h = book[name];
17  const TH1* const reco = book[base + "_reconstruction"];
18  const TH1* const field = book[base + "_field"];
19 
20  if (reco)
21  p.reco = std::make_pair<float, float>(reco->GetMean(), reco->GetMeanError());
22  if (field)
23  p.field = field->GetMean();
24  if (h) {
25  switch (m) {
26  case WIDTH: {
27  const TF1* const f = h->GetFunction("LA_profile_fit");
28  if (!f)
29  break;
30  p.measured = std::make_pair<float, float>(f->GetParameter(0), f->GetParError(0));
31  p.chi2 = f->GetChisquare();
32  p.ndof = f->GetNDF();
33  p.entries = (unsigned)(h->GetEntries());
34  break;
35  }
36  case PROB1:
37  case AVGV2:
38  case AVGV3:
39  case RMSV2:
40  case RMSV3: /*case MULTI:*/ {
41  const TF1* const f = h->GetFunction("SymmetryFit");
42  if (!f)
43  break;
44  p.measured = std::make_pair<float, float>(p.reco.first + f->GetParameter(0), f->GetParameter(1));
45  p.chi2 = f->GetParameter(2);
46  p.ndof = (unsigned)(f->GetParameter(3));
47  p.entries = (m & PROB1)
48  ? (unsigned)book[base + "_w1"]->GetEntries()
49  : (m & (AVGV2 | RMSV2))
50  ? (unsigned)book[base + method(AVGV2, false)]->GetEntries()
51  : (m & (AVGV3 | RMSV3)) ? (unsigned)book[base + method(AVGV3, false)]->GetEntries() : 0;
52  break;
53  }
54  default:
55  break;
56  }
57  }
58  return p;
59 }

References AVGV2, AVGV3, newFWLiteAna::base, f, visualization-live-secondInstance_cfg::m, method(), Skims_PA_cff::name, AlCaHLTBitMon_ParallelJobs::p, PROB1, RMSV2, RMSV3, AlCaHLTBitMon_QueryRunRegistry::string, tools::TF1, and WIDTH.

Referenced by ensemble_results(), layer_results(), and module_results().

◆ rms_profile()

TH1 * LA_Filler_Fitter::rms_profile ( const std::string  name,
const TProfile * const  prof 
)
static

Definition at line 121 of file LA_Fitter.cc.

121  {
122  const int bins = prof->GetNbinsX();
123  TH1* const rms =
124  new TH1F(name.c_str(), "", bins, prof->GetBinLowEdge(1), prof->GetBinLowEdge(bins) + prof->GetBinWidth(bins));
125  for (int i = 1; i <= bins; i++) {
126  const double Me = prof->GetBinError(i);
127  const double neff = prof->GetBinEntries(
128  i); //Should be prof->GetBinEffectiveEntries(i);, not availible this version ROOT. This is only ok for unweighted fills
129  rms->SetBinContent(i, Me * sqrt(neff));
130  rms->SetBinError(i, Me / sqrt(2));
131  }
132  return rms;
133 }

References trigObjTnPSource_cfi::bins, mps_fire::i, Skims_PA_cff::name, SiStripPI::rms, and mathSSE::sqrt().

Referenced by make_and_fit_symmchi2().

◆ subdetLabel()

std::string LA_Filler_Fitter::subdetLabel ( const SiStripDetId  detid)
static

Definition at line 122 of file LA_Filler.cc.

122  {
123  return detid.subDetector() == SiStripDetId::TOB ? "TOB" : "TIB";
124 }

References SiStripDetId::subDetector(), and SiStripDetId::TOB.

Referenced by granularity(), layerLabel(), and moduleLabel().

◆ subset_probability()

TH1 * LA_Filler_Fitter::subset_probability ( const std::string  name,
const TH1 * const  subset,
const TH1 * const  total 
)
static

Definition at line 135 of file LA_Fitter.cc.

135  {
136  const int bins = subset->GetNbinsX();
137  TH1* const prob = new TH1F(
138  name.c_str(), "", bins, subset->GetBinLowEdge(1), subset->GetBinLowEdge(bins) + subset->GetBinWidth(bins));
139  for (int i = 1; i <= bins; i++) {
140  const double s = subset->GetBinContent(i);
141  const double T = total->GetBinContent(i);
142  const double B = T - s;
143 
144  const double p = T ? s / T : 0;
145  const double perr = T ? ((s && B) ? sqrt(s * s * B + B * B * s) / (T * T) : 1 / T) : 0;
146 
147  prob->SetBinContent(i, p);
148  prob->SetBinError(i, perr);
149  }
150  return prob;
151 }

References TtFullHadDaughter::B, trigObjTnPSource_cfi::bins, mps_fire::i, Skims_PA_cff::name, AlCaHLTBitMon_ParallelJobs::p, TtFullHadEvtBuilder_cfi::prob, alignCSCRings::s, mathSSE::sqrt(), and dqmMemoryStats::total.

Referenced by make_and_fit_symmchi2().

◆ summarize_ensembles()

void LA_Filler_Fitter::summarize_ensembles ( Book &  book) const

Definition at line 90 of file LA_Results.cc.

90  {
91  typedef std::map<std::string, std::vector<Result> > results_t;
92  results_t results;
93  for (int m = FIRST_METHOD; m <= LAST_METHOD; m <<= 1)
94  if (methods_ & m) {
95  results_t g = ensemble_results(book, (Method)(m));
96  results.insert(g.begin(), g.end());
97  }
98 
99  for (auto const& group : results) {
100  const std::string name = group.first;
101  for (auto const& p : group.second) {
102  const float pad = (ensembleUp_ - ensembleLow_) / 10;
103  book.fill(p.reco.first, name + "_ensembleReco", 12 * ensembleBins_, ensembleLow_ - pad, ensembleUp_ + pad);
104  book.fill(p.measured.first, name + "_measure", 12 * ensembleBins_, ensembleLow_ - pad, ensembleUp_ + pad);
105  book.fill(p.measured.second, name + "_merr", 500, 0, 0.01);
106  book.fill((p.measured.first - p.reco.first) / p.measured.second, name + "_pull", 500, -10, 10);
107  }
108  //Need our own copy for thread safety
109  TF1 gaus("mygaus", "gaus");
110  book[name + "_measure"]->Fit(&gaus, "LLQ");
111  book[name + "_merr"]->Fit(&gaus, "LLQ");
112  book[name + "_pull"]->Fit(&gaus, "LLQ");
113  }
114 }

References ensemble_results(), ensembleBins_, ensembleLow_, ensembleUp_, Book::fill(), FIRST_METHOD, g, watchdog::group, LAST_METHOD, visualization-live-secondInstance_cfg::m, methods_, Skims_PA_cff::name, AlCaHLTBitMon_ParallelJobs::p, bookConverter::results, AlCaHLTBitMon_QueryRunRegistry::string, and tools::TF1.

Referenced by sistrip::EnsembleCalibrationLA::endJob().

◆ varWidth()

poly< std::string > LA_Filler_Fitter::varWidth ( const unsigned  width) const

Definition at line 85 of file LA_Filler.cc.

85  {
86  poly<std::string> vw;
87  vw++;
88  if (width == 2 && methods_ & (AVGV2 | RMSV2))
89  vw *= method(AVGV2, false);
90  if (width == 3 && methods_ & (AVGV3 | RMSV3))
91  vw *= method(AVGV3, false);
92  return vw;
93 }

References AVGV2, AVGV3, method(), methods_, RMSV2, RMSV3, and ApeEstimator_cff::width.

Referenced by fill_one_cluster().

Member Data Documentation

◆ byLayer_

const bool LA_Filler_Fitter::byLayer_
private

Definition at line 133 of file LA_Filler_Fitter.h.

Referenced by granularity().

◆ byModule_

const bool LA_Filler_Fitter::byModule_
private

Definition at line 133 of file LA_Filler_Fitter.h.

Referenced by granularity().

◆ ensembleBins_

const int LA_Filler_Fitter::ensembleBins_
private

Definition at line 131 of file LA_Filler_Fitter.h.

Referenced by fill(), fill_one_cluster(), granularity(), and summarize_ensembles().

◆ ensembleLow_

const double LA_Filler_Fitter::ensembleLow_
private

Definition at line 132 of file LA_Filler_Fitter.h.

Referenced by granularity(), and summarize_ensembles().

◆ ensembleSize_

const int LA_Filler_Fitter::ensembleSize_
private

Definition at line 131 of file LA_Filler_Fitter.h.

Referenced by granularity().

◆ ensembleUp_

const double LA_Filler_Fitter::ensembleUp_
private

Definition at line 132 of file LA_Filler_Fitter.h.

Referenced by granularity(), and summarize_ensembles().

◆ localYbin_

const float LA_Filler_Fitter::localYbin_
private

Definition at line 134 of file LA_Filler_Fitter.h.

Referenced by fill(), and granularity().

◆ maxEvents_

const Long64_t LA_Filler_Fitter::maxEvents_
private

Definition at line 136 of file LA_Filler_Fitter.h.

Referenced by fill().

◆ methods_

const int LA_Filler_Fitter::methods_
private

Definition at line 137 of file LA_Filler_Fitter.h.

Referenced by fill_one_cluster(), summarize_ensembles(), and varWidth().

◆ stripsPerBin_

const unsigned LA_Filler_Fitter::stripsPerBin_
private

Definition at line 135 of file LA_Filler_Fitter.h.

Referenced by fill(), and granularity().

◆ tTopo_

const TrackerTopology* LA_Filler_Fitter::tTopo_
private

Definition at line 138 of file LA_Filler_Fitter.h.

Referenced by layerLabel().

LA_Filler_Fitter::Method
Method
Definition: LA_Filler_Fitter.h:50
LA_Filler_Fitter::ensembleSize_
const int ensembleSize_
Definition: LA_Filler_Fitter.h:131
SectorBuilder_Tib_cff.TIB
TIB
Definition: SectorBuilder_Tib_cff.py:16
FastTimerService_cff.range
range
Definition: FastTimerService_cff.py:34
DDAxes::y
ApeEstimator_cff.width
width
Definition: ApeEstimator_cff.py:24
mps_fire.i
i
Definition: mps_fire.py:355
SiStripPI::mean
Definition: SiStripPayloadInspectorHelper.h:169
LA_Filler_Fitter::find_rebin
static unsigned find_rebin(const TH1 *const)
Definition: LA_Fitter.cc:104
LA_Filler_Fitter::PROB1
Definition: LA_Filler_Fitter.h:53
f
double f[11][100]
Definition: MuScleFitUtils.cc:78
LA_Filler_Fitter::methods_
const int methods_
Definition: LA_Filler_Fitter.h:137
min
T min(T a, T b)
Definition: MathUtil.h:58
LA_Filler_Fitter::byModule_
const bool byModule_
Definition: LA_Filler_Fitter.h:133
AlCaHLTBitMon_ParallelJobs.p
p
Definition: AlCaHLTBitMon_ParallelJobs.py:153
tree
Definition: tree.py:1
LA_Filler_Fitter::rms_profile
static TH1 * rms_profile(const std::string, const TProfile *const)
Definition: LA_Fitter.cc:121
edm::swap
void swap(Association< C > &lhs, Association< C > &rhs)
Definition: Association.h:117
LA_Filler_Fitter::allAndOne
poly< std::string > allAndOne(const unsigned width) const
Definition: LA_Filler.cc:78
LA_Filler_Fitter::tTopo_
const TrackerTopology * tTopo_
Definition: LA_Filler_Fitter.h:138
Book::const_iterator
boost::filter_iterator< match_name, book_t::const_iterator > const_iterator
Definition: Book.h:56
Validation_hcalonly_cfi.sign
sign
Definition: Validation_hcalonly_cfi.py:32
bookConverter.results
results
Definition: bookConverter.py:144
Book::book
TH1 * book(string_t name, TH1 *const hist)
Definition: Book.h:42
DDAxes::x
SymmetryFit::name
static std::string name(std::string base)
Definition: SymmetryFit.h:12
reco
fixed size matrix
Definition: AlignmentAlgorithmBase.h:45
hltPixelTracks_cff.chi2
chi2
Definition: hltPixelTracks_cff.py:25
tools.TF1
TF1
Definition: tools.py:23
SiStripPI::rms
Definition: SiStripPayloadInspectorHelper.h:169
python.cmstools.all
def all(container)
workaround iterator generators for ROOT classes
Definition: cmstools.py:26
end
#define end
Definition: vmac.h:39
TrackerTopology::tobStereo
uint32_t tobStereo(const DetId &id) const
Definition: TrackerTopology.h:275
edm::Exception
Definition: EDMException.h:77
LA_Filler_Fitter::AVGV2
Definition: LA_Filler_Fitter.h:54
LA_Filler_Fitter::subdetLabel
static std::string subdetLabel(const SiStripDetId)
Definition: LA_Filler.cc:122
LA_Filler_Fitter::layerLabel
std::string layerLabel(const SiStripDetId) const
Definition: LA_Filler.cc:128
LA_Filler_Fitter::fill_one_cluster
void fill_one_cluster(Book &, const poly< std::string > &, const unsigned, const float, const float, const float, const float) const
Definition: LA_Filler.cc:60
alignCSCRings.s
s
Definition: alignCSCRings.py:92
LA_Filler_Fitter::RMSV3
Definition: LA_Filler_Fitter.h:57
funct::cos
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
SymmetryFit::symmetryChi2
static TH1 * symmetryChi2(std::string, const std::vector< TH1 * > &, const std::pair< unsigned, unsigned >)
Definition: SymmetryFit.cc:6
compare.hist
hist
Definition: compare.py:376
w
const double w
Definition: UKUtility.cc:23
visualization-live-secondInstance_cfg.m
m
Definition: visualization-live-secondInstance_cfg.py:72
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:19
LA_Filler_Fitter::byLayer_
const bool byLayer_
Definition: LA_Filler_Fitter.h:133
LA_Filler_Fitter::method
static std::string method(Method m, bool fit=true)
Definition: LA_Filler_Fitter.h:61
N
#define N
Definition: blowfish.cc:9
LA_Filler_Fitter::varWidth
poly< std::string > varWidth(const unsigned width) const
Definition: LA_Filler.cc:85
Book::end
iterator end(string_t re=".*")
Definition: Book.h:65
LA_Filler_Fitter::ensembleBins_
const int ensembleBins_
Definition: LA_Filler_Fitter.h:131
StripSubdetector::TIB
static constexpr auto TIB
Definition: StripSubdetector.h:16
LA_Filler_Fitter::WIDTH
Definition: LA_Filler_Fitter.h:51
LA_Filler_Fitter::ensemble_results
static std::map< std::string, std::vector< Result > > ensemble_results(const Book &, const Method)
Definition: LA_Results.cc:80
h
LA_Filler_Fitter::moduleLabel
static std::string moduleLabel(const SiStripDetId)
Definition: LA_Filler.cc:125
Book::find
iterator find(string_t name, string_t re=".*")
Definition: Book.h:73
Book::erase
void erase(string_t name)
Definition: Book.h:84
DetId::subdetId
constexpr int subdetId() const
get the contents of the subdetector field (not cast into any detector's numbering enum)
Definition: DetId.h:48
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
LA_Filler_Fitter::LAST_METHOD
Definition: LA_Filler_Fitter.h:58
TrackerTopology::tibStereo
uint32_t tibStereo(const DetId &id) const
Definition: TrackerTopology.h:279
LA_Filler_Fitter::ensembleUp_
const double ensembleUp_
Definition: LA_Filler_Fitter.h:132
Book::fill
void fill(double_t X, const char *name, uint_t NbinsX, double_t Xlow, double_t Xup, double_t W=1)
Definition: Book.h:96
LA_Filler_Fitter::make_and_fit_symmchi2
static void make_and_fit_symmchi2(Book &)
Definition: LA_Fitter.cc:40
SiStripPI::max
Definition: SiStripPayloadInspectorHelper.h:169
beam_dqm_sourceclient-live_cfg.cerr
cerr
Definition: beam_dqm_sourceclient-live_cfg.py:17
funct::tan
Tan< T >::type tan(const T &t)
Definition: Tan.h:22
LA_Filler_Fitter::stripsPerBin_
const unsigned stripsPerBin_
Definition: LA_Filler_Fitter.h:135
edmLumisInFiles.summary
summary
Definition: edmLumisInFiles.py:39
LA_Filler_Fitter::subset_probability
static TH1 * subset_probability(const std::string name, const TH1 *const, const TH1 *const)
Definition: LA_Fitter.cc:135
SiStripDetId::TOB
static constexpr auto TOB
Definition: SiStripDetId.h:39
LA_Filler_Fitter::pull
static float pull(const std::vector< EnsembleSummary > &)
Definition: LA_Results.cc:169
module
Definition: vlib.h:198
LA_Filler_Fitter::granularity
poly< std::string > granularity(const SiStripDetId, const float, const Long64_t, const float, const unsigned) const
Definition: LA_Filler.cc:95
LA_Filler_Fitter::AVGV3
Definition: LA_Filler_Fitter.h:55
Book::iterator
boost::filter_iterator< match_name, book_t::iterator > iterator
Definition: Book.h:55
LA_Filler_Fitter::ensembleLow_
const double ensembleLow_
Definition: LA_Filler_Fitter.h:132
Book::begin
iterator begin(string_t re=".*")
Definition: Book.h:57
TtFullHadDaughter::B
static const std::string B
Definition: TtFullHadronicEvent.h:9
TrackerTopology::tobLayer
unsigned int tobLayer(const DetId &id) const
Definition: TrackerTopology.h:147
poly
Definition: poly.h:11
T
long double T
Definition: Basic3DVectorLD.h:48
LA_Filler_Fitter::fit_width_profile
static void fit_width_profile(Book &)
Definition: LA_Fitter.cc:8
LA_Filler_Fitter::maxEvents_
const Long64_t maxEvents_
Definition: LA_Filler_Fitter.h:136
SiStripDetId::subDetector
SubDetector subDetector() const
Definition: SiStripDetId.h:105
LA_Filler_Fitter::result
static Result result(Method, const std::string name, const Book &)
Definition: LA_Results.cc:12
SiStripDetId::TIB
static constexpr auto TIB
Definition: SiStripDetId.h:37
Skims_PA_cff.name
name
Definition: Skims_PA_cff.py:17
LA_Filler_Fitter::RMSV2
Definition: LA_Filler_Fitter.h:56
dqmMemoryStats.total
total
Definition: dqmMemoryStats.py:152
funct::pow
Power< A, B >::type pow(const A &a, const B &b)
Definition: Power.h:30
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
trigObjTnPSource_cfi.bins
bins
Definition: trigObjTnPSource_cfi.py:20
mps_check.guess
string guess
Definition: mps_check.py:394
SiStripDetId
Detector identifier class for the strip tracker.
Definition: SiStripDetId.h:18
LA_Filler_Fitter::FIRST_METHOD
Definition: LA_Filler_Fitter.h:52
newFWLiteAna.base
base
Definition: newFWLiteAna.py:92
TtFullHadEvtBuilder_cfi.prob
prob
Definition: TtFullHadEvtBuilder_cfi.py:33
LaserClient_cfi.low
low
Definition: LaserClient_cfi.py:52
begin
#define begin
Definition: vmac.h:32
fit
Definition: CombinedChiSquaredLikelihood.h:6
up
Definition: BitonicSort.h:7
LA_Filler_Fitter::localYbin_
const float localYbin_
Definition: LA_Filler_Fitter.h:134
g
The Signals That Services Can Subscribe To This is based on ActivityRegistry and is current per Services can connect to the signals distributed by the ActivityRegistry in order to monitor the activity of the application Each possible callback has some defined which we here list in angle e g
Definition: Activities.doc:4
MillePedeFileConverter_cfg.e
e
Definition: MillePedeFileConverter_cfg.py:37
TrackerTopology::tibLayer
unsigned int tibLayer(const DetId &id) const
Definition: TrackerTopology.h:150
watchdog.group
group
Definition: watchdog.py:82
Generated for CMSSW Reference Manual by doxygen 1.8.16