1 /**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
7 * Permission to use, copy, modify and distribute this software and its *
8 * documentation strictly for non-commercial purposes is hereby granted *
9 * without fee, provided that the above copyright notice appears in all *
10 * copies and that both the copyright notice and this permission notice *
11 * appear in the supporting documentation. The authors make no claims *
12 * about the suitability of this software for any purpose. It is *
13 * provided "as is" without express or implied warranty. *
14 **************************************************************************/
18 ////////////////////////////////////////////////////////////////////////////
20 // This class is a first implementation for the TRD. //
21 // It takes data from HLT and computes the parameters //
22 // and stores both reference data and online calibration //
23 // parameters in the CDB //
24 // It alsotakes DCS data, does spline fits //
25 // and stores both reference data and spline fits results //
29 // R. Bailhache (R.Bailhache@gsi.de) //
30 // W. Monange (w.monange@gsi.de) //
32 ////////////////////////////////////////////////////////////////////////////
35 #include <TProfile2D.h>
36 #include <TStopwatch.h>
37 #include <TObjString.h>
40 #include <TCollection.h>
41 //#include <TSAXParser.h>
43 #include "AliCDBMetaData.h"
46 #include "AliTRDPreprocessor.h"
47 #include "AliTRDSensorArray.h"
48 #include "AliTRDCalibraFit.h"
49 #include "AliTRDCalibraMode.h"
50 #include "AliTRDCalibPadStatus.h"
51 //#include "AliTRDSaxHandler.h"
53 #include "Cal/AliTRDCalDet.h"
54 #include "Cal/AliTRDCalPadStatus.h"
55 #include "Cal/AliTRDCalDCS.h"
57 ClassImp(AliTRDPreprocessor)
59 //______________________________________________________________________________________________
60 AliTRDPreprocessor::AliTRDPreprocessor(AliShuttleInterface *shuttle)
61 :AliPreprocessor("TRD", shuttle)
68 AddRunType("PHYSICS");
69 AddRunType("STANDALONE");
70 AddRunType("PEDESTAL");
74 //______________________________________________________________________________________________
75 AliTRDPreprocessor::~AliTRDPreprocessor()
83 //______________________________________________________________________________________________
84 void AliTRDPreprocessor::Initialize(Int_t run, UInt_t startTime, UInt_t endTime)
87 // Initialization routine for the TRD preprocessor
90 AliPreprocessor::Initialize(run,startTime,endTime);
94 //______________________________________________________________________________________________
95 UInt_t AliTRDPreprocessor::Process(TMap* dcsAliasMap)
98 // Process DCS and calibration part for HLT
101 TString runType = GetRunType();
102 Log(Form("runtype %s\n",runType.Data()));
104 // always process the configuration data
105 Int_t resultDCSC = ProcessDCSConfigData();
106 // if there was an error, return with its code
107 if (resultDCSC != 0) return resultDCSC;
109 if (runType=="PEDESTAL"){
110 if(ExtractPedestals()) return 1;
114 if ((runType=="PHYSICS") || (runType=="STANDALONE")){
116 if(ProcessDCS(dcsAliasMap)) return 1;
118 //TString runPar = GetRunParameter("HLTStatus");
121 if(ExtractHLT()) return 1;
125 if(ExtractDriftVelocityDAQ()) return 1;
133 //______________________________________________________________________________
134 Bool_t AliTRDPreprocessor::ProcessDCS()
137 // Default process DCS method
140 TString runType = GetRunType();
141 if ((runType == "PHYSICS") || (runType == "STANDALONE")) {
148 //______________________________________________________________________________
149 Bool_t AliTRDPreprocessor::ProcessDCS(TMap *dcsAliasMap)
152 // Process DCS method
155 Bool_t error = kFALSE;
157 AliCDBMetaData metaData;
158 metaData.SetBeamPeriod(0);
159 metaData.SetResponsible("Wilfried Monange/Raphaelle Bailhache");
160 metaData.SetComment("TRD calib test");
162 Log("****** DCS ******\n");
164 TObjArray * list=AliTRDSensorArray::GetList ();
167 Log ("Error during AliTRDSensorArray::GetList");
168 Log ("DCS will not be processing");
172 Int_t nEntries = list->GetEntries ();
173 Log (Form ("%d alias loaded", nEntries));
175 Bool_t * results=new Bool_t [nEntries];
176 Int_t * nGraph=new Int_t [nEntries];
178 for (Int_t iAlias = 0; iAlias < nEntries; iAlias++) {
180 AliTRDSensorArray * oneTRDDCS = (AliTRDSensorArray *)list->At (iAlias);
182 oneTRDDCS->SetStartTime (TTimeStamp (fStartTime));
183 oneTRDDCS->SetEndTime (TTimeStamp (fEndTime));
185 Log(Form("Processing DCS : \"%s\"", oneTRDDCS->GetStoreName ().Data ()));
189 map=oneTRDDCS->ExtractDCS (dcsAliasMap);
191 nGraph [iAlias] = map->GetEntries ();
193 if (nGraph [iAlias] == 0) {
194 Log("No TGraph for this dcsDatapointAlias : not stored");
195 results [iAlias] = kFALSE;
200 oneTRDDCS->SetGraph(map);
201 results[iAlias]=Store("Calib", oneTRDDCS->GetStoreName().Data(), oneTRDDCS, &metaData, 0, kTRUE);
204 //results [iAlias] = StoreReferenceData("Calib", oneTRDDCS->GetStoreName ().Data (), oneTRDDCS, &metaData);
206 if (!results[iAlias]) {
207 AliError("Problem during StoreRef DCS");
211 //BEGIN TEST (should not be removed ...)
213 oneTRDDCS->ClearGraph();
214 oneTRDDCS->ClearFit();
215 oneTRDDCS->SetDiffCut2 (0.1);
216 map=oneTRDDCS->ExtractDCS (dcsAliasMap);
217 oneTRDDCS->SetGraph (map);
218 Store("Calib", ("cut_"+oneTRDDCS->GetStoreName()).Data(), oneTRDDCS, &metaData, 0, kTRUE);
222 if(iAlias==1 || iAlias==19) continue;
224 oneTRDDCS->ClearGraph();
225 oneTRDDCS->ClearFit();
226 oneTRDDCS->SetDiffCut2(0);
227 map=oneTRDDCS->ExtractDCS(dcsAliasMap);
228 oneTRDDCS->MakeSplineFit(map);
229 Store("Calib", ("fit_"+oneTRDDCS->GetStoreName()).Data() , oneTRDDCS, &metaData, 0, kTRUE);
233 oneTRDDCS->ClearGraph();
234 oneTRDDCS->ClearFit();
235 oneTRDDCS->SetDiffCut2 (0.1);
236 map=oneTRDDCS->ExtractDCS (dcsAliasMap);
237 oneTRDDCS->MakeSplineFit(map);
238 Store("Calib", ("cutfit_"+oneTRDDCS->GetStoreName()).Data() , oneTRDDCS, &metaData, 0, kTRUE);
245 Log (" Summury of DCS :\n");
246 Log (Form("%30s %10s %10s", "dcsDatapointAlias", "Stored ?", "# graph"));
247 for (Int_t iAlias = 0; iAlias < nEntries; iAlias++) {
248 AliTRDSensorArray * oneTRDDCS = (AliTRDSensorArray *)list->At (iAlias);
249 Log (Form ("%30s %10s %4d",
250 oneTRDDCS->GetStoreName ().Data (),
251 results[iAlias] ? "ok" : "X",
254 Log ("*********** End of DCS **********");
263 //______________________________________________________________________________________________
264 Bool_t AliTRDPreprocessor::ExtractPedestals()
267 // Pedestal running on LDCs at the DAQ
270 Bool_t error = kFALSE;
272 // Init a AliTRDCalibPadStatus
273 AliTRDCalibPadStatus calPedSum = AliTRDCalibPadStatus();
275 AliCDBMetaData metaData;
276 metaData.SetBeamPeriod(0);
277 metaData.SetResponsible("Raphaelle Bailhache");
278 metaData.SetComment("TRD calib test");
280 // Sum the contributions of the LDCs
281 TList * listpad = GetFileSources(kDAQ,"PADSTATUS");
283 Log("No list found for the PEDESTRAL Run");
287 // loop through all files from LDCs
289 while (listpad->At(index)!=NULL) {
290 TObjString* fileNameEntry = (TObjString*) listpad->At(index);
291 if (fileNameEntry != NULL)
293 TString fileName = GetFile(kDAQ, "PADSTATUS",
294 fileNameEntry->GetString().Data());
295 if(fileName.Length() ==0){
296 Log(Form("Error by retrieving the file %d for the pedestal",(Int_t)index));
301 TFile *f = TFile::Open(fileName);
302 AliTRDCalibPadStatus *calPed;
303 f->GetObject("calibpadstatus",calPed);
310 //calPed->AnalyseHisto();
312 // Add to the calPedSum
313 for (Int_t idet=0; idet<540; idet++) {
314 AliTRDCalROC *rocMean = calPed->GetCalRocMean(idet, kFALSE);
316 calPedSum.SetCalRocMean(rocMean,idet);
317 sm = (Int_t) (idet / 30);
319 AliTRDCalROC *rocRMS = calPed->GetCalRocRMS(idet, kFALSE);
321 calPedSum.SetCalRocRMS(rocRMS,idet);
323 AliTRDCalROC *rocMeand = calPed->GetCalRocMeand(idet, kFALSE);
325 calPedSum.SetCalRocMeand(rocMeand,idet);
327 AliTRDCalROC *rocRMSd = calPed->GetCalRocRMSd(idet, kFALSE);
329 calPedSum.SetCalRocRMSd(rocRMSd,idet);
333 // store as reference data
334 TString name("PadStatus");
336 if(!StoreReferenceData("DAQData",(const char *)name,(TObject *) calPed,&metaData)){
337 Log(Form("Error storing AliTRDCalibPadStatus object %d as reference data",(Int_t)index));
346 Log(Form("%d elements found in the list for the pedestal",(Int_t)index));
353 // Store pedestal entry to OCDB
357 AliTRDCalPadStatus *calPadStatus = calPedSum.CreateCalPadStatus();
359 md3.SetObjectClassName("AliTRDCalPadStatus");
360 md3.SetResponsible("Raphaelle Bailhache");
361 md3.SetBeamPeriod(1);
362 md3.SetComment("TRD calib test");
363 if(!Store("Calib","PadStatus" ,(TObject *)calPadStatus, &md3, 0, kTRUE)){
364 Log("Error storing the pedestal");
370 //Make the AliTRDCalPad
371 AliTRDCalPad *calPad2 = calPedSum.CreateCalPad();
373 md4.SetObjectClassName("AliTRDCalPad");
374 md4.SetResponsible("Raphaelle Bailhache");
375 md4.SetBeamPeriod(1);
376 md4.SetComment("TRD calib test");
377 if(!Store("Calib","PadNoise" ,(TObject *)calPad2, &md4, 0, kTRUE)){
378 Log("Error storing the pedestal");
382 //Make the AliTRDCalDet correspondant
383 AliTRDCalDet *calDet = calPedSum.CreateCalDet();
385 md5.SetObjectClassName("AliTRDCalDet");
386 md5.SetResponsible("Raphaelle Bailhache");
387 md5.SetBeamPeriod(1);
388 md5.SetComment("TRD calib test");
389 if(!Store("Calib","DetNoise" ,(TObject *)calDet, &md5, 0, kTRUE)){
390 Log("Error storing the pedestal");
400 //______________________________________________________________________________________________
401 Bool_t AliTRDPreprocessor::ExtractDriftVelocityDAQ()
404 // Drift velocity DA running on monitoring servers at the DAQ
407 Bool_t error = kFALSE;
409 // Objects for HLT and DAQ zusammen
410 AliTRDCalibraFit *calibra = AliTRDCalibraFit::Instance();
411 AliCDBMetaData metaData;
412 metaData.SetBeamPeriod(0);
413 metaData.SetResponsible("Raphaelle Bailhache");
414 metaData.SetComment("TRD calib test");
415 // Store the infos for the detector
417 md1.SetObjectClassName("AliTRDCalDet");
418 md1.SetResponsible("Raphaelle Bailhache");
419 md1.SetBeamPeriod(0);
420 md1.SetComment("TRD calib test");
421 // Store the infos for the pads
423 md2.SetObjectClassName("AliTRDCalPad");
424 md2.SetResponsible("Raphaelle Bailhache");
425 md2.SetBeamPeriod(0);
426 md2.SetComment("TRD calib test");
428 // Take the file from the DAQ file exchange server
429 TList *listdaq = GetFileSources(kDAQ,"VDRIFT");
431 Log("No list found for vdrift (DAQ)");
435 if(listdaq->GetSize() !=1){
436 Log(Form("Problem on the size of the list: %d (DAQ)",listdaq->GetSize()));
441 TObjString* fileNameEntry = (TObjString*) listdaq->At(0);
442 if(fileNameEntry != NULL){
443 TString fileName = GetFile(kDAQ, "VDRIFT",
444 fileNameEntry->GetString().Data());
445 if(fileName.Length() ==0){
446 Log("Error retrieving the file vdrift (DAQ)");
450 TFile *filedaq = TFile::Open(fileName);
451 TProfile2D *histodriftvelocity = (TProfile2D *) filedaq->Get("PH2d");
452 if (histodriftvelocity) {
454 // store as reference data
455 if(!StoreReferenceData("DAQData","VdriftT0",(TObject *) histodriftvelocity,&metaData)){
456 Log("Error storing 2D Profile for vdrift from the DAQ");
462 Log("Take the PH reference data. Now we will try to fit\n");
463 calibra->SetMinEntries(2000); // If there is less than 2000
464 calibra->AnalysePH(histodriftvelocity);
466 Int_t nbtg = 6*4*18*((Int_t) ((AliTRDCalibraMode *)calibra->GetCalibraMode())->GetDetChamb0(1))
467 + 6* 18*((Int_t) ((AliTRDCalibraMode *)calibra->GetCalibraMode())->GetDetChamb2(1));
468 Int_t nbfit = calibra->GetNumberFit();
469 Int_t nbE = calibra->GetNumberEnt();
471 // if enough statistics store the results
473 (nbfit >= 0.95*nbE)) {
474 // create the cal objects
475 TObjArray object = calibra->GetVectorFit();
476 AliTRDCalDet *objdriftvelocitydet = calibra->CreateDetObjectVdrift(&object,kTRUE);
477 TObject *objdriftvelocitypad = calibra->CreatePadObjectVdrift();
478 object = calibra->GetVectorFit2();
479 AliTRDCalDet *objtime0det = calibra->CreateDetObjectT0(&object,kTRUE);
480 TObject *objtime0pad = calibra->CreatePadObjectT0();
481 calibra->ResetVectorFit();
483 if(!Store("Calib","ChamberVdrift" ,(TObject *) objdriftvelocitydet,&md1,0,kTRUE)){
484 Log("Error storing the calibration object for the chamber vdrift (DAQ)");
487 if(!Store("Calib","ChamberT0" ,(TObject *) objtime0det ,&md1,0,kTRUE)){
488 Log("Error storing the calibration object for the chamber t0 (DAQ)");
491 if(!Store("Calib","LocalVdrift" ,(TObject *) objdriftvelocitypad,&md2,0,kTRUE)){
492 Log("Error storing the calibration object for the local drift velocity (DAQ)");
495 if(!Store("Calib","LocalT0" ,(TObject *) objtime0pad ,&md2,0,kTRUE)){
496 Log("Error storing the calibration object for the local time0 (DAQ)");
501 Log("Not enough statistics for the average pulse height (DAQ)");
511 //______________________________________________________________________________________________
512 Bool_t AliTRDPreprocessor::ExtractHLT()
515 // Gain, vdrift and PRF calibration running on HLT
516 // return kTRUE if NULL pointer to the list
519 Bool_t error = kFALSE;
521 // Objects for HLT and DAQ zusammen
522 AliTRDCalibraFit *calibra = AliTRDCalibraFit::Instance();
523 AliCDBMetaData metaData;
524 metaData.SetBeamPeriod(0);
525 metaData.SetResponsible("Raphaelle Bailhache");
526 metaData.SetComment("TRD calib test");
527 // Store the infos for the detector
529 md1.SetObjectClassName("AliTRDCalDet");
530 md1.SetResponsible("Raphaelle Bailhache");
531 md1.SetBeamPeriod(0);
532 md1.SetComment("TRD calib test");
533 // Store the infos for the pads
535 md2.SetObjectClassName("AliTRDCalPad");
536 md2.SetResponsible("Raphaelle Bailhache");
537 md2.SetBeamPeriod(0);
538 md2.SetComment("TRD calib test");
540 // Take the file from the HLT file exchange server
541 TList *listhlt = GetFileSources(kHLT,"GAINDRIFTPRF");
543 Log("No list found for the HLT");
547 if(listhlt->GetSize() != 1) {
548 Log(Form("Problem on the size of the list: %d (HLT)",listhlt->GetSize()));
553 TObjString* fileNameEntry = (TObjString*) listhlt->At(0);
554 if(fileNameEntry != NULL){
555 TString fileName = GetFile(kHLT, "GAINDRIFTPRF",
556 fileNameEntry->GetString().Data());
557 if(fileName.Length() ==0){
558 Log("Error retrieving the file (HLT)");
563 TFile *filehlt = TFile::Open(fileName);
566 TH2I *histogain = (TH2I *) filehlt->Get("CH2d");
567 histogain->SetDirectory(0);
569 // store the reference data
570 if(!StoreReferenceData("HLTData","Gain",(TObject *) histogain,&metaData)){
571 Log("Error storing 2D histos for gain");
575 Log("Take the CH reference data. Now we will try to fit\n");
576 calibra->SetMinEntries(1000); // If there is less than 1000 entries in the histo: no fit
577 calibra->AnalyseCH(histogain);
578 Int_t nbtg = 6*4*18*((Int_t) ((AliTRDCalibraMode *)calibra->GetCalibraMode())->GetDetChamb0(0))
579 + 6* 18*((Int_t) ((AliTRDCalibraMode *)calibra->GetCalibraMode())->GetDetChamb2(0));
580 Int_t nbfit = calibra->GetNumberFit();
581 Int_t nbE = calibra->GetNumberEnt();
584 (nbfit >= 0.95*nbE)) {
585 // create the cal objects
586 TObjArray object = calibra->GetVectorFit();
587 AliTRDCalDet *objgaindet = calibra->CreateDetObjectGain(&object);
588 TObject *objgainpad = calibra->CreatePadObjectGain();
590 if(!Store("Calib","ChamberGainFactor",(TObject *) objgaindet ,&md1,0,kTRUE)){
591 Log("Error storing the calibration object for the chamber gain");
594 if(!Store("Calib","LocalGainFactor" ,(TObject *) objgainpad ,&md2,0,kTRUE)){
595 Log("Error storing the calibration object for the local gain factor");
599 calibra->ResetVectorFit();
604 TProfile2D *histodriftvelocity = (TProfile2D *) filehlt->Get("PH2d");
605 histodriftvelocity->SetDirectory(0);
606 if (histodriftvelocity) {
607 // store the reference data
608 if(!StoreReferenceData("HLTData","VdriftT0",(TObject *) histodriftvelocity,&metaData)){
609 Log("Error storing 2D Profile for average pulse height (HLT)");
613 Log("Take the PH reference data. Now we will try to fit\n");
614 calibra->SetMinEntries(1000*20); // If there is less than 20000
615 calibra->AnalysePH(histodriftvelocity);
616 Int_t nbtg = 6*4*18*((Int_t) ((AliTRDCalibraMode *)calibra->GetCalibraMode())->GetDetChamb0(1))
617 + 6* 18*((Int_t) ((AliTRDCalibraMode *)calibra->GetCalibraMode())->GetDetChamb2(1));
618 Int_t nbfit = calibra->GetNumberFit();
619 Int_t nbE = calibra->GetNumberEnt();
622 (nbfit >= 0.95*nbE)) {
623 // create the cal objects
624 TObjArray object = calibra->GetVectorFit();
625 AliTRDCalDet *objdriftvelocitydet = calibra->CreateDetObjectVdrift(&object,kTRUE);
626 TObject *objdriftvelocitypad = calibra->CreatePadObjectVdrift();
627 object = calibra->GetVectorFit2();
628 AliTRDCalDet *objtime0det = calibra->CreateDetObjectT0(&object,kTRUE);
629 TObject *objtime0pad = calibra->CreatePadObjectT0();
631 if(!Store("Calib","ChamberVdrift" ,(TObject *) objdriftvelocitydet,&md1,0,kTRUE)){
632 Log("Error storing the calibration object for the chamber vdrift (HLT)");
635 if(!Store("Calib","ChamberT0" ,(TObject *) objtime0det ,&md1,0,kTRUE)){
636 Log("Error storing the calibration object for the chamber t0 (HLT)");
639 if(!Store("Calib","LocalVdrift" ,(TObject *) objdriftvelocitypad,&md2,0,kTRUE)){
640 Log("Error storing the calibration object for the local drift velocity (HLT)");
643 if(!Store("Calib","LocalT0" ,(TObject *) objtime0pad ,&md2,0,kTRUE)){
644 Log("Error storing the calibration object for the local time0 (HLT)");
649 calibra->ResetVectorFit();
653 TProfile2D *histoprf = (TProfile2D *) filehlt->Get("PRF2d");
654 histoprf->SetDirectory(0);
656 // store reference data
657 if(!StoreReferenceData("HLTData","PRF",(TObject *) histoprf,&metaData)){
658 Log("Error storing the 2D Profile for Pad Response Function");
662 Log("Take the PRF reference data. Now we will try to fit\n");
663 calibra->SetMinEntries(600); // If there is less than 20000
664 calibra->AnalysePRFMarianFit(histoprf);
665 Int_t nbtg = 6*4*18*((Int_t) ((AliTRDCalibraMode *)calibra->GetCalibraMode())->GetDetChamb0(2))
666 + 6* 18*((Int_t) ((AliTRDCalibraMode *)calibra->GetCalibraMode())->GetDetChamb2(2));
667 Int_t nbfit = calibra->GetNumberFit();
668 Int_t nbE = calibra->GetNumberEnt();
671 (nbfit >= 0.95*nbE)) {
672 // create cal pad objects
673 TObjArray object = calibra->GetVectorFit();
674 TObject *objPRFpad = calibra->CreatePadObjectPRF(&object);
676 if(!Store("Calib","PRFWidth" ,(TObject *) objPRFpad ,&md2,0,kTRUE)){
677 Log("Error storing the calibration object for the Pad Response Function");
681 calibra->ResetVectorFit();
690 //_____________________________________________________________________________
691 UInt_t AliTRDPreprocessor::ProcessDCSConfigData()
694 // process the configuration of FEE, PTR and GTU
695 // reteive XML file from the DCS FXS
696 // parse it and store TObjArrays in the CDB
697 // return 0 for success, otherwise:
698 // 5: could not get the file from the FXS
699 // 6: ERROR in XML validation: something wrong with the file
700 // 7: ERROR while creating calibration objects in the handler
701 // 8: error while storing data in the CDB
704 Log("Processing the DCS config summary file.");
707 const char * nameFile = GetFile(kDCS,"CONFIGSUMMARY","");
708 if (nameFile == NULL) {
710 Log(Form("File %s not found!",nameFile));
713 // // create parser and parse
714 // TSAXParser saxParser;
715 // AliTRDSaxHandler saxHandler;
716 // saxParser.ConnectToHandler("AliTRDSaxHandler", &saxHandler);
717 // saxParser.ParseFile(nameFile);
719 // // report errors if present
720 // if (saxParser.GetParseCode() == 0) {
721 // Log("XML file validation OK");
723 // Log(Form("ERROR in XML file validation. Parsecode: %s", saxParser.GetParseCode()));
726 // if (saxHandler.GetHandlerStatus() != 0) {
727 // Log(Form("ERROR while creating calibration objects. Error code: %s", saxHandler.GetHandlerStatus()));
731 // // get the calibration object storing the data from the handler
732 // AliTRDCalDCS* fCalDCSObj = saxHandler.GetCalDCSObj();
734 // // store the DCS calib data in the CDB
735 // AliCDBMetaData metaData1;
736 // metaData1.SetBeamPeriod(0);
737 // metaData1.SetResponsible("Frederick Kramer");
738 // metaData1.SetComment("DCS configuration data in one AliTRDCalDCS object.");
739 // if (!Store("Calib", "DCSCONFIG", fCalDCSObj, &metaData1, 0, kTRUE)) {
740 // Log("problems while storing DCS config data object");