|
[Rivet-svn] r4090 - in branches/2012-06-aidarivet: . data/anainfo data/plotinfo data/refdata src/Analysesblackhole at projects.hepforge.org blackhole at projects.hepforge.orgMon Dec 17 17:45:07 GMT 2012
Author: hoeth Date: Mon Dec 17 17:45:07 2012 New Revision: 4090 Log: Added CMS_2012_I1193338 analysis Added: branches/2012-06-aidarivet/data/anainfo/CMS_2012_I1193338.info branches/2012-06-aidarivet/data/plotinfo/CMS_2012_I1193338.plot branches/2012-06-aidarivet/data/refdata/CMS_2012_I1193338.aida branches/2012-06-aidarivet/src/Analyses/CMS_2012_I1193338.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 Mon Dec 17 17:34:17 2012 (r4089) +++ branches/2012-06-aidarivet/ChangeLog Mon Dec 17 17:45:07 2012 (r4090) @@ -1,3 +1,9 @@ +2012-12-17 Hendrik Hoeth <hendrik.hoeth at cern.ch> + + * Added CMS_2012_I1193338 analysis + + * Fixed xi cut in ATLAS_2011_I894867 + 2012-12-17 Andy Buckley <andy.buckley at cern.ch> * Adding analysis descriptions to the HTML analysis page ToC. Added: branches/2012-06-aidarivet/data/anainfo/CMS_2012_I1193338.info ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ branches/2012-06-aidarivet/data/anainfo/CMS_2012_I1193338.info Mon Dec 17 17:45:07 2012 (r4090) @@ -0,0 +1,38 @@ +Name: CMS_2012_I1193338 +Year: 2012 +Summary: Measurement of the inelastic proton-proton cross section at sqrt(s)=7 TeV +Experiment: CMS +Collider: LHC +SpireID: +InspireID: 1193338 +Status: VALIDATED +Authors: + - Sercan Sen <ssen at cern.ch> +References: + - arXiv:1210.6718 [hep-ex] +RunInfo: + Inelastic events (non-diffractive and inelastic diffractive). +NumEvents: 100000 +NeedCrossSection: no +Beams: [p+, p+] +Energies: [7000] +PtCuts: [0] +Description: + 'The inelastic cross-section is measured through two independent methods + based on information from (i) forward calorimetry (for pseudorapidity + $3 < |\eta| < 5$), in collisions where at least one proton loses more than + $\xi > 5 \cdot 10^{-6}$ of its longitudinal momentum, and (ii) the central + tracker ($|\eta| < 2.4$), in collisions containing an interaction vertex + with more than 1, 2, or 3 tracks with $p_\perp > 200$ MeV/c.' +BibKey: :2012 +BibTeX: '@article{:2012, + author = "Chatrchyan, Serguei and others", + title = "{Measurement of the inelastic proton-proton cross section + at sqrt(s) = 7 TeV}", + collaboration = " CMS Collaboration", + year = "2012", + eprint = "1210.6718", + archivePrefix = "arXiv", + primaryClass = "hep-ex", + reportNumber = "CMS-FWD-11-001, CERN-PH-EP-2012-293", + SLACcitation = "%%CITATION = ARXIV:1210.6718;%%",}' Modified: branches/2012-06-aidarivet/data/anainfo/Makefile.am ============================================================================== --- branches/2012-06-aidarivet/data/anainfo/Makefile.am Mon Dec 17 17:34:17 2012 (r4089) +++ branches/2012-06-aidarivet/data/anainfo/Makefile.am Mon Dec 17 17:45:07 2012 (r4090) @@ -131,6 +131,7 @@ CMS_2012_I1102908.info \ CMS_2012_I1107658.info \ CMS_2012_I1184941.info \ + CMS_2012_I1193338.info \ CMS_2012_PAS_FWD_11_003.info \ CMS_2012_PAS_QCD_11_010.info \ CMS_QCD_10_024.info \ Added: branches/2012-06-aidarivet/data/plotinfo/CMS_2012_I1193338.plot ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ branches/2012-06-aidarivet/data/plotinfo/CMS_2012_I1193338.plot Mon Dec 17 17:45:07 2012 (r4090) @@ -0,0 +1,7 @@ +# BEGIN PLOT /CMS_2012_I1193338/d01-x01-y01 +Title=$\sigma_\text{inel}$ at $\sqrt{s}=7$ TeV +XCustomMajorTicks=0.5 $\xi>5x10^{-6}$ 1.5 $n_\text{ch} > 1$ 2.5 $n_\text{ch} > 2$ 3.5 $n_\text{ch} > 3$ +YLabel=$\sigma_\text{inel}$ [mb] +LegendYPos=0.4 +LogY=0 +# END PLOT Modified: branches/2012-06-aidarivet/data/plotinfo/Makefile.am ============================================================================== --- branches/2012-06-aidarivet/data/plotinfo/Makefile.am Mon Dec 17 17:34:17 2012 (r4089) +++ branches/2012-06-aidarivet/data/plotinfo/Makefile.am Mon Dec 17 17:45:07 2012 (r4090) @@ -128,6 +128,7 @@ CMS_2012_I1102908.plot \ CMS_2012_I1107658.plot \ CMS_2012_I1184941.plot \ + CMS_2012_I1193338.plot \ CMS_2012_PAS_FWD_11_003.plot \ CMS_2012_PAS_QCD_11_010.plot \ CMS_QCD_10_024.plot \ Added: branches/2012-06-aidarivet/data/refdata/CMS_2012_I1193338.aida ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ branches/2012-06-aidarivet/data/refdata/CMS_2012_I1193338.aida Mon Dec 17 17:45:07 2012 (r4090) @@ -0,0 +1,26 @@ +<?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.0" package="LWH"/> + <dataPointSet name="d01-x01-y01" + title="" path="/REF/CMS_2012_I1193338" dimension="2"> + <dimension dim="0" title="" /> + <dimension dim="1" title="" /> + <dataPoint> + <measurement value="0.5" errorPlus="0.5" errorMinus="0.5"/> + <measurement value="60.2" errorPlus="2.6476404589747453" errorMinus="2.6476404589747453"/> + </dataPoint> + <dataPoint> + <measurement value="1.5" errorPlus="0.5" errorMinus="0.5"/> + <measurement value="58.7" errorPlus="3.1240998703626617" errorMinus="3.1240998703626617"/> + </dataPoint> + <dataPoint> + <measurement value="2.5" errorPlus="0.5" errorMinus="0.5"/> + <measurement value="57.2" errorPlus="3.1240998703626617" errorMinus="3.1240998703626617"/> + </dataPoint> + <dataPoint> + <measurement value="3.5" errorPlus="0.5" errorMinus="0.5"/> + <measurement value="55.4" errorPlus="3.1240998703626617" errorMinus="3.1240998703626617"/> + </dataPoint> + </dataPointSet> +</aida> Modified: branches/2012-06-aidarivet/data/refdata/Makefile.am ============================================================================== --- branches/2012-06-aidarivet/data/refdata/Makefile.am Mon Dec 17 17:34:17 2012 (r4089) +++ branches/2012-06-aidarivet/data/refdata/Makefile.am Mon Dec 17 17:45:07 2012 (r4090) @@ -86,6 +86,7 @@ CMS_2012_I1102908.aida \ CMS_2012_I1107658.aida \ CMS_2012_I1184941.aida \ + CMS_2012_I1193338.aida \ CMS_2012_PAS_FWD_11_003.aida \ CMS_2012_PAS_QCD_11_010.aida \ CMS_QCD_10_024.aida \ Added: branches/2012-06-aidarivet/src/Analyses/CMS_2012_I1193338.cc ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ branches/2012-06-aidarivet/src/Analyses/CMS_2012_I1193338.cc Mon Dec 17 17:45:07 2012 (r4090) @@ -0,0 +1,88 @@ +// -*- C++ -*- +#include "Rivet/Analysis.hh" +#include "Rivet/RivetAIDA.hh" +#include "Rivet/Projections/ChargedFinalState.hh" +#include "Rivet/Projections/FinalState.hh" + +namespace Rivet { + + + class CMS_2012_I1193338 : public Analysis { + public: + + CMS_2012_I1193338() + : Analysis("CMS_2012_I1193338") + { } + + public: + + void init() { + addProjection(ChargedFinalState(-2.4, 2.4, 0.2*GeV), "CFS"); + addProjection(FinalState(), "FS"); + + _h_sigma = bookHistogram1D(1, 1, 1); + } + + + void analyze(const Event& event) { + const double weight = event.weight(); + + const ChargedFinalState& cfs = applyProjection<ChargedFinalState>(event, "CFS"); + if (cfs.size() > 1) {_h_sigma->fill(1.5, weight);} + if (cfs.size() > 2) {_h_sigma->fill(2.5, weight);} + if (cfs.size() > 3) {_h_sigma->fill(3.5, weight);} + + const FinalState& fs = applyProjection<FinalState>(event, "FS"); + if (fs.size() < 2) vetoEvent; // need at least two particles to calculate gaps + + double gapcenter = 0.; + double LRG = 0.; + double etapre = 0.; + bool first = true; + + foreach(const Particle& p, fs.particlesByEta()) { // sorted from minus to plus + if (first) { // First particle + first = false; + etapre = p.momentum().eta(); + } else { + double gap = fabs(p.momentum().eta()-etapre); + if (gap > LRG) { + LRG = gap; // largest gap + gapcenter = (p.momentum().eta()+etapre)/2.; // find the center of the gap to separate the X and Y systems. + } + etapre = p.momentum().eta(); + } + } + + + FourMomentum mxFourVector, myFourVector; + foreach(const Particle& p, fs.particlesByEta()) { + if (p.momentum().eta() > gapcenter) { + mxFourVector += p.momentum(); + } else { + myFourVector += p.momentum(); + } + } + const double M2 = max(mxFourVector.mass2(), myFourVector.mass2()); + const double xi = M2/sqr(sqrtS()); // sqrt(s)=7000 GeV, note that units cancel + if (xi < 5e-6) vetoEvent; + + _h_sigma->fill(0.5, weight); + } + + + void finalize() { + scale(_h_sigma, crossSection()/millibarn/sumOfWeights()); + } + + private: + + AIDA::IHistogram1D* _h_sigma; + + }; + + + // The hook for the plugin system + DECLARE_RIVET_PLUGIN(CMS_2012_I1193338); + +} Modified: branches/2012-06-aidarivet/src/Analyses/Makefile.am ============================================================================== --- branches/2012-06-aidarivet/src/Analyses/Makefile.am Mon Dec 17 17:34:17 2012 (r4089) +++ branches/2012-06-aidarivet/src/Analyses/Makefile.am Mon Dec 17 17:45:07 2012 (r4090) @@ -129,7 +129,8 @@ CMS_2012_I1087342.cc \ CMS_2012_I1102908.cc \ CMS_2012_I1107658.cc \ - CMS_2012_I1184941.cc + CMS_2012_I1184941.cc \ + CMS_2012_I1193338.cc if ENABLE_PRELIMINARY RivetCMSAnalyses_la_SOURCES += \ CMS_QCD_10_024.cc \
More information about the Rivet-svn mailing list |