CMS 3D CMS Logo

Public Member Functions | Private Member Functions | Private Attributes

BPhysicsOniaDQM Class Reference

#include <BPhysicsOniaDQM.h>

Inheritance diagram for BPhysicsOniaDQM:
edm::EDAnalyzer

List of all members.

Public Member Functions

void analyze (const edm::Event &, const edm::EventSetup &)
 Get the analysis.
void beginJob ()
 Inizialize parameters for histo binning.
void beginLuminosityBlock (const edm::LuminosityBlock &lumiBlock, const edm::EventSetup &iSetup)
void beginRun (const edm::Run &iRun, const edm::EventSetup &iSetup)
 BPhysicsOniaDQM (const edm::ParameterSet &)
 Constructor.
void endJob (void)
void endLuminosityBlock (const edm::LuminosityBlock &lumiBlock, const edm::EventSetup &iSetup)
void endRun (const edm::Run &iRun, const edm::EventSetup &iSetup)
virtual ~BPhysicsOniaDQM ()
 Destructor.

Private Member Functions

float computeMass (const math::XYZVector &vec1, const math::XYZVector &vec2)
bool isMuonInAccept (const reco::Muon &recoMu)
bool selGlobalMuon (const reco::Muon &recoMu)
bool selTrackerMuon (const reco::Muon &recoMu)

Private Attributes

MonitorElementdiMuonMass_global
MonitorElementdiMuonMass_standalone
MonitorElementdiMuonMass_tracker
MonitorElementglbBkgNoCut
MonitorElementglbSigCut
MonitorElementglbSigNoCut
MonitorElementglobal_background
std::map< int, int > jpsiGlbSig
int jpsiGlbSigPerLS
MonitorElementJPsiGlbYdLumi
std::map< int, int > jpsiStaSig
int jpsiStaSigPerLS
MonitorElementJPsiStaYdLumi
std::map< int, int > jpsiTrkSig
int jpsiTrkSigPerLS
MonitorElementJPsiTrkYdLumi
std::string metname
math::XYZPoint RefVtx
MonitorElementstaBkgNoCut
MonitorElementstandalone_background
MonitorElementstaSigCut
MonitorElementstaSigNoCut
DQMStoretheDbe
edm::InputTag theMuonCollectionLabel
MonitorElementtracker_background
MonitorElementtrkBkgNoCut
MonitorElementtrkSigCut
MonitorElementtrkSigNoCut
edm::InputTag vertex

Detailed Description

DQM offline for quarkonia

Date:
2010/11/11 17:33:03
Revision:
1.6
Author:
S. Bolognesi, Eric - CERN

Definition at line 28 of file BPhysicsOniaDQM.h.


Constructor & Destructor Documentation

BPhysicsOniaDQM::BPhysicsOniaDQM ( const edm::ParameterSet parameters)
BPhysicsOniaDQM::~BPhysicsOniaDQM ( ) [virtual]

Destructor.

Definition at line 60 of file BPhysicsOniaDQM.cc.

                                  { 
}

Member Function Documentation

void BPhysicsOniaDQM::analyze ( const edm::Event iEvent,
const edm::EventSetup iSetup 
) [virtual]

Get the analysis.

Implements edm::EDAnalyzer.

Definition at line 92 of file BPhysicsOniaDQM.cc.

References edm::Event::getByLabel(), muon::isGoodMuon(), edm::HandleBase::isValid(), LogTrace, metname, patZpeak::muons, NULL, and muon::TrackerMuonArbitrated.

                                                                           {

  LogTrace(metname)<<"[BPhysicsOniaDQM] Analysis of event # ";
  
  // Take the STA muon container
  Handle<MuonCollection> muons;
  iEvent.getByLabel(theMuonCollectionLabel,muons);

  Handle<reco::VertexCollection> privtxs;
  iEvent.getByLabel(vertex,privtxs);
  VertexCollection::const_iterator privtx;

  if(privtxs->begin() != privtxs->end()){
    privtx = privtxs->begin();
    RefVtx = privtx->position();
  } else {
    RefVtx.SetXYZ(0.,0.,0.);
  }

  if(muons.isValid()){
    for (MuonCollection::const_iterator recoMu1 = muons->begin(); recoMu1!=muons->end(); ++recoMu1){

      // only loop over the remaining muons if recoMu1 is one of the following
      if(recoMu1->isGlobalMuon() || recoMu1->isTrackerMuon() || recoMu1->isStandAloneMuon()){
        for (MuonCollection::const_iterator recoMu2 = recoMu1+1; recoMu2!=muons->end(); ++recoMu2){

          // fill the relevant histograms if recoMu2 satisfies one of the following
          if (recoMu1->isGlobalMuon() && recoMu2->isGlobalMuon()){
            math::XYZVector vec1 = recoMu1->globalTrack()->momentum();
            math::XYZVector vec2 = recoMu2->globalTrack()->momentum();
            float massJPsi = computeMass(vec1,vec2);

            // if opposite charges, fill glbSig, else fill glbBkg
            if (((*recoMu1).charge()*(*recoMu2).charge())<0) {
              if(diMuonMass_global!=NULL){  // BPhysicsOniaDQM original one
                diMuonMass_global->Fill(massJPsi);
              }

              if(glbSigNoCut!=NULL){
                glbSigNoCut->Fill(massJPsi);
                if (selGlobalMuon(*recoMu1) && selGlobalMuon(*recoMu2)) {
                  if (glbSigCut!=NULL) glbSigCut->Fill(massJPsi);
                  if (massJPsi >= 3.0 && massJPsi <= 3.2) jpsiGlbSigPerLS++;
                }
              }
            } else {
              if(global_background!=NULL){  // BPhysicsOniaDQM original one
                global_background->Fill (massJPsi);
              }

              if(glbBkgNoCut!=NULL){
                glbBkgNoCut->Fill(massJPsi);
              }
            }
          }
          
          if(recoMu1->isStandAloneMuon() && recoMu2->isStandAloneMuon() &&
            fabs(recoMu1->outerTrack()->d0()) < 5 && fabs(recoMu1->outerTrack()->dz()) < 30 &&
            fabs(recoMu2->outerTrack()->d0()) < 5 && fabs(recoMu2->outerTrack()->dz()) < 30){
            math::XYZVector vec1 = recoMu1->outerTrack()->momentum();
            math::XYZVector vec2 = recoMu2->outerTrack()->momentum();
            float massJPsi = computeMass(vec1,vec2);

            // if opposite charges, fill staSig, else fill staBkg
            if (((*recoMu1).charge()*(*recoMu2).charge())<0) {
              if(diMuonMass_standalone!=NULL){
                diMuonMass_standalone->Fill(massJPsi);
              }

              if(staSigNoCut!=NULL){
                staSigNoCut->Fill(massJPsi);
                /*if (selStandaloneMuon(*recoMu1) && selStandaloneMuon(*recoMu2)) {
                  if (staSigCut!=NULL) staSigCut->Fill(massJPsi);
                  if (massJPsi >= 3.0 && massJPsi <= 3.2) jpsiStaSigPerLS++;
                }*/
              }
            } else {
              if(standalone_background!=NULL){
                standalone_background->Fill (massJPsi);
              }

              if(staBkgNoCut!=NULL){
                staBkgNoCut->Fill(massJPsi);
              }
            }
          }

          if(recoMu1->isTrackerMuon() && recoMu2->isTrackerMuon() &&
            muon::isGoodMuon(*recoMu1, muon::TrackerMuonArbitrated) &&
            muon::isGoodMuon(*recoMu2, muon::TrackerMuonArbitrated)){
            math::XYZVector vec1 = recoMu1->innerTrack()->momentum();
            math::XYZVector vec2 = recoMu2->innerTrack()->momentum();
            float massJPsi = computeMass(vec1,vec2);

            // if opposite charges, fill trkSig, else fill trkBkg
            if (((*recoMu1).charge()*(*recoMu2).charge())<0) {
              if(diMuonMass_tracker!=NULL){
                diMuonMass_tracker->Fill(massJPsi);
              }

              if(trkSigNoCut!=NULL){
                trkSigNoCut->Fill(massJPsi);
                if (selTrackerMuon(*recoMu1) && selTrackerMuon(*recoMu2)) {
                  if (trkSigCut!=NULL) trkSigCut->Fill(massJPsi);
                  if(massJPsi >= 3.0 && massJPsi <= 3.2) jpsiTrkSigPerLS++;
                }
              }
            } else {
              if(tracker_background!=NULL){
                tracker_background->Fill (massJPsi);
              }

              if(trkBkgNoCut!=NULL){
                trkBkgNoCut->Fill(massJPsi);
              }
            }
          }

        }//end of 2nd MuonCollection
      }//end of GLB,STA,TRK muon check
    }//end of 1st MuonCollection
  }//Is this MuonCollection vaild?

}
void BPhysicsOniaDQM::beginJob ( void  ) [virtual]

Inizialize parameters for histo binning.

Reimplemented from edm::EDAnalyzer.

Definition at line 63 of file BPhysicsOniaDQM.cc.

References LogTrace, metname, NULL, and cmsCodeRules::cppFunctionSkipper::operator.

                               {
  // the services
  theDbe = Service<DQMStore>().operator->();

  metname = "oniaAnalyzer";
  LogTrace(metname)<<"[BPhysicsOniaDQM] Parameters initialization";

  if(theDbe!=NULL){
    theDbe->setCurrentFolder("Physics/BPhysics");  // Use folder with name of PAG
    global_background = theDbe->book1D("global_background", "Same-sign global-global dimuon mass", 750, 0, 15);
    diMuonMass_global = theDbe->book1D("diMuonMass_global", "Opposite-sign global-global dimuon mass", 750, 0, 15);
    tracker_background = theDbe->book1D("tracker_background", "Same-sign tracker-tracker (arbitrated) dimuon mass", 750, 0, 15);
    diMuonMass_tracker = theDbe->book1D("diMuonMass_tracker", "Opposite-sign tracker-tracker (arbitrated) dimuon mass", 750, 0, 15);
    standalone_background = theDbe->book1D("standalone_background", "Same-sign standalone-standalone dimuon mass", 500, 0, 15);
    diMuonMass_standalone = theDbe->book1D("diMuonMass_standalone", "Opposite-sign standalone-standalone dimuon mass", 500, 0, 15);

    glbSigCut = theDbe->book1D("glbSigCut", "Opposite-sign glb-glb dimuon mass", 650, 0, 130);
    glbSigNoCut = theDbe->book1D("glbSigNoCut", "Opposite-sign glb-glb dimuon mass (no cut)", 650, 0, 130);
    glbBkgNoCut = theDbe->book1D("glbBkgNoCut", "Same-sign glb-glb dimuon mass (no cut)", 650, 0, 130);
    staSigCut = theDbe->book1D("staSigCut", "Opposite-sign sta-sta dimuon mass", 430, 0, 129);
    staSigNoCut = theDbe->book1D("staSigNoCut", "Opposite-sign sta-sta dimuon mass (no cut)", 430, 0, 129);
    staBkgNoCut  = theDbe->book1D("staBkgNoCut", "Same-sign sta-sta dimuon mass (no cut)", 430, 0, 129);
    trkSigCut = theDbe->book1D("trkSigCut", "Opposite-sign trk-trk dimuon mass", 650, 0, 130);
    trkSigNoCut = theDbe->book1D("trkSigNoCut", "Opposite-sign trk-trk dimuon mass (no cut)", 650, 0, 130);
    trkBkgNoCut = theDbe->book1D("trkBkgNoCutt", "Same-sign trk-trk dimuon mass (no cut)", 650, 0, 130);
  }

}
void BPhysicsOniaDQM::beginLuminosityBlock ( const edm::LuminosityBlock lumiBlock,
const edm::EventSetup iSetup 
) [virtual]

Reimplemented from edm::EDAnalyzer.

Definition at line 221 of file BPhysicsOniaDQM.cc.

References LogTrace, and metname.

{
  LogTrace(metname)<<"[BPhysicsOniaDQM] Start of a LuminosityBlock";
  
  jpsiGlbSigPerLS = 0;
  jpsiStaSigPerLS = 0;
  jpsiTrkSigPerLS = 0;
}
void BPhysicsOniaDQM::beginRun ( const edm::Run iRun,
const edm::EventSetup iSetup 
) [virtual]

Reimplemented from edm::EDAnalyzer.

Definition at line 279 of file BPhysicsOniaDQM.cc.

References LogTrace, and metname.

{
  LogTrace(metname)<<"[BPhysicsOniaDQM] Start of a Run";
}
float BPhysicsOniaDQM::computeMass ( const math::XYZVector vec1,
const math::XYZVector vec2 
) [private]

Definition at line 295 of file BPhysicsOniaDQM.cc.

References mathSSE::sqrt().

                                                                                     {
  // mass of muon
  float massMu = 0.10566;
  float eMu1 = -999;
  if(massMu*massMu + vec1.Mag2()>0)
    eMu1 = sqrt(massMu*massMu + vec1.Mag2());
  float eMu2 = -999;
  if(massMu*massMu + vec2.Mag2()>0)
    eMu2 = sqrt(massMu*massMu + vec2.Mag2());

  float pJPsi = -999;
  if((vec1+vec2).Mag2()>0)
    pJPsi = sqrt((vec1+vec2).Mag2());
  float eJPsi = eMu1 + eMu2;

  float massJPsi = -999;
  if((eJPsi*eJPsi - pJPsi*pJPsi) > 0)
    massJPsi = sqrt(eJPsi*eJPsi - pJPsi*pJPsi);
 
 return massJPsi;
}
void BPhysicsOniaDQM::endJob ( void  ) [virtual]

Reimplemented from edm::EDAnalyzer.

Definition at line 217 of file BPhysicsOniaDQM.cc.

References LogTrace, and metname.

                                 {
  LogTrace(metname)<<"[BPhysicsOniaDQM] EndJob";
}
void BPhysicsOniaDQM::endLuminosityBlock ( const edm::LuminosityBlock lumiBlock,
const edm::EventSetup iSetup 
) [virtual]

Reimplemented from edm::EDAnalyzer.

Definition at line 230 of file BPhysicsOniaDQM.cc.

References newFWLiteAna::bin, first, edm::LuminosityBlock::getByLabel(), glb, edm::LuminosityBlockBase::id(), LogTrace, edm::LuminosityBlockID::luminosityBlock(), runregparse::lumiSummary, metname, and NULL.

{
  LogTrace(metname)<<"[BPhysicsOniaDQM] Start of a LuminosityBlock";

  edm::Handle<LumiSummary> lumiSummary;
  lumiBlock.getByLabel("lumiProducer",lumiSummary);

  int LBlockNum = lumiBlock.id().luminosityBlock();
  
  jpsiGlbSig.insert( pair<int,int>(LBlockNum, jpsiGlbSigPerLS) );
  jpsiStaSig.insert( pair<int,int>(LBlockNum, jpsiStaSigPerLS) );
  jpsiTrkSig.insert( pair<int,int>(LBlockNum, jpsiTrkSigPerLS) );
//  cout << "lumi: " << LBlockNum << "\t" << jpsiGlbSig[LBlockNum] << "\t" << jpsiStaSig[LBlockNum] << "\t" << jpsiTrkSig[LBlockNum] << endl;

  if (jpsiGlbSig.size()%5 != 0) return;

  theDbe->setCurrentFolder("Physics/BPhysics");
  if(JPsiGlbYdLumi!=NULL) {
    theDbe->removeElement("JPsiGlbYdLumi");   // Remove histograms from previous run
    theDbe->removeElement("JPsiStaYdLumi");
    theDbe->removeElement("JPsiTrkYdLumi");
  }

  int xmin = (*jpsiGlbSig.begin()).first;
  int xmax = (*jpsiGlbSig.rbegin()).first;
  int nx   = (xmax - xmin + 1)/5 + 1; // Merge 5 lumisections into 1 bin
//  cout << "x-axis " << xmin << " " << xmax << endl;

  JPsiGlbYdLumi = theDbe->book1D("JPsiGlbYdLumi", "JPsi yield from global-global dimuon", nx, xmin, xmax);
  JPsiStaYdLumi = theDbe->book1D("JPsiStaYdLumi", "JPsi yield from standalone-standalone dimuon", nx, xmin, xmax);
  JPsiTrkYdLumi = theDbe->book1D("JPsiTrkYdLumi", "JPsi yield from tracker-tracker dimuon", nx, xmin, xmax);

  map<int,int>::iterator glb;
  map<int,int>::iterator sta;
  map<int,int>::iterator trk;
  for (glb = jpsiGlbSig.begin(); glb != jpsiGlbSig.end(); ++glb)
  {
    int bin = ((*glb).first - xmin + 1)/5 + 1;  //X-axis bin #
    sta = jpsiStaSig.find((*glb).first);
    trk = jpsiTrkSig.find((*glb).first);
    JPsiGlbYdLumi->setBinContent(bin,JPsiGlbYdLumi->getBinContent(bin)+(*glb).second);
    JPsiStaYdLumi->setBinContent(bin,JPsiStaYdLumi->getBinContent(bin)+(*sta).second);
    JPsiTrkYdLumi->setBinContent(bin,JPsiTrkYdLumi->getBinContent(bin)+(*trk).second);
//    cout << "glb: " << bin << "\t" << (*glb).first << "\t" << (*glb).second << endl;
//    cout << "sta: " << bin << "\t" << (*sta).first << "\t" << (*sta).second << endl;
//    cout << "trk: " << bin << "\t" << (*trk).first << "\t" << (*trk).second << endl;
  }
}
void BPhysicsOniaDQM::endRun ( const edm::Run iRun,
const edm::EventSetup iSetup 
) [virtual]

Reimplemented from edm::EDAnalyzer.

Definition at line 284 of file BPhysicsOniaDQM.cc.

References LogTrace, and metname.

{
  LogTrace(metname)<<"[BPhysicsOniaDQM] End of a Run";
  
  if (!jpsiGlbSig.empty()) {
    jpsiGlbSig.clear();
    jpsiStaSig.clear();
    jpsiTrkSig.clear();
  }
}
bool BPhysicsOniaDQM::isMuonInAccept ( const reco::Muon recoMu) [private]

Definition at line 317 of file BPhysicsOniaDQM.cc.

References reco::LeafCandidate::eta(), reco::LeafCandidate::p(), and reco::LeafCandidate::pt().

{
  return (fabs(recoMu.eta()) < 2.4 &&
         ((fabs(recoMu.eta()) < 1.3 && recoMu.pt() > 3.3) ||
          (fabs(recoMu.eta()) > 1.3 && fabs(recoMu.eta()) < 2.2 && recoMu.p() > 2.9) ||
          (fabs(recoMu.eta()) > 2.2 && recoMu.pt() > 0.8)));
}
bool BPhysicsOniaDQM::selGlobalMuon ( const reco::Muon recoMu) [private]

Definition at line 325 of file BPhysicsOniaDQM.cc.

References reco::Muon::globalTrack(), reco::Muon::innerTrack(), reco::HitPattern::numberOfValidMuonHits(), AlCaHLTBitMon_ParallelJobs::p, reco::HitPattern::pixelLayersWithMeasurement(), and lumiQueryAPI::q.

{
  TrackRef iTrack = recoMu.innerTrack();
  const reco::HitPattern &p = iTrack->hitPattern();
  
  TrackRef gTrack = recoMu.globalTrack();
  const reco::HitPattern &q = gTrack->hitPattern();

  return (isMuonInAccept(recoMu) &&
          iTrack->found() > 11 &&
          gTrack->chi2()/gTrack->ndof() < 20.0 &&
          q.numberOfValidMuonHits() > 0 &&
          iTrack->chi2()/iTrack->ndof() < 4.0 &&
          //recoMu.muonID("TrackerMuonArbitrated") &&
          //recoMu.muonID("TMLastStationAngTight") &&
          p.pixelLayersWithMeasurement() > 1 &&
          fabs(iTrack->dxy(RefVtx)) < 3.0 &&
          fabs(iTrack->dz(RefVtx)) < 15.0 );
}
bool BPhysicsOniaDQM::selTrackerMuon ( const reco::Muon recoMu) [private]

Definition at line 345 of file BPhysicsOniaDQM.cc.

References reco::Muon::innerTrack(), AlCaHLTBitMon_ParallelJobs::p, and reco::HitPattern::pixelLayersWithMeasurement().

{
  TrackRef iTrack = recoMu.innerTrack();
  const reco::HitPattern &p = iTrack->hitPattern();

  return (isMuonInAccept(recoMu) &&
          iTrack->found() > 11 &&
          iTrack->chi2()/iTrack->ndof() < 4.0 &&
          //recoMu.muonID("TrackerMuonArbitrated") &&
          //recoMu.muonID("TMLastStationAngTight") &&
          p.pixelLayersWithMeasurement() > 1 &&
          fabs(iTrack->dxy(RefVtx)) < 3.0 &&
          fabs(iTrack->dz(RefVtx)) < 15.0 );
}

Member Data Documentation

Definition at line 67 of file BPhysicsOniaDQM.h.

Definition at line 69 of file BPhysicsOniaDQM.h.

Definition at line 68 of file BPhysicsOniaDQM.h.

Definition at line 76 of file BPhysicsOniaDQM.h.

Definition at line 74 of file BPhysicsOniaDQM.h.

Definition at line 75 of file BPhysicsOniaDQM.h.

Definition at line 70 of file BPhysicsOniaDQM.h.

std::map<int,int> BPhysicsOniaDQM::jpsiGlbSig [private]

Definition at line 92 of file BPhysicsOniaDQM.h.

Definition at line 89 of file BPhysicsOniaDQM.h.

Definition at line 84 of file BPhysicsOniaDQM.h.

std::map<int,int> BPhysicsOniaDQM::jpsiStaSig [private]

Definition at line 93 of file BPhysicsOniaDQM.h.

Definition at line 90 of file BPhysicsOniaDQM.h.

Definition at line 85 of file BPhysicsOniaDQM.h.

std::map<int,int> BPhysicsOniaDQM::jpsiTrkSig [private]

Definition at line 94 of file BPhysicsOniaDQM.h.

Definition at line 91 of file BPhysicsOniaDQM.h.

Definition at line 86 of file BPhysicsOniaDQM.h.

std::string BPhysicsOniaDQM::metname [private]

Definition at line 61 of file BPhysicsOniaDQM.h.

Definition at line 96 of file BPhysicsOniaDQM.h.

Definition at line 79 of file BPhysicsOniaDQM.h.

Definition at line 72 of file BPhysicsOniaDQM.h.

Definition at line 77 of file BPhysicsOniaDQM.h.

Definition at line 78 of file BPhysicsOniaDQM.h.

Definition at line 57 of file BPhysicsOniaDQM.h.

Definition at line 64 of file BPhysicsOniaDQM.h.

Definition at line 71 of file BPhysicsOniaDQM.h.

Definition at line 82 of file BPhysicsOniaDQM.h.

Definition at line 80 of file BPhysicsOniaDQM.h.

Definition at line 81 of file BPhysicsOniaDQM.h.

Definition at line 59 of file BPhysicsOniaDQM.h.