|
|
[Rivet-svn] r4007 - in trunk: data/anainfo data/plotinfo src/Analyses
blackhole at projects.hepforge.org
blackhole at projects.hepforge.org
Wed Nov 14 17:49:45 GMT 2012
Author: buckley
Date: Wed Nov 14 17:49:45 2012
New Revision: 4007
Log:
Merge changeset [3970] and [3971] from AIDA branch and YODAifying
Added:
trunk/data/anainfo/ATLAS_2012_I1190891.info
- copied unchanged from r3970, branches/2012-06-aidarivet/data/anainfo/ATLAS_2012_I1190891.info
trunk/data/plotinfo/ATLAS_2012_I1190891.plot
- copied unchanged from r3971, branches/2012-06-aidarivet/data/plotinfo/ATLAS_2012_I1190891.plot
trunk/src/Analyses/ATLAS_2012_I1190891.cc
- copied, changed from r3970, branches/2012-06-aidarivet/src/Analyses/ATLAS_2012_I1190891.cc
Modified:
trunk/data/anainfo/Makefile.am
trunk/data/plotinfo/Makefile.am
trunk/src/Analyses/Makefile.am
Copied: trunk/data/anainfo/ATLAS_2012_I1190891.info (from r3970, branches/2012-06-aidarivet/data/anainfo/ATLAS_2012_I1190891.info)
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ trunk/data/anainfo/ATLAS_2012_I1190891.info Wed Nov 14 17:49:45 2012 (r4007, copy of r3970, branches/2012-06-aidarivet/data/anainfo/ATLAS_2012_I1190891.info)
@@ -0,0 +1,37 @@
+Name: ATLAS_2012_I1190891
+Year: 2012
+Summary: 4 or more lepton plus missing transverse energy SUSY search
+Experiment: ATLAS
+Collider: LHC
+InspireID: 1190891
+Status: UNVALIDATED
+Authors:
+ - Peter Richardson <peter.richardson at durham.ac.uk>
+References:
+ - ATLAS-CONF-2012-001
+ - ATLAS-CONF-2012-035
+ - arXiv:1210.4457
+RunInfo:
+ BSM signal events at 7000 GeV.
+NumEvents: 25000 for BSM signals
+Beams: [p+, p+]
+Energies: [7000]
+NeedCrossSection: yes
+Description:
+ 'Search for R-parity violating SUSY using events with 4 or more leptons in association with
+ missing transverse energy in proton-proton
+ collisions at a centre-of-mass energy of 7 TeV. The data sample has a total
+ integrated luminosity of 4.7 fb$^{-1}$.'
+BibTeX: '@article{Aad:2012kr,
+ author = "Aad, Georges and others",
+ title = "{Search for R-parity-violating supersymmetry in events
+ with four or more leptons in sqrt(s) = 7 TeV pp collisions
+ with the ATLAS detector}",
+ collaboration = "ATLAS Collaboration",
+ year = "2012",
+ eprint = "1210.4457",
+ archivePrefix = "arXiv",
+ primaryClass = "hep-ex",
+ reportNumber = "CERN-PH-EP-2012-276",
+ SLACcitation = "%%CITATION = ARXIV:1210.4457;%%",
+}'
Modified: trunk/data/anainfo/Makefile.am
==============================================================================
--- trunk/data/anainfo/Makefile.am Wed Nov 14 17:47:23 2012 (r4006)
+++ trunk/data/anainfo/Makefile.am Wed Nov 14 17:49:45 2012 (r4007)
@@ -61,6 +61,7 @@
ATLAS_2012_I1117704.info \
ATLAS_2012_I1180197.info \
ATLAS_2012_I1186556.info \
+ ATLAS_2012_I1190891.info \
ATLAS_2012_CONF_2012_103.info \
ATLAS_2012_CONF_2012_104.info \
ATLAS_2012_CONF_2012_105.info \
Copied: trunk/data/plotinfo/ATLAS_2012_I1190891.plot (from r3971, branches/2012-06-aidarivet/data/plotinfo/ATLAS_2012_I1190891.plot)
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ trunk/data/plotinfo/ATLAS_2012_I1190891.plot Wed Nov 14 17:49:45 2012 (r4007, copy of r3971, branches/2012-06-aidarivet/data/plotinfo/ATLAS_2012_I1190891.plot)
@@ -0,0 +1,24 @@
+# BEGIN PLOT /ATLAS_2012_I1190891/count_SR1
+Title=Number of events in the SR1 signal region
+XLabel=
+YLabel=Number of Events
+LogY=0
+# END PLOT
+# BEGIN PLOT /ATLAS_2012_I1190891/count_SR2
+Title=Number of events in the SR2 signal region
+XLabel=
+YLabel=Number of Events
+LogY=0
+# END PLOT
+# BEGIN PLOT /ATLAS_2012_I1190891/hist_etmiss
+Title=Missing Transverse Energy
+XLabel=$E^{\rm miss}_{\mathrm{T}}$ [GeV]
+YLabel=$\mathrm{Events}/50$\,GeV
+LogY=1
+# END PLOT
+# BEGIN PLOT /ATLAS_2012_I1190891/hist_m_eff
+Title=Effective Mass
+XLabel=$M_{\mathrm{eff}}$ [GeV]
+YLabel=$\mathrm{Events}/150$\,GeV
+LogY=1
+# END PLOT
Modified: trunk/data/plotinfo/Makefile.am
==============================================================================
--- trunk/data/plotinfo/Makefile.am Wed Nov 14 17:47:23 2012 (r4006)
+++ trunk/data/plotinfo/Makefile.am Wed Nov 14 17:49:45 2012 (r4007)
@@ -58,6 +58,7 @@
ATLAS_2012_I1082009.plot \
ATLAS_2012_I1183818.plot \
ATLAS_2012_I1186556.plot \
+ ATLAS_2012_I1190891.plot \
ATLAS_2012_CONF_2012_103.plot \
ATLAS_2012_CONF_2012_104.plot \
ATLAS_2012_CONF_2012_105.plot \
Copied and modified: trunk/src/Analyses/ATLAS_2012_I1190891.cc (from r3970, branches/2012-06-aidarivet/src/Analyses/ATLAS_2012_I1190891.cc)
==============================================================================
--- branches/2012-06-aidarivet/src/Analyses/ATLAS_2012_I1190891.cc Mon Oct 22 12:33:05 2012 (r3970, copy source)
+++ trunk/src/Analyses/ATLAS_2012_I1190891.cc Wed Nov 14 17:49:45 2012 (r4007)
@@ -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"
@@ -37,14 +37,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, 10.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);
@@ -63,27 +63,11 @@
// all tracks (to do deltaR with leptons)
addProjection(ChargedFinalState(-3.0,3.0),"cfs");
- // // Book histograms
- // _hist_leptonpT.push_back(bookHistogram1D(1,1,1));
- // _hist_leptonpT.push_back(bookHistogram1D(2,1,1));
- // _hist_leptonpT.push_back(bookHistogram1D(3,1,1));
- // _hist_leptonpT.push_back(bookHistogram1D(4,1,1));
- // _hist_njet = bookHistogram1D(5,1,1);
- // _hist_etmiss = bookHistogram1D(6,1,1);
- // _hist_mSFOS = bookHistogram1D(7,1,1);
- // _hist_meff = bookHistogram1D(8,1,1);
-
- // _hist_leptonpT_MC.push_back(bookHistogram1D("hist_lepton_pT_1", 26, 0., 260));
- // _hist_leptonpT_MC.push_back(bookHistogram1D("hist_lepton_pT_2", 15, 0., 150));
- // _hist_leptonpT_MC.push_back(bookHistogram1D("hist_lepton_pT_3", 20, 0., 100));
- // _hist_leptonpT_MC.push_back(bookHistogram1D("hist_lepton_pT_4", 20, 0., 100));
- // _hist_njet_MC = bookHistogram1D("hist_njet", 7, -0.5, 6.5);
- // _hist_etmiss_MC = bookHistogram1D("hist_etmiss",11,0.,220.);
- // _hist_mSFOS_MC = bookHistogram1D("hist_m_SFOS",13,0.,260.);
- // _hist_meff_MC = bookHistogram1D("hist_m_eff",19,0.,950.);
-
- // _count_SR1 = bookHistogram1D("count_SR1", 1, 0., 1.);
- // _count_SR2 = bookHistogram1D("count_SR2", 1, 0., 1.);
+ // Book histograms
+ _hist_etmiss = bookHisto1D("hist_etmiss",10,0.,500.);
+ _hist_meff = bookHisto1D("hist_m_eff",7,0.,1050.);
+ _count_SR1 = bookHisto1D("count_SR1", 1, 0., 1.);
+ _count_SR2 = bookHisto1D("count_SR2", 1, 0., 1.);
}
@@ -166,30 +150,30 @@
// and 0.1 from electrons ( keep higher energy)
for(unsigned int ie2=0;ie2<cand_e.size();++ie2) {
if(ie==ie2) continue;
- if ( deltaR(e.momentum(),cand_e[ie2].momentum()) < 0.1 &&
- e.momentum().E() < cand_e[ie2].momentum().E() ) {
- away = false;
- break;
- }
+ if ( deltaR(e.momentum(),cand_e[ie2].momentum()) < 0.1 &&
+ e.momentum().E() < cand_e[ie2].momentum().E() ) {
+ away = false;
+ break;
+ }
}
// if isolated keep it
if ( away )
- cand2_e.push_back( e );
+ cand2_e.push_back( e );
}
// remove e+e- pairs with mass < 20.
ParticleVector recon_e;
for(unsigned int ie=0;ie<cand2_e.size();++ie) {
- bool pass = true;
- for(unsigned int ie2=0;ie2<cand2_e.size();++ie2) {
- if(cand2_e[ie].pdgId()*cand2_e[ie2].pdgId()>0) continue;
- double mtest = (cand2_e[ie].momentum()+cand2_e[ie2].momentum()).mass();
- if(mtest<=20.) {
- pass = false;
- break;
- }
- }
- if(pass) recon_e.push_back(cand2_e[ie]);
+ bool pass = true;
+ for(unsigned int ie2=0;ie2<cand2_e.size();++ie2) {
+ if(cand2_e[ie].pdgId()*cand2_e[ie2].pdgId()>0) continue;
+ double mtest = (cand2_e[ie].momentum()+cand2_e[ie2].momentum()).mass();
+ if(mtest<=20.) {
+ pass = false;
+ break;
+ }
+ }
+ if(pass) recon_e.push_back(cand2_e[ie]);
}
// only keep muons more than R=0.4 from jets
@@ -220,16 +204,16 @@
// remove mu+mu- pairs with mass < 20.
ParticleVector recon_mu;
for(unsigned int imu=0;imu<cand2_mu.size();++imu) {
- bool pass = true;
- for(unsigned int imu2=0;imu2<cand2_mu.size();++imu2) {
- if(cand2_mu[imu].pdgId()*cand2_mu[imu2].pdgId()>0) continue;
- double mtest = (cand2_mu[imu].momentum()+cand2_mu[imu2].momentum()).mass();
- if(mtest<=20.) {
- pass = false;
- break;
- }
- }
- if(pass) recon_mu.push_back(cand2_mu[imu]);
+ bool pass = true;
+ for(unsigned int imu2=0;imu2<cand2_mu.size();++imu2) {
+ if(cand2_mu[imu].pdgId()*cand2_mu[imu2].pdgId()>0) continue;
+ double mtest = (cand2_mu[imu].momentum()+cand2_mu[imu2].momentum()).mass();
+ if(mtest<=20.) {
+ pass = false;
+ break;
+ }
+ }
+ if(pass) recon_mu.push_back(cand2_mu[imu]);
}
// pTmiss
@@ -250,16 +234,16 @@
}
// check if passes single lepton trigger
- bool passSingle =
- ( !recon_e .empty() && recon_e[0] .momentum().perp()>25. )||
- ( !recon_mu.empty() && recon_mu[0].momentum().perp()>20.);
+ bool passSingle =
+ ( !recon_e .empty() && recon_e[0] .momentum().perp()>25. )||
+ ( !recon_mu.empty() && recon_mu[0].momentum().perp()>20.);
// or two lepton trigger
- bool passDouble =
- ( recon_mu.size()>=2 && recon_mu[1].momentum().perp()>12.) ||
- ( recon_e .size()>=2 && recon_e [1].momentum().perp()>17.) ||
- ( !recon_e.empty() && !recon_mu.empty() &&
- recon_e[0].momentum().perp()>15. && recon_mu[0].momentum().perp()>10.);
+ bool passDouble =
+ ( recon_mu.size()>=2 && recon_mu[1].momentum().perp()>12.) ||
+ ( recon_e .size()>=2 && recon_e [1].momentum().perp()>17.) ||
+ ( !recon_e.empty() && !recon_mu.empty() &&
+ recon_e[0].momentum().perp()>15. && recon_mu[0].momentum().perp()>10.);
// must pass a trigger
if( !passSingle && !passDouble ) {
@@ -283,15 +267,15 @@
for(unsigned int iy=ix+1;iy<recon_e.size();++iy) {
if(recon_e[ix].pdgId()*recon_e[iy].pdgId()>0) continue;
double mtest = (recon_e[ix].momentum()+recon_e[iy].momentum()).mass();
- if(mtest>81.2 && mtest<101.2) vetoEvent;
+ if(mtest>81.2 && mtest<101.2) vetoEvent;
}
}
for(unsigned int ix=0;ix<recon_mu.size();++ix) {
for(unsigned int iy=ix+1;iy<recon_mu.size();++iy) {
if(recon_mu[ix].pdgId()*recon_mu[iy].pdgId()>0) continue;
- double mtest = (recon_mu[ix].momentum()+recon_mu[iy].momentum()).mass();
- if(mtest>81.2 && mtest<101.2) vetoEvent;
- }
+ double mtest = (recon_mu[ix].momentum()+recon_mu[iy].momentum()).mass();
+ if(mtest>81.2 && mtest<101.2) vetoEvent;
+ }
}
// make the control plots
@@ -316,10 +300,10 @@
/// @name Histograms
//@{
- AIDA::IHistogram1D* _hist_etmiss;
- AIDA::IHistogram1D* _hist_meff;
- AIDA::IHistogram1D* _count_SR1;
- AIDA::IHistogram1D* _count_SR2;
+ Histo1DPtr _hist_etmiss;
+ Histo1DPtr _hist_meff;
+ Histo1DPtr _count_SR1;
+ Histo1DPtr _count_SR2;
//@}
};
Modified: trunk/src/Analyses/Makefile.am
==============================================================================
--- trunk/src/Analyses/Makefile.am Wed Nov 14 17:47:23 2012 (r4006)
+++ trunk/src/Analyses/Makefile.am Wed Nov 14 17:49:45 2012 (r4007)
@@ -97,6 +97,7 @@
ATLAS_2012_I946427.cc \
ATLAS_2012_I1180197.cc \
ATLAS_2012_I1186556.cc \
+ ATLAS_2012_I1190891.cc \
ATLAS_2012_CONF_2012_103.cc \
ATLAS_2012_CONF_2012_104.cc \
ATLAS_2012_CONF_2012_105.cc \
More information about the Rivet-svn
mailing list
|