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, int M, int N, double low, double up, unsigned max, const TrackerTopology *tTopo)
	LA_Filler_Fitter (int methods, bool layer, bool module, float localybin, unsigned stripbin, 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

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.

              {
    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
  };

Constructor & Destructor Documentation

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.

      : ensembleSize_(M),
        ensembleBins_(N),
        ensembleLow_(low),
        ensembleUp_(up),
        byLayer_(true),
        byModule_(false),
        localYbin_(0),
        stripsPerBin_(0),
        maxEvents_(max),
        methods_(methods),
        tTopo_(tTopo) {}

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.

      : ensembleSize_(0),
        ensembleBins_(0),
        ensembleLow_(0),
        ensembleUp_(0),
        byLayer_(layer),
        byModule_(module),
        localYbin_(localybin),
        stripsPerBin_(stripbin),
        maxEvents_(max),
        methods_(methods),
        tTopo_(tTopo) {}

Member Function Documentation

poly< std::string > LA_Filler_Fitter::allAndOne

(

const unsigned

width

)

const

Definition at line 78 of file LA_Filler.cc.

Referenced by LA_Filler_Fitter::EnsembleSummary::EnsembleSummary(), and fill_one_cluster().

                                                                      {
  poly<std::string> a1("_all");
  if (width == 1)
    a1 *= "_w1";
  return a1;
}

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.

References Book::begin(), Book::end(), method(), Skims_PA_cff::name, result(), bookConverter::results, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by fit(), and summarize_ensembles().

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

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.

References newFWLiteAna::base, Book::begin(), Book::end(), LA_Filler_Fitter::EnsembleSummary::meanMeasured, LA_Filler_Fitter::EnsembleSummary::meanUncertainty, Skims_PA_cff::name, LA_Filler_Fitter::EnsembleSummary::pull, pull(), alignCSCRings::s, LA_Filler_Fitter::EnsembleSummary::samples, LA_Filler_Fitter::EnsembleSummary::sigmaMeasured, AlCaHLTBitMon_QueryRunRegistry::string, edmLumisInFiles::summary, and LA_Filler_Fitter::EnsembleSummary::truth.

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

                      {
  std::map<std::string, std::vector<EnsembleSummary> > summary;
  for (Book::const_iterator it = book.begin(".*_ensembleReco"); it != book.end(); ++it) {
    const std::string base = boost::erase_all_copy(it->first, "_ensembleReco");

    const TH1* const reco = it->second;
    const TH1* const measure = book[base + "_measure"];
    const TH1* const merr = book[base + "_merr"];
    const TH1* const pull = book[base + "_pull"];

    EnsembleSummary s;
    s.samples = reco->GetEntries();
    s.truth = reco->GetMean();
    s.meanMeasured = std::make_pair<float, float>(measure->GetFunction("gaus")->GetParameter(1),
                                                  measure->GetFunction("gaus")->GetParError(1));
    s.sigmaMeasured = std::make_pair<float, float>(measure->GetFunction("gaus")->GetParameter(2),
                                                   measure->GetFunction("gaus")->GetParError(2));
    s.meanUncertainty = std::make_pair<float, float>(merr->GetFunction("gaus")->GetParameter(1),
                                                     merr->GetFunction("gaus")->GetParError(1));
    s.pull = std::make_pair<float, float>(pull->GetFunction("gaus")->GetParameter(2),
                                          pull->GetFunction("gaus")->GetParError(2));

    const std::string name = boost::regex_replace(base, boost::regex("ensembleBin\\d*_"), "");
    summary[name].push_back(s);
  }
  return summary;
}

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

Definition at line 8 of file LA_Filler.cc.

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(), LA_Filler_Fitter::EnsembleSummary::EnsembleSummary(), and sistrip::MeasureLA::MeasureLA().

                                                         {
  TTREE_FOREACH_ENTRY(tree) {
    TFE_MAX(maxEvents_);
    TFE_PRINTSTATUS;
    std::vector<unsigned> PLEAF(tsostrackmulti, tree);
    std::vector<unsigned> PLEAF(clusterdetid, tree);
    std::vector<unsigned> PLEAF(clusterwidth, tree);
    std::vector<float> PLEAF(clustervariance, tree);
    std::vector<float> PLEAF(tsosdriftx, tree);
    std::vector<float> PLEAF(tsosdriftz, tree);
    std::vector<float> PLEAF(tsoslocaltheta, tree);
    std::vector<float> PLEAF(tsoslocalphi, tree);
    std::vector<float> PLEAF(tsosglobalZofunitlocalY, tree);

    const unsigned N(clusterdetid.size());
    std::vector<float> BdotY(N, 0);
    if (!ensembleBins_) {
      std::vector<float> PLEAF(tsosBdotY, tree);
      swap(BdotY, tsosBdotY);
    }
    std::vector<float> localy(N, 0);
    if (localYbin_) {
      std::vector<float> PLEAF(tsoslocaly, tree);
      swap(localy, tsoslocaly);
    }
    std::vector<unsigned> seedstrip(N, 0);
    if (stripsPerBin_) {
      std::vector<unsigned> PLEAF(clusterseedstrip, tree);
      swap(seedstrip, clusterseedstrip);
    }

    for (unsigned i = 0; i < N; i++) {
      const SiStripDetId detid(clusterdetid[i]);
      if (tsostrackmulti[i] != 1 ||
          (detid.subDetector() != SiStripDetId::TIB && detid.subDetector() != SiStripDetId::TOB))
        continue;

      const int sign = tsosglobalZofunitlocalY[i] < 0 ? -1 : 1;
      const float tthetaL = sign * tsosdriftx[i] / tsosdriftz[i];
      const float tthetaT = sign * tan(tsoslocaltheta[i]) * cos(tsoslocalphi[i]);

      fill_one_cluster(book,
                       granularity(detid, tthetaL, TFE_index, localy[i], seedstrip[i] % 128),
                       clusterwidth[i],
                       clustervariance[i],
                       tthetaL,
                       tthetaT,
                       fabs(BdotY[i]));
    }
  }
}

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.

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

Referenced by LA_Filler_Fitter::EnsembleSummary::EnsembleSummary(), and fill().

                                                                 {
  book.fill(tthetaL, gran + "_reconstruction", 360, -1.0, 1.0);
  book.fill(tthetaT - tthetaL, gran + allAndOne(width), 360, -1.0, 1.0);
  book.fill(tthetaT - tthetaL, variance, gran + varWidth(width), 360, -1.0, 1.0);
  if (methods_ & WIDTH)
    book.fill(tthetaT, width, gran + method(WIDTH), 81, -0.6, 0.6);
  if (!ensembleBins_) {
    book.fill(BdotY, gran + "_field", 101, 1, 5);
    book.fill(width, gran + "_width", 10, 0, 10);
  }
}

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

static

Definition at line 104 of file LA_Fitter.cc.

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

Referenced by fit(), and make_and_fit_symmchi2().

                                                           {
  const double mean = hist->GetMean();
  const double rms = hist->GetRMS();
  const int begin = std::min(1, hist->GetXaxis()->FindFixBin(mean - rms));
  const int end = std::max(hist->GetNbinsX(), hist->GetXaxis()->FindFixBin(mean + rms)) + 1;
  unsigned current_hole(0), max_hole(0);
  for (int i = begin; i < end; i++) {
    if (!hist->GetBinError(i))
      current_hole++;
    else if (current_hole) {
      max_hole = std::max(current_hole, max_hole);
      current_hole = 0;
    }
  }
  return max_hole + 1;
}

static void LA_Filler_Fitter::fit ( Book &  book )

inlinestatic

Definition at line 110 of file LA_Filler_Fitter.h.

References ensemble_results(), ensemble_summary(), find_rebin(), fit_width_profile(), layer_results(), make_and_fit_symmchi2(), module_results(), Skims_PA_cff::name, offset_slope(), pull(), result(), rms_profile(), AlCaHLTBitMon_QueryRunRegistry::string, subset_probability(), and summarize_ensembles().

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

                              {
    make_and_fit_symmchi2(book);
    fit_width_profile(book);
  }

void LA_Filler_Fitter::fit_width_profile ( Book &  book )

static

Definition at line 8 of file LA_Fitter.cc.

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

Referenced by fit().

                                                   {
  for (Book::iterator it = book.begin(".*" + method(WIDTH)); it != book.end(); ++it) {
    it->second->SetTitle("Mean Cluster Width;tan#theta_{t}");
    TH1* const p = it->second;
    if (p->GetEntries() < 400) {
      delete p;
      book[it->first] = nullptr;
      continue;
    }
    p->SetTitle(";tan#theta_{t};");
    const float min = p->GetMinimum();
    const float max = p->GetMaximum();
    float xofmin = p->GetBinCenter(p->GetMinimumBin());
    if (xofmin > 0.0 || xofmin < -0.15)
      xofmin = -0.05;
    const float xofmax = p->GetBinCenter(p->GetMaximumBin());

    TF1* const fit = new TF1("LA_profile_fit", "[2]*(TMath::Abs(x-[0]))+[1]", -1, 1);
    fit->SetParLimits(0, -0.15, 0.01);
    fit->SetParLimits(1, 0.6 * min, 1.25 * min);
    fit->SetParLimits(2, 0.1, 10);
    fit->SetParameters(xofmin, min, (max - min) / fabs(xofmax - xofmin));

    int badfit = p->Fit(fit, "IEQ", "", -.5, .3);
    if (badfit)
      badfit = p->Fit(fit, "IEQ", "", -.46, .26);
    if (badfit) {
      book.erase(it);
    }
  }
}

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.

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

Referenced by LA_Filler_Fitter::EnsembleSummary::EnsembleSummary(), and fill().

                                                                               {
  poly<std::string> gran;
  gran += subdetLabel(detid);
  if (byLayer_)
    gran *= layerLabel(detid);
  if (byModule_)
    gran *= moduleLabel(detid);
  if (localYbin_)
    gran += (localy < 0 ? "_yM" : "_yP") + std::to_string(abs((int)(localy / localYbin_ + (localy < 0 ? -1 : 0))));
  if (stripsPerBin_)
    gran += "_strip" +
            std::to_string(
                (unsigned)((0.5 + ((apvstrip / 64) ? (127 - apvstrip) : apvstrip) / stripsPerBin_) * stripsPerBin_));
  if (ensembleBins_) {
    gran +=
        "_ensembleBin" + std::to_string((int)(ensembleBins_ * (tthetaL - ensembleLow_) / (ensembleUp_ - ensembleLow_)));
    gran += "";
    if (ensembleSize_)
      gran *= "_sample" + std::to_string(TFE_index % ensembleSize_);
  }
  return gran;
}

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

inlinestatic

Definition at line 105 of file LA_Filler_Fitter.h.

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

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

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.

References Book::begin(), Book::end(), method(), Skims_PA_cff::name, result(), bookConverter::results, and AlCaHLTBitMon_QueryRunRegistry::string.

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

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

std::string LA_Filler_Fitter::layerLabel

(

const SiStripDetId

detid

)

const

Definition at line 128 of file LA_Filler.cc.

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

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

                                                                     {
  const bool isTIB = detid.subdetId() == StripSubdetector::TIB;
  unsigned layer = isTIB ? tTopo_->tibLayer(detid) : tTopo_->tobLayer(detid);
  bool stereo = isTIB ? tTopo_->tibStereo(detid) : tTopo_->tobStereo(detid);

  return subdetLabel(detid) + "_layer" + std::to_string(layer) + (stereo ? "s" : "a");
}

void LA_Filler_Fitter::make_and_fit_symmchi2 ( Book &  book )

static

Definition at line 40 of file LA_Fitter.cc.

References Vispa.Plugins.EdmBrowser.EdmDataAccessor::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().

                                                       {
  for (Book::iterator it = book.begin(".*_all"); it != book.end(); ++it) {
    const std::string base = boost::erase_all_copy(it->first, "_all");

    std::vector<Book::iterator> rebin_hists;
    const Book::iterator& all = it;
    rebin_hists.push_back(all);
    Book::iterator w1 = book.find(base + "_w1");
    rebin_hists.push_back(w1);
    Book::iterator var_w2 = book.find(base + method(AVGV2, false));
    rebin_hists.push_back(var_w2);
    Book::iterator var_w3 = book.find(base + method(AVGV3, false));
    rebin_hists.push_back(var_w3);

    const unsigned rebin = std::max(var_w2 == book.end() ? 0 : find_rebin(var_w2->second),
                                    var_w3 == book.end() ? 0 : find_rebin(var_w3->second));
    for (const auto& it : rebin_hists)
      if (it != book.end())
        it->second->Rebin(rebin > 1 ? rebin < 7 ? rebin : 6 : 1);

    TH1* const prob_w1 =
        w1 == book.end() ? nullptr : subset_probability(base + method(PROB1, false), w1->second, all->second);
    TH1* const rmsv_w2 =
        var_w2 == book.end() ? nullptr : rms_profile(base + method(RMSV2, false), (TProfile* const)var_w2->second);
    TH1* const rmsv_w3 =
        var_w3 == book.end() ? nullptr : rms_profile(base + method(RMSV3, false), (TProfile* const)var_w3->second);

    std::vector<TH1*> fit_hists;
    if (prob_w1) {
      book.book(base + method(PROB1, false), prob_w1);
      fit_hists.push_back(prob_w1);
      prob_w1->SetTitle("Width==1 Probability;tan#theta_{t}-(dx/dz)_{reco}");
    }
    if (var_w2 != book.end()) {
      book.book(base + method(RMSV2, false), rmsv_w2);
      fit_hists.push_back(var_w2->second);
      var_w2->second->SetTitle("Width==2 Mean Variance;tan#theta_{t}-(dx/dz)_{reco}");
      fit_hists.push_back(rmsv_w2);
      rmsv_w2->SetTitle("Width==2 RMS Variance;tan#theta_{t}-(dx/dz)_{reco}");
    }
    if (var_w3 != book.end()) {
      book.book(base + method(RMSV3, false), rmsv_w3);
      fit_hists.push_back(var_w3->second);
      var_w3->second->SetTitle("Width==3 Mean Variance;tan#theta_{t}-(dx/dz)_{reco}");
      fit_hists.push_back(rmsv_w3);
      rmsv_w3->SetTitle("Width==3 RMS Variance;tan#theta_{t}-(dx/dz)_{reco}");
    }

    if (fit_hists.empty())
      continue;
    const unsigned bins = fit_hists[0]->GetNbinsX();
    const unsigned guess = fit_hists[0]->FindBin(0);
    const std::pair<unsigned, unsigned> range(guess - bins / 30, guess + bins / 30);

    for (auto const& hist : fit_hists) {
      TH1* const chi2 = SymmetryFit::symmetryChi2(hist, range);
      if (chi2) {
        book.book(chi2->GetName(), chi2);
        chi2->SetTitle("Symmetry #chi^{2};tan#theta_{t}-(dx/dz)_{reco}");
      }
    }
  }
}

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

inlinestatic

Definition at line 61 of file LA_Filler_Fitter.h.

References AVGV2, AVGV3, 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().

                                                     {
    switch (m) {
      case WIDTH:
        return "_width_prof";
      case PROB1:
        return fit ? SymmetryFit::name(method(m, false)) : "_prob_w1";
      case AVGV2:
        return fit ? SymmetryFit::name(method(m, false)) : "_avg_var_w2";
      case AVGV3:
        return fit ? SymmetryFit::name(method(m, false)) : "_avg_var_w3";
      case RMSV2:
        return fit ? SymmetryFit::name(method(m, false)) : "_rms_var_w2";
      case RMSV3:
        return fit ? SymmetryFit::name(method(m, false)) : "_rms_var_w3";
      default:
        return "_UNKNOWN";
    }
  }

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.

References Book::begin(), Book::end(), method(), result(), bookConverter::results, and AlCaHLTBitMon_QueryRunRegistry::string.

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

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

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

static

Definition at line 125 of file LA_Filler.cc.

References subdetLabel().

Referenced by LA_Filler_Fitter::EnsembleSummary::EnsembleSummary(), and granularity().

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

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.

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

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

                                                                 {
  try {
    std::vector<float> x, y, xerr, yerr;
    for (auto const& ensemble : ensembles) {
      x.push_back(ensemble.truth);
      xerr.push_back(0);
      y.push_back(ensemble.meanMeasured.first);
      yerr.push_back(ensemble.meanMeasured.second);
    }
    TGraphErrors graph(x.size(), &(x[0]), &(y[0]), &(xerr[0]), &(yerr[0]));
    //Need our own copy for thread safety
    TF1 pol1("mypol1", "pol1");
    graph.Fit(&pol1);
    const TF1* const fit = graph.GetFunction("pol1");

    return std::make_pair(std::make_pair(fit->GetParameter(0), fit->GetParError(0)),
                          std::make_pair(fit->GetParameter(1), fit->GetParError(1)));
  } catch (edm::Exception const& e) {
    std::cerr << "Fitting Line Failed " << std::endl << e << std::endl;
    return std::make_pair(std::make_pair(0, 0), std::make_pair(0, 0));
  }
}

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

static

Definition at line 169 of file LA_Results.cc.

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

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

                                                                                        {
  float p(0), w(0);
  for (auto const& ensemble : ensembles) {
    const float unc2 = pow(ensemble.pull.second, 2);
    p += ensemble.pull.first / unc2;
    w += 1 / unc2;
  }
  return p / w;
}

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.

References AVGV2, AVGV3, newFWLiteAna::base, LA_Filler_Fitter::Result::chi2, LA_Filler_Fitter::Result::entries, f, LA_Filler_Fitter::Result::field, h, LA_Filler_Fitter::Result::measured, method(), Skims_PA_cff::name, LA_Filler_Fitter::Result::ndof, AlCaHLTBitMon_ParallelJobs::p, PROB1, LA_Filler_Fitter::Result::reco, RMSV2, RMSV3, AlCaHLTBitMon_QueryRunRegistry::string, and WIDTH.

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

                                                                                                {
  Result p;
  const std::string base = boost::erase_all_copy(name, method(m));

  const TH1* const h = book[name];
  const TH1* const reco = book[base + "_reconstruction"];
  const TH1* const field = book[base + "_field"];

  if (reco)
    p.reco = std::make_pair<float, float>(reco->GetMean(), reco->GetMeanError());
  if (field)
    p.field = field->GetMean();
  if (h) {
    switch (m) {
      case WIDTH: {
        const TF1* const f = h->GetFunction("LA_profile_fit");
        if (!f)
          break;
        p.measured = std::make_pair<float, float>(f->GetParameter(0), f->GetParError(0));
        p.chi2 = f->GetChisquare();
        p.ndof = f->GetNDF();
        p.entries = (unsigned)(h->GetEntries());
        break;
      }
      case PROB1:
      case AVGV2:
      case AVGV3:
      case RMSV2:
      case RMSV3: /*case MULTI:*/ {
        const TF1* const f = h->GetFunction("SymmetryFit");
        if (!f)
          break;
        p.measured = std::make_pair<float, float>(p.reco.first + f->GetParameter(0), f->GetParameter(1));
        p.chi2 = f->GetParameter(2);
        p.ndof = (unsigned)(f->GetParameter(3));
        p.entries = (m & PROB1)
                        ? (unsigned)book[base + "_w1"]->GetEntries()
                        : (m & (AVGV2 | RMSV2))
                              ? (unsigned)book[base + method(AVGV2, false)]->GetEntries()
                              : (m & (AVGV3 | RMSV3)) ? (unsigned)book[base + method(AVGV3, false)]->GetEntries() : 0;
        break;
      }
      default:
        break;
    }
  }
  return p;
}

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

static

Definition at line 121 of file LA_Fitter.cc.

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

Referenced by fit(), and make_and_fit_symmchi2().

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

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

static

Definition at line 122 of file LA_Filler.cc.

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

Referenced by LA_Filler_Fitter::EnsembleSummary::EnsembleSummary(), granularity(), layerLabel(), and moduleLabel().

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

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.

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

Referenced by fit(), and make_and_fit_symmchi2().

                                                                                                               {
  const int bins = subset->GetNbinsX();
  TH1* const prob = new TH1F(
      name.c_str(), "", bins, subset->GetBinLowEdge(1), subset->GetBinLowEdge(bins) + subset->GetBinWidth(bins));
  for (int i = 1; i <= bins; i++) {
    const double s = subset->GetBinContent(i);
    const double T = total->GetBinContent(i);
    const double B = T - s;

    const double p = T ? s / T : 0;
    const double perr = T ? ((s && B) ? sqrt(s * s * B + B * B * s) / (T * T) : 1 / T) : 0;

    prob->SetBinContent(i, p);
    prob->SetBinError(i, perr);
  }
  return prob;
}

void LA_Filler_Fitter::summarize_ensembles

(

Book &

book

)

const

Definition at line 90 of file LA_Results.cc.

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, and AlCaHLTBitMon_QueryRunRegistry::string.

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

                                                           {
  typedef std::map<std::string, std::vector<Result> > results_t;
  results_t results;
  for (int m = FIRST_METHOD; m <= LAST_METHOD; m <<= 1)
    if (methods_ & m) {
      results_t g = ensemble_results(book, (Method)(m));
      results.insert(g.begin(), g.end());
    }

  for (auto const& group : results) {
    const std::string name = group.first;
    for (auto const& p : group.second) {
      const float pad = (ensembleUp_ - ensembleLow_) / 10;
      book.fill(p.reco.first, name + "_ensembleReco", 12 * ensembleBins_, ensembleLow_ - pad, ensembleUp_ + pad);
      book.fill(p.measured.first, name + "_measure", 12 * ensembleBins_, ensembleLow_ - pad, ensembleUp_ + pad);
      book.fill(p.measured.second, name + "_merr", 500, 0, 0.01);
      book.fill((p.measured.first - p.reco.first) / p.measured.second, name + "_pull", 500, -10, 10);
    }
    //Need our own copy for thread safety
    TF1 gaus("mygaus", "gaus");
    book[name + "_measure"]->Fit(&gaus, "LLQ");
    book[name + "_merr"]->Fit(&gaus, "LLQ");
    book[name + "_pull"]->Fit(&gaus, "LLQ");
  }
}

poly< std::string > LA_Filler_Fitter::varWidth

(

const unsigned

width

)

const

Definition at line 85 of file LA_Filler.cc.

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

Referenced by LA_Filler_Fitter::EnsembleSummary::EnsembleSummary(), and fill_one_cluster().

                                                                     {
  poly<std::string> vw;
  vw++;
  if (width == 2 && methods_ & (AVGV2 | RMSV2))
    vw *= method(AVGV2, false);
  if (width == 3 && methods_ & (AVGV3 | RMSV3))
    vw *= method(AVGV3, false);
  return vw;
}

Member Data Documentation

const bool LA_Filler_Fitter::byLayer_

private

Definition at line 133 of file LA_Filler_Fitter.h.

Referenced by granularity().

const bool LA_Filler_Fitter::byModule_

private

Definition at line 133 of file LA_Filler_Fitter.h.

Referenced by granularity().

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().

const double LA_Filler_Fitter::ensembleLow_

private

Definition at line 132 of file LA_Filler_Fitter.h.

Referenced by granularity(), and summarize_ensembles().

const int LA_Filler_Fitter::ensembleSize_

private

Definition at line 131 of file LA_Filler_Fitter.h.

Referenced by granularity().

const double LA_Filler_Fitter::ensembleUp_

private

Definition at line 132 of file LA_Filler_Fitter.h.

Referenced by granularity(), and summarize_ensembles().

const float LA_Filler_Fitter::localYbin_

private

Definition at line 134 of file LA_Filler_Fitter.h.

Referenced by fill(), and granularity().

const Long64_t LA_Filler_Fitter::maxEvents_

private

Definition at line 136 of file LA_Filler_Fitter.h.

Referenced by fill().

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().

const unsigned LA_Filler_Fitter::stripsPerBin_

private

Definition at line 135 of file LA_Filler_Fitter.h.

Referenced by fill(), and granularity().

const TrackerTopology* LA_Filler_Fitter::tTopo_

private

Definition at line 138 of file LA_Filler_Fitter.h.

Referenced by layerLabel().