MuScleFitPlotter.h

Go to the documentation of this file.

#ifndef MuScleFitPlotter_H
#define MuScleFitPlotter_H

// #include "SimDataFormats/HepMCProduct/interface/HepMCProduct.h"
#include "SimDataFormats/GeneratorProducts/interface/HepMCProduct.h"
#include "DataFormats/HepMCCandidate/interface/GenParticleFwd.h"
#include "SimDataFormats/Track/interface/SimTrackContainer.h"
#include "Histograms.h"
#include "MuonAnalysis/MomentumScaleCalibration/interface/Muon.h"

namespace edm {
  class ParameterSet;
  class Event;
  class EventSetup;
}  // namespace edm

class TFile;
class LeafCandidate;

class MuScleFitPlotter {
public:
  // Constructor
  // -----------
  MuScleFitPlotter(std::string);

  // Destructor
  // ----------
  virtual ~MuScleFitPlotter();

  // Operations
  // ----------
  void fillGen(const reco::GenParticleCollection& genParticles, bool = false);
  void fillGen(const edm::HepMCProduct& evtMC, bool sherpaFlag_);
  void fillSim(edm::Handle<edm::SimTrackContainer> simTracks);
  void fillGenSim(edm::Handle<edm::HepMCProduct> evtMC, edm::Handle<edm::SimTrackContainer> simTracks);
  //void fillRec(std::vector<reco::LeafCandidate>& muons);
  void fillRec(std::vector<MuScleFitMuon>& muons);

  // Root tree specific
  void fillTreeRec(
      const std::vector<std::pair<reco::Particle::LorentzVector, reco::Particle::LorentzVector> >& savedPairs);
  void fillTreeGen(
      const std::vector<std::pair<reco::Particle::LorentzVector, reco::Particle::LorentzVector> >& genPairs);

  void fillHistoMap();
  void writeHistoMap();

  bool debug;

protected:
private:
  // The map of histograms
  // ---------------------
  std::map<std::string, Histograms*> mapHisto;
  TFile* outputFile;
};
#endif