get_file_list.pl -keys runnumber,events,filename -delim '/' -limit 0 -cond \
  'storage=hpss,sname2~st_laser,filetype=online_daq,runnumber>16000000' >! listdb1
sed 's/_adc//' listdb1 | cut -d '_' -f 3 | sort -n -u >! listdb2
cat listdb2 | colrm 1 2 | colrm 4 | sed 's/^0//' >! listdb3

mysql -h dbx.star.bnl.gov -P3316 -C Calibrations_tpc -N -e \
 'select dayofyear(addtime(beginTime,"-5:00:00")),count(1) from tpcDriftVelocity where beginTime>"2015-01-01" and deactive=0 group by dayofyear(addtime(beginTime,"-5:00:00"));' \
  >! listfc

root -l

gStyle->SetPalette(1);gStyle->SetGridColor(kGray);gStyle->SetOptDate(0);gStyle->SetOptStat(0);

TNtuple db("db","db","day:count"); // database
db.ReadFile("listfc");
TNtuple fc("fc","fc","day"); // filecatalog
fc.ReadFile("listdb3");

db.Draw("day>>aa1(80,39.5,119.5)","count");
fc.Draw("day>>aa2(80,39.5,119.5)");
aa1->SetLineColor(2);
aa2->SetXTitle("day");
aa2->SetTitle("Laser runs per day (red have entries in database)");
aa2->Draw();
aa1->Draw("same");
c1->SaveAs("Run15_LaserRunsPerDay.png");

db.Draw("day>>aa3(80,39.5,119.5)","count");
paa2 = aa2;
aa3->Divide(paa2);
aa3->SetXTitle("day");
aa3->SetTitle("Fraction of laser runs with drift velocities in database");
aa3->Draw();
c1->SaveAs("Run15_LaserRunsFraction.png");