D0_2010_S8671338.cc

Go to the documentation of this file.
00001 // -*- C++ -*-
00002 #include "Rivet/Analysis.hh"
00003 #include "Rivet/RivetAIDA.hh"
00004 #include "Rivet/Tools/Logging.hh"
00005 #include "Rivet/Projections/FinalState.hh"
00006 #include "Rivet/Projections/ZFinder.hh"
00007 
00008 
00009 namespace Rivet {
00010 
00011   /// @brief Measurement of Z(->muon muon) pT differential cross-section
00012   /// @author Flavia Dias
00013   class D0_2010_S8671338 : public Analysis { 
00014 
00015   public:
00016 
00017     /// @name Construction
00018     //@{
00019     /// Constructor
00020     D0_2010_S8671338() : Analysis("D0_2010_S8671338")
00021     {
00022       setBeams(PROTON, ANTIPROTON);
00023     }
00024     
00025     //@}
00026 
00027     ///@name Analysis methods
00028     //@{
00029 
00030     ///Add Projections, Book histograms
00031     void init() {
00032       ZFinder zfinder(-1.7, 1.7, 15.0*GeV, MUON, 65.0*GeV, 115.0*GeV, 0.2);
00033       addProjection(zfinder, "ZFinder");
00034   
00035       _h_Z_pT_cross_section = bookHistogram1D(1, 1, 1);
00036     }
00037 
00038     // Do the analysis
00039     void analyze(const Event& e) {
00040       const double weight = e.weight();
00041       
00042       const ZFinder& zfinder = applyProjection<ZFinder>(e, "ZFinder");
00043       if (zfinder.particles().size()==1) {
00044      
00045         const FourMomentum Zmom = zfinder.particles()[0].momentum();
00046     // In Z pT
00047         _h_Z_pT_cross_section->fill(Zmom.pT(), weight);      
00048 
00049       }
00050 
00051     }
00052 
00053 
00054  /// Finalize
00055     void finalize() {
00056       normalize(_h_Z_pT_cross_section, 1.0);
00057     }
00058     //@}
00059 
00060 
00061   private:
00062 
00063     /// @name Histograms
00064     //@{
00065     AIDA::IHistogram1D * _h_Z_pT_cross_section;
00066     //@}
00067 
00068 
00069 
00070   };
00071 
00072   // This global object acts as a hook for the plugin system
00073   AnalysisBuilder<D0_2010_S8671338> plugin_D0_2010_S8671338;
00074 
00075 
00076 }