[Rivet-svn] r1988 - in trunk: data/anainfo src/Analyses

[blackhole at projects.hepforge.org]

Author: fsiegert
Date: Tue Nov  3 16:49:32 2009
New Revision: 1988

Log:
Minor finishing modifications to ALEPH photon fragmentation analysis.
Setting to VALIDATED, because Pythia looks ok in a first look.

Modified:
   trunk/data/anainfo/ALEPH_1996_S3196992.info
   trunk/src/Analyses/ALEPH_1996_S3196992.cc

Modified: trunk/data/anainfo/ALEPH_1996_S3196992.info
==============================================================================
--- trunk/data/anainfo/ALEPH_1996_S3196992.info	Tue Nov  3 15:38:41 2009	(r1987)
+++ trunk/data/anainfo/ALEPH_1996_S3196992.info	Tue Nov  3 16:49:32 2009	(r1988)
@@ -4,15 +4,14 @@
 Experiment: ALEPH
 Collider: LEP Run 1
 SpiresID: 3196992
-Status: UNVALIDATED
+Status: VALIDATED
 Authors:
  - Frank Siegert <frank.siegert at durham.ac.uk>
 References:
  - Z.Phys.C69:365-378,1996
  - doi:10.1007/s002880050037
 RunInfo:
-  e+ e- -> jets at shower level. The experimental results have been corrected
-  for hadronisation
+  e+ e- -> jets with hadron decays turned off.
 NumEvents: 1000000
 PtCuts: [0]
 Description:

Modified: trunk/src/Analyses/ALEPH_1996_S3196992.cc
==============================================================================
--- trunk/src/Analyses/ALEPH_1996_S3196992.cc	Tue Nov  3 15:38:41 2009	(r1987)
+++ trunk/src/Analyses/ALEPH_1996_S3196992.cc	Tue Nov  3 16:49:32 2009	(r1988)
@@ -25,9 +25,9 @@
 
     void init() {
       // Set up projections
-      FinalState fs;      
+      FinalState fs;
       addProjection(FastJets(fs, FastJets::DURHAM, 0.7), "DurhamJets");
-      IdentifiedFinalState ifs(-MAXRAPIDITY, +MAXRAPIDITY, 5.0*GeV);
+      IdentifiedFinalState ifs(-MAXRAPIDITY, +MAXRAPIDITY, 0.0);
       ifs.acceptId(PHOTON);
       addProjection(ifs, "Photons");
       addProjection(Thrust(fs), "Thrust");
@@ -51,7 +51,7 @@
       const ParticleVector allphotons = applyProjection<IdentifiedFinalState>(event, "Photons").particles();
       ParticleVector photons;
       foreach (const Particle& photon, allphotons) {
-        if (fabs(cos(photon.momentum().theta()))<0.95) {
+        if (fabs(cos(photon.momentum().theta()))<0.95 && photon.momentum().E()>5.0*GeV) {
           photons.push_back(photon);
         }
       }
@@ -74,7 +74,7 @@
             double zgamma = photon.momentum().E()/jet.E();
             if (jets_001.size() == 2) _h_z_2jet_001->fill(zgamma, weight);
             else if (jets_001.size() == 3) _h_z_3jet_001->fill(zgamma, weight);
-            else if (jets_001.size() == 4) _h_z_4jet_001->fill(zgamma, weight);
+            else if (jets_001.size() > 3) _h_z_4jet_001->fill(zgamma, weight);
             break;
           }
         }