#include "AliTRDQADataMakerSim.h"
#include "AliTRDdigitsManager.h"
#include "AliTRDgeometry.h"
-#include "AliTRDdataArrayS.h"
-//#include "AliTRDrawStreamTB.h"
+#include "AliTRDarrayADC.h"
+#include "AliTRDarraySignal.h"
+//#include "AliTRDrawStream.h"
#include "AliQAChecker.h"
//____________________________________________________________________________
AliTRDQADataMakerSim::AliTRDQADataMakerSim() :
- AliQADataMakerSim(AliQA::GetDetName(AliQA::kTRD), "TRD Quality Assurance Data Maker")
+ AliQADataMakerSim(AliQAv1::GetDetName(AliQAv1::kTRD), "TRD Quality Assurance Data Maker")
{
//
// Default constructor
}
//____________________________________________________________________________
-void AliTRDQADataMakerSim::EndOfDetectorCycle(AliQA::TASKINDEX task, TObjArray * list)
+void AliTRDQADataMakerSim::EndOfDetectorCycle(AliQAv1::TASKINDEX_t task, TObjArray ** list)
{
//
// Detector specific actions at end of cycle
//AliInfo(Form("EndOfCycle", "Fitting RecPoints %d", task));
// call the checker
- AliQAChecker::Instance()->Run(AliQA::kTRD, task, list) ;
+ AliQAChecker::Instance()->Run(AliQAv1::kTRD, task, list) ;
}
hist[0] = new TH1D("qaTRD_digits_det", ";Detector Id of the digit", 540, -0.5, 539.5);
hist[1] = new TH1D("qaTRD_digits_time", ";Time bin", 40, -0.5, 39.5);
- hist[2] = new TH1D("qaTRD_digits_amp", ";Amplitude", 100, 0, 100.);
+ hist[2] = new TH1D("qaTRD_digits_amp", ";Amplitude", 100, -5.5, 94.5);
for(Int_t i=0; i<kNhist; i++) {
hist[i]->Sumw2();
tmp->Delete();
delete tmp;
MakeHits(hits);
+ hits->Delete();
+ delete hits;
}
TIter next(digits) ;
AliTRDdigit * digit ;
+
+ // Info("Make digits", "From the arrya");
+
while ( (digit = dynamic_cast<AliTRDdigit *>(next())) ) {
+ if (digit->GetAmp() < 1) continue;
GetDigitsData(0)->Fill(digit->GetDetector());
GetDigitsData(1)->Fill(digit->GetTime());
GetDigitsData(2)->Fill(digit->GetAmp());
// Makes data from digits tree
//
+ // Info("Make digits", "From a tree");
+
AliTRDdigitsManager *digitsManager = new AliTRDdigitsManager();
digitsManager->CreateArrays();
digitsManager->ReadDigits(digits);
- for (Int_t i = 0; i < AliTRDgeometry::kNdet; i++) {
-
- AliTRDdataArrayS *digitsIn = (AliTRDdataArrayS *) digitsManager->GetDigits(i);
-
- // This is to take care of switched off super modules
- if (digitsIn->GetNtime() == 0) continue;
-
- digitsIn->Expand();
-
- //AliTRDSignalIndex* indexes = digitsManager->GetIndexes(i);
- //if (indexes->IsAllocated() == kFALSE) digitsManager->BuildIndexes(i);
-
- Int_t nRows = digitsIn->GetNrow();
- Int_t nCols = digitsIn->GetNcol();
- Int_t nTbins = digitsIn->GetNtime();
-
- for(Int_t row = 0; row < nRows; row++)
- for(Int_t col = 0; col < nCols; col++)
- for(Int_t time = 0; time < nTbins; time++) {
-
- Float_t signal = digitsIn->GetDataUnchecked(row,col,time);
- GetDigitsData(0)->Fill(i);
- GetDigitsData(1)->Fill(time);
- GetDigitsData(2)->Fill(signal);
- }
+ TH1D *histDet = (TH1D*)GetDigitsData(0);
+ TH1D *histTime = (TH1D*)GetDigitsData(1);
+ TH1D *histSignal = (TH1D*)GetDigitsData(2);
+
+ for (Int_t i = 0; i < AliTRDgeometry::kNdet; i++)
+ {
+ AliTRDarrayADC *digitsIn = (AliTRDarrayADC *) digitsManager->GetDigits(i);
+
+ // This is to take care of switched off super modules
+ if (digitsIn->GetNtime() == 0) continue;
+
+ digitsIn->Expand();
+
+ //AliTRDSignalIndex* indexes = digitsManager->GetIndexes(i);
+ //if (indexes->IsAllocated() == kFALSE) digitsManager->BuildIndexes(i);
+
+ Int_t nRows = digitsIn->GetNrow();
+ Int_t nCols = digitsIn->GetNcol();
+ Int_t nTbins = digitsIn->GetNtime();
+
+ for(Int_t row = 0; row < nRows; row++)
+ for(Int_t col = 0; col < nCols; col++)
+ for(Int_t time = 0; time < nTbins; time++)
+ {
+ Float_t signal = digitsIn->GetData(row,col,time);
+ if (signal < 1) continue;
+ histDet->Fill(i);
+ histTime->Fill(time);
+ histSignal->Fill(signal);
+ }
//delete digitsIn;
}
-
delete digitsManager;
-
}
//____________________________________________________________________________
//
AliTRDdigitsManager *digitsManager = new AliTRDdigitsManager();
+ digitsManager->SetSDigits(kTRUE);
digitsManager->CreateArrays();
digitsManager->ReadDigits(digits);
- for (Int_t i = 0; i < AliTRDgeometry::kNdet; i++) {
-
- AliTRDdataArrayS *digitsIn = (AliTRDdataArrayS *) digitsManager->GetDigits(i);
+ TH1D *histDet = (TH1D*)GetSDigitsData(0);
+ TH1D *histTime = (TH1D*)GetSDigitsData(1);
+ TH1D *histSignal = (TH1D*)GetSDigitsData(2);
+ for (Int_t i = 0; i < AliTRDgeometry::kNdet; i++)
+ {
+ AliTRDarraySignal *digitsIn = (AliTRDarraySignal *) digitsManager->GetSDigits(i);
+
// This is to take care of switched off super modules
- if (digitsIn->GetNtime() == 0) continue;
-
- digitsIn->Expand();
-
- //AliTRDSignalIndex* indexes = digitsManager->GetIndexes(i);
- //if (indexes->IsAllocated() == kFALSE) digitsManager->BuildIndexes(i);
-
- Int_t nRows = digitsIn->GetNrow();
- Int_t nCols = digitsIn->GetNcol();
- Int_t nTbins = digitsIn->GetNtime();
-
- for(Int_t row = 0; row < nRows; row++)
- for(Int_t col = 0; col < nCols; col++)
- for(Int_t time = 0; time < nTbins; time++) {
-
- Float_t signal = digitsIn->GetDataUnchecked(row,col,time);
- if (signal <= 0) continue;
- GetSDigitsData(0)->Fill(i);
- GetSDigitsData(1)->Fill(time);
- GetSDigitsData(2)->Fill(signal);
- }
-
- // delete digitsIn;
- }
-
+ if (digitsIn->GetNtime() == 0) continue;
+
+ digitsIn->Expand();
+
+ //AliTRDSignalIndex* indexes = digitsManager->GetIndexes(i);
+ //if (indexes->IsAllocated() == kFALSE) digitsManager->BuildIndexes(i);
+ Int_t nRows = digitsIn->GetNrow();
+ Int_t nCols = digitsIn->GetNcol();
+ Int_t nTbins = digitsIn->GetNtime();
+
+ for(Int_t row = 0; row < nRows; row++)
+ {
+ for(Int_t col = 0; col < nCols; col++)
+ {
+ for(Int_t time = 0; time < nTbins; time++)
+ {
+ Float_t signal = digitsIn->GetData(row,col,time);
+ if (signal < 1) continue;
+ histDet->Fill(i);
+ histTime->Fill(time);
+ histSignal->Fill(signal);
+ }
+ }
+ }
+ // delete digitsIn;
+ }
delete digitsManager;
-
}
//____________________________________________________________________________