/*
$Log$
+Revision 1.1 2000/07/10 13:58:01 fca
+New version of ZDC from E.Scomparin & C.Oppedisano
+
Revision 1.7 2000/01/19 17:17:40 fca
Revision 1.6 1999/09/29 09:24:35 fca
//
// Generate one trigger (n or p)
//
+ Int_t i;
+
Double_t mass, pLab[3], balp0, balp[3], ddp[3], dddp0, dddp[3];
Float_t ptot = fPMin;
Int_t nt;
pLab[1] = 0.;
pLab[2] = ptot*TMath::Cos(scang);
}
- for(Int_t i=0; i<=2; i++){
+ for(i=0; i<=2; i++){
fP[i] = pLab[i];
}
if(fIpart==kProton) {mass = 0.93956563;}
if(fIpart==kNeutron) {mass = 0.93827231;}
// printf(" pLABx = %f pLABy = %f pLABz = %f \n",pLab[0],pLab[1],pLab[2]);
- for(Int_t i=0; i<=2; i++){
+ for(i=0; i<=2; i++){
balp[i] = -pLab[i];
}
balp0 = TMath::Sqrt(pLab[0]*pLab[0]+pLab[1]*pLab[1]+pLab[2]*pLab[2]+mass*mass);
- for(Int_t i=0; i<=2; i++){
+ for(i=0; i<=2; i++){
dddp[i] = ddp[i];
}
dddp0 = TMath::Sqrt(dddp[0]*dddp[0]+dddp[1]*dddp[1]+dddp[2]*dddp[2]+mass*mass);
// printf(" Boosted momentum -> px = %f, py = %f, pz = %f\n",
// pFermi[0], pFermi[1], pFermi[2]);
- for(Int_t i=0; i<=2; i++){
+ for(i=0; i<=2; i++){
fBoostP[i] = pFermi[i];
}
{
Double_t tetpart, fipart, tetdiv, fidiv, angleSum[2], tetsum, fisum, dplab[3];
Double_t rvec;
+
+ Int_t i;
Double_t pmq = 0.;
- for(int i=0; i<=2; i++){
+ for(i=0; i<=2; i++){
dplab[i] = pLab[i];
pmq = pmq+pLab[i]*pLab[i];
}
pLab[2] = pmod*TMath::Cos(tetsum);
// printf(" pLab[0] = %f pLab[1] = %f pLab[2] = %f \n\n",
// pLab[0],pLab[1],pLab[2]);
- for(Int_t i=0; i<=2; i++){
+ for(i=0; i<=2; i++){
fDivP[i] = pLab[i];
}
}
/*
$Log$
+Revision 1.1 2000/07/10 13:58:01 fca
+New version of ZDC from E.Scomparin & C.Oppedisano
+
Revision 1.7 2000/01/19 17:17:40 fca
Revision 1.6 1999/09/29 09:24:35 fca
//_____________________________________________________________________________
void AliZDCv1::InitTables()
{
+ Int_t k, j;
//Initialize parameters for light tables and read them
fNalfan = 90;
fNalfap = 90;
return;
}
// printf(" --- Reading light tables for ZN \n");
- for(int k=0; k<fNalfan; k++){
- for(int j=0; j<fNben; j++){
+ for(k=0; k<fNalfan; k++){
+ for(j=0; j<fNben; j++){
fscanf(fp1,"%f",&fTablen[0][k][j]);
fscanf(fp2,"%f",&fTablen[1][k][j]);
fscanf(fp3,"%f",&fTablen[2][k][j]);
return;
}
// printf(" --- Reading light tables for ZP \n");
- for(int k=0; k<fNalfap; k++){
- for(int j=0; j<fNbep; j++){
+ for(k=0; k<fNalfap; k++){
+ for(j=0; j<fNbep; j++){
fscanf(fp5,"%f",&fTablep[0][k][j]);
fscanf(fp6,"%f",&fTablep[1][k][j]);
fscanf(fp7,"%f",&fTablep[2][k][j]);
// Routine called at every step in the Zero Degree Calorimeters
//
+ Int_t j;
+
Int_t vol[2], ibeta, ialfa, ibe;
Float_t x[3], xdet[3], destep, hits[9], m, ekin, um[3], ud[3], be, radius, out;
TLorentzVector s, p;
//Particle coordinates
gMC->TrackPosition(s);
- for(Int_t j=0; j<=2; j++){
+ for(j=0; j<=2; j++){
x[j] = s[j];
}
hits[0] = x[0];
// Distance between particle trajectory and fibre axis
gMC->TrackPosition(s);
- for(Int_t j=0; j<=2; j++){
+ for(j=0; j<=2; j++){
x[j] = s[j];
}
gMC->Gmtod(x,xdet,1);