<& /htmlheader, title => 'Themis TJD ' . ($nn + 1 - 50000) &> <& modelheader &>

Themis Family Asteroids

MOPS Night Number <% $nn %> (TJD <% $nn + 1 - 50000 %>)
<% mopslib_mjd2utctimestr($nn) %> to <% mopslib_mjd2utctimestr($nn + 1) %>

<& /lib/pager, rooturi => "?nn=$nn&minv=$minv&min_dets=$min_dets&sm=$sm&min_prob=$min_prob", limit => $limit, offset => $offset, found => $found &> % foreach $href (@rows) { % my $trkId = $href->{tracklet_id}; % my $trk = modct_retrieve($inst, trackletId => $trkId); % }
Tracklet ID digest2 (MOPS)
Sort 		Position Vtot (deg/day)
Sort 		Pos Ang (deg)
Sort 		Known As Known q
Sort 		GCR (arcsec)
Sort 		Probability
Sort 		Stamps
<% $trkId %>
MPCheck
Digest MOPS| MPC| OORB| All
Detections MPC | DES 		<% sprintf "nfuzz: %d", $href->{nfuzz} %>
<% sprintf "ffuzz: %.1f", $href->{ffuzz} %>
Stamp Info 		<% $href->{survey_mode} %>
α: <% sprintf "%.6f", $trk->extRA %>
<% turn2str(deg2turn($trk->extRA), 'H', 2, 'hms') %>		 δ: <% sprintf "%.6f", $trk->extDEC %>
<% turn2str(deg2turn($trk->extDEC), 'D', 2, 'deg') %>
 
S Elong: <% sprintf "%.2f", $href->{sol_elong} %> Ecl β: <% sprintf "%+.2f", $href->{ecl_lat} %>
<% sprintf "%.3f", $trk->vTot %> <% sprintf "%.1f", $trk->posAng %> <& /mpc_link, desig => $href->{known_name} &> {q} && $href->{q} <= 1.3) ? 'style="background:red"' : '' %> align="center"> <% $href->{q} ? sprintf("%.2f", $href->{q}) : 'N/A' %> <% sprintf "%.2f"", $href->{gcr_arcsec} %>
 <% sprintf "%.2f", $href->{probability} %> <& stamp4c_td, dets => $trk->detections, decider => 1, short => 1, ffuzzmap => $href->{ffuzzmap} &>
<& /lib/pager, rooturi => "?nn=$nn&minv=$minv&min_dets=$min_dets&sm=$sm&min_prob=$min_prob", limit => $limit, offset => $offset, found => $found &>
<& /htmlfooter &> <%args> $nn $limit => 100 $offset => 0 $found => undef $min_dets => 2 $min_prob => 0 $minv => 0 $sm => '' <%init> use Astro::Time; use PS::MOPS::DC::Orbit; use PS::MOPS::Lib qw(:all); use PS::MOPS::Constants qw(:all); use PS::MOPS::DC::Field; use PS::MOPS::DC::Tracklet; use PS::MOPS::DC::Known; # If $offset and $inc are defined, use them for the mysql query. # Otherwise use $limit, then as a last resort, all of em. my $limit_str = ''; my $found_str = ''; if (defined($offset)) { $limit = 100 unless $limit; $limit_str = sprintf 'limit %d,%d', $offset, $limit; } else { $limit_str = $limit ? "limit $limit" : ''; } my $prob_str = ''; if ($min_prob) { if ($min_prob < 0) { $prob_str = "and t.probability < " . (-$min_prob); } else { $prob_str = "and t.probability >= $min_prob"; } } if (!defined($found)) { $found_str = 'SQL_CALC_FOUND_ROWS'; } my $sm_str = ''; if ($sm) { $sm_str = "and survey_mode like '%$sm%'"; } my $inst = $m->notes('INST'); my $model = $m->notes('MODEL'); my $dbh = $inst->dbh(); my $sql = <<"SQL"; select $found_str t.tracklet_id tracklet_id, t.gcr_arcsec gcr_arcsec, t.probability probability, sum(1) ndet, group_concat(concat(det_id,'=',round(dr.moments_r1*dr.moments_r1/(dr.psf_major*dr.psf_minor),2))) ffuzzmap, sum(if(dr.psf_qf_perfect*dr.ratio_all*dr.f_pos > 0.5 and dr.moments_r1*dr.moments_r1/(dr.psf_major*dr.psf_minor) < 20,1,0)) nfuzz, /* Median of comet fuzz moment ratios */ ( substring_index( substring_index( group_concat( dr.moments_r1*dr.moments_r1/(dr.psf_major*dr.psf_minor) order by dr.moments_r1*dr.moments_r1/(dr.psf_major*dr.psf_minor) ) , ',' , ceiling(count(*)/2) ) , ',' , -1 ) + substring_index( substring_index( group_concat( dr.moments_r1*dr.moments_r1/(dr.psf_major*dr.psf_minor) order by dr.moments_r1*dr.moments_r1/(dr.psf_major*dr.psf_minor) ) , ',' , - ceiling(count(*)/2) ) , ',' , 1 ) ) / 2 as ffuzz, k.q q, k.known_name known_name, f.se_deg sol_elong, f.eb_deg ecl_lat, f.survey_mode from fields f join tracklets t using(field_id) join known k using(known_id) join custom.themis using (known_name) join tracklet_attrib ta using(tracklet_id) join detections d using(det_id) join det_rawattr_v2 dr using(det_id) where f.nn = $nn and t.classification='N' and t.v_tot > $minv $prob_str $sm_str group by tracklet_id having nfuzz > ndet / 2 and ndet >= $min_dets order by ffuzz desc $limit_str SQL my $sth = $dbh->prepare($sql) or die $dbh->errstr; $sth->execute() or die $sth->errstr; my @rows; my $href; while ($href = $sth->fetchrow_hashref()) { push @rows, $href; } $sth->finish(); if ($found_str) { ($found) = $dbh->selectrow_array("select FOUND_ROWS()"); }