|
[Rivet-svn] r4062 - in branches/2012-06-aidarivet: . data/anainfo data/plotinfo data/refdata src/Analysesblackhole at projects.hepforge.org blackhole at projects.hepforge.orgThu Dec 6 21:24:58 GMT 2012
Author: hoeth Date: Thu Dec 6 21:24:58 2012 New Revision: 4062 Log: Added CMS_2012_I1184941 analysis Added: branches/2012-06-aidarivet/data/anainfo/CMS_2012_I1184941.info branches/2012-06-aidarivet/data/plotinfo/CMS_2012_I1184941.plot branches/2012-06-aidarivet/data/refdata/CMS_2012_I1184941.aida branches/2012-06-aidarivet/src/Analyses/CMS_2012_I1184941.cc Modified: branches/2012-06-aidarivet/ChangeLog branches/2012-06-aidarivet/data/anainfo/Makefile.am branches/2012-06-aidarivet/data/plotinfo/Makefile.am branches/2012-06-aidarivet/data/refdata/Makefile.am branches/2012-06-aidarivet/src/Analyses/Makefile.am Modified: branches/2012-06-aidarivet/ChangeLog ============================================================================== --- branches/2012-06-aidarivet/ChangeLog Thu Dec 6 15:53:49 2012 (r4061) +++ branches/2012-06-aidarivet/ChangeLog Thu Dec 6 21:24:58 2012 (r4062) @@ -2,6 +2,7 @@ * Added ATLAS_2012_I1188891 analysis * Added ATLAS_2012_I1118269 analysis + * Added CMS_2012_I1184941 analysis * Added LHCB_2010_I867355 analysis * Added TGraphErrors support to root2flat Added: branches/2012-06-aidarivet/data/anainfo/CMS_2012_I1184941.info ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ branches/2012-06-aidarivet/data/anainfo/CMS_2012_I1184941.info Thu Dec 6 21:24:58 2012 (r4062) @@ -0,0 +1,40 @@ +Name: CMS_2012_I1184941 +Year: 2012 +Summary: Measurement of the differential cross section for inclusive dijet production as a function of $\xi$ in 7 TeV proton-proton collisions. +Experiment: CMS +Collider: LHC +SpireID: 9762035 +InspireID: 1184941 +Status: VALIDATED +Authors: + - Sercan Sen <ssen at cern.ch> + - Alexander Proskuryakov <aproskur at mail.cern.ch> +References: + - arXiv:1209.1805 [hep-ex] + - Submitted to Phys. Rev. D +RunInfo: + High statistics is needed to observe events in the lowest (xi) bin. + Distributions are presented for HardQCD events (i.e., in PYTHIA with ptHatMin=15GeV) and Diffractive events. +NumEvents: 100M +NeedCrossSection: no +Beams: [p+, p+] +Energies: [7000] +PtCuts: [15] +Description: + Measurement of the differential cross section for inclusive dijet production as a function of $\xi$ + which approximates the fractional momentum loss of the scattered proton in single-diffraction events. The + data used has a total integrated luminosity of 2.7 nb-1 collected during 2010 with low instantaneous luminosity. + Events are selected with at least two jets in $|\eta| < 4.4$ with $pT > 20$ GeV and all final states particles + are used for the reconstruction of $\xi$. +BibKey: :2012vc +BibTeX: '@article{:2012vc, + author = "Chatrchyan, Serguei and others", + title = "{Observation of a diffractive contribution to dijet + production in proton-proton collisions at sqrt(s)=7 TeV}", + collaboration = "CMS Collaboration", + year = "2012", + eprint = "1209.1805", + archivePrefix = "arXiv", + primaryClass = "hep-ex", + reportNumber = "CMS-FWD-10-004, CERN-PH-EP-2012-248", + SLACcitation = "%%CITATION = ARXIV:1209.1805;%%",}' Modified: branches/2012-06-aidarivet/data/anainfo/Makefile.am ============================================================================== --- branches/2012-06-aidarivet/data/anainfo/Makefile.am Thu Dec 6 15:53:49 2012 (r4061) +++ branches/2012-06-aidarivet/data/anainfo/Makefile.am Thu Dec 6 21:24:58 2012 (r4062) @@ -129,6 +129,7 @@ CMS_2012_I1087342.info \ CMS_2012_I1102908.info \ CMS_2012_I1107658.info \ + CMS_2012_I1184941.info \ CMS_QCD_10_024.info \ D0_1996_S3214044.info \ D0_1996_S3324664.info \ Added: branches/2012-06-aidarivet/data/plotinfo/CMS_2012_I1184941.plot ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ branches/2012-06-aidarivet/data/plotinfo/CMS_2012_I1184941.plot Thu Dec 6 21:24:58 2012 (r4062) @@ -0,0 +1,9 @@ +# BEGIN PLOT /CMS_2012_I1184941/d01-x01-y01 +Title=$\sqrt{s}=7$ TeV, $pp \to \text{jet}_{1}\text{jet}_{2}$, $|\eta^{j_1,j_2}|<4.4$, $p_{T}^{j_1,j_2} > 20$ GeV +XLabel=$\xi$ +YLabel=$\mathrm{d}{\sigma_{jj}}/\mathrm{d}{\xi}$ [$\mu$b] +FullRange=1 +LogX=1 +LogY=1 +LegendXPos=0.08 +# END PLOT Modified: branches/2012-06-aidarivet/data/plotinfo/Makefile.am ============================================================================== --- branches/2012-06-aidarivet/data/plotinfo/Makefile.am Thu Dec 6 15:53:49 2012 (r4061) +++ branches/2012-06-aidarivet/data/plotinfo/Makefile.am Thu Dec 6 21:24:58 2012 (r4062) @@ -126,6 +126,7 @@ CMS_2012_I1087342.plot \ CMS_2012_I1102908.plot \ CMS_2012_I1107658.plot \ + CMS_2012_I1184941.plot \ CMS_QCD_10_024.plot \ D0_1996_S3214044.plot \ D0_1996_S3324664.plot \ Added: branches/2012-06-aidarivet/data/refdata/CMS_2012_I1184941.aida ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ branches/2012-06-aidarivet/data/refdata/CMS_2012_I1184941.aida Thu Dec 6 21:24:58 2012 (r4062) @@ -0,0 +1,21 @@ +<?xml version="1.0" encoding="ISO-8859-1" ?> +<!DOCTYPE aida SYSTEM "http://aida.freehep.org/schemas/3.3/aida.dtd"> +<aida version="3.3"> + <implementation version="1.1" package="LWH"/> + <dataPointSet name="d01-x01-y01" title="" path="/REF/CMS_2012_I1184941" dimension="2"> + <dimension dim="0" title="" /> + <dimension dim="1" title="" /> + <dataPoint> + <measurement value="1.150000e-03" errorPlus="8.500000e-04" errorMinus="8.500000e-04"/> + <measurement value="5.000000e+00" errorPlus="1.749285e+00" errorMinus="1.581138e+00"/> + </dataPoint> + <dataPoint> + <measurement value="3.250000e-03" errorPlus="1.250000e-03" errorMinus="1.250000e-03"/> + <measurement value="8.200000e+00" errorPlus="2.376972e+00" errorMinus="2.563201e+00"/> + </dataPoint> + <dataPoint> + <measurement value="7.250000e-03" errorPlus="2.750000e-03" errorMinus="2.750000e-03"/> + <measurement value="1.350000e+01" errorPlus="4.589117e+00" errorMinus="3.228002e+00"/> + </dataPoint> + </dataPointSet> +</aida> Modified: branches/2012-06-aidarivet/data/refdata/Makefile.am ============================================================================== --- branches/2012-06-aidarivet/data/refdata/Makefile.am Thu Dec 6 15:53:49 2012 (r4061) +++ branches/2012-06-aidarivet/data/refdata/Makefile.am Thu Dec 6 21:24:58 2012 (r4062) @@ -83,6 +83,7 @@ CMS_2012_I1087342.aida \ CMS_2012_I1102908.aida \ CMS_2012_I1107658.aida \ + CMS_2012_I1184941.aida \ CMS_QCD_10_024.aida \ DELPHI_1994_S3021912.aida \ DELPHI_1995_S3137023.aida \ Added: branches/2012-06-aidarivet/src/Analyses/CMS_2012_I1184941.cc ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ branches/2012-06-aidarivet/src/Analyses/CMS_2012_I1184941.cc Thu Dec 6 21:24:58 2012 (r4062) @@ -0,0 +1,74 @@ +// -*- C++ -*- +#include "Rivet/Analysis.hh" +#include "Rivet/RivetAIDA.hh" +#include "Rivet/Projections/FastJets.hh" +#include "Rivet/Projections/FinalState.hh" +#include "Rivet/Tools/Logging.hh" + +namespace Rivet { + + + class CMS_2012_I1184941 : public Analysis { + public: + + CMS_2012_I1184941() + : Analysis("CMS_2012_I1184941") + { } + + public: + + void init() { + FinalState fs; + addProjection(fs, "FS"); + + const FastJets jets(FinalState(-4.9, 4.9, 0.0*GeV), FastJets::ANTIKT, 0.5); + addProjection(jets, "AntiKtJets05"); + + _h_xi = bookHistogram1D(1, 1, 1); + } + + + void analyze(const Event& event) { + const double weight = event.weight(); + + double xiM = 0.; + double xiP = 0.; + + const Jets jets = applyProjection<FastJets>(event, "AntiKtJets05").jetsByPt(20.*GeV, MAXDOUBLE, -4.4, 4.4); + if (jets.size() < 2) vetoEvent; // require a dijet system with a 20 GeV cut on both jets + + const FinalState& fsp = applyProjection<FinalState>(event, "FS"); + + foreach (const Particle& p, fsp.particlesByEta()) { + double eta = p.momentum().eta(); + double energy = p.momentum().E(); + double costheta = cos(p.momentum().theta()); + + // Yes, they really correct to +/- infinity, using Pythia 8 ... + if ( eta < 4.9 ) xiP += (energy + energy*costheta); + if ( -4.9 < eta ) xiM += (energy - energy*costheta); + } + + xiP = xiP / (sqrtS()/GeV); + xiM = xiM / (sqrtS()/GeV); + + _h_xi->fill( xiM, weight ); // Fill the histogram both with xiP and xiM, and get the average in the endjob. + _h_xi->fill( xiP, weight ); + + } + + + void finalize() { + scale( _h_xi, crossSection()/microbarn/sumOfWeights() / 2.); + } + + private: + + AIDA::IHistogram1D* _h_xi; + + }; + + // The hook for the plugin system + DECLARE_RIVET_PLUGIN(CMS_2012_I1184941); + +} Modified: branches/2012-06-aidarivet/src/Analyses/Makefile.am ============================================================================== --- branches/2012-06-aidarivet/src/Analyses/Makefile.am Thu Dec 6 15:53:49 2012 (r4061) +++ branches/2012-06-aidarivet/src/Analyses/Makefile.am Thu Dec 6 21:24:58 2012 (r4062) @@ -127,7 +127,8 @@ CMS_2011_S9215166.cc \ CMS_2012_I1087342.cc \ CMS_2012_I1102908.cc \ - CMS_2012_I1107658.cc + CMS_2012_I1107658.cc \ + CMS_2012_I1184941.cc if ENABLE_PRELIMINARY RivetCMSAnalyses_la_SOURCES += \ CMS_QCD_10_024.cc
More information about the Rivet-svn mailing list |