SimpleEvent.h

Go to the documentation of this file.

/*
 * Thermal-FIST package
 * 
 * Copyright (c) 2014-2018 Volodymyr Vovchenko
 *
 * GNU General Public License (GPLv3 or later)
 */
#ifndef SIMPLEEVENT_H
#define SIMPLEEVENT_H
 
#include <cmath>
#include <vector>
#include <fstream>
 
#include "HRGEventGenerator/SimpleParticle.h"
#include "HRGBase/ThermalParticleSystem.h"
 
namespace thermalfist {
  struct SimpleEvent {
    double weight;

    double logweight;

    std::vector<SimpleParticle> Particles;

    std::vector<SimpleParticle> AllParticles;

    std::vector<SimpleParticle> PhotonsLeptons;

    std::vector<int> DecayMap;

    std::vector<int> DecayMapFinal;

    SimpleEvent() { Particles.resize(0); AllParticles.resize(0); PhotonsLeptons.resize(0); weight = 1.; logweight = 0.; }

    void RapidityBoost(double dY);

    static SimpleEvent MergeEvents(const SimpleEvent &evt1, const SimpleEvent &evt2);

    struct EventOutputConfig {

      bool printEnergy;

      bool printMotherPdg;

      bool printPhotonsLeptons;

      bool printDecayEpoch;

      bool printCoordinates;

      bool printWeight;
 
      EventOutputConfig() :
        printEnergy(true), printMotherPdg(false), printPhotonsLeptons(false), printDecayEpoch(false), printCoordinates(false), printWeight(true) { }
    };

    void writeToFile(std::ofstream& fout, const EventOutputConfig& config = EventOutputConfig(), int eventnumber = 1) const;

    void writeToFile(std::ofstream& fout, int eventnumber = 1) const { writeToFile(fout, EventOutputConfig(), eventnumber); }

    void writeToFileForUrqmd(std::ofstream& fout) const;

    void writeToFileForSmash(std::ofstream& fout, const thermalfist::ThermalParticleSystem* TPS = NULL, const int eventnumber = 1) const;
  };
 
} // namespace thermalfist
 
#endif