+ //end of a QA cycle
+
+ AliDebug(AliQAv1::GetQADebugLevel(),"AliITSDM instantiates checker with Run(AliQAv1::kITS, task, list)\n");
+
+ Double_t entriescalibration1= 0.;
+ Double_t entriescalibrationL3=0.;
+ Double_t entriescalibrationL4=0.;
+ int offsRW = fGenRawsOffset[fAliITSQADataMakerRec->GetEventSpecie()];
+ int offsRP = fGenRecPointsOffset[fAliITSQADataMakerRec->GetEventSpecie()];
+ //
+ if(fHistoCalibration){
+ entriescalibration1= ((TH1F*)(fHistoCalibration->At(0)))->GetEntries();
+ entriescalibrationL3=((TH2F*)(fHistoCalibration->At(1)))->GetEntries();
+ entriescalibrationL4=((TH2F*)(fHistoCalibration->At(2)))->GetEntries();
+ }
+ else{ AliWarning("Calibration TObjArray is NULL! No Normalization and calibtaion plot will be filled\n");}
+ //
+ for (int trCl=-1;trCl<fAliITSQADataMakerRec->GetNTrigClasses();trCl++) { // RS Loop over all trigger classes (+ global counter, -1)
+ //
+ if(task==AliQAv1::kRAWS) {
+ TObjArray &harr = *fAliITSQADataMakerRec->GetRawsDataOfTrigClass(trCl);
+ int nEvent = GetNumberOfEvents(AliQAv1::kRAWS,trCl);
+ // printf("fNevent %d \n",nEvent);
+ TH1* h10 = (TH1*)harr[10+offsRW];
+ if(nEvent==0) if(h10) nEvent=h10->GetBinContent(1);
+ if (h10) h10->SetBinContent(1,nEvent);
+ //
+ if(nEvent!=0) {
+ if (h10 && harr[0+offsRW] && harr[1+offsRW] && harr[2]) { // were they all cloned?
+ //
+ Double_t entriesmodpattern = ((TH1*)harr[0+offsRW])->GetEntries();
+ Double_t entriesL3 = ((TH1*)harr[1+offsRW])->GetEntries();
+ Double_t entriesL4 = ((TH1*)harr[2+offsRW])->GetEntries();
+ //
+ Double_t normentriesmodpattern=0.;
+ Double_t normentriesL3= 0.;
+ Double_t normentriesL4= 0.;
+
+ if(entriescalibration1!=0.)normentriesmodpattern=entriesmodpattern/(entriescalibration1*nEvent);
+ if(entriesL3!=0.) normentriesL3= entriesL3/(entriescalibrationL3*nEvent);
+ if(entriesL4!=0.) normentriesL4= entriesL4/(entriescalibrationL4*nEvent);
+ //
+ h10->SetBinContent(2,normentriesmodpattern);
+ h10->SetBinContent(3,normentriesL3);
+ h10->SetBinContent(4,normentriesL4);
+ //
+ if(fHistoCalibration){
+ if (harr[3+offsRW] && harr[0+offsRW]) ((TH1*)harr[3+offsRW])->Divide( (TH1*)harr[0+offsRW],((TH1F*)(fHistoCalibration->At(0))),1.,nEvent);
+ if (harr[4+offsRW] && harr[1+offsRW]) ((TH1*)harr[4+offsRW])->Divide( (TH1*)harr[1+offsRW],((TH1F*)(fHistoCalibration->At(1))),1.,nEvent);
+ if (harr[5+offsRW] && harr[2+offsRW]) ((TH1*)harr[5+offsRW])->Divide( (TH1*)harr[2+offsRW],((TH1F*)(fHistoCalibration->At(2))),1.,nEvent);
+ }
+ }
+ }
+ //
+ TH1* htmp4 = (TH1*)harr[4+offsRW];
+ Int_t xbin3 = htmp4 ? ((TH1*)harr[4+offsRW])->GetNbinsX() : 0;
+ Int_t ybin3 = htmp4 ? ((TH1*)harr[4+offsRW])->GetNbinsY() : 0;
+ //
+ TH1* htmp6 = ((TH1*)harr[6+offsRW]);
+ if (htmp6 && htmp4) for(Int_t i=0; i<xbin3; i++) for(Int_t j=0; j<ybin3; j++) htmp6->Fill(htmp4->GetBinContent(i,j));
+ //
+ TH1* htmp5 = (TH1*)harr[5+offsRW];
+ Int_t xbin4 = htmp5 ? htmp5->GetNbinsX() : 0;
+ Int_t ybin4 = htmp5 ? htmp5->GetNbinsY() : 0;
+ //
+ TH1* htmp7 = (TH1*)harr[7+offsRW];
+ if (htmp7 && htmp5) for(Int_t i=0; i<xbin4; i++) for(Int_t j=0; j<ybin4; j++) htmp7->Fill(htmp5->GetBinContent(i,j));
+ //
+ if(fHistoCalibration){
+ for(Int_t inumb=8;inumb<10;inumb++) {
+ TH1* htmp = (TH1*)harr[inumb + offsRW];
+ if (!htmp) continue;
+ TH1* hcal = (TH1*)fHistoCalibration->At(inumb-7);
+ for(Int_t ii=1; ii<htmp->GetNbinsX()+1;ii++) {
+ for(Int_t jj=1; jj<htmp->GetNbinsY()+1;jj++) {
+ if( hcal->GetBinContent(ii,jj) != 0. ) htmp->SetBinContent(ii,jj,1);
+ else if( hcal->GetBinContent(ii,jj) == 0. ) htmp->SetBinContent(ii,jj,0);
+ //MM printf("x=%i \t y=%i => \t %f \n",ii,jj,hcal->GetBinContent(ii,jj));
+ }//end x axis
+ }//end y axis
+ }//end for inumb
+ }
+ //
+ }//end raws
+ //
+ //
+ if(task==AliQAv1::kRECPOINTS) {
+ //
+ TObjArray &harr = *fAliITSQADataMakerRec->GetRecPointsDataOfTrigClass(trCl);
+ //
+ TH1* htmp0 = (TH1*)harr[0+offsRP];
+ TH1* htmp1 = (TH1*)harr[1+offsRP];
+ TH1* htmp13 = (TH1*)harr[13+offsRP];
+ TH1* htmp14 = (TH1*)harr[14+offsRP];
+ //
+ int nEventRP = GetNumberOfEvents(AliQAv1::kRECPOINTS,trCl);
+ if (htmp0 && htmp1 && htmp13 && htmp14) {
+ Double_t chargeL3 = htmp0->GetMean();
+ Double_t errchargeL3 = htmp0->GetMeanError();
+ //
+ Double_t chargeL4 = htmp1->GetMean();
+ Double_t errchargeL4 = htmp1->GetMeanError();
+ //
+ Double_t radiusL3 = htmp13->GetMean();
+ Double_t errradiusL3 = htmp13->GetMeanError();
+
+ Double_t radiusL4= htmp14->GetMean();
+ Double_t errradiusL4= htmp14->GetMeanError();
+ //
+ // printf("fNeventRP %d \n",nEventRP);
+ TH1* htmp27 = (TH1*)harr[27+offsRP];
+ if(nEventRP==0) if(htmp27) nEventRP=htmp27->GetBinContent(1);
+
+ if (htmp27) {
+ htmp27->SetBinContent(1,nEventRP);
+
+ htmp27->SetBinContent(31,chargeL3);
+ htmp27->SetBinError(31,errchargeL3);
+
+ htmp27->SetBinContent(47,chargeL4);
+ htmp27->SetBinError(47,errchargeL4);
+
+ htmp27->SetBinContent(32,radiusL3);
+ htmp27->SetBinError(32,errradiusL3);
+
+ htmp27->SetBinContent(48,radiusL4);
+ htmp27->SetBinError(48,errradiusL4);
+ }
+ //
+ TH1 *htmp6=(TH1*)harr[6+offsRP],*htmp7=(TH1*)harr[7+offsRP],*htmp8=(TH1*)harr[8+offsRP];
+ Double_t entriesmodpattern = htmp6 ? htmp6->GetEntries() : 0;
+ Double_t entriesL3 = htmp7 ? htmp7->GetEntries() : 0;
+ Double_t entriesL4 = htmp8 ? htmp8->GetEntries() : 0;
+
+ //printf("entries modpatternrp = %02f \t L3= %02f \t L4 =%02f \n",entriesmodpattern,entriesL3,entriesL4 );
+
+ if(nEventRP!=0) {
+ Double_t normentriesmodpattern=0.;
+ Double_t normentriesL3= 0.;
+ Double_t normentriesL4= 0.;
+
+ if(entriescalibration1!=0.) normentriesmodpattern=entriesmodpattern/(entriescalibration1*nEventRP);
+ if(entriesL3!=0.) normentriesL3= entriesL3/(entriescalibrationL3*nEventRP);
+ if(entriesL4!=0.) normentriesL4= entriesL4/(entriescalibrationL4*nEventRP);
+
+ if (htmp27) {
+ htmp27->SetBinContent(2,normentriesmodpattern);
+ htmp27->SetBinContent(3,normentriesL3);
+ htmp27->SetBinContent(4,normentriesL4);
+ }
+ // printf("NORM entries modpatternrp = %02f \t L3= %02f \t L4 =%02f \n",normentriesmodpattern,normentriesL3,normentriesL4 );
+
+ if(fHistoCalibration) {
+ if (harr[9 +offsRP] && htmp6) ((TH1*)harr[9 +offsRP])->Divide(htmp6, (TH1*)fHistoCalibration->At(0),nEventRP);
+ if (harr[10+offsRP] && htmp7) ((TH1*)harr[10+offsRP])->Divide(htmp7, (TH1*)fHistoCalibration->At(1),nEventRP);
+ if (harr[11+offsRP] && htmp8) ((TH1*)harr[11+offsRP])->Divide(htmp8, (TH1*)fHistoCalibration->At(2),nEventRP);
+ //if (harr[21+offsRP] && htmp7) ((TH1*)harr[21+offsRP])->Divide(htmp7, (TH1*)fHistoCalibration->At(1),nEventRP);
+ //if (harr[22+offsRP] && htmp8) ((TH1*)harr[22+offsRP])->Divide(htmp8, (TH1*)fHistoCalibration->At(2),nEventRP);
+ }
+ }
+ TH1* htmp10 = (TH1*)harr[10+offsRP];
+ Int_t xbin3 = htmp10 ? htmp10->GetNbinsX() : 0;
+ Int_t ybin3 = htmp10 ? htmp10->GetNbinsY() : 0;
+ //
+ TH1* htmp19 = (TH1*)harr[19+offsRP];
+ if (htmp19&&htmp10) for(Int_t i=0; i<xbin3; i++) for(Int_t j=0; j<ybin3; j++) htmp19->Fill(htmp10->GetBinContent(i,j));
+ //
+ TH1* htmp11 = (TH1*)harr[11+offsRP];
+ Int_t xbin4 = htmp11 ? htmp11->GetNbinsX() : 0;
+ Int_t ybin4 = htmp11 ? htmp11->GetNbinsY() : 0;
+ //
+ TH1* htmp20 = (TH1*)harr[20+offsRP];
+ if (htmp11&&htmp20) for(Int_t i=0; i<xbin4; i++) for(Int_t j=0; j<ybin4; j++) htmp20->Fill(htmp11->GetBinContent(i,j));
+ }
+ //
+
+ // RecPoints 2 Raws Ratio
+ if(fAliITSQADataMakerRec->ListExists(AliQAv1::kRAWS)==kTRUE) {
+ //
+ TH1* htmp21 = (TH1*)harr[21+offsRP];
+ TH1* htmp23 = (TH1*)harr[23+offsRP];
+ TH1* htmp10 = (TH1*)harr[10+offsRP];
+ TH1* htmp4rw = fAliITSQADataMakerRec->GetRawsData(4 + offsRW,trCl);
+ if (htmp21 && htmp23 && htmp10 && htmp4rw) {
+ Int_t xbin3RP = htmp21->GetNbinsX();
+ Int_t ybin3RP = htmp21->GetNbinsY();
+ Int_t xbin3R = htmp4rw->GetNbinsX();
+ Int_t ybin3R = htmp4rw->GetNbinsY();
+ if((xbin3RP == xbin3R) && (ybin3RP == ybin3R)) {
+ htmp21->Divide(htmp10,htmp4rw);
+ for(Int_t i=0; i<xbin3R; i++) for(Int_t j=0; j<ybin3R; j++) htmp23->Fill(htmp21->GetBinContent(i,j));
+ }
+ else AliWarning("Number of bins for Raws and RecPoints (Layer 3) do not match\n");
+ }
+ //
+ TH1* htmp22 = (TH1*)harr[22+offsRP];
+ TH1* htmp24 = (TH1*)harr[24+offsRP];
+ TH1* htmp11 = (TH1*)harr[11+offsRP];
+ TH1* htmp5rw = fAliITSQADataMakerRec->GetRawsData(5 + offsRW,trCl);
+ if (htmp22 && htmp24 && htmp11 && htmp5rw) {
+ Int_t xbin4RP = htmp22->GetNbinsX();
+ Int_t ybin4RP = htmp22->GetNbinsY();
+ Int_t xbin4R = htmp5rw->GetNbinsX();
+ Int_t ybin4R = htmp5rw->GetNbinsY();
+ if((xbin4RP == xbin4R) && (ybin4RP == ybin4R)) {
+ htmp22->Divide(htmp11,htmp5rw);
+ for(Int_t i=0; i<xbin4R; i++) for(Int_t j=0; j<ybin4R; j++) htmp24->Fill(htmp22->GetBinContent(i,j));
+ }
+ else AliWarning("Number of bins for Raws and RecPoints (Layer 4) do not match\n");
+ }
+ }
+ else{AliWarning("Ratio between RecPoints and Raws not executed because the raw list has not been created\n");}
+ }//end recpoints
+ //
+ } // RS Loop over all trigger classes (+ global counter, -1)
+ //