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

blackhole at projects.hepforge.org blackhole at projects.hepforge.org
Tue May 18 06:52:05 BST 2010


Author: buckley
Date: Tue May 18 10:10:20 2010
New Revision: 2452

Log:
Removing unused sections rather than commenting

Modified:
   trunk/src/Analyses/CDF_1994_S2952106.cc

Modified: trunk/src/Analyses/CDF_1994_S2952106.cc
==============================================================================
--- trunk/src/Analyses/CDF_1994_S2952106.cc	Mon May 17 14:51:15 2010	(r2451)
+++ trunk/src/Analyses/CDF_1994_S2952106.cc	Tue May 18 10:10:20 2010	(r2452)
@@ -33,12 +33,9 @@
       // Veto neutrinos, and muons with pT above 1.0 GeV
       VetoedFinalState vfs(fs);
       vfs.vetoNeutrinos();
-      //VisibleFinalState vfs(fs);
       vfs.addVetoPairDetail(MUON, 1.0*GeV, MAXDOUBLE);
       addProjection(vfs, "VFS");
       addProjection(FastJets(vfs, FastJets::CDFJETCLU, 0.7), "ConeJets");
-      //addProjection(MissingMomentum(vfs), "CalMET");
-
 
       /// @todo Use histogram auto-booking
 
@@ -119,29 +116,17 @@
       const Jets jets = applyProjection<FastJets>(event, "ConeJets").jetsByPt();
       getLog() << Log::DEBUG << "Jet multiplicity before any cuts = " << jets.size() << endl;
 
-      // Check there isn't too much missing Et
-      ///const MissingMomentum& caloMissEt = applyProjection<MissingMomentum>(event, "CalMET");
-      ///getLog() << Log::DEBUG << "Missing ET = " << caloMissEt.visibleMomentum().Et()/GeV << " GeV" << endl;
-      /// @todo Looks dodgy to me... only difference is pT vs. Et. Really?
-      //if ((caloMissEt.visibleMomentum().pT()/GeV)/sqrt(caloMissEt.visibleMomentum().Et()/GeV) > 6.0) vetoEvent;
-      //if ((caloMissEt.visibleMomentum().ET()/GeV)/sqrt(caloMissEt.scalarET()/GeV) > 6.0) vetoEvent;
-      //if ((caloMissEt.visibleMomentum().Et()/GeV)/sqrt(caloMissEt.momentum().Et()/GeV) > 6.0) vetoEvent;
-
-
-      //Et's only from jets:
+      // ETs only from jets:
       _Et_sinphi = 0.;
       _Et_cosphi = 0.;
       _Et = 0.;
       for (int i=0; i< jets.size(); ++i) {
-	_Et_sinphi = jets[i].momentum().Et() * sin(jets[i].phi());
-	_Et_cosphi = jets[i].momentum().Et() * sin(jets[i].phi());
-	_Et = jets[i].momentum().Et() * sin(jets[i].phi());
-
+        _Et_sinphi = jets[i].momentum().Et() * sin(jets[i].phi());
+        _Et_cosphi = jets[i].momentum().Et() * sin(jets[i].phi());
+        _Et = jets[i].momentum().Et() * sin(jets[i].phi());
       }
       if (sqrt(_Et_sinphi*_Et_sinphi + _Et_cosphi*_Et_cosphi)/_Et > 6.0) vetoEvent;
 
-
-
       // Check jet requirements
       if (jets.size() < 3) vetoEvent;
       if (jets[0].momentum().pT() < 100*GeV) vetoEvent;
@@ -245,7 +230,7 @@
 
       vector<double> xval_R23, xerr_R23, yval_R23, yerr_R23;
       for (int ibin = 0;  ibin < 35; ++ibin) { // x-value + weight=err^2/y-value
-	xval_R23.push_back((ibin+0.5)*4.375/35.);
+        xval_R23.push_back((ibin+0.5)*4.375/35.);
         xerr_R23.push_back(0.5*4.375/35.);
         yval_R23.push_back(R23_CDF_sim[ibin]/R23_Ideal_sim[ibin]);
         yerr_R23.push_back(R23_CDF_sim_err[ibin]/R23_Ideal_sim[ibin]);


More information about the Rivet-svn mailing list