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

blackhole at projects.hepforge.org blackhole at projects.hepforge.org
Wed Sep 23 16:43:08 BST 2009


Author: holsch
Date: Wed Sep 23 16:43:08 2009
New Revision: 1829

Log:
Adding N, ptsum vs dphi histos. The normalisation is wrong, please check the filling methods.

Modified:
   trunk/src/Analyses/CDF_2001_S4751469.cc

Modified: trunk/src/Analyses/CDF_2001_S4751469.cc
==============================================================================
--- trunk/src/Analyses/CDF_2001_S4751469.cc	Wed Sep 23 15:55:49 2009	(r1828)
+++ trunk/src/Analyses/CDF_2001_S4751469.cc	Wed Sep 23 16:43:08 2009	(r1829)
@@ -54,6 +54,21 @@
     
     // Book histograms
     void init() {
+      // These histos bin N, pt in dphi
+      _hist_num_dphi_2  = bookProfile1D("dummy_num_2", 50, 0, 180);
+      _hist_num_dphi_5  = bookProfile1D("dummy_num_5", 50, 0, 180);
+      _hist_num_dphi_30 = bookProfile1D("dummy_num_30", 50, 0, 180);
+      _hist_pt_dphi_2   = bookProfile1D("dummy_pt_2", 50, 0, 180);
+      _hist_pt_dphi_5   = bookProfile1D("dummy_pt_5", 50, 0, 180);
+      _hist_pt_dphi_30  = bookProfile1D("dummy_pt_30", 50, 0, 180);
+      
+      _numvsDeltaPhi2 =  bookProfile1D(1, 1, 1);  
+      _numvsDeltaPhi5 =  bookProfile1D(1, 1, 2);  
+      _numvsDeltaPhi30 = bookProfile1D(1, 1, 3);  
+      _pTvsDeltaPhi2 =   bookProfile1D(2, 1, 1);  
+      _pTvsDeltaPhi5 =   bookProfile1D(2, 1, 2);  
+      _pTvsDeltaPhi30 =  bookProfile1D(2, 1, 3);
+      
       _numTowardMB = bookProfile1D(3, 1, 1);
       _numTransMB = bookProfile1D(3, 1, 2);
       _numAwayMB = bookProfile1D(3, 1, 3);
@@ -111,6 +126,14 @@
       // Run over tracks
       double ptSumToward(0.0), ptSumAway(0.0), ptSumTrans(0.0);
       size_t numToward(0), numTrans(0), numAway(0);
+      // Reset the histos that bin N, pt in dphi
+      _hist_num_dphi_2->reset();
+      _hist_num_dphi_5->reset();
+      _hist_num_dphi_30->reset();
+      _hist_pt_dphi_2->reset();
+      _hist_pt_dphi_5->reset();
+      _hist_pt_dphi_30->reset();
+
       foreach (const Jet& j, jets) {
         foreach (const FourMomentum& p, j) {
           // Calculate Delta(phi) from leading jet
@@ -145,6 +168,29 @@
             ++numAway;
           }
           
+          // Fill histos to bin pt, N in dphi
+          // dphi in degrees
+          const double dPhideg = 360*dPhi/TWOPI;
+          //
+          if (ptLead/GeV > 2.0) {
+            _hist_num_dphi_2->fill(dPhideg, 1, weight         );
+            _hist_pt_dphi_2->fill (dPhideg, p.pT()/GeV, weight);
+
+          }
+          if (ptLead/GeV > 5.0) { 
+            _hist_num_dphi_5->fill(dPhideg, 1, weight         );
+            _hist_pt_dphi_5->fill (dPhideg, p.pT()/GeV, weight);
+
+          }
+          if (ptLead/GeV > 30.0) {
+            _hist_num_dphi_30->fill(dPhideg, 1, weight         );
+            _hist_pt_dphi_30->fill (dPhideg, p.pT()/GeV, weight);
+          }
+
+
+
+
+
         }
       }
       
@@ -158,6 +204,23 @@
                << endl;
       
       // Update the pT profile histograms
+      //
+      // N, sumpt vs. dphi first
+      // TODO: normalisation
+      for (int i= 0; i < 50; i++) {
+        // pT > 2 GeV  
+        _numvsDeltaPhi2->fill(_hist_num_dphi_2->binMean(i), _hist_num_dphi_2->binHeight(i));///_sumWeightsPtLead2);
+        _pTvsDeltaPhi2->fill(_hist_pt_dphi_2->binMean(i), _hist_pt_dphi_2->binHeight(i));///_sumWeightsPtLead2);
+        
+        // pT > 5 GeV  
+        _numvsDeltaPhi5->fill(_hist_num_dphi_5->binMean(i),_hist_num_dphi_5->binHeight(i));///_sumWeightsPtLead5);
+        _pTvsDeltaPhi5->fill(_hist_pt_dphi_5->binMean(i), _hist_pt_dphi_5->binHeight(i));///_sumWeightsPtLead5);
+        
+        // pT > 30 GeV  
+        _numvsDeltaPhi30->fill(_hist_num_dphi_30->binMean(i),_hist_num_dphi_30->binHeight(i));///_sumWeightsPtLead30);
+        _pTvsDeltaPhi30->fill(_hist_pt_dphi_30->binMean(i), _hist_pt_dphi_30->binHeight(i));///_sumWeightsPtLead30);
+      }
+      
       _ptsumTowardMB->fill(ptLead/GeV, ptSumToward/GeV, weight);
       _ptsumTowardJ20->fill(ptLead/GeV, ptSumToward/GeV, weight);
       
@@ -208,6 +271,16 @@
 
     /// @name Histogram collections
     //@{
+    // These histos (binned in dphi) are filled per event and then reset
+    // TODO: use LWH
+    AIDA::IProfile1D *_hist_num_dphi_2, *_hist_num_dphi_5, *_hist_num_dphi_30;
+    AIDA::IProfile1D *_hist_pt_dphi_2, *_hist_pt_dphi_5, *_hist_pt_dphi_30;
+
+    // The sumpt vs. dphi and Nch vs. dphi histos
+    AIDA::IProfile1D *_numvsDeltaPhi2, *_numvsDeltaPhi5, *_numvsDeltaPhi30;
+    AIDA::IProfile1D *_pTvsDeltaPhi2, *_pTvsDeltaPhi5, *_pTvsDeltaPhi30;
+
+
     /// Profile histograms, binned in the \f$ p_T \f$ of the leading jet, for
     /// the \f$ p_T \f$ sum in the toward, transverse and away regions.
     AIDA::IProfile1D *_ptsumTowardMB,  *_ptsumTransMB,  *_ptsumAwayMB;


More information about the Rivet-svn mailing list