{
// Generate one event
// --------------------------------------------------------------------------
- fSpecn = 0;
- fSpecp = 0;
+ fProjectileSpecn = 0;
+ fProjectileSpecp = 0;
+ fTargetSpecn = 0;
+ fTargetSpecp = 0;
// --------------------------------------------------------------------------
fHijing->GenerateEvent();
fTrials++;
// --------------------------------------------------------------------------
// Count spectator neutrons and protons
- if(ks == 0 || ks == 1 || ks == 10 || ks == 11){
- if(kf == kNeutron) fSpecn += 1;
- if(kf == kProton) fSpecp += 1;
+ if(ks == 0 || ks == 1){
+ if(kf == kNeutron) fProjectileSpecn += 1;
+ if(kf == kProton) fProjectileSpecp += 1;
+ }
+ else if(ks == 10 || ks == 11){
+ if(kf == kNeutron) fTargetSpecn += 1;
+ if(kf == kProton) fTargetSpecp += 1;
}
// --------------------------------------------------------------------------
//
fHijing->GetN01(),
fHijing->GetN10(),
fHijing->GetN11());
- ((AliGenHijingEventHeader*) header)->SetSpectators(fSpecn, fSpecp);
+ ((AliGenHijingEventHeader*) header)->SetSpectators(fProjectileSpecn, fProjectileSpecp,
+ fTargetSpecn,fTargetSpecp);
// 4-momentum vectors of the triggered jets.
//
Int_t fNoGammas; // Don't write gammas if flag "on"
// ZDC proposal (by Chiara) to store num. of SPECTATORS protons and neutrons
- Int_t fSpecn; // Num. of spectator neutrons
- Int_t fSpecp; // Num. of spectator protons
+ Int_t fProjectileSpecn;// Num. of spectator neutrons from projectile nucleus
+ Int_t fProjectileSpecp;// Num. of spectator protons from projectile nucleus
+ Int_t fTargetSpecn; // Num. of spectator neutrons from target nucleus
+ Int_t fTargetSpecp; // Num. of spectator protons from target nucleus
Int_t fLHC; // Assume LHC as lab frame
private:
void Copy(TObject &rhs) const;