Skip to content

Commit

Permalink
Updating plotting scripts
Browse files Browse the repository at this point in the history
  • Loading branch information
@naodell
naodell committed Aug 21, 2014
1 parent dd16623 commit 7c7c73a
Show file tree
Hide file tree
Showing 6 changed files with 161 additions and 460 deletions.
48 changes: 24 additions & 24 deletions fcncAnalysis/batch_t3_cfg.py
View file
Original file line number Diff line number Diff line change
Expand Up @@ -65,10 +65,10 @@

bg.extend([
#cfg('ZJetsToTauTau_M-50', '{0}/DYJetsToLL_M-50'.format(mcDir), 60, 'DYJetsToTauTau {0} 2012'.format(mcTrigger)),
#cfg('ZJets_M-50', '{0}/DYJetsToLL_M-50'.format(mcDir), 30, 'ZJets_M-50 {0} 2012'.format(mcTrigger)),
#cfg('ZJets_M-10To50', '{0}/DYJetsToLL_M-10To50filter'.format(mcDir), 20, 'ZJets_M-10To50 {0} 2012'.format(mcTrigger)),
cfg('ZJets_M-50', '{0}/DYJetsToLL_M-50'.format(mcDir), 30, 'ZJets_M-50 {0} 2012'.format(mcTrigger)),
cfg('ZJets_M-10To50', '{0}/DYJetsToLL_M-10To50filter'.format(mcDir), 20, 'ZJets_M-10To50 {0} 2012'.format(mcTrigger)),
#cfg('WbbToLNu', '{0}/WbbJetsToLNu'.format(mcDir), 20, 'WbbToLNu {0} 2012'.format(mcTrigger)),
#cfg('WjetToLNu', '{0}/WJetsToLNu'.format(mcDir), 50, 'WJetsToLNu {0} 2012'.format(mcTrigger)),
cfg('WjetToLNu', '{0}/WJetsToLNu'.format(mcDir), 50, 'WJetsToLNu {0} 2012'.format(mcTrigger)),
cfg('WGStarLNu2Mu', '{0}/WGstarToLNu2Mu'.format(mcDir), 5, 'WGStarLNu2Mu {0} 2012'.format(mcTrigger)),
cfg('WGStarLNu2Tau', '{0}/WGstarToLNu2Tau'.format(mcDir), 5, 'WGStarLNu2Tau {0} 2012'.format(mcTrigger)),

Expand All @@ -77,26 +77,26 @@
#cfg('ZG', mcDir+'/ZGToLLG', 10, 'ZG muon 2012'),
cfg('WG', '{0}/WGToLNuG'.format('/tthome/naodell/storage/ntuples/MC_skimmed'), 10, 'WG {0} 2012'.format(mcTrigger)),

#cfg('ttbarHad', '{0}/TTJets'.format(mcDir), 30, 'ttbarHad {0} 2012'.format(mcTrigger)),
#cfg('ttbarLep', '{0}/TTJets'.format(mcDir), 30, 'ttbarLep {0} 2012'.format(mcTrigger)),
#cfg('tbarW', '{0}/Tbar_tW'.format(mcDir), 10, 'tbarW {0} 2012'.format(mcTrigger)),
#cfg('tW', '{0}/T_tW'.format(mcDir), 10, 'tW {0} 2012'.format(mcTrigger)),
cfg('ttbarHad', '{0}/TTJets'.format(mcDir), 30, 'ttbarHad {0} 2012'.format(mcTrigger)),
cfg('ttbarLep', '{0}/TTJets'.format(mcDir), 30, 'ttbarLep {0} 2012'.format(mcTrigger)),
cfg('tbarW', '{0}/Tbar_tW'.format(mcDir), 10, 'tbarW {0} 2012'.format(mcTrigger)),
cfg('tW', '{0}/T_tW'.format(mcDir), 10, 'tW {0} 2012'.format(mcTrigger)),
cfg('ttW', '{0}/TTWJets'.format(mcDir), 5, 'ttW {0} 2012'.format(mcTrigger)),
cfg('ttZ', '{0}/TTZJets'.format(mcDir), 5, 'ttZ {0} 2012'.format(mcTrigger)),
cfg('ttG', '{0}/TTGJets'.format(mcDir), 5, 'ttG {0} 2012'.format(mcTrigger)),
#cfg('t_t-channel', '{0}/T_t'.format(mcDir), 5, 't_t-channel {0} 2012'.format(mcTrigger)),
#cfg('tbar_t-channel', '{0}/Tbar_t'.format(mcDir), 5, 'tbar_t-channel {0} 2012'.format(mcTrigger)),
cfg('t_t-channel', '{0}/T_t'.format(mcDir), 5, 't_t-channel {0} 2012'.format(mcTrigger)),
cfg('tbar_t-channel', '{0}/Tbar_t'.format(mcDir), 5, 'tbar_t-channel {0} 2012'.format(mcTrigger)),

cfg('WWW', '{0}/WWWJets'.format(mcDir), 5, 'WWW {0} 2012'.format(mcTrigger)),
cfg('WWZ', '{0}/WWZNoGstarJets'.format(mcDir), 5, 'WWZ {0} 2012'.format(mcTrigger)),
cfg('WZZ', '{0}/WZZNoGstarJets'.format(mcDir), 5, 'WZZ {0} 2012'.format(mcTrigger)),
cfg('ZZZ', '{0}/ZZZNoGstarJets'.format(mcDir), 5, 'ZZZ {0} 2012'.format(mcTrigger)),
cfg('WWG', '{0}/WWGJets'.format(mcDir), 5, 'WWG {0} 2012'.format(mcTrigger)),

#cfg('WWJets2L2Nu', '{0}/WWJetsTo2L2Nu'.format(mcDir), 10, 'WWJets2L2Nu {0} 2012'.format(mcTrigger)),
#cfg('ZZJets2L2Nu', '{0}/ZZJetsTo2L2Nu'.format(mcDir), 10, 'ZZJets2L2Nu {0} 2012'.format(mcTrigger)),
#cfg('WZJets2L2Q', '{0}/WZJetsTo2L2Q'.format(mcDir), 10, 'WZJets2L2Q {0} 2012'.format(mcTrigger)),
#cfg('ZZJets2L2Q', '{0}/ZZJetsTo2L2Q'.format(mcDir), 10, 'ZZJets2L2Q {0} 2012'.format(mcTrigger)),
cfg('WWJets2L2Nu', '{0}/WWJetsTo2L2Nu'.format(mcDir), 10, 'WWJets2L2Nu {0} 2012'.format(mcTrigger)),
cfg('ZZJets2L2Nu', '{0}/ZZJetsTo2L2Nu'.format(mcDir), 10, 'ZZJets2L2Nu {0} 2012'.format(mcTrigger)),
cfg('WZJets2L2Q', '{0}/WZJetsTo2L2Q'.format(mcDir), 10, 'WZJets2L2Q {0} 2012'.format(mcTrigger)),
cfg('ZZJets2L2Q', '{0}/ZZJetsTo2L2Q'.format(mcDir), 10, 'ZZJets2L2Q {0} 2012'.format(mcTrigger)),
cfg('WZJets3LNu', '{0}/WZJetsTo3LNu'.format(mcDir), 10, 'WZJets3LNu {0} 2012'.format(mcTrigger)),
cfg('ZZTo4e', '{0}/ZZTo4e'.format(mcDir), 10, 'ZZ4e {0} 2012'.format(mcTrigger)),
cfg('ZZTo4mu', '{0}/ZZTo4mu'.format(mcDir), 10, 'ZZ4mu {0} 2012'.format(mcTrigger)),
Expand All @@ -105,17 +105,17 @@
cfg('ZZTo2e2tau', '{0}/ZZTo2e2tau'.format(mcDir), 10, 'ZZ2e2tau {0} 2012'.format(mcTrigger)),
cfg('ZZTo2mu2tau', '{0}/ZZTo2mu2tau'.format(mcDir), 10, 'ZZ2mu2tau {0} 2012'.format(mcTrigger)),

#cfg('WmWmqq', '{0}/WmWmqq'.format('/tthome/naodell/storage/ntuples/MC_skimmed'), 10, 'WmWmqq {0} 2012'.format(mcTrigger)),
#cfg('WpWpqq', '{0}/WpWpqq'.format('/tthome/naodell/storage/ntuples/MC_skimmed'), 10, 'WpWpqq {0} 2012'.format(mcTrigger)),
#cfg('WWDPS', '{0}/WW_DPS'.format('/tthome/naodell/storage/ntuples/MC_skimmed'), 10, 'WWDPS {0} 2012'.format(mcTrigger)),
cfg('WmWmqq', '{0}/WmWmqq'.format('/tthome/naodell/storage/ntuples/MC_skimmed'), 10, 'WmWmqq {0} 2012'.format(mcTrigger)),
cfg('WpWpqq', '{0}/WpWpqq'.format('/tthome/naodell/storage/ntuples/MC_skimmed'), 10, 'WpWpqq {0} 2012'.format(mcTrigger)),
cfg('WWDPS', '{0}/WW_DPS'.format('/tthome/naodell/storage/ntuples/MC_skimmed'), 10, 'WWDPS {0} 2012'.format(mcTrigger)),

#cfg('QCD_20-30_EM', '{0}/QCD_Pt_20_30_EMEnriched'.format(mcDir), 40, 'QCD_20-30_EM {0} 2012'.format(mcTrigger)),
#cfg('QCD_30-80_EM', '{0}/QCD_Pt_30_80_EMEnriched'.format(mcDir), 40, 'QCD_30-80_EM {0} 2012'.format(mcTrigger)),
#cfg('QCD_80-170_EM', '{0}/QCD_Pt_80_170_EMEnriched'.format(mcDir), 40, 'QCD_80-170_EM {0} 2012'.format(mcTrigger)),
#cfg('QCD_170-250_EM', '{0}/QCD_Pt_170_250_EMEnriched'.format(mcDir), 40, 'QCD_170-250_EM {0} 2012'.format(mcTrigger)),
#cfg('QCD_250-350_EM', '{0}/QCD_Pt_250_350_EMEnriched'.format(mcDir), 40, 'QCD_250-350_EM {0} 2012'.format(mcTrigger)),
#cfg('QCD_350_EM', '{0}/QCD_Pt_350_EMEnriched'.format(mcDir), 40, 'QCD_350_EM {0} 2012'.format(mcTrigger)),
#cfg('QCD_20_MU', '{0}/QCD_Pt_20_MuEnrichedPt_15'.format(mcDir), 40, 'QCD_20_MU {0} 2012'.format(mcTrigger)),
cfg('QCD_20-30_EM', '{0}/QCD_Pt_20_30_EMEnriched'.format(mcDir), 40, 'QCD_20-30_EM {0} 2012'.format(mcTrigger)),
cfg('QCD_30-80_EM', '{0}/QCD_Pt_30_80_EMEnriched'.format(mcDir), 40, 'QCD_30-80_EM {0} 2012'.format(mcTrigger)),
cfg('QCD_80-170_EM', '{0}/QCD_Pt_80_170_EMEnriched'.format(mcDir), 40, 'QCD_80-170_EM {0} 2012'.format(mcTrigger)),
cfg('QCD_170-250_EM', '{0}/QCD_Pt_170_250_EMEnriched'.format(mcDir), 40, 'QCD_170-250_EM {0} 2012'.format(mcTrigger)),
cfg('QCD_250-350_EM', '{0}/QCD_Pt_250_350_EMEnriched'.format(mcDir), 40, 'QCD_250-350_EM {0} 2012'.format(mcTrigger)),
cfg('QCD_350_EM', '{0}/QCD_Pt_350_EMEnriched'.format(mcDir), 40, 'QCD_350_EM {0} 2012'.format(mcTrigger)),
cfg('QCD_20_MU', '{0}/QCD_Pt_20_MuEnrichedPt_15'.format(mcDir), 40, 'QCD_20_MU {0} 2012'.format(mcTrigger)),

#cfg('ggHToZZ4L_M-125', mcDir+'/GluGluToHToZZTo4L_M-125', 5, 'ggHToZZ4L_M-125 muon 2012'),
#cfg('ggHToWW2L2Nu_M-125', mcDir+'/GluGluToHToWWTo2LAndTau2Nu_M-125', 5, 'ggHToWW2L2Nu_M-125 muon 2012'),
Expand All @@ -128,7 +128,7 @@
#cfg('FCNC_M125_tHj', mcDir+'/TToFCNHToWWTo2l2nuPlusTop_M125', 1, 'FCNC_M125_t mc 2012'),
#cfg('FCNC_M125_tbarHj', mcDir+'/TbarToFCNHToWWTo2l2nuPlusTop_M125/', 1, 'FCNC_M125_tbar mc 2012')
cfg('FCNC_M125_tHj', '{0}/FCNH_semihadronic'.format(signalDir), 1, 'FCNC_M125_t mc 2012'),
cfg('FCNC_M125_tbarHj', '{0}/FCNH_semihadronic_Minus'.format(signalDir), 1, 'FCNC_M125_tbar mc 2012'),
cfg('FCNC_M125_tbarHj', '{0}/FCNH_WW_semileptonic_Minus'.format(signalDir), 1, 'FCNC_M125_tbar mc 2012'),
cfg('FCNC_M125_ZZ_Plus', '{0}/FCNHToZZ_Plus'.format(signalDir), 1, 'FCNC_ZZ_t mc 2012'),
cfg('FCNC_M125_ZZ_Minus', '{0}/FCNHToZZ_Minus'.format(signalDir), 1, 'FCNC_ZZ_tbar mc 2012'),
cfg('FCNC_M125_TauTau_Plus', '{0}/FCNHToTauTau_Plus'.format(signalDir), 1, 'FCNC_TauTau_t mc 2012'),
Expand Down
6 changes: 5 additions & 1 deletion fcncAnalysis/fcncAnalyzer.C
View file
Original file line number Diff line number Diff line change
Expand Up @@ -425,7 +425,7 @@ bool fcncAnalyzer::Process(Long64_t entry)
// Vector bosons
for (unsigned i = 0; i < dubyas.size(); ++i)
if (dubyas[i].GetStatus() == 3 && dubyas[i].Mother() != 0)
cout << "\t status = " << dubyas[i].GetStatus() << ", mass = " << dubyas[i].M() << ", pt = " << dubyas[i].Pt() << ", dubyas" << endl;
cout << "\t status = " << dubyas[i].GetStatus() << ", mass = " << dubyas[i].M() << ", pt = " << dubyas[i].Pt() << ", charge = " << dubyas[i].Charge() << ", dubyas" << endl;

for (unsigned i = 0; i < Zeds.size(); ++i)
if (Zeds[i].GetStatus() == 3 && dubyas[i].Mother() != 0)
Expand Down Expand Up @@ -1991,6 +1991,8 @@ void fcncAnalyzer::MakeQMisIDPlots(vObj& electrons, vector<TCGenParticle>& gElec
"h2_TrailingElecQMisIDNumer", "trailing e charge misID (numerator);p_{T};#eta", nPtBins, elePtBins, nEtaBins, eleEtaBins);
histManager->Fill1DHist(nJets,
"h1_EleQMisIDNumerJets", "N_{jets};N_{jets};Entries", 5, -0.5, 4.5);
histManager->Fill1DHist(primaryVtx->GetSize(),
"h1_EleQMisIDNumerPVs", "N_{PV};N_{PV};Entries", 51, -0.5, 50.5);

histManager->Fill2DHist(nPtBins*iEta1 + iPt1, nPtBins*iEta2 + iPt2,
"h2_DileptonQMisIDNumer", "e charge misID (numerator);e_{leading};e_{trailing}", 9, 0.5, 9.5, 9, 0.5, 9.5);
Expand All @@ -2010,6 +2012,8 @@ void fcncAnalyzer::MakeQMisIDPlots(vObj& electrons, vector<TCGenParticle>& gElec
"h2_TrailingElecQMisIDDenom", "trailing e charge misID (denominator);p_{T};#eta", nPtBins, elePtBins, nEtaBins, eleEtaBins);
histManager->Fill1DHist(nJets,
"h1_EleQMisIDDenomJets", "N_{jets};N_{jets};Entries", 5, -0.5, 4.5);
histManager->Fill1DHist(primaryVtx->GetSize(),
"h1_EleQMisIDDenomPVs", "N_{PV};N_{PV};Entries", 51, -0.5, 50.5);

histManager->Fill2DHist(nPtBins*iEta1 + iPt1, nPtBins*iEta2 + iPt2,
"h2_DileptonQMisIDDenom", "e charge misID (denominator);e_{leading};e_{trailing}", 9, 0.5, 9.5, 9, 0.5, 9.5);
Expand Down
Loading

0 comments on commit 7c7c73a

Please sign in to comment.