fMomentum(0.),
fBeta(0.),
fPmax (0.),
- fATarget (0.),
- fZTarget (0.),
fCharge(0),
- fProtonDirection(0.),
+ fProtonDirection(1.),
fTemperatureG(0.),
fBetaSourceG(0.),
fTemperatureB(0.),
fMomentum(0.),
fBeta(0.),
fPmax (10.),
- fATarget (208.),
- fZTarget (82.),
fCharge(1),
- fProtonDirection(0.),
- fTemperatureG(0.04),
+ fProtonDirection(1.),
+ fTemperatureG(0.05),
fBetaSourceG(0.05),
- fTemperatureB(0.004),
+ fTemperatureB(0.005),
fBetaSourceB(0.),
fNgp(0),
fNgn(0),
// Initialization
//
Float_t kMass = TDatabasePDG::Instance()->GetParticle(kProton)->Mass();
- fMomentum = fCMS/2. * fZTarget / fATarget;
+ fMomentum = fCMS/2. * Float_t(fZTarget) / Float_t(fATarget);
fBeta = fMomentum / TMath::Sqrt(kMass * kMass + fMomentum * fMomentum);
if (fDebug) {
fDebugHist1 = new TH2F("DebugHist1", "nu vs N_slow", 100, 0., 100., 20, 0., 20.);
//
if (fCollisionGeometry) {
Float_t b = fCollisionGeometry->ImpactParameter();
- Int_t nn = fCollisionGeometry->NN();
- Int_t nwn = fCollisionGeometry->NwN();
- Int_t nnw = fCollisionGeometry->NNw();
- Int_t nwnw = fCollisionGeometry->NwNw();
+ // Int_t nn = fCollisionGeometry->NN();
+ // Int_t nwn = fCollisionGeometry->NwN();
+ // Int_t nnw = fCollisionGeometry->NNw();
+ // Int_t nwnw = fCollisionGeometry->NwNw();
fSlowNucleonModel->GetNumberOfSlowNucleons(fCollisionGeometry, fNgp, fNgn, fNbp, fNbn);
if (fDebug) {
- printf("Nucleons %d %d %d %d \n", fNgp, fNgn, fNbp, fNbn);
+ printf("Slow nucleons: %d grayp %d grayn %d blackp %d blackn \n", fNgp, fNgn, fNbp, fNbn);
fDebugHist1->Fill(Float_t(fNgp + fNgn + fNbp + fNbn), fCollisionGeometry->NwN(), 1.);
fDebugHist2->Fill(Float_t(fNgp + fNgn + fNbp + fNbn), b, 1.);
- printf("AliGenSlowNucleons: Impact parameter from Collision Geometry %f %d %d %d %d\n",
- b, nn, nwn, nnw, nwnw);
+ //printf("AliGenSlowNucleons: Impact parameter from Collision Geometry %f %d %d %d %d\n", b, nn, nwn, nnw, nwnw);
}
}
Three-momentum [GeV/c] is given back in q[3]
*/
+ //printf("Generating slow nuc. with: charge %d. temp. %1.3f, beta %1.2f\n",charge,T,beta);
+
Double_t m, pmax, p, f, phi;
TDatabasePDG * pdg = TDatabasePDG::Instance();
const Double_t kMassProton = pdg->GetParticle(kProton) ->Mass();
{
/* Lorentz transform in the direction of q[2] */
- Double_t gamma = 1/sqrt((1.-beta)*(1.+beta));
+ Double_t gamma = 1./sqrt((1.-beta)*(1.+beta));
Double_t energy = sqrt(m*m + q[0]*q[0] + q[1]*q[1] + q[2]*q[2]);
q[2] = gamma * (q[2] + beta*energy);
}
+