Skip to content

Commit

Permalink
Merge pull request #9010 from gzevi/dqmoffline_egamma_updates
Browse files Browse the repository at this point in the history 
Changes to Egamma/Electrons plots: separate +/- endcap for alignment,…
  • Loading branch information
@cmsbuild
cmsbuild committed May 11, 2015
2 parents bf6ed20 + c6f0d07 commit 0c17f04
Show file tree
Hide file tree
Showing 2 changed files with 51 additions and 14 deletions.
48 changes: 38 additions & 10 deletions DQMOffline/EGamma/plugins/ElectronAnalyzer.cc
View file
Expand Up @@ -192,20 +192,26 @@ void ElectronAnalyzer::bookHistograms( DQMStore::IBooker & iBooker, edm::Run con
h1_EeleOPout_barrel = bookH1(iBooker, "EeleOPout_barrel","ele E_{ele}/P_{out} in barrel",nbineop,0.,eopmax,"E_{ele}/P_{out}");
h1_EeleOPout_endcaps = bookH1(iBooker, "EeleOPout_endcaps","ele E_{ele}/P_{out} in endcaps",nbineop,0.,eopmax,"E_{ele}/P_{out}");
h1_dEtaSc_propVtx_barrel = bookH1(iBooker, "dEtaSc_propVtx_barrel","ele #eta_{sc} - #eta_{tr}, prop from vertex, in barrel",nbindetamatch,detamatchmin,detamatchmax,"#eta_{sc} - #eta_{tr}");
h1_dEtaSc_propVtx_endcaps = bookH1(iBooker, "dEtaSc_propVtx_endcaps","ele #eta_{sc} - #eta_{tr}, prop from vertex, in endcaps",nbindetamatch,detamatchmin,detamatchmax,"#eta_{sc} - #eta_{tr}");
h1_dEtaSc_propVtx_endcapsPos = bookH1(iBooker, "dEtaSc_propVtx_endcapsPos","ele #eta_{sc} - #eta_{tr}, prop from vertex, in positive endcap",nbindetamatch,detamatchmin,detamatchmax,"#eta_{sc} - #eta_{tr}");
h1_dEtaSc_propVtx_endcapsNeg = bookH1(iBooker, "dEtaSc_propVtx_endcapsNeg","ele #eta_{sc} - #eta_{tr}, prop from vertex, in negative endcap",nbindetamatch,detamatchmin,detamatchmax,"#eta_{sc} - #eta_{tr}");
py_dEtaSc_propVtxVsPhi = bookP1(iBooker, "dEtaSc_propVtxVsPhi","ele #eta_{sc} - #eta_{tr}, prop from vertex vs #phi",nbinphi2D,phimin,phimax,detamatchmin,detamatchmax,"#phi (rad)","<#eta_{sc} - #eta_{tr}>");
h1_dEtaEleCl_propOut_barrel = bookH1(iBooker, "dEtaEleCl_propOut_barrel","ele #eta_{EleCl} - #eta_{tr}, prop from outermost, in barrel",nbindetamatch,detamatchmin,detamatchmax,"#eta_{elecl} - #eta_{tr}");
h1_dEtaEleCl_propOut_endcaps = bookH1(iBooker, "dEtaEleCl_propOut_endcaps","ele #eta_{EleCl} - #eta_{tr}, prop from outermost, in endcaps",nbindetamatch,detamatchmin,detamatchmax,"#eta_{elecl} - #eta_{tr}");
h1_dEtaEleCl_propOut_endcapsPos = bookH1(iBooker, "dEtaEleCl_propOut_endcapsPos","ele #eta_{EleCl} - #eta_{tr}, prop from outermost, in positive endcap",nbindetamatch,detamatchmin,detamatchmax,"#eta_{elecl} - #eta_{tr}");
h1_dEtaEleCl_propOut_endcapsNeg = bookH1(iBooker, "dEtaEleCl_propOut_endcapsNeg","ele #eta_{EleCl} - #eta_{tr}, prop from outermost, in negative endcap",nbindetamatch,detamatchmin,detamatchmax,"#eta_{elecl} - #eta_{tr}");
h1_dPhiSc_propVtx_barrel = bookH1(iBooker, "dPhiSc_propVtx_barrel","ele #phi_{sc} - #phi_{tr}, prop from vertex, in barrel",nbindphimatch,dphimatchmin,dphimatchmax,"#phi_{sc} - #phi_{tr} (rad)");
h1_dPhiSc_propVtx_endcaps = bookH1(iBooker, "dPhiSc_propVtx_endcaps","ele #phi_{sc} - #phi_{tr}, prop from vertex, in endcaps",nbindphimatch,dphimatchmin,dphimatchmax,"#phi_{sc} - #phi_{tr} (rad)");
h1_dPhiSc_propVtx_endcapsPos = bookH1(iBooker, "dPhiSc_propVtx_endcapsPos","ele #phi_{sc} - #phi_{tr}, prop from vertex, in positive endcap",nbindphimatch,dphimatchmin,dphimatchmax,"#phi_{sc} - #phi_{tr} (rad)");
h1_dPhiSc_propVtx_endcapsNeg = bookH1(iBooker, "dPhiSc_propVtx_endcapsNeg","ele #phi_{sc} - #phi_{tr}, prop from vertex, in negative endcap",nbindphimatch,dphimatchmin,dphimatchmax,"#phi_{sc} - #phi_{tr} (rad)");
py_dPhiSc_propVtxVsPhi = bookP1(iBooker, "dPhiSc_propVtxVsPhi","ele #phi_{sc} - #phi_{tr}, prop from vertex vs #phi",nbinphi2D,phimin,phimax,dphimatchmin,dphimatchmax,"#phi (rad)","<#phi_{sc} - #phi_{tr}> (rad)");
h1_dPhiEleCl_propOut_barrel = bookH1(iBooker, "dPhiEleCl_propOut_barrel","ele #phi_{EleCl} - #phi_{tr}, prop from outermost, in barrel",nbindphimatch,dphimatchmin,dphimatchmax,"#phi_{elecl} - #phi_{tr} (rad)");
h1_dPhiEleCl_propOut_endcaps = bookH1(iBooker, "dPhiEleCl_propOut_endcaps","ele #phi_{EleCl} - #phi_{tr}, prop from outermost, in endcaps",nbindphimatch,dphimatchmin,dphimatchmax,"#phi_{elecl} - #phi_{tr} (rad)");
h1_dPhiEleCl_propOut_endcapsPos = bookH1(iBooker, "dPhiEleCl_propOut_endcapsPos","ele #phi_{EleCl} - #phi_{tr}, prop from outermost, in positive endcap",nbindphimatch,dphimatchmin,dphimatchmax,"#phi_{elecl} - #phi_{tr} (rad)");
h1_dPhiEleCl_propOut_endcapsNeg = bookH1(iBooker, "dPhiEleCl_propOut_endcapsNeg","ele #phi_{EleCl} - #phi_{tr}, prop from outermost, in negative endcap",nbindphimatch,dphimatchmin,dphimatchmax,"#phi_{elecl} - #phi_{tr} (rad)");
h1_Hoe_barrel = bookH1(iBooker, "Hoe_barrel","ele hadronic energy / em energy, in barrel", nbinhoe, hoemin, hoemax,"H/E","Events","ELE_LOGY E1 P") ;
h1_Hoe_endcaps = bookH1(iBooker, "Hoe_endcaps","ele hadronic energy / em energy, in endcaps", nbinhoe, hoemin, hoemax,"H/E","Events","ELE_LOGY E1 P") ;
py_HoeVsPhi = bookP1(iBooker, "HoeVsPhi","ele hadronic energy / em energy vs #phi",nbinphi2D,phimin,phimax,hoemin,hoemax,"#phi (rad)","<H/E>","E1 P") ;
h1_sclSigEtaEta_barrel = bookH1(iBooker, "sclSigEtaEta_barrel","ele supercluster sigma ieta ieta in barrel",100,0.,0.05,"sietaieta");
h1_sclSigEtaEta_endcaps = bookH1(iBooker, "sclSigEtaEta_endcaps","ele supercluster sigma ieta ieta in endcaps",100,0.,0.05,"sietaieta");
h1_sclSigEtaEta_barrel = bookH1(iBooker, "sclSigEtaEta_barrel","ele sigma eta eta in barrel",100,0.,0.05,"sietaieta");
h1_sclSigEtaEta_endcaps = bookH1(iBooker, "sclSigEtaEta_endcaps","ele sigma eta eta in endcaps",100,0.,0.05,"sietaieta");
h1_sigIEtaIEta5x5_barrel = bookH1(iBooker, "sigIEtaIEta5x5_barrel","ele sigma ieta ieta 5x5 in barrel",100,0.,0.05,"sietaieta5x5");
h1_sigIEtaIEta5x5_endcaps = bookH1(iBooker, "sigIEtaIEta5x5_endcaps","ele sigma ieta ieta 5x5 in endcaps",100,0.,0.05,"sietaieta5x5");

// fbrem
h1_fbrem = bookH1(iBooker, "fbrem","ele brem fraction",100,0.,1.,"P_{in} - P_{out} / P_{in}") ;
Expand All @@ -222,6 +228,11 @@ void ElectronAnalyzer::bookHistograms( DQMStore::IBooker & iBooker, edm::Run con
h1_ecalRecHitSumEt_dr03 = bookH1(iBooker, "ecalRecHitSumEt_dr03","ecal isolation sum, dR=0.3",100,0.0,20.,"EcalIsoSum (GeV)","Events","ELE_LOGY E1 P");
h1_hcalTowerSumEt_dr03 = bookH1(iBooker, "hcalTowerSumEt_dr03","hcal isolation sum, dR=0.3",100,0.0,20.,"HcalIsoSum (GeV)","Events","ELE_LOGY E1 P");

// pf isolation
h1_PFch_dr03 = bookH1(iBooker, "PFch_dr03","Charged PF candidate sum, dR=0.3",100,0.0,20.,"PF charged (GeV/c)","Events","ELE_LOGY E1 P");
h1_PFem_dr03 = bookH1(iBooker, "PFem_dr03","Neutral EM PF candidate sum, dR=0.3",100,0.0,20.,"PF neutral EM (GeV)","Events","ELE_LOGY E1 P");
h1_PFnh_dr03 = bookH1(iBooker, "PFnh_dr03","Neutral Had PF candidate sum, dR=0.3",100,0.0,20.,"PF neutral Had (GeV)","Events","ELE_LOGY E1 P");

// di-electron mass
setBookIndex(200) ;
h1_mee = bookH1(iBooker, "mee","ele pairs invariant mass", nbinmee, meemin, meemax,"m_{ee} (GeV/c^{2})");
Expand Down Expand Up @@ -339,17 +350,27 @@ void ElectronAnalyzer::analyze( const edm::Event& iEvent, const edm::EventSetup
h1_dPhiEleCl_propOut_barrel->Fill(gsfIter->deltaPhiEleClusterTrackAtCalo());
h1_Hoe_barrel->Fill(gsfIter->hadronicOverEm());
h1_sclSigEtaEta_barrel->Fill( gsfIter->scSigmaEtaEta() );
h1_sigIEtaIEta5x5_barrel->Fill( gsfIter->full5x5_sigmaIetaIeta() );
}
if (gsfIter->isEE())
{
h1_Eop_endcaps->Fill( gsfIter->eSuperClusterOverP() );
h1_EeleOPout_endcaps->Fill( gsfIter->eEleClusterOverPout() );
h1_dEtaSc_propVtx_endcaps->Fill(gsfIter->deltaEtaSuperClusterTrackAtVtx());
h1_dEtaEleCl_propOut_endcaps->Fill(gsfIter->deltaEtaEleClusterTrackAtCalo());
h1_dPhiSc_propVtx_endcaps->Fill(gsfIter->deltaPhiSuperClusterTrackAtVtx());
h1_dPhiEleCl_propOut_endcaps->Fill(gsfIter->deltaPhiEleClusterTrackAtCalo());
if (gsfIter->eta() > 0) {
h1_dEtaSc_propVtx_endcapsPos->Fill(gsfIter->deltaEtaSuperClusterTrackAtVtx());
h1_dEtaEleCl_propOut_endcapsPos->Fill(gsfIter->deltaEtaEleClusterTrackAtCalo());
h1_dPhiSc_propVtx_endcapsPos->Fill(gsfIter->deltaPhiSuperClusterTrackAtVtx());
h1_dPhiEleCl_propOut_endcapsPos->Fill(gsfIter->deltaPhiEleClusterTrackAtCalo());
}
else {
h1_dEtaSc_propVtx_endcapsNeg->Fill(gsfIter->deltaEtaSuperClusterTrackAtVtx());
h1_dEtaEleCl_propOut_endcapsNeg->Fill(gsfIter->deltaEtaEleClusterTrackAtCalo());
h1_dPhiSc_propVtx_endcapsNeg->Fill(gsfIter->deltaPhiSuperClusterTrackAtVtx());
h1_dPhiEleCl_propOut_endcapsNeg->Fill(gsfIter->deltaPhiEleClusterTrackAtCalo());
}
h1_Hoe_endcaps->Fill(gsfIter->hadronicOverEm());
h1_sclSigEtaEta_endcaps->Fill( gsfIter->scSigmaEtaEta() );
h1_sigIEtaIEta5x5_endcaps->Fill( gsfIter->full5x5_sigmaIetaIeta() );
}
py_EopVsPhi->Fill( gsfIter->phi(), gsfIter->eSuperClusterOverP() );
py_dEtaSc_propVtxVsPhi->Fill(gsfIter->phi(), gsfIter->deltaEtaSuperClusterTrackAtVtx());
Expand Down Expand Up @@ -378,6 +399,13 @@ void ElectronAnalyzer::analyze( const edm::Event& iEvent, const edm::EventSetup
h1_ecalRecHitSumEt_dr03->Fill(gsfIter->dr03EcalRecHitSumEt());
h1_hcalTowerSumEt_dr03->Fill(gsfIter->dr03HcalTowerSumEt());

// PF isolation
GsfElectron::PflowIsolationVariables pfIso = gsfIter->pfIsolationVariables();
h1_PFch_dr03->Fill( pfIso.sumChargedHadronPt );
h1_PFem_dr03->Fill( pfIso.sumPhotonEt );
h1_PFnh_dr03->Fill( pfIso.sumNeutralHadronEt );


}

// association matching object-reco electrons
Expand Down
17 changes: 13 additions & 4 deletions DQMOffline/EGamma/plugins/ElectronAnalyzer.h
View file
Expand Up @@ -229,25 +229,29 @@ class ElectronAnalyzer : public ElectronDqmAnalyzerBase
// MonitorElement * h1_dEtaSc_propVtx ;
// MonitorElement * h2_dEtaSc_propVtxVsEta ;
MonitorElement * h1_dEtaSc_propVtx_barrel ;
MonitorElement * h1_dEtaSc_propVtx_endcaps ;
MonitorElement * h1_dEtaSc_propVtx_endcapsPos ;
MonitorElement * h1_dEtaSc_propVtx_endcapsNeg ;
MonitorElement * py_dEtaSc_propVtxVsPhi ;
// MonitorElement * h2_dEtaSc_propVtxVsPt ;
// MonitorElement * h1_dEtaEleCl_propOut ;
// MonitorElement * h2_dEtaEleCl_propOutVsEta ;
MonitorElement * h1_dEtaEleCl_propOut_barrel ;
MonitorElement * h1_dEtaEleCl_propOut_endcaps ;
MonitorElement * h1_dEtaEleCl_propOut_endcapsPos ;
MonitorElement * h1_dEtaEleCl_propOut_endcapsNeg ;
// MonitorElement * h2_dEtaEleCl_propOutVsPhi ;
// MonitorElement * h2_dEtaEleCl_propOutVsPt ;
// MonitorElement * h1_dPhiSc_propVtx ;
// MonitorElement * h2_dPhiSc_propVtxVsEta ;
MonitorElement * h1_dPhiSc_propVtx_barrel ;
MonitorElement * h1_dPhiSc_propVtx_endcaps ;
MonitorElement * h1_dPhiSc_propVtx_endcapsPos ;
MonitorElement * h1_dPhiSc_propVtx_endcapsNeg ;
MonitorElement * py_dPhiSc_propVtxVsPhi ;
// MonitorElement * h2_dPhiSc_propVtxVsPt ;
// MonitorElement * h1_dPhiEleCl_propOut ;
// MonitorElement * h2_dPhiEleCl_propOutVsEta ;
MonitorElement * h1_dPhiEleCl_propOut_barrel ;
MonitorElement * h1_dPhiEleCl_propOut_endcaps ;
MonitorElement * h1_dPhiEleCl_propOut_endcapsPos ;
MonitorElement * h1_dPhiEleCl_propOut_endcapsNeg ;
// MonitorElement * h2_dPhiEleCl_propOutVsPhi ;
// MonitorElement * h2_dPhiEleCl_propOutVsPt ;
// MonitorElement * h1_Hoe ;
Expand All @@ -258,6 +262,8 @@ class ElectronAnalyzer : public ElectronDqmAnalyzerBase
// MonitorElement * h2_HoeVsPt ;
MonitorElement * h1_sclSigEtaEta_barrel ;
MonitorElement * h1_sclSigEtaEta_endcaps ;
MonitorElement * h1_sigIEtaIEta5x5_barrel ;
MonitorElement * h1_sigIEtaIEta5x5_endcaps ;

// fbrem related variables
//MonitorElement * h_outerP ;
Expand All @@ -281,6 +287,9 @@ class ElectronAnalyzer : public ElectronDqmAnalyzerBase
MonitorElement * h1_tkSumPt_dr03 ;
MonitorElement * h1_ecalRecHitSumEt_dr03 ;
MonitorElement * h1_hcalTowerSumEt_dr03 ;
MonitorElement * h1_PFch_dr03 ;
MonitorElement * h1_PFem_dr03 ;
MonitorElement * h1_PFnh_dr03 ;
// MonitorElement * h1_hcalDepth1TowerSumEt_dr03 ;
// MonitorElement * h1_hcalDepth2TowerSumEt_dr03 ;
// MonitorElement * h1_tkSumPt_dr04 ;
Expand Down

0 comments on commit 0c17f04

Please sign in to comment.