Changeset 38707

Timestamp:

Sep 3, 2015, 6:34:56 AM (11 years ago)

Author:

eugene

Message:

add some flexibility to 2mass

Location:

branches/eam_branches/ipp-20150625/Ohana/src/fakeastro

Files:

• include/fakeastro.h (modified) (1 diff)
• src/args.c (modified) (1 diff)
• src/make_2mass_measures.c (modified) (2 diffs)

branches/eam_branches/ipp-20150625/Ohana/src/fakeastro/include/fakeastro.h

@@ -103 +103 @@
 
 float  RADIUS;
+float  MAX_MAG_2MASS;
 
 SkyRegion UserPatch;

branches/eam_branches/ipp-20150625/Ohana/src/fakeastro/src/args.c

@@ -104 +104 @@
   if ((N = get_argument (*argc, argv, "-uniform-radec"))) {
     UNIFORM_RADEC = TRUE;
+    remove_argument (N, argc, argv);
+  }
+
+  MAX_MAG_2MASS = 16.0;
+  if ((N = get_argument (*argc, argv, "-2mass-limit"))) {
+    remove_argument (N, argc, argv);
+    MAX_MAG_2MASS = atof(argv[N]);
     remove_argument (N, argc, argv);
   }

branches/eam_branches/ipp-20150625/Ohana/src/fakeastro/src/make_2mass_measures.c

@@ -24 +24 @@
   // PhotCode *codeK = GetPhotcodebyName ("2MASS_K");
 
-  float ZP = 21.0; // this is chosen to give 16.0 mag stars a flux of 100.0 counts (crude, yes)
-  float SkyCts = 300.0; // this is chosen to make SN @ 16.0 = 5.0 (
+  float ZP = MAX_MAG_2MASS + 6.5; // this is chosen to give stars at the limit a flux of 400.0 counts (crude, yes)
+  float SkyCts = 1000.0; // this is chosen to make SN @ limit ~ 10.0 (
 
   int Nmeasure = catalog->Nmeasure;
@@ -45 +45 @@
     // make a crude JHK color for now:
     double J_PS1 = secfilt[i*Nsecfilt + Nsec].M - 1.0; // make all stars have y-J = 1.0
-    if (J_PS1 > 16.0) continue; // only generate 2MASS detections for objects with J_PS1 < 16.0
+    if (J_PS1 > MAX_MAG_2MASS) continue; // only generate 2MASS detections for objects with J_PS1 < 16.0
 
     if (isnan(secfilt[i*Nsecfilt + Nsec].M)) {