Added debugging printings. Fix in computation of statistical values.

[u/mrichter/AliRoot.git] / STEER / AliGlobalQADataMaker.cxx

diff --git a/STEER/AliGlobalQADataMaker.cxx b/STEER/AliGlobalQADataMaker.cxx
index 59f1ba43478b799e9f2f9b046efbdc8ab4075032..dfe17c3f68f300887f1e136f9e99c91fa7a72a9d 100644 (file)
--- a/STEER/AliGlobalQADataMaker.cxx
+++ b/STEER/AliGlobalQADataMaker.cxx
@@ -8,13 +8,30 @@
 #include <TPDGCode.h>
 #include <TH1F.h>
 
+#include "AliQAChecker.h"
 #include "AliGlobalQADataMaker.h"
 #include "AliGeomManager.h"
 #include "AliESDEvent.h"
 #include "AliESDv0.h"
+#include "AliRawReader.h"
 
 ClassImp(AliGlobalQADataMaker)
  
+//____________________________________________________________________________ 
+void AliGlobalQADataMaker::EndOfDetectorCycle(AliQA::TASKINDEX_t task, TObjArray * list)
+{
+  //Detector specific actions at end of cycle
+  // do the QA checking
+  AliQAChecker::Instance()->Run(AliQA::kGLOBAL, task, list) ;  
+}
+
+//____________________________________________________________________________ 
+void AliGlobalQADataMaker::InitRaws()
+{
+  // create Raws histograms in Raws subdir
+}
+
+//____________________________________________________________________________ 
 void AliGlobalQADataMaker::InitRecPoints() {
   //------------------------------------------------------
   // This function books the histograms of *track*residuals*
@@ -57,8 +74,37 @@ void AliGlobalQADataMaker::InitRecPoints() {
     h=new TH1F(name[i+1],name[i+1],100,-5.,5.);
     Add2RecPointsList(h,i+1);    
   }
+
+  Add2RecPointsList(
+  new TH1F("SSD1 absolute residuals Y for Z<0 (cm)",
+           "SSD1 absolute residuals Y for Z<0 (cm)",100,-2.,2.),40);
+  Add2RecPointsList(
+  new TH1F("SSD1 absolute residuals Z for Z<0 (cm)",
+           "SSD1 absolute residuals Z for Z<0 (cm)",100,-2.,2.),41);
+  Add2RecPointsList(
+  new TH1F("SSD1 absolute residuals Y for Z>0 (cm)",
+           "SSD1 absolute residuals Y for Z>0 (cm)",100,-2.,2.),42);
+  Add2RecPointsList(
+  new TH1F("SSD1 absolute residuals Z for Z>0 (cm)",
+           "SSD1 absolute residuals Z for Z>0 (cm)",100,-2.,2.),43);
+
+
+  Add2RecPointsList(
+  new TH1F("SSD2 absolute residuals Y for Z<0 (cm)",
+           "SSD2 absolute residuals Y for Z<0 (cm)",100,-3.,3.),44);
+  Add2RecPointsList(
+  new TH1F("SSD2 absolute residuals Z for Z<0 (cm)",
+           "SSD2 absolute residuals Z for Z<0 (cm)",100,-3.,3.),45);
+  Add2RecPointsList(
+  new TH1F("SSD2 absolute residuals Y for Z>0 (cm)",
+           "SSD2 absolute residuals Y for Z>0 (cm)",100,-3.,3.),46);
+  Add2RecPointsList(
+  new TH1F("SSD2 absolute residuals Z for Z>0 (cm)",
+           "SSD2 absolute residuals Z for Z>0 (cm)",100,-3.,3.),47);
+  
 }
 
+//____________________________________________________________________________ 
 void AliGlobalQADataMaker::InitESDs() {
   //------------------------------------------------------
   // This function books the ESD QA histograms
@@ -117,6 +163,15 @@ void AliGlobalQADataMaker::InitESDs() {
 
 }
 
+//____________________________________________________________________________
+void AliGlobalQADataMaker::MakeRaws(AliRawReader* rawReader)
+{
+  //Fill prepared histograms with Raw digit properties
+  rawReader->Reset() ;
+
+}
+
+//____________________________________________________________________________ 
 void AliGlobalQADataMaker::MakeESDs(AliESDEvent * event) {
   //-----------------------------------------------------------
   // This function fills the ESD QA histograms
@@ -137,7 +192,7 @@ void AliGlobalQADataMaker::MakeESDs(AliESDEvent * event) {
     if (track->IsOn(AliESDtrack::kTPCrefit)) {
       Int_t n =track->GetTPCNcls();
       Int_t nf=track->GetTPCNclsF();      // number of crossed TPC pad rows
-      GetESDsData(kClr1)->Fill(Float_t(n)/nf);
+      if (nf>0) GetESDsData(kClr1)->Fill(Float_t(n)/nf);
     }
 
     if (track->IsOn(AliESDtrack::kTRDrefit)) {