- Int_t iDetElt = 0;
- if (i<16) iDetElt = 100*(i/4+1) + i - int(i/4)*4;
- else if (i<52) iDetElt = 100*((i-16)/18+5) + i - (int((i-16)/18)*18 + 16);
- if (i>51) iDetElt = 100*((i-52)/26+7) + i - (int((i-52)/26)*26 + 52);
-
- Char_t histName[255];
- Char_t histTitle[255];
-
- sprintf(histName,"Eff_detEltNbr%d",iDetElt);
- sprintf(histTitle,"detEltNbr %d",iDetElt);
- if(i<16) new((*fDetEltEffHistList)[i]) TH2F(histName,histTitle,12,-10,110,12,-10,110);//!<Stations 1 & 2.
- else new((*fDetEltEffHistList)[i]) TH2F(histName,histTitle,28,-140,140,8,-40,40); //!<Stations 3 -> 5.
-
- sprintf(histName,"TD_detEltNbr%d",iDetElt);
- if(i<16) new((*fDetEltTDHistList)[i]) TH2F(histName,histTitle,12,-10,110,12,-10,110);//!<Stations 1 & 2.
- else new((*fDetEltTDHistList)[i]) TH2F(histName,histTitle,28,-140,140,8,-40,40); //!<Stations 3 -> 5.
-
- sprintf(histName,"TT_detEltNbr%d",iDetElt);
- if(i<16) new((*fDetEltTTHistList)[i]) TH2F(histName,histTitle,12,-10,110,12,-10,110); //!<Stations 1 & 2.
- else new((*fDetEltTTHistList)[i]) TH2F(histName,histTitle,28,-140,140,8,-40,40); //!<Stations 3 -> 5.
-
- ((TH2F*) fDetEltEffHistList->UncheckedAt(i)) -> GetXaxis() -> SetTitle("X (cm)");
- ((TH2F*) fDetEltEffHistList->UncheckedAt(i)) -> GetYaxis() -> SetTitle("Y (cm)");
- ((TH2F*) fDetEltEffHistList->UncheckedAt(i)) -> GetZaxis() -> SetTitle("Efficiency (%)");
- ((TH2F*) fDetEltEffHistList->UncheckedAt(i)) -> GetXaxis() -> SetTitleOffset(1.8);
- ((TH2F*) fDetEltEffHistList->UncheckedAt(i)) -> GetYaxis() -> SetTitleOffset(1.8);
- ((TH2F*) fDetEltEffHistList->UncheckedAt(i)) -> GetZaxis() -> SetTitleOffset(1.2);
- ((TH2F*) fDetEltEffHistList->UncheckedAt(i)) -> SetOption("LEGO");
- }
-
- for (Int_t j = 0; j < fTotNbrOfChamber; j++)
- {
- Char_t histName[255];
- Char_t histTitle[255];
-
-
- sprintf(histName,"Eff_ChamberNbr%d",j+1);
- sprintf(histTitle,"ChamberNbr %d",j+1);
- if (j<4) new ((*fChamberEffHistList)[j]) TH1F(histName, histTitle, 4, 0.0, 4.0);
- else if (j<6) new ((*fChamberEffHistList)[j]) TH1F(histName, histTitle, 18, 0.0, 18.0);
- if (j>=6) new ((*fChamberEffHistList)[j]) TH1F(histName, histTitle, 26, 0.0, 26.0);
-
- sprintf(histName,"TD_ChamberNbr%d",j+1);
- if (j<4) new ((*fChamberTDHistList)[j]) TH1F(histName, histTitle, 4, 0.0, 4.0);
- else if (j<6) new ((*fChamberTDHistList)[j]) TH1F(histName, histTitle, 18, 0.0, 18.0);
- if (j>=6) new ((*fChamberTDHistList)[j]) TH1F(histName, histTitle, 26, 0.0, 26.0);
-
- sprintf(histName,"TT_ChamberNbr%d",j+1);
- if (j<4) new ((*fChamberTTHistList)[j]) TH1F(histName, histTitle, 4, 0.0, 4.0);
- else if (j<6) new ((*fChamberTTHistList)[j]) TH1F(histName, histTitle, 18, 0.0, 18.0);
- if (j>=6) new ((*fChamberTTHistList)[j]) TH1F(histName, histTitle, 26, 0.0, 26.0);
-
- ((TH1F*) fChamberEffHistList->UncheckedAt(j)) -> GetXaxis() -> SetTitle("DetElement");
- ((TH1F*) fChamberEffHistList->UncheckedAt(j)) -> GetYaxis() -> SetTitle("Efficiency (%)");
- ((TH1F*) fChamberEffHistList->UncheckedAt(j)) -> GetXaxis() -> SetTitleOffset(1.8);
- ((TH1F*) fChamberEffHistList->UncheckedAt(j)) -> GetYaxis() -> SetTitleOffset(1.8);
- ((TH1F*) fChamberEffHistList->UncheckedAt(j)) -> Sumw2();
+ Int_t deId = FromLocalId2DetElt(iCh, iDE);
+ histTitle.Form("detEltNbr %d",deId);
+ if(iCh < 4)
+ {// chambers 1 -> 4
+ histName.Form("TD_detEltNbr%d",deId);
+ h3 = new TH3F(histName, histTitle, 12, -10.0 , 110.0, 12, -10.0, 110.0, nCentBins, centRange[0], centRange[1]);
+ fDetEltTDHistList->AddAt(h3, iDEGlobal);
+ histName.Form("TT_detEltNbr%d",deId);
+ h3 = new TH3F(histName, histTitle, 12, -10.0 , 110.0, 12, -10.0, 110.0, nCentBins, centRange[0], centRange[1]);
+ fDetEltTTHistList->AddAt(h3, iDEGlobal);
+ histName.Form("SD_detEltNbr%d",deId);
+ h3 = new TH3F(histName, histTitle, 12, -10.0 , 110.0, 12, -10.0, 110.0, nCentBins, centRange[0], centRange[1]);
+ fDetEltSDHistList->AddAt(h3, iDEGlobal);