src/framework/PreequilibriumDynamics.cc

0001 /*******************************************************************************
0002  * Copyright (c) The JETSCAPE Collaboration, 2018
0003  *
0004  * Modular, task-based framework for simulating all aspects of heavy-ion collisions
0005  *
0006  * For the list of contributors see AUTHORS.
0007  *
0008  * Report issues at https://github.com/JETSCAPE/JETSCAPE/issues
0009  *
0010  * or via email to bugs.jetscape@gmail.com
0011  *
0012  * Distributed under the GNU General Public License 3.0 (GPLv3 or later).
0013  * See COPYING for details.
0014  ******************************************************************************/
0015
0016 #include "PreequilibriumDynamics.h"
0017 #include "JetScapeLogger.h"
0018 #include "JetScapeXML.h"
0019 #include "JetScapeSignalManager.h"
0020 #include <string>
0021
0022 #include <iostream>
0023
0024 using namespace std;
0025
0026 #define MAGENTA "\033[35m"
0027
0028 namespace Jetscape {
0029
0030 PreequilibriumDynamics::PreequilibriumDynamics() {
0031   VERBOSE(8);
0032   SetId("PreequilibriumDynamics");
0033 }
0034
0035 PreequilibriumDynamics::~PreequilibriumDynamics() {
0036   VERBOSE(8);
0037   disconnect_all();
0038 }
0039
0040 void PreequilibriumDynamics::Init() {
0041   JetScapeModuleBase::Init();
0042
0043   JSINFO << "Initialize PreequilibriumDynamics : " << GetId() << " ...";
0044
0045   VERBOSE(8);
0046
0047   // this is grabbing the initial entropy density ?
0048   ini = JetScapeSignalManager::Instance()->GetInitialStatePointer().lock();
0049   if (!ini) {
0050     JSWARN << "No initialization module, try: "
0051            << "auto trento = make_shared<TrentoInitial>(); "
0052            << "jetscape->Add(trento);";
0053   }
0054
0055   InitializePreequilibrium(parameter_list_);
0056
0057   InitTask();
0058
0059   JetScapeTask::InitTasks();
0060 }
0061
0062 void PreequilibriumDynamics::Exec() {
0063   VERBOSE(2) << "Run Preequilibrium : " << GetId() << " ...";
0064   VERBOSE(8) << "Current Event #" << GetCurrentEvent();
0065
0066   if (ini) {
0067     VERBOSE(3) << "length of entropy density vector="
0068                << ini->GetEntropyDensityDistribution().size();
0069   }
0070
0071   EvolvePreequilibrium();
0072
0073   JetScapeTask::ExecuteTasks();
0074 }
0075
0076 void PreequilibriumDynamics::Clear() {
0077   e_.clear();
0078   P_.clear();
0079   utau_.clear();
0080   ux_.clear();
0081   uy_.clear();
0082   ueta_.clear();
0083   pi00_.clear();
0084   pi01_.clear();
0085   pi02_.clear();
0086   pi03_.clear();
0087   pi11_.clear();
0088   pi12_.clear();
0089   pi13_.clear();
0090   pi22_.clear();
0091   pi23_.clear();
0092   pi33_.clear();
0093   bulk_Pi_.clear();
0094 }
0095
0096 } // end namespace Jetscape