69 Float_t t, cost, sint, cosp, sinp;
72 const Float_t g2rad = TMath::Pi() / 180.;
80 }
else if (
mMode == 1) {
81 Float_t theta = 2. * TMath::ATan(TMath::Exp(-c1));
82 sint = TMath::Sin(theta);
83 cost = TMath::Cos(theta);
85 cost = TMath::Cos(c1 * g2rad);
86 sint = TMath::Sin(c1 * g2rad);
89 cosp = TMath::Cos(c2 * g2rad);
90 sinp = TMath::Sin(c2 * g2rad);
92 pmom[0] = cosp * sint;
93 pmom[1] = sinp * sint;
97 orig[0] = orig[1] = orig[2] = 0;
104 orig[0] = pmom[0] * t;
105 orig[1] = pmom[1] * t;
106 orig[2] = pmom[2] * t;
107 if (TMath::Abs(orig[2]) >
mZMax) {
111 Float_t polar[3] = {0., 0., 0.};
112 primGen->AddTrack(0, pmom[0], pmom[1], pmom[2], orig[0], orig[1], orig[2], -1,
true, 0, 0, 1.);