From: Christian Holm Christensen (cholm@hehi03.nbi.dk)
Date: Mon Jun 02 2003 - 10:55:56 EDT
Hi Steve, [changed to `dev' list] "Stephen J. Sanders" <ssanders@ku.edu> wrote concerning brhijing [Wed, 21 May 2003 16:35:23 -0500] ---------------------------------------------------------------------- > Hi, > I just checked out breg with the hope of running Hijing on > the rcas machines. However, the make fails miserably with > numerous invalid declarations and inconsistent common block > length errors for hijing_main.f The problem is, that Hiro introduced additional variables in the common block `eventinfo', which unfortunately has the same name as existing variables in Hijing itself. Hence, there's a conflict. I would suggest changing the variable names in `event.inc' from `breg/base'. > Is there a trick for building hijing on rcas? If you want to use plain Hijing, I would suggest getting `THijing' from my web-site [1]. It's a sub-class of `TGenerator', which makes it easy to use in ROOT. Then, you can output the relevant information to a plain text file: int testHijing(Int_t nEvents=10) { // Define your main function here gROOT->LoadClass("THijing", ".libs/libHijing.so"); // gDebug = 1; TClonesArray* particles = new TClonesArray("TParticle"); THijing* hijing = new THijing("hijing", "HIJING Simulation"); hijing->SetParameter("ihpr2",12,1); hijing->Initialize(200, "CMS", "A", "A", 197, 79, 197, 79); std::ostream output("data.dat"); if (!output) return 1; Int_t i; for (i = 0; i < nEvents; i++) { cout << "Event # " << i << " ... " << flush; Int_t npart = hijing->GenerateEvent(0, 0, "final", particles); output << npart << "\t" << i << std::endl; TIter next(particles); TParticle* particle = 0; while ((particle = (TParticle*)next())) { pid = 0; switch (particle->GetPdgCode()) { case 211: pid = ... ... } output << pid << "\t" << particle->Px() << "\t" << particle->Py() << "\t" << particle->Pz() << "\t" << std::endl; } } output.close(); return 0; } (see also the `testHijing.C' example that comes with the tar-ball.) Then, if you're concerned with the disk-space, you can `bzip2' the files after you're done with them (a cool thing to do with `.dat', and `.cdat' files). Yours, ___ | Christian Holm Christensen |_| | ------------------------------------------------------------- | | Address: Sankt Hansgade 23, 1. th. Phone: (+45) 35 35 96 91 _| DK-2200 Copenhagen N Cell: (+45) 24 61 85 91 _| Denmark Office: (+45) 353 25 305 ____| Email: cholm@nbi.dk Web: www.nbi.dk/~cholm | | [1] http://cholm.home.cern.ch/cholm/root/#thijing
This archive was generated by hypermail 2.1.5 : Mon Jun 02 2003 - 10:57:23 EDT