|
[Rivet-svn] r4061 - in branches/2012-06-aidarivet: . data/anainfo data/plotinfo data/refdata src/Analysesblackhole at projects.hepforge.org blackhole at projects.hepforge.orgThu Dec 6 15:53:49 GMT 2012
Author: hoeth Date: Thu Dec 6 15:53:49 2012 New Revision: 4061 Log: Added ATLAS_2012_I1118269 analysis Added: branches/2012-06-aidarivet/data/anainfo/ATLAS_2012_I1118269.info branches/2012-06-aidarivet/data/plotinfo/ATLAS_2012_I1118269.plot branches/2012-06-aidarivet/data/refdata/ATLAS_2012_I1118269.aida branches/2012-06-aidarivet/src/Analyses/ATLAS_2012_I1118269.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:26:38 2012 (r4060) +++ branches/2012-06-aidarivet/ChangeLog Thu Dec 6 15:53:49 2012 (r4061) @@ -1,6 +1,7 @@ 2012-12-06 Hendrik Hoeth <hendrik.hoeth at cern.ch> * Added ATLAS_2012_I1188891 analysis + * Added ATLAS_2012_I1118269 analysis * Added LHCB_2010_I867355 analysis * Added TGraphErrors support to root2flat Added: branches/2012-06-aidarivet/data/anainfo/ATLAS_2012_I1118269.info ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ branches/2012-06-aidarivet/data/anainfo/ATLAS_2012_I1118269.info Thu Dec 6 15:53:49 2012 (r4061) @@ -0,0 +1,42 @@ +Name: ATLAS_2012_I1118269 +Year: 2012 +Summary: Measurement of the b-hadron production cross section using decays to $D^{*}\mu^-X$ final states in pp collisions at sqrt(s) = 7 TeV with the ATLAS detector +Experiment: ATLAS +Collider: LHC +SpiresID: 9626190 +InspireID: 1118269 +Status: VALIDATED +Authors: + - Andy Buckley <andy.buckley at ed.ac.uk> + - Sercan Sen <sercan.sen at cern.ch> + - Peter Skands <Peter.Skands at cern.ch> +References: + - arXiv:1206.3122 [hep-ex] +RunInfo: + pp to b-hadron + X at 7 TeV, i.e. switch on "HardQCD:gg2bbbar" and "HardQCD:qqbar2bbbar" flags in Pythia8. +NumEvents: 1000000 +Beams: [p+, p+] +Energies: [7000] +PtCuts: [0] +Description: + 'Measurement of $b$-hadron production cross section using 3.3 pb$^{-1}$ of integrated luminosity, collected during the 2010 LHC run. + The $b$-hadrons are selected by partially reconstructing $D^{*}\mu^-X$ final states using only direct semileptonic decays of $b$ to $D^{*}\mu^-X$. + Differential cross sections as functions of $p_\perp$ and $|\eta|$.' +BibKey: 2012fc +BibTeX: '@article{:2012fc, + author = "Aad, Georges and others", + title = "{Measurement of the b-hadron production cross section + using decays to $D^{*}\mu^-X$ final states in pp + collisions at sqrt(s) = 7 TeV with the ATLAS detector}", + collaboration = "ATLAS Collaboration", + journal = "Nucl.Phys.", + volume = "B864", + pages = "341-381", + doi = "10.1016/j.nuclphysb.2012.07.009", + year = "2012", + eprint = "1206.3122", + archivePrefix = "arXiv", + primaryClass = "hep-ex", + reportNumber = "CERN-PH-EP-2012-121", + SLACcitation = "%%CITATION = ARXIV:1206.3122;%%", +}' Modified: branches/2012-06-aidarivet/data/anainfo/Makefile.am ============================================================================== --- branches/2012-06-aidarivet/data/anainfo/Makefile.am Thu Dec 6 15:26:38 2012 (r4060) +++ branches/2012-06-aidarivet/data/anainfo/Makefile.am Thu Dec 6 15:53:49 2012 (r4061) @@ -62,6 +62,7 @@ ATLAS_2012_I1125961.info \ ATLAS_2012_I1126136.info \ ATLAS_2012_I1117704.info \ + ATLAS_2012_I1118269.info \ ATLAS_2012_I1180197.info \ ATLAS_2012_I1186556.info \ ATLAS_2012_I1190891.info \ Added: branches/2012-06-aidarivet/data/plotinfo/ATLAS_2012_I1118269.plot ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ branches/2012-06-aidarivet/data/plotinfo/ATLAS_2012_I1118269.plot Thu Dec 6 15:53:49 2012 (r4061) @@ -0,0 +1,13 @@ +# BEGIN PLOT /ATLAS_2012_I1118269/d01-x01-y01 +Title=$b$ hadron $p_\perp$ at $\sqrt{s}=7$ TeV +XLabel=$p_\perp(H_b)$ [GeV] +YLabel=$\mathrm{d}\sigma(H_bX)/\mathrm{d} p_\perp(H_b)$ [nb/GeV] +# END PLOT + +# BEGIN PLOT /ATLAS_2012_I1118269/d02-x01-y01 +Title=$b$ hadron $\eta$ at $\sqrt{s}=7$ TeV +XLabel=$|\eta|(H_b)$ +YLabel=$\mathrm{d}\sigma(H_bX) / \mathrm{d}|\eta|(H_b)$ [$\mu$b] +LogY=0 +LegendYPos=0.4 +# END PLOT Modified: branches/2012-06-aidarivet/data/plotinfo/Makefile.am ============================================================================== --- branches/2012-06-aidarivet/data/plotinfo/Makefile.am Thu Dec 6 15:26:38 2012 (r4060) +++ branches/2012-06-aidarivet/data/plotinfo/Makefile.am Thu Dec 6 15:53:49 2012 (r4061) @@ -56,6 +56,7 @@ ATLAS_2012_I1125961.plot \ ATLAS_2012_I1126136.plot \ ATLAS_2012_I1117704.plot \ + ATLAS_2012_I1118269.plot \ ATLAS_2012_I1180197.plot \ ATLAS_2012_I1082009.plot \ ATLAS_2012_I1183818.plot \ Added: branches/2012-06-aidarivet/data/refdata/ATLAS_2012_I1118269.aida ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ branches/2012-06-aidarivet/data/refdata/ATLAS_2012_I1118269.aida Thu Dec 6 15:53:49 2012 (r4061) @@ -0,0 +1,59 @@ +<?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_2012_I1118269" dimension="2"> + <dimension dim="0" title="" /> + <dimension dim="1" title="" /> + <dataPoint> + <measurement value="11.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="5800.0" errorPlus="2.2" errorMinus="2.2"/> + </dataPoint> + <dataPoint> + <measurement value="13.5" errorPlus="1.5" errorMinus="1.5"/> + <measurement value="2300.0" errorPlus="2.2" errorMinus="2.2"/> + </dataPoint> + <dataPoint> + <measurement value="17.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="900.0" errorPlus="2.2" errorMinus="2.2"/> + </dataPoint> + <dataPoint> + <measurement value="25.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="220.0" errorPlus="2.2" errorMinus="2.2"/> + </dataPoint> + <dataPoint> + <measurement value="37.5" errorPlus="7.5" errorMinus="7.5"/> + <measurement value="30.0" errorPlus="2.2" errorMinus="2.2"/> + </dataPoint> + <dataPoint> + <measurement value="62.5" errorPlus="17.5" errorMinus="17.5"/> + <measurement value="2.8" errorPlus="0.2" errorMinus="0.2"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d02-x01-y01" + title="" path="/REF/ATLAS_2012_I1118269" dimension="2"> + <dimension dim="0" title="" /> + <dimension dim="1" title="" /> + <dataPoint> + <measurement value="0.25" errorPlus="0.25" errorMinus="0.25"/> + <measurement value="14.3" errorPlus="1.0" errorMinus="1.0"/> + </dataPoint> + <dataPoint> + <measurement value="0.75" errorPlus="0.25" errorMinus="0.25"/> + <measurement value="13.2" errorPlus="1.0" errorMinus="1.0"/> + </dataPoint> + <dataPoint> + <measurement value="1.25" errorPlus="0.25" errorMinus="0.25"/> + <measurement value="13.0" errorPlus="1.2" errorMinus="1.2"/> + </dataPoint> + <dataPoint> + <measurement value="1.75" errorPlus="0.25" errorMinus="0.25"/> + <measurement value="16.0" errorPlus="1.4" errorMinus="1.4"/> + </dataPoint> + <dataPoint> + <measurement value="2.25" errorPlus="0.25" errorMinus="0.25"/> + <measurement value="13.2" errorPlus="1.5" errorMinus="1.5"/> + </dataPoint> + </dataPointSet> +</aida> Modified: branches/2012-06-aidarivet/data/refdata/Makefile.am ============================================================================== --- branches/2012-06-aidarivet/data/refdata/Makefile.am Thu Dec 6 15:26:38 2012 (r4060) +++ branches/2012-06-aidarivet/data/refdata/Makefile.am Thu Dec 6 15:53:49 2012 (r4061) @@ -54,6 +54,7 @@ ATLAS_2012_I1094568.aida \ ATLAS_2012_I943401.aida \ ATLAS_2012_I1082009.aida \ + ATLAS_2012_I1118269.aida \ ATLAS_2012_I1183818.aida \ ATLAS_2012_I1188891.aida \ ATLAS_2012_CONF_2012_001.aida \ Added: branches/2012-06-aidarivet/src/Analyses/ATLAS_2012_I1118269.cc ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ branches/2012-06-aidarivet/src/Analyses/ATLAS_2012_I1118269.cc Thu Dec 6 15:53:49 2012 (r4061) @@ -0,0 +1,80 @@ +// -*- C++ -*- +#include "Rivet/Analysis.hh" +#include "Rivet/RivetAIDA.hh" +#include "Rivet/Tools/Logging.hh" +#include "Rivet/Particle.hh" +#include "Rivet/Tools/ParticleIdUtils.hh" + +namespace Rivet { + + class ATLAS_2012_I1118269 : public Analysis { + public: + + ATLAS_2012_I1118269() : Analysis("ATLAS_2012_I1118269") + { } + + void init() { + _h_sigma_vs_pt = bookHistogram1D(1, 1, 1); + _h_sigma_vs_eta = bookHistogram1D(2, 1, 1); + } + + /// Perform the per-event analysis + void analyze(const Event& event) { + double weight = event.weight(); + + ParticleVector bhadrons; + foreach (const GenParticle* p, particles(event.genEvent())) { + + if (!( PID::isHadron( p->pdg_id() ) && PID::hasBottom( p->pdg_id() )) ) continue; + + const GenVertex* dv = p->end_vertex(); + + /// @todo In future, convert to use built-in 'last B hadron' function + bool hasBdaughter = false; + if ( PID::isHadron( p->pdg_id() ) && PID::hasBottom( p->pdg_id() )) { // b-hadron selection + if (dv) { + for (GenVertex::particles_out_const_iterator pp = dv->particles_out_const_begin() ; pp != dv->particles_out_const_end() ; ++pp) { + if ( PID::isHadron( (*pp)->pdg_id() ) && PID::hasBottom( (*pp)->pdg_id()) ) { + hasBdaughter = true; + } + } + } + } + if (hasBdaughter) continue; + + bhadrons += Particle(*p); + } + + foreach (const Particle& particle, bhadrons) { + double eta = particle.momentum().eta(); + double pt = particle.momentum().pT(); + + if (!(inRange(eta, -2.5, 2.5))) continue; + if (pt < 9.*GeV) continue; + + _h_sigma_vs_pt->fill(pt, weight); + _h_sigma_vs_eta->fill(fabs(eta), weight); + + } + + } + + + void finalize() { + scale(_h_sigma_vs_pt, crossSection()/nanobarn/sumOfWeights()); + scale(_h_sigma_vs_eta, crossSection()/microbarn/sumOfWeights()); + } + + + private: + + AIDA::IHistogram1D *_h_sigma_vs_pt; + AIDA::IHistogram1D *_h_sigma_vs_eta; + + }; + + + // Hook for the plugin system + DECLARE_RIVET_PLUGIN(ATLAS_2012_I1118269); + +} Modified: branches/2012-06-aidarivet/src/Analyses/Makefile.am ============================================================================== --- branches/2012-06-aidarivet/src/Analyses/Makefile.am Thu Dec 6 15:26:38 2012 (r4060) +++ branches/2012-06-aidarivet/src/Analyses/Makefile.am Thu Dec 6 15:53:49 2012 (r4061) @@ -78,6 +78,7 @@ ATLAS_2012_I1094568.cc \ ATLAS_2012_I1112263.cc \ ATLAS_2012_I1117704.cc \ + ATLAS_2012_I1118269.cc \ ATLAS_2012_I1125961.cc \ ATLAS_2012_I1183818.cc \ ATLAS_2012_I1188891.cc
More information about the Rivet-svn mailing list |