#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 40 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 19 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 30 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 61 of file LA_Filler.cc.

References varWidth().

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

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

poly

Definition: poly.h:11

ApeEstimator_cff.width

width

Definition: ApeEstimator_cff.py:24

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

References Book::begin(), Book::end(), method(), dataset::name, result(), mps_update::results, AlCaHLTBitMon_QueryRunRegistry::string, and summarize_ensembles().

Referenced by fit(), layer_results(), 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 107 of file LA_Results.cc.

References runEdmFileComparison::base, Book::begin(), Book::end(), LA_Filler_Fitter::EnsembleSummary::meanMeasured, LA_Filler_Fitter::EnsembleSummary::meanUncertainty, dataset::name, offset_slope(), LA_Filler_Fitter::EnsembleSummary::pull, pull(), alignCSCRings::s, LA_Filler_Fitter::EnsembleSummary::samples, LA_Filler_Fitter::EnsembleSummary::sigmaMeasured, AlCaHLTBitMon_QueryRunRegistry::string, heppy_report::summary, and LA_Filler_Fitter::EnsembleSummary::truth.

Referenced by fit(), summarize_ensembles(), 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 9 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 46 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 81 of file LA_Fitter.cc.

References begin, end, mps_fire::i, hpstanc_transforms::max, RecoTauDiscriminantConfiguration::mean, min(), RecoTauDiscriminantConfiguration::rms, and rms_profile().

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 89 of file LA_Filler_Fitter.h.

References ensemble_results(), ensemble_summary(), find_rebin(), fit_width_profile(), layer_results(), make_and_fit_symmchi2(), module_results(), dataset::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().

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

LA_Filler_Fitter::fit_width_profile

static void fit_width_profile(Book &)

Definition: LA_Fitter.cc:10

LA_Filler_Fitter::make_and_fit_symmchi2

static void make_and_fit_symmchi2(Book &)

Definition: LA_Fitter.cc:34

void LA_Filler_Fitter::fit_width_profile ( Book & book )

static

Definition at line 10 of file LA_Fitter.cc.

References Book::begin(), Book::end(), Book::erase(), make_and_fit_symmchi2(), hpstanc_transforms::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]=0; 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 73 of file LA_Filler.cc.

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

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

                                                                                                                                         {
   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 86 of file LA_Filler_Fitter.h.

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

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

SectorBuilder_Tib_cff.TIB

TIB

Definition: SectorBuilder_Tib_cff.py:16

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

static

Definition at line 63 of file LA_Results.cc.

References Book::begin(), Book::end(), ensemble_results(), method(), dataset::name, result(), mps_update::results, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by fit(), module_results(), 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 93 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(), moduleLabel(), 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 34 of file LA_Fitter.cc.

References Vispa.Plugins.EdmBrowser.EdmDataAccessor::all(), AVGV2, AVGV3, runEdmFileComparison::base, Book::begin(), create_public_pileup_plots::bins, Book::book(), HiEvtPlane_cfi::chi2, Book::end(), Book::find(), find_rebin(), create_public_lumi_plots::hist, hpstanc_transforms::max, method(), PROB1, rms_profile(), RMSV2, RMSV3, AlCaHLTBitMon_QueryRunRegistry::string, subset_probability(), and SymmetryFit::symmetryChi2().

Referenced by fit(), and fit_width_profile().

                                   {
   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;              
     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,0)); rebin_hists.push_back(var_w2);
     Book::iterator var_w3 = book.find(base+method(AVGV3,0)); 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) );
     BOOST_FOREACH(Book::iterator it, rebin_hists) if(it!=book.end()) it->second->Rebin( rebin>1 ? rebin<7 ? rebin : 6 : 1);
 
     TH1* const prob_w1 = w1==book.end()     ? 0 : subset_probability( base+method(PROB1,0) ,w1->second,all->second);
     TH1* const rmsv_w2 = var_w2==book.end() ? 0 :        rms_profile( base+method(RMSV2,0), (TProfile*const)var_w2->second);
     TH1* const rmsv_w3 = var_w3==book.end() ? 0 :        rms_profile( base+method(RMSV3,0), (TProfile*const)var_w3->second);
     
     std::vector<TH1*> fit_hists;
     if(prob_w1) {
       book.book(base+method(PROB1,0),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,0),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,0),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.size()) 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);
 
     BOOST_FOREACH(TH1*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 47 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,0)): "_prob_w1"   ;  
     case AVGV2: return fit? SymmetryFit::name(method(m,0)): "_avg_var_w2";  
     case AVGV3: return fit? SymmetryFit::name(method(m,0)): "_avg_var_w3";  
     case RMSV2: return fit? SymmetryFit::name(method(m,0)): "_rms_var_w2";  
     case RMSV3: return fit? SymmetryFit::name(method(m,0)): "_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 51 of file LA_Results.cc.

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

Referenced by fit(), sistrip::MeasureLA::process_reports(), result(), 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 91 of file LA_Filler.cc.

References layerLabel(), AlCaHLTBitMon_QueryRunRegistry::string, and subdetLabel().

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

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

LA_Filler_Fitter::subdetLabel

static std::string subdetLabel(const SiStripDetId)

Definition: LA_Filler.cc:89

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

static

Definition at line 132 of file LA_Results.cc.

References MessageLogger_cfi::cerr, MillePedeFileConverter_cfg::e, LA_Filler_Fitter::EnsembleSummary::meanMeasured, pull(), LA_Filler_Fitter::EnsembleSummary::truth, x, and y.

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

                                                                           { 
   try {
     std::vector<float> x,y,xerr,yerr;
     BOOST_FOREACH(EnsembleSummary 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 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 156 of file LA_Results.cc.

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

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

                                                                   {
   float p(0),w(0);
   BOOST_FOREACH(EnsembleSummary 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 14 of file LA_Results.cc.

References AVGV2, AVGV3, runEdmFileComparison::base, LA_Filler_Fitter::Result::chi2, LA_Filler_Fitter::Result::entries, f, LA_Filler_Fitter::Result::field, LA_Filler_Fitter::Result::measured, method(), module_results(), dataset::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,0)]->GetEntries() : 
     (m&(AVGV3|RMSV3)) ? (unsigned) book[base+method(AVGV3,0)]->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 95 of file LA_Fitter.cc.

References create_public_pileup_plots::bins, mps_fire::i, RecoTauDiscriminantConfiguration::rms, mathSSE::sqrt(), and subset_probability().

Referenced by find_rebin(), 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 89 of file LA_Filler.cc.

References moduleLabel(), AlCaHLTBitMon_QueryRunRegistry::string, SiStripDetId::subDetector(), and SiStripDetId::TOB.

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

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

SiStripDetId::TOB

Definition: SiStripDetId.h:36

SiStripDetId::subDetector

SubDetector subDetector() const

Definition: SiStripDetId.h:114

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

static

Definition at line 108 of file LA_Fitter.cc.

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

Referenced by fit(), make_and_fit_symmchi2(), and rms_profile().

                                                                                           {
   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 83 of file LA_Results.cc.

References ensemble_results(), ensemble_summary(), ensembleBins_, ensembleLow_, ensembleUp_, Book::fill(), FIRST_METHOD, g, edmIntegrityCheck::group, LAST_METHOD, funct::m, LA_Filler_Fitter::Result::measured, methods_, dataset::name, AlCaHLTBitMon_ParallelJobs::p, LA_Filler_Fitter::Result::reco, mps_update::results, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by sistrip::EnsembleCalibrationLA::endJob(), ensemble_results(), 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());}
   
   BOOST_FOREACH(const results_t::value_type group, results) {
     const std::string name = group.first;
     BOOST_FOREACH(const Result 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 65 of file LA_Filler.cc.

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

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

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

Member Data Documentation

const bool LA_Filler_Fitter::byLayer_

private

Definition at line 110 of file LA_Filler_Fitter.h.

Referenced by granularity().

const bool LA_Filler_Fitter::byModule_

private

Definition at line 110 of file LA_Filler_Fitter.h.

Referenced by granularity().

const int LA_Filler_Fitter::ensembleBins_

private

Definition at line 108 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 109 of file LA_Filler_Fitter.h.

Referenced by granularity(), and summarize_ensembles().

const int LA_Filler_Fitter::ensembleSize_

private

Definition at line 108 of file LA_Filler_Fitter.h.

Referenced by granularity().

const double LA_Filler_Fitter::ensembleUp_

private

Definition at line 109 of file LA_Filler_Fitter.h.

Referenced by granularity(), and summarize_ensembles().

const float LA_Filler_Fitter::localYbin_

private

Definition at line 111 of file LA_Filler_Fitter.h.

Referenced by fill(), and granularity().

const Long64_t LA_Filler_Fitter::maxEvents_

private

Definition at line 113 of file LA_Filler_Fitter.h.

Referenced by fill().

const int LA_Filler_Fitter::methods_

private

Definition at line 114 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 112 of file LA_Filler_Fitter.h.

Referenced by fill(), and granularity().

const TrackerTopology* LA_Filler_Fitter::tTopo_

private

Definition at line 115 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

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation