|
|
[Rivet-svn] r2115 - in trunk: data/plotinfo src/Analyses
blackhole at projects.hepforge.org
blackhole at projects.hepforge.org
Mon Nov 30 20:29:30 GMT 2009
Author: buckley
Date: Mon Nov 30 20:29:30 2009
New Revision: 2115
Log:
W analysis plot name changes
Modified:
trunk/data/plotinfo/MC_LHC_WANALYSIS.plot
trunk/src/Analyses/MC_LHC_WANALYSIS.cc
Modified: trunk/data/plotinfo/MC_LHC_WANALYSIS.plot
==============================================================================
--- trunk/data/plotinfo/MC_LHC_WANALYSIS.plot Mon Nov 30 19:40:49 2009 (r2114)
+++ trunk/data/plotinfo/MC_LHC_WANALYSIS.plot Mon Nov 30 20:29:30 2009 (r2115)
@@ -1,144 +1,144 @@
-# BEGIN PLOT /MC_LHC_WANALYSIS/n-ch
+# BEGIN PLOT /MC_LHC_WANALYSIS/track-n
Title=Charged multiplicity
XLabel=$n_\text{ch}$
YLabel=$\mathrm{d}N_\text{evt}/\mathrm{d}n_\text{ch}$
LogY=0
# END PLOT
-# BEGIN PLOT /MC_LHC_WANALYSIS/n-w
+# BEGIN PLOT /MC_LHC_WANALYSIS/w-n
Title=W candidate multiplicity
XLabel=$n_\text{W}^\pm$
YLabel=$\mathrm{d}N_\text{evt}/\mathrm{d}n_\text{W}$
LogY=0
# END PLOT
-# BEGIN PLOT /MC_LHC_WANALYSIS/n-wplus
+# BEGIN PLOT /MC_LHC_WANALYSIS/wplus-n
Title=\PWplus candidate multiplicity
XLabel=$n_\text{W}^+$
YLabel=$\mathrm{d}N_\text{evt}/\mathrm{d}n_\text{W}$
LogY=0
# END PLOT
-# BEGIN PLOT /MC_LHC_WANALYSIS/n-wminus
+# BEGIN PLOT /MC_LHC_WANALYSIS/wminus-n
Title=\PWminus candidate multiplicity
XLabel=$n_\text{W}^-$
YLabel=$\mathrm{d}N_\text{evt}/\mathrm{d}n_\text{W}$
LogY=0
# END PLOT
-# BEGIN PLOT /MC_LHC_WANALYSIS/n-jet
+# BEGIN PLOT /MC_LHC_WANALYSIS/jet-n
Title=Jet multiplicity
XLabel=$n_\text{jet}$
YLabel=$\mathrm{d}N_\text{evt}/\mathrm{d}n_\text{jet}$
LogY=0
# END PLOT
-# BEGIN PLOT /MC_LHC_WANALYSIS/pt-ch
+# BEGIN PLOT /MC_LHC_WANALYSIS/track-pt
Title=Charged track $p_\perp$
XLabel=$p_\perp$ / GeV
YLabel=$\mathrm{d}N_\text{ch}/\mathrm{d}p_\perp^\text{ch}$ / GeV$^{-1}$
LogY=1
# END PLOT
-# BEGIN PLOT /MC_LHC_WANALYSIS/ptavg-ch
+# BEGIN PLOT /MC_LHC_WANALYSIS/track-ptavg
Title=Average charged track $p_\perp$
XLabel=$\langle p_\perp^\text{ch} \rangle$ / GeV
YLabel=$\mathrm{d}N_\text{ch}/\mathrm{d}\langle p_\perp^\text{ch} \langle$ / GeV$^{-1}$
LogY=1
# END PLOT
-# BEGIN PLOT /MC_LHC_WANALYSIS/ptrms-ch
+# BEGIN PLOT /MC_LHC_WANALYSIS/track-ptrms
Title=Width of charged track $p_\perp$
XLabel=$\hat{\sigma}_{p_\perp^\text{ch}}$ / GeV
YLabel=$\mathrm{d}N_\text{ch}/\mathrm{d}\hat{\sigma}_{p_\perp^\text{ch}}$ / GeV$^{-1}$
LogY=1
# END PLOT
-# BEGIN PLOT /MC_LHC_WANALYSIS/pt-w
+# BEGIN PLOT /MC_LHC_WANALYSIS/w-pt
Title=W candidate $p_\perp$
XLabel=$p_\perp(\text{W})$ / GeV
YLabel=$\mathrm{d}N_\text{W}/\mathrm{d}p_\perp(\text{W})$ / GeV$^{-1}$
LogY=1
# END PLOT
-# BEGIN PLOT /MC_LHC_WANALYSIS/pt-wplus
+# BEGIN PLOT /MC_LHC_WANALYSIS/wplus-pt
Title=\PWplus candidate $p_\perp$
XLabel=$p_\perp(\text{\PWplus})$ / GeV
YLabel=$\mathrm{d}N_\text{W}/\mathrm{d}p_\perp(\text{W})$ / GeV$^{-1}$
LogY=1
# END PLOT
-# BEGIN PLOT /MC_LHC_WANALYSIS/pt-wminus
+# BEGIN PLOT /MC_LHC_WANALYSIS/wminus-pt
Title=\PWminus candidate $p_\perp$
XLabel=$p_\perp(\text{\PWminus})$ / GeV
YLabel=$\mathrm{d}N_\text{W}/\mathrm{d}p_\perp(\text{W})$ / GeV$^{-1}$
LogY=1
# END PLOT
-# BEGIN PLOT /MC_LHC_WANALYSIS/pt-jet
+# BEGIN PLOT /MC_LHC_WANALYSIS/jet-pt
Title=Jet $p_\perp$
XLabel=$p_\perp(\text{jet})$ / GeV
YLabel=$\mathrm{d}N_\text{jet}/\mathrm{d}p_\perp(\text{jet})$ / GeV$^{-1}$
LogY=1
# END PLOT
-# BEGIN PLOT /MC_LHC_WANALYSIS/eta-w
+# BEGIN PLOT /MC_LHC_WANALYSIS/w-eta
Title=W candidate $\eta$
XLabel=$\eta(\text{W})$
YLabel=$\mathrm{d}N_\text{W}/\mathrm{d}\eta(\text{W})$
LogY=0
# END PLOT
-# BEGIN PLOT /MC_LHC_WANALYSIS/eta-wplus
+# BEGIN PLOT /MC_LHC_WANALYSIS/wplus-eta
Title=\PWplus candidate $\eta$
XLabel=$\eta(\text{W})$
YLabel=$\mathrm{d}N_\text{W}/\mathrm{d}\eta(\text{W})$
LogY=0
# END PLOT
-# BEGIN PLOT /MC_LHC_WANALYSIS/eta-wminus
+# BEGIN PLOT /MC_LHC_WANALYSIS/wminus-eta
Title=\PWminus candidate $\eta$
XLabel=$\eta(\text{W})$
YLabel=$\mathrm{d}N_\text{W}/\mathrm{d}\eta(\text{W})$
LogY=0
# END PLOT
-# BEGIN PLOT /MC_LHC_WANALYSIS/phi-w
+# BEGIN PLOT /MC_LHC_WANALYSIS/w-phi
Title=W candidate $\phi$
XLabel=$\phi(\text{W})$
YLabel=$\mathrm{d}N_\text{W}/\mathrm{d}\phi(\text{W})$
LogY=0
# END PLOT
-# BEGIN PLOT /MC_LHC_WANALYSIS/phi-wplus
+# BEGIN PLOT /MC_LHC_WANALYSIS/wplus-phi
Title=\PWplus candidate $\phi$
XLabel=$\phi(\text{W})$
YLabel=$\mathrm{d}N_\text{W}/\mathrm{d}\phi(\text{W})$
LogY=0
# END PLOT
-# BEGIN PLOT /MC_LHC_WANALYSIS/phi-wminus
+# BEGIN PLOT /MC_LHC_WANALYSIS/wminus-phi
Title=\PWminus candidate $\phi$
XLabel=$\phi(\text{W})$
YLabel=$\mathrm{d}N_\text{W}/\mathrm{d}\phi(\text{W})$
LogY=0
# END PLOT
-# BEGIN PLOT /MC_LHC_WANALYSIS/m-w
+# BEGIN PLOT /MC_LHC_WANALYSIS/w-m
Title=W candidate mass
XLabel=$m(\text{W})$ / GeV
YLabel=$\mathrm{d}N_\text{W}/\mathrm{d}m(\text{W})$ / GeV$^{-1}$
LogY=1
# END PLOT
-# BEGIN PLOT /MC_LHC_WANALYSIS/m-wplus
+# BEGIN PLOT /MC_LHC_WANALYSIS/wplus-m
Title=\PWplus candidate mass
XLabel=$m(\text{W})$ / GeV
YLabel=$\mathrm{d}N_\text{W}/\mathrm{d}m(\text{W})$ / GeV$^{-1}$
LogY=1
# END PLOT
-# BEGIN PLOT /MC_LHC_WANALYSIS/m-wminus
+# BEGIN PLOT /MC_LHC_WANALYSIS/wminus-m
Title=\PWminus candidate mass
XLabel=$m(\text{W})$ / GeV
YLabel=$\mathrm{d}N_\text{W}/\mathrm{d}m(\text{W})$ / GeV$^{-1}$
Modified: trunk/src/Analyses/MC_LHC_WANALYSIS.cc
==============================================================================
--- trunk/src/Analyses/MC_LHC_WANALYSIS.cc Mon Nov 30 19:40:49 2009 (r2114)
+++ trunk/src/Analyses/MC_LHC_WANALYSIS.cc Mon Nov 30 20:29:30 2009 (r2115)
@@ -32,42 +32,37 @@
FastJets fastjets(wf.remainingFinalState(), FastJets::KT, 0.7);
addProjection(fastjets, "Jets");
- _hist_chargemulti = bookHistogram1D("n-ch", 50, -0.5, 399.5);
- _hist_chargept = bookHistogram1D("pt-ch", 25, 0, 25);
+ _hist_chargemulti = bookHistogram1D("track-n", 50, -0.5, 399.5);
+ _hist_chargept = bookHistogram1D("track-pt", 25, 0, 25);
/// @todo Use profile plots instead:
- _hist_chargemeanpt = bookHistogram1D("ptavg-ch", 20, 0, 10);
+ _hist_chargemeanpt = bookHistogram1D("track-ptavg", 20, 0, 10);
/// @todo Use profile plots instead (and this isn't really an "RMS")
- _hist_chargermspt = bookHistogram1D("ptrms-ch", 20, 0, 10);
+ _hist_chargermspt = bookHistogram1D("track-ptrms", 20, 0, 10);
//
- _hist_wcount = bookHistogram1D("n-w", 6, -0.5, 5.5);
- _hist_wpluscount = bookHistogram1D("n-wplus", 6, -0.5, 5.5);
- _hist_wminuscount = bookHistogram1D("n-wminus", 6, -0.5, 5.5);
- _hist_wpt = bookHistogram1D("pt-w", 25, 0, 100);
- _hist_wpluspt = bookHistogram1D("pt-wplus", 25, 0, 100);
- _hist_wminuspt = bookHistogram1D("pt-wminus", 25, 0, 100);
- _hist_weta = bookHistogram1D("eta-w", 36, -6, 6);
- _hist_wpluseta = bookHistogram1D("eta-wplus", 36, -6, 6);
- _hist_wminuseta = bookHistogram1D("eta-wminus", 36, -6, 6);
- _hist_wphi = bookHistogram1D("phi-w", 25, 0, TWOPI);
- _hist_wplusphi = bookHistogram1D("phi-wplus", 25, 0, TWOPI);
- _hist_wminusphi = bookHistogram1D("phi-wminus", 25, 0, TWOPI);
- _hist_wmass = bookHistogram1D("m-w", 40, 60, 100);
- _hist_wplusmass = bookHistogram1D("m-wplus", 40, 60, 100);
- _hist_wminusmass = bookHistogram1D("m-wminus", 40, 60, 100);
+ _hist_wcount = bookHistogram1D("w-n", 6, -0.5, 5.5);
+ _hist_wpluscount = bookHistogram1D("wplus-n", 6, -0.5, 5.5);
+ _hist_wminuscount = bookHistogram1D("wminus-n", 6, -0.5, 5.5);
+ _hist_wpt = bookHistogram1D("w-pt", 25, 0, 100);
+ _hist_wpluspt = bookHistogram1D("wplus-pt", 25, 0, 100);
+ _hist_wminuspt = bookHistogram1D("wminus-pt", 25, 0, 100);
+ _hist_weta = bookHistogram1D("w-eta", 36, -6, 6);
+ _hist_wpluseta = bookHistogram1D("wplus-eta", 36, -6, 6);
+ _hist_wminuseta = bookHistogram1D("wminus-eta", 36, -6, 6);
+ _hist_wphi = bookHistogram1D("w-phi", 25, 0, TWOPI);
+ _hist_wplusphi = bookHistogram1D("wplus-phi", 25, 0, TWOPI);
+ _hist_wminusphi = bookHistogram1D("wminus-phi", 25, 0, TWOPI);
+ _hist_wmass = bookHistogram1D("w-m", 40, 60, 100);
+ _hist_wplusmass = bookHistogram1D("wplus-m", 40, 60, 100);
+ _hist_wminusmass = bookHistogram1D("wminus-m", 40, 60, 100);
//_hist_weta_asymm = bookProfile1D("asymm-eta-w", 20, -5.0, 5.0);
//
- _hist_jetcount = bookHistogram1D("n-jet", 11, -0.5, 10.5);
- _hist_jetpt = bookHistogram1D("pt-jet", 50, 20, 100);
+ _hist_jetcount = bookHistogram1D("jet-n", 6, -0.5, 5.5);
+ _hist_jetpt = bookHistogram1D("jet-pt", 50, 20, 100);
}
void analyze(const Event& event) {
- const FinalState& cfs = applyProjection<FinalState>(event, "CFS");
const WFinder& wf = applyProjection<WFinder>(event, "WF");
- const FastJets& fastjets = applyProjection<FastJets>(event, "Jets");
- const Jets jets = fastjets.jetsByPt();
-
- // Veto if no Ws found
if (wf.size() == 0) {
getLog() << Log::DEBUG << "No W candidates found: vetoing" << endl;
vetoEvent;
@@ -76,6 +71,7 @@
_sumWPass += weight;
// Charged particles part
+ const FinalState& cfs = applyProjection<FinalState>(event, "CFS");
_hist_chargemulti->fill(cfs.particles().size(), weight);
double meanpt(0), rmspt(0);
foreach (const Particle& p, cfs.particles()) {
@@ -123,8 +119,10 @@
_hist_wminuscount->fill(n_wminus, weight);
// Jet part
- _hist_jetcount->fill(fastjets.jets().size(), weight);
- foreach(const Jet& j, fastjets.jetsByPt()) {
+ const FastJets& fastjets = applyProjection<FastJets>(event, "Jets");
+ const Jets jets = fastjets.jetsByPt();
+ _hist_jetcount->fill(jets.size(), weight);
+ foreach (const Jet& j, jets) {
const double pT = j.momentum().pT();
_hist_jetpt->fill(pT/GeV, weight);
}
More information about the Rivet-svn
mailing list
|