|
|
[Rivet-svn] r4048 - in trunk: data/anainfo data/plotinfo src/Analyses
blackhole at projects.hepforge.org
blackhole at projects.hepforge.org
Tue Nov 27 12:01:47 GMT 2012
Author: hoeth
Date: Tue Nov 27 12:01:41 2012
New Revision: 4048
Log:
merge r4039-4040 from branches/2012-06-aidarivet
Added:
trunk/data/anainfo/ATLAS_2012_CONF_2012_153.info
- copied unchanged from r4039, branches/2012-06-aidarivet/data/anainfo/ATLAS_2012_CONF_2012_153.info
trunk/data/plotinfo/ATLAS_2012_CONF_2012_153.plot
- copied unchanged from r4039, branches/2012-06-aidarivet/data/plotinfo/ATLAS_2012_CONF_2012_153.plot
trunk/src/Analyses/ATLAS_2012_CONF_2012_153.cc
- copied, changed from r4039, branches/2012-06-aidarivet/src/Analyses/ATLAS_2012_CONF_2012_153.cc
Modified:
trunk/data/anainfo/Makefile.am
trunk/data/plotinfo/Makefile.am
trunk/src/Analyses/Makefile.am
Copied: trunk/data/anainfo/ATLAS_2012_CONF_2012_153.info (from r4039, branches/2012-06-aidarivet/data/anainfo/ATLAS_2012_CONF_2012_153.info)
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ trunk/data/anainfo/ATLAS_2012_CONF_2012_153.info Tue Nov 27 12:01:41 2012 (r4048, copy of r4039, branches/2012-06-aidarivet/data/anainfo/ATLAS_2012_CONF_2012_153.info)
@@ -0,0 +1,25 @@
+Name: ATLAS_2012_CONF_2012_153
+Year: 2012
+Summary: 4 or more lepton plus missing transverse energy SUSY search
+Experiment: ATLAS
+Collider: LHC
+Status: PRELIMINARY
+Authors:
+ - Peter Richardson <peter.richardson at durham.ac.uk>
+References:
+ - ATLAS-CONF-2012-153
+RunInfo:
+ BSM signal events at 8000 GeV.
+NumEvents: 25000 for BSM signals
+Beams: [p+, p+]
+Energies: [8000]
+NeedCrossSection: yes
+Description:
+ 'Search for 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 8 TeV. The data sample has a total
+ integrated luminosity of 13.0 fb$^{-1}$. There is no reference data and
+ in addition to the control plots from the paper the number of events in
+ the two signal regions,
+ correctly normalized to an integrated luminosity 13.0 fb$^{-1}$,
+ are calculated.'
Modified: trunk/data/anainfo/Makefile.am
==============================================================================
--- trunk/data/anainfo/Makefile.am Tue Nov 27 11:54:33 2012 (r4047)
+++ trunk/data/anainfo/Makefile.am Tue Nov 27 12:01:41 2012 (r4048)
@@ -68,6 +68,7 @@
ATLAS_2012_CONF_2012_104.info \
ATLAS_2012_CONF_2012_105.info \
ATLAS_2012_CONF_2012_109.info \
+ ATLAS_2012_CONF_2012_153.info \
BABAR_2005_S6181155.info \
BABAR_2007_S6895344.info \
BABAR_2007_S7266081.info \
Copied: trunk/data/plotinfo/ATLAS_2012_CONF_2012_153.plot (from r4039, branches/2012-06-aidarivet/data/plotinfo/ATLAS_2012_CONF_2012_153.plot)
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ trunk/data/plotinfo/ATLAS_2012_CONF_2012_153.plot Tue Nov 27 12:01:41 2012 (r4048, copy of r4039, branches/2012-06-aidarivet/data/plotinfo/ATLAS_2012_CONF_2012_153.plot)
@@ -0,0 +1,24 @@
+# BEGIN PLOT /ATLAS_2012_CONF_2012_153/count_SR1
+Title=Number of events in the SR1 signal region
+XLabel=
+YLabel=Number of Events
+LogY=0
+# END PLOT
+# BEGIN PLOT /ATLAS_2012_CONF_2012_153/count_SR2
+Title=Number of events in the SR2 signal region
+XLabel=
+YLabel=Number of Events
+LogY=0
+# END PLOT
+# BEGIN PLOT /ATLAS_2012_CONF_2012_153/hist_etmiss
+Title=Missing Transverse Energy
+XLabel=$E^{\rm miss}_{\mathrm{T}}$ [GeV]
+YLabel=$\mathrm{Events}/20$\,GeV
+LogY=1
+# END PLOT
+# BEGIN PLOT /ATLAS_2012_CONF_2012_153/hist_m_eff
+Title=Effective Mass
+XLabel=$M_{\mathrm{eff}}$ [GeV]
+YLabel=$\mathrm{Events}/20$\,GeV
+LogY=1
+# END PLOT
Modified: trunk/data/plotinfo/Makefile.am
==============================================================================
--- trunk/data/plotinfo/Makefile.am Tue Nov 27 11:54:33 2012 (r4047)
+++ trunk/data/plotinfo/Makefile.am Tue Nov 27 12:01:41 2012 (r4048)
@@ -65,6 +65,7 @@
ATLAS_2012_CONF_2012_104.plot \
ATLAS_2012_CONF_2012_105.plot \
ATLAS_2012_CONF_2012_109.plot \
+ ATLAS_2012_CONF_2012_153.plot \
BABAR_2005_S6181155.plot \
BABAR_2007_S6895344.plot \
BABAR_2007_S7266081.plot \
Copied and modified: trunk/src/Analyses/ATLAS_2012_CONF_2012_153.cc (from r4039, branches/2012-06-aidarivet/src/Analyses/ATLAS_2012_CONF_2012_153.cc)
==============================================================================
--- branches/2012-06-aidarivet/src/Analyses/ATLAS_2012_CONF_2012_153.cc Fri Nov 23 21:23:17 2012 (r4039, copy source)
+++ trunk/src/Analyses/ATLAS_2012_CONF_2012_153.cc Tue Nov 27 12:01:41 2012 (r4048)
@@ -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"
@@ -81,17 +81,17 @@
edges_eT.push_back(150);
edges_eT.push_back(300);
edges_eT.push_back(500);
-
+
// Book histograms
- _hist_electrons = bookHistogram1D("hist_electrons_before", 11, -0.5,10.5);
- _hist_muons = bookHistogram1D("hist_muons_before" , 11, -0.5,10.5);
- _hist_leptons = bookHistogram1D("hist_leptons_before" , 11, -0.5,10.5);
- _hist_4leptons = bookHistogram1D("hist_4leptons", 1, 0.,1.);
- _hist_veto = bookHistogram1D("hist_veto", 1, 0., 1.);
- _hist_etmiss = bookHistogram1D("hist_etmiss",edges_eT);
- _hist_meff = bookHistogram1D("hist_m_eff",edges_meff);
- _count_SR1 = bookHistogram1D("count_SR1", 1, 0., 1.);
- _count_SR2 = bookHistogram1D("count_SR2", 1, 0., 1.);
+ _hist_electrons = bookHisto1D("hist_electrons_before", 11, -0.5,10.5);
+ _hist_muons = bookHisto1D("hist_muons_before" , 11, -0.5,10.5);
+ _hist_leptons = bookHisto1D("hist_leptons_before" , 11, -0.5,10.5);
+ _hist_4leptons = bookHisto1D("hist_4leptons", 1, 0.,1.);
+ _hist_veto = bookHisto1D("hist_veto", 1, 0., 1.);
+ _hist_etmiss = bookHisto1D("hist_etmiss",edges_eT);
+ _hist_meff = bookHisto1D("hist_m_eff",edges_meff);
+ _count_SR1 = bookHisto1D("count_SR1", 1, 0., 1.);
+ _count_SR2 = bookHisto1D("count_SR2", 1, 0., 1.);
}
@@ -117,13 +117,13 @@
vector<bool> vetoed(temp.size(),false);
ParticleVector cand_e;
for (unsigned int ix=0; ix<temp.size(); ++ix) {
- if(vetoed[ix]) continue;
- for (unsigned int iy=ix+1; iy<temp.size(); ++iy) {
- if( deltaR(temp[ix].momentum(),temp[iy].momentum()) < 0.1 ) {
- vetoed[iy] = true;
+ if(vetoed[ix]) continue;
+ for (unsigned int iy=ix+1; iy<temp.size(); ++iy) {
+ if( deltaR(temp[ix].momentum(),temp[iy].momentum()) < 0.1 ) {
+ vetoed[iy] = true;
}
- }
- if(!vetoed[ix]) cand_e.push_back(temp[ix]);
+ }
+ if(!vetoed[ix]) cand_e.push_back(temp[ix]);
}
// Sort by transverse momentum
@@ -156,7 +156,7 @@
}
// if isolated keep it
if ( away )
- cand2_e.push_back( e );
+ cand2_e.push_back( e );
}
// only keep muons more than R=0.4 from jets
@@ -177,25 +177,25 @@
// electron and muon more than 0.1 apart
ParticleVector cand3_e;
foreach ( const Particle & e, cand2_e ) {
- bool away = true;
- foreach( const Particle & mu, cand2_mu ) {
- if( deltaR(e.momentum(),mu.momentum()) < 0.1) {
- away = false;
- break;
- }
- }
- if(away) cand3_e.push_back(e);
+ bool away = true;
+ foreach( const Particle & mu, cand2_mu ) {
+ if( deltaR(e.momentum(),mu.momentum()) < 0.1) {
+ away = false;
+ break;
+ }
+ }
+ if(away) cand3_e.push_back(e);
}
ParticleVector cand3_mu;
foreach( const Particle & mu, cand2_mu ) {
- bool away = true;
- foreach ( const Particle & e, cand2_e ) {
- if( deltaR(e.momentum(),mu.momentum()) < 0.1) {
- away = false;
- break;
- }
- }
- if(away) cand3_mu.push_back(mu);
+ bool away = true;
+ foreach ( const Particle & e, cand2_e ) {
+ if( deltaR(e.momentum(),mu.momentum()) < 0.1) {
+ away = false;
+ break;
+ }
+ }
+ if(away) cand3_mu.push_back(mu);
}
// pTmiss
@@ -206,30 +206,30 @@
pTmiss -= p.momentum();
}
double eTmiss = pTmiss.pT();
-
+
// apply electron isolation
ParticleVector chg_tracks =
applyProjection<ChargedFinalState>(event, "cfs").particles();
ParticleVector cand4_e;
foreach ( const Particle & e, cand3_e ) {
- // charge isolation
+ // charge isolation
double pTinCone = -e.momentum().perp();
foreach ( const Particle & track, chg_tracks ) {
- if(track.momentum().perp()>0.4 &&
- deltaR(e.momentum(),track.momentum()) <= 0.3 )
+ if(track.momentum().perp()>0.4 &&
+ deltaR(e.momentum(),track.momentum()) <= 0.3 )
pTinCone += track.momentum().pT();
}
- if (pTinCone/e.momentum().perp()>0.16) continue;
- // all particles isolation
+ if (pTinCone/e.momentum().perp()>0.16) continue;
+ // all particles isolation
pTinCone = -e.momentum().perp();
- foreach ( const Particle & p, vfs_particles ) {
- if(abs(p.pdgId())!=MUON &&
- deltaR(e.momentum(),p.momentum()) <= 0.3 )
+ foreach ( const Particle & p, vfs_particles ) {
+ if(abs(p.pdgId())!=MUON &&
+ deltaR(e.momentum(),p.momentum()) <= 0.3 )
pTinCone += p.momentum().pT();
- }
+ }
if (pTinCone/e.momentum().perp()<0.18) {
cand4_e.push_back(e);
- }
+ }
}
// apply muon isolation
@@ -237,37 +237,37 @@
foreach ( const Particle & mu, cand3_mu ) {
double pTinCone = -mu.momentum().perp();
foreach ( const Particle & track, chg_tracks ) {
- if(track.momentum().perp()>1.0 &&
- deltaR(mu.momentum(),track.momentum()) <= 0.3 )
+ if(track.momentum().perp()>1.0 &&
+ deltaR(mu.momentum(),track.momentum()) <= 0.3 )
pTinCone += track.momentum().pT();
}
if (pTinCone/mu.momentum().perp()<0.12) {
cand4_mu.push_back(mu);
- }
+ }
}
// same SOSF pairs m>12.
ParticleVector recon_e;
foreach(const Particle & e, cand4_e) {
- bool veto=false;
- foreach(const Particle & e2, cand4_e) {
- if(e.pdgId()*e2.pdgId()<0&&(e.momentum()+e2.momentum()).mass()<12.) {
- veto=true;
- break;
- }
- }
- if(!veto) recon_e.push_back(e);
+ bool veto=false;
+ foreach(const Particle & e2, cand4_e) {
+ if(e.pdgId()*e2.pdgId()<0&&(e.momentum()+e2.momentum()).mass()<12.) {
+ veto=true;
+ break;
+ }
+ }
+ if(!veto) recon_e.push_back(e);
}
ParticleVector recon_mu;
foreach(const Particle & mu, cand4_mu) {
- bool veto=false;
- foreach(const Particle & mu2, cand4_mu) {
- if(mu.pdgId()*mu2.pdgId()<0&&(mu.momentum()+mu2.momentum()).mass()<12.) {
- veto=true;
- break;
- }
- }
- if(!veto) recon_mu.push_back(mu);
+ bool veto=false;
+ foreach(const Particle & mu2, cand4_mu) {
+ if(mu.pdgId()*mu2.pdgId()<0&&(mu.momentum()+mu2.momentum()).mass()<12.) {
+ veto=true;
+ break;
+ }
+ }
+ if(!veto) recon_mu.push_back(mu);
}
// now only use recon_jets, recon_mu, recon_e
@@ -275,7 +275,7 @@
_hist_muons->fill(recon_mu.size(), weight);
_hist_leptons->fill(recon_mu.size() + recon_e.size(), weight);
if( recon_mu.size() + recon_e.size() > 3) {
- _hist_4leptons->fill(0.5, weight);
+ _hist_4leptons->fill(0.5, weight);
}
// reject events with less than 4 electrons and muons
@@ -286,14 +286,14 @@
// or two lepton trigger
- bool passDouble =
- (recon_mu.size()>=2 && ( (recon_mu[1].momentum().perp()>14.) ||
- (recon_mu[0].momentum().perp()>18. && recon_mu[1].momentum().perp()>10.) )) ||
- (recon_e.size() >=2 && ( (recon_e [1].momentum().perp()>14.) ||
- (recon_e [0].momentum().perp()>25. && recon_e [1].momentum().perp()>10.) )) ||
- (!recon_e.empty() && !recon_mu.empty() &&
- ( (recon_e[0].momentum().perp()>14. && recon_mu[0].momentum().perp()>10.)||
- (recon_e[0].momentum().perp()>10. && recon_mu[0].momentum().perp()>18.) ));
+ bool passDouble =
+ (recon_mu.size()>=2 && ( (recon_mu[1].momentum().perp()>14.) ||
+ (recon_mu[0].momentum().perp()>18. && recon_mu[1].momentum().perp()>10.) )) ||
+ (recon_e.size() >=2 && ( (recon_e [1].momentum().perp()>14.) ||
+ (recon_e [0].momentum().perp()>25. && recon_e [1].momentum().perp()>10.) )) ||
+ (!recon_e.empty() && !recon_mu.empty() &&
+ ( (recon_e[0].momentum().perp()>14. && recon_mu[0].momentum().perp()>10.)||
+ (recon_e[0].momentum().perp()>10. && recon_mu[0].momentum().perp()>18.) ));
// must pass a trigger
if(!passDouble ) {
@@ -322,14 +322,14 @@
double mtest = ppair.mass();
if(mtest>81.2 && mtest<101.2) vetoEvent;
- // check triplets with electron
+ // check triplets with electron
for(unsigned int iz=0;iz<recon_e.size();++iz) {
if(iz==ix||iz==iy) continue;
mtest = (ppair+recon_e[iz].momentum()).mass();
if(mtest>81.2 && mtest<101.2) vetoEvent;
- }
+ }
- // check triplets with muon
+ // check triplets with muon
for(unsigned int iz=0;iz<recon_mu.size();++iz) {
mtest = (ppair+recon_mu[iz].momentum()).mass();
if(mtest>81.2 && mtest<101.2) vetoEvent;
@@ -362,39 +362,39 @@
FourMomentum ppair = recon_mu[ix].momentum()+recon_mu[iy].momentum();
double mtest = ppair.mass();
if(mtest>81.2 && mtest<101.2) vetoEvent;
-
- // check triplets with muon
+
+ // check triplets with muon
for(unsigned int iz=0;iz<recon_mu.size();++iz) {
if(iz==ix||iz==iy) continue;
mtest = (ppair+recon_mu[iz].momentum()).mass();
if(mtest>81.2 && mtest<101.2) vetoEvent;
- }
+ }
- // check triplets with electron
+ // check triplets with electron
for(unsigned int iz=0;iz<recon_e.size();++iz) {
mtest = (ppair+recon_e[iz].momentum()).mass();
if(mtest>81.2 && mtest<101.2) vetoEvent;
- }
+ }
- // check muon quadruplets
+ // check muon quadruplets
for(unsigned int iz=0;iz<recon_mu.size();++iz) {
for(unsigned int iw=iz+1;iy<recon_mu.size();++iy) {
- if(iz==ix||iz==iy||iw==ix||iw==iy) continue;
+ if(iz==ix||iz==iy||iw==ix||iw==iy) continue;
if(recon_mu[iz].pdgId()*recon_mu[iw].pdgId()>0) continue;
mtest = (ppair+recon_mu[iz].momentum()+recon_mu[iw].momentum()).mass();
if(mtest>81.2 && mtest<101.2) vetoEvent;
- }
- }
+ }
+ }
}
}
-
+
//make the control plots
_hist_etmiss ->fill(eTmiss,weight);
_hist_meff ->fill(meff ,weight);
// finally the counts
if(eTmiss>50.) _count_SR1->fill(0.5,weight);
if(meff >0. ) _count_SR2->fill(0.5,weight);
-
+
}
//@}
@@ -411,15 +411,15 @@
/// @name Histograms
//@{
- AIDA::IHistogram1D* _hist_electrons;
- AIDA::IHistogram1D* _hist_muons;
- AIDA::IHistogram1D* _hist_leptons;
- AIDA::IHistogram1D* _hist_4leptons;
- AIDA::IHistogram1D* _hist_veto;
- AIDA::IHistogram1D* _hist_etmiss;
- AIDA::IHistogram1D* _hist_meff;
- AIDA::IHistogram1D* _count_SR1;
- AIDA::IHistogram1D* _count_SR2;
+ Histo1DPtr _hist_electrons;
+ Histo1DPtr _hist_muons;
+ Histo1DPtr _hist_leptons;
+ Histo1DPtr _hist_4leptons;
+ Histo1DPtr _hist_veto;
+ Histo1DPtr _hist_etmiss;
+ Histo1DPtr _hist_meff;
+ Histo1DPtr _count_SR1;
+ Histo1DPtr _count_SR2;
//@}
};
Modified: trunk/src/Analyses/Makefile.am
==============================================================================
--- trunk/src/Analyses/Makefile.am Tue Nov 27 11:54:33 2012 (r4047)
+++ trunk/src/Analyses/Makefile.am Tue Nov 27 12:01:41 2012 (r4048)
@@ -103,7 +103,8 @@
ATLAS_2012_CONF_2012_103.cc \
ATLAS_2012_CONF_2012_104.cc \
ATLAS_2012_CONF_2012_105.cc \
- ATLAS_2012_CONF_2012_109.cc
+ ATLAS_2012_CONF_2012_109.cc \
+ ATLAS_2012_CONF_2012_153.cc
endif
More information about the Rivet-svn
mailing list
|