Updated DAs
authorcoppedis <coppedis@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 3 Jun 2009 12:14:01 +0000 (12:14 +0000)
committercoppedis <coppedis@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 3 Jun 2009 12:14:01 +0000 (12:14 +0000)
ZDC/AliZDCRawStream.cxx
ZDC/ZDCLASERda.cxx
ZDC/ZDCPEDESTALda.cxx

index 2d3c8f9..58a54b0 100644 (file)
@@ -179,13 +179,13 @@ void AliZDCRawStream::ReadCDHHeader()
     //printf("\t AliZDCRawStream::ReadCDHHeader -> Attributes %x\n",message);
     
     if(message & 0x10){ // COSMIC RUN
-       printf("\t STANDALONE_COSMIC RUN raw data found\n");
+       //printf("\t STANDALONE_COSMIC RUN raw data found\n");
     }
     else if(message & 0x20){ // PEDESTAL RUN
-       printf("\t STANDALONE_PEDESTAL RUN raw data found\n");
+       //printf("\t STANDALONE_PEDESTAL RUN raw data found\n");
     }
     else if(message & 0x30){ // LASER RUN
-       printf("\t STANDALONE_LASER RUN raw data found\n");
+       //printf("\t STANDALONE_LASER RUN raw data found\n");
     }
     
     if(header->GetL1TriggerMessage() & 0x1){ // Calibration bit set in CDH
index 7e9961f..f986f2c 100644 (file)
@@ -23,6 +23,7 @@ Trigger Types Used: Standalone Trigger
 */
 #define PEDDATA_FILE  "ZDCPedestal.dat"
 #define MAPDATA_FILE  "ZDCChMapping.dat"
+#define LASHISTO_FILE "ZDCLaserHisto.root"
 #define LASDATA_FILE  "ZDCLaserCalib.dat"
 
 #include <stdio.h>
@@ -35,11 +36,13 @@ Trigger Types Used: Standalone Trigger
 #include <daqDA.h>
 
 //ROOT
-#include <TRandom.h>
+#include <TROOT.h>
+#include <TPluginManager.h>
 #include <TH1F.h>
 #include <TF1.h>
 #include <TFile.h>
 #include <TFitter.h>
+#include "TMinuitMinimizer.h"
 
 //AliRoot
 #include <AliRawReaderDate.h>
@@ -52,7 +55,16 @@ Trigger Types Used: Standalone Trigger
 */
 int main(int argc, char **argv) {
   
-//  TVirtualFitter::SetDefaultFitter("Minuit2");
+  gROOT->GetPluginManager()->AddHandler("TVirtualStreamerInfo",
+                                       "*",
+                                       "TStreamerInfo",
+                                       "RIO",
+                                       "TStreamerInfo()"); 
+
+  TMinuitMinimizer m; 
+  gROOT->GetPluginManager()->AddHandler("ROOT::Math::Minimizer", "Minuit","TMinuitMinimizer",
+      "Minuit", "TMinuitMinimizer(const char *)");
+  TVirtualFitter::SetDefaultFitter("Minuit");
 
   int status = 0;
   int const kNChannels = 24;
@@ -398,6 +410,12 @@ int main(int argc, char **argv) {
     if(nBin[k]!=0) maxXval[k] = maxBin[k]*xMax[k]/nBin[k];
     if(maxXval[k]-150.<0.) xlow[k]=0.;
     else xlow[k] = maxXval[k]-150.;
+    // checking if histos are empty
+    if(hZNChg[k]->GetEntries() == 0){
+      printf("\n WARNING! Empty LASER histos -> ending DA\n\n");
+      return -1;
+    } 
+    //
     hZNChg[k]->Fit("gaus","Q","",xlow[k],maxXval[k]+150.);
     fun[k] = hZNChg[k]->GetFunction("gaus");
     mean[k]  = (Float_t) (fun[k]->GetParameter(1));
@@ -589,6 +607,31 @@ int main(int argc, char **argv) {
   }
   //                                                  
   fclose(fileShuttle);
+  /* report progress */
+  daqDA_progressReport(80);
+  //
+  TFile *histofile = new TFile(LASHISTO_FILE,"RECREATE");
+  histofile->cd();
+  for(int j=0; j<5; j++){
+     hZNChg[j]->Write();
+     hZPChg[j]->Write();
+     hZNAhg[j]->Write();
+     hZPAhg[j]->Write();
+     hZNClg[j]->Write();
+     hZPClg[j]->Write();
+     hZNAlg[j]->Write();
+     hZPAlg[j]->Write();  
+     if(j<2){
+       hZEMhg[j]->Write();
+       hZEMlg[j]->Write();
+    }
+  }
+  hPMRefChg->Write();
+  hPMRefAhg->Write();
+  hPMRefClg->Write();
+  hPMRefAlg->Write();  
+  //
+  histofile->Close();
   //
   for(Int_t j=0; j<5; j++){
     delete hZNChg[j];
@@ -619,17 +662,25 @@ int main(int argc, char **argv) {
   daqDA_progressReport(90);
   
   /* store the result file on FES */
+  // [1] File with mapping
   status = daqDA_FES_storeFile(MAPDATA_FILE,MAPDATA_FILE);
   if(status){
     printf("Failed to export file : %d\n",status);
     return -1;
   }
   //
+  // [2] File with laser data
   status = daqDA_FES_storeFile(LASDATA_FILE,LASDATA_FILE);
   if(status){
     printf("Failed to export file : %d\n",status);
     return -1;
   }
+  // [3] File with laser histos
+  status = daqDA_FES_storeFile(LASHISTO_FILE,LASHISTO_FILE);
+  if(status){
+    printf("Failed to export pedestal histos file to DAQ FES\n");
+    return -1;
+  }
 
   /* report progress */
   daqDA_progressReport(100);
index a3dcd8c..f96bd5e 100644 (file)
@@ -20,6 +20,7 @@ Trigger Types Used: Standalone Trigger
 
 */
 #define PEDDATA_FILE  "ZDCPedestal.dat"
+#define PEDHISTO_FILE "ZDCPedHisto.root"
 #define MAPDATA_FILE  "ZDCChMapping.dat"
 
 #include <stdio.h>
@@ -32,13 +33,15 @@ Trigger Types Used: Standalone Trigger
 #include <daqDA.h>
 
 //ROOT
-#include <TRandom.h>
+#include <TROOT.h>
+#include <TPluginManager.h>
 #include <TH1F.h>
 #include <TH2F.h>
 #include <TProfile.h>
 #include <TF1.h>
 #include <TFile.h>
 #include <TFitter.h>
+#include "TMinuitMinimizer.h"
 
 //AliRoot
 #include <AliRawReaderDate.h>
@@ -51,7 +54,17 @@ Trigger Types Used: Standalone Trigger
 */
 int main(int argc, char **argv) {
   
-//  TVirtualFitter::SetDefaultFitter("Minuit2");
+
+  gROOT->GetPluginManager()->AddHandler("TVirtualStreamerInfo",
+                                       "*",
+                                       "TStreamerInfo",
+                                       "RIO",
+                                       "TStreamerInfo()"); 
+
+  TMinuitMinimizer m; 
+  gROOT->GetPluginManager()->AddHandler("ROOT::Math::Minimizer", "Minuit","TMinuitMinimizer",
+      "Minuit", "TMinuitMinimizer(const char *)");
+  TVirtualFitter::SetDefaultFitter("Minuit");
 
   int status = 0;
   int const kNChannels = 24;
@@ -60,7 +73,7 @@ int main(int argc, char **argv) {
   printf("\n ZDC PEDESTAL program started\n");  
 
   /* check that we got some arguments = list of files */
-  if (argc<2) {
+  if (argc<2){
     printf("Wrong number of arguments\n");
     return -1;
   }
@@ -405,7 +418,7 @@ int main(int argc, char **argv) {
   
   // --- Correlations
 
-  Float_t CorrCoeff0[2*kNChannels], CorrCoeff1[2*kNChannels];
+/*  Float_t CorrCoeff0[2*kNChannels], CorrCoeff1[2*kNChannels];
   TProfile *hPedCorrProfhg[kNChannels], *hPedCorrProflg[kNChannels];
   TF1 *ffunchg[kNChannels], *ffunclg[kNChannels];
   char namhist4[50];
@@ -432,10 +445,26 @@ int main(int argc, char **argv) {
      //printf("\t CorrCoeff0[%d] = %f, CorrCoeff1[%d] = %f\n",
      //                i+kNChannels, CorrCoeff0[i+kNChannels], i+kNChannels, CorrCoeff1[i+kNChannels]);
   }    
-
+*/
   //                                                  
   fclose(fileShuttle);
   //
+  /* report progress */
+  daqDA_progressReport(80);
+  //
+  TFile *histofile = new TFile(PEDHISTO_FILE,"RECREATE");
+  histofile->cd();
+  for(int k=0; k<kNChannels; k++){
+     hPedhg[k]->Write();
+     hPedOutOfTimehg[k]->Write();
+     hPedCorrhg[k]->Write();
+     hPedlg[k]->Write();
+     hPedOutOfTimelg[k]->Write();
+     hPedCorrlg[k]->Write();
+  }
+  //                                                  
+  histofile->Close();
+  //
   for(Int_t j=0; j<kNChannels; j++){
      delete hPedhg[j];
      delete hPedOutOfTimehg[j];
@@ -458,20 +487,26 @@ int main(int argc, char **argv) {
   // [1] File with mapping
   status = daqDA_FES_storeFile(MAPDATA_FILE,MAPDATA_FILE);
   if(status){
-    printf("Failed to export file %d to DAQ FES\n",status);
+    printf("Failed to export mapping data file to DAQ FES\n");
     return -1;
   }
   // [2] File with pedestal data
   status = daqDA_FES_storeFile(PEDDATA_FILE,PEDDATA_FILE);
   if(status){
-    printf("Failed to export file %d to DAQ FES\n",status);
+    printf("Failed to export pedestal data file to DAQ FES\n");
+    return -1;
+  }
+  // [3] File with pedestal histos
+  status = daqDA_FES_storeFile(PEDHISTO_FILE,PEDHISTO_FILE);
+  if(status){
+    printf("Failed to export pedestal histos file to DAQ FES\n");
     return -1;
   }
   
   /* store the result files on DB */
   status = daqDA_DB_storeFile(PEDDATA_FILE,PEDDATA_FILE);  
   if(status){
-    printf("Failed to export file %d to DAQ DB\n",status);
+    printf("Failed to store pedestal data file to DAQ DB\n");
     return -1;
   }