#if !defined(__CLING__) || defined(__ROOTCLING__)

#include <iostream>

#include "Rtypes.h"

#include "TSystem.h"

#include "TMath.h"

#include "TString.h"

#include "TStopwatch.h"

#include "TGeoManager.h"

#include "FairRunSim.h"

#include "FairRuntimeDb.h"

#include "FairPrimaryGenerator.h"

#include "FairBoxGenerator.h"

#include "FairParRootFileIo.h"

#include "FairSystemInfo.h"

#include "Field/MagneticField.h"

#include "DetectorsPassive/Cave.h"

#include "TPCSimulation/Detector.h"

#include "Generators/Pythia8Generator.h"

#endif

void run_sim_pythia8hi(Int_t nEvents = 10, TString mcEngine = "TGeant3")

{

TString dir = getenv("VMCWORKDIR");

TString geom_dir = dir + "/Detectors/Geometry/";

gSystem->Setenv("GEOMPATH",geom_dir.Data());

TString tut_configdir = dir + "/Detectors/gconfig";

gSystem->Setenv("CONFIG_DIR",tut_configdir.Data());

// Output file name

char fileout[100];

sprintf(fileout, "AliceO2_%s.aligen.mc_%i_event.root", mcEngine.Data(), nEvents);

TString outFile = fileout;

// Parameter file name

char filepar[100];

sprintf(filepar, "AliceO2_%s.aligen.params_%i.root", mcEngine.Data(), nEvents);

TString parFile = filepar;

// In general, the following parts need not be touched

// Debug option

gDebug = 0;

// Timer

TStopwatch timer;

timer.Start();

// Create simulation run

FairRunSim* run = new FairRunSim();

// enable usage of the fair link mechanism

run->SetUseFairLinks(kTRUE);

run->SetName(mcEngine); // Transport engine

run->SetOutputFile(outFile); // Output file

FairRuntimeDb* rtdb = run->GetRuntimeDb();

// Create media

run->SetMaterials("media.geo"); // Materials

// Create geometry

o2::passive::Cave* cave = new o2::passive::Cave("CAVE");

cave->SetGeometryFileName("cave.geo");

run->AddModule(cave);

auto magField = std::make_unique<o2::field::MagneticField>("Maps","Maps", -1., -1., o2::field::MagFieldParam::k5kG);

run->SetField(magField.get());

// ===| Add TPC |============================================================

o2::tpc::Detector* tpc = new o2::tpc::Detector(kTRUE);

tpc->SetGeoFileName("TPCGeometry.root");

run->AddModule(tpc);

// Create TGenerator interface

auto gen = new o2::eventgen::Pythia8Generator();

gen->SetParameters("Beams:idA 1000822080"); // Pb ion

gen->SetParameters("Beams:idB 1000822080"); // Pb ion

gen->SetParameters("Beams:eCM 5520.0"); // [GeV]

gen->SetParameters("HeavyIon:SigFitNGen 0"); // valid for Pb-Pb 5520 only

gen->SetParameters("HeavyIon:SigFitDefPar 14.82,1.82,0.25,0.0,0.0,0.0,0.0,0.0"); // valid for Pb-Pb 5520 only

gen->SetParameters("HeavyIon:bWidth 1"); // impact parameter from 0-x [fm]

// Create PrimaryGenerator

auto primGen = new FairPrimaryGenerator();

primGen->AddGenerator(gen);

run->SetGenerator(primGen);

// store track trajectories

// run->SetStoreTraj(kTRUE);

// Initialize simulation run

run->Init();

// Runtime database

Bool_t kParameterMerged = kTRUE;

FairParRootFileIo* parOut = new FairParRootFileIo(kParameterMerged);

parOut->open(parFile.Data());

rtdb->setOutput(parOut);

rtdb->saveOutput();

rtdb->print();

// Start run

run->Run(nEvents);

delete run;

// run->CreateGeometryFile("geofile_full.root");

// Finish

timer.Stop();

Double_t rtime = timer.RealTime();

Double_t ctime = timer.CpuTime();

// extract max memory usage

FairSystemInfo sysinfo;

std::cout << std::endl << std::endl;

std::cout << "Macro finished succesfully." << std::endl;

std::cout << "Output file is " << outFile << std::endl;

std::cout << "Parameter file is " << parFile << std::endl;

std::cout << "Real time " << rtime << " s, CPU time " << ctime << "s" << std::endl << std::endl;

std::cout << "Memory used " << sysinfo.GetMaxMemory() << "\n";

}