53 using namespace Jetscape;
60 int main(
int argc,
char** argv)
62 clock_t
t; t = clock();
69 JetScapeLogger::Instance()->SetInfo(
true);
70 JetScapeLogger::Instance()->SetDebug(
false);
71 JetScapeLogger::Instance()->SetRemark(
false);
74 JetScapeLogger::Instance()->SetVerboseLevel(0);
79 auto jetscape = make_shared<JetScape>(argv[1],10000);
80 jetscape->SetId(
"primary");
83 auto trento = make_shared<TrentoInitial>();
84 auto pythiaGun= make_shared<PythiaGun> ();
85 auto hydro = make_shared<Brick> ();
86 jetscape->SetReuseHydro(
false);
87 jetscape->Add(trento);
88 jetscape->Add(pythiaGun);
93 auto jlossmanager = make_shared<JetEnergyLossManager> ();
94 auto jloss = make_shared<JetEnergyLoss> ();
96 auto matter = make_shared<Matter> ();
108 jlossmanager->Add(jloss);
109 jetscape->Add(jlossmanager);
111 auto printer = make_shared<PartonPrinter> () ;
112 jetscape->Add(printer);
115 auto hadroMgr = make_shared<HadronizationManager> ();
116 auto hadro = make_shared<Hadronization> ();
119 auto colorless = make_shared<ColorlessHadronization> ();
120 hadro->Add(colorless);
121 hadroMgr->Add(hadro);
122 jetscape->Add(hadroMgr);
126 auto writer= make_shared<JetScapeWriterAscii> (argv[2]);
153 printf (
"CPU time: %f seconds.\n",((
float)t)/CLOCKS_PER_SEC);
154 printf (
"Real time: %f seconds.\n",difftime(end,start));
160 Pythia8::Info& info = pythiaGun->info;
167 ofstream WriteSigmaHard;
169 WriteSigmaHard<<info.sigmaGen()<<
"\t"<<info.sigmaErr()<<endl;
170 WriteSigmaHard.close();
179 INFO_NICE<<
"------------------------------------";
180 INFO_NICE<<
"| Brick Test JetScape Framework ... |";
181 INFO_NICE<<
"------------------------------------";