Shower.h

// This file collects top-level classes for shower evolution. 

#ifndef Shower_H
#define Shower_H

#include <iostream>
#include <vector>
#include <cmath>
#include <ctime>
#include "Pythia7/Config/Pythia7.h"

namespace Pythia7 {
namespace Shower {

using namespace std;

class TimeShower;
class SpaceShower;

//**************************************************************************


class Event {
    
public:
  Event(long capacity = 100) {entry.reserve(capacity); maxColIndx = 0;}


  long append(Particle entryin) {    
    entry.push_back(entryin); 
    if (entryin.col() > maxColIndx) maxColIndx = entryin.col();   
    if (entryin.anticol() > maxColIndx) maxColIndx = entryin.anticol();
    return entry.size() - 1;
  }


  void remove(long index) {entry.erase(entry.begin() + index);}


  Particle& operator[](long i) {return entry[i];}
  const Particle & operator[](long i) const {return entry[i];}
  bool isDescendent(long, long) const;

  long size() const {return entry.size();}
  void zero() {entry.resize(0); maxColIndx = 0;}
  void colIndx(long indx) {maxColIndx = indx;}
  long colIndx() const {return maxColIndx;}


  friend ostream& operator<<(ostream&, const Event&) ;  


private: 
  vector<Particle> entry;
  long maxColIndx;
};
 
//**************************************************************************


class Shower {

public:
  static long ISR;

  static long FSR;

  static long FSRONISR;

  static TimeShower * theTimeShower;

  static SpaceShower * theSpaceShower;

  void shower(Event&, BeamParticle&, BeamParticle&, bool = false);
  void shower(Event&); 


private: 

  long in1Old, in2Old, in1New, in2New, sizeOld,sizeMid, sizeNew;
  bool hadISR;
  void doSpaceShower(Event&, BeamParticle&, BeamParticle&);
  void doTimeOnSpaceShower(Event&);
  void copyFinalState(Event&);
  void doTimeShower(Event&);
};
 
//**************************************************************************

}
}

#endif // Shower_H

---