[Rivet-svn] r3001 - in trunk: data/plotinfo src/Analyses

[blackhole at projects.hepforge.org]

Author: fsiegert
Date: Fri Mar  4 09:06:20 2011
New Revision: 3001

Log:
Add charge asymmetry vs. pT(W) boson to MC_WJETS analysis.

Modified:
   trunk/data/plotinfo/MC_WJETS.plot
   trunk/src/Analyses/MC_WJETS.cc

Modified: trunk/data/plotinfo/MC_WJETS.plot
==============================================================================
--- trunk/data/plotinfo/MC_WJETS.plot	Wed Mar  2 08:07:07 2011	(r3000)
+++ trunk/data/plotinfo/MC_WJETS.plot	Fri Mar  4 09:06:20 2011	(r3001)
@@ -1,3 +1,22 @@
+# BEGIN PLOT /MC_WJETS/W_chargeasymm_eta
+Title=Charge asymmetry vs. $\eta(\ell)$
+XLabel=$\eta(\ell)$ [GeV]
+YLabel=$(\mathrm{d}\sigma(\ell^+)-\mathrm{d}\sigma(\ell^-))\;/\;(\mathrm{d}\sigma(\ell^+)+\mathrm{d}\sigma(\ell^-))$
+LogY=0
+Scale=1.0
+# END PLOT
+
+# BEGIN PLOT /MC_WJETS/W_chargeasymm_pT
+Title=Charge asymmetry vs. $p_\perp(W)$
+XLabel=$p_\perp(W)$ [GeV]
+YLabel=$\mathrm{d}\sigma(W^+)\;/\;\mathrm{d}\sigma(W^-)$
+LogX=1
+LogY=0
+Scale=1.0
+YMin=0.0
+YMax=3.0
+# END PLOT
+
 # BEGIN PLOT /MC_WJETS/W_jet1_dR
 Title=Separation between W boson and leading jet
 XLabel=$\Delta R(\text{W, 1st jet})$
@@ -19,7 +38,7 @@
 # END PLOT
 
 # BEGIN PLOT /MC_WJETS/W_pT$
-Title=W $p_\perp$
+Title=W boson $p_\perp$
 XLabel=$p_\perp^{\text{W}}$ [GeV]
 YLabel=$\text{d}\sigma/\text{d}p_\perp^{\text{W}}$ [pb/GeV]
 LogX=1
@@ -28,7 +47,7 @@
 # END PLOT
 
 # BEGIN PLOT /MC_WJETS/W_pT_peak
-Title=W $p_\perp$ in peak region
+Title=W boson $p_\perp$ in peak region
 XLabel=$p_\perp^{\text{W}}$ [GeV]
 YLabel=$\text{d}\sigma/\text{d}p_\perp^{\text{W}}$ [pb/GeV]
 LogY=0
@@ -38,7 +57,7 @@
 # END PLOT
 
 # BEGIN PLOT /MC_WJETS/W_y
-Title=W rapidity
+Title=W boson rapidity
 XLabel=$y_{\text{W}}$
 YLabel=$\text{d}\sigma/\text{d}y_{\text{W}}$ [pb]
 LegendXPos=0.35
@@ -51,6 +70,24 @@
 YLabel=$\text{d}\sigma/\text{d}\phi_{\text{W}}$ [pb]
 # END PLOT
 
+# BEGIN PLOT /MC_WJETS/Wplus_pT$
+Title=$p_\perp$ of $W^+$ boson
+XLabel=$p_\perp^{W^+}$ [GeV]
+YLabel=$\text{d}\sigma/\text{d}p_\perp^{W^+}$ [pb/GeV]
+LogX=1
+LegendXPos=0.05
+LegendYPos=0.5
+# END PLOT
+
+# BEGIN PLOT /MC_WJETS/Wminus_pT$
+Title=$p_\perp$ of $W^-$ boson
+XLabel=$p_\perp^{W^-}$ [GeV]
+YLabel=$\text{d}\sigma/\text{d}p_\perp^{W^-}$ [pb/GeV]
+LogX=1
+LegendXPos=0.05
+LegendYPos=0.5
+# END PLOT
+
 # BEGIN PLOT /MC_WJETS/lepton_pT
 Title=Lepton $p_\perp$
 XLabel=$p_\perp^{\text{l}}$ [GeV]

Modified: trunk/src/Analyses/MC_WJETS.cc
==============================================================================
--- trunk/src/Analyses/MC_WJETS.cc	Wed Mar  2 08:07:07 2011	(r3000)
+++ trunk/src/Analyses/MC_WJETS.cc	Fri Mar  4 09:06:20 2011	(r3001)
@@ -38,6 +38,8 @@
       _h_W_phi = bookHistogram1D("W_phi", 25, 0.0, TWOPI);
       _h_W_jet1_deta = bookHistogram1D("W_jet1_deta", 50, -5.0, 5.0);
       _h_W_jet1_dR = bookHistogram1D("W_jet1_dR", 25, 0.5, 7.0);
+      _h_Wplus_pT = bookHistogram1D("Wplus_pT", logBinEdges(100, 1.0, 0.5*sqrtS()));
+      _h_Wminus_pT = bookHistogram1D("Wminus_pT", logBinEdges(100, 1.0, 0.5*sqrtS()));
       _h_lepton_pT = bookHistogram1D("lepton_pT", logBinEdges(100, 10.0, 0.25*sqrtS()));
       _h_lepton_eta = bookHistogram1D("lepton_eta", 40, -4.0, 4.0);
       _htmp_dsigminus_deta = bookHistogram1D("lepton_dsigminus_deta", 20, 0.0, 4.0);
@@ -80,6 +82,11 @@
           _htmp_dsigplus_deta->fill(emom.eta(), weight);
         }
       }
+      if (charge3_x_eta < 0) {
+        _h_Wminus_pT->fill(wmom.pT(), weight);
+      } else {
+        _h_Wplus_pT->fill(wmom.pT(), weight);
+      }
 
       const FastJets& jetpro = applyProjection<FastJets>(e, "Jets");
       const Jets& jets = jetpro.jetsByPt(20.0*GeV);
@@ -115,6 +122,11 @@
       hf.destroy(_htmp_dsigminus_deta);
       hf.destroy(_htmp_dsigplus_deta);
 
+      // W charge asymmetry vs. pTW: dsig+/dpT / dsig-/dpT
+      hf.divide(histoDir() + "/W_chargeasymm_pT", *_h_Wplus_pT, *_h_Wminus_pT);
+      scale(_h_Wplus_pT, crossSection()/sumOfWeights());
+      scale(_h_Wminus_pT, crossSection()/sumOfWeights());
+
       MC_JetAnalysis::finalize();
     }
 
@@ -132,6 +144,8 @@
     AIDA::IHistogram1D * _h_W_phi;
     AIDA::IHistogram1D * _h_W_jet1_deta;
     AIDA::IHistogram1D * _h_W_jet1_dR;
+    AIDA::IHistogram1D * _h_Wplus_pT;
+    AIDA::IHistogram1D * _h_Wminus_pT;
     AIDA::IHistogram1D * _h_lepton_pT;
     AIDA::IHistogram1D * _h_lepton_eta;