|
|
[Rivet-svn] r3732 - in trunk: data/anainfo data/plotinfo src/Analyses
blackhole at projects.hepforge.org
blackhole at projects.hepforge.org
Sun May 20 10:04:54 BST 2012
Author: richardn
Date: Sun May 20 10:04:53 2012
New Revision: 3732
Log:
added ATLAS_2012_I946427 diphoton production + MET
Added:
trunk/data/anainfo/ATLAS_2012_I946427.info
- copied, changed from r3731, trunk/data/anainfo/ATLAS_2012_I1095236.info
trunk/data/plotinfo/ATLAS_2012_I946427.plot
- copied, changed from r3731, trunk/data/plotinfo/ATLAS_2012_CONF_2012_001.plot
trunk/src/Analyses/ATLAS_2012_I946427.cc
- copied, changed from r3731, trunk/src/Analyses/ATLAS_2012_I1095236.cc
Modified:
trunk/data/anainfo/Makefile.am
trunk/data/plotinfo/Makefile.am
trunk/src/Analyses/Makefile.am
Copied and modified: trunk/data/anainfo/ATLAS_2012_I946427.info (from r3731, trunk/data/anainfo/ATLAS_2012_I1095236.info)
==============================================================================
--- trunk/data/anainfo/ATLAS_2012_I1095236.info Thu May 17 20:38:03 2012 (r3731, copy source)
+++ trunk/data/anainfo/ATLAS_2012_I946427.info Sun May 20 10:04:53 2012 (r3732)
@@ -1,40 +1,38 @@
-Name: ATLAS_2012_I1095236
-Year: 2011
-Summary: b-jets search for supersymmetry with 0- and 1-leptons
+Name: ATLAS_2012_I946427
+Year: 2012
+Summary: Search for supersymmetry with diphotons and mising Transverse Momentum
Experiment: ATLAS
Collider: LHC
-SpiresID: 1095236
+InspireID: 946427
Status: UNVALIDATED
Authors:
- Peter Richardson <Peter.Richardson at durham.ac.uk>
References:
- - arXiv:1203.6193
+ - arXiv:1111.4116
+ - Phys. Lett. B710 (2012) 519-537
RunInfo:
BSM signal events at 7000 GeV.
NumEvents: 25000 for BSM signals
Beams: [p+, p+]
Energies: [7000]
Description:
- 'Search for supersymmmetric particles by ATLAS at 7 TeV in events with b-jets,
- large missing
- energy, and zero or one leptons.
- Event counts in six zero lepton and two one lepton signal regions are implemented
- as one-bin histograms.
- Histograms for missing transverse energy, and effective mass are also
- implemented for some signal regions.'
-ibKey: ATLAS2012:1095236
-BibTeX: '@article{ATLAS2012:1095236,
+ 'Search for diphoton events with large missing transverse momentum
+ with integrated luminosity 1.07$\text{fb}^{-1}$ at $\sqrt{s}=7$.
+ No excess of events was observed.'
+BibKey: Aad:2011zj
+BibTeX: '@article{Aad:2011zj,
author = "Aad, Georges and others",
- title = "{Search for supersymmetry in pp collisions at sqrt(s) = 7
- TeV in final states with missing transverse momentum and
- b-jets with the ATLAS detector}",
+ title = "{Search for Diphoton Events with Large Missing Transverse
+ Momentum in $1\,\text{fb}^{-1}$ of 7 TeV Proton-Proton Collision Data
+ with the ATLAS Detector}",
collaboration = "ATLAS Collaboration",
+ journal = "Phys.Lett.",
+ volume = "B710",
+ pages = "519-537",
year = "2012",
- note = "15 pages plus author list (28 pages total), 11 figures, 8
- tables, submitted to Physical Review D",
- eprint = "1203.6193",
+ eprint = "1111.4116",
archivePrefix = "arXiv",
primaryClass = "hep-ex",
- reportNumber = "CERN-PH-EP-2012-052",
- SLACcitation = "%%CITATION = ARXIV:1203.6193;%%",
+ reportNumber = "CERN-PH-EP-2011-160",
+ SLACcitation = "%%CITATION = ARXIV:1111.4116;%%",
}'
Modified: trunk/data/anainfo/Makefile.am
==============================================================================
--- trunk/data/anainfo/Makefile.am Thu May 17 20:38:03 2012 (r3731)
+++ trunk/data/anainfo/Makefile.am Sun May 20 10:04:53 2012 (r3732)
@@ -49,6 +49,7 @@
ATLAS_2012_I1094568.info \
ATLAS_2012_I1095236.info \
ATLAS_2012_I1082009.info \
+ ATLAS_2012_I946427.info \
ATLAS_2012_CONF_2012_001.info \
ATLAS_2012_I1112263.info \
ATLAS_2012_CONF_2012_033.info \
Copied and modified: trunk/data/plotinfo/ATLAS_2012_I946427.plot (from r3731, trunk/data/plotinfo/ATLAS_2012_CONF_2012_001.plot)
==============================================================================
--- trunk/data/plotinfo/ATLAS_2012_CONF_2012_001.plot Thu May 17 20:38:03 2012 (r3731, copy source)
+++ trunk/data/plotinfo/ATLAS_2012_I946427.plot Sun May 20 10:04:53 2012 (r3732)
@@ -1,116 +1,19 @@
-# BEGIN PLOT /ATLAS_2012_CONF_2012_001/count_SR1
-Title=Number of events in the SR1 signal region
+# BEGIN PLOT /ATLAS_2012_I946427/count_SR
+Title=Number of signal events
XLabel=
YLabel=Number of Events
LogY=0
# END PLOT
-# BEGIN PLOT /ATLAS_2012_CONF_2012_001/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_001/d01-x01-y01
-Title=$E^e_{\mathrm{T}}(p^\mu_{\mathrm{T}})$ for the leading lepton
-XLabel=$E^e_{\mathrm{T}}(p^\mu_{\mathrm{T}})$ [GeV]
-YLabel=$\mathrm{Events}/10 \mathrm{GeV}$
-LogY=1
-# END PLOT
-# BEGIN PLOT /ATLAS_2012_CONF_2012_001/d02-x01-y01
-Title=$E^e_{\mathrm{T}}(p^\mu_{\mathrm{T}})$ for the $2^{\rm nd}$ leading lepton
-XLabel=$E^e_{\mathrm{T}}(p^\mu_{\mathrm{T}})$ [GeV]
-YLabel=$\mathrm{Events}/10 \mathrm{GeV}$
-LogY=1
-# END PLOT
-# BEGIN PLOT /ATLAS_2012_CONF_2012_001/d03-x01-y01
-Title=$E^e_{\mathrm{T}}(p^\mu_{\mathrm{T}})$ for the $3^{\rm rd}$ leading lepton
-XLabel=$E^e_{\mathrm{T}}(p^\mu_{\mathrm{T}})$ [GeV]
-YLabel=$\mathrm{Events}/5 \mathrm{GeV}$
-LogY=1
-# END PLOT
-# BEGIN PLOT /ATLAS_2012_CONF_2012_001/d04-x01-y01
-Title=$E^e_{\mathrm{T}}(p^\mu_{\mathrm{T}})$ for the $4^{\rm th}$ leading lepton
-XLabel=$E^e_{\mathrm{T}}(p^\mu_{\mathrm{T}})$ [GeV]
-YLabel=$\mathrm{Events}/5 \mathrm{GeV}$
-LogY=1
-# END PLOT
-
-# BEGIN PLOT /ATLAS_2012_CONF_2012_001/d05-x01-y01
-Title=Number of Jets
-XLabel=Number of Jets
-YLabel=Events
-LogY=1
-# END PLOT
-# BEGIN PLOT /ATLAS_2012_CONF_2012_001/d06-x01-y01
-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_001/d07-x01-y01
-Title=Mass of SFOF lepton pair closest to the $Z^0$ mass
-XLabel=$M_{\mathrm{SFOS}}$ [GeV]
-YLabel=$\mathrm{Events}/20$\,GeV
-LogY=1
-# END PLOT
-# BEGIN PLOT /ATLAS_2012_CONF_2012_001/d08-x01-y01
-Title=Effective Mass
-XLabel=$M_{\mathrm{eff}}$ [GeV]
-YLabel=$\mathrm{Events}/50$\,GeV
-LogY=1
-# END PLOT
-
-
-
-
-# BEGIN PLOT /ATLAS_2012_CONF_2012_001/hist_lepton_pT_1
-Title=$E^e_{\mathrm{T}}(p^\mu_{\mathrm{T}})$ for the leading lepton
-XLabel=$E^e_{\mathrm{T}}(p^\mu_{\mathrm{T}})$ [GeV]
-YLabel=$\mathrm{Events}/10 \mathrm{GeV}$
-LogY=1
-# END PLOT
-# BEGIN PLOT /ATLAS_2012_CONF_2012_001/hist_lepton_pT_2
-Title=$E^e_{\mathrm{T}}(p^\mu_{\mathrm{T}})$ for the $2^{\rm nd}$ leading lepton
-XLabel=$E^e_{\mathrm{T}}(p^\mu_{\mathrm{T}})$ [GeV]
+# BEGIN PLOT /ATLAS_2012_I946427/hist_ET_photon
+Title=$E^\gamma_{\mathrm{T}}$ for the leading photon
+XLabel=$E^\gamma_{\mathrm{T}}$ [GeV]
YLabel=$\mathrm{Events}/10 \mathrm{GeV}$
LogY=1
# END PLOT
-# BEGIN PLOT /ATLAS_2012_CONF_2012_001/hist_lepton_pT_3
-Title=$E^e_{\mathrm{T}}(p^\mu_{\mathrm{T}})$ for the $3^{\rm rd}$ leading lepton
-XLabel=$E^e_{\mathrm{T}}(p^\mu_{\mathrm{T}})$ [GeV]
-YLabel=$\mathrm{Events}/5 \mathrm{GeV}$
-LogY=1
-# END PLOT
-# BEGIN PLOT /ATLAS_2012_CONF_2012_001/hist_lepton_pT_4
-Title=$E^e_{\mathrm{T}}(p^\mu_{\mathrm{T}})$ for the $4^{\rm th}$ leading lepton
-XLabel=$E^e_{\mathrm{T}}(p^\mu_{\mathrm{T}})$ [GeV]
-YLabel=$\mathrm{Events}/5 \mathrm{GeV}$
-LogY=1
-# END PLOT
-
-# BEGIN PLOT /ATLAS_2012_CONF_2012_001/hist_njet
-Title=Number of Jets
-XLabel=Number of Jets
-YLabel=Events
-LogY=1
-# END PLOT
-# BEGIN PLOT /ATLAS_2012_CONF_2012_001/hist_etmiss
+# BEGIN PLOT /ATLAS_2012_I946427/hist_met
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_001/hist_m_SFOS
-Title=Mass of SFOF lepton pair closest to the $Z^0$ mass
-XLabel=$M_{\mathrm{SFOS}}$ [GeV]
-YLabel=$\mathrm{Events}/20$\,GeV
-LogY=1
-# END PLOT
-# BEGIN PLOT /ATLAS_2012_CONF_2012_001/hist_m_eff
-Title=Effective Mass
-XLabel=$M_{\mathrm{eff}}$ [GeV]
-YLabel=$\mathrm{Events}/50$\,GeV
+YLabel=$\mathrm{Events}/5$\,GeV
LogY=1
# END PLOT
Modified: trunk/data/plotinfo/Makefile.am
==============================================================================
--- trunk/data/plotinfo/Makefile.am Thu May 17 20:38:03 2012 (r3731)
+++ trunk/data/plotinfo/Makefile.am Sun May 20 10:04:53 2012 (r3732)
@@ -45,6 +45,7 @@
ATLAS_2012_I1093738.plot \
ATLAS_2012_I1094568.plot \
ATLAS_2012_I1095236.plot \
+ ATLAS_2012_I946427.plot \
ATLAS_2012_CONF_2012_001.plot \
ATLAS_2012_I1112263.plot \
ATLAS_2012_CONF_2012_033.plot \
Copied and modified: trunk/src/Analyses/ATLAS_2012_I946427.cc (from r3731, trunk/src/Analyses/ATLAS_2012_I1095236.cc)
==============================================================================
--- trunk/src/Analyses/ATLAS_2012_I1095236.cc Thu May 17 20:38:03 2012 (r3731, copy source)
+++ trunk/src/Analyses/ATLAS_2012_I946427.cc Sun May 20 10:04:53 2012 (r3732)
@@ -14,15 +14,15 @@
namespace Rivet {
/// @author Peter Richardson
- class ATLAS_2012_I1095236 : public Analysis {
+ class ATLAS_2012_I946427 : public Analysis {
public:
/// @name Constructors etc.
//@{
/// Constructor
- ATLAS_2012_I1095236()
- : Analysis("ATLAS_2012_I1095236")
+ ATLAS_2012_I946427()
+ : Analysis("ATLAS_2012_I946427")
{ }
//@}
@@ -36,130 +36,66 @@
/// Book histograms and initialise projections before the run
void init() {
- // Projection to find the electrons
- std::vector<std::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;
- eta_m.push_back(make_pair(-2.4,2.4));
- IdentifiedFinalState muons(eta_m, 10.0*GeV);
- muons.acceptIdPair(MUON);
- addProjection(muons, "muons");
-
- // Jet finder
- VetoedFinalState vfs;
- vfs.addVetoPairId(MUON);
- addProjection(FastJets(vfs, FastJets::ANTIKT, 0.4), "AntiKtJets04");
-
- // All tracks (to do deltaR with leptons)
- addProjection(ChargedFinalState(-3.0,3.0),"cfs");
+ // photons
+ IdentifiedFinalState photonfs(-1.81, 1.81, 25.0*GeV);
+ photonfs.acceptId(PHOTON);
+ addProjection(photonfs, "Photon");
+
+ //
+ FinalState fs;
+ addProjection(fs, "FS");
// Used for pTmiss
addProjection(VisibleFinalState(-4.9,4.9),"vfs");
// Book histograms
- _count_SR0_A1 = bookHistogram1D("count_SR0_A1", 1, 0., 1.);
- _count_SR0_B1 = bookHistogram1D("count_SR0_B1", 1, 0., 1.);
- _count_SR0_C1 = bookHistogram1D("count_SR0_C1", 1, 0., 1.);
- _count_SR0_A2 = bookHistogram1D("count_SR0_A2", 1, 0., 1.);
- _count_SR0_B2 = bookHistogram1D("count_SR0_B2", 1, 0., 1.);
- _count_SR0_C2 = bookHistogram1D("count_SR0_C2", 1, 0., 1.);
- _count_SR1_D = bookHistogram1D("count_SR1_D" , 1, 0., 1.);
- _count_SR1_E = bookHistogram1D("count_SR1_E" , 1, 0., 1.);
-
- _hist_meff_SR0_A1 = bookHistogram1D("hist_m_eff_SR0_A1", 14, 400., 1800.);
- _hist_meff_SR0_A2 = bookHistogram1D("hist_m_eff_SR0_A2", 14, 400., 1800.);
- _hist_meff_SR1_D_e = bookHistogram1D("hist_meff_SR1_D_e" , 16, 600., 2200.);
- _hist_meff_SR1_D_mu = bookHistogram1D("hist_meff_SR1_D_mu", 16, 600., 2200.);
-
- _hist_met_SR0_A1 = bookHistogram1D("hist_met_SR0_A1", 14, 0., 700.);
- _hist_met_SR0_A2 = bookHistogram1D("hist_met_SR0_A2", 14, 0., 700.);
- _hist_met_SR0_D_e = bookHistogram1D("hist_met_SR1_D_e" , 15, 0., 600.);
- _hist_met_SR0_D_mu = bookHistogram1D("hist_met_SR1_D_mu", 15, 0., 600.);
+ _count_SR = bookHistogram1D("count_SR", 1, 0., 1.);
+
+ _hist_ET_photon = bookHistogram1D("hist_ET_photon", 48 , 20., 500.);
+ _hist_met = bookHistogram1D("hist_met" , 100, 0., 500.);
}
/// Perform the per-event analysis
void analyze(const Event& event) {
- const double weight = event.weight();
-
- Jets cand_jets;
- const Jets jets = applyProjection<FastJets>(event, "AntiKtJets04").jetsByPt(20.0*GeV);
- foreach (const Jet& jet, jets) {
- if ( fabs( jet.momentum().eta() ) < 2.8 ) {
- cand_jets.push_back(jet);
- }
- }
- const ParticleVector cand_e = applyProjection<IdentifiedFinalState>(event, "elecs").particlesByPt();
+ const double weight = event.weight();
- const ParticleVector cand_mu = applyProjection<IdentifiedFinalState>(event, "muons").particlesByPt();
- // Resolve jet-lepton overlap for jets with |eta| < 2.8
- Jets recon_jets;
- foreach ( const Jet& jet, cand_jets ) {
- if ( fabs( jet.momentum().eta() ) >= 2.8 ) continue;
- 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 );
- }
-
- // get the loose leptons used to define the 0 lepton channel
- ParticleVector loose_e, loose_mu;
- foreach ( const Particle & e, cand_e ) {
- bool away = true;
- foreach ( const Jet& jet, recon_jets ) {
- if ( deltaR(e.momentum(),jet.momentum()) < 0.4 ) {
- away = false;
- break;
- }
+ // require at least 2 photons in final state
+ ParticleVector photons =
+ applyProjection<IdentifiedFinalState>(event, "Photon").particlesByPt();
+ if (photons.size() < 2) {
+ vetoEvent;
+ }
+
+ // Loop over photons and fill vector of isolated ones
+ ParticleVector fs = applyProjection<FinalState>(event, "FS").particles();
+ ParticleVector isolated_photons;
+ foreach (const Particle& photon, photons) {
+ // remove photons in crack
+ double eta_P = photon.momentum().eta();
+ if (fabs(eta_P)>=1.37 && fabs(eta_P)<1.52) continue;
+
+ double phi_P = photon.momentum().phi();
+
+ FourMomentum mom_in_EtCone = -photon.momentum();
+ foreach (const Particle& p, fs) {
+ // check if it's in the cone of .2
+ if (deltaR(eta_P, phi_P, p.momentum().eta(),
+ p.momentum().phi()) >= 0.2) continue;
+ mom_in_EtCone += p.momentum();
}
- if ( away ) loose_e.push_back( e );
- }
- foreach ( const Particle & mu, cand_mu ) {
- bool away = true;
- foreach ( const Jet& jet, recon_jets ) {
- if ( deltaR(mu.momentum(),jet.momentum()) < 0.4 ) {
- away = false;
- break;
- }
- }
- if ( away ) loose_mu.push_back( mu );
- }
- // tight leptons for the 1-lepton channel
- ParticleVector tight_mu;
- ParticleVector chg_tracks =
- applyProjection<ChargedFinalState>(event, "cfs").particles();
- foreach ( const Particle & mu, loose_mu) {
- if(mu.momentum().perp()<20.) continue;
- double pTinCone = -mu.momentum().pT();
- foreach ( const Particle & track, chg_tracks ) {
- if ( deltaR(mu.momentum(),track.momentum()) <= 0.2 )
- pTinCone += track.momentum().pT();
- }
- if ( pTinCone < 1.8*GeV )
- tight_mu.push_back(mu);
+ // apply isolation
+ if(mom_in_EtCone.Et()>5.) continue;
+
+ // add photon to list of isolated ones
+ isolated_photons.push_back(photon);
}
- ParticleVector tight_e;
- foreach ( const Particle & e, loose_e ) {
- if(e.momentum().perp()<25.) 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) {
- tight_e.push_back(e);
- }
+
+ // need two isolated photons
+ if(isolated_photons.size() < 2 ) {
+ vetoEvent;
}
// pTmiss
@@ -171,141 +107,22 @@
}
double eTmiss = pTmiss.pT();
- // get the number of b-tagged jets
- unsigned int ntagged=0;
- foreach (const Jet & jet, recon_jets ) {
- if(jet.momentum().perp()>50. && abs(jet.momentum().eta())<2.5 &&
- jet.containsBottom() && rand()/static_cast<double>(RAND_MAX)<=0.60)
- ++ntagged;
- }
-
- // ATLAS calo problem
- if(rand()/static_cast<double>(RAND_MAX)<=0.42) {
- foreach ( const Jet & jet, recon_jets ) {
- double eta = jet.momentum().rapidity();
- double phi = jet.momentum().azimuthalAngle(MINUSPI_PLUSPI);
- if(jet.momentum().perp()>50 && eta>-0.1&&eta<1.5&&phi>-0.9&&phi<-0.5)
- vetoEvent;
- }
- }
+ _hist_ET_photon->fill(isolated_photons[0].momentum().Et(),weight);
+ _hist_met ->fill(eTmiss ,weight);
- // at least 1 b tag
- if(ntagged==0) vetoEvent;
-
- // minumum Et miss
- if(eTmiss<80.) vetoEvent;
-
- // at least 3 jets pT > 50
- if(recon_jets.size()<3 || recon_jets[2].momentum().perp()<50.)
- vetoEvent;
-
- // m_eff
- double m_eff = eTmiss;
- for(unsigned int ix=0;ix<3;++ix)
- m_eff += recon_jets[ix].momentum().perp();
-
- // delta Phi
- double min_dPhi = 999.999;
- double pTmiss_phi = pTmiss.phi();
- for(unsigned int ix=0;ix<3;++ix) {
- min_dPhi = min( min_dPhi, deltaPhi( pTmiss_phi, recon_jets[ix].momentum().phi() ) );
- }
-
- // 0-lepton channels
- if(loose_e.empty() && loose_mu.empty() &&
- recon_jets[0].momentum().perp()>130. && eTmiss>130. &&
- eTmiss/m_eff>0.25 && min_dPhi>0.4) {
- // jet charge cut
- bool jetCharge = true;
- for(unsigned int ix=0;ix<3;++ix) {
- if(fabs(recon_jets[ix].momentum().eta())>2.) continue;
- double trackpT=0;
- foreach(const Particle & p, recon_jets[ix].particles()) {
- if(PID::threeCharge(p.pdgId())==0) continue;
- trackpT += p.momentum().perp();
- }
- if(trackpT/recon_jets[ix].momentum().perp()<0.05)
- jetCharge = false;
- }
- if(jetCharge) {
- // SR0-A region
- if(m_eff>500.) {
- _count_SR0_A1->fill(0.5,weight);
- _hist_meff_SR0_A1->fill(m_eff,weight);
- _hist_met_SR0_A1 ->fill(eTmiss,weight);
- if(ntagged>=2) {
- _count_SR0_A2->fill(0.5,weight);
- _hist_meff_SR0_A2->fill(m_eff,weight);
- _hist_met_SR0_A2 ->fill(eTmiss,weight);
- }
- }
- // SR0-B
- if(m_eff>700.) {
- _count_SR0_B1->fill(0.5,weight);
- if(ntagged>=2) _count_SR0_B2->fill(0.5,weight);
- }
- // SR0-C
- if(m_eff>900.) {
- _count_SR0_C1->fill(0.5,weight);
- if(ntagged>=2) _count_SR0_C2->fill(0.5,weight);
- }
- }
- }
-
- // 1-lepton channels
- if(tight_e.size() + tight_mu.size() == 1 &&
- recon_jets.size()>=4 && recon_jets[3].momentum().perp()>50.&&
- recon_jets[0].momentum().perp()>60.) {
- Particle lepton = tight_e.empty() ? tight_mu[0] : tight_e[0];
- m_eff += lepton.momentum().perp() + recon_jets[3].momentum().perp();
- // transverse mass cut
- double mT = 2.*(lepton.momentum().perp()*eTmiss-
- lepton.momentum().x()*pTmiss.x()-
- lepton.momentum().y()*pTmiss.y());
- mT = sqrt(mT);
- if(mT>100.&&m_eff>700.) {
- // D region
- _count_SR1_D->fill(0.5,weight);
- if(abs(lepton.pdgId())==ELECTRON) {
- _hist_meff_SR1_D_e->fill(m_eff,weight);
- _hist_met_SR0_D_e->fill(eTmiss,weight);
- }
- else {
- _hist_meff_SR1_D_mu->fill(m_eff,weight);
- _hist_met_SR0_D_mu->fill(eTmiss,weight);
- }
- // E region
- if(eTmiss>200.) {
- _count_SR1_E->fill(0.5,weight);
- }
- }
- }
+ if(eTmiss>125.) _count_SR->fill(0.5,weight);
}
void finalize() {
- double norm = crossSection()/femtobarn*2.05/sumOfWeights();
- // these are number of events at 2.05fb^-1 per 100 GeV
- scale( _hist_meff_SR0_A1 , 100. * norm );
- scale( _hist_meff_SR0_A2 , 100. * norm );
- scale( _hist_meff_SR1_D_e , 100. * norm );
- scale( _hist_meff_SR1_D_mu , 100. * norm );
- // these are number of events at 2.05fb^-1 per 50 GeV
- scale( _hist_met_SR0_A1, 50. * norm );
- scale( _hist_met_SR0_A2, 40. * norm );
- // these are number of events at 2.05fb^-1 per 40 GeV
- scale( _hist_met_SR0_D_e , 40. * norm );
- scale( _hist_met_SR0_D_mu, 40. * norm );
- // these are number of events at 2.05fb^-1
- scale(_count_SR0_A1,norm);
- scale(_count_SR0_B1,norm);
- scale(_count_SR0_C1,norm);
- scale(_count_SR0_A2,norm);
- scale(_count_SR0_B2,norm);
- scale(_count_SR0_C2,norm);
- scale(_count_SR1_D ,norm);
- scale(_count_SR1_E ,norm);
+ double norm = crossSection()/femtobarn*1.07/sumOfWeights();
+ // these are number of events at 1.07fb^-1 per 10 GeV
+ scale( _hist_ET_photon, 10. * norm );
+ // these are number of events at 1.07fb^-1 per 5 GeV
+ scale( _hist_met, 5. * norm );
+ // these are number of events at 1.07fb^-1
+ scale(_count_SR,norm);
}
//@}
@@ -313,28 +130,14 @@
private:
- AIDA::IHistogram1D* _count_SR0_A1;
- AIDA::IHistogram1D* _count_SR0_B1;
- AIDA::IHistogram1D* _count_SR0_C1;
- AIDA::IHistogram1D* _count_SR0_A2;
- AIDA::IHistogram1D* _count_SR0_B2;
- AIDA::IHistogram1D* _count_SR0_C2;
- AIDA::IHistogram1D* _count_SR1_D;
- AIDA::IHistogram1D* _count_SR1_E;
-
- AIDA::IHistogram1D* _hist_meff_SR0_A1;
- AIDA::IHistogram1D* _hist_meff_SR0_A2;
- AIDA::IHistogram1D* _hist_meff_SR1_D_e;
- AIDA::IHistogram1D* _hist_meff_SR1_D_mu;
- AIDA::IHistogram1D* _hist_met_SR0_A1;
- AIDA::IHistogram1D* _hist_met_SR0_A2;
- AIDA::IHistogram1D* _hist_met_SR0_D_e;
- AIDA::IHistogram1D* _hist_met_SR0_D_mu;
+ AIDA::IHistogram1D* _count_SR;
+ AIDA::IHistogram1D* _hist_ET_photon;
+ AIDA::IHistogram1D* _hist_met;
};
// This global object acts as a hook for the plugin system
- DECLARE_RIVET_PLUGIN(ATLAS_2012_I1095236);
+ DECLARE_RIVET_PLUGIN(ATLAS_2012_I946427);
}
Modified: trunk/src/Analyses/Makefile.am
==============================================================================
--- trunk/src/Analyses/Makefile.am Thu May 17 20:38:03 2012 (r3731)
+++ trunk/src/Analyses/Makefile.am Sun May 20 10:04:53 2012 (r3732)
@@ -89,6 +89,7 @@
ATLAS_2012_I1083318.cc \
ATLAS_2012_I1095236.cc \
ATLAS_2012_I1082009.cc \
+ ATLAS_2012_I946427.cc \
ATLAS_2012_CONF_2012_041.cc
endif
More information about the Rivet-svn
mailing list
|