[Rivet-svn] r3967 - in branches/2012-06-aidarivet: data/anainfo data/plotinfo src/Analyses

blackhole at projects.hepforge.org blackhole at projects.hepforge.org
Wed Oct 17 15:10:40 BST 2012


Author: richardn
Date: Wed Oct 17 15:10:40 2012
New Revision: 3967

Log:
added ATLAS stop/top partner search paper

Added:
   branches/2012-06-aidarivet/data/anainfo/ATLAS_2012_I1186556.info
      - copied, changed from r3954, branches/2012-06-aidarivet/data/anainfo/ATLAS_2012_I1180197.info
   branches/2012-06-aidarivet/data/plotinfo/ATLAS_2012_I1186556.plot
      - copied, changed from r3954, branches/2012-06-aidarivet/data/plotinfo/ATLAS_2012_I1180197.plot
   branches/2012-06-aidarivet/src/Analyses/ATLAS_2012_I1186556.cc
      - copied, changed from r3954, branches/2012-06-aidarivet/src/Analyses/ATLAS_2012_I1180197.cc
Modified:
   branches/2012-06-aidarivet/data/anainfo/Makefile.am
   branches/2012-06-aidarivet/data/plotinfo/Makefile.am
   branches/2012-06-aidarivet/src/Analyses/Makefile.am

Copied and modified: branches/2012-06-aidarivet/data/anainfo/ATLAS_2012_I1186556.info (from r3954, branches/2012-06-aidarivet/data/anainfo/ATLAS_2012_I1180197.info)
==============================================================================
--- branches/2012-06-aidarivet/data/anainfo/ATLAS_2012_I1180197.info	Fri Oct  5 10:47:32 2012	(r3954, copy source)
+++ branches/2012-06-aidarivet/data/anainfo/ATLAS_2012_I1186556.info	Wed Oct 17 15:10:40 2012	(r3967)
@@ -1,24 +1,30 @@
-Name: ATLAS_2012_I1180197
+Name: ATLAS_2012_I1186556
 Year: 2012
-Summary: Search for supersymmetry at 7 TeV in final states with jets, missing transverse momentum and isolated leptons with the ATLAS detector.
+Summary: Search for a heavy top-quark partner in final states with two leptons.
 Experiment: ATLAS
 Collider: LHC
-SpiresID: 
+InspireID: 1186556
 Status: UNVALIDATED
 Authors:
  - Peter Richardson <Peter.Richardson at durham.ac.uk>
 References:
- - ATLAS-CONF-2012-041
- - arXiv:1208.4688 
+ - arXiv:1209.4186
 RunInfo:
   BSM signal events at 7000 GeV.
 NumEvents: 25000 for BSM signals
 Beams: [p+, p+]
 Energies: [7000]
 Description:
-  'One and two lepton search for supersymmmetric particles by ATLAS at 7 TeV.
-  Event counts in the signal regions are implemented as one-bin histograms.
-  Histograms for effective mass are implemented for the two
-  signal hard lepton signal regions and the ratio of missing transverse energy
-  to effective mass for the soft lepton region. Only the one lepton plots are
-  currently implemented as taken from a conf note originally.'
+  'Search for direct pair production of heavy top-quark partners with 4.7 fb$^{-1}$ integrated luminosity at $\sqrt{s} = 7 TeV$ by the ATLAS experiment. Heavy top-quark partners decaying into a top quark and a neutral non-interacting particle are searched for in events with two leptons in the final state.'
+BibTeX: '@article{Aad:2012uu,
+      author         = "Aad, Georges and others",
+      title          = "{Search for a heavy top-quark partner in final states
+                        with two leptons with the ATLAS detector at the LHC}",
+      collaboration  = "ATLAS Collaboration",
+      year           = "2012",
+      eprint         = "1209.4186",
+      archivePrefix  = "arXiv",
+      primaryClass   = "hep-ex",
+      reportNumber   = "CERN-PH-EP-2012-199",
+      SLACcitation   = "%%CITATION = ARXIV:1209.4186;%%",
+}'

Modified: branches/2012-06-aidarivet/data/anainfo/Makefile.am
==============================================================================
--- branches/2012-06-aidarivet/data/anainfo/Makefile.am	Wed Oct 17 15:06:36 2012	(r3966)
+++ branches/2012-06-aidarivet/data/anainfo/Makefile.am	Wed Oct 17 15:10:40 2012	(r3967)
@@ -60,6 +60,7 @@
   ATLAS_2012_I1125961.info \
   ATLAS_2012_I1117704.info \
   ATLAS_2012_I1180197.info \
+  ATLAS_2012_I1186556.info \
   ATLAS_2012_CONF_2012_103.info \
   ATLAS_2012_CONF_2012_104.info \
   ATLAS_2012_CONF_2012_105.info \

Copied and modified: branches/2012-06-aidarivet/data/plotinfo/ATLAS_2012_I1186556.plot (from r3954, branches/2012-06-aidarivet/data/plotinfo/ATLAS_2012_I1180197.plot)
==============================================================================
--- branches/2012-06-aidarivet/data/plotinfo/ATLAS_2012_I1180197.plot	Fri Oct  5 10:47:32 2012	(r3954, copy source)
+++ branches/2012-06-aidarivet/data/plotinfo/ATLAS_2012_I1186556.plot	Wed Oct 17 15:10:40 2012	(r3967)
@@ -1,39 +1,39 @@
-# BEGIN PLOT /ATLAS_2012_I1180197/count_3jet_channel
-Title=Number of events in the 3 jet channel
+# BEGIN PLOT /ATLAS_2012_I1186556/count_SR_SF
+Title=Number of events in the same-flavour region
 XLabel=
 YLabel=Number of Events
 LogY=0
 # END PLOT
-# BEGIN PLOT /ATLAS_2012_I1180197/count_4jet_channel
-Title=Number of events in the 4 jet channel
-XLabel=
-YLabel=Number of Events
-LogY=0
-# END PLOT
-# BEGIN PLOT /ATLAS_2012_I1180197/count_soft_channel
-Title=Number of events in the soft lepton channel
+# BEGIN PLOT /ATLAS_2012_I1186556/count_SR_OF
+Title=Number of events in the opposite-flavour region
 XLabel=
 YLabel=Number of Events
 LogY=0
 # END PLOT
 
-# BEGIN PLOT /ATLAS_2012_I1180197/hist_m_eff_3jet
-Title=Effective Mass in the 3 jet channe
-XLabel=$m_{\mathrm{eff}}$ [GeV]
-YLabel=$\mathrm{Events}/200 \mathrm{GeV}$
+# BEGIN PLOT /ATLAS_2012_I1186556/hist_mT2_SF_exp
+Title=$m_{T2}$ for the same-flavour region
+XLabel=$m_{T2}$ [GeV]
+YLabel=$\mathrm{Events}/5 \mathrm{GeV}$
 LogY=1
 # END PLOT
-
-# BEGIN PLOT /ATLAS_2012_I1180197/hist_m_eff_4jet
-Title=Effective Mass in the 4 jet channe
-XLabel=$m_{\mathrm{eff}}$ [GeV]
-YLabel=$\mathrm{Events}/200 \mathrm{GeV}$
+# BEGIN PLOT /ATLAS_2012_I1186556/hist_mT2_OF_exp
+Title=$m_{T2}$ for the opposite-flavour region
+XLabel=$m_{T2}$ [GeV]
+YLabel=$\mathrm{Events}/5 \mathrm{GeV}$
 LogY=1
 # END PLOT
 
-# BEGIN PLOT /ATLAS_2012_I1180197/hist_eTmiss_m_eff_soft
-Title=Ratio $E_{\mathrm{T}}^{\rm miss}/m_{\mathrm{eff}}$ in the soft lepton channel
-XLabel=$E_{\mathrm{T}}^{\rm miss}/m_{\mathrm{eff}}$
-YLabel=$\mathrm{Events}/0.1$
+# BEGIN PLOT /ATLAS_2012_I1186556/hist_mT2_SF_MC
+Title=$m_{T2}$ for the same-flavour region
+XLabel=$m_{T2}$ [GeV]
+YLabel=$\mathrm{d}{\sigma}/\mathrm{d}{m_{T2}} [\mathrm{fb/GeV}] $
+
+LogY=1
+# END PLOT
+# BEGIN PLOT /ATLAS_2012_I1186556/hist_mT2_OF_MC
+Title=$m_{T2}$ for the opposite-flavour region
+XLabel=$m_{T2}$ [GeV]
+YLabel=$\mathrm{d}{\sigma}/\mathrm{d}{m_{T2}} [\mathrm{fb/GeV}] $
 LogY=1
 # END PLOT

Modified: branches/2012-06-aidarivet/data/plotinfo/Makefile.am
==============================================================================
--- branches/2012-06-aidarivet/data/plotinfo/Makefile.am	Wed Oct 17 15:06:36 2012	(r3966)
+++ branches/2012-06-aidarivet/data/plotinfo/Makefile.am	Wed Oct 17 15:10:40 2012	(r3967)
@@ -57,6 +57,7 @@
   ATLAS_2012_I1180197.plot \
   ATLAS_2012_I1082009.plot \
   ATLAS_2012_I1183818.plot \
+  ATLAS_2012_I1186556.plot \
   ATLAS_2012_CONF_2012_103.plot \
   ATLAS_2012_CONF_2012_104.plot \
   ATLAS_2012_CONF_2012_105.plot \

Copied and modified: branches/2012-06-aidarivet/src/Analyses/ATLAS_2012_I1186556.cc (from r3954, branches/2012-06-aidarivet/src/Analyses/ATLAS_2012_I1180197.cc)
==============================================================================
--- branches/2012-06-aidarivet/src/Analyses/ATLAS_2012_I1180197.cc	Fri Oct  5 10:47:32 2012	(r3954, copy source)
+++ branches/2012-06-aidarivet/src/Analyses/ATLAS_2012_I1186556.cc	Wed Oct 17 15:10:40 2012	(r3967)
@@ -9,11 +9,12 @@
 #include "Rivet/Projections/IdentifiedFinalState.hh"
 #include "Rivet/Projections/VetoedFinalState.hh"
 #include "Rivet/Projections/FastJets.hh"
+#include "Rivet/Tools/RivetMT2.hh"
 
 namespace Rivet {
 
 
-  class ATLAS_2012_I1180197 : public Analysis {
+  class ATLAS_2012_I1186556 : public Analysis {
   public:
 
     /// @name Constructors etc.
@@ -21,8 +22,8 @@
 
     /// Constructor
 
-    ATLAS_2012_I1180197()
-      : Analysis("ATLAS_2012_I1180197")
+    ATLAS_2012_I1186556()
+      : Analysis("ATLAS_2012_I1186556")
     {    }
 
     //@}
@@ -39,14 +40,14 @@
       // 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, 7.0*GeV);
+      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, 6.0*GeV);
+      IdentifiedFinalState muons(eta_m, 10.0*GeV);
       muons.acceptIdPair(MUON);
       addProjection(muons, "muons");
 
@@ -57,19 +58,19 @@
                    "AntiKtJets04");
 
       // all tracks (to do deltaR with leptons)
-      addProjection(ChargedFinalState(-3.0,3.0,0.5*GeV),"cfs");
+      addProjection(ChargedFinalState(-3.0,3.0,1.*GeV),"cfs");
 
       // for pTmiss
       addProjection(VisibleFinalState(-4.9,4.9),"vfs");
 
       // Book histograms
-      _count_3jet_channel = bookHistogram1D("count_3jet_channel", 1, 0., 1.);
-      _count_4jet_channel = bookHistogram1D("count_4jet_channel", 1, 0., 1.);
-      _count_soft_channel = bookHistogram1D("count_soft_channel", 1, 0., 1.);
-
-      _hist_m_eff_3jet        = bookHistogram1D("hist_m_eff_3jet"       ,  6, 400., 1600.);
-      _hist_m_eff_4jet        = bookHistogram1D("hist_m_eff_4jet"       ,  6, 400., 1600.);
-      _hist_eTmiss_m_eff_soft = bookHistogram1D("hist_eTmiss_m_eff_soft",  6, 0.1 , 0.7  );
+      _count_SR_SF     = bookHistogram1D("count_SR_SF"    , 1, 0., 1.);
+      _count_SR_OF     = bookHistogram1D("count_SR_OF"    , 1, 0., 1.);
+      					                 
+      _hist_mT2_SF_exp = bookHistogram1D("hist_mT2_SF_exp", 40 , 0., 200. );
+      _hist_mT2_OF_exp = bookHistogram1D("hist_mT2_OF_exp", 40 , 0., 200. );
+      _hist_mT2_SF_MC  = bookHistogram1D("hist_mT2_SF_MC" , 500, 0., 1000.);
+      _hist_mT2_OF_MC  = bookHistogram1D("hist_mT2_OF_MC" , 500, 0., 1000.);
 
     }
 
@@ -85,13 +86,13 @@
           cand_jets.push_back(jet);
         }
       }
-
-      // get the candidate "medium" leptons without isolation
-      ParticleVector cand_soft_e,cand_hard_e;
+      // charged tracks for isolation
+      ParticleVector chg_tracks =
+        applyProjection<ChargedFinalState>(event, "cfs").particles();
+      // find the electrons
+      ParticleVector cand_e;
       foreach( const Particle & e,
                applyProjection<IdentifiedFinalState>(event, "elecs").particlesByPt()) {
-        double pT  = e.momentum().perp();
-        double eta = e.momentum().eta();
         // remove any leptons within 0.4 of any candidate jets
         bool e_near_jet = false;
         foreach ( const Jet& jet, cand_jets ) {
@@ -101,19 +102,12 @@
             break;
           }
         }
-        if ( e_near_jet ) continue;
-        // soft selection
-        if(pT>7.&&!(fabs(eta)>1.37&&fabs(eta)<1.52)) {
-          cand_soft_e.push_back(e);
-        }
-        // hard selection
-        if(pT>10.) cand_hard_e.push_back(e);
+	if ( e_near_jet ) continue;
+	cand_e.push_back(e);
       }
-      ParticleVector cand_soft_mu,cand_hard_mu;
+      ParticleVector cand_mu;
       foreach( const Particle & mu,
                applyProjection<IdentifiedFinalState>(event, "muons").particlesByPt()) {
-        double pT  = mu.momentum().perp();
-        double eta = mu.momentum().eta();
         // remove any leptons within 0.4 of any candidate jets
         bool mu_near_jet = false;
         foreach ( const Jet& jet, cand_jets ) {
@@ -123,74 +117,27 @@
           }
         }
         if ( mu_near_jet ) continue;
-        // soft selection
-        if(pT>6.&&!(fabs(eta)>1.37&&fabs(eta)<1.52)) {
-          cand_soft_mu.push_back(mu);
-        }
-        // hard selection
-        if(pT>10.) cand_hard_mu.push_back(mu);
+	cand_mu.push_back(mu);
       }
-      // apply the isolation
-      ParticleVector chg_tracks =
-        applyProjection<ChargedFinalState>(event, "cfs").particles();
-      // pTcone around muon track (hard)
-      ParticleVector recon_hard_mu;
-      foreach ( const Particle & mu, cand_hard_mu ) {
+      // pTcone around muon track
+      ParticleVector recon_mu;
+      foreach ( const Particle & mu, cand_mu ) {
         double pTinCone = -mu.momentum().pT();
-        if(-pTinCone<20.) continue;
         foreach ( const Particle & track, chg_tracks ) {
           if ( deltaR(mu.momentum(),track.momentum()) < 0.2 )
             pTinCone += track.momentum().pT();
         }
-        if ( pTinCone < 1.8*GeV ) recon_hard_mu.push_back(mu);
+        if ( pTinCone < 1.8*GeV ) recon_mu.push_back(mu);
       }
-      // pTcone around muon track (soft)
-      ParticleVector recon_soft_mu;
-      foreach ( const Particle & mu, cand_soft_mu ) {
-        double pTinCone = -mu.momentum().pT();
-        if(-pTinCone>20.) continue;
-        foreach ( const Particle & track, chg_tracks ) {
-          if ( deltaR(mu.momentum(),track.momentum()) < 0.2 )
-            pTinCone += track.momentum().pT();
-        }
-        if ( pTinCone < 1.8*GeV ) recon_soft_mu.push_back(mu);
-      }
-      // pTcone around electron track (hard)
-      ParticleVector recon_hard_e;
-      foreach ( const Particle & e, cand_hard_e ) {
-        double pTinCone = -e.momentum().pT();
-        if(-pTinCone<25.) continue;
-        foreach ( const Particle & track, chg_tracks ) {
-          if ( deltaR(e.momentum(),track.momentum()) < 0.2 )
-            pTinCone += track.momentum().pT();
-        }
-        if ( pTinCone < 0.1 * e.momentum().pT() ) recon_hard_e.push_back(e);
-      }
-      // pTcone around electron track (soft)
-      ParticleVector recon_soft_e;
-      foreach ( const Particle & e, cand_soft_e ) {
+      // pTcone around electron track
+      ParticleVector recon_e;
+      foreach ( const Particle & e, cand_e ) {
         double pTinCone = -e.momentum().pT();
-        if(-pTinCone>25.) continue;
         foreach ( const Particle & track, chg_tracks ) {
           if ( deltaR(e.momentum(),track.momentum()) < 0.2 )
             pTinCone += track.momentum().pT();
         }
-        if ( pTinCone < 0.1 * e.momentum().pT() ) recon_soft_e.push_back(e);
-      }
-
-      // discard jets that overlap with electrons
-      Jets recon_jets;
-      foreach ( const Jet& jet, cand_jets ) {
-        if(fabs(jet.momentum().eta())>2.5||
-           jet.momentum().perp()<25.) continue;
-        bool away_from_e = true;
-        foreach ( const Particle & e, cand_hard_e ) {
-          if ( deltaR(e.momentum(),jet.momentum()) < 0.2 ) {
-            away_from_e = false;
-            break;
-          }
-        }
-        if ( away_from_e ) recon_jets.push_back( jet );
+        if ( pTinCone < 0.1 * e.momentum().pT() ) recon_e.push_back(e);
       }
 
       // pTmiss
@@ -199,78 +146,70 @@
                 applyProjection<VisibleFinalState>(event, "vfs").particles() ) {
         pTmiss -= p.momentum();
       }
-      double eTmiss = pTmiss.pT();
-
-      // both selections require at least 2 jets
-      if(recon_jets.size()<2) vetoEvent;
 
-      // start of meff calculation
-      double HT=0.;
-      foreach( const Jet & jet, recon_jets) {
-        HT += jet.momentum().perp();
-      }
-      double m_eff_inc  = HT+eTmiss;
-
-      // hard selection exactly one candidate
-      // and 1 recon and at least 3 jets
-      if( cand_hard_e.size()  +  cand_hard_mu.size() == 1 &&
-          recon_hard_e.size() + recon_hard_mu.size() == 1 &&
-          recon_jets.size() >= 3 ) {
-        // get the lepton
-        Particle lepton = recon_hard_e.empty() ?
-          recon_hard_mu[0] : recon_hard_e[0];
-        // lepton variables
-        double pT = lepton.momentum().perp();
-        double mT  = 2.*(pT*eTmiss -
-                         lepton.momentum().x()*pTmiss.x() -
-                         lepton.momentum().y()*pTmiss.y());
-        mT = sqrt(mT);
-        HT += pT;
-        m_eff_inc += pT;
-        double m_eff = pT+eTmiss+recon_jets[0].momentum().perp()+
-          recon_jets[1].momentum().perp()+recon_jets[2].momentum().perp();
-        // three jet selection
-        if(recon_jets[0].momentum().perp()>100. &&
-           (recon_jets.size() == 3 ||
-            recon_jets[3].momentum().perp() < 80. ) &&
-           mT>100. && eTmiss>250. && eTmiss/m_eff>0.3) {
-          if(m_eff_inc>1200.) _count_3jet_channel->fill(0.5,weight);
-          _hist_m_eff_3jet->fill(min(1599.,m_eff_inc),weight);
-        }
-        // four jet selecton
-        if(recon_jets.size() >= 4) {
-          m_eff += recon_jets[3].momentum().perp();
-          if(recon_jets[3].momentum().perp() > 80.  &&
-             mT>100. && eTmiss>250. && eTmiss/m_eff>0.2) {
-            if(m_eff_inc>800.) _count_4jet_channel->fill(0.5,weight);
-            _hist_m_eff_4jet->fill(min(1599.,m_eff_inc),weight);
-          }
-        }
-      }
-
-      // soft selection exactly one candidate
-      // and 1 recon and 4 jets
-      if( cand_soft_e.size()  +  cand_soft_mu.size() == 1 &&
-          recon_soft_e.size() + recon_soft_mu.size() == 1 &&
-          recon_jets.size() >= 2 &&
-          recon_jets[0].momentum().perp()>130.) {
-        // get the lepton
-        Particle lepton = recon_soft_e.empty() ?
-          recon_soft_mu[0] : recon_soft_e[0];
-        // lepton variables
-        double pT = lepton.momentum().perp();
-        double mT  = 2.*(pT*eTmiss -
-                         lepton.momentum().x()*pTmiss.x() -
-                         lepton.momentum().y()*pTmiss.y());
-        mT = sqrt(mT);
-        HT += pT;
-        m_eff_inc += pT;
-        double m_eff = pT+eTmiss+recon_jets[0].momentum().perp()+
-          recon_jets[1].momentum().perp();
-        if (mT>100. && eTmiss>250.) {
-          if( eTmiss/m_eff>0.3 ) _count_soft_channel->fill(0.5,weight);
-          _hist_eTmiss_m_eff_soft->fill( eTmiss/m_eff_inc,weight);
-        }
+      // discard jets that overlap with electrons
+      Jets recon_jets;
+      foreach ( const Jet& jet, cand_jets ) {
+        if(fabs(jet.momentum().eta())>2.5||
+           jet.momentum().perp()<20.) 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 );
+      }
+
+      // put leptons into 1 vector and order by pT
+      ParticleVector leptons(recon_e.begin(),recon_e.end());
+      leptons.insert(leptons.begin(),recon_mu.begin(),recon_mu.end());
+      std::sort(leptons.begin(),leptons.end(),cmpParticleByPt);
+
+      // exactly two leptons
+      if(leptons.size() !=2) vetoEvent;
+
+      // hardest lepton pT greater the 25 (20) e(mu)
+      if( (abs(leptons[0].pdgId())==ELECTRON && leptons[0].momentum().perp()<25.) ||
+	  (abs(leptons[0].pdgId())==ELECTRON && leptons[0].momentum().perp()<20.)) 
+	vetoEvent;
+
+      // require opposite sign
+      if(leptons[0].pdgId()*leptons[1].pdgId()>0) vetoEvent;
+
+      // and invariant mass > 20
+      double mll = (leptons[0].momentum()+leptons[1].momentum()).mass();
+      if(mll<20.) vetoEvent;
+
+      // two jets 1st pT > 50 and second pT> 25
+      if(recon_jets.size()<2 || recon_jets[0].momentum().perp()<50. ||
+	 recon_jets[1].momentum().perp()<25.) vetoEvent;
+
+      // calculate mT2
+      double m_T2 = mT2::mT2( leptons[0].momentum(),leptons[1].momentum(),
+			      pTmiss,0.0 ); // zero mass invisibles
+
+      // same flavour region
+      if(leptons[0].pdgId()==-leptons[1].pdgId()) {
+	// remove Z region
+	if(mll>71.&&mll<111.) vetoEvent;
+	// 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)
+	     ++n_b;
+	}
+	if(n_b==0) vetoEvent;
+	_hist_mT2_SF_exp->fill(m_T2,weight);
+	_hist_mT2_SF_MC ->fill(m_T2,weight);
+	if(m_T2>120.) _count_SR_SF->fill(0.5,weight);
+      }
+      // opposite flavour region
+      else {
+	_hist_mT2_OF_exp->fill(m_T2,weight);
+	_hist_mT2_OF_MC ->fill(m_T2,weight);
+	if(m_T2>120.) _count_SR_OF->fill(0.5,weight);
       }
     }
     //@}
@@ -279,12 +218,12 @@
     void finalize() {
 
       double norm = 4.7* crossSection()/sumOfWeights()/femtobarn;
-      scale(_count_3jet_channel    ,norm);
-      scale(_count_4jet_channel    ,norm);
-      scale(_count_soft_channel    ,norm);
-      scale(_hist_m_eff_3jet       ,200.*norm);
-      scale(_hist_m_eff_4jet       ,200.*norm);
-      scale(_hist_eTmiss_m_eff_soft,0.1*norm);
+      scale(_count_SR_SF    ,   norm);
+      scale(_count_SR_OF    ,   norm);
+      scale(_hist_mT2_SF_exp,5.*norm);
+      scale(_hist_mT2_OF_exp,5.*norm);
+      scale(_hist_mT2_SF_MC ,   norm/4.7);
+      scale(_hist_mT2_OF_MC ,   norm/4.7);
 
     }
 
@@ -292,18 +231,18 @@
 
     /// @name Histograms
     //@{
-    AIDA::IHistogram1D* _count_3jet_channel;
-    AIDA::IHistogram1D* _count_4jet_channel;
-    AIDA::IHistogram1D* _count_soft_channel;
-
-    AIDA::IHistogram1D* _hist_m_eff_3jet;
-    AIDA::IHistogram1D* _hist_m_eff_4jet;
-    AIDA::IHistogram1D* _hist_eTmiss_m_eff_soft;
+    AIDA::IHistogram1D* _count_SR_SF;
+    AIDA::IHistogram1D* _count_SR_OF;
+
+    AIDA::IHistogram1D* _hist_mT2_SF_exp;
+    AIDA::IHistogram1D* _hist_mT2_OF_exp;
+    AIDA::IHistogram1D* _hist_mT2_SF_MC;
+    AIDA::IHistogram1D* _hist_mT2_OF_MC;
     //@}
 
   };
 
   // The hook for the plugin system
-  DECLARE_RIVET_PLUGIN(ATLAS_2012_I1180197);
+  DECLARE_RIVET_PLUGIN(ATLAS_2012_I1186556);
 
 }

Modified: branches/2012-06-aidarivet/src/Analyses/Makefile.am
==============================================================================
--- branches/2012-06-aidarivet/src/Analyses/Makefile.am	Wed Oct 17 15:06:36 2012	(r3966)
+++ branches/2012-06-aidarivet/src/Analyses/Makefile.am	Wed Oct 17 15:10:40 2012	(r3967)
@@ -96,6 +96,7 @@
     ATLAS_2012_I1082009.cc \
     ATLAS_2012_I946427.cc \
     ATLAS_2012_I1180197.cc \
+    ATLAS_2012_I1186556.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