#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

◆ 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.

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

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

       : 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

◆ allAndOne()

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

Definition at line 77 of file LA_Filler.cc.

References isotrackTrainRegressor::a1, and ApeEstimator_cff::width.

Referenced by fill_one_cluster().

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

◆ 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 78 of file LA_Results.cc.

References Book::begin(), Book::end(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, visualization-live-secondInstance_cfg::m, method(), mergeVDriftHistosByStation::name, result(), mysort::results, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by 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 = std::regex_replace(it->first, std::regex("sample\\d*_"), "");
     results[name].push_back(result(m, it->first, book));
   }
   return results;
 }

◆ ensemble_summary()

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

static

Definition at line 114 of file LA_Results.cc.

References edmMakeDummyCfis::base, Book::begin(), Book::end(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, mergeVDriftHistosByStation::name, pull(), alignCSCRings::s, AlCaHLTBitMon_QueryRunRegistry::string, and edmLumisInFiles::summary.

Referenced by 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 = std::regex_replace(base, std::regex("ensembleBin\\d*_"), "");
     summary[name].push_back(s);
   }
   return summary;
 }

◆ fill()

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

Definition at line 7 of file LA_Filler.cc.

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

Referenced by sistrip::EnsembleCalibrationLA::endJob(), 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]));
     }
   }
 }

◆ 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 59 of file LA_Filler.cc.

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

Referenced by 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);
   }
 }

◆ find_rebin()

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

static

Definition at line 104 of file LA_Fitter.cc.

References SplitLinear::begin, mps_fire::end, compareTotals::hist, mps_fire::i, WZElectronSkims53X_cff::max, SiStripPI::mean, SiStripPI::min, and SiStripPI::rms.

Referenced by 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;
 }

◆ fit()

static void LA_Filler_Fitter::fit ( Book & book )

inlinestatic

Definition at line 110 of file LA_Filler_Fitter.h.

References fit_width_profile(), and make_and_fit_symmchi2().

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

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

◆ fit_width_profile()

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(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, WZElectronSkims53X_cff::max, method(), SiStripPI::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);
     }
   }
 }

◆ 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 94 of file LA_Filler.cc.

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

Referenced by 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;
 }

◆ 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.

References nano_mu_digi_cff::layer, and SectorBuilder_Tib_cff::TIB.

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

◆ 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 69 of file LA_Results.cc.

References Book::begin(), Book::end(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, visualization-live-secondInstance_cfg::m, method(), mergeVDriftHistosByStation::name, result(), mysort::results, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by 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;
 }

◆ layerLabel()

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

Definition at line 127 of file LA_Filler.cc.

References ALCARECOPPSCalTrackBasedSel_cff::detid, nano_mu_digi_cff::layer, subdetLabel(), StripSubdetector::TIB, TrackerTopology::tibLayer(), TrackerTopology::tibStereo(), to_string(), TrackerTopology::tobLayer(), TrackerTopology::tobStereo(), and tTopo_.

Referenced by 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");
 }

◆ make_and_fit_symmchi2()

void LA_Filler_Fitter::make_and_fit_symmchi2 ( Book & book )

static

Definition at line 40 of file LA_Fitter.cc.

References python.cmstools::all(), AVGV2, AVGV3, edmMakeDummyCfis::base, Book::begin(), isotrackApplyRegressor::bins, Book::book(), isoTrack_cff::chi2, Book::end(), Book::find(), find_rebin(), mps_check::guess, compareTotals::hist, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, WZElectronSkims53X_cff::max, method(), PROB1, isotrackApplyRegressor::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}");
       }
     }
   }
 }

◆ method()

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

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

◆ 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 59 of file LA_Results.cc.

References Book::begin(), ALCARECOPPSCalTrackBasedSel_cff::detid, Book::end(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, visualization-live-secondInstance_cfg::m, method(), result(), mysort::results, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by 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 =
         std::stoul(std::regex_replace(it->first, std::regex(".*_module(\\d*)_.*"), std::string("\\1")));
     results[detid] = result(m, it->first, book);
   }
   return results;
 }

◆ moduleLabel()

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

static

Definition at line 124 of file LA_Filler.cc.

References ALCARECOPPSCalTrackBasedSel_cff::detid, subdetLabel(), and to_string().

Referenced by granularity().

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

◆ 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 143 of file LA_Results.cc.

References DMR_cfg::cerr, MillePedeFileConverter_cfg::e, x, and y.

Referenced by 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));
   }
 }

◆ pull()

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

static

Definition at line 167 of file LA_Results.cc.

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

Referenced by ensemble_summary(), 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;
 }

◆ result()

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

static

Definition at line 11 of file LA_Results.cc.

References AVGV2, AVGV3, edmMakeDummyCfis::base, f, h, visualization-live-secondInstance_cfg::m, method(), mergeVDriftHistosByStation::name, AlCaHLTBitMon_ParallelJobs::p, PROB1, RMSV2, RMSV3, AlCaHLTBitMon_QueryRunRegistry::string, and WIDTH.

Referenced by ensemble_results(), 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;
 }

◆ 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.

References isotrackApplyRegressor::bins, mps_fire::i, mergeVDriftHistosByStation::name, SiStripPI::rms, and mathSSE::sqrt().

Referenced by 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;
 }

◆ subdetLabel()

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

static

Definition at line 121 of file LA_Filler.cc.

References ALCARECOPPSCalTrackBasedSel_cff::detid, and SiStripDetId::TOB.

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

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

◆ 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.

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

Referenced by 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;
 }

◆ summarize_ensembles()

void LA_Filler_Fitter::summarize_ensembles ( Book & book ) const

Definition at line 88 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_, mergeVDriftHistosByStation::name, AlCaHLTBitMon_ParallelJobs::p, mysort::results, and AlCaHLTBitMon_QueryRunRegistry::string.

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

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

◆ varWidth()

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

Definition at line 84 of file LA_Filler.cc.

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

Referenced by 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

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

Classes

Public Types

Public Member Functions

Static Public Member Functions

Private Attributes

Detailed Description

Member Enumeration Documentation

◆ Method

Constructor & Destructor Documentation

◆ LA_Filler_Fitter() [1/2]

◆ LA_Filler_Fitter() [2/2]

Member Function Documentation

◆ allAndOne()

◆ ensemble_results()

◆ ensemble_summary()

◆ fill()

◆ fill_one_cluster()

◆ find_rebin()

◆ fit()

◆ fit_width_profile()

◆ granularity()

◆ layer_index()

◆ layer_results()

◆ layerLabel()

◆ make_and_fit_symmchi2()

◆ method()

◆ module_results()

◆ moduleLabel()

◆ offset_slope()

◆ pull()

◆ result()

◆ rms_profile()

◆ subdetLabel()

◆ subset_probability()

◆ summarize_ensembles()

◆ varWidth()

Member Data Documentation

◆ byLayer_

◆ byModule_

◆ ensembleBins_

◆ ensembleLow_

◆ ensembleSize_

◆ ensembleUp_

◆ localYbin_

◆ maxEvents_

◆ methods_

◆ stripsPerBin_

◆ tTopo_