[u/mrichter/AliRoot.git] / ANALYSIS / TagMacros / ParseRCTDump.C

// Macro to parse the text dump of the RCT and convert it into
// a RunTag object. This object needs to be then merged with the
// proper RunTag from the Run.
//
// Parameters:
// Name of the text file with the dump from MonaLisa RCT page
//
// Author: Adam.Kisiel@cern.ch

int getposition(const char *token, int poscount)
{
  int position = -1;

  if (strstr(token, "#raw_run")) {
    position = 0;
    cout << "Run number at " << poscount << endl;
  }
  else if (strstr(token, "partition")) {
    position = 1;
    cout << "LHC Period at " << poscount << endl;
  }
  else if (strstr(token, "fillno")) {
    position = 4;
    cout << "Fill number at " << poscount << endl;
  }
  else if (strstr(token, "energy")) {
    position = 5;
    cout << "Beam Energy at " << poscount << endl;
  }
  else if (strstr(token, "intensity_per_bunch")) {
    position = 6;
    cout << "Bunch Intensity at " << poscount << endl;
  }
  else if (strstr(token, "cbeamb_abce_nopf_all")) {
    position = 7;
    cout << "Beam Triggers at " << poscount << endl;
  }
  else if (strstr(token, "cint1b_abce_nopf_all")) {
    position = 8;
    cout << "Collision Triggers at " << poscount << endl;
  }
  else if (strstr(token, "rate")) {
    position = 9;
    cout << "Collision rate at " << poscount << endl;
  }
  else if (strstr(token, "cint1_e_nopf_all")) {
    position = 10;
    cout << "Empty Triggers at " << poscount << endl;
  }
  else if (strstr(token, "cint1a_abce_nopf_all")) {
    position = 11;
    cout << "ASide Triggers at " << poscount << endl;
  }
  else if (strstr(token, "cint1c_abce_nopf_all")) {
    position = 12;
    cout << "CSide Triggers at " << poscount << endl;
  }
  else if (strstr(token, "mean_vertex_xyz")) {
    position = 13;
    cout << "Mean Vertex at " << poscount << endl;
  }
  else if (strstr(token, "vertex_quality")) {
    position = 14;
    cout << "Vertex Quality at " << poscount << endl;
  }
  else if (strstr(token, "field")) {
    position = 15;
    cout << "Magnetic Field at " << poscount << endl;
  }
  else if (strstr(token, "det_aco")) {
    position = 16;
    cout << "ACORDE at " << poscount << endl;
  }
  else if (strstr(token, "det_emc")) {
    position = 17;
    cout << "EMC at " << poscount << endl;
  }
  else if (strstr(token, "det_fmd")) {
    position = 18;
    cout << "FMD at " << poscount << endl;
  }
  else if (strstr(token, "det_hlt")) {
    position = 19;
    cout << "HLT at " << poscount << endl;
  }
  else if (strstr(token, "det_hmp")) {
    position = 20;
    cout << "HMPID at " << poscount << endl;
  }
  else if (strstr(token, "det_mch")) {
    position = 21;
    cout << "MCH at " << poscount << endl;
  }
  else if (strstr(token, "det_mtr")) {
    position = 22;
    cout << "MTR at " << poscount << endl;
  }
  else if (strstr(token, "det_phs")) {
    position = 23;
    cout << "PHOS at " << poscount << endl;
  }
  else if (strstr(token, "det_pmd")) {
    position = 24;
    cout << "pmd at " << poscount << endl;
  }
  else if (strstr(token, "det_spd")) {
    position = 25;
    cout << "SPD at " << poscount << endl;
  }
  else if (strstr(token, "det_sdd")) {
    position = 26;
    cout << "SDD at " << poscount << endl;
  }
  else if (strstr(token, "det_ssd")) {
    position = 27;
    cout << "SSD at " << poscount << endl;
  }
  else if (strstr(token, "det_tof")) {
    position = 28;
    cout << "TOF at " << poscount << endl;
  }
  else if (strstr(token, "det_tpc")) {
    position = 29;
    cout << "TPC at " << poscount << endl;
  }
  else if (strstr(token, "det_trd")) {
    position = 30;
    cout << "TRD at " << poscount << endl;
  }
  else if (strstr(token, "det_t00")) {
    position = 31;
    cout << "TZERO at " << poscount << endl;
  }
  else if (strstr(token, "det_v00")) {
    position = 32;
    cout << "VZERO at " << poscount << endl;
  }
  else if (strstr(token, "det_zdc")) {
    position = 33;
    cout << "ZDC at " << poscount << endl;
  }
  else if (strstr(token, "filling_config")) {
    position = 39;
    cout << "Filling Scheme at " << poscount << endl;
  }
  else if (strstr(token, "quality")) {
    position = 40;
    cout << "Quality at " << poscount << endl;
  }
  else if (strstr(token, "muon_trigger")) {
    position = 41;
    cout << "Muon Trigger at " << poscount << endl;
  }
  else if (strstr(token, "high_multiplicity_trigger")) {
    position = 42;
    cout << "High Multiplicity trigger at " << poscount << endl;
  }

  return position;
}

void ParseRCTDump(const char *insum = "rctdump.txt")
{
  ifstream *istr = new ifstream(insum);

  char buf[10000];
  istr->getline(buf, 10000);
  
  if (strstr(buf, "#raw_run"))
    cout << "Found raw run " << endl;

  char *nval;
  char *nlin = strstr(buf, "|");
  char *nbeg = buf-1;
  char ntok[1000];
  char nstr[1000];

  // Parse field name list

  int fieldpos[200];
  int poscount = 0;
  
  do {
    if (nlin == NULL) {
      strcpy(ntok, nbeg+1);
      //      ntok[nlin-nbeg-1] = '\0';
    }
    else {
      strncpy(ntok, nbeg+1, nlin-nbeg-1);
      ntok[nlin-nbeg-1] = '\0';
    }

    cout << "Token is " << ntok << endl;

    fieldpos[poscount] = getposition(ntok, poscount);
    
    nbeg = nlin;
    nlin = strstr(nlin+1, "|");
    poscount++;
  }
  while (nlin != NULL);

  strcpy(ntok, nbeg+1);

  cout << "Token is " << ntok << endl;
  
  fieldpos[poscount] = getposition(ntok, poscount);
  
  TFile* ftag = TFile::Open("RunTableTagSummary.root", "recreate");

  AliRunTag *tag = new AliRunTag();
  TTree * ttag = new TTree("T","A Tree with event tags");
  TBranch * btag = ttag->Branch("AliTAG", &tag);
  btag->SetCompressionLevel(9);

  int nfield = 0;


  while (!istr->eof()) {
    istr->getline(buf, 10000);
    nlin = strstr(buf, "|");
    nbeg = buf-1;

    nfield = 0;

    cout << "Line " << nlin << endl;
    
    int endline = 0;

    while (!endline) {
      
      if (nlin == NULL) {
	strcpy(ntok, nbeg+1);
	endline = 1;
      }
      else {
	strncpy(ntok, nbeg+1, nlin-nbeg-1);
	ntok[nlin-nbeg-1] = '\0';
      }
      if (ntok[0] == '"') {
	strncpy(nstr, ntok+1, strlen(ntok)-2); 
	nstr[strlen(ntok)-2] = '\0';
      }
      else if (ntok[0] == '(') {
	char *com = strstr(ntok, ";");
	cout << ntok << " "<< com << endl;
	strncpy(nstr, ntok+1, com-ntok-1);
	nstr[com-ntok-1] = '\0';
      }
      else {
	nstr[0] = '\0';
	nstr[1] = '\0';
	nstr[2] = '\0';
      }

      if (nlin != NULL) {
	nbeg = nlin;
	nlin = strstr(nlin+1, "|");
      }

      switch (fieldpos[nfield]) {
      case 0:
	tag->SetRunId(atoi(ntok));
	break;
      case 1:
	tag->SetLHCPeriod(nstr);
	break;
      case 2:
	break;
      case 3:
	break;
      case 4:
	tag->GetLHCTag()->SetFillNo(atoi(ntok));
	break;
      case 5:
	tag->GetLHCTag()->SetBeamEnergy(atof(ntok));
	break;
      case 6:
	tag->GetLHCTag()->SetBunchIntensity(atof(ntok));
	break;
      case 7:
	tag->SetBeamTriggers(atol(ntok));
	break;
      case 8:
	tag->SetCollisionTriggers(atol(ntok));
	break;
      case 9:
	tag->SetCollisionRate(atof(ntok));
	break;
      case 10:
	tag->SetEmptyTriggers(atol(ntok));
	break;
      case 11:
	tag->SetASideTriggers(atol(ntok));
	break;
      case 12:
	tag->SetCSideTriggers(atol(ntok));
	break;
      case 13:
	tag->SetMeanVertex(atof(ntok));
	break;
      case 14:
	tag->SetVertexQuality(atof(ntok));
	break;
      case 15:
	tag->SetMagneticField(atoi(ntok+1));
	break;
      case 16:
	tag->GetDetectorTags()->SetDetectorStatus(17, nstr);
	break;
      case 17:
	tag->GetDetectorTags()->SetDetectorStatus(18, nstr);
	break;
      case 18:
	tag->GetDetectorTags()->SetDetectorStatus(12, nstr);
	break;
      case 19:
	tag->GetDetectorTags()->SetDetectorStatus(21, nstr);
	break;
      case 20:
	tag->GetDetectorTags()->SetDetectorStatus(6, nstr);
	break;
      case 21:
	tag->GetDetectorTags()->SetDetectorStatus(10, nstr);
	break;
      case 22:
	tag->GetDetectorTags()->SetDetectorStatus(11, nstr);
	break;
      case 23:
	tag->GetDetectorTags()->SetDetectorStatus(7, nstr);
	break;
      case 24:
	tag->GetDetectorTags()->SetDetectorStatus(9, nstr);
	break;
      case 25:
	tag->GetDetectorTags()->SetDetectorStatus(0, nstr);
	break;
      case 26:
	tag->GetDetectorTags()->SetDetectorStatus(1, nstr);
	break;
      case 27:
	tag->GetDetectorTags()->SetDetectorStatus(2, nstr);
	break;
      case 28:
	tag->GetDetectorTags()->SetDetectorStatus(5, nstr);
	break;
      case 29:
	tag->GetDetectorTags()->SetDetectorStatus(3, nstr);
	break;
      case 30:
	tag->GetDetectorTags()->SetDetectorStatus(4, nstr);
	break;
      case 31:
	tag->GetDetectorTags()->SetDetectorStatus(13, nstr);
	break;
      case 32:
	tag->GetDetectorTags()->SetDetectorStatus(14, nstr);
	break;
      case 33:
	tag->GetDetectorTags()->SetDetectorStatus(15, nstr);
	break;
      case 34:
	// comment
	break;
      case 35:
	// detectors
	break;
      case 36:
	// changedby
	break;
      case 37:
	// changedon
	break;
      case 38:
	// run duration
	break;
      case 39:
	tag->GetLHCTag()->SetFillingScheme(nstr);
	break;
      case 40:
	cout << "Setting Quality ]" << nstr << "[" << atoi(nstr) << endl;
	if (strlen(nstr) == 0)
	  tag->SetRunQuality(0);
	else
	  tag->SetRunQuality(atoi(nstr));
	break;
      case 41:
	tag->SetMuonTriggers(atol(ntok));
	break;
      case 42:
	tag->SetHMTriggers(atol(ntok));
	break;
      }

      nfield++;
    }


    if (nfield > 1) {
      ttag->Fill();
      tag->Clear("");
    }
  }

  ftag->cd();
  tag->Clear();
  ttag->Write();
  ftag->Close();

}
Commit	Line	Data
b6461bb1	1	// Macro to parse the text dump of the RCT and convert it into
	2	// a RunTag object. This object needs to be then merged with the
	3	// proper RunTag from the Run.
	4	//
	5	// Parameters:
	6	// Name of the text file with the dump from MonaLisa RCT page
	7	//
	8	// Author: Adam.Kisiel@cern.ch
	9
	10	int getposition(const char *token, int poscount)
	11	{
	12	int position = -1;
	13
	14	if (strstr(token, "#raw_run")) {
	15	position = 0;
	16	cout << "Run number at " << poscount << endl;
	17	}
	18	else if (strstr(token, "partition")) {
	19	position = 1;
	20	cout << "LHC Period at " << poscount << endl;
	21	}
	22	else if (strstr(token, "fillno")) {
	23	position = 4;
	24	cout << "Fill number at " << poscount << endl;
	25	}
	26	else if (strstr(token, "energy")) {
	27	position = 5;
	28	cout << "Beam Energy at " << poscount << endl;
	29	}
	30	else if (strstr(token, "intensity_per_bunch")) {
	31	position = 6;
	32	cout << "Bunch Intensity at " << poscount << endl;
	33	}
	34	else if (strstr(token, "cbeamb_abce_nopf_all")) {
	35	position = 7;
	36	cout << "Beam Triggers at " << poscount << endl;
	37	}
	38	else if (strstr(token, "cint1b_abce_nopf_all")) {
	39	position = 8;
	40	cout << "Collision Triggers at " << poscount << endl;
	41	}
	42	else if (strstr(token, "rate")) {
	43	position = 9;
	44	cout << "Collision rate at " << poscount << endl;
	45	}
	46	else if (strstr(token, "cint1_e_nopf_all")) {
	47	position = 10;
	48	cout << "Empty Triggers at " << poscount << endl;
	49	}
	50	else if (strstr(token, "cint1a_abce_nopf_all")) {
	51	position = 11;
	52	cout << "ASide Triggers at " << poscount << endl;
	53	}
	54	else if (strstr(token, "cint1c_abce_nopf_all")) {
	55	position = 12;
	56	cout << "CSide Triggers at " << poscount << endl;
	57	}
	58	else if (strstr(token, "mean_vertex_xyz")) {
	59	position = 13;
	60	cout << "Mean Vertex at " << poscount << endl;
	61	}
	62	else if (strstr(token, "vertex_quality")) {
	63	position = 14;
	64	cout << "Vertex Quality at " << poscount << endl;
65	}
66	else if (strstr(token, "field")) {
67	position = 15;
68	cout << "Magnetic Field at " << poscount << endl;
69	}
70	else if (strstr(token, "det_aco")) {
71	position = 16;
72	cout << "ACORDE at " << poscount << endl;
73	}
74	else if (strstr(token, "det_emc")) {
75	position = 17;
76	cout << "EMC at " << poscount << endl;
77	}
78	else if (strstr(token, "det_fmd")) {
79	position = 18;
80	cout << "FMD at " << poscount << endl;
81	}
82	else if (strstr(token, "det_hlt")) {
83	position = 19;
84	cout << "HLT at " << poscount << endl;
85	}
86	else if (strstr(token, "det_hmp")) {
87	position = 20;
88	cout << "HMPID at " << poscount << endl;
89	}
90	else if (strstr(token, "det_mch")) {
91	position = 21;
92	cout << "MCH at " << poscount << endl;
93	}
94	else if (strstr(token, "det_mtr")) {
95	position = 22;
96	cout << "MTR at " << poscount << endl;
97	}
98	else if (strstr(token, "det_phs")) {
99	position = 23;
100	cout << "PHOS at " << poscount << endl;
101	}
102	else if (strstr(token, "det_pmd")) {
103	position = 24;
104	cout << "pmd at " << poscount << endl;
105	}
106	else if (strstr(token, "det_spd")) {
107	position = 25;
108	cout << "SPD at " << poscount << endl;
109	}
110	else if (strstr(token, "det_sdd")) {
111	position = 26;
112	cout << "SDD at " << poscount << endl;
113	}
114	else if (strstr(token, "det_ssd")) {
115	position = 27;
116	cout << "SSD at " << poscount << endl;
117	}
118	else if (strstr(token, "det_tof")) {
119	position = 28;
120	cout << "TOF at " << poscount << endl;
121	}
122	else if (strstr(token, "det_tpc")) {
123	position = 29;
124	cout << "TPC at " << poscount << endl;
125	}
126	else if (strstr(token, "det_trd")) {
127	position = 30;
128	cout << "TRD at " << poscount << endl;
129	}
130	else if (strstr(token, "det_t00")) {
131	position = 31;
132	cout << "TZERO at " << poscount << endl;
133	}
134	else if (strstr(token, "det_v00")) {
135	position = 32;
136	cout << "VZERO at " << poscount << endl;
137	}
138	else if (strstr(token, "det_zdc")) {
139	position = 33;
140	cout << "ZDC at " << poscount << endl;
141	}
142	else if (strstr(token, "filling_config")) {
143	position = 39;
144	cout << "Filling Scheme at " << poscount << endl;
145	}
146	else if (strstr(token, "quality")) {
147	position = 40;
148	cout << "Quality at " << poscount << endl;
149	}
150	else if (strstr(token, "muon_trigger")) {
151	position = 41;
152	cout << "Muon Trigger at " << poscount << endl;
153	}
154	else if (strstr(token, "high_multiplicity_trigger")) {
155	position = 42;
156	cout << "High Multiplicity trigger at " << poscount << endl;
157	}
158
159	return position;
160	}
161
162	void ParseRCTDump(const char *insum = "rctdump.txt")
163	{
164	ifstream *istr = new ifstream(insum);
165
166	char buf[10000];
167	istr->getline(buf, 10000);
168
169	if (strstr(buf, "#raw_run"))
170	cout << "Found raw run " << endl;
171
172	char *nval;
173	char *nlin = strstr(buf, "\|");
174	char *nbeg = buf-1;
175	char ntok[1000];
176	char nstr[1000];
177
178	// Parse field name list
179
180	int fieldpos[200];
181	int poscount = 0;
182
183	do {
184	if (nlin == NULL) {
185	strcpy(ntok, nbeg+1);
186	// ntok[nlin-nbeg-1] = '\0';
187	}
188	else {
189	strncpy(ntok, nbeg+1, nlin-nbeg-1);
190	ntok[nlin-nbeg-1] = '\0';
191	}
192
193	cout << "Token is " << ntok << endl;
194
195	fieldpos[poscount] = getposition(ntok, poscount);
196
197	nbeg = nlin;
198	nlin = strstr(nlin+1, "\|");
199	poscount++;
200	}
201	while (nlin != NULL);
202
203	strcpy(ntok, nbeg+1);
204
205	cout << "Token is " << ntok << endl;
206
207	fieldpos[poscount] = getposition(ntok, poscount);
208
209	TFile* ftag = TFile::Open("RunTableTagSummary.root", "recreate");
210
211	AliRunTag *tag = new AliRunTag();
212	TTree * ttag = new TTree("T","A Tree with event tags");
213	TBranch * btag = ttag->Branch("AliTAG", &tag);
214	btag->SetCompressionLevel(9);
215
216	int nfield = 0;
217
218
219	while (!istr->eof()) {
220	istr->getline(buf, 10000);
221	nlin = strstr(buf, "\|");
222	nbeg = buf-1;
223
224	nfield = 0;
225
226	cout << "Line " << nlin << endl;
227
228	int endline = 0;
229
230	while (!endline) {
231
232	if (nlin == NULL) {
233	strcpy(ntok, nbeg+1);
234	endline = 1;
235	}
236	else {
237	strncpy(ntok, nbeg+1, nlin-nbeg-1);
238	ntok[nlin-nbeg-1] = '\0';
239	}
240	if (ntok[0] == '"') {
241	strncpy(nstr, ntok+1, strlen(ntok)-2);
242	nstr[strlen(ntok)-2] = '\0';
243	}
244	else if (ntok[0] == '(') {
245	char *com = strstr(ntok, ";");
246	cout << ntok << " "<< com << endl;
247	strncpy(nstr, ntok+1, com-ntok-1);
248	nstr[com-ntok-1] = '\0';
249	}
250	else {
251	nstr[0] = '\0';
252	nstr[1] = '\0';
253	nstr[2] = '\0';
254	}
255
256	if (nlin != NULL) {
257	nbeg = nlin;
258	nlin = strstr(nlin+1, "\|");
259	}
260
261	switch (fieldpos[nfield]) {
262	case 0:
263	tag->SetRunId(atoi(ntok));
264	break;
265	case 1:
266	tag->SetLHCPeriod(nstr);
267	break;
268	case 2:
269	break;
270	case 3:
271	break;
272	case 4:
273	tag->GetLHCTag()->SetFillNo(atoi(ntok));
274	break;
275	case 5:
276	tag->GetLHCTag()->SetBeamEnergy(atof(ntok));
277	break;
278	case 6:
279	tag->GetLHCTag()->SetBunchIntensity(atof(ntok));
280	break;
281	case 7:
282	tag->SetBeamTriggers(atol(ntok));
283	break;
284	case 8:
285	tag->SetCollisionTriggers(atol(ntok));
286	break;
287	case 9:
288	tag->SetCollisionRate(atof(ntok));
289	break;
290	case 10:
291	tag->SetEmptyTriggers(atol(ntok));
292	break;
293	case 11:
294	tag->SetASideTriggers(atol(ntok));
295	break;
296	case 12:
297	tag->SetCSideTriggers(atol(ntok));
298	break;
299	case 13:
300	tag->SetMeanVertex(atof(ntok));
301	break;
302	case 14:
303	tag->SetVertexQuality(atof(ntok));
304	break;
305	case 15:
306	tag->SetMagneticField(atoi(ntok+1));
307	break;
308	case 16:
309	tag->GetDetectorTags()->SetDetectorStatus(17, nstr);
310	break;
311	case 17:
312	tag->GetDetectorTags()->SetDetectorStatus(18, nstr);
313	break;
314	case 18:
315	tag->GetDetectorTags()->SetDetectorStatus(12, nstr);
316	break;
317	case 19:
318	tag->GetDetectorTags()->SetDetectorStatus(21, nstr);
319	break;
320	case 20:
321	tag->GetDetectorTags()->SetDetectorStatus(6, nstr);
322	break;
323	case 21:
324	tag->GetDetectorTags()->SetDetectorStatus(10, nstr);
325	break;
326	case 22:
327	tag->GetDetectorTags()->SetDetectorStatus(11, nstr);
328	break;
329	case 23:
330	tag->GetDetectorTags()->SetDetectorStatus(7, nstr);
331	break;
332	case 24:
333	tag->GetDetectorTags()->SetDetectorStatus(9, nstr);
334	break;
335	case 25:
336	tag->GetDetectorTags()->SetDetectorStatus(0, nstr);
337	break;
338	case 26:
339	tag->GetDetectorTags()->SetDetectorStatus(1, nstr);
340	break;
341	case 27:
342	tag->GetDetectorTags()->SetDetectorStatus(2, nstr);
343	break;
344	case 28:
345	tag->GetDetectorTags()->SetDetectorStatus(5, nstr);
346	break;
347	case 29:
348	tag->GetDetectorTags()->SetDetectorStatus(3, nstr);
349	break;
350	case 30:
351	tag->GetDetectorTags()->SetDetectorStatus(4, nstr);
352	break;
353	case 31:
354	tag->GetDetectorTags()->SetDetectorStatus(13, nstr);
355	break;
356	case 32:
357	tag->GetDetectorTags()->SetDetectorStatus(14, nstr);
358	break;
359	case 33:
360	tag->GetDetectorTags()->SetDetectorStatus(15, nstr);
361	break;
362	case 34:
363	// comment
364	break;
365	case 35:
366	// detectors
367	break;
368	case 36:
369	// changedby
370	break;
371	case 37:
372	// changedon
373	break;
374	case 38:
375	// run duration
376	break;
377	case 39:
378	tag->GetLHCTag()->SetFillingScheme(nstr);
379	break;
380	case 40:
381	cout << "Setting Quality ]" << nstr << "[" << atoi(nstr) << endl;
382	if (strlen(nstr) == 0)
383	tag->SetRunQuality(0);
384	else
385	tag->SetRunQuality(atoi(nstr));
386	break;
387	case 41:
388	tag->SetMuonTriggers(atol(ntok));
389	break;
390	case 42:
391	tag->SetHMTriggers(atol(ntok));
392	break;
393	}
394
395	nfield++;
396	}
397
398
399	if (nfield > 1) {
400	ttag->Fill();
401	tag->Clear("");
402	}
403	}
404
405	ftag->cd();
406	tag->Clear();
407	ttag->Write();
408	ftag->Close();
409
410	}