2 // $MpId: testIndicesLimits.C,v 1.4 2005/09/26 16:05:25 ivana Exp $
4 // Test macro for indices limits.
6 void testIndicesLimits(AliMpStationType station = kStation1,
7 AliMpPlaneType plane = kBendingPlane,
8 Bool_t rootInput = false)
10 AliMpSector *sector = 0;
12 AliMpSectorReader r(station, plane);
13 sector=r.BuildSector();
16 TString filePath = AliMpFiles::Instance()->SectorFilePath(station,plane);
17 filePath.ReplaceAll("zones.dat", "sector.root");
19 TFile f(filePath.Data(), "READ");
20 sector = (AliMpSector*)f.Get("Sector");
24 for (Int_t i=0; i<sector->GetNofRows(); i++) {
25 AliMpRow* row = sector->GetRow(i);
28 << row->GetLowIndicesLimit() << " "
29 << row->GetHighIndicesLimit() << endl;
31 // Loop over row segments
32 for (Int_t j=0; j<row->GetNofRowSegments(); j++) {
33 AliMpVRowSegment* rowSegment = row->GetRowSegment(j);
36 << "th row segment limits: "
37 << rowSegment->GetLowIndicesLimit() << " "
38 << rowSegment->GetHighIndicesLimit() << endl;
40 // Loop over motif positions
41 for (Int_t k=0; k<rowSegment->GetNofMotifs(); k++) {
42 Int_t mposID = rowSegment->GetMotifPositionId(k);
43 AliMpMotifPosition* mPos =
44 sector->GetMotifMap()->FindMotifPosition(mposID);
48 << " motif position limits: "
49 << mPos->GetLowIndicesLimit() << " "
50 << mPos->GetHighIndicesLimit() << endl;
53 cerr << "Motif position "
54 << mposID << " not found in the map" << endl;