Bug fixes (Melinda)
authormasera <masera@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 23 Jun 2009 17:15:38 +0000 (17:15 +0000)
committermasera <masera@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 23 Jun 2009 17:15:38 +0000 (17:15 +0000)
ITS/AliITSQADataMakerRec.cxx
ITS/AliITSQASDDDataMakerRec.cxx

index 25ebf9c..3d2fb79 100644 (file)
@@ -161,6 +161,7 @@ void AliITSQADataMakerRec::EndOfDetectorCycle(const char * /*fgDataName*/)
 //____________________________________________________________________________ 
 void AliITSQADataMakerRec::InitRaws()
 {  
+
   // Initialization for RAW data 
        if(fSubDetector == 0 || fSubDetector == 1) {
          AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SPD InitRaws\n");
@@ -180,18 +181,18 @@ void AliITSQADataMakerRec::InitRaws()
 void AliITSQADataMakerRec::MakeRaws(AliRawReader* rawReader)
 { 
   // Fill QA for RAW   
-  return ; 
+  //return ; 
 
   if(fSubDetector == 0 || fSubDetector == 1)  {
     Int_t rv = fSPDDataMaker->MakeRaws(rawReader) ; 
     if ( rv != 0 )
-      fSDDDataMaker->SetOffset(AliQAv1::kRAWS, fRawsQAList [AliRecoParam::AConvert(fEventSpecie)]->GetEntries());
+      fSDDDataMaker->SetOffset(AliQAv1::kRAWS, fRawsQAList[AliRecoParam::AConvert(fEventSpecie)]->GetEntries());
   }
   
   if(fSubDetector == 0 || fSubDetector == 2) {
-    Int_t rv = fSPDDataMaker->MakeRaws(rawReader) ; 
+    Int_t rv = fSDDDataMaker->MakeRaws(rawReader) ; 
     if ( rv != 0 )
-      fSSDDataMaker->SetOffset(AliQAv1::kRAWS, fRawsQAList [AliRecoParam::AConvert(fEventSpecie)]->GetEntries());
+      fSSDDataMaker->SetOffset(AliQAv1::kRAWS, fRawsQAList[AliRecoParam::AConvert(fEventSpecie)]->GetEntries());
   }
 
   if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeRaws(rawReader);
@@ -200,6 +201,7 @@ void AliITSQADataMakerRec::MakeRaws(AliRawReader* rawReader)
 //____________________________________________________________________________ 
 void AliITSQADataMakerRec::InitDigits()
 {
+
   // Initialization for DIGITS
   if(fSubDetector == 0 || fSubDetector == 1) {
     AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SPD InitDigitss\n");
@@ -218,8 +220,8 @@ void AliITSQADataMakerRec::InitDigits()
 //____________________________________________________________________________ 
 void AliITSQADataMakerRec::MakeDigits(TTree * digitsTree)
 {
+
   // Fill QA for recpoints
-  return ; 
   if(fSubDetector == 0 || fSubDetector == 1) {
     Int_t rv = fSPDDataMaker->MakeDigits(digitsTree) ; 
     if ( rv != 0 )
@@ -227,7 +229,7 @@ void AliITSQADataMakerRec::MakeDigits(TTree * digitsTree)
   }
   
   if(fSubDetector == 0 || fSubDetector == 2) {
-    Int_t rv = fSPDDataMaker->MakeDigits(digitsTree) ; 
+    Int_t rv = fSDDDataMaker->MakeDigits(digitsTree) ; 
     if ( rv != 0 )
       fSSDDataMaker->SetOffset(AliQAv1::kDIGITSR, fDigitsQAList [AliRecoParam::AConvert(fEventSpecie)]->GetEntries());
   }
@@ -238,6 +240,7 @@ void AliITSQADataMakerRec::MakeDigits(TTree * digitsTree)
 //____________________________________________________________________________ 
 void AliITSQADataMakerRec::InitRecPoints()
 {
+
   // Initialization for RECPOINTS
   if(fSubDetector == 0 || fSubDetector == 1) {
     AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SPD InitRecPoints\n");
@@ -256,18 +259,18 @@ void AliITSQADataMakerRec::InitRecPoints()
 //____________________________________________________________________________ 
 void AliITSQADataMakerRec::MakeRecPoints(TTree * clustersTree)
 {
-  return ; 
   // Fill QA for recpoints
   if(fSubDetector == 0 || fSubDetector == 1) {
     Int_t rv = fSPDDataMaker->MakeRecPoints(clustersTree) ; 
     if ( rv != 0 )
-      fSDDDataMaker->SetOffset(AliQAv1::kRECPOINTS, fDigitsQAList [AliRecoParam::AConvert(fEventSpecie)]->GetEntries());
+      fSDDDataMaker->SetOffset(AliQAv1::kRECPOINTS, fRecPointsQAList[AliRecoParam::AConvert(fEventSpecie)]->GetEntries());
   }
   
   if(fSubDetector == 0 || fSubDetector == 2) {
-    Int_t rv = fSPDDataMaker->MakeRecPoints(clustersTree) ; 
+    Int_t rv = fSDDDataMaker->MakeRecPoints(clustersTree) ; 
     if ( rv != 0 )
-      fSSDDataMaker->SetOffset(AliQAv1::kRECPOINTS, fDigitsQAList [AliRecoParam::AConvert(fEventSpecie)]->GetEntries());
+      fSSDDataMaker->SetOffset(AliQAv1::kRECPOINTS, fRecPointsQAList [AliRecoParam::AConvert(fEventSpecie)]->GetEntries());
   }
   
   if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeRecPoints(clustersTree);
@@ -277,6 +280,7 @@ void AliITSQADataMakerRec::MakeRecPoints(TTree * clustersTree)
 //____________________________________________________________________________ 
 void AliITSQADataMakerRec::InitESDs()
 {
+
   // Create ESDs histograms in ESDs subdir
 
   Bool_t expertHistogram = kTRUE;
index 3667c88..0ca00e6 100644 (file)
@@ -283,9 +283,13 @@ Int_t AliITSQASDDDataMakerRec::InitRaws()
 Int_t AliITSQASDDDataMakerRec::MakeRaws(AliRawReader* rawReader)
 { 
   // Fill QA for RAW - SDD -
-  
+
   Int_t rv = 0 ; 
 
+  // Check id histograms already created for this Event Specie
+  if ( ! fAliITSQADataMakerRec->GetRawsData(fGenRawsOffset) )
+    rv = InitRaws () ;
+
   if(!fDDLModuleMap){
     AliError("SDD DDL module map not available - skipping SDD QA");
     return rv;
@@ -293,9 +297,7 @@ Int_t AliITSQASDDDataMakerRec::MakeRaws(AliRawReader* rawReader)
   if(rawReader->GetType() != 7) return rv;  // skips non physical triggers
   AliDebug(AliQAv1::GetQADebugLevel(),"entering MakeRaws\n");                 
   
-  // Check id histograms already created for this Event Specie
-  if ( ! fAliITSQADataMakerRec->GetRawsData(fGenRawsOffset) )
-    rv = InitRaws () ;
+
   
   rawReader->Reset();       
   AliITSRawStream *stream;
@@ -419,7 +421,7 @@ Int_t AliITSQASDDDataMakerRec::MakeRaws(AliRawReader* rawReader)
 Int_t AliITSQASDDDataMakerRec::InitDigits()
 { 
 
-  //printf(" ======================================================> Init digits\n " );
+
   // Initialization for DIGIT data - SDD -  
   const Bool_t expert   = kTRUE ; 
   const Bool_t image    = kTRUE ;
@@ -453,7 +455,7 @@ Int_t AliITSQASDDDataMakerRec::InitDigits()
 //____________________________________________________________________________
 Int_t AliITSQASDDDataMakerRec::MakeDigits(TTree * digits)
 { 
-  //printf(" ======================================================> make digits\n " );
+
   // Fill QA for DIGIT - SDD -
   //AliITS *fITS  = (AliITS*)gAlice->GetModule("ITS");
   //fITS->SetTreeAddress();
@@ -481,7 +483,7 @@ Int_t AliITSQASDDDataMakerRec::MakeDigits(TTree * digits)
     digits->GetEvent(nmod);
     Int_t ndigits = iITSdigits->GetEntries();
     fAliITSQADataMakerRec->GetDigitsData(fGenDigitsOffset)->Fill(nmod,ndigits);
-    //printf(" Filled:  =======================================> %s \t %i \t %i \n",fAliITSQADataMakerRec->GetDigitsData(fGenDigitsOffset)->GetName(), nmod, ndigits );
+
     for (Int_t idig=0; idig<ndigits; idig++) {
       AliITSdigit *dig=(AliITSdigit*)iITSdigits->UncheckedAt(idig);
       Int_t iz=dig->GetCoord1();  // cell number z
@@ -498,6 +500,8 @@ Int_t AliITSQASDDDataMakerRec::MakeDigits(TTree * digits)
 //____________________________________________________________________________ 
 Int_t AliITSQASDDDataMakerRec::InitRecPoints()
 {
+
+
   // Initialization for RECPOINTS - SDD -
   const Bool_t expert   = kTRUE ; 
   const Bool_t image    = kTRUE ; 
@@ -631,9 +635,6 @@ Int_t AliITSQASDDDataMakerRec::InitRecPoints()
       
     }//online
 
-  //printf("%d SDD Recs histograms booked\n",fSDDhRecPointsTask);
-
-
   AliDebug(AliQAv1::GetQADebugLevel(),Form("%d SDD Recs histograms booked\n",fSDDhRecPointsTask));
 
   return rv ; 
@@ -644,15 +645,16 @@ Int_t AliITSQASDDDataMakerRec::MakeRecPoints(TTree * clustersTree)
 {
  // Fill QA for RecPoints - SDD -
   Int_t rv = 0 ; 
+
+  // Check id histograms already created for this Event Specie
+  if ( ! fAliITSQADataMakerRec->GetRecPointsData(fGenRecPointsOffset) )
+    rv = InitRecPoints() ;
   Int_t lay, lad, det; 
   TBranch *branchRecP = clustersTree->GetBranch("ITSRecPoints");
   if (!branchRecP) {
     AliError("can't get the branch with the ITS clusters !");
     return rv ;
   }
-  // Check id histograms already created for this Event Specie
-  if ( ! fAliITSQADataMakerRec->GetRecPointsData(fGenRecPointsOffset) )
-    rv = InitRecPoints() ;
 
   static TClonesArray statRecpoints("AliITSRecPoint") ;
   TClonesArray *recpoints = &statRecpoints;
@@ -677,6 +679,7 @@ Int_t AliITSQASDDDataMakerRec::MakeRecPoints(TTree * clustersTree)
       recp->GetGlobalXYZ(cluglo);
       Float_t rad=TMath::Sqrt(cluglo[0]*cluglo[0]+cluglo[1]*cluglo[1]); 
       Float_t phi=TMath::ATan2(cluglo[1],cluglo[0]);
       fAliITSQADataMakerRec->GetRecPointsData(9 +fGenRecPointsOffset)->Fill(recp->GetDetLocalX(),recp->GetDetLocalZ());//local distribution
       fAliITSQADataMakerRec->GetRecPointsData(2 +fGenRecPointsOffset)->Fill(cluglo[0],cluglo[1]);//global distribution YX
       fAliITSQADataMakerRec->GetRecPointsData(3 +fGenRecPointsOffset)->Fill(cluglo[2],rad);//global distribution rz
@@ -685,14 +688,14 @@ Int_t AliITSQASDDDataMakerRec::MakeRecPoints(TTree * clustersTree)
          fAliITSQADataMakerRec->GetRecPointsData(14 +fGenRecPointsOffset)->Fill(cluglo[0],cluglo[1]);//global distribution YX FSE
          fAliITSQADataMakerRec->GetRecPointsData(15 +fGenRecPointsOffset)->Fill(cluglo[2],rad);//global distribution rz FSE
        }
-      if(recp->GetLayer() ==2) {
+      if(recp->GetLayer() ==3) {
        fAliITSQADataMakerRec->GetRecPointsData(0  +fGenRecPointsOffset)->Fill(recp->GetQ()) ;//total charge of layer 3
        fAliITSQADataMakerRec->GetRecPointsData(7  +fGenRecPointsOffset)->Fill(lad);//lad pattern layer 3
        fAliITSQADataMakerRec->GetRecPointsData(10 +fGenRecPointsOffset)->Fill(rad);//r distribution layer 3
        fAliITSQADataMakerRec->GetRecPointsData(12 +fGenRecPointsOffset)->Fill(phi);// phi distribution layer 3
        fAliITSQADataMakerRec->GetRecPointsData(4  +fGenRecPointsOffset)->Fill(cluglo[2],phi);// phi distribution layer 3
       }
-      else if(recp->GetLayer() ==3) {
+      else if(recp->GetLayer() ==4) {
        fAliITSQADataMakerRec->GetRecPointsData(1  +fGenRecPointsOffset)->Fill(recp->GetQ()) ;//total charge layer 4
        fAliITSQADataMakerRec->GetRecPointsData(8  +fGenRecPointsOffset)->Fill(lad);//ladpatternlayer4
        fAliITSQADataMakerRec->GetRecPointsData(11 +fGenRecPointsOffset)->Fill(rad);//r distribution