HMassResolutionVSPart Class Reference

#include <Histograms.h>

Inherits Histograms.

Public Member Functions
void	Clear () override
void	Fill (const reco::Particle::LorentzVector &recoP1, const int charge1, const reco::Particle::LorentzVector &genP1, const reco::Particle::LorentzVector &recoP2, const int charge2, const reco::Particle::LorentzVector &genP2, const double &recoMass, const double &genMass) override
void	Fill (const reco::Particle::LorentzVector &recoP1, const int charge1, const reco::Particle::LorentzVector &recoP2, const int charge2, const double &recoMass, const double &genMass) override
	HMassResolutionVSPart (TFile *outputFile, const TString &name)
void	Write () override
	~HMassResolutionVSPart () override

Protected Attributes
std::map< TString, TH1 * >	mapHisto_
std::unique_ptr< HDelta >	muMinus
std::unique_ptr< HDelta >	muPlus
TString	nameSuffix_ [2]

Detailed Description

A set of histograms for resolution. The fill method requires the two selected muons, their charges and the reconstructed and generated masses. It evaluates the resolution on the mass vs: Pt of the pair DeltaEta of the pair DeltaPhi of the pair pt, eta and phi of the plus and minus muon separately

Definition at line 2265 of file Histograms.h.

Constructor & Destructor Documentation

HMassResolutionVSPart::HMassResolutionVSPart ( TFile *  outputFile, const TString &  name  )

inline

Definition at line 2267 of file Histograms.h.

References mps_fire::i, and Skims_PA_cff::name.

                                                                : Histograms(outputFile, name) {
    // Kinematical variables
    nameSuffix_[0] = "Plus";
    nameSuffix_[1] = "Minus";
    TString titleSuffix[] = {" for mu+", " for mu-"};

    mapHisto_[name] = new TH1F(name, "#Delta M/M", 4000, -1, 1);
    mapHisto_[name + "VSPairPt"] =
        new HResolution(name + "VSPairPt", "resolution VS pt of the pair", 100, 0, 200, -1, 1, histoDir_);
    mapHisto_[name + "VSPairDeltaEta"] = new HResolution(
        name + "VSPairDeltaEta", "resolution VS #Delta#eta of the pair", 100, -0.1, 6.2, -1, 1, histoDir_);
    mapHisto_[name + "VSPairDeltaPhi"] = new HResolution(
        name + "VSPairDeltaPhi", "resolution VS #Delta#phi of the pair", 100, -0.1, 3.2, -1, 1, histoDir_);

    for (int i = 0; i < 2; ++i) {
      mapHisto_[name + "VSPt" + nameSuffix_[i]] = new HResolution(
          name + "VSPt" + nameSuffix_[i], "resolution VS pt" + titleSuffix[i], 100, 0, 200, -1, 1, histoDir_);
      mapHisto_[name + "VSEta" + nameSuffix_[i]] = new HResolution(
          name + "VSEta" + nameSuffix_[i], "resolution VS #eta" + titleSuffix[i], 100, -3, 3, -1, 1, histoDir_);
      mapHisto_[name + "VSPhi" + nameSuffix_[i]] = new HResolution(
          name + "VSPhi" + nameSuffix_[i], "resolution VS #phi" + titleSuffix[i], 100, -3.2, 3.2, -1, 1, histoDir_);
    }

    // single particles histograms
    muMinus.reset(new HDelta("muMinus"));
    muPlus.reset(new HDelta("muPlus"));
  }

HMassResolutionVSPart::~HMassResolutionVSPart ( )

inlineoverride

Definition at line 2295 of file Histograms.h.

References timingPdfMaker::histo.

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

Member Function Documentation

void HMassResolutionVSPart::Clear ( )

inlineoverride

Definition at line 2377 of file Histograms.h.

References timingPdfMaker::histo.

                        {
    for (std::map<TString, TH1*>::const_iterator histo = mapHisto_.begin(); histo != mapHisto_.end(); histo++) {
      (*histo).second->Clear();
    }
    muMinus->Clear();
    muPlus->Clear();
  }

void HMassResolutionVSPart::Fill ( const reco::Particle::LorentzVector &  recoP1, const int  charge1, const reco::Particle::LorentzVector &  genP1, const reco::Particle::LorentzVector &  recoP2, const int  charge2, const reco::Particle::LorentzVector &  genP2, const double &  recoMass, const double &  genMass  )

inlineoverride

Definition at line 2301 of file Histograms.h.

References HcalObjRepresent::Fill().

                                            {
    muMinus->Fill(recoP1, genP1);
    muPlus->Fill(recoP2, genP2);

    Fill(recoP1, charge1, recoP2, charge2, recoMass, genMass);
  }

void HMassResolutionVSPart::Fill ( const reco::Particle::LorentzVector &  recoP1, const int  charge1, const reco::Particle::LorentzVector &  recoP2, const int  charge2, const double &  recoMass, const double &  genMass  )

inlineoverride

Definition at line 2315 of file Histograms.h.

References gather_cfg::cout, MuScleFitUtils::deltaPhi(), and mps_fire::i.

                                            {
    if (charge1 == charge2)
      std::cout << "Error: must get two opposite charge particles" << std::endl;

    double massRes = (recoMass - genMass) / genMass;

    reco::Particle::LorentzVector recoPair(recoP1 + recoP2);
    double pairPt = recoPair.Pt();

    double recoPt[2] = {recoP1.Pt(), recoP2.Pt()};
    double recoEta[2] = {recoP1.Eta(), recoP2.Eta()};
    double recoPhi[2] = {recoP1.Phi(), recoP2.Phi()};

    // std::cout << "pairPt = " << pairPt << ", massRes = ("<<recoMass<<" - "<<genMass<<")/"<<genMass<<" = " << massRes
    //      << ", recoPt[0] = " << recoPt[0] << ", recoPt[1] = " << recoPt[1] << std::endl;

    // Index of the histogram. If the muons have charge1 = -1 and charge2 = 1, they already have the
    // correct histogram indeces. Otherwise, swap the indeces.
    // In any case the negative muon has index = 0  and the positive muon has index = 1.
    int id[2] = {0, 1};
    if (charge1 > 0) {
      id[0] = 1;
      id[1] = 0;
    }

    double pairDeltaEta = fabs(recoEta[0] - recoEta[1]);
    double pairDeltaPhi = MuScleFitUtils::deltaPhi(recoPhi[0], recoPhi[1]);

    mapHisto_[name_]->Fill(massRes);
    mapHisto_[name_ + "VSPairPt"]->Fill(pairPt, massRes);
    mapHisto_[name_ + "VSPairDeltaEta"]->Fill(pairDeltaEta, massRes);
    mapHisto_[name_ + "VSPairDeltaPhi"]->Fill(pairDeltaPhi, massRes);

    // Resolution vs single muon quantities
    // ------------------------------------
    int index = 0;
    for (int i = 0; i < 2; ++i) {
      index = id[i];
      mapHisto_[name_ + "VSPt" + nameSuffix_[i]]->Fill(recoPt[index], massRes);  // EM [index] or [i]???
      mapHisto_[name_ + "VSEta" + nameSuffix_[i]]->Fill(recoEta[index], massRes);
      mapHisto_[name_ + "VSPhi" + nameSuffix_[i]]->Fill(recoPhi[index], massRes);
    }
  }

void HMassResolutionVSPart::Write ( )

inlineoverride

Definition at line 2366 of file Histograms.h.

References hippyaddtobaddatafiles::cd(), and timingPdfMaker::histo.

                        {
    histoDir_->cd();
    for (std::map<TString, TH1*>::const_iterator histo = mapHisto_.begin(); histo != mapHisto_.end(); histo++) {
      (*histo).second->Write();
    }
    // Create the new dir and cd into it
    (histoDir_->mkdir("singleMuonsVSgen"))->cd();
    muMinus->Write();
    muPlus->Write();
  }

Member Data Documentation

std::map<TString, TH1*> HMassResolutionVSPart::mapHisto_

protected

Definition at line 2400 of file Histograms.h.

std::unique_ptr<HDelta> HMassResolutionVSPart::muMinus

protected

Definition at line 2402 of file Histograms.h.

std::unique_ptr<HDelta> HMassResolutionVSPart::muPlus

protected

Definition at line 2403 of file Histograms.h.

TString HMassResolutionVSPart::nameSuffix_[2]

protected

Definition at line 2401 of file Histograms.h.