]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - ITS/AliITSQASSDDataMakerRec.cxx
Changed error estimate in AliBWTools::GetMeanDataAndExtrapolation + removed debug...
[u/mrichter/AliRoot.git] / ITS / AliITSQASSDDataMakerRec.cxx
index 1a1987a2debcb71d9d847bf0f4041b32634772ad..83d466b52b3baa0589bef8f52c3a42d1e1ab0a85 100644 (file)
@@ -42,6 +42,7 @@
 #include "AliITSgeomTGeo.h"
 #include "AliRawEventHeaderBase.h"
 #include "AliITSRecPoint.h"
+#include "AliITSRecPointContainer.h"
 #include "AliITSdigitSSD.h"
 #include "AliITSBadChannelsSSDv2.h"
 
@@ -78,12 +79,14 @@ fCDBManager(0) {
        }
   if(fkOnline) {
     fCDBManager = AliCDBManager::Instance();
-    fCDBManager->SetDefaultStorage("local://$ALICE_ROOT/OCDB");
-    //fCDBManager->SetDefaultStorage(gSystem->Getenv("AMORE_CDB_URI"));
-    Int_t runNumber = atoi(gSystem->Getenv("DATE_RUN_NUMBER"));
-    if(!runNumber) 
-      AliWarning("DATE_RUN_NUMBER not defined!!!\n");
-    
+    //fCDBManager->SetDefaultStorage("local://$ALICE_ROOT/OCDB");
+    fCDBManager->SetDefaultStorage(gSystem->Getenv("AMORE_CDB_URI"));
+
+    Int_t runNumber = 1;
+    if(!gSystem->Getenv("DATE_RUN_NUMBER")) {
+      AliWarning("DATE_RUN_NUMBER not defined!!!\n"); }
+    else {
+      runNumber = atoi(gSystem->Getenv("DATE_RUN_NUMBER"));}
     fCDBManager->SetRun(runNumber);
     AliCDBEntry *geomGRP = fCDBManager->Get("GRP/Geometry/Data");
     if(!geomGRP) AliWarning("GRP geometry not found!!!\n");
@@ -185,8 +188,10 @@ AliITSQASSDDataMakerRec::~AliITSQASSDDataMakerRec() {
 void AliITSQASSDDataMakerRec::StartOfDetectorCycle()
 {
 
- if ( fAliITSQADataMakerRec->GetRawsData(0) == NULL ) // Raws not defined
-       return ;
+  if(fAliITSQADataMakerRec->ListExists(AliQAv1::kRAWS)==kFALSE)return;
+
+  //if ( fAliITSQADataMakerRec->GetRawsData(0) == NULL ) // Raws not defined
+  //return ;
  
  //for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
  //if (!AliQAv1::Instance()->IsEventSpecieSet(specie)) continue;
@@ -419,37 +424,36 @@ Int_t AliITSQASSDDataMakerRec::InitRaws() {
   TH1D *fHistSSDEventType = new TH1D("fHistSSDEventType",
                                     "SSD Event Type;Event type;Events",
                                     31,-1,30);
-  rv = fAliITSQADataMakerRec->Add2RawsList((new TH1D(*fHistSSDEventType)), 
+  fHistSSDEventType->SetStats(kFALSE);
+  rv = fAliITSQADataMakerRec->Add2RawsList(fHistSSDEventType, 
                                           fGenRawsOffset[specie]+fSSDRawsOffset, expert, !image, !saveCorr);
   fSSDRawsOffset += 1;
-  delete fHistSSDEventType;
   //cout<<"(AliITSQASSDDataMakerRec::InitRaws): "<<fGenRawsOffset[specie]+fSSDRawsOffset-1<<" - Name: "<<(fAliITSQADataMakerRec->GetRawsData(fGenRawsOffset[specie]))->GetName()<<endl;
   TH1D *fHistSSDDataSize = new TH1D("fHistSSDDataSize",
                                    "SSD Data Size;(SSD data size) [KB];Events",
                                    1000,0,500);
-  rv = fAliITSQADataMakerRec->Add2RawsList((new TH1D(*fHistSSDDataSize)), 
-                                          fGenRawsOffset[specie]+fSSDRawsOffset, expert, !image, !saveCorr);
-  delete fHistSSDDataSize;
+  rv = fAliITSQADataMakerRec->Add2RawsList(fHistSSDDataSize, 
+                                          fGenRawsOffset[specie]+fSSDRawsOffset, !expert, !image, !saveCorr);
   fSSDRawsOffset += 1;
   TH1D *fHistSSDDataSizePercentage = new TH1D("fHistSSDDataSizePercentage",
                                              "SSD Data Size Percentage;SSD data size [%];Events",
                                              1000,0,100);
-  rv = fAliITSQADataMakerRec->Add2RawsList((new TH1D(*fHistSSDDataSizePercentage))
+  rv = fAliITSQADataMakerRec->Add2RawsList(fHistSSDDataSizePercentage
                                           fGenRawsOffset[specie]+fSSDRawsOffset, expert, !image, !saveCorr);
-  delete fHistSSDDataSizePercentage;
   fSSDRawsOffset += 1;
   TH1D *fHistSSDDDLId = new TH1D("fHistSSDDDLId",
                                 "SSD DDL Id;DDL id;Events",16,511.5,527.5);
-  rv = fAliITSQADataMakerRec->Add2RawsList((new TH1D(*fHistSSDDDLId)), 
+  fHistSSDDDLId->SetStats(kFALSE);
+  rv = fAliITSQADataMakerRec->Add2RawsList(fHistSSDDDLId, 
                                           fGenRawsOffset[specie]+fSSDRawsOffset, expert, !image, !saveCorr);
-  delete fHistSSDDDLId;
   fSSDRawsOffset += 1;
   TH1D *fHistSSDDataSizePerDDL = new TH1D("fHistSSDDataSizePerDDL",
                                          "SSD Data Size Per DDL;DDL id;<SSD data size> [KB]",
                                          16,511.5,527.5);
-  rv = fAliITSQADataMakerRec->Add2RawsList((new TH1D(*fHistSSDDataSizePerDDL)), 
+  fHistSSDDataSizePerDDL->SetStats(kFALSE);
+  rv = fAliITSQADataMakerRec->Add2RawsList(fHistSSDDataSizePerDDL, 
                                           fGenRawsOffset[specie]+fSSDRawsOffset, !expert, image, !saveCorr);
-  delete fHistSSDDataSizePerDDL;
+
   fSSDRawsOffset += 1;
   TH1D *fHistSSDDataSizeDDL[fgkNumOfDDLs];
   for(Int_t i = 1; i < fgkNumOfDDLs+1; i++) {
@@ -458,24 +462,22 @@ Int_t AliITSQASSDDataMakerRec::InitRaws() {
     fHistSSDDataSizeDDL[i-1] = new TH1D(gName.Data(),
                                        Form("%s;(SSD data size) [KB];Events", gTitle.Data()),
                                        100,0,50);
-    rv = fAliITSQADataMakerRec->Add2RawsList((new TH1D(*fHistSSDDataSizeDDL[i-1]))
+    rv = fAliITSQADataMakerRec->Add2RawsList(fHistSSDDataSizeDDL[i-1]
                                             fGenRawsOffset[specie]+fSSDRawsOffset, expert, !image, !saveCorr);
     fSSDRawsOffset += 1;
   }
-  for(Int_t i = 1; i < fgkNumOfDDLs+1; i++) delete fHistSSDDataSizeDDL[i-1];
   
   TH1D *fHistSSDLDCId = new TH1D("fHistSSDLDCId","SSD LDC Id;LDC id;Events",8,169.5,177.5);
-  rv = fAliITSQADataMakerRec->Add2RawsList((new TH1D(*fHistSSDLDCId)), 
+  fHistSSDLDCId->SetStats(kFALSE);
+  rv = fAliITSQADataMakerRec->Add2RawsList(fHistSSDLDCId, 
                                           fGenRawsOffset[specie]+fSSDRawsOffset, expert, !image, !saveCorr);
-  delete fHistSSDLDCId;
   fSSDRawsOffset += 1;
   TH1D *fHistSSDDataSizePerLDC = new TH1D("fHistSSDDataSizePerLDC",
                                          "SSD Data Size Per LDC;LDC id;<SSD data size> [KB]",
                                          8,169.5,177.5);
-  rv = fAliITSQADataMakerRec->Add2RawsList((new TH1D(*fHistSSDDataSizePerLDC)), 
-                                          fGenRawsOffset[specie]+fSSDRawsOffset, !expert, image, !saveCorr);
-  delete fHistSSDDataSizePerLDC;
-  fSSDRawsOffset += 1;
+  fHistSSDDataSizePerLDC->SetStats(kFALSE);
+  rv = fAliITSQADataMakerRec->Add2RawsList(fHistSSDDataSizePerLDC, 
+                                          fGenRawsOffset[specie]+fSSDRawsOffset, expert, image, !saveCorr);  fSSDRawsOffset += 1;
   TH1D *fHistSSDDataSizeLDC[fgkNumOfLDCs];
   for(Int_t i = 1; i < fgkNumOfLDCs+1; i++) {
     gName = "fHistSSDDataSizeLDC"; 
@@ -492,11 +494,10 @@ Int_t AliITSQASSDDataMakerRec::InitRaws() {
     fHistSSDDataSizeLDC[i-1] = new TH1D(gName.Data(),
                                        Form("%s;SSD data size [KB];Events", gTitle.Data()),
                                        1000,0,100);
-    rv = fAliITSQADataMakerRec->Add2RawsList((new TH1D(*fHistSSDDataSizeLDC[i-1]))
+    rv = fAliITSQADataMakerRec->Add2RawsList(fHistSSDDataSizeLDC[i-1]
                                             fGenRawsOffset[specie]+fSSDRawsOffset, expert, !image, !saveCorr);
     fSSDRawsOffset += 1;
   }
-  for(Int_t i = 1; i < fgkNumOfLDCs+1; i++) delete fHistSSDDataSizeLDC[i-1];
   fSSDRawsCommonLevelOffset = fSSDRawsOffset;
   
   if(fkOnline) {
@@ -522,7 +523,7 @@ Int_t AliITSQASSDDataMakerRec::InitRaws() {
       gName += "_Module"; gName += gModule; 
       gTitle += "_Module"; gTitle += gModule; 
       
-      fHistSSDOccupancyModule[i-500] = new TH1D(gName.Data(),Form("%s;N_{strip};Occupancy [%]", gTitle.Data()),
+      fHistSSDOccupancyModule[i-500] = new TH1D(gName.Data(),Form("%s;N_{strip};Occupancy [%%]", gTitle.Data()),
                                                2*fgkNumberOfPSideStrips,0,2*fgkNumberOfPSideStrips);
       fHistSSDOccupancyModule[i-500]->GetXaxis()->SetTitleColor(1);
       rv = fAliITSQADataMakerRec->Add2RawsList(fHistSSDOccupancyModule[i-500], 
@@ -553,7 +554,7 @@ Int_t AliITSQASSDDataMakerRec::InitRaws() {
         gName += "_PSide";
         gTitle += "_PSide";
         fHistSSDOccupancyLadder[occupancyCounter] = new TH1D(gName.Data(),
-                                                             Form("%s;Module number;Occupancy [%]", gTitle.Data()),
+                                                             Form("%s;Module number;Occupancy [%%]", gTitle.Data()),
                                                              AliITSgeomTGeo::GetNDetectors(iLayer),
                                                              0.5,AliITSgeomTGeo::GetNDetectors(iLayer)+0.5);
         fHistSSDOccupancyLadder[occupancyCounter]->GetXaxis()->SetTitleColor(1);
@@ -578,7 +579,7 @@ Int_t AliITSQASSDDataMakerRec::InitRaws() {
         gName += "_NSide";
         gTitle += "_NSide";
         fHistSSDOccupancyLadder[occupancyCounter] = new TH1D(gName.Data(),
-                                                             Form("%s;Module number;Occupancy [%]", gTitle.Data()),
+                                                             Form("%s;Module number;Occupancy [%%]", gTitle.Data()),
                                                              AliITSgeomTGeo::GetNDetectors(iLayer),
                                                              0.5,AliITSgeomTGeo::GetNDetectors(iLayer)+0.5);
         fHistSSDOccupancyLadder[occupancyCounter]->GetXaxis()->SetTitleColor(1);
@@ -602,6 +603,7 @@ Int_t AliITSQASSDDataMakerRec::InitRaws() {
       sprintf(fLabel,"%d",iBin);
       fHistSSDOccupancyLayer5->GetXaxis()->SetBinLabel(iBin,fLabel);
     }
+    fHistSSDOccupancyLayer5->SetStats(kFALSE);
     rv = fAliITSQADataMakerRec->Add2RawsList(fHistSSDOccupancyLayer5, 
                                        fGenRawsOffset[specie]+fSSDRawsOffset, expert, !image, !saveCorr);
     fSSDRawsOffset += 1;
@@ -616,6 +618,7 @@ Int_t AliITSQASSDDataMakerRec::InitRaws() {
       sprintf(fLabel,"%d",iBin);
       fHistSSDOccupancyLayer6->GetXaxis()->SetBinLabel(iBin,fLabel);
     }
+    fHistSSDOccupancyLayer6->SetStats(kFALSE);
     rv = fAliITSQADataMakerRec->Add2RawsList(fHistSSDOccupancyLayer6, 
                                         fGenRawsOffset[specie]+fSSDRawsOffset, expert, !image, !saveCorr);
     fSSDRawsOffset += 1;
@@ -632,6 +635,7 @@ Int_t AliITSQASSDDataMakerRec::InitRaws() {
       sprintf(fLabel,"%d",iBin);
       fHistSSDOccupancyThresholdLayer5->GetXaxis()->SetBinLabel(iBin,fLabel);
     }
+    fHistSSDOccupancyThresholdLayer5->SetStats(kFALSE);
     rv = fAliITSQADataMakerRec->Add2RawsList(fHistSSDOccupancyThresholdLayer5, 
                                         fGenRawsOffset[specie]+fSSDRawsOffset, expert, !image, !saveCorr);
     fSSDRawsOffset += 1;
@@ -646,6 +650,7 @@ Int_t AliITSQASSDDataMakerRec::InitRaws() {
       sprintf(fLabel,"%d",iBin);
       fHistSSDOccupancyThresholdLayer6->GetXaxis()->SetBinLabel(iBin,fLabel);
     }
+    fHistSSDOccupancyThresholdLayer6->SetStats(kFALSE);
     rv = fAliITSQADataMakerRec->Add2RawsList(fHistSSDOccupancyThresholdLayer6, 
                                        fGenRawsOffset[specie]+fSSDRawsOffset, expert, !image, !saveCorr);
     fSSDRawsOffset += 1;
@@ -662,6 +667,7 @@ Int_t AliITSQASSDDataMakerRec::InitRaws() {
       sprintf(fLabel,"%d",iBin);
       fHistSSDAverageOccupancyLayer5->GetXaxis()->SetBinLabel(iBin,fLabel);
     }
+    fHistSSDAverageOccupancyLayer5->SetStats(kFALSE);
     rv = fAliITSQADataMakerRec->Add2RawsList(fHistSSDAverageOccupancyLayer5, 
                                             fGenRawsOffset[specie]+fSSDRawsOffset, !expert, image, !saveCorr);
     fSSDRawsOffset += 1;
@@ -676,6 +682,7 @@ Int_t AliITSQASSDDataMakerRec::InitRaws() {
       sprintf(fLabel,"%d",iBin);
       fHistSSDAverageOccupancyLayer6->GetXaxis()->SetBinLabel(iBin,fLabel);
     }
+    fHistSSDAverageOccupancyLayer6->SetStats(kFALSE);
     rv = fAliITSQADataMakerRec->Add2RawsList(fHistSSDAverageOccupancyLayer6, 
                                             fGenRawsOffset[specie]+fSSDRawsOffset, !expert, image, !saveCorr);
     fSSDRawsOffset += 1;
@@ -1534,23 +1541,23 @@ Int_t AliITSQASSDDataMakerRec::MakeRecPoints(TTree *clustersTree)
   Int_t rv = 0 ; 
   Int_t gLayer = 0, gLadder = 0, gModule = 0;
   Int_t lLadderLocationY = 0;
-  TBranch *branchRecP = clustersTree->GetBranch("ITSRecPoints");
-  if (!branchRecP) { 
-    AliError("can't get the branch with the ITS clusters !");
+  AliITSRecPointContainer* rpcont=AliITSRecPointContainer::Instance();
+  TClonesArray *recpoints = rpcont->FetchClusters(0,clustersTree); 
+  if(!rpcont->GetStatusOK() || !rpcont->IsSSDActive()){
+    AliError("can't get SSD clusters !");
     return rv;
   }
  
   //AliInfo(Form("fAliITSQADataMakerRec->GetEventSpecie() %d\n",fAliITSQADataMakerRec->GetEventSpecie()));
   //AliInfo(Form("fGenRecPointsOffset[fAliITSQADataMakerRec->GetEventSpecie()] %d\n",fGenRecPointsOffset[fAliITSQADataMakerRec->GetEventSpecie()]));
-  static TClonesArray statRecpoints("AliITSRecPoint");
-  TClonesArray *recpoints = &statRecpoints;
-  branchRecP->SetAddress(&recpoints);
   Int_t nClustersLayer5 = 0, nClustersLayer6 = 0;
   Int_t npoints = 0;      
   Float_t cluglo[3]={0.,0.,0.}; 
   //printf("*-*-*-*-*-*-*---*-*-*-------*-*-*-*-*-*-***************AliITSQASSDataMakerRec::MakeRecpoints STEP1 \n");
-  for(Int_t module = 0; module < clustersTree->GetEntries(); module++){
-    branchRecP->GetEvent(module);
+  Int_t firMod = AliITSgeomTGeo::GetModuleIndex(5,1,1);
+  Int_t lasMod =  AliITSgeomTGeo::GetNModules();
+  for(Int_t module = firMod; module < lasMod; module++){
+    recpoints = rpcont->UncheckedGetClusters(module);
     npoints += recpoints->GetEntries();
     AliITSgeomTGeo::GetModuleId(module,gLayer,gLadder,gModule);
     //printf("SSDDataMAkerRec:::::::::::::::::::::::gLayer ========== %d \n\n",gLayer);
@@ -1624,22 +1631,21 @@ Int_t AliITSQASSDDataMakerRec::MakeRecPoints(TTree *clustersTree)
   fAliITSQADataMakerRec->GetRecPointsData(fGenRecPointsOffset[fAliITSQADataMakerRec->GetEventSpecie()] + 2)->Fill(nClustersLayer5);
   fAliITSQADataMakerRec->GetRecPointsData(fGenRecPointsOffset[fAliITSQADataMakerRec->GetEventSpecie()] + 3)->Fill(nClustersLayer6);
 
-  statRecpoints.Clear();
   return rv ; 
 }
 
 //____________________________________________________________________________ 
-Int_t AliITSQASSDDataMakerRec::GetOffset(AliQAv1::TASKINDEX_t task) {
+Int_t AliITSQASSDDataMakerRec::GetOffset(AliQAv1::TASKINDEX_t task,Int_t specie) {
   // Returns offset number according to the specified task 
   Int_t offset=0;
   if( task == AliQAv1::kRAWS ) {
-    offset=fGenRawsOffset[fAliITSQADataMakerRec->GetEventSpecie()];  
+    offset=fGenRawsOffset[specie];  
   }
   else if( task == AliQAv1::kDIGITSR ) {
-    offset=fGenDigitsOffset[fAliITSQADataMakerRec->GetEventSpecie()];   
+    offset=fGenDigitsOffset[specie];   
   }
   else if( task == AliQAv1::kRECPOINTS ) {
-    offset=fGenRecPointsOffset[fAliITSQADataMakerRec->GetEventSpecie()];   
+    offset=fGenRecPointsOffset[specie];   
   }
 
   return offset;
@@ -1680,3 +1686,11 @@ Int_t AliITSQASSDDataMakerRec::GetTaskHisto(AliQAv1::TASKINDEX_t task) {
 
   return histotot;
 }
+
+//____________________________________________________________________________ 
+void AliITSQASSDDataMakerRec::ResetDetector(AliQAv1::TASKINDEX_t task)
+{
+  if(task==AliQAv1::kRAWS)ResetRawsMonitoredObjects();
+
+
+}