[Rivet-svn] r4009 - in trunk: data/plotinfo src/Analyses

[blackhole at projects.hepforge.org]

Author: buckley
Date: Wed Nov 14 17:52:36 2012
New Revision: 4009

Log:
Merge changesets [3973], [3974] from AIDA branch

Modified:
   trunk/data/plotinfo/ATLAS_2012_I1126136.plot
   trunk/src/Analyses/ATLAS_2012_CONF_2012_001.cc

Modified: trunk/data/plotinfo/ATLAS_2012_I1126136.plot
==============================================================================
--- trunk/data/plotinfo/ATLAS_2012_I1126136.plot	Wed Nov 14 17:52:02 2012	(r4008)
+++ trunk/data/plotinfo/ATLAS_2012_I1126136.plot	Wed Nov 14 17:52:36 2012	(r4009)
@@ -1,34 +1,34 @@
-# BEGIN PLOT /ATLAS_2012_I1186556/count_SR_A
+# BEGIN PLOT /ATLAS_2012_I1126136/count_SR_A
 Title=Number of events in the $E_{\mathrm{T}}^{\rm miss}>130\,$GeV
 XLabel=
 YLabel=Number of Events
 LogY=0
 # END PLOT
-# BEGIN PLOT /ATLAS_2012_I1186556/count_SR_B
+# BEGIN PLOT /ATLAS_2012_I1126136/count_SR_B
 Title=Number of events in the $E_{\mathrm{T}}^{\rm miss}>260\,$GeV
 XLabel=
 YLabel=Number of Events
 LogY=0
 # END PLOT
-# BEGIN PLOT /ATLAS_2012_I1186556/hist_mjjj1
+# BEGIN PLOT /ATLAS_2012_I1126136/hist_mjjj1
 Title=$m_{jjj}$ for closest-triplet
 XLabel=$m_{jjj}$ [GeV]
 YLabel=$\mathrm{Events}/20 \mathrm{GeV}$
 LogY=0
 # END PLOT
-# BEGIN PLOT /ATLAS_2012_I1186556/hist_mjjj2
+# BEGIN PLOT /ATLAS_2012_I1126136/hist_mjjj2
 Title=$m_{jjj}$ for second-triplet
 XLabel=$m_{jjj}$ [GeV]
 YLabel=$\mathrm{Events}/20 \mathrm{GeV}$
 LogY=0
 # END PLOT
-# BEGIN PLOT /ATLAS_2012_I1186556/hist_ETmiss
+# BEGIN PLOT /ATLAS_2012_I1126136/hist_ETmiss
 Title=$E_{\mathrm{T}}^{\rm miss}$
 XLabel=$E_{\mathrm{T}}^{\rm miss}$ [GeV]
 YLabel=$\mathrm{Events}/50 \mathrm{GeV}$
 LogY=0
 # END PLOT
-# BEGIN PLOT /ATLAS_2012_I1186556/hist_mT2
+# BEGIN PLOT /ATLAS_2012_I1126136/hist_mT2
 Title=$m_{T2}$
 XLabel=$m_{T2}$ [GeV]
 YLabel=$\mathrm{Events}/\mathrm{GeV}$

Modified: trunk/src/Analyses/ATLAS_2012_CONF_2012_001.cc
==============================================================================
--- trunk/src/Analyses/ATLAS_2012_CONF_2012_001.cc	Wed Nov 14 17:52:02 2012	(r4008)
+++ trunk/src/Analyses/ATLAS_2012_CONF_2012_001.cc	Wed Nov 14 17:52:36 2012	(r4009)
@@ -147,7 +147,7 @@
       }
 
       // only keep electrons more than R=0.4 from jets
-      ParticleVector recon_e;
+      ParticleVector cand2_e;
       for(unsigned int ie=0;ie<cand_e.size();++ie) {
         const Particle & e = cand_e[ie];
         // at least 0.4 from any jets
@@ -176,12 +176,25 @@
           }
         }
         // if isolated keep it
-        if ( away )
-          recon_e.push_back( e );
+        if ( away ) 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]);
       }
 
       // only keep muons more than R=0.4 from jets
-      ParticleVector recon_mu;
+      ParticleVector cand2_mu;
       for(unsigned int imu=0;imu<cand_mu.size();++imu) {
         const Particle & mu = cand_mu[imu];
         bool away = true;
@@ -210,7 +223,22 @@
           }
         }
         if ( away )
-          recon_mu.push_back( mu );
+          cand2_mu.push_back( mu );
+      }
+
+      // 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]);
       }
 
       // pTmiss
@@ -264,7 +292,7 @@
         if(pT>40.) meff += pT;
       }
 
-      double mSFOS=1e30, mdiff=1e30,mMin=1e30;
+      double mSFOS=1e30, mdiff=1e30;
       // mass of SFOS pairs closest to the Z mass
       for(unsigned int ix=0;ix<recon_e.size();++ix) {
         for(unsigned int iy=ix+1;iy<recon_e.size();++iy) {
@@ -274,9 +302,6 @@
             mSFOS = mtest;
             mdiff = fabs(mtest-90.);
           }
-          else if(mMin>mtest) {
-            mMin = mtest;
-          }
         }
       }
       for(unsigned int ix=0;ix<recon_mu.size();++ix) {
@@ -287,15 +312,8 @@
             mSFOS = mtest;
             mdiff = fabs(mtest-91.118);
           }
-          else if(mMin>mtest) {
-            mMin = mtest;
-          }
         }
       }
-      // cut to reject low mass Drell-Yan
-      if(mMin<=20.) {
-        vetoEvent;
-      }
 
       // make the control plots
       // lepton pT
@@ -356,7 +374,7 @@
       scale(_hist_etmiss_MC   ,norm*20.);
       scale(_hist_mSFOS_MC    ,norm*20.);
       // these are number of events at 2.06fb^-1 per 50 GeV
-      scale(_hist_meff_MC     ,norm*50.);
+      scale(_hist_meff        ,norm*50.);
       scale(_hist_meff_MC     ,norm*50.);
       // these are number of events at 2.06fb^-1
       scale(_hist_njet        ,norm);