return kFALSE;
}
}
-
+
Int_t iCh = (chamber > AliMpConstants::NofTrackingChambers()) ? chamber - 11 : chamber;
Int_t iCath = cathode;
-
+
TArrayD xAxisStrip;
TArrayD yAxisStrip;
TArrayD xAxisBoard;
TArrayD yAxisBoard;
-
+
Float_t yOffsetLine, xOffsetLine = 0.;
-
+
const Float_t kResetValue=1234567.;
-
+
if(!inputObject){
xAxisBoard.Set(55);
xAxisBoard.Reset(kResetValue);
yAxisBoard.Set(50);
yAxisBoard.Reset(kResetValue);
-
+
xAxisStrip.Set(420);
xAxisStrip.Reset(kResetValue);
yAxisStrip.Set(710);
yAxisStrip.Reset(kResetValue);
}
else if(!displayHisto){
- AliWarning("Display histogram not initialized. Please initialize it first!");
- return kFALSE;
+ AliWarning("Display histogram not initialized. Please initialize it first!");
+ return kFALSE;
}
else {
- TH1* inputHisto = dynamic_cast<TH1*>(inputObject);
+ TH1* inputHisto = dynamic_cast<TH1*>(inputObject);
if ( inputHisto ) {
if ( inputHisto->GetEntries() == 0 ) {
- return kTRUE;
+ return kTRUE;
}
}
else {
TGraph* inputGraph = dynamic_cast<TGraph*>(inputObject);
if ( inputGraph ) {
if ( inputGraph->GetN() == 0 ){
- return kTRUE;
+ return kTRUE;
}
- }
+ }
else {
AliWarning("The object should inherit from TH1 or TGraph!");
return kFALSE;
}
}
}
-
+
Float_t xWidth, yWidth, yWidthSlat=0., xWidthCol=0.;
Float_t x1,x2,y1,y2;
Float_t x1b=0., x2b=0., y1b=0., y2b=0.;
Float_t xcPad, ycPad;
Int_t line=0, slat;
Float_t sign = 1.;
-
+
const Float_t kShiftB = 0.5;
const Float_t kShiftS = 0.1;
-
+
const Float_t kShiftX = (iCath==0) ? kShiftB : kShiftS;
const Float_t kShiftY = (iCath==0) ? kShiftS : kShiftB;
const Float_t kShiftEl = (displayType==kDisplaySlats) ? 0.01 : kShiftB;
-
+
Int_t iChamber = iCh + AliMpConstants::NofTrackingChambers();
- for(Int_t iLoc = 0; iLoc < AliMpConstants::NofLocalBoards(); iLoc++) {
+ for(Int_t iLoc = 0; iLoc < AliMpConstants::NofLocalBoards(); iLoc++) {
Int_t iBoard = iLoc+1;
Int_t detElemId = AliMpDDLStore::Instance()->GetDEfromLocalBoard(iBoard, iChamber);
-
+
if (!detElemId) continue;
-
+
AliMpLocalBoard* localBoard = AliMpDDLStore::Instance()->GetLocalBoard(iBoard, kFALSE);
-
+
// skip copy cards
- if( !localBoard->IsNotified())
+ if( !localBoard->IsNotified())
continue;
-
+
// get segmentation
const AliMpVSegmentation* seg[2] = {
AliMpSegmentation::Instance()->GetMpSegmentation(detElemId, AliMp::GetCathodType(0)),
AliMpSegmentation::Instance()->GetMpSegmentation(detElemId, AliMp::GetCathodType(1))};
-
+
if(iLoc==0){
AliMpPad pad1 = seg[1]->PadByLocation(iBoard,0,kFALSE);
yWidthSlat = pad1.GetDimensionY();
xOffsetLine = TMath::Abs(pad0.GetPositionX()) + pad0.GetDimensionX();
xWidthCol = 2.* pad0.GetDimensionX();
}
-
+
// Get ideal global position of DetElemId center
slat = detElemId%100;
line = (4 + slat)%18;
yOffsetLine = (Float_t)(line - 4) * 2. * yWidthSlat;
Int_t nLocations = 1;
-
+
for(Int_t cath=0; cath<AliMpConstants::NofCathodes(); cath++){
- // Loop on cathodes:
+ // Loop on cathodes:
// necessary because strip info is read for each cathode
// board info is read only from cathode 0
-
+
// loop over strips
for (Int_t ibitxy = 0; ibitxy < 16; ++ibitxy) {
- // get pad from electronics
-
- Int_t offset = 0;
- if (cath && localBoard->GetSwitch(AliMpLocalBoard::kZeroAllYLSB)) offset = -8;
-
- AliMpPad pad = seg[cath]->PadByLocation(iBoard,ibitxy+offset,kFALSE);
-
- if (!pad.IsValid()) continue;
+ // get pad from electronics
+
+ Int_t offset = 0;
+ if (cath && localBoard->GetSwitch(AliMpLocalBoard::kZeroAllYLSB)) offset = -8;
- // For non-bending plane fill only the first board covered by the strip
- // i.e. avoide filling many times the same information
- if ( cath == 1 && pad.GetLocalBoardId(0) != iBoard ) continue;
-
- xWidth = pad.GetDimensionX();
- yWidth = pad.GetDimensionY();
- xcPad = sign * (pad.GetPositionX() + xOffsetLine);
- if(line==4) xcPad += 0.75 * sign * xWidthCol;
- ycPad = pad.GetPositionY() + yOffsetLine;
- nLocations = pad.GetNofLocations();
- Int_t iStrip = pad.GetLocalBoardChannel(0);
-
- if(cath==iCath){
- x1 = xcPad - xWidth + kShiftX;
- y1 = ycPad - yWidth + kShiftY;
- x2 = xcPad + xWidth - kShiftX;
- y2 = ycPad + yWidth - kShiftY;
-
- if(!inputObject){
- AddSortedPoint(x1, xAxisStrip, kResetValue);
- AddSortedPoint(x2, xAxisStrip, kResetValue);
-
- AddSortedPoint(y1, yAxisStrip, kResetValue);
- AddSortedPoint(y2, yAxisStrip, kResetValue);
- }
- else if(displayType==kDisplayStrips)
- FillBins(inputObject, displayHisto, iBoard, iStrip, x1, x2, y1, y2, kShiftX, kShiftY, displayOpt);
- }
-
- if(cath==0){
- if(iStrip==0) {
- x1b = xcPad - xWidth + kShiftEl;
- y1b = ycPad - yWidth + kShiftEl;
- }
- x2b = xcPad + xWidth - kShiftEl;
- y2b = ycPad + yWidth - kShiftEl;
- }
+ AliMpPad pad = seg[cath]->PadByLocation(iBoard,ibitxy+offset,kFALSE);
+
+ if (!pad.IsValid()) continue;
+
+ xWidth = pad.GetDimensionX();
+ yWidth = pad.GetDimensionY();
+ xcPad = sign * (pad.GetPositionX() + xOffsetLine);
+ if(line==4) xcPad += 0.75 * sign * xWidthCol;
+ ycPad = pad.GetPositionY() + yOffsetLine;
+ nLocations = pad.GetNofLocations();
+ Int_t iStrip = pad.GetLocalBoardChannel(0);
+
+ if(cath==0){
+ if(iStrip==0) {
+ x1b = xcPad - xWidth + kShiftEl;
+ y1b = ycPad - yWidth + kShiftEl;
+ }
+ x2b = xcPad + xWidth - kShiftEl;
+ y2b = ycPad + yWidth - kShiftEl;
+ }
+
+ // For non-bending plane fill only the first board covered by the strip
+ // i.e. avoid filling many times the same information
+ if ( cath == 1 && pad.GetLocalBoardId(0) != iBoard ) continue;
+
+ if(cath==iCath){
+ x1 = xcPad - xWidth + kShiftX;
+ y1 = ycPad - yWidth + kShiftY;
+ x2 = xcPad + xWidth - kShiftX;
+ y2 = ycPad + yWidth - kShiftY;
+
+ if(!inputObject){
+ AddSortedPoint(x1, xAxisStrip, kResetValue);
+ AddSortedPoint(x2, xAxisStrip, kResetValue);
+
+ AddSortedPoint(y1, yAxisStrip, kResetValue);
+ AddSortedPoint(y2, yAxisStrip, kResetValue);
+ }
+ else if(displayType==kDisplayStrips)
+ FillBins(inputObject, displayHisto, iBoard, iStrip, x1, x2, y1, y2, kShiftX, kShiftY, displayOpt);
+ }
+
} // loop on strips
-
+
// if iCath==0 strip info and board info are both filled -> break!
// if iCath==1 board info is filled at cath==0. Strip info to be filled at cath==1
if(iCath==0) break;
} // loop on cathodes
-
+
if(!inputObject){
// Per board
AddSortedPoint(x1b, xAxisBoard, kResetValue);
AddSortedPoint(x2b, xAxisBoard, kResetValue);
-
+
AddSortedPoint(y1b, yAxisBoard, kResetValue);
AddSortedPoint(y2b, yAxisBoard, kResetValue);
}
- else if(displayType==kDisplayBoards)
+ else if(displayType==kDisplayBoards)
FillBins(inputObject, displayHisto, iBoard, -nLocations, x1b, x2b, y1b, y2b, kShiftEl, kShiftEl, displayOpt);
- else if(displayType==kDisplaySlats)
+ else if(displayType==kDisplaySlats)
FillBins(inputObject, displayHisto, slat, -1, x1b, x2b, y1b, y2b, kShiftEl, kShiftEl, displayOpt);
} // loop on local boards
-
+
if ( inputObject ) return kTRUE;
-
+
displayHisto->Reset();
-
+
// Book histos
const Float_t kMinDiff = 0.1;
-
- TArrayD* currArray[4] = {&xAxisStrip, &yAxisStrip,
- &xAxisBoard, &yAxisBoard};
+
+ TArrayD* currArray[4] = {&xAxisStrip, &yAxisStrip,
+ &xAxisBoard, &yAxisBoard};
for(Int_t iaxis=0; iaxis<4; iaxis++){
Int_t ipoint=0;
- while(TMath::Abs((*currArray[iaxis])[ipoint]-kResetValue)>kMinDiff) {
+ while(TMath::Abs((*currArray[iaxis])[ipoint]-kResetValue)>kMinDiff) {
ipoint++;
}
- if(ipoint>currArray[iaxis]->GetSize()-2)
+ if(ipoint>currArray[iaxis]->GetSize()-2)
AliWarning(Form("Array size (%i) lower than the number of points!", currArray[iaxis]->GetSize()));
currArray[iaxis]->Set(ipoint);
}
-
+
switch(displayType){
- case kDisplayStrips:
- displayHisto->SetBins(xAxisStrip.GetSize()-1, xAxisStrip.GetArray(),
- yAxisStrip.GetSize()-1, yAxisStrip.GetArray());
- break;
- case kDisplayBoards:
- case kDisplaySlats:
- displayHisto->SetBins(xAxisBoard.GetSize()-1, xAxisBoard.GetArray(),
- yAxisBoard.GetSize()-1, yAxisBoard.GetArray());
- break;
+ case kDisplayStrips:
+ displayHisto->SetBins(xAxisStrip.GetSize()-1, xAxisStrip.GetArray(),
+ yAxisStrip.GetSize()-1, yAxisStrip.GetArray());
+ break;
+ case kDisplayBoards:
+ case kDisplaySlats:
+ displayHisto->SetBins(xAxisBoard.GetSize()-1, xAxisBoard.GetArray(),
+ yAxisBoard.GetSize()-1, yAxisBoard.GetArray());
+ break;
}
-
+
displayHisto->SetName(displayHistoName.Data());
displayHisto->SetTitle(displayHistoTitle.Data());
displayHisto->SetXTitle("X (cm)");
displayHisto->SetYTitle("Y (cm)");
//displayHisto->SetStats(kFALSE);
-
+
return kTRUE;
}
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff