Changeset 17622 for trunk/ippScripts/scripts/ipp_simulation_data.pl
- Timestamp:
- May 9, 2008, 3:28:25 PM (18 years ago)
- File:
-
- 1 edited
-
trunk/ippScripts/scripts/ipp_simulation_data.pl (modified) (9 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/ippScripts/scripts/ipp_simulation_data.pl
r16489 r17622 13 13 use Pod::Usage qw( pod2usage ); 14 14 use PS::IPP::Config qw( 15 16 17 18 19 20 21 22 23 15 $PS_EXIT_SUCCESS 16 $PS_EXIT_UNKNOWN_ERROR 17 $PS_EXIT_SYS_ERROR 18 $PS_EXIT_CONFIG_ERROR 19 $PS_EXIT_PROG_ERROR 20 $PS_EXIT_DATA_ERROR 21 $PS_EXIT_TIMEOUT_ERROR 22 caturi 23 ); 24 24 25 25 my $ipprc = PS::IPP::Config->new(); # IPP configuration 26 26 27 my ($name, # Base name for output images 28 $camera, # Name of camera to use 29 $telescope, # Telescope name 30 $dbname, # Database name 31 $path, # Path to data 32 $workdir, # Working directory for data 33 $no_cal, # Don't produce calibration files 34 $no_update # Don't update the database 27 my ($name, # Base name for output images 28 $camera, # Name of camera to use 29 $telescope, # Telescope name 30 $tess_id, # Tessellation identifier 31 $dbname, # Database name 32 $path, # Path to data 33 $workdir, # Working directory for data 34 $no_cal, # Don't produce calibration files 35 $no_update # Don't update the database 35 36 ); 36 37 37 38 GetOptions( 38 'name=s' => \$name, 39 'camera=s' => \$camera, 40 'telescope=s' => \$telescope, 41 'dbname=s' => \$dbname, 42 'path=s' => \$path, 43 'workdir=s' => \$workdir, 44 'no-cal' => \$no_cal, 45 'no-update' => \$no_update, 46 ) or pod2usage( 2 ); 39 'name=s' => \$name, 40 'camera=s' => \$camera, 41 'telescope=s' => \$telescope, 42 'tess_id=s' => \$tess_id, 43 'dbname=s' => \$dbname, 44 'path=s' => \$path, 45 'workdir=s' => \$workdir, 46 'no-cal' => \$no_cal, 47 'no-update' => \$no_update, 48 ) or pod2usage( 2 ); 47 49 48 50 pod2usage( -msg => "Unknown option: @ARGV", -exitval => 2 ) if @ARGV; … … 50 52 pod2usage( 51 53 -msg => "Required options: --name --path --camera --telescope --dbname", 52 53 54 -exitval => 3, 55 ) unless 54 56 defined $name and 55 57 defined $path and … … 67 69 or (warn "Can't find ipp_serial_inject_mosaic.pl" and $missing_tools = 1); 68 70 69 if ($missing_tools) { 71 if ($missing_tools) { 70 72 warn ("Can't find required tools"); 71 exit($PS_EXIT_CONFIG_ERROR); 73 exit($PS_EXIT_CONFIG_ERROR); 72 74 } 73 75 … … 76 78 # Dark exposure times 77 79 use constant DARK => [ 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 78 80 300, 300, 300, 300, 300, 300, 300, 300, 300, 300 ]; 79 81 # Flat-field image characteristics 80 82 use constant FLAT => [ 81 82 83 84 85 86 87 88 89 83 { 84 filter => 'r', 85 exptime => [ 0.1, 0.1, 0.1, 0.1, 0.1, 0.5, 1, 2, 5, 10, 20, 20, 20, 20, 20 ], 86 }, 87 { 88 filter => 'i', 89 exptime => [ 20, 20, 20, 20, 20 ], 90 } 91 ]; 90 92 # Object image characteristics 91 93 use constant OBJECT => [ 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 use constant SCALE => 0.2; 94 { 95 filter => 'r', 96 exptime => [ 5, 10, 10, 10, 10, 10, 240 ], 97 seeing => [ 0.8, 0.4, 0.6, 0.8, 1.2, 1.5, 0.7 ], 98 ra => 150.119167, 99 dec => 2.205833, 100 pa => 0, 101 zp => 25.15, 102 sky => 20.86, 103 dither => 40, 104 }, 105 { 106 filter => 'i', 107 exptime => [ 5, 30, 30, 30, 30, 30, 240 ], 108 seeing => [ 0.8, 0.4, 0.6, 0.8, 1.2, 1.5, 0.7 ], 109 ra => 150.119167, 110 dec => 2.205833, 111 pa => 0, 112 zp => 25.00, 113 sky => 20.15, 114 dither => 40, 115 }, 116 ]; 117 use constant SCALE => 0.2; # Plate scale 116 118 117 119 ############################################################################################################# … … 120 122 121 123 my $counter = 0; 124 my $tessellation = (defined $tess_id) ? "--tess_id $tess_id" : ""; # Additional switch for tessellation 122 125 123 126 # Generate bias images 124 127 for (my $i = 0; $i < BIAS; $i++) { 125 my $basename; 128 my $basename; # Output base filename 126 129 ( $basename, $counter ) = filename( $name, $counter ); 127 130 my $filename = caturi( $path, $basename ); 128 131 unless ($no_cal) { 129 130 131 132 133 "--workdir $workdir --dbname $dbname $basename",134 135 132 run( command => "$ppSim -camera $camera -type BIAS $filename", 133 verbose => 1 ) or die "Unable to run ppSim"; 134 unless ($no_update) { 135 run( command => "$inject --camera $camera --telescope $telescope --path $path " . 136 "$tessellation --workdir $workdir --dbname $dbname $basename", 137 verbose => 1 ) or die "Unable to inject file."; 138 } 136 139 } 137 140 } … … 139 142 # Generate dark images 140 143 foreach my $exptime ( @{DARK()} ) { 141 my $basename; 144 my $basename; # Output base filename 142 145 ( $basename, $counter ) = filename( $name, $counter ); 143 146 my $filename = caturi( $path, $basename ); 144 147 unless ($no_cal) { 145 146 147 148 149 "--workdir $workdir --dbname $dbname $basename",150 151 148 run ( command => "$ppSim -camera $camera -type DARK -exptime $exptime $filename", 149 verbose => 1 ) or die "Unable to run ppSim"; 150 unless ($no_update) { 151 run( command => "$inject --camera $camera --telescope $telescope --path $path " . 152 "$tessellation --workdir $workdir --dbname $dbname $basename", 153 verbose => 1 ) or die "Unable to inject file."; 154 } 152 155 } 153 156 } … … 157 160 my $filter = $set->{filter}; # Name of filter 158 161 foreach my $exptime ( @{$set->{exptime}} ) { 159 my $basename;# Output base filename160 161 162 163 164 165 166 167 "--workdir $workdir --dbname $dbname $basename",168 169 170 162 my $basename; # Output base filename 163 ( $basename, $counter ) = filename( $name, $counter ); 164 my $filename = caturi( $path, $basename ); 165 unless ($no_cal) { 166 run( command => "$ppSim -camera $camera -type FLAT -filter $filter -exptime $exptime $filename", 167 verbose => 1 ) or die "Unable to run ppSim"; 168 unless ($no_update) { 169 run( command => "$inject --camera $camera --telescope $telescope --path $path " . 170 "$tessellation --workdir $workdir --dbname $dbname $basename", 171 verbose => 1 ) or die "Unable to inject file."; 172 } 173 } 171 174 } 172 175 } … … 175 178 foreach my $set ( @{OBJECT()} ) { 176 179 my $filter = $set->{filter}; # Name of filter 177 my $ra0 = $set->{ra}; 178 my $dec0 = $set->{dec}; 179 my $pa = $set->{pa}; 180 my $zp = $set->{zp}; 181 my $scale = SCALE(); 182 my $sky = 10**( -0.4 * ( $set->{sky} - $zp ) ) * $scale**2; 183 my $dither = $set->{dither} / 3600; 180 my $ra0 = $set->{ra}; # Base Right Ascension (deg) 181 my $dec0 = $set->{dec}; # Base Declination (deg) 182 my $pa = $set->{pa}; # Position angle (deg) 183 my $zp = $set->{zp}; # Zero point 184 my $scale = SCALE(); # Plate scale (arcsec/pix) 185 my $sky = 10**( -0.4 * ( $set->{sky} - $zp ) ) * $scale**2; # Sky background (counts/s) 186 my $dither = $set->{dither} / 3600; # Dither size (deg) 184 187 185 188 for (my $i = 0; $i < scalar @{$set->{exptime}}; $i++) { 186 187 188 189 190 191 my $basename;# Output base filename192 193 194 195 196 197 198 199 "--workdir $workdir --dbname $dbname $basename",200 201 189 my $exptime = ${$set->{exptime}}[$i]; # Exposure time 190 my $seeing = ${$set->{seeing}}[$i]; # Seeing (pix) 191 my $ra = $ra0 + (2*rand() - 1) * $dither * cos(deg2rad($dec0)); # RA with dither 192 my $dec = $dec0 + (2*rand() - 1) * $dither; # Dec with dither 193 194 my $basename; # Output base filename 195 ( $basename, $counter ) = filename( $name, $counter ); 196 my $filename = caturi( $path, $basename ); 197 run( command => "$ppSim -camera $camera -type OBJECT -filter $filter -exptime $exptime " . 198 "-skyrate $sky -ra $ra -dec $dec -pa $pa -scale $scale -zp $zp -seeing $seeing $filename", 199 verbose => 1 ) or die "Unable to run ppSim"; 200 unless ($no_update) { 201 run( command => "$inject --camera $camera --telescope $telescope --path $path " . 202 "$tessellation --workdir $workdir --dbname $dbname $basename", 203 verbose => 1 ) or die "Unable to inject file."; 204 } 202 205 } 203 206 } … … 210 213 sub filename 211 214 { 212 my $base = shift; 213 my $num = shift; 214 my $workdir = shift; 215 my $base = shift; # Base name 216 my $num = shift; # Number 217 my $workdir = shift; # Working directory 215 218 my $name = sprintf("$base%04d", $num); 216 219 $num++;
Note:
See TracChangeset
for help on using the changeset viewer.
