Changeset 40939 for branches/eam_branches/ipp-20191011/ppSim/src/ppSimMakeGalaxies.c

Timestamp:

Oct 13, 2019, 1:00:44 PM (7 years ago)

Author:

eugene

Message:

add option to insert a gaussian-profile globular cluster

File:

• branches/eam_branches/ipp-20191011/ppSim/src/ppSimMakeGalaxies.c (modified) (3 diffs)

branches/eam_branches/ipp-20191011/ppSim/src/ppSimMakeGalaxies.c

@@ -45 +45 @@
     float zp              = psMetadataLookupF32(&mdok, recipe, "ZEROPOINT"); // Photometric zero point
     float seeing          = psMetadataLookupF32(&mdok, recipe, "SEEING"); // Seeing sigma (pix)
-    float scale           = psMetadataLookupF32(&mdok, recipe, "PIXEL.SCALE") * M_PI / 3600.0 / 180.0; // Plate scale (radians/pixel)
+    float scale           = psMetadataLookupF32(&mdok, recipe, "PIXEL.SCALE"); // Plate scale (arcsec/pixel)
     float skyRate         = psMetadataLookupF32(&mdok, recipe, "SKY.RATE"); // Sky rate
     if (isnan(skyRate)) {
         float skyMags = psMetadataLookupF32(&mdok, recipe, "SKY.MAGS");  assert (mdok);
         skyRate = scale * scale * ppSimMagToFlux (skyMags, zp);
+        // skyMags is in mags / square arcsec so scale must be arcsec / pixel
     }
 
@@ -94 +95 @@
 
     // Normalisation, set by the specified stellar density at the specified bright magnitude
-    float refSum = galaxyDensity * xSize * ySize * PS_SQR(scale * 180.0 / M_PI);
+    float refSum = galaxyDensity * xSize * ySize * PS_SQR(scale / 3600.0);
     float normLum =  refSum / (galaxyAlpha * logf(10.0) * powf(10.0, (galaxyAlpha * brightMag)));
     float normScale = normLum * galaxyAlpha * logf(10.0);
@@ -143 +144 @@
                 rndValue = galaxyGridRandom ? drand48() : i / (float) nTotal;
                 float scale = (galaxyRmajorMin + rndValue * galaxyRmajorSlope);
-
                 galaxy->Rmaj  = scale;