Commit 881626f4 authored by Andrzej Olszewski's avatar Andrzej Olszewski Committed by Graeme Stewart
Browse files

Hijing_i-00-01-03 (Hijing_i-00-01-03)

	* src/Hijing.cxx
	        Add filling HeavyIon component in GenEvent with collision parameters
	* Tagged as Hijing_i-00-01-03

2015-04-09  David Quarrie <David.Quarrie@cern.ch>
	* cmt/requirements
		Add cmake-specific dependency against Generators/GeneratorFortranCommon (transparent to CMT)
	* Tagged as Hijing_i-00-01-02
parent 90f47dc2
......@@ -15,6 +15,9 @@ use GeneratorUtils GeneratorUtils-* Generators
use TruthUtils TruthUtils-* Generators
use StoreGate StoreGate-* Control
# Add cmake-specific dependency (transparent to CMT)
apply_pattern cmake_add_dependency package=Generators/GeneratorFortranCommon
#public
# new to store hijing event parameters
# ------------------------------------
......
......@@ -265,6 +265,39 @@ Hijing::fillEvt(HepMC::GenEvent* evt)
// Set the generator id
evt->set_signal_process_id(HIJING + m_iap);
// Store collision parameters
int np = m_himain1.np();
int nt = m_himain1.nt();
int n0 = m_himain1.n0();
int n01 = m_himain1.n01();
int n10 = m_himain1.n10();
int n11 = m_himain1.n11();
//int natt = m_himain1.natt();
int jatt = m_himain1.jatt();
float b = m_hiparnt.hint1(19);
float bphi = m_hiparnt.hint1(20);
float sigmainel = m_hiparnt.hint1(12);
HepMC::HeavyIon ion
(
static_cast<int>(jatt), // Ncoll_hard
static_cast<int>(np), // Npart_proj
static_cast<int>(nt), // Npart_targ
static_cast<int>(n0+n10+n01+n11), // Ncoll
static_cast<int>(-1), // spectator_neutrons
static_cast<int>(-1), // spectator_protons
static_cast<int>(n01), // N_Nwounded_collisions
static_cast<int>(n10), // Nwounded_N_collisions
static_cast<int>(n11), // Nwounded_Nwounded_collisions
b, // impact_parameter
bphi, // event_plane_angle
-1, // eccentricity
sigmainel ); // sigma_inel_NN
evt->set_heavy_ion(ion);
std::cout << " heavy ion " << evt->heavy_ion() << std::endl;
// Did we keep decay history?
//
bool keptHistory = (m_hiparnt.ihpr2(21) == 1);
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment