#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
- fHits(0),\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
- fHits(0),\r
- fDigit(0) \r
+ AliQADataMakerSim()\r
{\r
//copy ctor \r
SetName((const char*)qadm.GetName()); \r
//\r
const Bool_t expert = kTRUE ; \r
const Bool_t image = kTRUE ; \r
+\r
+ TH1F * hHitsZNCTot = new TH1F("hHitsZNCTot", "Signal in ZNC; N_{phe}", 100, 0., 6000.);\r
+ TH1F * hHitsZNATot = new TH1F("hHitsZNATot", "Signal in ZNA; N_{phe}", 100, 0., 6000.);\r
+ TH1F * hHitsZPCTot = new TH1F("hHitsZPCTot", "Signal in ZPC; N_{phe}", 100, 0., 6000.);\r
+ TH1F * hHitsZPATot = new TH1F("hHitsZPATot", "Signal in ZPA; N_{phe}", 100, 0., 6000.);\r
+ Add2HitsList(hHitsZNCTot, 0, !expert, image);\r
+ Add2HitsList(hHitsZNATot, 1, !expert, image);\r
+ Add2HitsList(hHitsZPCTot, 2, !expert, image);\r
+ Add2HitsList(hHitsZPATot, 3, !expert, image);\r
+ //\r
+ TH1F * hHitsSumQZNC = new TH1F("hHitsSumQZNC", "Signal in 4 ZNC PMQ; N_{phe}",100, 0., 4000.);\r
+ TH1F * hHitsSumQZNA = new TH1F("hHitsSumQZNA", "Signal in 4 ZNA PMQ; N_{phe}",100, 0., 4000.);\r
+ TH1F * hHitsSumQZPC = new TH1F("hHitsSumQZPC", "Signal in 4 ZPC PMQ; N_{phe}",100, 0., 4000.);\r
+ TH1F * hHitsSumQZPA = new TH1F("hHitsSumQZPA", "Signal in 4 ZPA PMQ; N_{phe}",100, 0., 4000.);\r
+ Add2HitsList(hHitsSumQZNC, 4, expert, !image);\r
+ Add2HitsList(hHitsSumQZNA, 5, expert, !image);\r
+ Add2HitsList(hHitsSumQZPC, 6, expert, !image);\r
+ Add2HitsList(hHitsSumQZPA, 7, expert, !image);\r
+ //\r
+ TH1F * hHitsPMCZNC = new TH1F("hHitsPMCZNC", "Signal in ZNC PMC; N_{phe}",100, 0., 4000.);\r
+ TH1F * hHitsPMCZNA = new TH1F("hHitsPMCZNA", "Signal in ZNA PMC; N_{phe}",100, 0., 4000.);\r
+ TH1F * hHitsPMCZPC = new TH1F("hHitsPMCZPC", "Signal in ZPC PMC; N_{phe}",100, 0., 4000.);\r
+ TH1F * hHitsPMCZPA = new TH1F("hHitsPMCZPA", "Signal in ZPA PMC; N_{phe}",100, 0., 4000.);\r
+ Add2HitsList(hHitsPMCZNC, 8, expert, !image);\r
+ Add2HitsList(hHitsPMCZNA, 9, expert, !image);\r
+ Add2HitsList(hHitsPMCZPC, 10, expert, !image);\r
+ Add2HitsList(hHitsPMCZPA, 11, expert, !image);\r
\r
- TH2F * hZNCh = new TH2F("hZNCh", "Hits centroid in ZNC;Centroid position [cm];Counts", 100, -5.,5.,100,-5.,5.);\r
- TH2F * hZNAh = new TH2F("hZNAh", "Hits centroid in ZNA;Centroid position [cm];Counts", 100, -5.,5.,100,-5.,5.);\r
+ TH2F * hHitsZNCh = new TH2F("hHitsZNCh", "Hits centroid in ZNC;Centroid position [cm];Counts", 100, -5.,5.,100,-5.,5.);\r
+ TH2F * hHitsZNAh = new TH2F("hHitsZNAh", "Hits centroid in ZNA;Centroid position [cm];Counts", 100, -5.,5.,100,-5.,5.);\r
+ Add2HitsList(hHitsZNCh, 12, !expert, image);\r
+ Add2HitsList(hHitsZNAh, 13, !expert, image);\r
// NB -> For the moment no check is performesd on ZP centroids\r
-// TH2F * hZPCh = new TH2F("hZPCh", "Hits centroid in ZPC", 100,-12.,12.,100,-12.,12.); \r
-// TH2F * hZPAh = new TH2F("hZPAh", "Hits centroid in ZPA", 100,-12.,12.,100,-12.,12.); \r
- Add2HitsList(hZNCh, 0, !expert, image);\r
- Add2HitsList(hZNAh, 1, !expert, image);\r
-// Add2HitsList(hZPCh, 2, !expert, image);\r
-// Add2HitsList(hZPAh, 3, !expert, image);\r
+ TH2F * hHitsZPCh = new TH2F("hHitsZPCh", "Hits centroid in ZPC;Centroid position [cm];Counts", 100,-12.,12.,100,-12.,12.); \r
+ TH2F * hHitsZPAh = new TH2F("hHitsZPAh", "Hits centroid in ZPA;Centroid position [cm];Counts", 100,-12.,12.,100,-12.,12.); \r
+ Add2HitsList(hHitsZPCh, 14, !expert, image);\r
+ Add2HitsList(hHitsZPAh, 15, !expert, image);\r
}\r
\r
\r
Add2DigitsList(hDigPMCZPA, 11, expert, !image);\r
// \r
// ------------------- LOW GAIN CHAIN ---------------------------\r
-/* TH1F * hDigZNCTotlg = new TH1F("hDigZNCTotlg", "Digit lg signal in ZNC", 100, 0., 6000.);\r
+ TH1F * hDigZNCTotlg = new TH1F("hDigZNCTotlg", "Digit lg signal in ZNC", 100, 0., 6000.);\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
Add2DigitsList(hDigPMCZNAlg, 21, expert, !image);\r
Add2DigitsList(hDigPMCZPClg, 22, expert, !image);\r
Add2DigitsList(hDigPMCZPAlg, 23, expert, !image);\r
-*/\r
+\r
}\r
\r
//____________________________________________________________________________\r
-void AliZDCQADataMakerSim::MakeHits(TClonesArray * /*data*/)\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(fHits); \r
- AliZDCHit * hit; \r
- while((hit = dynamic_cast<AliZDCHit *>(next()))){\r
- if(hit->GetVolume(0)==1) GetHitsData(0)->Fill(hit->GetXImpact(),hit->GetYImpact());\r
- else if(hit->GetVolume(0)==4) GetHitsData(1)->Fill(hit->GetXImpact(), hit->GetYImpact());\r
-// else if(hit->GetVolume(0)==2) GetHitsData(1)->Fill(hit->GetXImpact(), hit->GetYImpact());\r
-// else if(hit->GetVolume(0)==5) GetHitsData(3)->Fill(hit->GetXImpact(), hit->GetYImpact());\r
+ TIter next(fHitsArray); \r
+ AliZDCHit * hit; \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
+ while((hit = dynamic_cast<AliZDCHit *>(next()))){\r
+ if(hit->GetVolume(0)==1){\r
+ adcSumQ_ZNC += hit->GetLightPMQ();\r
+ adcSum_ZNC += hit->GetLightPMC() + hit->GetLightPMQ();\r
+ //\r
+ GetHitsData(8)->Fill(hit->GetLightPMQ());\r
+ //\r
+ GetHitsData(12)->Fill(hit->GetXImpact(),hit->GetYImpact()); \r
+ }\r
+ else if(hit->GetVolume(0)==4){\r
+ adcSumQ_ZNA += hit->GetLightPMQ();\r
+ adcSum_ZNA += hit->GetLightPMC() + hit->GetLightPMQ();\r
+ //\r
+ GetHitsData(9)->Fill(hit->GetLightPMQ());\r
+ //\r
+ GetHitsData(13)->Fill(hit->GetXImpact(), hit->GetYImpact());\r
+ }\r
+ else if(hit->GetVolume(0)==2){\r
+ adcSumQ_ZNC += hit->GetLightPMQ();\r
+ adcSum_ZNC += hit->GetLightPMC() + hit->GetLightPMQ();\r
+ //\r
+ GetHitsData(10)->Fill(hit->GetLightPMQ());\r
+ //\r
+ GetHitsData(14)->Fill(hit->GetXImpact(), hit->GetYImpact());\r
}\r
+ else if(hit->GetVolume(0)==5){\r
+ adcSumQ_ZNC += hit->GetLightPMQ();\r
+ adcSum_ZNC += hit->GetLightPMC() + hit->GetLightPMQ();\r
+ //\r
+ GetHitsData(11)->Fill(hit->GetLightPMQ());\r
+ //\r
+ GetHitsData(15)->Fill(hit->GetXImpact(), hit->GetYImpact());\r
+ }\r
+ //\r
+ GetHitsData(0)->Fill(adcSum_ZNC);\r
+ GetHitsData(1)->Fill(adcSum_ZNA);\r
+ GetHitsData(2)->Fill(adcSum_ZPC);\r
+ GetHitsData(3)->Fill(adcSum_ZPA);\r
+ //\r
+ GetHitsData(4)->Fill(adcSumQ_ZNC);\r
+ GetHitsData(5)->Fill(adcSumQ_ZNA);\r
+ GetHitsData(6)->Fill(adcSumQ_ZPC);\r
+ GetHitsData(7)->Fill(adcSumQ_ZPA);\r
+ }\r
}\r
\r
//___________________________________________________________________________\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
- char** add = (char**) (branch->GetAddress());\r
- if(add){\r
- fHits = (TClonesArray*)(*add);\r
- } \r
- else{\r
- if(!fHits) fHits = new TClonesArray("AliZDCHit", 1000);\r
- branch->SetAddress(&fHits);\r
- }\r
- Int_t ntracks = (Int_t) hitTree->GetEntries();\r
- //printf("\n\t *** no.track %d\n",ntracks);\r
- if (ntracks<=0) return;\r
- //\r
- for(Int_t itrack=0; itrack<ntracks; itrack++){\r
- \r
- branch->GetEntry(itrack);\r
- //\r
- //printf("\t *** track %d",itrack);\r
- //hits->Print("");\r
- //printf("\n");\r
- //\r
- MakeHits(); \r
- fHits->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
+ 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
+ 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(6)->Fill(adcSumQ_ZPC);\r
GetDigitsData(7)->Fill(adcSumQ_ZPA);\r
//\r
- /*GetDigitsData(12)->Fill(adcSum_ZNC_lg);\r
+ GetDigitsData(12)->Fill(adcSum_ZNC_lg);\r
GetDigitsData(13)->Fill(adcSum_ZNA_lg);\r
GetDigitsData(14)->Fill(adcSum_ZPC_lg);\r
GetDigitsData(15)->Fill(adcSum_ZPA_lg);\r
GetDigitsData(16)->Fill(adcSumQ_ZNC_lg);\r
GetDigitsData(17)->Fill(adcSumQ_ZNA_lg);\r
GetDigitsData(18)->Fill(adcSumQ_ZPC_lg);\r
- GetDigitsData(19)->Fill(adcSumQ_ZPA_lg);*/\r
+ GetDigitsData(19)->Fill(adcSumQ_ZPA_lg);\r
+\r
}\r
\r
//____________________________________________________________________________\r
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff