Logo coherent WaveBurst  
Library Reference Guide
Logo
TestCreateFrList.C
Go to the documentation of this file.
1 //
2 // Test Create FrList
3 // Author : Gabriele Vedovato
4 
5 #define SORT_FILE_NAME "root/GHLTV-HBRST14_S6D_R1-Sorted.root"
6 #define CHANNEL_NAME "L1:GW-H"
7 
8 {
9 
10  int nDQF=12;
11  dqfile DQF[12]={
12  {0 ,"SEGMENTS/S6D_R10_segments/S6D_OFFLINE_L1SCIENCE.txt" , 0, 0., false, false},
13  {0 ,"SEGMENTS/S6D_R10_segments/S6D_OFFLINE_L1_DQCAT1SEGMENTS.txt", 1, 0., true , false},
14  {0 ,"SEGMENTS/S6D_R10_segments/S6D_OFFLINE_L1_DQCAT2SEGMENTS.txt", 2, 0., true , false},
15  {0 ,"SEGMENTS/S6D_R10_segments/S6D_OFFLINE_L1_DQCAT4SEGMENTS.txt", 1, 0., true , false},
16  {1 ,"SEGMENTS/S6D_R10_segments/S6D_OFFLINE_H1SCIENCE.txt" , 0, 0., false, false},
17  {1 ,"SEGMENTS/S6D_R10_segments/S6D_OFFLINE_H1_DQCAT1SEGMENTS.txt", 1, 0., true , false},
18  {1 ,"SEGMENTS/S6D_R10_segments/S6D_OFFLINE_H1_DQCAT2SEGMENTS.txt", 2, 0., true , false},
19  {1 ,"SEGMENTS/S6D_R10_segments/S6D_OFFLINE_H1_DQCAT4SEGMENTS.txt", 1, 0., true , false},
20  {2 ,"SEGMENTS/S6D_R10_segments/S6D_OFFLINE_V1SCIENCE.txt" , 0, 0., false, false},
21  {2 ,"SEGMENTS/S6D_R10_segments/S6D_OFFLINE_V1_DQCAT1SEGMENTS.txt", 1, 0., true , false},
22  {2 ,"SEGMENTS/S6D_R10_segments/S6D_OFFLINE_V1_DQCAT2SEGMENTS.txt", 2, 0., true , false},
23  {2 ,"SEGMENTS/S6D_R10_segments/S6D_OFFLINE_V1_DQCAT4SEGMENTS.txt", 1, 0., true , false},
24  };
25 
26  cwbtb tb;
27 
28  // create merged dq file list
29  int dqcat = 1;
30  std::vector<waveSegment> olist;
31  olist=tb.readSegList(nDQF, DQF, dqcat);
32  double ctime=tb.getSegTime(olist);
33  cout << "ctime : " << int(ctime) << " sec " << ctime/3600. << " h " << ctime/86400. << " day" << endl;
34 
35  // create job file list
36  std::vector<waveSegment> jlist;
37  jlist=tb.getJobList(olist, 300, 600, 8);
38 
40  int start,stop;
41  int jobId=0;
42  char ofile[256];
45  cout << "nJobs : " << jlist.size() << endl;
46  //for(int j=0;j<jlist.size();j++) {
47  //for(int j=4030;j<jlist.size();j++) {
48  for(int j=0;j<2;j++) {
49  start=jlist[j].start;
50  stop=jlist[j].stop;
51  jobId=j+1;
52  cout << "jobId : " << jobId << " " << "start " << start << " stop " << stop << endl;
53  frf=tb.getFrList(sort_file_name, start, stop, 8);
54  sprintf(ofile, "lists/HBRST14_S6D_2.lst.%d",jobId);
55  tb.dumpFrList(frf, ofile);
56  tb.readFrames(ofile, CHANNEL_NAME, w);
57  cout << "Data : " << w.size() << " " << w.rate() << endl;
58  }
59 
60  olist.clear();
61 
62  exit(0);
63 }
virtual void rate(double r)
Definition: wavearray.hh:141
TString("c")
int stop
std::vector< waveSegment > olist
char ofile[256]
int j
Definition: cwb_net.C:28
#define CHANNEL_NAME
cout<< "ctime : "<< int(ctime)<< " sec "<< ctime/3600.<< " h "<< ctime/86400.<< " day"<< endl;std::vector< waveSegment > jlist
nDQF
Definition: cwb_eced.C:109
virtual size_t size() const
Definition: wavearray.hh:145
int dqcat
i() int(T_cor *100))
frfile frf
#define SORT_FILE_NAME
wavearray< double > w
cwbtb tb
int start
sprintf(tfres,"(1/%g)x(%g) (sec)x(Hz)", 2 *df, df)
dqfile DQF[12]
int jobId
double ctime
TString sort_file_name
exit(0)