Make coherent ouput
authorschutz <schutz@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 4 May 2007 17:41:11 +0000 (17:41 +0000)
committerschutz <schutz@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 4 May 2007 17:41:11 +0000 (17:41 +0000)
ESDCheck/AliAnalysisTaskPt.cxx
ESDCheck/AliEMCALQATask.cxx
ESDCheck/AliHMPIDQATask.cxx
ESDCheck/AliMUONQATask.cxx
ESDCheck/AliPHOSQATask.cxx
ESDCheck/AliPMDQATask.cxx
ESDCheck/AliT0QATask.cxx
ESDCheck/AliTOFQATask.cxx
ESDCheck/AliTRDQATask.cxx
ESDCheck/AliVZEROQATask.cxx
ESDCheck/ana.C

index fe3f16d..ab95280 100644 (file)
@@ -120,14 +120,17 @@ void AliAnalysisTaskPt::Exec(Option_t *)
 void AliAnalysisTaskPt::Terminate(Option_t *) 
 {
   // Processing when the event loop is ended
-  
+
+  AliInfo(Form(" *** %s Report:", GetName())) ; 
+
   TCanvas *c1 = new TCanvas("c1","Pt",10,10,310,310);
   c1->SetFillColor(10);
   c1->SetHighLightColor(10);
   
   c1->cd(1)->SetLeftMargin(0.15);
   c1->cd(1)->SetBottomMargin(0.15);  
-  c1->cd(1)->SetLogy();
+  if (fhPt->GetMaximum() > 0 ) 
+    c1->cd(1)->SetLogy();
   fOutputContainer = (TObjArray*)GetOutputData(0);
   fhPt = (TH1F*)fOutputContainer->At(0);
   if (fhPt) fhPt->DrawCopy("E");
@@ -135,7 +138,7 @@ void AliAnalysisTaskPt::Terminate(Option_t *)
   c1->Print("TracksPt.eps");
 
   char line[1024] ; 
-  sprintf(line, ".!tar -zcvf %s.tar.gz *.eps", GetName()) ; 
+  sprintf(line, ".!tar -zcf %s.tar.gz *.eps", GetName()) ; 
   gROOT->ProcessLine(line);
   sprintf(line, ".!rm -fR *.eps"); 
   gROOT->ProcessLine(line);
index b3d9c8e..045aa7f 100644 (file)
@@ -210,36 +210,42 @@ void AliEMCALQATask::Terminate(Option_t *)
   fhEMCALPhotons = (TH1I*)fOutputContainer->At(5);
   fhEMCALInvariantMass = (TH1D*)fOutputContainer->At(6);
   fhEMCALDigitsEvent = (TH1I*)fOutputContainer->At(7);
-  printf("EMCALEnergy Mean        : %5.3f , RMS : %5.3f \n", fhEMCALEnergy->GetMean(),        fhEMCALEnergy->GetRMS()        ) ;
-  printf("EMCALDigits Mean        : %5.3f , RMS : %5.3f \n", fhEMCALDigits->GetMean(),        fhEMCALDigits->GetRMS()        ) ;
-  printf("EMCALRecParticles Mean  : %5.3f , RMS : %5.3f \n", fhEMCALRecParticles->GetMean(),  fhEMCALRecParticles->GetRMS()  ) ;
-  printf("EMCALPhotons Mean       : %5.3f , RMS : %5.3f \n", fhEMCALPhotons->GetMean(),       fhEMCALPhotons->GetRMS()       ) ;
-  printf("EMCALInvariantMass Mean : %5.3f , RMS : %5.3f \n", fhEMCALInvariantMass->GetMean(), fhEMCALInvariantMass->GetRMS() ) ;
-  printf("EMCALDigitsEvent Mean   : %5.3f , RMS : %5.3f \n", fhEMCALDigitsEvent->GetMean(),   fhEMCALDigitsEvent->GetRMS()   ) ;
+  AliInfo(Form(" *** %s Report:", GetName())) ; 
+  printf("        EMCALEnergy Mean        : %5.3f , RMS : %5.3f \n", fhEMCALEnergy->GetMean(),        fhEMCALEnergy->GetRMS()        ) ;
+  printf("        EMCALDigits Mean        : %5.3f , RMS : %5.3f \n", fhEMCALDigits->GetMean(),        fhEMCALDigits->GetRMS()        ) ;
+  printf("        EMCALRecParticles Mean  : %5.3f , RMS : %5.3f \n", fhEMCALRecParticles->GetMean(),  fhEMCALRecParticles->GetRMS()  ) ;
+  printf("        EMCALPhotons Mean       : %5.3f , RMS : %5.3f \n", fhEMCALPhotons->GetMean(),       fhEMCALPhotons->GetRMS()       ) ;
+  printf("        EMCALInvariantMass Mean : %5.3f , RMS : %5.3f \n", fhEMCALInvariantMass->GetMean(), fhEMCALInvariantMass->GetRMS() ) ;
+  printf("        EMCALDigitsEvent Mean   : %5.3f , RMS : %5.3f \n", fhEMCALDigitsEvent->GetMean(),   fhEMCALDigitsEvent->GetRMS()   ) ;
 
   TCanvas  * cEMCAL = new TCanvas("EMCAL", "EMCAL ESD Test", 400, 10, 600, 700);
   cEMCAL->Divide(3, 2) ; 
 
   cEMCAL->cd(1) ; 
-  gPad->SetLogy();
+  if ( fhEMCALEnergy->GetMaximum() > 0. ) 
+    gPad->SetLogy();
   fhEMCALEnergy->SetAxisRange(0, 25.);
   fhEMCALEnergy->SetXTitle("Energy (GeV)");
   fhEMCALEnergy->Draw();
   
   cEMCAL->cd(2) ; 
-  gPad->SetLogy();
+  if ( fhEMCALDigits->GetMaximum() > 0. ) 
+    gPad->SetLogy();
   fhEMCALDigits->SetAxisRange(0, 25.);
   fhEMCALDigits->SetXTitle("DigitsPerCluster");
   fhEMCALDigits->Draw();
  
   cEMCAL->cd(3) ; 
-  gPad->SetLogy();
+  if ( fhEMCALRecParticles->GetMaximum() > 0. ) 
+     gPad->SetLogy();
   fhEMCALRecParticles->SetAxisRange(0, 25.);
   fhEMCALRecParticles->SetXTitle("RecParticles");
   fhEMCALRecParticles->Draw();
  
   cEMCAL->cd(4) ; 
-  gPad->SetLogy();
+  if ( fhEMCALPhotons->GetMaximum() > 0. ) 
+    gPad->SetLogy();
   fhEMCALPhotons->SetAxisRange(0, 25.);
   fhEMCALPhotons->SetXTitle("Photons");
   fhEMCALPhotons->Draw();
@@ -249,7 +255,8 @@ void AliEMCALQATask::Terminate(Option_t *)
   fhEMCALInvariantMass->Draw();
  
   cEMCAL->cd(6) ; 
-  gPad->SetLogy();
+  if ( fhEMCALDigitsEvent->GetMaximum() > 0. ) 
+    gPad->SetLogy();
   fhEMCALDigitsEvent->SetAxisRange(0, 40.);
   fhEMCALDigitsEvent->SetXTitle("DigitsPerEvent");
   fhEMCALDigitsEvent->Draw();
@@ -257,7 +264,7 @@ void AliEMCALQATask::Terminate(Option_t *)
   cEMCAL->Print("EMCAL.eps");
  
   char line[1024] ; 
-  sprintf(line, ".!tar -zcvf %s.tar.gz *.eps", GetName()) ; 
+  sprintf(line, ".!tar -zcf %s.tar.gz *.eps", GetName()) ; 
   gROOT->ProcessLine(line);
   sprintf(line, ".!rm -fR *.eps"); 
   gROOT->ProcessLine(line);
index 839966c..2e0b6e2 100644 (file)
@@ -183,7 +183,9 @@ void AliHMPIDQATask::Terminate(Option_t *)
   fhHMPIDProb[2] = (TH1F*)fOutputContainer->At(6);
   fhHMPIDProb[3] = (TH1F*)fOutputContainer->At(7);
   fhHMPIDProb[4] = (TH1F*)fOutputContainer->At(8);
-  
+
+  AliInfo(Form(" *** %s Report:", GetName())) ; 
+
   Float_t n = 1.292 ; //mean freon ref idx 
   TF1 * hHMPIDpPi = new TF1("RiPiTheo", "acos(sqrt(x*x+[0]*[0])/(x*[1]))", 1.2, 7) ; 
   hHMPIDpPi->SetLineWidth(1) ; 
@@ -233,7 +235,7 @@ void AliHMPIDQATask::Terminate(Option_t *)
   cHMPID->Print("HMPID.eps");
   
   char line[1024] ; 
-  sprintf(line, ".!tar -zcvf %s.tar.gz *.eps", GetName()) ; 
+  sprintf(line, ".!tar -zcf %s.tar.gz *.eps", GetName()) ; 
   gROOT->ProcessLine(line);
   sprintf(line, ".!rm -fR *.eps"); 
   gROOT->ProcessLine(line);
index 949d2e3..0d9f4fe 100644 (file)
@@ -145,7 +145,7 @@ void AliMUONQATask::CreateOutputObjects()
   fhChi2match = new TH1F("hChi2match","Chi2 of trig/track matching",100, 0.,20.);
   // create output container
   
-  fOutputContainer = new TObjArray(10) ; 
+  fOutputContainer = new TObjArray(12) ; 
   fOutputContainer->SetName(GetName()) ; 
   fOutputContainer->AddAt(fhMUONVertex,             0) ; 
   fOutputContainer->AddAt(fhMUONMult,               1) ; 
@@ -279,7 +279,8 @@ void AliMUONQATask::Terminate(Option_t *)
       effMatch = 100 * fnTrackTrig / ftracktot ;
       fheffMatchT->Fill(effMatch);}
 
-  AliInfo(Form("Terminate %s:", GetName())) ;
+  AliInfo(Form(" *** %s Report:", GetName())) ; 
+
   fOutputContainer = (TObjArray*)GetOutputData(0);
   fhMUONVertex = (TH1F*)fOutputContainer->At(0);
   fhMUONMult   = (TH1F*)fOutputContainer->At(1); 
@@ -291,30 +292,24 @@ void AliMUONQATask::Terminate(Option_t *)
   fhUSHighpt=(TH1F*)fOutputContainer->At(7);
   fhLSLowpt=(TH1F*)fOutputContainer->At(8); 
   fhLSHighpt=(TH1F*)fOutputContainer->At(9);
- //______________________________________________________________________
-  printf("===================================================\n") ;
-  printf("================  %s ESD SUMMARY    ==============\n", GetName()) ;
-  printf("                                                   \n") ;
+
   printf("         Total number of processed events  %d      \n", fnevents) ;
-  printf("\n")  ;
-  printf("\n")  ;
-  printf("Table 1:                                         \n") ;
-   printf("===================================================\n") ;
- printf(" Global Trigger output       Low pt  High pt \n") ;
-  printf(" number of Single      :\t");
-  printf("%i\t", fSLowpt) ;
+  printf("     \n")  ;
+  printf("     \n")  ;
+  printf("     Table 1:                                         \n") ;
+   printf("    ===================================================\n") ;
+ printf("      Global Trigger output       Low pt  High pt \n") ;
+  printf("     number of Single      :\t");
+  printf("     %i\t", fSLowpt) ;
   printf("\n");
-  printf(" number of UnlikeSign pair  :\t"); 
-  printf("%i\t%i\t", fUSLowpt, fUSHighpt) ;
+  printf("     number of UnlikeSign pair  :\t"); 
+  printf("     %i\t%i\t", fUSLowpt, fUSHighpt) ;
   printf("\n");
-  printf(" number of LikeSign pair    :\t");  
-  printf("%i\t%i\t", fLSLowpt, fLSHighpt) ;
+  printf("     number of LikeSign pair    :\t");  
+  printf("     %i\t%i\t", fLSLowpt, fLSHighpt) ;
   printf("\n");
-  printf("===================================================\n") ;
-  printf("\n") ;
-  printf("matching efficiency with the trigger for single tracks = %2d %% \n", effMatch);
+  printf("     matching efficiency with the trigger for single tracks = %2d %% \n", effMatch);
   printf("\n") ;
-  printf("===================================================\n") ;
   
 TCanvas * cMUON1 = new TCanvas("cMUON1", "MUON ESD Vert & Mult", 400, 10, 600, 700) ;
   cMUON1->Divide(1,2) ;
@@ -369,7 +364,7 @@ TCanvas * cMUON1 = new TCanvas("cMUON1", "MUON ESD Vert & Mult", 400, 10, 600, 7
   cMUON4->Print("MUON4.eps") ;
  
   char line[1024] ; 
-  sprintf(line, ".!tar -zcvf %s.tar.gz *.eps", GetName()) ; 
+  sprintf(line, ".!tar -zcf %s.tar.gz *.eps", GetName()) ; 
   gROOT->ProcessLine(line);
   sprintf(line, ".!rm -fR *.eps"); 
   gROOT->ProcessLine(line);
index 94e3382..a2dfdd9 100644 (file)
@@ -212,18 +212,20 @@ void AliPHOSQATask::Terminate(Option_t *)
   fhPHOSInvariantMass  = (TH1D*)fOutputContainer->At(6);
   fhPHOSDigitsEvent    = (TH1I*)fOutputContainer->At(7);
 
-  printf("PHOSEnergy Mean         : %5.3f , RMS : %5.3f \n", fhPHOSEnergy->GetMean(),         fhPHOSEnergy->GetRMS()         ) ;
-  printf("PHOSDigits Mean         : %5.3f , RMS : %5.3f \n", fhPHOSDigits->GetMean(),         fhPHOSDigits->GetRMS()         ) ;
-  printf("PHOSRecParticles Mean   : %5.3f , RMS : %5.3f \n", fhPHOSRecParticles->GetMean(),   fhPHOSRecParticles->GetRMS()   ) ;
-  printf("PHOSPhotons Mean        : %5.3f , RMS : %5.3f \n", fhPHOSPhotons->GetMean(),        fhPHOSPhotons->GetRMS()        ) ;
-  printf("PHOSInvariantMass Mean  : %5.3f , RMS : %5.3f \n", fhPHOSInvariantMass->GetMean(),  fhPHOSInvariantMass->GetRMS()  ) ;
-  printf("PHOSDigitsEvent Mean    : %5.3f , RMS : %5.3f \n", fhPHOSDigitsEvent->GetMean(),    fhPHOSDigitsEvent->GetRMS()    ) ;
+  AliInfo(Form(" *** %s Report:", GetName())) ; 
+  printf("        PHOSEnergy Mean         : %5.3f , RMS : %5.3f \n", fhPHOSEnergy->GetMean(),         fhPHOSEnergy->GetRMS()         ) ;
+  printf("        PHOSDigits Mean         : %5.3f , RMS : %5.3f \n", fhPHOSDigits->GetMean(),         fhPHOSDigits->GetRMS()         ) ;
+  printf("        PHOSRecParticles Mean   : %5.3f , RMS : %5.3f \n", fhPHOSRecParticles->GetMean(),   fhPHOSRecParticles->GetRMS()   ) ;
+  printf("        PHOSPhotons Mean        : %5.3f , RMS : %5.3f \n", fhPHOSPhotons->GetMean(),        fhPHOSPhotons->GetRMS()        ) ;
+  printf("        PHOSInvariantMass Mean  : %5.3f , RMS : %5.3f \n", fhPHOSInvariantMass->GetMean(),  fhPHOSInvariantMass->GetRMS()  ) ;
+  printf("        PHOSDigitsEvent Mean    : %5.3f , RMS : %5.3f \n", fhPHOSDigitsEvent->GetMean(),    fhPHOSDigitsEvent->GetRMS()    ) ;
 
   TCanvas  * cPHOS = new TCanvas("cPHOS", "PHOS ESD Test", 400, 10, 600, 700) ;
   cPHOS->Divide(3, 2);
 
   cPHOS->cd(1) ; 
-  gPad->SetLogy();
+  if ( fhPHOSEnergy->GetMaximum() > 0. ) 
+    gPad->SetLogy();
   fhPHOSEnergy->SetAxisRange(0, 25.);
   fhPHOSEnergy->SetLineColor(2);
   fhPHOSEnergy->Draw();
@@ -234,13 +236,15 @@ void AliPHOSQATask::Terminate(Option_t *)
   fhPHOSDigits->Draw();
 
   cPHOS->cd(3) ; 
-  gPad->SetLogy();
+  if ( fhPHOSRecParticles->GetMaximum() > 0. ) 
+    gPad->SetLogy();
   fhPHOSRecParticles->SetAxisRange(0, 25.);
   fhPHOSRecParticles->SetLineColor(2);
   fhPHOSRecParticles->Draw();
 
   cPHOS->cd(4) ; 
-  gPad->SetLogy();
+  if ( fhPHOSPhotons->GetMaximum() > 0. ) 
+    gPad->SetLogy();
   fhPHOSPhotons->SetAxisRange(0,25.);
   fhPHOSPhotons->SetLineColor(2);
   fhPHOSPhotons->Draw();
@@ -250,7 +254,8 @@ void AliPHOSQATask::Terminate(Option_t *)
   fhPHOSInvariantMass->Draw();
  
   cPHOS->cd(6) ; 
-  gPad->SetLogy();
+  if ( fhPHOSDigitsEvent->GetMaximum() > 0. ) 
+    gPad->SetLogy();
   fhPHOSDigitsEvent->SetAxisRange(0,40.);
   fhPHOSDigitsEvent->SetLineColor(2);
   fhPHOSDigitsEvent->Draw();
@@ -258,7 +263,7 @@ void AliPHOSQATask::Terminate(Option_t *)
   cPHOS->Print("PHOS.eps");
  
   char line[1024] ; 
-  sprintf(line, ".!tar -zcvf %s.tar.gz *.eps", GetName()) ; 
+  sprintf(line, ".!tar -zcf %s.tar.gz *.eps", GetName()) ; 
   gROOT->ProcessLine(line);
   sprintf(line, ".!rm -fR *.eps"); 
   gROOT->ProcessLine(line);
index ecceb3f..91a4c93 100644 (file)
@@ -474,6 +474,8 @@ void AliPMDQATask::Terminate(Option_t *)
   fhPMDSM24 = (TH2F*)fOutputContainer->At(36);
   fhPMDSM   = (TH1F*)fOutputContainer->At(37);
 
+  AliInfo(Form(" *** %s Report:", GetName())) ; 
+
   gStyle->SetOptStat(110000);
   gStyle->SetOptFit(1);
 
@@ -581,7 +583,7 @@ void AliPMDQATask::Terminate(Option_t *)
   cPMD4->Print("CPVPREAdc.eps");
 
   char line[1024] ; 
-  sprintf(line, ".!tar -zcvf %s.tar.gz *.eps", GetName()) ; 
+  sprintf(line, ".!tar -zcf %s.tar.gz *.eps", GetName()) ; 
   gROOT->ProcessLine(line);
   
   AliInfo(Form("!!! All the eps files are in %s.tar.gz !!! \n", GetName())) ;
index a0353cf..451ef5a 100644 (file)
@@ -142,6 +142,8 @@ void AliT0QATask::Terminate(Option_t *)
   fhT02 = (TH1F*)fOutputContainer->At(1);
   fhT03 = (TH1F*)fOutputContainer->At(2);
 
+  AliInfo(Form(" *** %s Report:", GetName())) ; 
+
   Float_t mean = fhT02->GetMean();
 
   printf ("mean time T0 ps %f\n", mean) ;
@@ -165,7 +167,7 @@ void AliT0QATask::Terminate(Option_t *)
   cTO1->Print("T0.eps");
 
   char line[1024] ; 
-  sprintf(line, ".!tar -zcvf %s.tar.gz *.eps", GetName()) ; 
+  sprintf(line, ".!tar -zcf %s.tar.gz *.eps", GetName()) ; 
   gROOT->ProcessLine(line);
   sprintf(line, ".!rm -fR *.eps"); 
   gROOT->ProcessLine(line);
index 60c7f92..5a1ff1d 100644 (file)
@@ -326,7 +326,7 @@ void AliTOFQATask::Exec(Option_t *)
   fhTOFMatch->Fill(fracM);
 
   PostData(0, fOutputContainer);  
-  AliInfo("Finishing event processing...") ; 
+  // AliInfo("Finishing event processing...") ; 
 
 }
 //______________________________________________________________________________
@@ -468,7 +468,6 @@ void AliTOFQATask::DrawHistos()
   fhTOFeffMom->Draw();
   fhTOFeffMom->Draw("histo,same");
   
-  cTOFeff->Print("TOF_eff.gif");
   cTOFeff->Print("TOF_eff.eps");
 
 
@@ -476,31 +475,34 @@ void AliTOFQATask::DrawHistos()
   TCanvas * cTOFtime = new TCanvas("cTOFtime", "TOF measured Times ", 400, 30, 550, 630) ;  
   cTOFtime->Divide(1,2) ;
   cTOFtime->cd(1);
-  cTOFtime->GetPad(1)->SetLogy(1);
+  if ( fhTOFTime->GetMaximum() > 0 )  
+    cTOFtime->GetPad(1)->SetLogy(1);
   fhTOFTime->GetXaxis()->SetTitle("TOF time (ns)");
   fhTOFTime->GetYaxis()->SetTitle("Entries");
   fhTOFTime->SetFillColor(4);
   fhTOFTime->Draw();
   cTOFtime->cd(2);
-  cTOFtime->GetPad(2)->SetLogy(1);
+  if ( fhTOFDeltaTime->GetMaximum() > 0 )  
+    cTOFtime->GetPad(2)->SetLogy(1);
   fhTOFDeltaTime->GetXaxis()->SetTitle("t^{TOF}-t^{exp}_{#pi} (ns)");
   fhTOFDeltaTime->GetYaxis()->SetTitle("Entries");
   fhTOFDeltaTime->SetFillColor(4);
   fhTOFDeltaTime->Draw();
 
-  cTOFtime->Print("TOF_time.gif");
   cTOFtime->Print("TOF_time.eps");
 
 
   TCanvas * cTOFpid = new TCanvas("cTOFpid", "TOF PID ", 400, 30, 550, 630) ;  
   cTOFpid->Divide(1,3) ;
   cTOFpid->cd(1);
-  cTOFpid->SetLogy(1);
+  if ( fhTOFMass->GetMaximum() > 0 )  
+    cTOFpid->SetLogy(1);
   fhTOFMass->GetXaxis()->SetTitle("Reconstructed Mass (GeV/c^{2})");
   fhTOFMass->GetYaxis()->SetTitle("Entries");
   fhTOFMass->SetFillColor(4);
   fhTOFMass->Draw();
-  cTOFpid->SetLogy(0);
+  if ( fhTOFMassVsMom->GetMaximum() > 0 )  
+    cTOFpid->SetLogy(0);
   cTOFpid->cd(2);
   fhTOFMassVsMom->GetYaxis()->SetRange(0,400);
   fhTOFMassVsMom->GetXaxis()->SetTitle("Reconstructed Mass (GeV/c^{2})");
@@ -555,7 +557,6 @@ void AliTOFQATask::DrawHistos()
   texpr->Draw();
 
 
-  cTOFpid->Print("TOF_pid.gif");
   cTOFpid->Print("TOF_pid.eps");
 
   // draw all 
@@ -572,7 +573,6 @@ void AliTOFQATask::Terminate(Option_t *)
   
   // some plots
 
-  AliInfo("TOF QA Task: End of events loop");
   fOutputContainer = (TObjArray*)GetOutputData(0);
   fhTOFMatch       = (TH1F*)fOutputContainer->At(0);
   fhESDeffPhi      = (TH1F*)fOutputContainer->At(1);
@@ -592,7 +592,16 @@ void AliTOFQATask::Terminate(Option_t *)
   fhTOFIDSpecies   = (TH1F*)fOutputContainer->At(15);
   fhTOFMassVsMom   = (TH2F*)fOutputContainer->At(16);
   fhTOFMass        = (TH1F*)fOutputContainer->At(17);
+
+  AliInfo(Form(" *** %s Report:", GetName())) ; 
+
   GetEfficiency();
 //  PostData(0, fOutputContainer);
   DrawHistos() ; 
+
+  char line[1024] ; 
+  sprintf(line, ".!tar -zcf %s.tar.gz *.eps", GetName()) ; 
+  gROOT->ProcessLine(line);
+  
+  AliInfo(Form("!!! All the eps files are in %s.tar.gz !!! \n", GetName())) ;
 }
index 196a314..e0798e9 100644 (file)
@@ -212,7 +212,7 @@ void AliTRDQATask::CreateOutputObjects()
      fOutputContainer->AddAt(fTrdSigMomPID[i], counter++);
   }
 
-  AliInfo(Form("Number of histograms = %d", counter));
+  //AliInfo(Form("Number of histograms = %d", counter));
 
  }
 
@@ -352,7 +352,6 @@ void AliTRDQATask::Exec(Option_t *)
 void AliTRDQATask::Terminate(Option_t *)
 {
   // Processing when the event loop is ended
-  AliInfo("TRD QA module");
   fOutputContainer = (TObjArray*)GetOutputData(0);
   int counter = 0;
   fNTracks     = (TH1D*)fOutputContainer->At(counter++);
@@ -399,14 +398,21 @@ void AliTRDQATask::Terminate(Option_t *)
   }
 
   // create efficiency histograms
+  AliInfo(Form(" *** %s Report:", GetName())) ; 
   
   CalculateEff();
-//  PostData(0, fOutputContainer);
 
   DrawESD() ; 
   DrawGeoESD() ; 
   //DrawConvESD() ; 
   DrawPidESD() ; 
+
+  char line[1024] ; 
+  sprintf(line, ".!tar -zcf %s.tar.gz *.eps", GetName()) ; 
+  gROOT->ProcessLine(line);
+  
+  AliInfo(Form("!!! All the eps files are in %s.tar.gz !!! \n", GetName())) ;
+
 }
 
 //______________________________________________________________________________
@@ -457,8 +463,6 @@ void AliTRDQATask::DrawESD()
 {
   // Makes a few plots
 
-  AliInfo("Plotting....") ; 
-  
   TCanvas * cTRD = new TCanvas("cTRD", "TRD ESD Test", 400, 10, 600, 700) ;
   cTRD->Divide(6,3) ;
 
@@ -500,12 +504,12 @@ void AliTRDQATask::DrawESD()
   cTRD->cd(i+1) ;
   
   //  new TCanvas(names[i], names[nhist], 500, 300);
-    gPad->SetLogy(logy[i]);
       
     for(int j=0; j<nover[i]; j++) {
       TH1D *hist = dynamic_cast<TH1D*>(gDirectory->FindObject(names[nhist++]));
       if (!hist) continue;
-      
+      if (hist->GetMaximum() > 0 ) 
+       gPad->SetLogy(logy[i]);
       if (strstr(hist->GetName(), "eff")) {
        hist->SetMarkerStyle(20);
        hist->SetMinimum(0);
@@ -516,7 +520,7 @@ void AliTRDQATask::DrawESD()
       else hist->Draw("SAME");
     }
   }
-  cTRD->Print("TRD_ESD.gif");
+  cTRD->Print("TRD_ESD.eps");
 }
 
 //______________________________________________________________________________
@@ -524,8 +528,6 @@ void AliTRDQATask::DrawGeoESD()
 {
   // Makes a few plots
 
-  AliInfo("Plotting....") ; 
-
   TCanvas * cTRDGeo = new TCanvas("cTRDGeo", "TRD ESDGeo Test", 400, 10, 600, 700) ;
   cTRDGeo->Divide(4,2) ;
 
@@ -563,14 +565,15 @@ void AliTRDQATask::DrawGeoESD()
     //if (i<2) new TCanvas(names[i], names[i], 500, 300);
     //else new TCanvas(names[i], names[i], 300, 900);
    
-    gPad->SetLogy(logy[i]);
+    if (hist->GetMaximum() > 0 ) 
+      gPad->SetLogy(logy[i]);
     if (strstr(opt[i],"colz")) gPad->SetRightMargin(0.1);
     
     hist->Draw(opt[i]);    
     AliInfo(Form("%s\t%d", names[i], hist->GetEntries()));
   }
   
-  cTRDGeo->Print("TRD_Geo.gif");
+  cTRDGeo->Print("TRD_Geo.eps");
 }
 
 //______________________________________________________________________________
@@ -613,11 +616,12 @@ void AliTRDQATask::DrawConvESD()
   for(int i=0; i<nplots; i++) {
     cTRDConv->cd(i+1) ;
     //new TCanvas(names[i], names[i], 500, 300);
-    gPad->SetLogy(logy[i]);
     if (strstr(opt[i],"colz")) gPad->SetRightMargin(0.1);
    
     for(int j=0; j<nover[i]; j++) {
       TH1D *hist = dynamic_cast<TH1D*>(gDirectory->FindObject(names[nhist++]));
+      if ( hist->GetMaximum() > 0 ) 
+       gPad->SetLogy(logy[i]);
       if (!j) hist->Draw(opt[i]);
       else hist->Draw("same");
     }
@@ -631,7 +635,6 @@ void AliTRDQATask::DrawPidESD()
 {
   // Makes a few plots
 
-  AliInfo("Plotting....") ; 
   TCanvas * cTRDPid = new TCanvas("cTRDPid", "TRD ESDPid Test", 400, 10, 600, 700) ;
   cTRDPid->Divide(9,3) ;
 
@@ -692,7 +695,8 @@ void AliTRDQATask::DrawPidESD()
     if (!hist) continue;
     
     //new TCanvas(names[i], names[i], 500, 300);
-    gPad->SetLogy(logy[i]);
+    if ( hist->GetMaximum() > 0  ) 
+      gPad->SetLogy(logy[i]);
     if (strstr(opt[i],"colz")) gPad->SetRightMargin(0.1);
     
     if (strstr(names[i],"sigMom")) gPad->SetLogz(1);
@@ -701,5 +705,5 @@ void AliTRDQATask::DrawPidESD()
     hist->Draw(opt[i]);    
     AliInfo(Form("%s\t%d", names[i], hist->GetEntries()));
   }
-   cTRDPid->Print("TRD_Pid.gif");
+   cTRDPid->Print("TRD_Pid.eps");
 }
index 995e619..cbce7da 100644 (file)
@@ -150,8 +150,9 @@ void AliVZEROQATask::Terminate(Option_t *)
   fhVZEROMultA     = (TH1I*)fOutputContainer->At(2);
   fhVZEROMultC     = (TH1I*)fOutputContainer->At(3);
   
-  printf("V0A Multiplicity Mean : %5.3f , RMS : %5.3f \n",fhVZEROMultA->GetMean(),fhVZEROMultA->GetRMS());
-  printf("V0C Multiplicity Mean : %5.3f , RMS : %5.3f \n",fhVZEROMultC->GetMean(),fhVZEROMultC->GetRMS());
+  AliInfo(Form(" *** %s Report:", GetName())) ; 
+  printf("        V0A Multiplicity Mean : %5.3f , RMS : %5.3f \n",fhVZEROMultA->GetMean(),fhVZEROMultA->GetRMS());
+  printf("        V0C Multiplicity Mean : %5.3f , RMS : %5.3f \n",fhVZEROMultC->GetMean(),fhVZEROMultC->GetRMS());
 
   TCanvas  * c1 = new TCanvas("Number of PMs fired in V0A", "Number of PMs fired in V0A", 1);
   fhVZERONbPMA->SetAxisRange(0, 99);
@@ -183,7 +184,7 @@ void AliVZEROQATask::Terminate(Option_t *)
   c4->Print("NumberV0CPMs.eps");
   
   char line[1024] ; 
-  sprintf(line, ".!tar -zcvf %s.tar.gz *.eps", GetName()) ; 
+  sprintf(line, ".!tar -zcf %s.tar.gz *.eps", GetName()) ; 
   gROOT->ProcessLine(line);
   sprintf(line, ".!rm -fR *.eps"); 
   gROOT->ProcessLine(line);
index acbfc23..d6f22ee 100644 (file)
@@ -36,10 +36,10 @@ Bool_t LoadLib( const char* pararchivename)
   }
 
   if ( strstr(pararchivename, "ESD") ) {
-    gSystem->Load("libVMC.so");
+    //gSystem->Load("libVMC.so");
     gSystem->Load("libESD.so");
-    gSystem->Load("libRAliEn.so") ;
-    gSystem->Load("libProof.so") ;
+    //gSystem->Load("libRAliEn.so") ;
+    //gSystem->Load("libProof.so") ;
   }
 
   printf("*** %s library loaded *** %s **\n", pararchivename);
@@ -56,7 +56,6 @@ void ana()
   if (! gIsAnalysisLoaded ) {
     LoadLib("ESD") ; 
     LoadLib("ANALYSIS") ; 
-    printf("Include path = %s\n", gSystem->GetIncludePath()) ; 
     LoadLib("AnalysisCheck") ; 
   }
   
@@ -120,7 +119,11 @@ void ana()
 
   // definition of Tag cuts 
   const char * runCuts = 0x0 ; 
-  const char * evtCuts = 0x0 ; //"fEventTag.fNPHOSClustersMin == 1 && fEventTag.fNEMCALClustersMin == 1" ; 
+  const char * evtCuts = 0x0 ; 
+  const char * lhcCuts = 0x0 ; 
+  const char * detCuts = 0x0 ; 
+  
+//"fEventTag.fNPHOSClustersMin == 1 && fEventTag.fNEMCALClustersMin == 1" ; 
 
   
   TString input = gSystem->Getenv("ANA_INPUT") ; 
@@ -130,7 +133,7 @@ void ana()
       //create the ESD collection from the tag collection 
       input.ReplaceAll("tag?", "") ; 
       const char * collESD = "esdCollection.xml" ;
-      ag->MakeEsdCollectionFromTagCollection(runCuts, evtCuts, input.Data(), collESD) ;
+      ag->MakeEsdCollectionFromTagCollection(runCuts, lhcCuts, detCuts, evtCuts, input.Data(), collESD) ;
       sprintf(argument, "esd?%s", collESD) ; 
     } else if ( input.Contains("esd?") ) 
       sprintf(argument, "%s", input.Data()) ; 
@@ -141,7 +144,7 @@ void ana()
     TChain* analysisChain = new TChain("esdTree") ;
     //   input = "alien:///alice/cern.ch/user/a/aliprod/prod2006_2/output_pp/105/411/AliESDs.root" ; 
     //   analysisChain->AddFile(input);
-    input = "/home/schutz/group/schutz/work/analysis/QA/HEAD/RunV0/AliESDs.root" ; 
+    input = "AliESDs.root" ; 
     analysisChain->AddFile(input);
     ag->Process(analysisChain) ; 
   }