|
[Rivet-svn] r4023 - in branches/2012-06-aidarivet: . data/anainfo data/plotinfo data/refdata src/Analysesblackhole at projects.hepforge.org blackhole at projects.hepforge.orgThu Nov 15 15:51:48 GMT 2012
Author: hoeth Date: Thu Nov 15 15:51:47 2012 New Revision: 4023 Log: Added ATLAS_2011_I894867 analysis Added: branches/2012-06-aidarivet/data/anainfo/ATLAS_2011_I894867.info branches/2012-06-aidarivet/data/plotinfo/ATLAS_2011_I894867.plot branches/2012-06-aidarivet/data/refdata/ATLAS_2011_I894867.aida branches/2012-06-aidarivet/src/Analyses/ATLAS_2011_I894867.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 Nov 15 15:17:03 2012 (r4022) +++ branches/2012-06-aidarivet/ChangeLog Thu Nov 15 15:51:47 2012 (r4023) @@ -1,6 +1,7 @@ 2012-11-15 Hendrik Hoeth <hendrik.hoeth at cern.ch> * Added LHCF_2012_I1115479 analysis + * Added ATLAS_2011_I894867 analysis 2012-11-14 Hendrik Hoeth <hendrik.hoeth at cern.ch> Added: branches/2012-06-aidarivet/data/anainfo/ATLAS_2011_I894867.info ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ branches/2012-06-aidarivet/data/anainfo/ATLAS_2011_I894867.info Thu Nov 15 15:51:47 2012 (r4023) @@ -0,0 +1,36 @@ +Name: ATLAS_2011_I894867 +Year: 2011 +Summary: Measurement of the inelastic proton-proton cross-section at sqrt{s}=7 TeV. +Experiment: ATLAS +Collider: LHC +SpireID: 9021310 +InspireID: 894867 +Status: VALIDATED +Authors: + - Sercan Sen <ssen at cern.ch> +References: + - arXiv:1104.0326 [hep-ex] +RunInfo: + Inelastic events (non-diffractive and inelastic diffractive). +NumEvents: 100K +Beams: [p+, p+] +Energies: [7000] +Description: + Inelastic cross-section is measured for $\xi > 5x10^{-6}$, where xi=M_X^2/s is calculated from the invariant mass, M_X, + of hadrons selected using the largest rapidity gap in the event. +BibKey: 2011eu +BibTeX: '@Article{Aad:2011eu, + author = "Aad, Georges and others", + title = "{Measurement of the Inelastic Proton-Proton Cross-Section + at $\sqrt{s}=7$ TeV with the ATLAS Detector}", + collaboration = "ATLAS Collaboration", + journal = "Nature Commun.", + volume = "2", + pages = "463", + doi = "10.1038/ncomms1472", + year = "2011", + eprint = "1104.0326", + archivePrefix = "arXiv", + primaryClass = "hep-ex", + reportNumber = "CERN-PH-EP-2011-047", + SLACcitation = "%%CITATION = ARXIV:1104.0326;%%",}' Modified: branches/2012-06-aidarivet/data/anainfo/Makefile.am ============================================================================== --- branches/2012-06-aidarivet/data/anainfo/Makefile.am Thu Nov 15 15:17:03 2012 (r4022) +++ branches/2012-06-aidarivet/data/anainfo/Makefile.am Thu Nov 15 15:51:47 2012 (r4023) @@ -36,6 +36,7 @@ ATLAS_2011_S9131140.info \ ATLAS_2011_S9108483.info \ ATLAS_2011_S9212183.info \ + ATLAS_2011_I894867.info \ ATLAS_2011_I919017.info \ ATLAS_2011_I925932.info \ ATLAS_2011_I926145.info \ Added: branches/2012-06-aidarivet/data/plotinfo/ATLAS_2011_I894867.plot ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ branches/2012-06-aidarivet/data/plotinfo/ATLAS_2011_I894867.plot Thu Nov 15 15:51:47 2012 (r4023) @@ -0,0 +1,7 @@ +# BEGIN PLOT /ATLAS_2011_I894867/d01-x01-y01 +Title=$\sigma_\text{inel}$ for $\xi > 5 \cdot 10^{-6}$ at $\sqrt{s}=7$ TeV +XLabel=$\sqrt{s}$ [GeV] +YLabel=$\sigma_\text{inel}$ [mb] +YMin=45 +XCustomMajorTicks=7000 7000 +# END PLOT Modified: branches/2012-06-aidarivet/data/plotinfo/Makefile.am ============================================================================== --- branches/2012-06-aidarivet/data/plotinfo/Makefile.am Thu Nov 15 15:17:03 2012 (r4022) +++ branches/2012-06-aidarivet/data/plotinfo/Makefile.am Thu Nov 15 15:51:47 2012 (r4023) @@ -29,6 +29,7 @@ ATLAS_2011_S9126244.plot \ ATLAS_2011_S9128077.plot \ ATLAS_2011_S9131140.plot \ + ATLAS_2011_I894867.plot \ ATLAS_2011_I919017.plot \ ATLAS_2011_I925932.plot \ ATLAS_2011_I926145.plot \ Added: branches/2012-06-aidarivet/data/refdata/ATLAS_2011_I894867.aida ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ branches/2012-06-aidarivet/data/refdata/ATLAS_2011_I894867.aida Thu Nov 15 15:51:47 2012 (r4023) @@ -0,0 +1,14 @@ +<?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/ATLAS_2011_I894867" dimension="2"> + <dimension dim="0" title="" /> + <dimension dim="1" title="" /> + <dataPoint> + <measurement value="7000.0" errorPlus="0.5" errorMinus="0.5"/> + <measurement value="60.33" errorPlus="2.10" errorMinus="2.10"/> + </dataPoint> + </dataPointSet> +</aida> Modified: branches/2012-06-aidarivet/data/refdata/Makefile.am ============================================================================== --- branches/2012-06-aidarivet/data/refdata/Makefile.am Thu Nov 15 15:17:03 2012 (r4022) +++ branches/2012-06-aidarivet/data/refdata/Makefile.am Thu Nov 15 15:51:47 2012 (r4023) @@ -36,6 +36,7 @@ ATLAS_2011_S9126244.aida \ ATLAS_2011_S9128077.aida \ ATLAS_2011_S9131140.aida \ + ATLAS_2011_I894867.aida \ ATLAS_2011_I919017.aida \ ATLAS_2011_I925932.aida \ ATLAS_2011_I926145.aida \ Added: branches/2012-06-aidarivet/src/Analyses/ATLAS_2011_I894867.cc ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ branches/2012-06-aidarivet/src/Analyses/ATLAS_2011_I894867.cc Thu Nov 15 15:51:47 2012 (r4023) @@ -0,0 +1,86 @@ +// -*- C++ -*- +#include "Rivet/Analysis.hh" +#include "Rivet/RivetAIDA.hh" +#include "Rivet/Projections/FinalState.hh" + +namespace Rivet { + + + class ATLAS_2011_I894867 : public Analysis { + public: + + ATLAS_2011_I894867() + : Analysis("ATLAS_2011_I894867") + { } + + public: + + void init() { + addProjection(FinalState(),"FS"); + _h_sigma = bookHistogram1D(1, 1, 1); + } + + + void analyze(const Event& event) { + const double weight = event.weight(); + + const FinalState& fs = applyProjection<FinalState>(event, "FS"); + + 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(0.,0.,0.,0.); + FourMomentum MyFourVector(0.,0.,0.,0.); + + foreach(const Particle& p, fs.particlesByEta()) { + if (p.momentum().eta() > gapcenter) { + MxFourVector += p.momentum(); + } else { + MyFourVector += p.momentum(); + } + } + + double Mx2 = FourMomentum(MxFourVector).mass2(); + double My2 = FourMomentum(MyFourVector).mass2(); + + const double M2 = (Mx2 > My2 ? Mx2 : My2); + const double xi = M2/(7000*7000); // sqrt(s)=7000 GeV + + if (xi < 5*10e-6) vetoEvent; + + _h_sigma->fill(7000/GeV, weight); + + } + + + void finalize() { + scale(_h_sigma, crossSection()/millibarn/sumOfWeights()); + } + + private: + + AIDA::IHistogram1D* _h_sigma; + + }; + + + // The hook for the plugin system + DECLARE_RIVET_PLUGIN(ATLAS_2011_I894867); + +} Modified: branches/2012-06-aidarivet/src/Analyses/Makefile.am ============================================================================== --- branches/2012-06-aidarivet/src/Analyses/Makefile.am Thu Nov 15 15:17:03 2012 (r4022) +++ branches/2012-06-aidarivet/src/Analyses/Makefile.am Thu Nov 15 15:51:47 2012 (r4023) @@ -63,6 +63,7 @@ ATLAS_2011_S9212183.cc \ ATLAS_2011_S9225137.cc \ ATLAS_2011_S9019561.cc \ + ATLAS_2011_I894867.cc \ ATLAS_2011_I919017.cc \ ATLAS_2011_I925932.cc \ ATLAS_2011_I926145.cc \
More information about the Rivet-svn mailing list |