]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - PYTHIA8/pythia8170/examples/main12.cc
PYTHIA8: removing legacy pythia8170
[u/mrichter/AliRoot.git] / PYTHIA8 / pythia8170 / examples / main12.cc
diff --git a/PYTHIA8/pythia8170/examples/main12.cc b/PYTHIA8/pythia8170/examples/main12.cc
deleted file mode 100644 (file)
index b749d17..0000000
+++ /dev/null
@@ -1,111 +0,0 @@
-// main12.cc is a part of the PYTHIA event generator.
-// Copyright (C) 2012 Torbjorn Sjostrand.
-// PYTHIA is licenced under the GNU GPL version 2, see COPYING for details.
-// Please respect the MCnet Guidelines, see GUIDELINES for details.
-
-// This is a simple test program. 
-// It illustrates how Les Houches Event File input can be used in PYTHIA.
-// It uses two LHE files, ttbar.lhe and ttbar2.lhe, which are combined
-// using Beams:newLHEFsameInit = on to skip new initialization second time.
-// Then the second file is viewed as a simple continuation of the first,
-// just split for practical reasons, rather than as a separate new run 
-// with a new set of processes.
-// In the first file top decays have been performed, in the second not,
-// and are instead handled by the internal PYTHIA resonance-decay machinery. 
-// Furthermore the internal top production processes are switched on and
-// mixed in, giving an unrealistic "double up" total top cross section.
-// Much of this of course is not intended to be realistic, 
-// but rather illustrates several tricks that can be useful.
-
-#include "Pythia.h"
-using namespace Pythia8; 
-int main() {
-
-  //  Number of listed events. Allow for possibility of a few faulty events.
-  int nPrintLHA  = 1;             
-  int nPrintRest = 0;             
-  int nAbort     = 10;
-
-  // Generator           
-  Pythia pythia;                            
-
-  // Switch on internal ttbar production.
-  pythia.readString("Top:gg2ttbar = on");  
-  pythia.readString("Top:qqbar2ttbar = on"); 
-
-  // Use same top mass as in Pythia 6.4 to simplify comparison.
-  pythia.readString("6:m0 = 175.");    
-
-  // No automatic event listings - do it manually below.
-  pythia.readString("Next:numberShowLHA = 0"); 
-  pythia.readString("Next:numberShowInfo = 0"); 
-  pythia.readString("Next:numberShowProcess = 0"); 
-  pythia.readString("Next:numberShowEvent = 0"); 
-
-  // Initialize Les Houches Event File run.
-  pythia.readString("Beams:frameType = 4"); 
-  pythia.readString("Beams:LHEF = ttbar.lhe");   
-  pythia.init();   
-
-  // Book histogram.
-  Hist nCharged("charged particle multiplicity",100,-0.5,399.5); 
-
-  // Set counters.
-  int iPrintLHA  = 0;             
-  int iPrintRest = 0;             
-  int iAbort     = 0;
-  int iFile      = 1;
-
-  // Begin event loop   
-  for (int iEvent = 0; ; ++iEvent) {
-
-    // Generate until none left in input file.  
-    if (!pythia.next()) {
-      if (pythia.info.atEndOfFile()) {
-
-        // First time open next file, second time stop event loop.
-        if (iFile == 1) {
-          pythia.readString("Beams:newLHEFsameInit = on"); 
-          pythia.readString("Beams:LHEF = ttbar2.lhe");   
-          pythia.init();   
-          ++iFile;
-          continue;
-       } else break;          
-      }
-
-      // First few failures write off as "acceptable" errors, then quit.
-      if (++iAbort < nAbort) continue;
-      break;
-    }
-  
-    // List first few Les Houches and other events.
-    if (pythia.info.isLHA() && iPrintLHA < nPrintLHA) {     
-      pythia.LHAeventList();               
-      pythia.info.list();          
-      pythia.process.list();          
-      pythia.event.list();  
-      ++iPrintLHA;         
-    } else if (!pythia.info.isLHA() && iPrintRest < nPrintRest) {     
-      pythia.info.list();          
-      pythia.process.list();          
-      pythia.event.list();           
-      ++iPrintRest;         
-    }                 
-
-    // Sum up final charged multiplicity and fill in histogram.
-    int nChg = 0;                 
-    for (int i = 0; i < pythia.event.size(); ++i) 
-    if (pythia.event[i].isFinal() && pythia.event[i].isCharged()) 
-      ++nChg;
-    nCharged.fill(nChg);               
-
-  // End of event loop.        
-  }                                           
-
-  // Give statistics. Print histogram.
-  pythia.stat();
-  cout << nCharged;  
-
-  // Done.                           
-  return 0;
-}