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 792 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 795 of file Histograms.h.

References PFRecoTauDiscriminationByHPSSelection_cfi::maxMass, MuonErrorMatrixAnalyzer_cfi::maxPt, PFRecoTauDiscriminationByHPSSelection_cfi::minMass, and dataset::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 837 of file Histograms.h.

References dataset::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 860 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 1012 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 878 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 884 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 896 of file Histograms.h.

References RPCpg::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 961 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 986 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 1042 of file Histograms.h.

TProfile* HMassVSPart::hMassVSCosThetaCS_prof

protected

Definition at line 1056 of file Histograms.h.

TH2F* HMassVSPart::hMassVSEta_

protected

Definition at line 1039 of file Histograms.h.

TH2F* HMassVSPart::hMassVSEtaMinus_

protected

Definition at line 1048 of file Histograms.h.

TH3F* HMassVSPart::hMassVSEtaPhiMinus_

protected

Definition at line 1046 of file Histograms.h.

TH3F* HMassVSPart::hMassVSEtaPhiPlus_

protected

Definition at line 1045 of file Histograms.h.

TH2F* HMassVSPart::hMassVSEtaPlus_

protected

Definition at line 1047 of file Histograms.h.

TH3F* HMassVSPart::hMassVSEtaPlusEtaMinus_

protected

Definition at line 1054 of file Histograms.h.

TH2F* HMassVSPart::hMassVSEtaPlusMinusDiff_

protected

Definition at line 1051 of file Histograms.h.

TH2F* HMassVSPart::hMassVSPhiCS_

protected

Definition at line 1043 of file Histograms.h.

TH2F* HMassVSPart::hMassVSPhiMinus_

protected

Definition at line 1041 of file Histograms.h.

TH2F* HMassVSPart::hMassVSPhiPlus_

protected

Definition at line 1040 of file Histograms.h.

TH2F* HMassVSPart::hMassVSPhiPlusMinusDiff_

protected

Definition at line 1050 of file Histograms.h.

TH3F* HMassVSPart::hMassVSPhiPlusPhiMinus_

protected

Definition at line 1053 of file Histograms.h.

TH2F* HMassVSPart::hMassVSPt_

protected

Definition at line 1038 of file Histograms.h.