[Rivet-svn] r3386 - trunk/src/Analyses

[blackhole at projects.hepforge.org]

Author: richardn
Date: Thu Sep 22 15:35:48 2011
New Revision: 3386

Log:
removed testing printout#]

Modified:
   trunk/src/Analyses/ATLAS_2011_CONF_2011_090.cc

Modified: trunk/src/Analyses/ATLAS_2011_CONF_2011_090.cc
==============================================================================
--- trunk/src/Analyses/ATLAS_2011_CONF_2011_090.cc	Thu Sep 22 15:17:24 2011	(r3385)
+++ trunk/src/Analyses/ATLAS_2011_CONF_2011_090.cc	Thu Sep 22 15:35:48 2011	(r3386)
@@ -63,7 +63,7 @@
 
       // Jet finder
       VetoedFinalState vfs;
-      vfs.addVetoPairId(MUON);
+      vfs.addVetoPair(MUON);
       addProjection(FastJets(vfs, FastJets::ANTIKT, 0.4),
                    "AntiKtJets04");
 
@@ -94,7 +94,7 @@
 
     /// Perform the per-event analysis
     void analyze(const Event& event) {
-      cerr << __LINE__ << "\n";
+
       const double weight = event.weight();
 
 
@@ -104,7 +104,6 @@
        	MSG_DEBUG("electrons in veto region");
        	vetoEvent;
       }
-      cerr << __LINE__ << "\n";
 
       Jets cand_jets;
       foreach ( const Jet& jet,
@@ -113,7 +112,6 @@
           cand_jets.push_back(jet);
         }
       }
-      cerr << __LINE__ << "\n";
 
       ParticleVector candtemp_e =
 	applyProjection<IdentifiedFinalState>(event, "elecs").particlesByPt();
@@ -123,7 +121,6 @@
 	applyProjection<ChargedFinalState>(event, "cfs").particles();
       ParticleVector cand_mu;
       ParticleVector cand_e;
-      cerr << __LINE__ << "\n";
 
 
       // pTcone around muon track
@@ -137,7 +134,6 @@
 	  cand_mu.push_back(mu);
       }
 
-      cerr << __LINE__ << "\n";
 
       // pTcone around electron
       foreach ( const Particle e, candtemp_e ) {
@@ -150,7 +146,6 @@
 	  cand_e.push_back(e);
       }
 
-      cerr << __LINE__ << "\n";
 
 
       // discard jets that overlap with electrons
@@ -166,7 +161,6 @@
 	  if ( away_from_e )
 	    cand_jets_2.push_back( jet );
       }
-      cerr << __LINE__ << "\n";
 
       // only consider leptons far from jet
       ParticleVector recon_e, recon_mu;
@@ -180,7 +174,6 @@
         if ( ! e_near_jet )
           recon_e.push_back( e );
       }
-      cerr << __LINE__ << "\n";
 
       foreach ( const Particle & mu, cand_mu ) {
 	bool mu_near_jet = false;
@@ -191,7 +184,6 @@
 	if ( ! mu_near_jet )
 	  recon_mu.push_back( mu );
       }
-      cerr << __LINE__ << "\n";
       
       // pTmiss
       ParticleVector vfs_particles
@@ -202,7 +194,6 @@
       }
       double eTmiss = pTmiss.pT();
 
-      cerr << __LINE__ << "\n";
 
       // final jet filter
       Jets recon_jets;
@@ -212,7 +203,6 @@
 
 
 
-      cerr << __LINE__ << "\n";
 
       // ==================== observables ====================
 
@@ -225,24 +215,20 @@
 	if ( fabs(jet.momentum().eta()) < 2.8 )
 	  Njets+=1;
       }
-      cerr << __LINE__ << "\n";
       if ( Njets < 3 ) {
 	MSG_DEBUG("Only " << Njets << " jets w/ eta<2.8 left");
 	vetoEvent;
       }
 
-      cerr << __LINE__ << "\n";
       if ( recon_jets[0].momentum().pT() <= 60.0 * GeV ) {
 	MSG_DEBUG("No hard leading jet in " << recon_jets.size() << " jets");
 	vetoEvent;
       }
-      cerr << __LINE__ << "\n";
       for ( int i = 1; i < 3; ++i ) {
 	if ( recon_jets[i].momentum().pT() <= 25*GeV ) {
 	  vetoEvent;
 	}
       }
-      cerr << __LINE__ << "\n";
 
       for ( int i = 0; i < 3; ++i ) {
 	double dPhi = deltaPhi( pTmiss_phi, recon_jets[i].momentum().phi() );
@@ -252,7 +238,6 @@
 	  break;
 	}
       }
-      cerr << __LINE__ << "\n";
 
 
       ParticleVector lepton;
@@ -266,7 +251,6 @@
         foreach ( const Particle & e, recon_e )
 	    lepton.push_back(e);
       }
-      cerr << __LINE__ << "\n";
 
       std::sort(lepton.begin(), lepton.end(), cmpParticleByPt);
 
@@ -282,7 +266,6 @@
                 lepton[0].momentum().pT() <= 20*GeV ) {
 	vetoEvent;
       }
-      cerr << __LINE__ << "\n";
 
       // exactly one hard leading lepton cut
       if ( fabs(lepton[1].pdgId()) == e_id &&
@@ -294,7 +277,6 @@
 	  vetoEvent;
       }
 
-      cerr << __LINE__ << "\n";
 
 
     // ==================== FILL ====================
@@ -313,7 +295,6 @@
 	+ recon_jets[1].momentum().pT()
 	+ recon_jets[2].momentum().pT();
 
-      cerr << __LINE__ << "\n";
 
       // Electron channel signal region
 
@@ -345,7 +326,6 @@
         }
 
       }
-      cerr << __LINE__ << "\n";
 
 
     }