sPhenix code displayed by LXR master/​JETSCAPE/​examples/​custom_examples/​freestream-milneTest.cc	Source navigation Diff markup Identifier search General search
	Version: master Revert   Change

File indexing completed on 2025-08-05 08:18:30

 /*******************************************************************************
  * Copyright (c) The JETSCAPE Collaboration, 2018
  *
  * Modular, task-based framework for simulating all aspects of heavy-ion collisions
  * 
  * For the list of contributors see AUTHORS.
  *
  * Report issues at https://github.com/JETSCAPE/JETSCAPE/issues
  *
  * or via email to bugs.jetscape@gmail.com
  *
  * Distributed under the GNU General Public License 3.0 (GPLv3 or later).
  * See COPYING for details.
  ******************************************************************************/
 
 #include <iostream>
 #include <time.h>
 
 // JetScape Framework includes ...
 #include "JetScape.h"
 #include "JetEnergyLoss.h"
 #include "JetEnergyLossManager.h"
 #include "JetScapeWriterStream.h"
 #ifdef USE_HEPMC
 #include "JetScapeWriterHepMC.h"
 #endif
 
 //User modules derived from jetscape framework clasess
 //to be used to run Jetscape ...
 #include "AdSCFT.h"
 #include "Matter.h"
 #include "Martini.h"
 #include "FreestreamMilneWrapper.h"
 #include "MusicWrapper.h"
 #include "iSpectraSamplerWrapper.h"
 #include "TrentoInitial.h"
 #include "PGun.h"
 #include "PythiaGun.h"
 //#include "HadronizationManager.h"
 //#include "Hadronization.h"
 //#include "ColoredHadronization.h"
 
 #include <chrono>
 #include <thread>
 
 using namespace std;
 
 using namespace Jetscape;
 
 // Forward declaration
 void Show();
 
 // -------------------------------------
 
 int main(int argc, char** argv)
 {
   clock_t t; t = clock();
   time_t start, end; time(&start);
 
   cout<<endl;
 
   //DEBUG=true by default and REMARK=false
   //can be also set also via XML file (at least partially)
   JetScapeLogger::Instance()->SetInfo(true);
   JetScapeLogger::Instance()->SetDebug(true);
   JetScapeLogger::Instance()->SetRemark(false);
   //SetVerboseLevel (9 a lot of additional debug output ...)
   //If you want to suppress it: use SetVerboseLevel(0) or max  SetVerboseLevel(9) or 10
   JetScapeLogger::Instance()->SetVerboseLevel(0);
 
   Show();
 
   auto jetscape = make_shared<JetScape>();
   const char* mainXMLName = "../config/jetscape_main.xml";
   const char* userXMLName = "../config/jetscape_user.xml";
 
   jetscape->SetXMLMainFileName(mainXMLName);
   jetscape->SetXMLUserFileName(userXMLName);
 
   jetscape->SetNumberOfEvents(1);
   jetscape->SetReuseHydro (false);
   jetscape->SetNReuseHydro (0);
   // jetscape->SetNumberOfEvents(10);
   //jetscape->SetReuseHydro (true);
   //jetscape->SetNReuseHydro (5);
 
   auto trento = make_shared<TrentoInitial> ();
   auto freestream = make_shared<FreestreamMilneWrapper> ();
   auto hydro = make_shared<MpiMusic> ();
   //auto hydro = make_shared<GubserHydro> ();
   
   // surface sampler
   auto iSS = make_shared<iSpectraSamplerWrapper> ();
 
   auto jlossmanager = make_shared<JetEnergyLossManager> ();
   auto jloss = make_shared<JetEnergyLoss> ();
   auto matter = make_shared<Matter> ();
   //auto martini = make_shared<Martini> ();
   //auto adscft = make_shared<AdSCFT> ();
   //DBEUG: Remark:
   //does not matter unfortunately since not called recursively, done by JetEnergyLoss class ...
   //matter->SetActive(false);
   //martini->SetActive(false);
   // This works ... (check with above logic ...)
   //jloss->SetActive(false);
 
   // auto pGun= make_shared<PGun> ();
   auto pythiaGun= make_shared<PythiaGun> ();
 
   //auto hadroMgr = make_shared<HadronizationManager> ();
   //auto hadro = make_shared<Hadronization> ();
   //auto hadroModule = make_shared<ColoredHadronization> ();
 
   // only pure Ascii writer implemented and working with graph output ...
   auto writer= make_shared<JetScapeWriterAscii> ("test_out.dat");
   //auto writer= make_shared<JetScapeWriterAsciiGZ> ("test_out.dat.gz");
 #ifdef USE_HEPMC
   //auto writer= make_shared<JetScapeWriterHepMC> ("test_out.dat");
 #endif
   //writer->SetActive(false);
 
   //Remark: For now modules have to be added
   //in proper "workflow" order (can be defined via xml and sorted if necessary)
 
   jetscape->Add(trento);
 
   // jetscape->Add(pGun);
   jetscape->Add(pythiaGun);
 
   jetscape->Add(freestream);
 
   //Some modifications will be needed for reusing hydro events, so far
   //simple test hydros always executed "on the fly" ...
   jetscape->Add(hydro);
   jetscape->Add(iSS);
 
   // Matter with silly "toy shower (no physics)
   // and Martini dummy ...
   // Switching Q2 (or whatever variable used
   // hardcoded at 5 to be changed to xml)
   jloss->Add(matter);
   //jloss->Add(martini);
   //jloss->Add(adscft);
 
   jlossmanager->Add(jloss);
 
   jetscape->Add(jlossmanager);
 
 
   //hadro->Add(hadroModule);
   //hadroMgr->Add(hadro);
   //jetscape->Add(hadroMgr);
 
 
 
   jetscape->Add(writer);
 
   // Initialize all modules tasks
   jetscape->Init();
 
   // Run JetScape with all task/modules as specified ...
   jetscape->Exec();
 
   // "dummy" so far ...
   // Most thinkgs done in write and clear ...
   jetscape->Finish();
 
   INFO_NICE<<"Finished!";
   cout<<endl;
 
   // wait for 5s
   //std::this_thread::sleep_for(std::chrono::milliseconds(500000));
 
   t = clock() - t;
   time(&end);
   printf ("CPU time: %f seconds.\n",((float)t)/CLOCKS_PER_SEC);
   printf ("Real time: %f seconds.\n",difftime(end,start));
   //printf ("Real time: %f seconds.\n",(start-end));
   return 0;
 }
 
 // -------------------------------------
 
 void Show()
 {
   INFO_NICE<<"-----------------------------------------------";
   INFO_NICE<<"| freestream-milne Test JetScape Framework ... |";
   INFO_NICE<<"-----------------------------------------------";
   INFO_NICE;
 }

This page was automatically generated by the 2.3.7 LXR engine. The LXR team