Changeset 14153 for trunk


Ignore:
Timestamp:
06/11/12 10:54:24 (12 years ago)
Author:
tbretz
Message:
Replaced my own homogeneous number rgenerator with gRandom->Circle
File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/Mars/msimreflector/MSimRays.cc

    r10054 r14153  
    164164        MPhotonData &dat = *static_cast<MPhotonData*>(arr.UncheckedAt(idx));
    165165
     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
    166174        // 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  */
    173184        // The is the incident direction of the photon
    174185        // h==0 means infinitiy
Note: See TracChangeset for help on using the changeset viewer.