- // loop over the trigger DDL (Right: 20, Left: 21)
- for (Int_t iDDL = 20; iDDL <= 21; ++iDDL) {
- AliMpDDL* ddl = ddlStore->GetDDL(iDDL);
- Int_t nCrate = ddl->GetNofTriggerCrates();
- // loop over the number of crates in DDL
- for (Int_t index = 0; index < nCrate; ++index) {
- // get crate object
- AliMpTriggerCrate* crate = ddlStore->GetTriggerCrate(iDDL, index);
- Int_t nLocal = crate->GetNofLocalBoards();
- for (Int_t iLocal = 0; iLocal < nLocal; ++iLocal) {
- // get local board Id from crate object
- localId = crate->GetLocalBoardId(iLocal);
- // get local board object
- AliMpLocalBoard* localBoard = ddlStore->GetLocalBoard(localId);
-
- if (!localBoard->IsNotified()) continue;
-
- // get detection element connected to this board
- for (Int_t ch = 0; ch < fgkNchambers; ++ch) {
- Int_t iCh = ch + AliMpConstants::NofTrackingChambers();
- Int_t detElemId = ddlStore->GetDEfromLocalBoard(localId, iCh);
-
- if (!detElemId) continue;
-
- // get segmentation
- for (Int_t cath = 0; cath < fgkNcathodes; ++cath) {
- const AliMpVSegmentation* seg = AliMpSegmentation::Instance()
- ->GetMpSegmentation(detElemId,
- AliMp::GetCathodType(cath));
-
- Int_t chCath = fgkNchambers*cath + ch;
- slat = detElemId%100;
- Int_t nStrips=0;
-
- // loop over strips
- for (Int_t ibitxy = 0; ibitxy < 16; ++ibitxy) {
- // get pad from electronics
- AliMpPad pad = seg->PadByLocation(AliMpIntPair(localId,
- ibitxy),kFALSE);
-
- if (!pad.IsValid()) continue;
- if(cath==0){ // Geometry info from bending plane only
- if(ibitxy==0) {
- xLocal1 = pad.Position().X();
- yLocal1 = pad.Position().Y();
- xWidth = pad.Dimensions().X();
- yWidth = pad.Dimensions().Y();
- }
- xLocal2 = pad.Position().X();
- yLocal2 = pad.Position().Y();
- }
-
- // Check fired pads (when available)
- if(padFired[cath]) {
- Int_t padX = pad.GetIndices().GetFirst();
- Int_t padY = pad.GetIndices().GetSecond();
- Int_t currPair = padX*maxY[cath] + padY;
- nStrips++;
- //printf("cath %i board = %i (%2i, %2i) -> %i\n",
- //cath, localId, padX, padY, currPair);
- Int_t chBin = padFired[cath]->GetXaxis()->FindBin(ch);
- Int_t slatBin = padFired[cath]->GetYaxis()->FindBin(slat);
- Int_t pairBin = padFired[cath]->GetZaxis()->FindBin(currPair);
- Float_t nFired = padFired[cath]->GetBinContent(chBin, slatBin, pairBin);;
-
- Float_t dimX = pad.Dimensions().X();
- Float_t dimY = pad.Dimensions().Y();
-
- Float_t stripX1 = pad.Position().X();
- Float_t stripY1 = pad.Position().Y();
- Float_t stripX2 = pad.Position().X();
- Float_t stripY2 = pad.Position().Y();
-
- transform->Local2Global(detElemId, stripX1, stripY1, 0, xg1, yg1, zg1);
- transform->Local2Global(detElemId, stripX2, stripY2, 0, xg2, yg2, zg2);
-
- Float_t x1Float = TMath::Min(xg1,xg2) - dimX;
- Float_t y1Float = TMath::Min(yg1,yg2) - dimY;
- Float_t x2Float = TMath::Max(xg1,xg2) + dimX;
- Float_t y2Float = TMath::Max(yg1,yg2) + dimY;
-
- Int_t x1Int = histoFiredDisplay[chCath]->GetXaxis()->FindBin(x1Float)+1;
- Int_t y1Int = histoFiredDisplay[chCath]->GetYaxis()->FindBin(y1Float)+1;
- Int_t x2Int = histoFiredDisplay[chCath]->GetXaxis()->FindBin(x2Float)-1;
- Int_t y2Int = histoFiredDisplay[chCath]->GetYaxis()->FindBin(y2Float)-1;
-
- for(Int_t binX=x1Int; binX<=x2Int; binX++){
- for(Int_t binY=y1Int; binY<=y2Int; binY++){
- histoFiredDisplay[chCath]->SetBinContent(binX, binY, nFired);
- }
- }