|
|
[Rivet-svn] r4008 - in trunk: data/anainfo data/plotinfo src/Analyses
blackhole at projects.hepforge.org
blackhole at projects.hepforge.org
Wed Nov 14 17:52:02 GMT 2012
Author: buckley
Date: Wed Nov 14 17:52:02 2012
New Revision: 4008
Log:
Merge changeset [3972] from AIDA branch and YODAifying
Added:
trunk/data/anainfo/ATLAS_2012_I1126136.info
- copied unchanged from r3972, branches/2012-06-aidarivet/data/anainfo/ATLAS_2012_I1126136.info
trunk/data/plotinfo/ATLAS_2012_I1126136.plot
- copied unchanged from r3972, branches/2012-06-aidarivet/data/plotinfo/ATLAS_2012_I1126136.plot
trunk/src/Analyses/ATLAS_2012_I1126136.cc
- copied, changed from r3972, branches/2012-06-aidarivet/src/Analyses/ATLAS_2012_I1126136.cc
Modified:
trunk/data/anainfo/Makefile.am
trunk/data/plotinfo/Makefile.am
trunk/src/Analyses/ATLAS_2012_I1186556.cc
trunk/src/Analyses/Makefile.am
Copied: trunk/data/anainfo/ATLAS_2012_I1126136.info (from r3972, branches/2012-06-aidarivet/data/anainfo/ATLAS_2012_I1126136.info)
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ trunk/data/anainfo/ATLAS_2012_I1126136.info Wed Nov 14 17:52:02 2012 (r4008, copy of r3972, branches/2012-06-aidarivet/data/anainfo/ATLAS_2012_I1126136.info)
@@ -0,0 +1,35 @@
+Name: ATLAS_2012_I1126136
+Year: 2012
+Summary: SUSY Top partner search in jets with missing transverse momentum
+Experiment: ATLAS
+Collider: LHC
+InspireID: 1126136
+Status: UNVALIDATED
+Authors:
+ - Peter Richardson <peter.richardson at durham.ac.uk>
+References:
+ - arXiv:1208.1447
+RunInfo:
+ BSM signal events at 7000 GeV.
+NumEvents: 25000 for BSM signals
+Beams: [p+, p+]
+Energies: [7000]
+NeedCrossSection: yes
+Description:
+ 'Search for direct pair production of supersymmetric top squarks, assuming the stop_1 decays
+ into a top quark and the lightest supersymmetric particle, and that both top quarks decay
+ to purely hadronic final states. This search has an integrated luminosity
+ of 4.7 fb$^{-1}$ at $\sqrt{s}=8$\,TeV.'
+BibTeX: '@article{:2012si,
+ author = "Aad, Georges and others",
+ title = "{Search for a supersymmetric partner to the top quark in
+ final states with jets and missing transverse momentum at
+ sqrt(s) = 7 TeV with the ATLAS detector}",
+ collaboration = "ATLAS Collaboration",
+ year = "2012",
+ eprint = "1208.1447",
+ archivePrefix = "arXiv",
+ primaryClass = "hep-ex",
+ reportNumber = "CERN-PH-EP-2012-201",
+ SLACcitation = "%%CITATION = ARXIV:1208.1447;%%",
+}'
Modified: trunk/data/anainfo/Makefile.am
==============================================================================
--- trunk/data/anainfo/Makefile.am Wed Nov 14 17:49:45 2012 (r4007)
+++ trunk/data/anainfo/Makefile.am Wed Nov 14 17:52:02 2012 (r4008)
@@ -58,6 +58,7 @@
ATLAS_2012_CONF_2012_001.info \
ATLAS_2012_I1112263.info \
ATLAS_2012_I1125961.info \
+ ATLAS_2012_I1126136.info \
ATLAS_2012_I1117704.info \
ATLAS_2012_I1180197.info \
ATLAS_2012_I1186556.info \
@@ -215,3 +216,4 @@
UA5_1987_S1640666.info \
UA5_1988_S1867512.info \
UA5_1989_S1926373.info
+
Copied: trunk/data/plotinfo/ATLAS_2012_I1126136.plot (from r3972, branches/2012-06-aidarivet/data/plotinfo/ATLAS_2012_I1126136.plot)
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ trunk/data/plotinfo/ATLAS_2012_I1126136.plot Wed Nov 14 17:52:02 2012 (r4008, copy of r3972, branches/2012-06-aidarivet/data/plotinfo/ATLAS_2012_I1126136.plot)
@@ -0,0 +1,36 @@
+# BEGIN PLOT /ATLAS_2012_I1186556/count_SR_A
+Title=Number of events in the $E_{\mathrm{T}}^{\rm miss}>130\,$GeV
+XLabel=
+YLabel=Number of Events
+LogY=0
+# END PLOT
+# BEGIN PLOT /ATLAS_2012_I1186556/count_SR_B
+Title=Number of events in the $E_{\mathrm{T}}^{\rm miss}>260\,$GeV
+XLabel=
+YLabel=Number of Events
+LogY=0
+# END PLOT
+# BEGIN PLOT /ATLAS_2012_I1186556/hist_mjjj1
+Title=$m_{jjj}$ for closest-triplet
+XLabel=$m_{jjj}$ [GeV]
+YLabel=$\mathrm{Events}/20 \mathrm{GeV}$
+LogY=0
+# END PLOT
+# BEGIN PLOT /ATLAS_2012_I1186556/hist_mjjj2
+Title=$m_{jjj}$ for second-triplet
+XLabel=$m_{jjj}$ [GeV]
+YLabel=$\mathrm{Events}/20 \mathrm{GeV}$
+LogY=0
+# END PLOT
+# BEGIN PLOT /ATLAS_2012_I1186556/hist_ETmiss
+Title=$E_{\mathrm{T}}^{\rm miss}$
+XLabel=$E_{\mathrm{T}}^{\rm miss}$ [GeV]
+YLabel=$\mathrm{Events}/50 \mathrm{GeV}$
+LogY=0
+# END PLOT
+# BEGIN PLOT /ATLAS_2012_I1186556/hist_mT2
+Title=$m_{T2}$
+XLabel=$m_{T2}$ [GeV]
+YLabel=$\mathrm{Events}/\mathrm{GeV}$
+LogY=0
+# END PLOT
Modified: trunk/data/plotinfo/Makefile.am
==============================================================================
--- trunk/data/plotinfo/Makefile.am Wed Nov 14 17:49:45 2012 (r4007)
+++ trunk/data/plotinfo/Makefile.am Wed Nov 14 17:52:02 2012 (r4008)
@@ -53,6 +53,7 @@
ATLAS_2012_CONF_2012_001.plot \
ATLAS_2012_I1112263.plot \
ATLAS_2012_I1125961.plot \
+ ATLAS_2012_I1126136.plot \
ATLAS_2012_I1117704.plot \
ATLAS_2012_I1180197.plot \
ATLAS_2012_I1082009.plot \
Copied and modified: trunk/src/Analyses/ATLAS_2012_I1126136.cc (from r3972, branches/2012-06-aidarivet/src/Analyses/ATLAS_2012_I1126136.cc)
==============================================================================
--- branches/2012-06-aidarivet/src/Analyses/ATLAS_2012_I1126136.cc Thu Oct 25 12:25:01 2012 (r3972, copy source)
+++ trunk/src/Analyses/ATLAS_2012_I1126136.cc Wed Nov 14 17:52:02 2012 (r4008)
@@ -1,7 +1,7 @@
// -*- C++ -*-
#include "Rivet/Analysis.hh"
#include "Rivet/Tools/BinnedHistogram.hh"
-#include "Rivet/RivetAIDA.hh"
+#include "Rivet/RivetYODA.hh"
#include "Rivet/Tools/Logging.hh"
#include "Rivet/Projections/FinalState.hh"
#include "Rivet/Projections/ChargedFinalState.hh"
@@ -39,14 +39,14 @@
void init() {
// projection to find the electrons
- std::vector<std::pair<double, double> > eta_e;
+ vector<pair<double, double> > eta_e;
eta_e.push_back(make_pair(-2.47,2.47));
IdentifiedFinalState elecs(eta_e, 20.0*GeV);
elecs.acceptIdPair(ELECTRON);
addProjection(elecs, "elecs");
// projection to find the muons
- std::vector<std::pair<double, double> > eta_m;
+ vector<pair<double, double> > eta_m;
eta_m.push_back(make_pair(-2.4,2.4));
IdentifiedFinalState muons(eta_m, 10.0*GeV);
muons.acceptIdPair(MUON);
@@ -62,13 +62,13 @@
addProjection(VisibleFinalState(-4.9,4.9),"vfs");
// Book histograms
- _count_SR_A = bookHistogram1D("count_SR_A" , 1, 0., 1.);
- _count_SR_B = bookHistogram1D("count_SR_B" , 1, 0., 1.);
-
- _hist_mjjj1 = bookHistogram1D("hist_mjjj1" , 30 , 0. , 600. );
- _hist_mjjj2 = bookHistogram1D("hist_mjjj2" , 30 , 0. , 600. );
- _hist_ETmiss = bookHistogram1D("hist_ETmiss", 20 , 100. , 600. );
- _hist_mT2 = bookHistogram1D("hist_mT2" , 200, 0. , 1000. );
+ _count_SR_A = bookHisto1D("count_SR_A" , 1, 0., 1.);
+ _count_SR_B = bookHisto1D("count_SR_B" , 1, 0., 1.);
+
+ _hist_mjjj1 = bookHisto1D("hist_mjjj1" , 30 , 0. , 600. );
+ _hist_mjjj2 = bookHisto1D("hist_mjjj2" , 30 , 0. , 600. );
+ _hist_ETmiss = bookHisto1D("hist_ETmiss", 20 , 100. , 600. );
+ _hist_mT2 = bookHisto1D("hist_mT2" , 200, 0. , 1000. );
}
@@ -225,7 +225,7 @@
for(unsigned int i1=0;i1<recon_jets.size();++i1) {
for(unsigned int i2=i1+1;i2<recon_jets.size();++i2) {
for(unsigned int i3=i2+1;i3<recon_jets.size();++i3) {
- double delR2 =
+ double delR2 =
sqr(deltaR(recon_jets[i1].momentum(),recon_jets[i2].momentum())) +
sqr(deltaR(recon_jets[i1].momentum(),recon_jets[i3].momentum())) +
sqr(deltaR(recon_jets[i2].momentum(),recon_jets[i3].momentum()));
@@ -239,7 +239,7 @@
}
}
// 4-momentum and mass of first triplet
- FourMomentum pjjj1 = recon_jets[j1].momentum() +
+ FourMomentum pjjj1 = recon_jets[j1].momentum() +
recon_jets[j2].momentum()+ recon_jets[j3].momentum();
double mjjj1 = pjjj1.mass();
@@ -252,7 +252,7 @@
if(i2==j1||i2==j2||i2==j3) continue;
for(unsigned int i3=i2+1;i3<recon_jets.size();++i3) {
if(i3==j1||i3==j2||i3==j3) continue;
- double delR2 =
+ double delR2 =
sqr(deltaR(recon_jets[i1].momentum(),recon_jets[i2].momentum())) +
sqr(deltaR(recon_jets[i1].momentum(),recon_jets[i3].momentum())) +
sqr(deltaR(recon_jets[i2].momentum(),recon_jets[i3].momentum()));
@@ -267,7 +267,7 @@
}
// 4-momentum and mass of first triplet
- FourMomentum pjjj2 = recon_jets[j4].momentum() +
+ FourMomentum pjjj2 = recon_jets[j4].momentum() +
recon_jets[j5].momentum()+ recon_jets[j6].momentum();
double mjjj2 = pjjj2.mass();
@@ -305,13 +305,13 @@
/// @name Histograms
//@{
- AIDA::IHistogram1D* _count_SR_A;
- AIDA::IHistogram1D* _count_SR_B;
+ Histo1DPtr _count_SR_A;
+ Histo1DPtr _count_SR_B;
- AIDA::IHistogram1D* _hist_mjjj1;
- AIDA::IHistogram1D* _hist_mjjj2;
- AIDA::IHistogram1D* _hist_ETmiss;
- AIDA::IHistogram1D* _hist_mT2;
+ Histo1DPtr _hist_mjjj1;
+ Histo1DPtr _hist_mjjj2;
+ Histo1DPtr _hist_ETmiss;
+ Histo1DPtr _hist_mT2;
//@}
};
Modified: trunk/src/Analyses/ATLAS_2012_I1186556.cc
==============================================================================
--- trunk/src/Analyses/ATLAS_2012_I1186556.cc Wed Nov 14 17:49:45 2012 (r4007)
+++ trunk/src/Analyses/ATLAS_2012_I1186556.cc Wed Nov 14 17:52:02 2012 (r4008)
@@ -196,7 +196,7 @@
// require at least 1 b jet
unsigned int n_b=0;
for(unsigned int ix=0;ix<recon_jets.size();++ix) {
- if(recon_jets[ix].containsBottom() && rand()/static_cast<double>(RAND_MAX)>=0.60)
+ if(recon_jets[ix].containsBottom() && rand()/static_cast<double>(RAND_MAX)<=0.60)
++n_b;
}
if(n_b==0) vetoEvent;
Modified: trunk/src/Analyses/Makefile.am
==============================================================================
--- trunk/src/Analyses/Makefile.am Wed Nov 14 17:49:45 2012 (r4007)
+++ trunk/src/Analyses/Makefile.am Wed Nov 14 17:52:02 2012 (r4008)
@@ -95,6 +95,7 @@
ATLAS_2012_I1095236.cc \
ATLAS_2012_I1082009.cc \
ATLAS_2012_I946427.cc \
+ ATLAS_2012_I1126136.cc \
ATLAS_2012_I1180197.cc \
ATLAS_2012_I1186556.cc \
ATLAS_2012_I1190891.cc \
More information about the Rivet-svn
mailing list
|