|
|
[Rivet-svn] r3755 - in branches/2011-07-aida2yoda: data/anainfo data/plotinfo src/Analyses
blackhole at projects.hepforge.org
blackhole at projects.hepforge.org
Fri Jun 8 15:23:42 BST 2012
Author: hoeth
Date: Fri Jun 8 15:23:42 2012
New Revision: 3755
Log:
merge c3718, 3743 from trunk
Added:
branches/2011-07-aida2yoda/data/anainfo/ATLAS_2012_I1112263.info
- copied unchanged from r3718, trunk/data/anainfo/ATLAS_2012_I1112263.info
branches/2011-07-aida2yoda/data/plotinfo/ATLAS_2012_I1112263.plot
- copied unchanged from r3718, trunk/data/plotinfo/ATLAS_2012_I1112263.plot
branches/2011-07-aida2yoda/src/Analyses/ATLAS_2012_I1112263.cc
- copied, changed from r3718, trunk/src/Analyses/ATLAS_2012_I1112263.cc
Deleted:
branches/2011-07-aida2yoda/data/anainfo/ATLAS_2012_CONF_2012_023.info
branches/2011-07-aida2yoda/data/plotinfo/ATLAS_2012_CONF_2012_023.plot
branches/2011-07-aida2yoda/src/Analyses/ATLAS_2012_CONF_2012_023.cc
Modified:
branches/2011-07-aida2yoda/data/anainfo/Makefile.am
branches/2011-07-aida2yoda/data/plotinfo/Makefile.am
branches/2011-07-aida2yoda/src/Analyses/Makefile.am
Copied: branches/2011-07-aida2yoda/data/anainfo/ATLAS_2012_I1112263.info (from r3718, trunk/data/anainfo/ATLAS_2012_I1112263.info)
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ branches/2011-07-aida2yoda/data/anainfo/ATLAS_2012_I1112263.info Fri Jun 8 15:23:42 2012 (r3755, copy of r3718, trunk/data/anainfo/ATLAS_2012_I1112263.info)
@@ -0,0 +1,42 @@
+Name: ATLAS_2012_I1112263
+Year: 2012
+Summary: 3 lepton plus missing transverse energy SUSY search
+Experiment: ATLAS
+Collider: LHC
+Status: VALIDATED
+Authors:
+ - Peter Richardson <peter.richardson at durham.ac.uk>
+References:
+ - ATLAS-CONF-2012-023
+ - arXiv:1204.5638
+RunInfo:
+ BSM signal events at 7000 GeV.
+NumEvents: 25000 for BSM signals
+Beams: [p+, p+]
+Energies: [7000]
+NeedCrossSection: yes
+Description:
+ 'Search for SUSY using events with 3 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 2.06 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 2.06 fb$^{-1}$,
+ are calculated.'
+BibKey: Aad:2012I1112263
+BibTeX: '@article{Aad:2012I1112263,
+ author = "Aad, Georges and others",
+ title = "{Search for supersymmetry in events with three leptons
+ and missing transverse momentum in $\sqrt{s}$ = 7 TeV pp
+ collisions with the ATLAS detector}",
+ collaboration = "ATLAS Collaboration",
+ year = "2012",
+ note = "5 pages plus author list (18 pages total), 2 figures, 1
+ table, submitted to Physics Review Letters",
+ eprint = "1204.5638",
+ archivePrefix = "arXiv",
+ primaryClass = "hep-ex",
+ reportNumber = "CERN-PH-EP-2012-098",
+ SLACcitation = "%%CITATION = ARXIV:1204.5638;%%",
+}'
Modified: branches/2011-07-aida2yoda/data/anainfo/Makefile.am
==============================================================================
--- branches/2011-07-aida2yoda/data/anainfo/Makefile.am Fri Jun 8 15:17:44 2012 (r3754)
+++ branches/2011-07-aida2yoda/data/anainfo/Makefile.am Fri Jun 8 15:23:42 2012 (r3755)
@@ -48,7 +48,7 @@
ATLAS_2012_I1094568.info \
ATLAS_2012_I1095236.info \
ATLAS_2012_CONF_2012_001.info \
- ATLAS_2012_CONF_2012_023.info \
+ ATLAS_2012_I1112263.info \
ATLAS_2012_CONF_2012_033.info \
ATLAS_2012_CONF_2012_037.info \
ATLAS_2012_CONF_2012_041.info \
Copied: branches/2011-07-aida2yoda/data/plotinfo/ATLAS_2012_I1112263.plot (from r3718, trunk/data/plotinfo/ATLAS_2012_I1112263.plot)
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ branches/2011-07-aida2yoda/data/plotinfo/ATLAS_2012_I1112263.plot Fri Jun 8 15:23:42 2012 (r3755, copy of r3718, trunk/data/plotinfo/ATLAS_2012_I1112263.plot)
@@ -0,0 +1,82 @@
+# BEGIN PLOT /ATLAS_2012_I1112263/count_SR1
+Title=Number of events in the SR1 signal region
+XLabel=
+YLabel=Number of Events
+LogY=0
+# END PLOT
+# BEGIN PLOT /ATLAS_2012_I1112263/count_SR2
+Title=Number of events in the SR2 signal region
+XLabel=
+YLabel=Number of Events
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2012_I1112263/hist_lepton_pT_1_SR1
+Title=$E^e_{\mathrm{T}}(p^\mu_{\mathrm{T}})$ for the leading lepton in SR1
+XLabel=$E^e_{\mathrm{T}}(p^\mu_{\mathrm{T}})$ [GeV]
+YLabel=$\mathrm{Events}/20 \mathrm{GeV}$
+LogY=1
+# END PLOT
+# BEGIN PLOT /ATLAS_2012_I1112263/hist_lepton_pT_2_SR1
+Title=$E^e_{\mathrm{T}}(p^\mu_{\mathrm{T}})$ for the $2^{\rm nd}$ leading lepton in SR1
+XLabel=$E^e_{\mathrm{T}}(p^\mu_{\mathrm{T}})$ [GeV]
+YLabel=$\mathrm{Events}/20 \mathrm{GeV}$
+LogY=1
+# END PLOT
+# BEGIN PLOT /ATLAS_2012_I1112263/hist_lepton_pT_3_SR1
+Title=$E^e_{\mathrm{T}}(p^\mu_{\mathrm{T}})$ for the $3^{\rm rd}$ leading lepton in SR1
+XLabel=$E^e_{\mathrm{T}}(p^\mu_{\mathrm{T}})$ [GeV]
+YLabel=$\mathrm{Events}/20 \mathrm{GeV}$
+LogY=1
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2012_I1112263/hist_lepton_pT_1_SR2
+Title=$E^e_{\mathrm{T}}(p^\mu_{\mathrm{T}})$ for the leading lepton in SR2
+XLabel=$E^e_{\mathrm{T}}(p^\mu_{\mathrm{T}})$ [GeV]
+YLabel=$\mathrm{Events}/20 \mathrm{GeV}$
+LogY=1
+# END PLOT
+# BEGIN PLOT /ATLAS_2012_I1112263/hist_lepton_pT_2_SR2
+Title=$E^e_{\mathrm{T}}(p^\mu_{\mathrm{T}})$ for the $2^{\rm nd}$ leading lepton in SR2
+XLabel=$E^e_{\mathrm{T}}(p^\mu_{\mathrm{T}})$ [GeV]
+YLabel=$\mathrm{Events}/20 \mathrm{GeV}$
+LogY=1
+# END PLOT
+# BEGIN PLOT /ATLAS_2012_I1112263/hist_lepton_pT_3_SR2
+Title=$E^e_{\mathrm{T}}(p^\mu_{\mathrm{T}})$ for the $3^{\rm rd}$ leading lepton in SR2
+XLabel=$E^e_{\mathrm{T}}(p^\mu_{\mathrm{T}})$ [GeV]
+YLabel=$\mathrm{Events}/20 \mathrm{GeV}$
+LogY=1
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2012_I1112263/hist_etmiss_SR1_A
+Title=Missing Transverse Energy in SR1
+XLabel=$E^{\rm miss}_{\mathrm{T}}$ [GeV]
+YLabel=$\mathrm{Events}/20$\,GeV
+LogY=0
+# END PLOT
+# BEGIN PLOT /ATLAS_2012_I1112263/hist_etmiss_SR1_B
+Title=Missing Transverse Energy in SR1
+XLabel=$E^{\rm miss}_{\mathrm{T}}$ [GeV]
+YLabel=$\mathrm{Events}/20$\,GeV
+LogY=1
+# END PLOT
+# BEGIN PLOT /ATLAS_2012_I1112263/hist_etmiss_SR2_A
+Title=Missing Transverse Energy in SR2
+XLabel=$E^{\rm miss}_{\mathrm{T}}$ [GeV]
+YLabel=$\mathrm{Events}/20$\,GeV
+LogY=0
+# END PLOT
+# BEGIN PLOT /ATLAS_2012_I1112263/hist_etmiss_SR2_B
+Title=Missing Transverse Energy in SR2
+XLabel=$E^{\rm miss}_{\mathrm{T}}$ [GeV]
+YLabel=$\mathrm{Events}/20$\,GeV
+LogY=1
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2012_I1112263/hist_mSFOF
+Title=Mass of SFOF lepton pair closest to the $Z^0$ mass in SR1
+XLabel=$M_{\mathrm{SFOS}}$ [GeV]
+YLabel=$\mathrm{Events}/20$\,GeV
+LogY=1
+# END PLOT
Modified: branches/2011-07-aida2yoda/data/plotinfo/Makefile.am
==============================================================================
--- branches/2011-07-aida2yoda/data/plotinfo/Makefile.am Fri Jun 8 15:17:44 2012 (r3754)
+++ branches/2011-07-aida2yoda/data/plotinfo/Makefile.am Fri Jun 8 15:23:42 2012 (r3755)
@@ -45,7 +45,7 @@
ATLAS_2012_I1094568.plot \
ATLAS_2012_I1095236.plot \
ATLAS_2012_CONF_2012_001.plot \
- ATLAS_2012_CONF_2012_023.plot \
+ ATLAS_2012_I1112263.plot \
ATLAS_2012_CONF_2012_033.plot \
ATLAS_2012_CONF_2012_037.plot \
ATLAS_2012_CONF_2012_041.plot \
Copied and modified: branches/2011-07-aida2yoda/src/Analyses/ATLAS_2012_I1112263.cc (from r3718, trunk/src/Analyses/ATLAS_2012_I1112263.cc)
==============================================================================
--- trunk/src/Analyses/ATLAS_2012_I1112263.cc Wed May 2 10:24:37 2012 (r3718, copy source)
+++ branches/2011-07-aida2yoda/src/Analyses/ATLAS_2012_I1112263.cc Fri Jun 8 15:23:42 2012 (r3755)
@@ -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"
@@ -64,20 +64,20 @@
addProjection(ChargedFinalState(-3.0,3.0),"cfs");
// Book histograms
- _hist_leptonpT_SR1.push_back(bookHistogram1D("hist_lepton_pT_1_SR1",11,0.,220.));
- _hist_leptonpT_SR1.push_back(bookHistogram1D("hist_lepton_pT_2_SR1", 7,0.,140.));
- _hist_leptonpT_SR1.push_back(bookHistogram1D("hist_lepton_pT_3_SR1", 8,0.,160.));
- _hist_leptonpT_SR2.push_back(bookHistogram1D("hist_lepton_pT_1_SR2",11,0.,220.));
- _hist_leptonpT_SR2.push_back(bookHistogram1D("hist_lepton_pT_2_SR2", 7,0.,140.));
- _hist_leptonpT_SR2.push_back(bookHistogram1D("hist_lepton_pT_3_SR2", 8,0.,160.));
- _hist_etmiss_SR1_A = bookHistogram1D("hist_etmiss_SR1_A",15,10.,310.);
- _hist_etmiss_SR1_B = bookHistogram1D("hist_etmiss_SR1_B", 9,10.,190.);
- _hist_etmiss_SR2_A = bookHistogram1D("hist_etmiss_SR2_A",15,10.,310.);
- _hist_etmiss_SR2_B = bookHistogram1D("hist_etmiss_SR2_B", 9,10.,190.);
- _hist_mSFOS= bookHistogram1D("hist_mSFOF",9,0.,180.);
+ _hist_leptonpT_SR1.push_back(bookHisto1D("hist_lepton_pT_1_SR1",11,0.,220.));
+ _hist_leptonpT_SR1.push_back(bookHisto1D("hist_lepton_pT_2_SR1", 7,0.,140.));
+ _hist_leptonpT_SR1.push_back(bookHisto1D("hist_lepton_pT_3_SR1", 8,0.,160.));
+ _hist_leptonpT_SR2.push_back(bookHisto1D("hist_lepton_pT_1_SR2",11,0.,220.));
+ _hist_leptonpT_SR2.push_back(bookHisto1D("hist_lepton_pT_2_SR2", 7,0.,140.));
+ _hist_leptonpT_SR2.push_back(bookHisto1D("hist_lepton_pT_3_SR2", 8,0.,160.));
+ _hist_etmiss_SR1_A = bookHisto1D("hist_etmiss_SR1_A",15,10.,310.);
+ _hist_etmiss_SR1_B = bookHisto1D("hist_etmiss_SR1_B", 9,10.,190.);
+ _hist_etmiss_SR2_A = bookHisto1D("hist_etmiss_SR2_A",15,10.,310.);
+ _hist_etmiss_SR2_B = bookHisto1D("hist_etmiss_SR2_B", 9,10.,190.);
+ _hist_mSFOS= bookHisto1D("hist_mSFOF",9,0.,180.);
- _count_SR1 = bookHistogram1D("count_SR1", 1, 0., 1.);
- _count_SR2 = bookHistogram1D("count_SR2", 1, 0., 1.);
+ _count_SR1 = bookHisto1D("count_SR1", 1, 0., 1.);
+ _count_SR2 = bookHisto1D("count_SR2", 1, 0., 1.);
}
@@ -111,21 +111,21 @@
// candidate electrons
ParticleVector cand_e;
foreach ( const Particle & e,
- applyProjection<IdentifiedFinalState>(event, "elecs").particlesByPt() ) {
- double eta = e.momentum().eta();
- // remove electrons with pT<15 in old veto region
- // (NOT EXPLICIT IN THIS PAPER BUT IN SIMILAR 4 LEPTON PAPER and THIS DESCRPITION
- // IS MUCH WORSE SO ASSUME THIS IS DONE)
- if( fabs(eta)>1.37 && fabs(eta) < 1.52 && e.momentum().perp()< 15.*GeV)
- continue;
+ applyProjection<IdentifiedFinalState>(event, "elecs").particlesByPt() ) {
+ double eta = e.momentum().eta();
+ // remove electrons with pT<15 in old veto region
+ // (NOT EXPLICIT IN THIS PAPER BUT IN SIMILAR 4 LEPTON PAPER and THIS DESCRPITION
+ // IS MUCH WORSE SO ASSUME THIS IS DONE)
+ if( fabs(eta)>1.37 && fabs(eta) < 1.52 && e.momentum().perp()< 15.*GeV)
+ continue;
double pTinCone = -e.momentum().perp();
foreach ( const Particle & track, chg_tracks ) {
if ( deltaR(e.momentum(),track.momentum()) <= 0.2 )
pTinCone += track.momentum().pT();
}
- if (pTinCone/e.momentum().perp()<0.1) {
+ if (pTinCone/e.momentum().perp()<0.1) {
cand_e.push_back(e);
- }
+ }
}
// resolve jet/lepton ambiguity
@@ -133,15 +133,15 @@
// IS MUCH WORSE SO ASSUME THIS IS DONE)
Jets recon_jets;
foreach ( const Jet& jet, cand_jets ) {
- bool away_from_e = true;
- foreach ( const Particle & e, cand_e ) {
- if ( deltaR(e.momentum(),jet.momentum()) <= 0.2 ) {
- away_from_e = false;
- break;
- }
- }
- if ( away_from_e )
- recon_jets.push_back( jet );
+ bool away_from_e = true;
+ foreach ( const Particle & e, cand_e ) {
+ if ( deltaR(e.momentum(),jet.momentum()) <= 0.2 ) {
+ away_from_e = false;
+ break;
+ }
+ }
+ if ( away_from_e )
+ recon_jets.push_back( jet );
}
// only keep electrons more than R=0.4 from jets
@@ -154,15 +154,15 @@
break;
}
}
- // and 0.1 from any muons
- if ( ! away ) {
- foreach ( const Particle & mu, cand_e ) {
- if ( deltaR(mu.momentum(),e.momentum()) < 0.1 ) {
- away = false;
- break;
- }
- }
- }
+ // and 0.1 from any muons
+ if ( ! away ) {
+ foreach ( const Particle & mu, cand_e ) {
+ if ( deltaR(mu.momentum(),e.momentum()) < 0.1 ) {
+ away = false;
+ break;
+ }
+ }
+ }
if ( away )
recon_e.push_back( e );
}
@@ -176,15 +176,15 @@
break;
}
}
- // and 0.1 from any electrona
- if ( ! away ) {
- foreach ( const Particle & e, cand_e ) {
- if ( deltaR(mu.momentum(),e.momentum()) < 0.1 ) {
- away = false;
- break;
- }
- }
- }
+ // and 0.1 from any electrona
+ if ( ! away ) {
+ foreach ( const Particle & e, cand_e ) {
+ if ( deltaR(mu.momentum(),e.momentum()) < 0.1 ) {
+ away = false;
+ break;
+ }
+ }
+ }
if ( away )
recon_mu.push_back( mu );
}
@@ -208,25 +208,25 @@
// ATLAS calo problem
if(rand()/static_cast<double>(RAND_MAX)<=0.42) {
- foreach ( const Particle & e, recon_e ) {
- double eta = e.momentum().eta();
- double phi = e.momentum().azimuthalAngle(MINUSPI_PLUSPI);
- if(eta>-0.1&&eta<1.5&&phi>-0.9&&phi<-0.5)
- vetoEvent;
- }
- foreach ( const Jet & jet, recon_jets ) {
- double eta = jet.momentum().rapidity();
- double phi = jet.momentum().azimuthalAngle(MINUSPI_PLUSPI);
- if(jet.momentum().perp()>40 && eta>-0.1&&eta<1.5&&phi>-0.9&&phi<-0.5)
- vetoEvent;
- }
+ foreach ( const Particle & e, recon_e ) {
+ double eta = e.momentum().eta();
+ double phi = e.momentum().azimuthalAngle(MINUSPI_PLUSPI);
+ if(eta>-0.1&&eta<1.5&&phi>-0.9&&phi<-0.5)
+ vetoEvent;
+ }
+ foreach ( const Jet & jet, recon_jets ) {
+ double eta = jet.momentum().rapidity();
+ double phi = jet.momentum().azimuthalAngle(MINUSPI_PLUSPI);
+ if(jet.momentum().perp()>40 && eta>-0.1&&eta<1.5&&phi>-0.9&&phi<-0.5)
+ vetoEvent;
+ }
}
// check at least one e/mu passing trigger
if( !( !recon_e .empty() && recon_e[0] .momentum().perp()>25.) &&
- !( !recon_mu.empty() && recon_mu[0].momentum().perp()>20.) ) {
- MSG_DEBUG("Hardest lepton fails trigger");
- vetoEvent;
+ !( !recon_mu.empty() && recon_mu[0].momentum().perp()>20.) ) {
+ MSG_DEBUG("Hardest lepton fails trigger");
+ vetoEvent;
}
// eTmiss cut
@@ -236,72 +236,72 @@
double mSFOS=1e30, mdiff=1e30;
unsigned int nSFOS=0;
for(unsigned int ix=0;ix<recon_e.size();++ix) {
- for(unsigned int iy=ix+1;iy<recon_e.size();++iy) {
- if(recon_e[ix].pdgId()*recon_e[iy].pdgId()>0) continue;
- ++nSFOS;
- double mtest = (recon_e[ix].momentum()+recon_e[iy].momentum()).mass();
- // veto is mass<20
- if(mtest<20.) vetoEvent;
- if(fabs(mtest-90.)<mdiff) {
- mSFOS = mtest;
- mdiff = fabs(mtest-90.);
- }
- }
+ for(unsigned int iy=ix+1;iy<recon_e.size();++iy) {
+ if(recon_e[ix].pdgId()*recon_e[iy].pdgId()>0) continue;
+ ++nSFOS;
+ double mtest = (recon_e[ix].momentum()+recon_e[iy].momentum()).mass();
+ // veto is mass<20
+ if(mtest<20.) vetoEvent;
+ if(fabs(mtest-90.)<mdiff) {
+ mSFOS = mtest;
+ mdiff = fabs(mtest-90.);
+ }
+ }
}
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;
- ++nSFOS;
- double mtest = (recon_mu[ix].momentum()+recon_mu[iy].momentum()).mass();
- // veto is mass<20
- if(mtest<20.) vetoEvent;
- if(fabs(mtest-90.)<mdiff) {
- mSFOS = mtest;
- mdiff = fabs(mtest-90.);
- }
- }
+ for(unsigned int iy=ix+1;iy<recon_mu.size();++iy) {
+ if(recon_mu[ix].pdgId()*recon_mu[iy].pdgId()>0) continue;
+ ++nSFOS;
+ double mtest = (recon_mu[ix].momentum()+recon_mu[iy].momentum()).mass();
+ // veto is mass<20
+ if(mtest<20.) vetoEvent;
+ if(fabs(mtest-90.)<mdiff) {
+ mSFOS = mtest;
+ mdiff = fabs(mtest-90.);
+ }
+ }
}
// require at least 1 SFOS pair
if(nSFOS==0) vetoEvent;
// b-jet veto in SR!
if(mdiff>10.) {
- foreach (const Jet & jet, recon_jets ) {
- if(jet.containsBottom() &&
- rand()/static_cast<double>(RAND_MAX)<=0.60)
- vetoEvent;
- }
+ foreach (const Jet & jet, recon_jets ) {
+ if(jet.containsBottom() &&
+ rand()/static_cast<double>(RAND_MAX)<=0.60)
+ vetoEvent;
+ }
}
// region SR1, Z depleted
if(mdiff>10.) {
- _count_SR1->fill(0.5,weight);
- _hist_etmiss_SR1_A->fill(eTmiss,weight);
- _hist_etmiss_SR1_B->fill(eTmiss,weight);
- _hist_mSFOS->fill(mSFOS,weight);
+ _count_SR1->fill(0.5,weight);
+ _hist_etmiss_SR1_A->fill(eTmiss,weight);
+ _hist_etmiss_SR1_B->fill(eTmiss,weight);
+ _hist_mSFOS->fill(mSFOS,weight);
}
// region SR2, Z enriched
else {
- _count_SR2->fill(0.5,weight);
- _hist_etmiss_SR2_A->fill(eTmiss,weight);
- _hist_etmiss_SR2_B->fill(eTmiss,weight);
+ _count_SR2->fill(0.5,weight);
+ _hist_etmiss_SR2_A->fill(eTmiss,weight);
+ _hist_etmiss_SR2_B->fill(eTmiss,weight);
}
// make the control plots
// lepton pT
unsigned int ie=0,imu=0;
for(unsigned int ix=0;ix<3;++ix) {
- AIDA::IHistogram1D* hist = mdiff>10. ?
- _hist_leptonpT_SR1[ix] : _hist_leptonpT_SR2[ix];
- double pTe = ie <recon_e .size() ?
- recon_e [ie ].momentum().perp() : -1*GeV;
- double pTmu = imu<recon_mu.size() ?
- recon_mu[imu].momentum().perp() : -1*GeV;
- if(pTe>pTmu) {
- hist->fill(pTe ,weight);
- ++ie;
- }
- else {
- hist->fill(pTmu,weight);
- ++imu;
- }
+ Histo1DPtr hist = mdiff>10. ?
+ _hist_leptonpT_SR1[ix] : _hist_leptonpT_SR2[ix];
+ double pTe = ie <recon_e .size() ?
+ recon_e [ie ].momentum().perp() : -1*GeV;
+ double pTmu = imu<recon_mu.size() ?
+ recon_mu[imu].momentum().perp() : -1*GeV;
+ if(pTe>pTmu) {
+ hist->fill(pTe ,weight);
+ ++ie;
+ }
+ else {
+ hist->fill(pTmu,weight);
+ ++imu;
+ }
}
}
@@ -311,8 +311,8 @@
double norm = crossSection()/femtobarn*2.06/sumOfWeights();
// these are number of events at 2.06fb^-1 per 20 GeV
for(unsigned int ix=0;ix<3;++ix) {
- scale(_hist_leptonpT_SR1[ix],norm*20.);
- scale(_hist_leptonpT_SR2[ix],norm*20.);
+ scale(_hist_leptonpT_SR1[ix],norm*20.);
+ scale(_hist_leptonpT_SR2[ix],norm*20.);
}
scale(_hist_etmiss_SR1_A,norm*20.);
scale(_hist_etmiss_SR1_B,norm*20.);
@@ -328,15 +328,15 @@
/// @name Histograms
//@{
- vector<AIDA::IHistogram1D*> _hist_leptonpT_SR1;
- vector<AIDA::IHistogram1D*> _hist_leptonpT_SR2;
- AIDA::IHistogram1D* _hist_etmiss_SR1_A;
- AIDA::IHistogram1D* _hist_etmiss_SR1_B;
- AIDA::IHistogram1D* _hist_etmiss_SR2_A;
- AIDA::IHistogram1D* _hist_etmiss_SR2_B;
- AIDA::IHistogram1D* _hist_mSFOS;
- AIDA::IHistogram1D* _count_SR1;
- AIDA::IHistogram1D* _count_SR2;
+ vector<Histo1DPtr> _hist_leptonpT_SR1;
+ vector<Histo1DPtr> _hist_leptonpT_SR2;
+ Histo1DPtr _hist_etmiss_SR1_A;
+ Histo1DPtr _hist_etmiss_SR1_B;
+ Histo1DPtr _hist_etmiss_SR2_A;
+ Histo1DPtr _hist_etmiss_SR2_B;
+ Histo1DPtr _hist_mSFOS;
+ Histo1DPtr _count_SR1;
+ Histo1DPtr _count_SR2;
//@}
};
Modified: branches/2011-07-aida2yoda/src/Analyses/Makefile.am
==============================================================================
--- branches/2011-07-aida2yoda/src/Analyses/Makefile.am Fri Jun 8 15:17:44 2012 (r3754)
+++ branches/2011-07-aida2yoda/src/Analyses/Makefile.am Fri Jun 8 15:23:42 2012 (r3755)
@@ -68,7 +68,8 @@
ATLAS_2011_I944826.cc \
ATLAS_2012_I1084540.cc \
ATLAS_2012_I1093738.cc \
- ATLAS_2012_I1094568.cc
+ ATLAS_2012_I1094568.cc \
+ ATLAS_2012_I1112263.cc
if ENABLE_PRELIMINARY
RivetATLASAnalyses_la_SOURCES += \
ATLAS_2012_CONF_2012_033.cc \
@@ -86,7 +87,6 @@
ATLAS_2012_I1083318.cc \
ATLAS_2012_I1095236.cc \
ATLAS_2012_CONF_2012_001.cc \
- ATLAS_2012_CONF_2012_023.cc \
ATLAS_2012_CONF_2012_041.cc
endif
More information about the Rivet-svn
mailing list
|