X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=ZDC%2FAliZDCQADataMakerSim.cxx;h=949fbcc2f1e128b73197c961c93abe7ffc42e706;hb=62e9025c1ca2cce61f873a05630531ea451df623;hp=7892e9a1981600c7f3a9e884fa21b9ff42b0cdbf;hpb=6252ceeb59406d4f518422dfad532f414e0a23c4;p=u%2Fmrichter%2FAliRoot.git diff --git a/ZDC/AliZDCQADataMakerSim.cxx b/ZDC/AliZDCQADataMakerSim.cxx index 7892e9a1981..949fbcc2f1e 100644 --- a/ZDC/AliZDCQADataMakerSim.cxx +++ b/ZDC/AliZDCQADataMakerSim.cxx @@ -28,22 +28,19 @@ #include "AliZDCQADataMakerSim.h" #include "AliZDCHit.h" #include "AliZDCDigit.h" -#include "AliZDCRawStream.h" ClassImp(AliZDCQADataMakerSim) //____________________________________________________________________________ AliZDCQADataMakerSim::AliZDCQADataMakerSim() : - AliQADataMakerSim(AliQAv1::GetDetName(AliQAv1::kZDC), "ZDC Quality Assurance Data Maker"), - fDigit(0) + AliQADataMakerSim(AliQAv1::GetDetName(AliQAv1::kZDC), "ZDC Quality Assurance Data Maker") { // ctor } //____________________________________________________________________________ AliZDCQADataMakerSim::AliZDCQADataMakerSim(const AliZDCQADataMakerSim& qadm) : - AliQADataMakerSim(), - fDigit(0) + AliQADataMakerSim() { //copy ctor SetName((const char*)qadm.GetName()); @@ -66,16 +63,35 @@ void AliZDCQADataMakerSim::InitHits() // const Bool_t expert = kTRUE ; const Bool_t image = kTRUE ; + + TH1F * hHitsZNCTot = new TH1F("hHitsZNCTot", "Signal in ZNC; N_{phe}", 100, 0., 6000.); + TH1F * hHitsZNATot = new TH1F("hHitsZNATot", "Signal in ZNA; N_{phe}", 100, 0., 6000.); + TH1F * hHitsZPCTot = new TH1F("hHitsZPCTot", "Signal in ZPC; N_{phe}", 100, 0., 6000.); + TH1F * hHitsZPATot = new TH1F("hHitsZPATot", "Signal in ZPA; N_{phe}", 100, 0., 6000.); + Add2HitsList(hHitsZNCTot, 0, !expert, image); + Add2HitsList(hHitsZNATot, 1, !expert, image); + Add2HitsList(hHitsZPCTot, 2, !expert, image); + Add2HitsList(hHitsZPATot, 3, !expert, image); + // + TH1F * hHitsSumQZNC = new TH1F("hHitsSumQZNC", "Signal in 4 ZNC PMQ; N_{phe}",100, 0., 4000.); + TH1F * hHitsSumQZNA = new TH1F("hHitsSumQZNA", "Signal in 4 ZNA PMQ; N_{phe}",100, 0., 4000.); + TH1F * hHitsSumQZPC = new TH1F("hHitsSumQZPC", "Signal in 4 ZPC PMQ; N_{phe}",100, 0., 4000.); + TH1F * hHitsSumQZPA = new TH1F("hHitsSumQZPA", "Signal in 4 ZPA PMQ; N_{phe}",100, 0., 4000.); + Add2HitsList(hHitsSumQZNC, 4, expert, !image); + Add2HitsList(hHitsSumQZNA, 5, expert, !image); + Add2HitsList(hHitsSumQZPC, 6, expert, !image); + Add2HitsList(hHitsSumQZPA, 7, expert, !image); + // + TH1F * hHitsPMCZNC = new TH1F("hHitsPMCZNC", "Signal in ZNC PMC; N_{phe}",100, 0., 4000.); + TH1F * hHitsPMCZNA = new TH1F("hHitsPMCZNA", "Signal in ZNA PMC; N_{phe}",100, 0., 4000.); + TH1F * hHitsPMCZPC = new TH1F("hHitsPMCZPC", "Signal in ZPC PMC; N_{phe}",100, 0., 4000.); + TH1F * hHitsPMCZPA = new TH1F("hHitsPMCZPA", "Signal in ZPA PMC; N_{phe}",100, 0., 4000.); + Add2HitsList(hHitsPMCZNC, 8, expert, !image); + Add2HitsList(hHitsPMCZNA, 9, expert, !image); + Add2HitsList(hHitsPMCZPC, 10, expert, !image); + Add2HitsList(hHitsPMCZPA, 11, expert, !image); - TH2F * hZNCh = new TH2F("hZNCh", "Hits centroid in ZNC;Centroid position [cm];Counts", 100, -5.,5.,100,-5.,5.); - TH2F * hZNAh = new TH2F("hZNAh", "Hits centroid in ZNA;Centroid position [cm];Counts", 100, -5.,5.,100,-5.,5.); - // NB -> For the moment no check is performesd on ZP centroids -// TH2F * hZPCh = new TH2F("hZPCh", "Hits centroid in ZPC", 100,-12.,12.,100,-12.,12.); -// TH2F * hZPAh = new TH2F("hZPAh", "Hits centroid in ZPA", 100,-12.,12.,100,-12.,12.); - Add2HitsList(hZNCh, 0, !expert, image); - Add2HitsList(hZNAh, 1, !expert, image); -// Add2HitsList(hZPCh, 2, !expert, image); -// Add2HitsList(hZPAh, 3, !expert, image); + ClonePerTrigClass(AliQAv1::kHITS); // this should be the last line } @@ -115,66 +131,68 @@ void AliZDCQADataMakerSim::InitDigits() Add2DigitsList(hDigPMCZPC, 10, expert, !image); Add2DigitsList(hDigPMCZPA, 11, expert, !image); // - // ------------------- LOW GAIN CHAIN --------------------------- -/* TH1F * hDigZNCTotlg = new TH1F("hDigZNCTotlg", "Digit lg signal in ZNC", 100, 0., 6000.); - TH1F * hDigZNATotlg = new TH1F("hDigZNATotlg", "Digit lg signal in ZNA", 100, 0., 6000.); - TH1F * hDigZPCTotlg = new TH1F("hDigZPCTotlg", "Digit lg signal in ZPC", 100, 0., 6000.); - TH1F * hDigZPATotlg = new TH1F("hDigZPATotlg", "Digit lg signal in ZPA", 100, 0., 6000.); - Add2DigitsList(hDigZNCTotlg, 12, !expert, image); - Add2DigitsList(hDigZNATotlg, 13, !expert, image); - Add2DigitsList(hDigZPCTotlg, 14, !expert, image); - Add2DigitsList(hDigZPATotlg, 15, !expert, image); // - TH1F * hDigSumQZNClg = new TH1F("hDigSumQZNClg", "Signal in 4 ZNC PMQlg",100, 0., 4000.); - TH1F * hDigSumQZNAlg = new TH1F("hDigSumQZNAlg", "Signal in 4 ZNA PMQlg",100, 0., 4000.); - TH1F * hDigSumQZPClg = new TH1F("hDigSumQZPClg", "Signal in 4 ZPC PMQlg",100, 0., 4000.); - TH1F * hDigSumQZPAlg = new TH1F("hDigSumQZPAlg", "Signal in 4 ZPA PMQlg",100, 0., 4000.); - Add2DigitsList(hDigSumQZNClg, 16, expert, !image); - Add2DigitsList(hDigSumQZNAlg, 17, expert, !image); - Add2DigitsList(hDigSumQZPClg, 18, expert, !image); - Add2DigitsList(hDigSumQZPAlg, 19, expert, !image); - // - TH1F * hDigPMCZNClg = new TH1F("hDigPMCZNClg", "Signal in ZNC PMClg",100, 0., 4000.); - TH1F * hDigPMCZNAlg = new TH1F("hDigPMCZNAlg", "Signal in ZNA PMClg",100, 0., 4000.); - TH1F * hDigPMCZPClg = new TH1F("hDigPMCZPClg", "Signal in ZPC PMClg",100, 0., 4000.); - TH1F * hDigPMCZPAlg = new TH1F("hDigPMCZPAlg", "Signal in ZPA PMClg",100, 0., 4000.); - Add2DigitsList(hDigPMCZNClg, 20, expert, !image); - Add2DigitsList(hDigPMCZNAlg, 21, expert, !image); - Add2DigitsList(hDigPMCZPClg, 22, expert, !image); - Add2DigitsList(hDigPMCZPAlg, 23, expert, !image); -*/ + ClonePerTrigClass(AliQAv1::kDIGITS); // this should be the last line } //____________________________________________________________________________ void AliZDCQADataMakerSim::MakeHits() { //filling QA histos for Hits - // // Check id histograms already created for this Event Specie - if ( ! GetHitsData(0) ) - InitHits() ; + if( !GetHitsData(0) ) InitHits(); TIter next(fHitsArray); - AliZDCHit * hit; - while((hit = dynamic_cast(next()))){ - if(hit->GetVolume(0)==1) GetHitsData(0)->Fill(hit->GetXImpact(),hit->GetYImpact()); - else if(hit->GetVolume(0)==4) GetHitsData(1)->Fill(hit->GetXImpact(), hit->GetYImpact()); -// else if(hit->GetVolume(0)==2) GetHitsData(1)->Fill(hit->GetXImpact(), hit->GetYImpact()); -// else if(hit->GetVolume(0)==5) GetHitsData(3)->Fill(hit->GetXImpact(), hit->GetYImpact()); + AliZDCHit * hit; + Float_t adcSum_ZNC=0., adcSum_ZNA=0., adcSum_ZPC=0., adcSum_ZPA=0.; + Float_t adcSumQ_ZNC=0., adcSumQ_ZNA=0., adcSumQ_ZPC=0., adcSumQ_ZPA=0.; + while((hit = dynamic_cast(next()))){ + if(hit->GetVolume(0)==1){ + adcSumQ_ZNC += hit->GetLightPMQ(); + adcSum_ZNC += hit->GetLightPMC() + hit->GetLightPMQ(); + // + FillHitsData(8,hit->GetLightPMC()); + } + else if(hit->GetVolume(0)==4){ + adcSumQ_ZNA += hit->GetLightPMQ(); + adcSum_ZNA += hit->GetLightPMC() + hit->GetLightPMQ(); + // + FillHitsData(9,hit->GetLightPMC()); + } + else if(hit->GetVolume(0)==2){ + adcSumQ_ZNC += hit->GetLightPMQ(); + adcSum_ZNC += hit->GetLightPMC() + hit->GetLightPMQ(); + // + FillHitsData(10,hit->GetLightPMC()); } + else if(hit->GetVolume(0)==5){ + adcSumQ_ZNC += hit->GetLightPMQ(); + adcSum_ZNC += hit->GetLightPMC() + hit->GetLightPMQ(); + // + FillHitsData(11,hit->GetLightPMC()); + } + // + FillHitsData(0,adcSum_ZNC); + FillHitsData(1,adcSum_ZNA); + FillHitsData(2,adcSum_ZPC); + FillHitsData(3,adcSum_ZPA); + // + FillHitsData(4,adcSumQ_ZNC); + FillHitsData(5,adcSumQ_ZNA); + FillHitsData(6,adcSumQ_ZPC); + FillHitsData(7,adcSumQ_ZPA); + } } //___________________________________________________________________________ void AliZDCQADataMakerSim::MakeHits(TTree * hitTree) { // make QA data from Hit Tree - // if(!hitTree){ - AliError("Hit Tree not found!"); - return; - } - // + AliError("Can't get ZDC hit tree!!"); + return; + } TBranch * branch = hitTree->GetBranch("ZDC") ; @@ -182,133 +200,118 @@ void AliZDCQADataMakerSim::MakeHits(TTree * hitTree) AliError("ZDC branch in Hit Tree not found!"); return; } - else{ - if (fHitsArray) - fHitsArray->Clear() ; - char** add = (char**) (branch->GetAddress()); - if(add){ - fHitsArray = (TClonesArray*)(*add); - } - else{ - if(!fHitsArray) fHitsArray = new TClonesArray("AliZDCHit", 1000); - branch->SetAddress(&fHitsArray); - } - Int_t ntracks = (Int_t) hitTree->GetEntries(); - //printf("\n\t *** no.track %d\n",ntracks); - if (ntracks<=0) return; - // - for(Int_t itrack=0; itrackGetEntry(itrack); - // - //printf("\t *** track %d",itrack); - //hits->Print(""); - //printf("\n"); - // - MakeHits(); - fHitsArray->Clear(); - } - } + + if(fHitsArray) fHitsArray->Clear() ; + else fHitsArray = new TClonesArray("AliZDCHit", 1000); + + branch->SetAddress(&fHitsArray) ; + for (Int_t ientry = 0 ; ientry < branch->GetEntries() ; ientry++) { + branch->GetEntry(ientry) ; + MakeHits() ; + fHitsArray->Clear() ; + } + // + IncEvCountCycleHits(); + IncEvCountTotalHits(); + // } //___________________________________________________________________________ -void AliZDCQADataMakerSim::MakeDigits(TTree *digitTree ) +void AliZDCQADataMakerSim::MakeDigits(TTree *digitTree) { // makes data from Digit Tree + if( !GetDigitsData(0) ) InitDigits(); + + if(!digitTree){ + AliError("Can't get ZDC digit tree!!"); + return; + } + TBranch * branch = digitTree->GetBranch("ZDC"); if(!branch){ - AliError("ZDC branch in Digit Tree not found"); + AliError("ZDC branch in digit tree not found"); return; } - - // Check id histograms already created for this Event Specie - if ( ! GetDigitsData(0) ) - InitDigits() ; - - branch->SetAddress(&fDigit); - - Int_t ndig = digitTree->GetEntries(); - + + AliZDCDigit *digit = 0x0; + branch->SetAddress(&digit); + Float_t adcSum_ZNC=0., adcSum_ZNA=0., adcSum_ZPC=0., adcSum_ZPA=0.; Float_t adcSumQ_ZNC=0., adcSumQ_ZNA=0., adcSumQ_ZPC=0., adcSumQ_ZPA=0.; - //Float_t adcSum_ZNC_lg=0., adcSum_ZNA_lg=0., adcSum_ZPC_lg=0., adcSum_ZPA_lg=0.; - //Float_t adcSumQ_ZNC_lg=0., adcSumQ_ZNA_lg=0., adcSumQ_ZPC_lg=0., adcSumQ_ZPA_lg=0.; - // - for(Int_t i = 0; i < ndig; i++){ - digitTree->GetEntry(i); - if(fDigit->GetSector(0)==1){ - adcSum_ZNC += fDigit->GetADCValue(0); - //adcSum_ZNC_lg += fDigit->GetADCValue(1); + // Float_t adcSum_ZNC_lg=0., adcSum_ZNA_lg=0., adcSum_ZPC_lg=0., adcSum_ZPA_lg=0.; + // Float_t adcSumQ_ZNC_lg=0., adcSumQ_ZNA_lg=0., adcSumQ_ZPC_lg=0., adcSumQ_ZPA_lg=0.; + + Int_t ndig = digitTree->GetEntries(); + for(Int_t i=0; iGetEntry(i); + + if(digit->GetSector(0)==1 && digit->GetSector(1)!=5){ + adcSum_ZNC += digit->GetADCValue(0); + //adcSum_ZNC_lg += digit->GetADCValue(1); // - if(fDigit->GetSector(1)!=0){ - adcSumQ_ZNC += fDigit->GetADCValue(0); - //adcSumQ_ZNC_lg+= fDigit->GetADCValue(1); + if(digit->GetSector(1)!=0){ + adcSumQ_ZNC += digit->GetADCValue(0); + //adcSumQ_ZNC_lg+= digit->GetADCValue(1); } else{ - GetDigitsData(8)->Fill(fDigit->GetADCValue(0)); - //GetDigitsData(20)->Fill(fDigit->GetADCValue(1)); + FillDigitsData(8,digit->GetADCValue(0)); + //FillDigitsData(20,digit->GetADCValue(1)); } } - else if(fDigit->GetSector(0)==2){ - adcSum_ZPC += fDigit->GetADCValue(0); - //adcSum_ZPC_lg += fDigit->GetADCValue(1); + else if(digit->GetSector(0)==2){ + adcSum_ZPC += digit->GetADCValue(0); + //adcSum_ZPC_lg += digit->GetADCValue(1); // - if(fDigit->GetSector(1)!=0){ - adcSumQ_ZPC += fDigit->GetADCValue(0); - //adcSumQ_ZPC_lg+= fDigit->GetADCValue(1); + if(digit->GetSector(1)!=0){ + adcSumQ_ZPC += digit->GetADCValue(0); + //adcSumQ_ZPC_lg+= digit->GetADCValue(1); } else{ - GetDigitsData(10)->Fill(fDigit->GetADCValue(0)); - //GetDigitsData(22)->Fill(fDigit->GetADCValue(1)); + FillDigitsData(10,digit->GetADCValue(0)); + //FillDigitsData(22,digit->GetADCValue(1)); } } - else if(fDigit->GetSector(0)==4){ - adcSum_ZNA += fDigit->GetADCValue(0); - //adcSum_ZNA_lg += fDigit->GetADCValue(1); + else if(digit->GetSector(0)==4 && digit->GetSector(1)!=5){ + adcSum_ZNA += digit->GetADCValue(0); + //adcSum_ZNA_lg += digit->GetADCValue(1); // - if(fDigit->GetSector(1)!=0){ - adcSumQ_ZNA += fDigit->GetADCValue(0); - //adcSumQ_ZNA_lg+= fDigit->GetADCValue(1); + if(digit->GetSector(1)!=0){ + adcSumQ_ZNA += digit->GetADCValue(0); + //adcSumQ_ZNA_lg+= digit->GetADCValue(1); } else{ - GetDigitsData(9)->Fill(fDigit->GetADCValue(0)); - //GetDigitsData(21)->Fill(fDigit->GetADCValue(1)); + FillDigitsData(9,digit->GetADCValue(0)); + //FillDigitsData(21,digit->GetADCValue(1)); } } - else if(fDigit->GetSector(0)==5){ - adcSum_ZPA += fDigit->GetADCValue(0); - //adcSum_ZPA_lg += fDigit->GetADCValue(1); + else if(digit->GetSector(0)==5){ + adcSum_ZPA += digit->GetADCValue(0); + //adcSum_ZPA_lg += digit->GetADCValue(1); // - if(fDigit->GetSector(1)!=0){ - adcSumQ_ZPA += fDigit->GetADCValue(0); - //adcSumQ_ZPA_lg+= fDigit->GetADCValue(1); + if(digit->GetSector(1)!=0){ + adcSumQ_ZPA += digit->GetADCValue(0); + //adcSumQ_ZPA_lg+= digit->GetADCValue(1); } else{ - GetDigitsData(11)->Fill(fDigit->GetADCValue(0)); - //GetDigitsData(23)->Fill(fDigit->GetADCValue(1)); + FillDigitsData(11,digit->GetADCValue(0)); + //FillDigitsData(23,digit->GetADCValue(1)); } } } // - GetDigitsData(0)->Fill(adcSum_ZNC); - GetDigitsData(1)->Fill(adcSum_ZNA); - GetDigitsData(2)->Fill(adcSum_ZPC); - GetDigitsData(3)->Fill(adcSum_ZPA); + FillDigitsData(0,adcSum_ZNC); + FillDigitsData(1,adcSum_ZNA); + FillDigitsData(2,adcSum_ZPC); + FillDigitsData(3,adcSum_ZPA); // - GetDigitsData(4)->Fill(adcSumQ_ZNC); - GetDigitsData(5)->Fill(adcSumQ_ZNA); - GetDigitsData(6)->Fill(adcSumQ_ZPC); - GetDigitsData(7)->Fill(adcSumQ_ZPA); + FillDigitsData(4,adcSumQ_ZNC); + FillDigitsData(5,adcSumQ_ZNA); + FillDigitsData(6,adcSumQ_ZPC); + FillDigitsData(7,adcSumQ_ZPA); // - /*GetDigitsData(12)->Fill(adcSum_ZNC_lg); - GetDigitsData(13)->Fill(adcSum_ZNA_lg); - GetDigitsData(14)->Fill(adcSum_ZPC_lg); - GetDigitsData(15)->Fill(adcSum_ZPA_lg); + IncEvCountCycleDigits(); + IncEvCountTotalDigits(); // - GetDigitsData(16)->Fill(adcSumQ_ZNC_lg); - GetDigitsData(17)->Fill(adcSumQ_ZNA_lg); - GetDigitsData(18)->Fill(adcSumQ_ZPC_lg); - GetDigitsData(19)->Fill(adcSumQ_ZPA_lg);*/ } //____________________________________________________________________________ @@ -323,5 +326,6 @@ void AliZDCQADataMakerSim::EndOfDetectorCycle(AliQAv1::TASKINDEX_t task, TObjArr { // Detector specific actions at end of cycle // do the QA checking + ResetEventTrigClasses(); AliQAChecker::Instance()->Run(AliQAv1::kZDC, task, list); }