- Timestamp:
- 06/11/12 10:54:24 (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/Mars/msimreflector/MSimRays.cc
r10054 r14153 164 164 MPhotonData &dat = *static_cast<MPhotonData*>(arr.UncheckedAt(idx)); 165 165 166 Double_t x, y; 167 const Double_t r = gRandom->Uniform(); 168 gRandom->Circle(x, y, maxr*TMath::Sqrt(r)); 169 /* 170 Double_t ra = gRandom->Uniform(maxr); 171 Double_t ph = gRandom->Uniform(TMath::TwoPi()); 172 173 166 174 // Get radom incident point on the mirror plane. 167 const Double_t x = gRandom->Uniform(-maxr, maxr); 168 const Double_t y = gRandom->Uniform(-maxr, maxr); 169 170 if (x*x + y*y > maxr*maxr) 171 continue; 172 175 //const Double_t x = gRandom->Uniform(-maxr, maxr); 176 //const Double_t y = gRandom->Uniform(-maxr, maxr); 177 178 Double_t x = ra*sin(ph); 179 Double_t y = ra*cos(ph); 180 181 // if (x*x + y*y > maxr*maxr) 182 // continue; 183 */ 173 184 // The is the incident direction of the photon 174 185 // h==0 means infinitiy
Note:
See TracChangeset
for help on using the changeset viewer.