ErrorsPropagationAnalyzer.h

Go to the documentation of this file.

#ifndef ERRORSANALYZER_HH
#define ERRORSANALYZER_HH
 
// -*- C++ -*-
//
// Package:    ErrorsPropagationAnalyzer
// Class:      ErrorsPropagationAnalyzer
//
//
// Original Author:  Marco De Mattia
//         Created:  Thu Sep 11 12:16:00 CEST 2008
//
//
 
// system include files
#include <memory>
#include <string>
#include <vector>
 
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/EDAnalyzer.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
 
#include <TH1D.h>
#include <TProfile.h>
#include <TString.h>
#include <TCanvas.h>
#include <TGraphAsymmErrors.h>
#include <TROOT.h>
 
#include "MuonAnalysis/MomentumScaleCalibration/interface/Functions.h"
#include "MuonAnalysis/MomentumScaleCalibration/interface/RootTreeHandler.h"
#include "MuScleFitUtils.h"
#include "MuonAnalysis/MomentumScaleCalibration/interface/SigmaPtDiff.h"
 
//
// class declaration
//
 
class ErrorsPropagationAnalyzer : public edm::EDAnalyzer {
public:
  explicit ErrorsPropagationAnalyzer(const edm::ParameterSet&);
  ~ErrorsPropagationAnalyzer() override;
 
private:
  void analyze(const edm::Event&, const edm::EventSetup&) override;
  void fillHistograms();
  void drawHistograms(const TProfile* histo,
                      const TProfile* histoPlusErr,
                      const TProfile* histoMinusErr,
                      const TString& type,
                      const TString& yLabel);
  void fillValueError();
  void endJob() override{};
  double massResolution(const lorentzVector& mu1,
                        const lorentzVector& mu2,
                        const std::vector<double>& parval,
                        const double& sigmaPt1,
                        const double& sigmaPt2);
  double massResolution(const lorentzVector& mu1,
                        const lorentzVector& mu2,
                        double* parval,
                        const double& sigmaPt1,
                        const double& sigmaPt2);
 
  TString treeFileName_;
  int resolFitType_;
  uint32_t maxEvents_;
  TString outputFileName_;
  int ptBins_;
  double ptMin_;
  double ptMax_;
  int etaBins_;
  double etaMin_;
  double etaMax_;
  bool debug_;
 
  double ptMinCut_, ptMaxCut_, etaMinCut_, etaMaxCut_;
 
  std::vector<double> parameters_;
  std::vector<double> errors_;
  std::vector<int> errorFactors_;
 
  std::vector<double> valuePlusError_;
  std::vector<double> valueMinusError_;
 
  TProfile* sigmaPtVsEta_;
  TProfile* sigmaPtVsEtaPlusErr_;
  TProfile* sigmaPtVsEtaMinusErr_;
 
  TProfile* sigmaPtVsPt_;
  TProfile* sigmaPtVsPtPlusErr_;
  TProfile* sigmaPtVsPtMinusErr_;
 
  TProfile* sigmaPtVsEtaDiff_;
  TProfile* sigmaPtVsPtDiff_;
 
  // Mass resolution
  TProfile* sigmaMassVsEta_;
  TProfile* sigmaMassVsEtaPlusErr_;
  TProfile* sigmaMassVsEtaMinusErr_;
 
  TProfile* sigmaMassVsPt_;
  TProfile* sigmaMassVsPtPlusErr_;
  TProfile* sigmaMassVsPtMinusErr_;
 
  TProfile* sigmaMassOverMassVsEta_;
  TProfile* sigmaMassOverMassVsEtaPlusErr_;
  TProfile* sigmaMassOverMassVsEtaMinusErr_;
 
  TProfile* sigmaMassOverMassVsPt_;
  TProfile* sigmaMassOverMassVsPtPlusErr_;
  TProfile* sigmaMassOverMassVsPtMinusErr_;
};
 
#endif  // RESOLUTIONANALYZER_HH