MuScleFitBase Class Reference

#include <MuScleFitBase.h>

Inheritance diagram for MuScleFitBase:

Classes
class	ProbForIntegral
	Functor used to compute the normalization integral of probability functions. More...

Public Member Functions
	MuScleFitBase (const edm::ParameterSet &iConfig)
virtual	~MuScleFitBase ()(false)

Protected Member Functions
void	clearHistoMap ()
	Clean the histograms map. More...
void	fillHistoMap (TFile *outputFile, unsigned int iLoop)
	Create the histograms map. More...
void	readProbabilityDistributionsFromFile ()
	Read probability distributions from a local root file. More...
void	writeHistoMap (const unsigned int iLoop)
	Save the histograms map to file. More...

Protected Attributes
int	debug_
std::vector< GenMuonPair >	genMuonPairs_
	Stores the genMuon pairs and the motherId prior to the creation of the internal tree. More...
std::map< std::string, Histograms * >	mapHisto_
	The map of histograms. More...
std::vector< MuonPair >	muonPairs_
	Used to store the muon pairs plus run and event number prior to the creation of the internal tree. More...
std::string	probabilitiesFile_
std::string	probabilitiesFileInPath_
std::vector< TFile * >	theFiles_
	The files were the histograms are saved. More...
std::string	theGenInfoRootFileName_
edm::InputTag	theMuonLabel_
int	theMuonType_
std::string	theRootFileName_

Detailed Description

This class is used as a base for MuSlceFit. The reason for putting some of the methods inside this base class is that they are used also by the TestCorrection analyzer.

Definition at line 18 of file MuScleFitBase.h.

Constructor & Destructor Documentation

MuScleFitBase::MuScleFitBase ( const edm::ParameterSet &  iConfig )

inline

Definition at line 20 of file MuScleFitBase.h.

      : probabilitiesFileInPath_(iConfig.getUntrackedParameter<std::string>(
            "ProbabilitiesFileInPath", "MuonAnalysis/MomentumScaleCalibration/test/Probs_new_Horace_CTEQ_1000.root")),
        probabilitiesFile_(iConfig.getUntrackedParameter<std::string>("ProbabilitiesFile", "")),
        theMuonType_(iConfig.getParameter<int>("MuonType")),
        theMuonLabel_(iConfig.getParameter<edm::InputTag>("MuonLabel")),
        theRootFileName_(iConfig.getUntrackedParameter<std::string>("OutputFileName")),
        theGenInfoRootFileName_(
            iConfig.getUntrackedParameter<std::string>("OutputGenInfoFileName", "genSimRecoPlots.root")),
        debug_(iConfig.getUntrackedParameter<int>("debug", 0)) {}

virtual MuScleFitBase::~MuScleFitBase ( )

inlinevirtual

Definition at line 30 of file MuScleFitBase.h.

References clearHistoMap(), fillHistoMap(), download_sqlite_cfg::outputFile, readProbabilityDistributionsFromFile(), and writeHistoMap().

{}

Member Function Documentation

void MuScleFitBase::clearHistoMap ( )

protected

Clean the histograms map.

Definition at line 143 of file MuScleFitBase.cc.

References timingPdfMaker::histo, and mapHisto_.

Referenced by MuScleFit::endOfFastLoop(), and ~MuScleFitBase().

                                  {
  for (std::map<std::string, Histograms*>::const_iterator histo = mapHisto_.begin(); histo != mapHisto_.end();
       histo++) {
    delete (*histo).second;
  }
}

void MuScleFitBase::fillHistoMap ( TFile *  outputFile, unsigned int  iLoop  )

protected

Create the histograms map.

Definition at line 9 of file MuScleFitBase.cc.

References MuScleFitUtils::debugMassResol_, LogDebug, mapHisto_, B2GTnPMonitor_cfi::maxMass, MuonErrorMatrixAnalyzer_cfi::maxPt, B2GTnPMonitor_cfi::minMass, MuScleFitUtils::resfind, and theMuonType_.

Referenced by MuScleFit::startingNewLoop(), TestCorrection::TestCorrection(), and ~MuScleFitBase().

                                                                      {
  //Reconstructed muon kinematics
  //-----------------------------
  outputFile->cd();
  // If no Z is required, use a smaller mass range.
  double minMass = 0.;
  double maxMass = 200.;
  double maxPt = 100.;
  double yMaxEta = 4.;
  double yMaxPt = 2.;
  if (MuScleFitUtils::resfind[0] != 1) {
    maxMass = 20.;
    maxPt = 20.;
    yMaxEta = 0.2;
    yMaxPt = 0.2;
    // If running on standalone muons we need to expand the window range
    if (theMuonType_ == 2) {
      yMaxEta = 20.;
    }
  }

  LogDebug("MuScleFitBase") << "Creating new histograms" << std::endl;

  mapHisto_["hRecBestMu"] = new HParticle("hRecBestMu", minMass, maxMass, maxPt);
  mapHisto_["hRecBestMuVSEta"] = new HPartVSEta("hRecBestMuVSEta", minMass, maxMass, maxPt);
  //mapHisto_["hRecBestMuVSPhi"] = new HPartVSPhi ("hRecBestMuVSPhi");
  //mapHisto_["hRecBestMu_Acc"]  = new HParticle ("hRecBestMu_Acc", minMass, maxMass);
  mapHisto_["hDeltaRecBestMu"] = new HDelta("hDeltaRecBestMu");

  mapHisto_["hRecBestRes"] = new HParticle("hRecBestRes", minMass, maxMass, maxPt);
  mapHisto_["hRecBestResAllEvents"] = new HParticle("hRecBestResAllEvents", minMass, maxMass, maxPt);
  //mapHisto_["hRecBestRes_Acc"] = new HParticle   ("hRecBestRes_Acc", minMass, maxMass);
  // If not finding Z, use a smaller mass window
  mapHisto_["hRecBestResVSMu"] = new HMassVSPart("hRecBestResVSMu", minMass, maxMass, maxPt);
  mapHisto_["hRecBestResVSRes"] = new HMassVSPart("hRecBestResVSRes", minMass, maxMass, maxPt);
  //Generated Mass versus pt
  mapHisto_["hGenResVSMu"] = new HMassVSPart("hGenResVSMu", minMass, maxMass, maxPt);
  // Likelihood values VS muon variables
  // -------------------------------------
  mapHisto_["hLikeVSMu"] = new HLikelihoodVSPart("hLikeVSMu");
  mapHisto_["hLikeVSMuMinus"] = new HLikelihoodVSPart("hLikeVSMuMinus");
  mapHisto_["hLikeVSMuPlus"] = new HLikelihoodVSPart("hLikeVSMuPlus");

  //Resolution VS muon kinematic
  //----------------------------
  mapHisto_["hResolMassVSMu"] =
      new HResolutionVSPart(outputFile, "hResolMassVSMu", maxPt, 0., yMaxEta, 0., yMaxPt, true);
  mapHisto_["hFunctionResolMassVSMu"] =
      new HResolutionVSPart(outputFile, "hFunctionResolMassVSMu", maxPt, 0, 0.1, 0, 0.1, true);
  mapHisto_["hResolPtGenVSMu"] = new HResolutionVSPart(outputFile, "hResolPtGenVSMu", maxPt, -0.1, 0.1, -0.1, 0.1);
  mapHisto_["hResolPtSimVSMu"] = new HResolutionVSPart(outputFile, "hResolPtSimVSMu", maxPt, -0.1, 0.1, -0.1, 0.1);
  mapHisto_["hResolEtaGenVSMu"] =
      new HResolutionVSPart(outputFile, "hResolEtaGenVSMu", maxPt, -0.02, 0.02, -0.02, 0.02);
  mapHisto_["hResolEtaSimVSMu"] =
      new HResolutionVSPart(outputFile, "hResolEtaSimVSMu", maxPt, -0.02, 0.02, -0.02, 0.02);
  mapHisto_["hResolThetaGenVSMu"] =
      new HResolutionVSPart(outputFile, "hResolThetaGenVSMu", maxPt, -0.02, 0.02, -0.02, 0.02);
  mapHisto_["hResolThetaSimVSMu"] =
      new HResolutionVSPart(outputFile, "hResolThetaSimVSMu", maxPt, -0.02, 0.02, -0.02, 0.02);
  mapHisto_["hResolCotgThetaGenVSMu"] =
      new HResolutionVSPart(outputFile, "hResolCotgThetaGenVSMu", maxPt, -0.02, 0.02, -0.02, 0.02);
  mapHisto_["hResolCotgThetaSimVSMu"] =
      new HResolutionVSPart(outputFile, "hResolCotgThetaSimVSMu", maxPt, -0.02, 0.02, -0.02, 0.02);
  mapHisto_["hResolPhiGenVSMu"] =
      new HResolutionVSPart(outputFile, "hResolPhiGenVSMu", maxPt, -0.02, 0.02, -0.02, 0.02);
  mapHisto_["hResolPhiSimVSMu"] =
      new HResolutionVSPart(outputFile, "hResolPhiSimVSMu", maxPt, -0.02, 0.02, -0.02, 0.02);

  if (MuScleFitUtils::debugMassResol_) {
    mapHisto_["hdMdPt1"] = new HResolutionVSPart(outputFile, "hdMdPt1", maxPt, 0, 100, -3.2, 3.2, true);
    mapHisto_["hdMdPt2"] = new HResolutionVSPart(outputFile, "hdMdPt2", maxPt, 0, 100, -3.2, 3.2, true);
    mapHisto_["hdMdPhi1"] = new HResolutionVSPart(outputFile, "hdMdPhi1", maxPt, 0, 100, -3.2, 3.2, true);
    mapHisto_["hdMdPhi2"] = new HResolutionVSPart(outputFile, "hdMdPhi2", maxPt, 0, 100, -3.2, 3.2, true);
    mapHisto_["hdMdCotgTh1"] = new HResolutionVSPart(outputFile, "hdMdCotgTh1", maxPt, 0, 100, -3.2, 3.2, true);
    mapHisto_["hdMdCotgTh2"] = new HResolutionVSPart(outputFile, "hdMdCotgTh2", maxPt, 0, 100, -3.2, 3.2, true);
  }

  HTH2D* recoGenHisto =
      new HTH2D(outputFile, "hPtRecoVsPtGen", "Pt reco vs Pt gen", "hPtRecoVsPtGen", 120, 0., 120., 120, 0, 120.);
  (*recoGenHisto)->SetXTitle("Pt gen (GeV)");
  (*recoGenHisto)->SetYTitle("Pt reco (GeV)");
  mapHisto_["hPtRecoVsPtGen"] = recoGenHisto;
  HTH2D* recoSimHisto =
      new HTH2D(outputFile, "hPtRecoVsPtSim", "Pt reco vs Pt sim", "hPtRecoVsPtSim", 120, 0., 120., 120, 0, 120.);
  (*recoSimHisto)->SetXTitle("Pt sim (GeV)");
  (*recoSimHisto)->SetYTitle("Pt reco (GeV)");
  mapHisto_["hPtRecoVsPtSim"] = recoSimHisto;
  // Resolutions from resolution functions
  // -------------------------------------
  mapHisto_["hFunctionResolPt"] = new HFunctionResolution(outputFile, "hFunctionResolPt", maxPt);
  mapHisto_["hFunctionResolCotgTheta"] = new HFunctionResolution(outputFile, "hFunctionResolCotgTheta", maxPt);
  mapHisto_["hFunctionResolPhi"] = new HFunctionResolution(outputFile, "hFunctionResolPhi", maxPt);

  // Mass probability histograms
  // ---------------------------
  // The word "profile" is added to the title automatically
  mapHisto_["hMass_P"] = new HTProfile(outputFile, "Mass_P", "Mass probability", 4000, 0., 200., 0., 50.);
  mapHisto_["hMass_fine_P"] = new HTProfile(outputFile, "Mass_fine_P", "Mass probability", 4000, 0., 20., 0., 50.);
  mapHisto_["hMass_Probability"] = new HTH1D(outputFile, "Mass_Probability", "Mass probability", 4000, 0., 200.);
  mapHisto_["hMass_fine_Probability"] =
      new HTH1D(outputFile, "Mass_fine_Probability", "Mass probability", 4000, 0., 20.);
  mapHisto_["hMassProbVsMu"] = new HMassVSPartProfile("hMassProbVsMu", minMass, maxMass, maxPt);
  mapHisto_["hMassProbVsRes"] = new HMassVSPartProfile("hMassProbVsRes", minMass, maxMass, maxPt);
  mapHisto_["hMassProbVsMu_fine"] = new HMassVSPartProfile("hMassProbVsMu_fine", minMass, maxMass, maxPt);
  mapHisto_["hMassProbVsRes_fine"] = new HMassVSPartProfile("hMassProbVsRes_fine", minMass, maxMass, maxPt);

  // (M_reco-M_gen)/M_gen vs (pt, eta) of the muons from MC
  mapHisto_["hDeltaMassOverGenMassVsPt"] = new HTH2D(outputFile,
                                                     "DeltaMassOverGenMassVsPt",
                                                     "DeltaMassOverGenMassVsPt",
                                                     "DeltaMassOverGenMass",
                                                     200,
                                                     0,
                                                     maxPt,
                                                     200,
                                                     -0.2,
                                                     0.2);
  mapHisto_["hDeltaMassOverGenMassVsEta"] = new HTH2D(outputFile,
                                                      "DeltaMassOverGenMassVsEta",
                                                      "DeltaMassOverGenMassVsEta",
                                                      "DeltaMassOverGenMass",
                                                      200,
                                                      -3.,
                                                      3.,
                                                      200,
                                                      -0.2,
                                                      0.2);

  // Square of mass resolution vs (pt, eta) of the muons from MC
  // EM 2012.12.19  mapHisto_["hMassResolutionVsPtEta"] = new HCovarianceVSxy( "Mass", "Mass", 100, 0., maxPt, 60, -3, 3, outputFile->mkdir("MassCovariance") );
  // Mass resolution vs (pt, eta) from resolution function
  mapHisto_["hFunctionResolMass"] = new HFunctionResolution(outputFile, "hFunctionResolMass", maxPt);
}

void MuScleFitBase::readProbabilityDistributionsFromFile ( )

protected

Read probability distributions from a local root file.

Definition at line 159 of file MuScleFitBase.cc.

References gather_cfg::cout, beamvalidation::exit(), FrontierConditions_GlobalTag_cff::file, GL, MuScleFitUtils::GLNorm, MuScleFitUtils::GLValue, MuScleFitUtils::GLZNorm, MuScleFitUtils::GLZValue, mps_fire::i, ires, MuScleFitUtils::nbins, probabilitiesFile_, probabilitiesFileInPath_, MuScleFitUtils::rapidityBinsForZ_, MuScleFitUtils::resfind, MuScleFitUtils::ResHalfWidth, MuScleFitUtils::ResMaxSigma, MuScleFitUtils::ResMinMass, and theMuonType_.

Referenced by MuScleFit::beginOfJobInConstructor(), TestCorrection::initialize(), and ~MuScleFitBase().

                                                         {
  std::array<TH2D*, 6> GLZ = {{nullptr}};
  std::array<TH2D*, 6> GL = {{nullptr}};
  std::unique_ptr<TFile> ProbsFile;
  if (!probabilitiesFile_.empty()) {
    ProbsFile = std::make_unique<TFile>(probabilitiesFile_.c_str());
    std::cout << "[MuScleFit-Constructor]: Reading TH2D probabilities from " << probabilitiesFile_ << std::endl;
  } else {
    // edm::FileInPath file("MuonAnalysis/MomentumScaleCalibration/test/Probs_new_1000_CTEQ.root");
    // edm::FileInPath file("MuonAnalysis/MomentumScaleCalibration/test/Probs_new_Horace_CTEQ_1000.root");
    // edm::FileInPath file("MuonAnalysis/MomentumScaleCalibration/test/Probs_merge.root");
    edm::FileInPath file(probabilitiesFileInPath_.c_str());
    ProbsFile = std::make_unique<TFile>(file.fullPath().c_str());
    std::cout << "[MuScleFit-Constructor]: Reading TH2D probabilities from " << probabilitiesFileInPath_ << std::endl;
  }

  ProbsFile->cd();
  bool resfindEmpty = true;
  if (MuScleFitUtils::rapidityBinsForZ_ && MuScleFitUtils::resfind[0] && theMuonType_ != 2) {
    resfindEmpty = false;
    for (unsigned char i = 0; i < GLZ.size(); i++) {
      char nameh[6];
      snprintf(nameh, 6, "GLZ%hhu", i);
      GLZ[i] = dynamic_cast<TH2D*>(ProbsFile->Get(nameh));
    }
  } else if (MuScleFitUtils::resfind[0] && theMuonType_ == 2) {
    GL[0] = dynamic_cast<TH2D*>(ProbsFile->Get("GL0"));
    resfindEmpty = false;
  }
  for (unsigned char i = 1; i < GL.size(); ++i) {
    if (MuScleFitUtils::resfind[i]) {
      char nameh[6];
      snprintf(nameh, 6, "GL%hhu", i);
      GL[i] = dynamic_cast<TH2D*>(ProbsFile->Get(nameh));
      resfindEmpty = false;
    }
  }
  if (resfindEmpty) {
    std::cout << "[MuScleFit-Constructor]: No resonance selected, please fill the resfind array" << std::endl;
    exit(1);
  }

  // Read the limits for M and Sigma axis for each pdf
  // Note: we assume all the Z histograms to have the same limits
  // x is mass, y is sigma
  if (MuScleFitUtils::rapidityBinsForZ_ && MuScleFitUtils::resfind[0] && theMuonType_ != 2) {
    MuScleFitUtils::ResHalfWidth[0] = (GLZ[0]->GetXaxis()->GetXmax() - GLZ[0]->GetXaxis()->GetXmin()) / 2.;
    MuScleFitUtils::ResMaxSigma[0] = (GLZ[0]->GetYaxis()->GetXmax() - GLZ[0]->GetYaxis()->GetXmin());
    MuScleFitUtils::ResMinMass[0] = (GLZ[0]->GetXaxis()->GetXmin());
  }
  if (MuScleFitUtils::resfind[0] && theMuonType_ == 2) {
    MuScleFitUtils::ResHalfWidth[0] = (GL[0]->GetXaxis()->GetXmax() - GL[0]->GetXaxis()->GetXmin()) / 2.;
    MuScleFitUtils::ResMaxSigma[0] = (GL[0]->GetYaxis()->GetXmax() - GL[0]->GetYaxis()->GetXmin());
    MuScleFitUtils::ResMinMass[0] = (GL[0]->GetXaxis()->GetXmin());
  }
  for (unsigned int i = 1; i < GL.size(); ++i) {
    if (MuScleFitUtils::resfind[i]) {
      MuScleFitUtils::ResHalfWidth[i] = (GL[i]->GetXaxis()->GetXmax() - GL[i]->GetXaxis()->GetXmin()) / 2.;
      MuScleFitUtils::ResMaxSigma[i] = (GL[i]->GetYaxis()->GetXmax() - GL[i]->GetYaxis()->GetXmin());
      MuScleFitUtils::ResMinMass[i] = (GL[i]->GetXaxis()->GetXmin());
      // if( debug_>2 ) {
      std::cout << "MuScleFitUtils::ResHalfWidth[" << i << "] = " << MuScleFitUtils::ResHalfWidth[i] << std::endl;
      std::cout << "MuScleFitUtils::ResMaxSigma[" << i << "] = " << MuScleFitUtils::ResMaxSigma[i] << std::endl;
      // }
    }
  }

  // Extract normalization for mass slice in Y bins of Z
  // ---------------------------------------------------
  if (MuScleFitUtils::rapidityBinsForZ_ && MuScleFitUtils::resfind[0] && theMuonType_ != 2) {
    for (unsigned int iY = 0; iY < GLZ.size(); iY++) {
      int nBinsX = GLZ[iY]->GetNbinsX();
      int nBinsY = GLZ[iY]->GetNbinsY();
      if (nBinsX != MuScleFitUtils::nbins + 1 || nBinsY != MuScleFitUtils::nbins + 1) {
        std::cout << "Error: for histogram \"" << GLZ[iY]->GetName() << "\" bins are not " << MuScleFitUtils::nbins
                  << std::endl;
        std::cout << "nBinsX = " << nBinsX << ", nBinsY = " << nBinsY << std::endl;
        exit(1);
      }
      for (int iy = 0; iy <= MuScleFitUtils::nbins; iy++) {
        MuScleFitUtils::GLZNorm[iY][iy] = 0.;
        for (int ix = 0; ix <= MuScleFitUtils::nbins; ix++) {
          MuScleFitUtils::GLZValue[iY][ix][iy] = GLZ[iY]->GetBinContent(ix + 1, iy + 1);
          MuScleFitUtils::GLZNorm[iY][iy] +=
              MuScleFitUtils::GLZValue[iY][ix][iy] * (2 * MuScleFitUtils::ResHalfWidth[0]) / MuScleFitUtils::nbins;
        }
      }
    }
  }

  if (MuScleFitUtils::resfind[0] && theMuonType_ == 2) {
    int nBinsX = GL[0]->GetNbinsX();
    int nBinsY = GL[0]->GetNbinsY();
    if (nBinsX != MuScleFitUtils::nbins + 1 || nBinsY != MuScleFitUtils::nbins + 1) {
      std::cout << "Error: for histogram \"" << GL[0]->GetName() << "\" bins are not " << MuScleFitUtils::nbins
                << std::endl;
      std::cout << "nBinsX = " << nBinsX << ", nBinsY = " << nBinsY << std::endl;
      exit(1);
    }

    for (int iy = 0; iy <= MuScleFitUtils::nbins; iy++) {
      MuScleFitUtils::GLNorm[0][iy] = 0.;
      for (int ix = 0; ix <= MuScleFitUtils::nbins; ix++) {
        MuScleFitUtils::GLValue[0][ix][iy] = GL[0]->GetBinContent(ix + 1, iy + 1);
        // N.B. approximation: we should compute the integral of the function used to compute the probability (linear
        // interpolation of the mass points). This computation could be troublesome because the points have a steep
        // variation near the mass peak and the normal integral is not precise in these conditions.
        // Furthermore it is slow.
        MuScleFitUtils::GLNorm[0][iy] +=
            MuScleFitUtils::GLValue[0][ix][iy] * (2 * MuScleFitUtils::ResHalfWidth[0]) / MuScleFitUtils::nbins;
      }
    }
  }
  // Extract normalization for each mass slice
  // -----------------------------------------
  for (unsigned int ires = 1; ires < GL.size(); ires++) {
    if (MuScleFitUtils::resfind[ires]) {
      int nBinsX = GL[ires]->GetNbinsX();
      int nBinsY = GL[ires]->GetNbinsY();
      if (nBinsX != MuScleFitUtils::nbins + 1 || nBinsY != MuScleFitUtils::nbins + 1) {
        std::cout << "Error: for histogram \"" << GL[ires]->GetName() << "\" bins are not " << MuScleFitUtils::nbins
                  << std::endl;
        std::cout << "nBinsX = " << nBinsX << ", nBinsY = " << nBinsY << std::endl;
        exit(1);
      }

      for (int iy = 0; iy <= MuScleFitUtils::nbins; iy++) {
        MuScleFitUtils::GLNorm[ires][iy] = 0.;
        for (int ix = 0; ix <= MuScleFitUtils::nbins; ix++) {
          MuScleFitUtils::GLValue[ires][ix][iy] = GL[ires]->GetBinContent(ix + 1, iy + 1);
          // N.B. approximation: we should compute the integral of the function used to compute the probability (linear
          // interpolation of the mass points). This computation could be troublesome because the points have a steep
          // variation near the mass peak and the normal integral is not precise in these conditions.
          // Furthermore it is slow.
          MuScleFitUtils::GLNorm[ires][iy] +=
              MuScleFitUtils::GLValue[ires][ix][iy] * (2 * MuScleFitUtils::ResHalfWidth[ires]) / MuScleFitUtils::nbins;
        }
      }
    }
  }
  // Free all the memory for the probability histograms.
  if (MuScleFitUtils::rapidityBinsForZ_ && MuScleFitUtils::resfind[0] && theMuonType_ != 2) {
    for (unsigned int i = 0; i < GLZ.size(); i++) {
      delete GLZ[i];
    }
  }
  if (MuScleFitUtils::resfind[0] && theMuonType_ == 2)
    delete GL[0];
  for (unsigned int ires = 1; ires < GL.size(); ires++) {
    if (MuScleFitUtils::resfind[ires])
      delete GL[ires];
  }
}

void MuScleFitBase::writeHistoMap ( const unsigned int  iLoop )

protected

Save the histograms map to file.

Definition at line 150 of file MuScleFitBase.cc.

References timingPdfMaker::histo, mapHisto_, and theFiles_.

Referenced by MuScleFit::endOfFastLoop(), ~MuScleFitBase(), and TestCorrection::~TestCorrection().

                                                          {
  for (std::map<std::string, Histograms*>::const_iterator histo = mapHisto_.begin(); histo != mapHisto_.end();
       histo++) {
    // This is to avoid writing into subdirs. Need a workaround.
    theFiles_[iLoop]->cd();
    (*histo).second->Write();
  }
}

Member Data Documentation

int MuScleFitBase::debug_

protected

Definition at line 51 of file MuScleFitBase.h.

Referenced by MuScleFit::applyBias(), MuScleFit::applySmearing(), MuScleFit::beginOfJobInConstructor(), MuScleFit::checkDeltaR(), MuScleFit::duringFastLoop(), MuScleFit::duringLoop(), MuScleFit::endOfJob(), TestCorrection::fillMuonCollection(), MuScleFit::MuScleFit(), MuScleFit::selectMuons(), MuScleFit::startingNewLoop(), and MuScleFit::~MuScleFit().

std::vector<GenMuonPair> MuScleFitBase::genMuonPairs_

protected

Stores the genMuon pairs and the motherId prior to the creation of the internal tree.

Definition at line 82 of file MuScleFitBase.h.

Referenced by MuScleFit::selectMuons(), and MuScleFit::~MuScleFit().

std::map<std::string, Histograms*> MuScleFitBase::mapHisto_

protected

The map of histograms.

Definition at line 77 of file MuScleFitBase.h.

Referenced by TestCorrection::analyze(), clearHistoMap(), MuScleFit::duringFastLoop(), MuScleFit::fillComparisonHistograms(), fillHistoMap(), and writeHistoMap().

std::vector<MuonPair> MuScleFitBase::muonPairs_

protected

Used to store the muon pairs plus run and event number prior to the creation of the internal tree.

Definition at line 80 of file MuScleFitBase.h.

Referenced by MuScleFit::selectMuons(), and MuScleFit::~MuScleFit().

std::string MuScleFitBase::probabilitiesFile_

protected

Definition at line 44 of file MuScleFitBase.h.

Referenced by readProbabilityDistributionsFromFile().

std::string MuScleFitBase::probabilitiesFileInPath_

protected

Definition at line 43 of file MuScleFitBase.h.

Referenced by readProbabilityDistributionsFromFile().

std::vector<TFile*> MuScleFitBase::theFiles_

protected

The files were the histograms are saved.

Definition at line 74 of file MuScleFitBase.h.

Referenced by MuScleFit::beginOfJobInConstructor(), MuScleFit::endOfFastLoop(), MuScleFit::startingNewLoop(), TestCorrection::TestCorrection(), writeHistoMap(), and TestCorrection::~TestCorrection().

std::string MuScleFitBase::theGenInfoRootFileName_

protected

Definition at line 49 of file MuScleFitBase.h.

Referenced by MuScleFit::beginOfJobInConstructor().

edm::InputTag MuScleFitBase::theMuonLabel_

protected

Definition at line 47 of file MuScleFitBase.h.

Referenced by MuScleFit::MuScleFit().

int MuScleFitBase::theMuonType_

protected

Definition at line 46 of file MuScleFitBase.h.

Referenced by TestCorrection::analyze(), fillHistoMap(), MuScleFit::MuScleFit(), readProbabilityDistributionsFromFile(), MuScleFit::selectMuons(), and MuScleFit::~MuScleFit().

std::string MuScleFitBase::theRootFileName_

protected

Definition at line 48 of file MuScleFitBase.h.

Referenced by MuScleFit::beginOfJobInConstructor(), and TestCorrection::TestCorrection().