|
[Rivet-svn] r3238 - in trunk: data/plotinfo src/Analysesblackhole at projects.hepforge.org blackhole at projects.hepforge.orgFri Jul 22 10:35:03 BST 2011
Author: fsiegert Date: Fri Jul 22 10:35:03 2011 New Revision: 3238 Log: Add linear photon pT histo to MC_PHOTONJETS. Modified: trunk/data/plotinfo/MC_PHOTONJETS.plot trunk/src/Analyses/MC_PHOTONJETS.cc Modified: trunk/data/plotinfo/MC_PHOTONJETS.plot ============================================================================== --- trunk/data/plotinfo/MC_PHOTONJETS.plot Fri Jul 22 10:26:56 2011 (r3237) +++ trunk/data/plotinfo/MC_PHOTONJETS.plot Fri Jul 22 10:35:03 2011 (r3238) @@ -16,13 +16,19 @@ YLabel=$\text{d}\sigma/\text{d}\Delta{\phi}(\gamma\text{, 1st jet})$ [pb] # END PLOT -# BEGIN PLOT /MC_PHOTONJETS/photon_pT +# BEGIN PLOT /MC_PHOTONJETS/photon_pT$ Title=Photon $p_\perp$ XLabel=$p_\perp^{\gamma}$ [GeV] YLabel=$\text{d}\sigma/\text{d}p_\perp^{\gamma}$ [pb/GeV] LogX=1 # END PLOT +# BEGIN PLOT /MC_PHOTONJETS/photon_pT_lin$ +Title=Photon $p_\perp$ +XLabel=$p_\perp^{\gamma}$ [GeV] +YLabel=$\text{d}\sigma/\text{d}y_{\gamma}$ [pb/GeV] +# END PLOT + # BEGIN PLOT /MC_PHOTONJETS/photon_y Title=Photon rapidity XLabel=$y_{\gamma}$ Modified: trunk/src/Analyses/MC_PHOTONJETS.cc ============================================================================== --- trunk/src/Analyses/MC_PHOTONJETS.cc Fri Jul 22 10:26:56 2011 (r3237) +++ trunk/src/Analyses/MC_PHOTONJETS.cc Fri Jul 22 10:35:03 2011 (r3238) @@ -42,6 +42,7 @@ addProjection(jetpro, "Jets"); _h_photon_pT = bookHistogram1D("photon_pT", logBinEdges(50, 30.0, 0.5*sqrtS())); + _h_photon_pT_lin = bookHistogram1D("photon_pT_lin", 50, 0.0, 70.0); _h_photon_y = bookHistogram1D("photon_y", 50, -5.0, 5.0); _h_photon_jet1_deta = bookHistogram1D("photon_jet1_deta", 50, -5.0, 5.0); _h_photon_jet1_dphi = bookHistogram1D("photon_jet1_dphi", 20, 0.0, M_PI); @@ -83,6 +84,7 @@ } _h_photon_pT->fill(photon.pT(),weight); + _h_photon_pT_lin->fill(photon.pT(),weight); _h_photon_y->fill(photon.rapidity(),weight); const FastJets& jetpro = applyProjection<FastJets>(e, "Jets"); @@ -100,6 +102,7 @@ // Finalize void finalize() { scale(_h_photon_pT, crossSectionPerEvent()); + scale(_h_photon_pT_lin, crossSectionPerEvent()); scale(_h_photon_y, crossSectionPerEvent()); scale(_h_photon_jet1_deta, crossSectionPerEvent()); scale(_h_photon_jet1_dphi, crossSectionPerEvent()); @@ -116,6 +119,7 @@ /// @name Histograms //@{ AIDA::IHistogram1D * _h_photon_pT; + AIDA::IHistogram1D * _h_photon_pT_lin; AIDA::IHistogram1D * _h_photon_y; AIDA::IHistogram1D * _h_photon_jet1_deta; AIDA::IHistogram1D * _h_photon_jet1_dphi;
More information about the Rivet-svn mailing list |