HMassVSPart Class Reference

#include <Histograms.h>

Inherits Histograms.

Public Member Functions
void	Clear () override
void	Fill (const reco::Particle::LorentzVector &p41, const reco::Particle::LorentzVector &p42, const int charge, const double &weight=1.) override
void	Fill (const reco::Particle::LorentzVector &p41, const reco::Particle::LorentzVector &p42, const reco::Particle::LorentzVector &p4Res, const double &weight=1.) override
void	Fill (const CLHEP::HepLorentzVector &momentum1, const CLHEP::HepLorentzVector &momentum2, const CLHEP::HepLorentzVector &momentumRes, const double &weight=1.) override
	Used to fill 2D histograms for comparison of opposite charge muons quantities. More...
void	Fill (const CLHEP::HepLorentzVector &momentum1, const CLHEP::HepLorentzVector &momentum2, const int charge, const double &weight=1.) override
	HMassVSPart (const TString &name, const double &minMass=0., const double &maxMass=150., const double maxPt=100.)
	HMassVSPart (const TString &name, TFile *file)
void	Write () override
	~HMassVSPart () override

Protected Attributes
TH2F *	hMassVSCosThetaCS_
TProfile *	hMassVSCosThetaCS_prof
TH2F *	hMassVSEta_
TH2F *	hMassVSEtaMinus_
TH3F *	hMassVSEtaPhiMinus_
TH3F *	hMassVSEtaPhiPlus_
TH2F *	hMassVSEtaPlus_
TH3F *	hMassVSEtaPlusEtaMinus_
TH2F *	hMassVSEtaPlusMinusDiff_
TH2F *	hMassVSPhiCS_
TH2F *	hMassVSPhiMinus_
TH2F *	hMassVSPhiPlus_
TH2F *	hMassVSPhiPlusMinusDiff_
TH3F *	hMassVSPhiPlusPhiMinus_
TH2F *	hMassVSPt_

Detailed Description

Definition at line 776 of file Histograms.h.

Constructor & Destructor Documentation

HMassVSPart::HMassVSPart ( const TString &  name, const double &  minMass = 0., const double &  maxMass = 150., const double  maxPt = 100.  )

inline

Definition at line 778 of file Histograms.h.

References B2GTnPMonitor_cfi::maxMass, MuonErrorMatrixAnalyzer_cfi::maxPt, B2GTnPMonitor_cfi::minMass, and Skims_PA_cff::name.

                                                                                                                        {
    name_ = name;

    // Kinematical variables
    // ---------------------
    hMassVSPt_ = new TH2F(
        name + "_MassVSPt", "resonance mass vs muon transverse momentum", 200, 0., maxPt, 6000, minMass, maxMass);

    hMassVSEta_ =
        new TH2F(name + "_MassVSEta", "resonance mass vs muon pseudorapidity", 64, -6.4, 6.4, 6000, minMass, maxMass);
    hMassVSEtaPlus_ = new TH2F(
        name + "_MassVSEtaPlus", "resonance mass vs muon+ pseudorapidity", 64, -6.4, 6.4, 6000, minMass, maxMass);
    hMassVSEtaMinus_ = new TH2F(
        name + "_MassVSEtaMinus", "resonance mass vs muon- pseudorapidity", 64, -6.4, 6.4, 6000, minMass, maxMass);

    hMassVSPhiPlus_ =
        new TH2F(name + "_MassVSPhiPlus", "resonance mass vs muon+ phi angle", 64, -3.2, 3.2, 6000, minMass, maxMass);
    hMassVSPhiMinus_ =
        new TH2F(name + "_MassVSPhiMinus", "resonance mass vs muon- phi angle", 64, -3.2, 3.2, 6000, minMass, maxMass);

    // J/Psi mass -----
    //     hMassVSEtaPhiPlus_      = new TH3F( name+"_MassVSEtaPhiPlus", "resonance mass vs muon+ phi/pseudorapidity",6, -3.2, 3.2, 20, -2.5, 2.5, 6000, minMass, maxMass );
    //     hMassVSEtaPhiMinus_      = new TH3F( name+"_MassVSEtaPhiMinus", "resonance mass vs muon- phi/pseudorapidity", 6, -3.2, 3.2, 20, -2.5, 2.5, 6000, minMass, maxMass );

    //Z mass -----------
    hMassVSEtaPhiPlus_ = new TH3F(name + "_MassVSEtaPhiPlus",
                                  "resonance mass vs muon+ phi/pseudorapidity",
                                  16,
                                  -3.2,
                                  3.2,
                                  20,
                                  -2.4,
                                  2.4,
                                  300,
                                  minMass,
                                  maxMass);
    hMassVSEtaPhiMinus_ = new TH3F(name + "_MassVSEtaPhiMinus",
                                   "resonance mass vs muon- phi/pseudorapidity",
                                   16,
                                   -3.2,
                                   3.2,
                                   20,
                                   -2.4,
                                   2.4,
                                   300,
                                   minMass,
                                   maxMass);

    hMassVSCosThetaCS_ = new TH2F(
        name + "_MassVSCosThetaCS", "resonance mass vs cos(theta) (CS frame)", 40, -1., 1., 6000, minMass, maxMass);
    hMassVSPhiCS_ =
        new TH2F(name + "_MassVSPhiCS", "resonance mass vs phi (CS frame)", 64, -3.2, 3.2, 6000, minMass, maxMass);

    hMassVSPhiPlusPhiMinus_ = new TH3F(name + "_MassVSPhiPlusPhiMinus",
                                       "resonance mass vs muon+ phi/muon- phi",
                                       16,
                                       -3.2,
                                       3.2,
                                       16,
                                       -3.2,
                                       3.2,
                                       6000,
                                       minMass,
                                       maxMass);
    hMassVSEtaPlusEtaMinus_ = new TH3F(name + "_MassVSEtaPlusEtaMinus",
                                       "resonance mass vs muon+ eta/muon- eta",
                                       16,
                                       -3.2,
                                       3.2,
                                       16,
                                       -3.2,
                                       3.2,
                                       6000,
                                       minMass,
                                       maxMass);

    hMassVSPhiPlusMinusDiff_ = new TH2F(name + "_MassVSPhiPlusMinusDiff",
                                        "resonance mass vs delta phi between mu+/mu-",
                                        64,
                                        -6.4,
                                        6.4,
                                        6000,
                                        minMass,
                                        maxMass);
    hMassVSEtaPlusMinusDiff_ = new TH2F(name + "_MassVSEtaPlusMinusDiff",
                                        "resonance mass vs delta pseudorapidity between mu+/mu-",
                                        32,
                                        -4.4,
                                        4.4,
                                        6000,
                                        minMass,
                                        maxMass);
    hMassVSCosThetaCS_prof =
        new TProfile(name + "_MassVScosTheta_prof", "resonance mass vs cosTheta", 40, -1., 1., 85., 95.);

    //hMassVSPt_prof       = new TProfile (name+"_MassVSPt_prof", "resonance mass vs muon transverse momentum", 100, 0., 200., minMass, maxMass);
    //hMassVSEta_prof      = new TProfile (name+"_MassVSEta_prof", "resonance mass vs muon pseudorapidity", 30, -6., 6., minMass, maxMass);
    //hMassVSPhiPlus_prof  = new TProfile (name+"_MassVSPhiPlus_prof", "resonance mass vs muon+ phi angle", 32, -3.2, 3.2, minMass, maxMass);
    //hMassVSPhiMinus_prof = new TProfile (name+"_MassVSPhiMinus_prof", "resonance mass vs muon- phi angle", 32, -3.2, 3.2, minMass, maxMass);
  }

HMassVSPart::HMassVSPart ( const TString &  name, TFile *  file  )

inline

Definition at line 879 of file Histograms.h.

References Skims_PA_cff::name.

                                                {
    name_ = name;
    hMassVSPt_ = (TH2F*)file->Get(name + "_MassVSPt");
    hMassVSEta_ = (TH2F*)file->Get(name + "_MassVSEta");

    hMassVSEtaPhiPlus_ = (TH3F*)file->Get(name + "_MassVSEtaPlus");
    hMassVSEtaPhiMinus_ = (TH3F*)file->Get(name + "_MassVSEtaMinus");
    hMassVSEtaPlus_ = (TH2F*)file->Get(name + "_MassVSEtaPlus");
    hMassVSEtaMinus_ = (TH2F*)file->Get(name + "_MassVSEtaMinus");

    hMassVSPhiPlusMinusDiff_ = (TH2F*)file->Get(name + "_MassVSPhiPlusMinusDiff");
    hMassVSEtaPlusMinusDiff_ = (TH2F*)file->Get(name + "_MassVSEtaPlusMinusDiff");

    hMassVSPhiPlus_ = (TH2F*)file->Get(name + "_MassVSPhiPlus");
    hMassVSPhiMinus_ = (TH2F*)file->Get(name + "_MassVSPhiMinus");

    hMassVSCosThetaCS_prof = (TProfile*)file->Get(name + "_MassVScosTheta_prof");
    //hMassVSPt_prof       = (TProfile *) file->Get(name+"_MassVSPt_prof");
    //hMassVSEta_prof      = (TProfile *) file->Get(name+"_MassVSEta_prof");
    //hMassVSPhiPlus_prof  = (TProfile *) file->Get(name+"_MassVSPhiPlus_prof");
    //hMassVSPhiMinus_prof = (TProfile *) file->Get(name+"_MassVSPhiMinus_prof");
  }

HMassVSPart::~HMassVSPart ( )

inlineoverride

Definition at line 902 of file Histograms.h.

                          {
    delete hMassVSPt_;
    delete hMassVSEta_;
    delete hMassVSPhiPlus_;
    delete hMassVSPhiMinus_;
    delete hMassVSEtaPhiPlus_;
    delete hMassVSEtaPhiMinus_;
    delete hMassVSEtaPlus_;
    delete hMassVSEtaMinus_;
    delete hMassVSPhiPlusPhiMinus_;
    delete hMassVSEtaPlusEtaMinus_;
    delete hMassVSCosThetaCS_;
    delete hMassVSPhiCS_;
    delete hMassVSPhiPlusMinusDiff_;
    delete hMassVSEtaPlusMinusDiff_;
    delete hMassVSCosThetaCS_prof;
  }

Member Function Documentation

void HMassVSPart::Clear ( )

inlineoverride

Definition at line 1053 of file Histograms.h.

                        {
    hMassVSPt_->Clear();
    hMassVSEta_->Clear();
    hMassVSPhiPlus_->Clear();
    hMassVSPhiMinus_->Clear();

    hMassVSEtaPhiPlus_->Clear();
    hMassVSEtaPhiMinus_->Clear();
    hMassVSEtaPlus_->Clear();
    hMassVSEtaMinus_->Clear();

    hMassVSPhiPlusPhiMinus_->Clear();
    hMassVSEtaPlusEtaMinus_->Clear();
    hMassVSCosThetaCS_->Clear();
    hMassVSPhiCS_->Clear();
    hMassVSPhiPlusMinusDiff_->Clear();
    hMassVSEtaPlusMinusDiff_->Clear();
    hMassVSCosThetaCS_prof->Clear();

    //hMassVSPt_prof_->Clear();
    //hMassVSEta_prof_->Clear();
    //hMassVSPhiPlus_prof_->Clear();
    //hMassVSPhiMinus_prof_->Clear();
  }

void HMassVSPart::Fill ( const reco::Particle::LorentzVector &  p41, const reco::Particle::LorentzVector &  p42, const int  charge, const double &  weight = 1.  )

inlineoverride

Definition at line 920 of file Histograms.h.

References HcalObjRepresent::Fill().

                                                {
    Fill(CLHEP::HepLorentzVector(p41.x(), p41.y(), p41.z(), p41.t()),
         CLHEP::HepLorentzVector(p42.x(), p42.y(), p42.z(), p42.t()),
         charge,
         weight);
  }

void HMassVSPart::Fill ( const reco::Particle::LorentzVector &  p41, const reco::Particle::LorentzVector &  p42, const reco::Particle::LorentzVector &  p4Res, const double &  weight = 1.  )

inlineoverride

Definition at line 930 of file Histograms.h.

References HcalObjRepresent::Fill().

                                                {
    Fill(CLHEP::HepLorentzVector(p41.x(), p41.y(), p41.z(), p41.t()),
         CLHEP::HepLorentzVector(p42.x(), p42.y(), p42.z(), p42.t()),
         CLHEP::HepLorentzVector(p4Res.x(), p4Res.y(), p4Res.z(), p4Res.t()),
         weight);
  }

void HMassVSPart::Fill ( const CLHEP::HepLorentzVector &  momentum1, const CLHEP::HepLorentzVector &  momentum2, const CLHEP::HepLorentzVector &  momentumRes, const double &  weight = 1.  )

inlineoverride

Used to fill 2D histograms for comparison of opposite charge muons quantities.

Definition at line 941 of file Histograms.h.

References amptDefaultParameters_cff::mu, RPCpg::mubar, BPhysicsValidation_cfi::muminus, BPhysicsValidation_cfi::muplus, funct::pow(), class-composition::Q, dttmaxenums::R, mathSSE::sqrt(), and mps_merge::weight.

                                                {
    /************************************************************************
     *
     * Observable: cos(theta) = 2 Q^-1 (Q^2+Qt^2)^-(1/2) (mu^+ mubar^- - mu^- mubar^+)
     * (computed in Collins-Soper frame)
     *
     ************************************************************************/

    double costhetaCS, phiCS;

    const CLHEP::HepLorentzVector& mu = momentum1;
    const CLHEP::HepLorentzVector& mubar = momentum2;
    CLHEP::HepLorentzVector Q(mu + mubar);
    double muplus = 1.0 / sqrt(2.0) * (mu.e() + mu.z());
    double muminus = 1.0 / sqrt(2.0) * (mu.e() - mu.z());
    double mubarplus = 1.0 / sqrt(2.0) * (mubar.e() + mubar.z());
    double mubarminus = 1.0 / sqrt(2.0) * (mubar.e() - mubar.z());
    //double costheta = 2.0 / Q.Mag() / sqrt(pow(Q.Mag(), 2) + pow(Q.Pt(), 2)) * (muplus * mubarminus - muminus * mubarplus);
    costhetaCS = 2.0 / Q.mag() / sqrt(pow(Q.mag(), 2) + pow(Q.perp(), 2)) * (muplus * mubarminus - muminus * mubarplus);
    if (momentumRes.rapidity() < 0)
      costhetaCS = -costhetaCS;

    /************************************************************************
   *
   * 3) tanphi = (Q^2 + Qt^2)^1/2 / Q (Dt dot R unit) /(Dt dot Qt unit)
   *
   ************************************************************************/

    // unit vector on R direction
    CLHEP::HepLorentzVector Pbeam(0., 0., 3500., 3500.);
    CLHEP::Hep3Vector R = Pbeam.vect().cross(Q.vect());
    CLHEP::Hep3Vector Runit = R.unit();

    // unit vector on Qt
    CLHEP::Hep3Vector Qt = Q.vect();
    Qt.setZ(0);
    CLHEP::Hep3Vector Qtunit = Qt.unit();

    CLHEP::HepLorentzVector D(mu - mubar);
    CLHEP::Hep3Vector Dt = D.vect();
    Dt.setZ(0);
    double tanphi = sqrt(pow(Q.mag(), 2) + pow(Q.perp(), 2)) / Q.mag() * Dt.dot(Runit) / Dt.dot(Qtunit);
    if (momentumRes.rapidity() < 0)
      tanphi = -tanphi;
    phiCS = atan(tanphi);

    hMassVSPhiCS_->Fill(phiCS, momentumRes.m(), weight);
    hMassVSCosThetaCS_->Fill(costhetaCS, momentumRes.m(), weight);
    hMassVSCosThetaCS_prof->Fill(costhetaCS, momentumRes.m());
    /*************************************************************************
     *************************************************************************/

    hMassVSPhiPlusPhiMinus_->Fill(momentum1.phi(), momentum2.phi(), momentumRes.m(), weight);
    hMassVSEtaPlusEtaMinus_->Fill(momentum1.eta(), momentum2.eta(), momentumRes.m(), weight);

    hMassVSPhiPlusMinusDiff_->Fill((momentum1.phi() - momentum2.phi()), momentumRes.m(), weight);
    hMassVSEtaPlusMinusDiff_->Fill((momentum1.eta() - momentum2.eta()), momentumRes.m(), weight);
  }

void HMassVSPart::Fill ( const CLHEP::HepLorentzVector &  momentum1, const CLHEP::HepLorentzVector &  momentum2, const int  charge, const double &  weight = 1.  )

inlineoverride

Definition at line 1003 of file Histograms.h.

References LogDebug, and mps_merge::weight.

                                                {
    hMassVSPt_->Fill(momentum1.perp(), momentum2.m(), weight);
    //hMassVSPt_prof_->Fill(momentum1.perp(),momentum2.m());

    hMassVSEta_->Fill(momentum1.eta(), momentum2.m(), weight);
    //hMassVSEta_prof_->Fill(momentum1.eta(),momentum2.m());

    if (charge > 0) {
      hMassVSPhiPlus_->Fill(momentum1.phi(), momentum2.m(), weight);
      hMassVSEtaPlus_->Fill(momentum1.eta(), momentum2.m(), weight);
      hMassVSEtaPhiPlus_->Fill(momentum1.phi(), momentum1.eta(), momentum2.m(), weight);
    } else if (charge < 0) {
      hMassVSPhiMinus_->Fill(momentum1.phi(), momentum2.m(), weight);
      hMassVSEtaMinus_->Fill(momentum1.eta(), momentum2.m(), weight);
      hMassVSEtaPhiMinus_->Fill(momentum1.phi(), momentum1.eta(), momentum2.m(), weight);
    } else {
      LogDebug("Histograms") << "HMassVSPart: wrong charge value = " << charge << std::endl;
      abort();
    }
  }

void HMassVSPart::Write ( )

inlineoverride

Definition at line 1027 of file Histograms.h.

                        {
    hMassVSPt_->Write();
    hMassVSEta_->Write();
    hMassVSPhiPlus_->Write();
    hMassVSPhiMinus_->Write();

    hMassVSEtaPhiPlus_->Write();
    hMassVSEtaPhiMinus_->Write();
    hMassVSEtaPlus_->Write();
    hMassVSEtaMinus_->Write();

    hMassVSPhiPlusPhiMinus_->Write();
    hMassVSEtaPlusEtaMinus_->Write();
    hMassVSCosThetaCS_->Write();
    hMassVSPhiCS_->Write();

    hMassVSPhiPlusMinusDiff_->Write();
    hMassVSEtaPlusMinusDiff_->Write();
    hMassVSCosThetaCS_prof->Write();

    //hMassVSPt_prof_->Write();
    //hMassVSEta_prof_->Write();
    //hMassVSPhiPlus_prof_->Write();
    //hMassVSPhiMinus_prof_->Write();
  }

Member Data Documentation

TH2F* HMassVSPart::hMassVSCosThetaCS_

protected

Definition at line 1083 of file Histograms.h.

TProfile* HMassVSPart::hMassVSCosThetaCS_prof

protected

Definition at line 1097 of file Histograms.h.

TH2F* HMassVSPart::hMassVSEta_

protected

Definition at line 1080 of file Histograms.h.

TH2F* HMassVSPart::hMassVSEtaMinus_

protected

Definition at line 1089 of file Histograms.h.

TH3F* HMassVSPart::hMassVSEtaPhiMinus_

protected

Definition at line 1087 of file Histograms.h.

TH3F* HMassVSPart::hMassVSEtaPhiPlus_

protected

Definition at line 1086 of file Histograms.h.

TH2F* HMassVSPart::hMassVSEtaPlus_

protected

Definition at line 1088 of file Histograms.h.

TH3F* HMassVSPart::hMassVSEtaPlusEtaMinus_

protected

Definition at line 1095 of file Histograms.h.

TH2F* HMassVSPart::hMassVSEtaPlusMinusDiff_

protected

Definition at line 1092 of file Histograms.h.

TH2F* HMassVSPart::hMassVSPhiCS_

protected

Definition at line 1084 of file Histograms.h.

TH2F* HMassVSPart::hMassVSPhiMinus_

protected

Definition at line 1082 of file Histograms.h.

TH2F* HMassVSPart::hMassVSPhiPlus_

protected

Definition at line 1081 of file Histograms.h.

TH2F* HMassVSPart::hMassVSPhiPlusMinusDiff_

protected

Definition at line 1091 of file Histograms.h.

TH3F* HMassVSPart::hMassVSPhiPlusPhiMinus_

protected

Definition at line 1094 of file Histograms.h.

TH2F* HMassVSPart::hMassVSPt_

protected

Definition at line 1079 of file Histograms.h.