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"

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;

      EventOutputConfig() :
        printEnergy(false), printMotherPdg(false), printPhotonsLeptons(false), printDecayEpoch(false) { }
    };

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

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

} // namespace thermalfist

#endif