#include "AliZDCQADataMakerSim.h"\r
#include "AliZDCHit.h"\r
#include "AliZDCDigit.h"\r
-#include "AliZDCRawStream.h"\r
\r
ClassImp(AliZDCQADataMakerSim)\r
\r
//____________________________________________________________________________ \r
AliZDCQADataMakerSim::AliZDCQADataMakerSim() : \r
- AliQADataMakerSim(AliQAv1::GetDetName(AliQAv1::kZDC), "ZDC Quality Assurance Data Maker"),\r
- fDigit(0)\r
+ AliQADataMakerSim(AliQAv1::GetDetName(AliQAv1::kZDC), "ZDC Quality Assurance Data Maker")\r
{\r
// ctor\r
}\r
\r
//____________________________________________________________________________ \r
AliZDCQADataMakerSim::AliZDCQADataMakerSim(const AliZDCQADataMakerSim& qadm) :\r
- AliQADataMakerSim(), \r
- fDigit(0) \r
+ AliQADataMakerSim()\r
{\r
//copy ctor \r
SetName((const char*)qadm.GetName()); \r
TH1F * hDigZNATotlg = new TH1F("hDigZNATotlg", "Digit lg signal in ZNA", 100, 0., 6000.);\r
TH1F * hDigZPCTotlg = new TH1F("hDigZPCTotlg", "Digit lg signal in ZPC", 100, 0., 6000.);\r
TH1F * hDigZPATotlg = new TH1F("hDigZPATotlg", "Digit lg signal in ZPA", 100, 0., 6000.);\r
- Add2DigitsList(hDigZNCTotlg, 12, !expert, image);\r
- Add2DigitsList(hDigZNATotlg, 13, !expert, image);\r
- Add2DigitsList(hDigZPCTotlg, 14, !expert, image);\r
- Add2DigitsList(hDigZPATotlg, 15, !expert, image);\r
+ Add2DigitsList(hDigZNCTotlg, 12, expert, !image);\r
+ Add2DigitsList(hDigZNATotlg, 13, expert, !image);\r
+ Add2DigitsList(hDigZPCTotlg, 14, expert, !image);\r
+ Add2DigitsList(hDigZPATotlg, 15, expert, !image);\r
//\r
TH1F * hDigSumQZNClg = new TH1F("hDigSumQZNClg", "Signal in 4 ZNC PMQlg",100, 0., 4000.);\r
TH1F * hDigSumQZNAlg = new TH1F("hDigSumQZNAlg", "Signal in 4 ZNA PMQlg",100, 0., 4000.);\r
void AliZDCQADataMakerSim::MakeHits()\r
{\r
//filling QA histos for Hits\r
- //\r
\r
// Check id histograms already created for this Event Specie\r
- if ( ! GetHitsData(0) )\r
- InitHits() ;\r
+ if( !GetHitsData(0) ) InitHits();\r
\r
TIter next(fHitsArray); \r
AliZDCHit * hit; \r
void AliZDCQADataMakerSim::MakeHits(TTree * hitTree)\r
{\r
// make QA data from Hit Tree\r
- //\r
if(!hitTree){\r
- AliError("Hit Tree not found!"); \r
- return;\r
- }\r
- //\r
+ AliError("Can't get ZDC hit tree!!");\r
+ return; \r
+ } \r
\r
TBranch * branch = hitTree->GetBranch("ZDC") ;\r
\r
AliError("ZDC branch in Hit Tree not found!"); \r
return;\r
} \r
- else{\r
- if (fHitsArray) \r
- fHitsArray->Clear() ; \r
- char** add = (char**) (branch->GetAddress());\r
- if(add){\r
- fHitsArray = (TClonesArray*)(*add);\r
- } \r
- else{\r
- if(!fHitsArray) fHitsArray = new TClonesArray("AliZDCHit", 1000);\r
- branch->SetAddress(&fHitsArray);\r
- }\r
- Int_t ntracks = (Int_t) hitTree->GetEntries();\r
- if (ntracks<=0) return;\r
- //\r
- for(Int_t itrack=0; itrack<ntracks; itrack++){\r
- \r
- branch->GetEntry(itrack);\r
- //\r
- MakeHits(); \r
- fHitsArray->Clear();\r
- } \r
- }\r
+ \r
+ if(fHitsArray) fHitsArray->Clear() ; \r
+ else fHitsArray = new TClonesArray("AliZDCHit", 1000);\r
+ \r
+ branch->SetAddress(&fHitsArray) ;\r
+ for (Int_t ientry = 0 ; ientry < branch->GetEntries() ; ientry++) {\r
+ branch->GetEntry(ientry) ;\r
+ MakeHits() ; \r
+ fHitsArray->Clear() ; \r
+ } \r
+\r
}\r
\r
//___________________________________________________________________________\r
-void AliZDCQADataMakerSim::MakeDigits(TTree *digitTree )\r
+void AliZDCQADataMakerSim::MakeDigits(TTree *digitTree)\r
{\r
// makes data from Digit Tree\r
+ if( !GetDigitsData(0) ) InitDigits();\r
+\r
+ if(!digitTree){\r
+ AliError("Can't get ZDC digit tree!!");\r
+ return; \r
+ } \r
+ \r
TBranch * branch = digitTree->GetBranch("ZDC");\r
if(!branch){\r
- AliError("ZDC branch in Digit Tree not found"); \r
+ AliError("ZDC branch in digit tree not found"); \r
return;\r
} \r
- \r
- // Check id histograms already created for this Event Specie\r
- if ( ! GetDigitsData(0) )\r
- InitDigits() ;\r
- \r
- branch->SetAddress(&fDigit);\r
- \r
- Int_t ndig = digitTree->GetEntries();\r
- \r
+ \r
+ AliZDCDigit *digit = 0x0;\r
+ branch->SetAddress(&digit);\r
+ \r
Float_t adcSum_ZNC=0., adcSum_ZNA=0., adcSum_ZPC=0., adcSum_ZPA=0.;\r
Float_t adcSumQ_ZNC=0., adcSumQ_ZNA=0., adcSumQ_ZPC=0., adcSumQ_ZPA=0.;\r
Float_t adcSum_ZNC_lg=0., adcSum_ZNA_lg=0., adcSum_ZPC_lg=0., adcSum_ZPA_lg=0.;\r
Float_t adcSumQ_ZNC_lg=0., adcSumQ_ZNA_lg=0., adcSumQ_ZPC_lg=0., adcSumQ_ZPA_lg=0.;\r
- //\r
- for(Int_t i = 0; i < ndig; i++){\r
- digitTree->GetEntry(i);\r
- if(fDigit->GetSector(0)==1){\r
- adcSum_ZNC += fDigit->GetADCValue(0);\r
- adcSum_ZNC_lg += fDigit->GetADCValue(1);\r
+ \r
+ Int_t ndig = digitTree->GetEntries();\r
+ for(Int_t i=0; i<ndig; i++){\r
+ branch->GetEntry(i);\r
+ \r
+ if(digit->GetSector(0)==1 && digit->GetSector(1)!=5){\r
+ adcSum_ZNC += digit->GetADCValue(0);\r
+ adcSum_ZNC_lg += digit->GetADCValue(1);\r
//\r
- if(fDigit->GetSector(1)!=0){\r
- adcSumQ_ZNC += fDigit->GetADCValue(0);\r
- adcSumQ_ZNC_lg+= fDigit->GetADCValue(1);\r
+ if(digit->GetSector(1)!=0){\r
+ adcSumQ_ZNC += digit->GetADCValue(0);\r
+ adcSumQ_ZNC_lg+= digit->GetADCValue(1);\r
}\r
else{\r
- GetDigitsData(8)->Fill(fDigit->GetADCValue(0));\r
- GetDigitsData(20)->Fill(fDigit->GetADCValue(1));\r
+ GetDigitsData(8)->Fill(digit->GetADCValue(0));\r
+ GetDigitsData(20)->Fill(digit->GetADCValue(1));\r
}\r
}\r
- else if(fDigit->GetSector(0)==2){\r
- adcSum_ZPC += fDigit->GetADCValue(0);\r
- adcSum_ZPC_lg += fDigit->GetADCValue(1);\r
+ else if(digit->GetSector(0)==2){\r
+ adcSum_ZPC += digit->GetADCValue(0);\r
+ adcSum_ZPC_lg += digit->GetADCValue(1);\r
//\r
- if(fDigit->GetSector(1)!=0){\r
- adcSumQ_ZPC += fDigit->GetADCValue(0);\r
- adcSumQ_ZPC_lg+= fDigit->GetADCValue(1);\r
+ if(digit->GetSector(1)!=0){\r
+ adcSumQ_ZPC += digit->GetADCValue(0);\r
+ adcSumQ_ZPC_lg+= digit->GetADCValue(1);\r
}\r
else{\r
- GetDigitsData(10)->Fill(fDigit->GetADCValue(0));\r
- GetDigitsData(22)->Fill(fDigit->GetADCValue(1));\r
+ GetDigitsData(10)->Fill(digit->GetADCValue(0));\r
+ GetDigitsData(22)->Fill(digit->GetADCValue(1));\r
}\r
}\r
- else if(fDigit->GetSector(0)==4){\r
- adcSum_ZNA += fDigit->GetADCValue(0);\r
- adcSum_ZNA_lg += fDigit->GetADCValue(1);\r
+ else if(digit->GetSector(0)==4 && digit->GetSector(1)!=5){\r
+ adcSum_ZNA += digit->GetADCValue(0);\r
+ adcSum_ZNA_lg += digit->GetADCValue(1);\r
//\r
- if(fDigit->GetSector(1)!=0){\r
- adcSumQ_ZNA += fDigit->GetADCValue(0);\r
- adcSumQ_ZNA_lg+= fDigit->GetADCValue(1);\r
+ if(digit->GetSector(1)!=0){\r
+ adcSumQ_ZNA += digit->GetADCValue(0);\r
+ adcSumQ_ZNA_lg+= digit->GetADCValue(1);\r
}\r
else{\r
- GetDigitsData(9)->Fill(fDigit->GetADCValue(0));\r
- GetDigitsData(21)->Fill(fDigit->GetADCValue(1));\r
+ GetDigitsData(9)->Fill(digit->GetADCValue(0));\r
+ GetDigitsData(21)->Fill(digit->GetADCValue(1));\r
}\r
}\r
- else if(fDigit->GetSector(0)==5){\r
- adcSum_ZPA += fDigit->GetADCValue(0);\r
- adcSum_ZPA_lg += fDigit->GetADCValue(1);\r
+ else if(digit->GetSector(0)==5){\r
+ adcSum_ZPA += digit->GetADCValue(0);\r
+ adcSum_ZPA_lg += digit->GetADCValue(1);\r
//\r
- if(fDigit->GetSector(1)!=0){\r
- adcSumQ_ZPA += fDigit->GetADCValue(0);\r
- adcSumQ_ZPA_lg+= fDigit->GetADCValue(1);\r
+ if(digit->GetSector(1)!=0){\r
+ adcSumQ_ZPA += digit->GetADCValue(0);\r
+ adcSumQ_ZPA_lg+= digit->GetADCValue(1);\r
}\r
else{\r
- GetDigitsData(11)->Fill(fDigit->GetADCValue(0));\r
- GetDigitsData(23)->Fill(fDigit->GetADCValue(1));\r
+ GetDigitsData(11)->Fill(digit->GetADCValue(0));\r
+ GetDigitsData(23)->Fill(digit->GetADCValue(1));\r
}\r
}\r
}\r
GetDigitsData(17)->Fill(adcSumQ_ZNA_lg);\r
GetDigitsData(18)->Fill(adcSumQ_ZPC_lg);\r
GetDigitsData(19)->Fill(adcSumQ_ZPA_lg);\r
+\r
}\r
\r
//____________________________________________________________________________\r