CMS 3D CMS Logo
List of all members | Public Member Functions | Private Member Functions | Private Attributes
ErrorsAnalyzer Class Reference

#include <MuonAnalysis/MomentumScaleCalibration/plugins/ErrorsAnalyzer.cc>

Inheritance diagram for ErrorsAnalyzer:
edm::EDAnalyzer edm::EDConsumerBase

Public Member Functions

 ErrorsAnalyzer (const edm::ParameterSet &)
 
 ~ErrorsAnalyzer () override
 
- Public Member Functions inherited from edm::EDAnalyzer
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
 EDAnalyzer ()
 
SerialTaskQueue * globalLuminosityBlocksQueue ()
 
SerialTaskQueue * globalRunsQueue ()
 
ModuleDescription const & moduleDescription () const
 
std::string workerType () const
 
 ~EDAnalyzer () override
 
- Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo > consumesInfo () const
 
void convertCurrentProcessAlias (std::string const &processName)
 Convert "@currentProcess" in InputTag process names to the actual current process name. More...
 
 EDConsumerBase ()
 
 EDConsumerBase (EDConsumerBase &&)=default
 
 EDConsumerBase (EDConsumerBase const &)=delete
 
ESProxyIndex const * esGetTokenIndices (edm::Transition iTrans) const
 
std::vector< ESProxyIndex > const & esGetTokenIndicesVector (edm::Transition iTrans) const
 
std::vector< ESRecordIndex > const & esGetTokenRecordIndicesVector (edm::Transition iTrans) const
 
ProductResolverIndexAndSkipBit indexFrom (EDGetToken, BranchType, TypeID const &) const
 
void itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
void itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
std::vector< ProductResolverIndexAndSkipBit > const & itemsToGetFrom (BranchType iType) const
 
void labelsForToken (EDGetToken iToken, Labels &oLabels) const
 
void modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const * > *, NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
 
EDConsumerBase & operator= (EDConsumerBase &&)=default
 
EDConsumerBase const & operator= (EDConsumerBase const &)=delete
 
bool registeredToConsume (ProductResolverIndex, bool, BranchType) const
 
bool registeredToConsumeMany (TypeID const &, BranchType) const
 
void selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)
 
ProductResolverIndexAndSkipBit uncheckedIndexFrom (EDGetToken) const
 
void updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
 
void updateLookup (eventsetup::ESRecordsToProxyIndices const &)
 
virtual ~EDConsumerBase () noexcept(false)
 

Private Member Functions

void analyze (const edm::Event &, const edm::EventSetup &) override
 
void drawHistograms (const TProfile *histo, const TProfile *histoPlusErr, const TProfile *histoMinusErr, const TString &type)
 
void endJob () override
 
void fillHistograms ()
 
void fillValueError ()
 

Private Attributes

bool debug_
 
std::vector< int > errorFactors_
 
std::vector< double > errors_
 
int etaBins_
 
double etaMax_
 
double etaMin_
 
uint32_t maxEvents_
 
TString outputFileName_
 
std::vector< double > parameters_
 
int ptBins_
 
double ptMax_
 
double ptMin_
 
int resolFitType_
 
TProfile * sigmaMassVsEta_
 
TProfile * sigmaMassVsEtaMinusErr_
 
TProfile * sigmaMassVsEtaPlusErr_
 
TProfile * sigmaMassVsPt_
 
TProfile * sigmaMassVsPtMinusErr_
 
TProfile * sigmaMassVsPtPlusErr_
 
TProfile * sigmaPtVsEta_
 
TProfile * sigmaPtVsEtaMinusErr_
 
TProfile * sigmaPtVsEtaPlusErr_
 
TProfile * sigmaPtVsPt_
 
TProfile * sigmaPtVsPtMinusErr_
 
TProfile * sigmaPtVsPtPlusErr_
 
TString treeFileName_
 
std::vector< double > valueMinusError_
 
std::vector< double > valuePlusError_
 

Additional Inherited Members

- Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer ModuleType
 
- Public Types inherited from edm::EDConsumerBase
typedef ProductLabels Labels
 
- Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &)
 
static bool wantsGlobalLuminosityBlocks ()
 
static bool wantsGlobalRuns ()
 
static bool wantsInputProcessBlocks ()
 
static bool wantsProcessBlocks ()
 
static bool wantsStreamLuminosityBlocks ()
 
static bool wantsStreamRuns ()
 
- Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken consumes (const TypeToGet &id, edm::InputTag const &tag)
 
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > consumes (edm::InputTag const &tag)
 
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B > consumes (edm::InputTag tag) noexcept
 
template<BranchType B>
EDGetToken consumes (TypeToGet const &id, edm::InputTag const &tag)
 
ConsumesCollector consumesCollector ()
 Use a ConsumesCollector to gather consumes information from helper functions. More...
 
template<typename ProductType , BranchType B = InEvent>
void consumesMany ()
 
void consumesMany (const TypeToGet &id)
 
template<BranchType B>
void consumesMany (const TypeToGet &id)
 
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes ()
 
template<Transition Tr = Transition::Event>
constexpr auto esConsumes () noexcept
 
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes (ESInputTag const &tag)
 
template<Transition Tr = Transition::Event>
auto esConsumes (ESInputTag tag) noexcept
 
template<Transition Tr = Transition::Event>
ESGetTokenGeneric esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
 Used with EventSetupRecord::doGet. More...
 
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > mayConsume (edm::InputTag const &tag)
 

Detailed Description

Description: <one line="" class="" summary>="">

Implementation: <Notes on="" implementation>="">

Definition at line 48 of file ErrorsAnalyzer.h.

Constructor & Destructor Documentation

◆ ErrorsAnalyzer()

ErrorsAnalyzer::ErrorsAnalyzer ( const edm::ParameterSet &  iConfig)
explicit

Definition at line 6 of file ErrorsAnalyzer.cc.

7  : treeFileName_(iConfig.getParameter<std::string>("InputFileName")),
8  resolFitType_(iConfig.getParameter<int>("ResolFitType")),
9  maxEvents_(iConfig.getParameter<int>("MaxEvents")),
10  outputFileName_(iConfig.getParameter<std::string>("OutputFileName")),
11  ptBins_(iConfig.getParameter<int>("PtBins")),
12  ptMin_(iConfig.getParameter<double>("PtMin")),
13  ptMax_(iConfig.getParameter<double>("PtMax")),
14  etaBins_(iConfig.getParameter<int>("EtaBins")),
15  etaMin_(iConfig.getParameter<double>("EtaMin")),
16  etaMax_(iConfig.getParameter<double>("EtaMax")),
17  debug_(iConfig.getParameter<bool>("Debug")) {
18  parameters_ = iConfig.getParameter<std::vector<double> >("Parameters");
19  errors_ = iConfig.getParameter<std::vector<double> >("Errors");
20  errorFactors_ = iConfig.getParameter<std::vector<int> >("ErrorFactors");
21 
22  if ((parameters_.size() != errors_.size()) || (parameters_.size() != errorFactors_.size())) {
23  std::cout << "Error: parameters and errors have different number of values" << std::endl;
24  exit(1);
25  }
26 
27  fillValueError();
28 
29  sigmaPtVsPt_ = new TProfile("sigmaPtVsPtProfile", "sigmaPtVsPt", ptBins_, ptMin_, ptMax_);
30  sigmaPtVsPtPlusErr_ = new TProfile("sigmaPtVsPtPlusErrProfile", "sigmaPtVsPtPlusErr", ptBins_, ptMin_, ptMax_);
31  sigmaPtVsPtMinusErr_ = new TProfile("sigmaPtVsPtMinusErrProfile", "sigmaPtVsPtMinusErr", ptBins_, ptMin_, ptMax_);
32 
33  sigmaPtVsEta_ = new TProfile("sigmaPtVsEtaProfile", "sigmaPtVsEta", etaBins_, etaMin_, etaMax_);
34  sigmaPtVsEtaPlusErr_ = new TProfile("sigmaPtVsEtaPlusErrProfile", "sigmaPtVsEtaPlusErr", etaBins_, etaMin_, etaMax_);
35  sigmaPtVsEtaMinusErr_ =
36  new TProfile("sigmaPtVsEtaMinusErrProfile", "sigmaPtVsEtaMinusErr", etaBins_, etaMin_, etaMax_);
37 
38  sigmaMassVsPt_ = new TProfile("sigmaMassVsPtProfile", "sigmaMassVsPt", ptBins_, ptMin_, ptMax_);
39  sigmaMassVsPtPlusErr_ = new TProfile("sigmaMassVsPtPlusErrProfile", "sigmaMassVsPtPlusErr", ptBins_, ptMin_, ptMax_);
40  sigmaMassVsPtMinusErr_ =
41  new TProfile("sigmaMassVsPtMinusErrProfile", "sigmaMassVsPtMinusErr", ptBins_, ptMin_, ptMax_);
42 
43  sigmaMassVsEta_ = new TProfile("sigmaMassVsEtaProfile", "sigmaMassVsEta", etaBins_, etaMin_, etaMax_);
44  sigmaMassVsEtaPlusErr_ =
45  new TProfile("sigmaMassVsEtaPlusErrProfile", "sigmaMassVsEtaPlusErr", etaBins_, etaMin_, etaMax_);
46  sigmaMassVsEtaMinusErr_ =
47  new TProfile("sigmaMassVsEtaMinusErrProfile", "sigmaMassVsEtaMinusErr", etaBins_, etaMin_, etaMax_);
48 }

References gather_cfg::cout, errorFactors_, errors_, etaBins_, etaMax_, etaMin_, beamvalidation::exit(), fillValueError(), edm::ParameterSet::getParameter(), parameters_, ptBins_, ptMax_, ptMin_, sigmaMassVsEta_, sigmaMassVsEtaMinusErr_, sigmaMassVsEtaPlusErr_, sigmaMassVsPt_, sigmaMassVsPtMinusErr_, sigmaMassVsPtPlusErr_, sigmaPtVsEta_, sigmaPtVsEtaMinusErr_, sigmaPtVsEtaPlusErr_, sigmaPtVsPt_, sigmaPtVsPtMinusErr_, and sigmaPtVsPtPlusErr_.

◆ ~ErrorsAnalyzer()

ErrorsAnalyzer::~ErrorsAnalyzer ( )
override

Definition at line 64 of file ErrorsAnalyzer.cc.

64  {
65  gROOT->SetStyle("Plain");
66 
67  fillHistograms();
68 
69  TFile* outputFile = new TFile(outputFileName_, "RECREATE");
70 
71  drawHistograms(sigmaPtVsEta_, sigmaPtVsEtaPlusErr_, sigmaPtVsEtaMinusErr_, "sigmaPtVsEta");
72  drawHistograms(sigmaPtVsPt_, sigmaPtVsPtPlusErr_, sigmaPtVsPtMinusErr_, "sigmaPtVsPt");
73 
74  drawHistograms(sigmaMassVsEta_, sigmaMassVsEtaPlusErr_, sigmaMassVsEtaMinusErr_, "sigmaMassVsEta");
75  drawHistograms(sigmaMassVsPt_, sigmaMassVsPtPlusErr_, sigmaMassVsPtMinusErr_, "sigmaMassVsPt");
76 
77  outputFile->Write();
78  outputFile->Close();
79 }

References drawHistograms(), fillHistograms(), download_sqlite_cfg::outputFile, outputFileName_, sigmaMassVsEta_, sigmaMassVsEtaMinusErr_, sigmaMassVsEtaPlusErr_, sigmaMassVsPt_, sigmaMassVsPtMinusErr_, sigmaMassVsPtPlusErr_, sigmaPtVsEta_, sigmaPtVsEtaMinusErr_, sigmaPtVsEtaPlusErr_, sigmaPtVsPt_, sigmaPtVsPtMinusErr_, and sigmaPtVsPtPlusErr_.

Member Function Documentation

◆ analyze()

void ErrorsAnalyzer::analyze ( const edm::Event &  iEvent,
const edm::EventSetup &  iSetup 
)
overrideprivatevirtual

Implements edm::EDAnalyzer.

Definition at line 194 of file ErrorsAnalyzer.cc.

194 {}

◆ drawHistograms()

void ErrorsAnalyzer::drawHistograms ( const TProfile *  histo,
const TProfile *  histoPlusErr,
const TProfile *  histoMinusErr,
const TString &  type 
)
private

Definition at line 81 of file ErrorsAnalyzer.cc.

84  {
85  TH1D* sigmaPtVsEtaTH1D =
86  new TH1D(type, type, histo->GetNbinsX(), histo->GetXaxis()->GetXmin(), histo->GetXaxis()->GetXmax());
87  TH1D* sigmaPtVsEtaPlusErrTH1D = new TH1D(type + "PlusErr",
88  type + "PlusErr",
89  histo->GetNbinsX(),
90  histo->GetXaxis()->GetXmin(),
91  histo->GetXaxis()->GetXmax());
92  TH1D* sigmaPtVsEtaMinusErrTH1D = new TH1D(type + "MinusErr",
93  type + "MinusErr",
94  histo->GetNbinsX(),
95  histo->GetXaxis()->GetXmin(),
96  histo->GetXaxis()->GetXmax());
97 
98  TH1D* sigmaMassVsEtaTH1D =
99  new TH1D(type, type, histo->GetNbinsX(), histo->GetXaxis()->GetXmin(), histo->GetXaxis()->GetXmax());
100  TH1D* sigmaMassVsEtaPlusErrTH1D = new TH1D(type + "PlusErr",
101  type + "PlusErr",
102  histo->GetNbinsX(),
103  histo->GetXaxis()->GetXmin(),
104  histo->GetXaxis()->GetXmax());
105  TH1D* sigmaMassVsEtaMinusErrTH1D = new TH1D(type + "MinusErr",
106  type + "MinusErr",
107  histo->GetNbinsX(),
108  histo->GetXaxis()->GetXmin(),
109  histo->GetXaxis()->GetXmax());
110 
111  TCanvas* canvas = new TCanvas("canvas" + type, "canvas" + type, 1000, 800);
112  for (int iBin = 1; iBin <= histo->GetNbinsX(); ++iBin) {
113  sigmaPtVsEtaTH1D->SetBinContent(iBin, histo->GetBinContent(iBin));
114  sigmaPtVsEtaPlusErrTH1D->SetBinContent(iBin, histoPlusErr->GetBinContent(iBin));
115  sigmaPtVsEtaMinusErrTH1D->SetBinContent(iBin, histoMinusErr->GetBinContent(iBin));
116  }
117  int numBins = sigmaPtVsEtaTH1D->GetNbinsX();
118  // Draw TGraph with asymmetric errors
119  Double_t* values = sigmaPtVsEtaTH1D->GetArray();
120  Double_t* valuesPlus = sigmaPtVsEtaPlusErrTH1D->GetArray();
121  Double_t* valuesMinus = sigmaPtVsEtaMinusErrTH1D->GetArray();
122  double* posErrors = new double[numBins];
123  double* negErrors = new double[numBins];
124 
125  TGraphAsymmErrors* graphAsymmErrors = new TGraphAsymmErrors(sigmaPtVsEtaTH1D);
126  TGraph* graph = new TGraph(sigmaPtVsEtaTH1D);
127 
128  for (int i = 1; i <= numBins; ++i) {
129  // std::cout << "filling " << i << std::endl;
130  posErrors[i - 1] = valuesPlus[i] - values[i];
131  if (valuesMinus[i] < 0)
132  negErrors[i - 1] = values[i];
133  else
134  negErrors[i - 1] = values[i] - valuesMinus[i];
135 
136  graphAsymmErrors->SetPointEYlow(i - 1, negErrors[i - 1]);
137  graphAsymmErrors->SetPointEYhigh(i - 1, posErrors[i - 1]);
138  }
139 
140  canvas->Draw();
141 
142  delete[] posErrors;
143  delete[] negErrors;
144 
145  graphAsymmErrors->SetTitle("");
146  graphAsymmErrors->SetFillColor(kGray);
147  graphAsymmErrors->Draw("A2");
148  TString title(type);
149  if (type == "Eta")
150  title = "#eta";
151  graphAsymmErrors->GetXaxis()->SetTitle(title);
152  graphAsymmErrors->GetYaxis()->SetTitle("#sigmaPt/Pt");
153  graph->Draw();
154 
155  // graph->SetLineColor(kGray);
156  // graph->SetMarkerColor(kBlack);
157  // graph->Draw("AP");
158 
159  // if( debug_ ) {
160  // sigmaPtVsEtaPlusErrTH1D->SetLineColor(kRed);
161  // sigmaPtVsEtaMinusErrTH1D->SetLineColor(kBlue);
162  // }
163  // else {
164  // sigmaPtVsEtaPlusErrTH1D->SetFillColor(kGray);
165  // sigmaPtVsEtaPlusErrTH1D->SetLineColor(kWhite);
166  // sigmaPtVsEtaMinusErrTH1D->SetFillColor(kWhite);
167  // sigmaPtVsEtaMinusErrTH1D->SetLineColor(kWhite);
168  // }
169  // sigmaPtVsEtaPlusErrTH1D->Draw();
170  // sigmaPtVsEtaTH1D->Draw("SAME");
171  // sigmaPtVsEtaMinusErrTH1D->Draw("SAME");
172 
173  sigmaPtVsEtaPlusErrTH1D->Write();
174  sigmaPtVsEtaTH1D->Write();
175  sigmaPtVsEtaMinusErrTH1D->Write();
176 
177  sigmaPtVsEtaPlusErr_->Write();
178  sigmaPtVsEta_->Write();
179  sigmaPtVsEtaMinusErr_->Write();
180 
181  // Mass
182  sigmaMassVsEtaPlusErrTH1D->Write();
183  sigmaMassVsEtaTH1D->Write();
184  sigmaMassVsEtaMinusErrTH1D->Write();
185 
186  sigmaMassVsEtaPlusErr_->Write();
187  sigmaMassVsEta_->Write();
188  sigmaMassVsEtaMinusErr_->Write();
189 
190  canvas->Write();
191 }

References svgfig::canvas(), timingPdfMaker::histo, mps_fire::i, sigmaMassVsEta_, sigmaMassVsEtaMinusErr_, sigmaMassVsEtaPlusErr_, sigmaPtVsEta_, sigmaPtVsEtaMinusErr_, sigmaPtVsEtaPlusErr_, runGCPTkAlMap::title, and contentValuesCheck::values.

Referenced by ~ErrorsAnalyzer().

◆ endJob()

void ErrorsAnalyzer::endJob ( void  )
inlineoverrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 61 of file ErrorsAnalyzer.h.

61 {};

◆ fillHistograms()

void ErrorsAnalyzer::fillHistograms ( )
private

Definition at line 196 of file ErrorsAnalyzer.cc.

196  {
197  std::cout << "Reading muon pairs from Root Tree in " << treeFileName_ << std::endl;
198  RootTreeHandler rootTreeHandler;
199 
200  typedef std::vector<std::pair<lorentzVector, lorentzVector> > MuonPairVector;
201  MuonPairVector savedPair;
202  std::vector<std::pair<unsigned int, unsigned long long> > evtRun;
203  rootTreeHandler.readTree(maxEvents_, treeFileName_, &savedPair, 0, &evtRun);
204  // rootTreeHandler.readTree(maxEvents, inputRootTreeFileName_, &savedPair, &(MuScleFitUtils::genPair));
205 
206  resolutionFunctionBase<std::vector<double> >* resolutionFunctionForVec = resolutionFunctionVecService(resolFitType_);
207  MuScleFitUtils::resolutionFunction = resolutionFunctionService(resolFitType_);
208  MuScleFitUtils::debugMassResol_ = false;
209  MuScleFitUtils::debug = 0;
210  MuScleFitUtils::resfind = std::vector<int>(6, 0);
211 
212  // Loop on all the pairs
213  unsigned int i = 0;
214  MuonPairVector::iterator it = savedPair.begin();
215  std::cout << "Starting loop on " << savedPair.size() << " muons" << std::endl;
216  for (; it != savedPair.end(); ++it, ++i) {
217  double pt1 = it->first.pt();
218  double eta1 = it->first.eta();
219  double pt2 = it->second.pt();
220  double eta2 = it->second.eta();
221 
222  if (debug_) {
223  std::cout << "pt1 = " << pt1 << ", eta1 = " << eta1 << ", pt2 = " << pt2 << ", eta2 = " << eta2 << std::endl;
224  }
225 
226  if (pt1 == 0 && pt2 == 0 && eta1 == 0 && eta2 == 0)
227  continue;
228 
229  double sigmaPt1 = resolutionFunctionForVec->sigmaPt(pt1, eta1, parameters_);
230  double sigmaPt2 = resolutionFunctionForVec->sigmaPt(pt2, eta2, parameters_);
231  double sigmaPtPlusErr1 = resolutionFunctionForVec->sigmaPt(pt1, eta1, valuePlusError_);
232  double sigmaPtPlusErr2 = resolutionFunctionForVec->sigmaPt(pt2, eta2, valuePlusError_);
233  double sigmaPtMinusErr1 = resolutionFunctionForVec->sigmaPt(pt1, eta1, valueMinusError_);
234  double sigmaPtMinusErr2 = resolutionFunctionForVec->sigmaPt(pt2, eta2, valueMinusError_);
235 
236  double sigmaMass = MuScleFitUtils::massResolution(it->first, it->second, parameters_);
237  double sigmaMassPlusErr = MuScleFitUtils::massResolution(it->first, it->second, valuePlusError_);
238  double sigmaMassMinusErr = MuScleFitUtils::massResolution(it->first, it->second, valueMinusError_);
239 
240  if (debug_) {
241  std::cout << "sigmaPt1 = " << sigmaPt1 << " + " << sigmaPtPlusErr1 << " - " << sigmaPtMinusErr1 << std::endl;
242  std::cout << "sigmaPt2 = " << sigmaPt2 << " + " << sigmaPtPlusErr2 << " - " << sigmaPtMinusErr2 << std::endl;
243  std::cout << "sigmaMass = " << sigmaMass << " + " << sigmaMassPlusErr << " - " << sigmaMassMinusErr << std::endl;
244  }
245 
246  // Protections from nans
247  if (pt1 != pt1)
248  continue;
249  if (pt2 != pt2)
250  continue;
251  if (eta1 != eta1)
252  continue;
253  if (eta2 != eta2)
254  continue;
255  if (sigmaPt1 != sigmaPt1)
256  continue;
257  if (sigmaPt2 != sigmaPt2)
258  continue;
259  if (sigmaPtPlusErr1 != sigmaPtPlusErr1)
260  continue;
261  if (sigmaPtPlusErr2 != sigmaPtPlusErr2)
262  continue;
263  if (sigmaPtMinusErr1 != sigmaPtMinusErr1)
264  continue;
265  if (sigmaPtMinusErr2 != sigmaPtMinusErr2)
266  continue;
267  if (sigmaMass != sigmaMass)
268  continue;
269  if (sigmaMassPlusErr != sigmaMassPlusErr)
270  continue;
271  if (sigmaMassMinusErr != sigmaMassMinusErr)
272  continue;
273 
274  std::cout << "Muon pair number " << i << std::endl;
275 
276  // std::cout << "sigmaPt1 = " << sigmaPt1 << ", sigmaPt2 = " << sigmaPt2 << std::endl;
277  // std::cout << "sigmaPtPlusErr1 = " << sigmaPtPlusErr1 << ", sigmaPtMinusErr2 = " << sigmaPtMinusErr2 << std::endl;
278  // std::cout << "sigmaPtMinusErr1 = " << sigmaPtPlusErr1 << ", sigmaPtMinusErr2 = " << sigmaPtMinusErr2 << std::endl;
279 
280  sigmaPtVsPt_->Fill(pt1, sigmaPt1);
281  sigmaPtVsPt_->Fill(pt2, sigmaPt2);
282  sigmaPtVsPtPlusErr_->Fill(pt1, sigmaPtPlusErr1);
283  sigmaPtVsPtPlusErr_->Fill(pt2, sigmaPtPlusErr2);
284  sigmaPtVsPtMinusErr_->Fill(pt1, sigmaPtMinusErr1);
285  sigmaPtVsPtMinusErr_->Fill(pt2, sigmaPtMinusErr2);
286 
287  sigmaPtVsEta_->Fill(eta1, sigmaPt1);
288  sigmaPtVsEta_->Fill(eta2, sigmaPt2);
289  sigmaPtVsEtaPlusErr_->Fill(eta1, sigmaPtPlusErr1);
290  sigmaPtVsEtaPlusErr_->Fill(eta2, sigmaPtPlusErr2);
291  sigmaPtVsEtaMinusErr_->Fill(eta1, sigmaPtMinusErr1);
292  sigmaPtVsEtaMinusErr_->Fill(eta2, sigmaPtMinusErr2);
293 
294  sigmaMassVsPt_->Fill(pt1, sigmaMass);
295  sigmaMassVsPt_->Fill(pt2, sigmaMass);
296  sigmaMassVsPtPlusErr_->Fill(pt1, sigmaMassPlusErr);
297  sigmaMassVsPtPlusErr_->Fill(pt2, sigmaMassPlusErr);
298  sigmaMassVsPtMinusErr_->Fill(pt1, sigmaMassMinusErr);
299  sigmaMassVsPtMinusErr_->Fill(pt2, sigmaMassMinusErr);
300 
301  sigmaMassVsEta_->Fill(eta1, sigmaMass);
302  sigmaMassVsEta_->Fill(eta2, sigmaMass);
303  sigmaMassVsEtaPlusErr_->Fill(eta1, sigmaMassPlusErr);
304  sigmaMassVsEtaPlusErr_->Fill(eta2, sigmaMassPlusErr);
305  sigmaMassVsEtaMinusErr_->Fill(eta1, sigmaMassMinusErr);
306  sigmaMassVsEtaMinusErr_->Fill(eta2, sigmaMassMinusErr);
307  }
308 }

References gather_cfg::cout, MuScleFitUtils::debug, debug_, MuScleFitUtils::debugMassResol_, HLT_FULL_cff::eta1, HLT_FULL_cff::eta2, mps_fire::i, MuScleFitUtils::massResolution(), maxEvents_, parameters_, HLT_FULL_cff::pt1, HLT_FULL_cff::pt2, RootTreeHandler::readTree(), MuScleFitUtils::resfind, resolFitType_, MuScleFitUtils::resolutionFunction, resolutionFunctionService(), resolutionFunctionVecService(), sigmaMassVsEta_, sigmaMassVsEtaMinusErr_, sigmaMassVsEtaPlusErr_, sigmaMassVsPt_, sigmaMassVsPtMinusErr_, sigmaMassVsPtPlusErr_, resolutionFunctionBase< T >::sigmaPt(), sigmaPtVsEta_, sigmaPtVsEtaMinusErr_, sigmaPtVsEtaPlusErr_, sigmaPtVsPt_, sigmaPtVsPtMinusErr_, sigmaPtVsPtPlusErr_, treeFileName_, valueMinusError_, and valuePlusError_.

Referenced by ~ErrorsAnalyzer().

◆ fillValueError()

void ErrorsAnalyzer::fillValueError ( )
private

Definition at line 50 of file ErrorsAnalyzer.cc.

50  {
51  valuePlusError_.resize(parameters_.size());
52  valueMinusError_.resize(parameters_.size());
53 
54  std::vector<double>::const_iterator parIt = parameters_.begin();
55  std::vector<double>::const_iterator errIt = errors_.begin();
56  std::vector<int>::const_iterator errFactorIt = errorFactors_.begin();
57  int i = 0;
58  for (; parIt != parameters_.end(); ++parIt, ++errIt, ++errFactorIt, ++i) {
59  valuePlusError_[i] = *parIt + (*errIt) * (*errFactorIt);
60  valueMinusError_[i] = *parIt - (*errIt) * (*errFactorIt);
61  }
62 }

References errorFactors_, errors_, mps_fire::i, parameters_, valueMinusError_, and valuePlusError_.

Referenced by ErrorsAnalyzer().

Member Data Documentation

◆ debug_

bool ErrorsAnalyzer::debug_
private

Definition at line 73 of file ErrorsAnalyzer.h.

Referenced by fillHistograms().

◆ errorFactors_

std::vector<int> ErrorsAnalyzer::errorFactors_
private

Definition at line 77 of file ErrorsAnalyzer.h.

Referenced by ErrorsAnalyzer(), and fillValueError().

◆ errors_

std::vector<double> ErrorsAnalyzer::errors_
private

Definition at line 76 of file ErrorsAnalyzer.h.

Referenced by ErrorsAnalyzer(), and fillValueError().

◆ etaBins_

int ErrorsAnalyzer::etaBins_
private

Definition at line 70 of file ErrorsAnalyzer.h.

Referenced by ErrorsAnalyzer().

◆ etaMax_

double ErrorsAnalyzer::etaMax_
private

Definition at line 72 of file ErrorsAnalyzer.h.

Referenced by ErrorsAnalyzer().

◆ etaMin_

double ErrorsAnalyzer::etaMin_
private

Definition at line 71 of file ErrorsAnalyzer.h.

Referenced by ErrorsAnalyzer().

◆ maxEvents_

uint32_t ErrorsAnalyzer::maxEvents_
private

Definition at line 65 of file ErrorsAnalyzer.h.

Referenced by fillHistograms().

◆ outputFileName_

TString ErrorsAnalyzer::outputFileName_
private

Definition at line 66 of file ErrorsAnalyzer.h.

Referenced by ~ErrorsAnalyzer().

◆ parameters_

std::vector<double> ErrorsAnalyzer::parameters_
private

Definition at line 75 of file ErrorsAnalyzer.h.

Referenced by Mixins._TypedParameterizable::clone(), Types.PSet::clone(), Types.EDAlias::clone(), Mixins._TypedParameterizable::copy(), ErrorsAnalyzer(), fillHistograms(), and fillValueError().

◆ ptBins_

int ErrorsAnalyzer::ptBins_
private

Definition at line 67 of file ErrorsAnalyzer.h.

Referenced by ErrorsAnalyzer().

◆ ptMax_

double ErrorsAnalyzer::ptMax_
private

Definition at line 69 of file ErrorsAnalyzer.h.

Referenced by ErrorsAnalyzer().

◆ ptMin_

double ErrorsAnalyzer::ptMin_
private

Definition at line 68 of file ErrorsAnalyzer.h.

Referenced by ErrorsAnalyzer().

◆ resolFitType_

int ErrorsAnalyzer::resolFitType_
private

Definition at line 64 of file ErrorsAnalyzer.h.

Referenced by fillHistograms().

◆ sigmaMassVsEta_

TProfile* ErrorsAnalyzer::sigmaMassVsEta_
private

Definition at line 91 of file ErrorsAnalyzer.h.

Referenced by drawHistograms(), ErrorsAnalyzer(), fillHistograms(), and ~ErrorsAnalyzer().

◆ sigmaMassVsEtaMinusErr_

TProfile* ErrorsAnalyzer::sigmaMassVsEtaMinusErr_
private

Definition at line 93 of file ErrorsAnalyzer.h.

Referenced by drawHistograms(), ErrorsAnalyzer(), fillHistograms(), and ~ErrorsAnalyzer().

◆ sigmaMassVsEtaPlusErr_

TProfile* ErrorsAnalyzer::sigmaMassVsEtaPlusErr_
private

Definition at line 92 of file ErrorsAnalyzer.h.

Referenced by drawHistograms(), ErrorsAnalyzer(), fillHistograms(), and ~ErrorsAnalyzer().

◆ sigmaMassVsPt_

TProfile* ErrorsAnalyzer::sigmaMassVsPt_
private

Definition at line 95 of file ErrorsAnalyzer.h.

Referenced by ErrorsAnalyzer(), fillHistograms(), and ~ErrorsAnalyzer().

◆ sigmaMassVsPtMinusErr_

TProfile* ErrorsAnalyzer::sigmaMassVsPtMinusErr_
private

Definition at line 97 of file ErrorsAnalyzer.h.

Referenced by ErrorsAnalyzer(), fillHistograms(), and ~ErrorsAnalyzer().

◆ sigmaMassVsPtPlusErr_

TProfile* ErrorsAnalyzer::sigmaMassVsPtPlusErr_
private

Definition at line 96 of file ErrorsAnalyzer.h.

Referenced by ErrorsAnalyzer(), fillHistograms(), and ~ErrorsAnalyzer().

◆ sigmaPtVsEta_

TProfile* ErrorsAnalyzer::sigmaPtVsEta_
private

Definition at line 82 of file ErrorsAnalyzer.h.

Referenced by drawHistograms(), ErrorsAnalyzer(), fillHistograms(), and ~ErrorsAnalyzer().

◆ sigmaPtVsEtaMinusErr_

TProfile* ErrorsAnalyzer::sigmaPtVsEtaMinusErr_
private

Definition at line 84 of file ErrorsAnalyzer.h.

Referenced by drawHistograms(), ErrorsAnalyzer(), fillHistograms(), and ~ErrorsAnalyzer().

◆ sigmaPtVsEtaPlusErr_

TProfile* ErrorsAnalyzer::sigmaPtVsEtaPlusErr_
private

Definition at line 83 of file ErrorsAnalyzer.h.

Referenced by drawHistograms(), ErrorsAnalyzer(), fillHistograms(), and ~ErrorsAnalyzer().

◆ sigmaPtVsPt_

TProfile* ErrorsAnalyzer::sigmaPtVsPt_
private

Definition at line 86 of file ErrorsAnalyzer.h.

Referenced by ErrorsAnalyzer(), fillHistograms(), and ~ErrorsAnalyzer().

◆ sigmaPtVsPtMinusErr_

TProfile* ErrorsAnalyzer::sigmaPtVsPtMinusErr_
private

Definition at line 88 of file ErrorsAnalyzer.h.

Referenced by ErrorsAnalyzer(), fillHistograms(), and ~ErrorsAnalyzer().

◆ sigmaPtVsPtPlusErr_

TProfile* ErrorsAnalyzer::sigmaPtVsPtPlusErr_
private

Definition at line 87 of file ErrorsAnalyzer.h.

Referenced by ErrorsAnalyzer(), fillHistograms(), and ~ErrorsAnalyzer().

◆ treeFileName_

TString ErrorsAnalyzer::treeFileName_
private

Definition at line 61 of file ErrorsAnalyzer.h.

Referenced by fillHistograms().

◆ valueMinusError_

std::vector<double> ErrorsAnalyzer::valueMinusError_
private

Definition at line 80 of file ErrorsAnalyzer.h.

Referenced by fillHistograms(), and fillValueError().

◆ valuePlusError_

std::vector<double> ErrorsAnalyzer::valuePlusError_
private

Definition at line 79 of file ErrorsAnalyzer.h.

Referenced by fillHistograms(), and fillValueError().

ErrorsAnalyzer::ptMin_
double ptMin_
Definition: ErrorsAnalyzer.h:68
svgfig.canvas
def canvas(*sub, **attr)
Definition: svgfig.py:482
ErrorsAnalyzer::errors_
std::vector< double > errors_
Definition: ErrorsAnalyzer.h:76
runGCPTkAlMap.title
string title
Definition: runGCPTkAlMap.py:94
ErrorsAnalyzer::outputFileName_
TString outputFileName_
Definition: ErrorsAnalyzer.h:66
mps_fire.i
i
Definition: mps_fire.py:428
RootTreeHandler
Definition: RootTreeHandler.h:24
resolutionFunctionBase
Definition: Functions.h:828
ErrorsAnalyzer::sigmaMassVsPtPlusErr_
TProfile * sigmaMassVsPtPlusErr_
Definition: ErrorsAnalyzer.h:96
ErrorsAnalyzer::etaMin_
double etaMin_
Definition: ErrorsAnalyzer.h:71
gather_cfg.cout
cout
Definition: gather_cfg.py:144
ErrorsAnalyzer::valueMinusError_
std::vector< double > valueMinusError_
Definition: ErrorsAnalyzer.h:80
timingPdfMaker.histo
histo
Definition: timingPdfMaker.py:278
resolutionFunctionService
resolutionFunctionBase< double * > * resolutionFunctionService(const int identifier)
Service to build the resolution functor corresponding to the passed identifier.
Definition: Functions.cc:70
ErrorsAnalyzer::sigmaMassVsEta_
TProfile * sigmaMassVsEta_
Definition: ErrorsAnalyzer.h:91
MuScleFitUtils::resfind
static std::vector< int > resfind
Definition: MuScleFitUtils.h:210
ErrorsAnalyzer::drawHistograms
void drawHistograms(const TProfile *histo, const TProfile *histoPlusErr, const TProfile *histoMinusErr, const TString &type)
Definition: ErrorsAnalyzer.cc:81
ErrorsAnalyzer::resolFitType_
int resolFitType_
Definition: ErrorsAnalyzer.h:64
resolutionFunctionVecService
resolutionFunctionBase< std::vector< double > > * resolutionFunctionVecService(const int identifier)
Service to build the resolution functor corresponding to the passed identifier when receiving a std::...
Definition: Functions.cc:92
ErrorsAnalyzer::sigmaMassVsPtMinusErr_
TProfile * sigmaMassVsPtMinusErr_
Definition: ErrorsAnalyzer.h:97
MuonPairVector
std::vector< std::pair< lorentzVector, lorentzVector > > MuonPairVector
Definition: RootTreeHandler.h:13
download_sqlite_cfg.outputFile
outputFile
Definition: download_sqlite_cfg.py:5
ErrorsAnalyzer::etaBins_
int etaBins_
Definition: ErrorsAnalyzer.h:70
MuScleFitUtils::debugMassResol_
static bool debugMassResol_
Definition: MuScleFitUtils.h:268
ErrorsAnalyzer::sigmaMassVsEtaPlusErr_
TProfile * sigmaMassVsEtaPlusErr_
Definition: ErrorsAnalyzer.h:92
contentValuesCheck.values
values
Definition: contentValuesCheck.py:38
HLT_FULL_cff.eta2
eta2
Definition: HLT_FULL_cff.py:9551
ErrorsAnalyzer::errorFactors_
std::vector< int > errorFactors_
Definition: ErrorsAnalyzer.h:77
HLT_FULL_cff.pt1
pt1
Definition: HLT_FULL_cff.py:9893
ErrorsAnalyzer::maxEvents_
uint32_t maxEvents_
Definition: ErrorsAnalyzer.h:65
ErrorsAnalyzer::valuePlusError_
std::vector< double > valuePlusError_
Definition: ErrorsAnalyzer.h:79
MuScleFitUtils::massResolution
static double massResolution(const lorentzVector &mu1, const lorentzVector &mu2)
ErrorsAnalyzer::sigmaPtVsEtaMinusErr_
TProfile * sigmaPtVsEtaMinusErr_
Definition: ErrorsAnalyzer.h:84
HLT_FULL_cff.eta1
eta1
Definition: HLT_FULL_cff.py:9550
MuScleFitUtils::debug
static int debug
Definition: MuScleFitUtils.h:130
ErrorsAnalyzer::fillHistograms
void fillHistograms()
Definition: ErrorsAnalyzer.cc:196
type
type
Definition: SiPixelVCal_PayloadInspector.cc:39
ErrorsAnalyzer::sigmaPtVsEta_
TProfile * sigmaPtVsEta_
Definition: ErrorsAnalyzer.h:82
resolutionFunctionBase::sigmaPt
virtual double sigmaPt(const double &pt, const double &eta, const T &parval)=0
ErrorsAnalyzer::sigmaPtVsEtaPlusErr_
TProfile * sigmaPtVsEtaPlusErr_
Definition: ErrorsAnalyzer.h:83
ErrorsAnalyzer::parameters_
std::vector< double > parameters_
Definition: ErrorsAnalyzer.h:75
HLT_FULL_cff.pt2
pt2
Definition: HLT_FULL_cff.py:9895
ErrorsAnalyzer::etaMax_
double etaMax_
Definition: ErrorsAnalyzer.h:72
ErrorsAnalyzer::sigmaMassVsEtaMinusErr_
TProfile * sigmaMassVsEtaMinusErr_
Definition: ErrorsAnalyzer.h:93
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
ErrorsAnalyzer::fillValueError
void fillValueError()
Definition: ErrorsAnalyzer.cc:50
ErrorsAnalyzer::sigmaMassVsPt_
TProfile * sigmaMassVsPt_
Definition: ErrorsAnalyzer.h:95
ErrorsAnalyzer::debug_
bool debug_
Definition: ErrorsAnalyzer.h:73
ErrorsAnalyzer::ptMax_
double ptMax_
Definition: ErrorsAnalyzer.h:69
ErrorsAnalyzer::treeFileName_
TString treeFileName_
Definition: ErrorsAnalyzer.h:61
ErrorsAnalyzer::sigmaPtVsPtMinusErr_
TProfile * sigmaPtVsPtMinusErr_
Definition: ErrorsAnalyzer.h:88
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
ErrorsAnalyzer::sigmaPtVsPtPlusErr_
TProfile * sigmaPtVsPtPlusErr_
Definition: ErrorsAnalyzer.h:87
RootTreeHandler::readTree
void readTree(const int maxEvents, const TString &fileName, MuonPairVector *savedPair, const int muonType, std::vector< std::pair< unsigned int, unsigned long long > > *evtRun, MuonPairVector *genPair=nullptr)
Definition: RootTreeHandler.h:86
ErrorsAnalyzer::sigmaPtVsPt_
TProfile * sigmaPtVsPt_
Definition: ErrorsAnalyzer.h:86
beamvalidation.exit
def exit(msg="")
Definition: beamvalidation.py:52
MuScleFitUtils::resolutionFunction
static resolutionFunctionBase< double * > * resolutionFunction
Definition: MuScleFitUtils.h:153
ErrorsAnalyzer::ptBins_
int ptBins_
Definition: ErrorsAnalyzer.h:67
Generated for CMSSW Reference Manual by doxygen 1.8.16