Added Bool_t rootInput argument; if set to true, the sector
authorivana <ivana@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 27 Sep 2005 09:29:09 +0000 (09:29 +0000)
committerivana <ivana@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 27 Sep 2005 09:29:09 +0000 (09:29 +0000)
is loaded from the Root file instead of the ASCI data

14 files changed:
MUON/mapping/macros/testAnyPadIterators.C
MUON/mapping/macros/testExistingPads.C
MUON/mapping/macros/testGraphics.C
MUON/mapping/macros/testIndicesLimits.C
MUON/mapping/macros/testMotifTypeIterators.C
MUON/mapping/macros/testNeighboursPadIterator.C
MUON/mapping/macros/testPadDimensions.C
MUON/mapping/macros/testPadsUpEtc.C
MUON/mapping/macros/testPrintLimits.C
MUON/mapping/macros/testReadMotifType.C
MUON/mapping/macros/testReadSector.C
MUON/mapping/macros/testSectorAreaIterator.C
MUON/mapping/macros/testSectorFind.C
MUON/mapping/macros/testSectorPadIterators.C

index efa667b757bc622058ee52ed431edfff4c65152a..e666cbb2e830fe873b94e8a2a969c59f67f25ff3 100644 (file)
@@ -1,5 +1,5 @@
 // $Id$
-// $MpId: testAnyPadIterators.C,v 1.12 2005/08/24 08:53:27 ivana Exp $
+// $MpId: testAnyPadIterators.C,v 1.13 2005/09/26 16:05:24 ivana Exp $
 //
 // Test macro for reading  sector, and iterate over it
 
@@ -25,19 +25,30 @@ void MarkPads(AliMpVPadIterator& it,Int_t xmax,Int_t ymax,Bool_t print=kTRUE)
 
 void testAnyPadIterators(AliMpStationType station = kStation1,
                          AliMpPlaneType plane = kBendingPlane, 
+                        Bool_t rootInput = false,
                          Int_t i=50, Int_t j=50)
 {
-  AliMpSectorReader r(station, plane);
-  AliMpSector* sect = r.BuildSector();
+  AliMpSector *sector = 0;
+  if (!rootInput) {
+    AliMpSectorReader r(station, plane);
+    sector=r.BuildSector();
+  }
+  else  {
+    TString filePath = AliMpFiles::Instance()->SectorFilePath(station,plane);
+    filePath.ReplaceAll("zones.dat", "sector.root"); 
+
+    TFile f(filePath.Data(), "READ");
+    sector = (AliMpSector*)f.Get("Sector");
+  }  
     
   TCanvas *canv = new TCanvas("canv");
   canv->Divide(2,2);
   //canv_1->Divide(2);
   
   canv->cd(1);
-  MarkPads(AliMpSectorPadIterator(sect), 150, 250, kFALSE);
+  MarkPads(AliMpSectorPadIterator(sector), 150, 250, kFALSE);
   canv->cd(2);
-  AliMpVMotif* motif = sect->FindMotif(TVector2(300,30));
+  AliMpVMotif* motif = sector->FindMotif(TVector2(300,30));
 
   if (motif) {
     AliMpMotifType* motifType = motif->GetMotifType();
@@ -47,16 +58,16 @@ void testAnyPadIterators(AliMpStationType station = kStation1,
   } else cout<<"No motif found at given position..."<<endl;
   
   canv->cd(3);
-  //MarkPads(*AliMpPadIteratorPtr(AliMpSectorSegmentation(sect)->CreateIterator(AliMpIntPair(i,j)))
-  AliMpSectorSegmentation segm(sect);
+  //MarkPads(*AliMpPadIteratorPtr(AliMpSectorSegmentation(sector)->CreateIterator(AliMpIntPair(i,j)))
+  AliMpSectorSegmentation segm(sector);
   AliMpPad pad = segm.PadByIndices(AliMpIntPair(i,j));
-  MarkPads(AliMpNeighboursPadIterator(&AliMpSectorSegmentation(sect),pad)
+  MarkPads(AliMpNeighboursPadIterator(&AliMpSectorSegmentation(sector),pad)
            ,i+8,j+8);
   cout<<"________________ Neighbours __________________________"<<endl;
   canv->cd(4);
   Int_t motifPosId = 20; 
   if (plane == kNonBendingPlane) motifPosId = 19;
-  AliMpMotifPosition* motifPos = sect->GetMotifMap()->FindMotifPosition(motifPosId);
+  AliMpMotifPosition* motifPos = sector->GetMotifMap()->FindMotifPosition(motifPosId);
   if (motifPos){
     //MarkPads(*AliMpPadIteratorPtr(motifPos->CreateIterator()),15,15);
     MarkPads(AliMpMotifPositionPadIterator(motifPos),15,15);
index a940f057d4a907488586b3557c061c12be176caf..f215f9ba7b5154fd29f7381cc4c20dabc7454bcd 100644 (file)
@@ -1,13 +1,25 @@
 // $Id$
-// $MpId: testExistingPads.C,v 1.10 2005/08/24 08:53:27 ivana Exp $
+// $MpId: testExistingPads.C,v 1.11 2005/09/26 16:05:25 ivana Exp $
 //
 // Test macro for testing which pad is seen as "existing" by AliMpSector.
 
 void testExistingPads(AliMpStationType station = kStation1,
-                      AliMpPlaneType plane = kBendingPlane) 
+                      AliMpPlaneType plane = kBendingPlane, 
+                     Bool_t rootInput = false)
 {
-  AliMpSectorReader r(station, plane);
-  AliMpSector *sector=r.BuildSector();
+  AliMpSector *sector = 0;
+  if (!rootInput) {
+    AliMpSectorReader r(station, plane);
+    sector=r.BuildSector();
+  }
+  else  {
+    TString filePath = AliMpFiles::Instance()->SectorFilePath(station,plane);
+    filePath.ReplaceAll("zones.dat", "sector.root"); 
+
+    TFile f(filePath.Data(), "READ");
+    sector = (AliMpSector*)f.Get("Sector");
+  }  
+
   AliMpSectorSegmentation segmentation(sector);
   AliMpVPainter* painter = AliMpVPainter::CreatePainter(sector);
 
index a57ab9e3ce83a8d4c33b0d03d225cb1a190a50bd..28e0a9c02e4b7cfa140b46acf6216a053516de80 100644 (file)
@@ -1,13 +1,25 @@
 // $Id$
-// $MpId: testGraphics.C,v 1.11 2005/08/24 08:53:27 ivana Exp $
+// $MpId: testGraphics.C,v 1.12 2005/09/26 16:05:25 ivana Exp $
 //
 // Test macro for drawing sector data.
 
 void testGraphics(AliMpStationType station = kStation1,
-                  AliMpPlaneType plane = kBendingPlane) 
+                  AliMpPlaneType plane = kBendingPlane,
+                 Bool_t rootInput = false) 
 {
-  AliMpSectorReader r(station, plane);
-  AliMpSector *sector=r.BuildSector();
+  AliMpSector *sector = 0;
+  if (!rootInput) {
+    AliMpSectorReader r(station, plane);
+    sector=r.BuildSector();
+  }
+  else  {
+    TString filePath = AliMpFiles::Instance()->SectorFilePath(station,plane);
+    filePath.ReplaceAll("zones.dat", "sector.root"); 
+
+    TFile f(filePath.Data(), "READ");
+    sector = (AliMpSector*)f.Get("Sector");
+  }  
+    
   AliMpVPainter *painter=AliMpVPainter::CreatePainter(sector);
 
   TCanvas* canvas = new TCanvas();
index e69425cc96dae7d2336fb97b48bc05c7fff29647..41c3fa80ba1b11c2cddb83f0e83b76c5cd5991c2 100644 (file)
@@ -1,16 +1,24 @@
 // $Id$
-// $MpId: testIndicesLimits.C,v 1.3 2005/08/24 08:53:27 ivana Exp $
+// $MpId: testIndicesLimits.C,v 1.4 2005/09/26 16:05:25 ivana Exp $
 //
 // Test macro for indices limits.
 
 void testIndicesLimits(AliMpStationType station = kStation1,
-                       AliMpPlaneType plane = kBendingPlane) 
+                       AliMpPlaneType plane = kBendingPlane, 
+                      Bool_t rootInput = false)
 {
-  AliMpSectorReader reader(station, plane);  
-  //reader.SetVerboseLevel(1);
-  
-  // Read data 
-  AliMpSector* sector = reader.BuildSector();
+  AliMpSector *sector = 0;
+  if (!rootInput) {
+    AliMpSectorReader r(station, plane);
+    sector=r.BuildSector();
+  }
+  else  {
+    TString filePath = AliMpFiles::Instance()->SectorFilePath(station,plane);
+    filePath.ReplaceAll("zones.dat", "sector.root"); 
+
+    TFile f(filePath.Data(), "READ");
+    sector = (AliMpSector*)f.Get("Sector");
+  }  
 
   // Loop over rows
   for (Int_t i=0; i<sector->GetNofRows(); i++) {
index 4480f0fdd55331b5a33b17afd2d32760d66d0606..2c2d998d1e54882fd5f5bd6212730f78343f6ae0 100644 (file)
@@ -1,10 +1,11 @@
 // $Id$
-// $MpId: testMotifTypeIterators.C,v 1.10 2005/08/24 08:53:27 ivana Exp $
+// $MpId: testMotifTypeIterators.C,v 1.11 2005/09/26 16:05:25 ivana Exp $
 //
 // Test macro for reading motif type data and iterate over them.
 
 void testMotifTypeIterators(AliMpStationType station = kStation1,
-                            AliMpPlaneType plane = kBendingPlane)
+                            AliMpPlaneType plane = kBendingPlane,
+                           Bool_t rootInput = false)
 {
   TString names;
   TString names2;
index 44cd2a5e3d22b59ee8d029b2a44b5ea12ef3c268..7a8146e42257dbf53ce96dd81756df0cadcc4580 100644 (file)
@@ -1,19 +1,26 @@
 // $Id$
-// $MpId: testNeighboursPadIterator.C,v 1.8 2005/08/24 08:53:27 ivana Exp $
+// $MpId: testNeighboursPadIterator.C,v 1.9 2005/09/26 16:05:25 ivana Exp $
 //
 // Test macro for reading  sector, and iterate over it
 
 void testNeighboursPadIterator(AliMpStationType station = kStation1,
                                AliMpPlaneType plane = kBendingPlane, 
+                              Bool_t rootInput = false,
                                Int_t i=50, Int_t j=50)
 {
-  if (!gInterpreter->IsLoaded("mlibs.C")){ 
-    gROOT->LoadMacro("mlibs.C");
-    gInterpreter->ProcessLine("mlibs()");
+  AliMpSector *sector = 0;
+  if (!rootInput) {
+    AliMpSectorReader r(station, plane);
+    sector=r.BuildSector();
+  }
+  else  {
+    TString filePath = AliMpFiles::Instance()->SectorFilePath(station,plane);
+    filePath.ReplaceAll("zones.dat", "sector.root"); 
+
+    TFile f(filePath.Data(), "READ");
+    sector = (AliMpSector*)f.Get("Sector");
   }  
 
-  AliMpSectorReader r(station, plane);
-  AliMpSector* sect = r.BuildSector();
   AliMpSectorSegmentation segm(sect);  
   
   TCanvas *can = new TCanvas("canv");
index 1c7dba9691cc2c9c35f7628adb99290e4ccae99f..31638d34e41e2abd854fe82271334c217dc22d7d 100644 (file)
@@ -1,14 +1,26 @@
 // $Id$
-// $MpId: testPadDimensions.C,v 1.4 2005/08/24 08:53:27 ivana Exp $
+// $MpId: testPadDimensions.C,v 1.5 2005/09/26 16:05:25 ivana Exp $
 //
 // Test macro for testing retrieving of pad dimensions from
 // the map in AliMpSectorSegmentation.
 
 void testPadDimensions(AliMpStationType station = kStation1,
-                       AliMpPlaneType plane = kBendingPlane) 
+                       AliMpPlaneType plane = kBendingPlane,
+                      Bool_t rootInput = false)
 {
-  AliMpSectorReader r(station, plane);
-  AliMpSector* sector=r.BuildSector();
+  AliMpSector *sector = 0;
+  if (!rootInput) {
+    AliMpSectorReader r(station, plane);
+    sector=r.BuildSector();
+  }
+  else  {
+    TString filePath = AliMpFiles::Instance()->SectorFilePath(station,plane);
+    filePath.ReplaceAll("zones.dat", "sector.root"); 
+
+    TFile f(filePath.Data(), "READ");
+    sector = (AliMpSector*)f.Get("Sector");
+  }  
+
   AliMpSectorSegmentation segmentation(sector);  
   segmentation.PrintZones(); 
   
index ac19b54340f2e0e42acc0396c373f599ec85379a..b8b071e2652257d44b4446e98fa9e218caaaaac2 100644 (file)
@@ -1,19 +1,31 @@
 // $Id$
-// $MpId: testPadsUpEtc.C,v 1.6 2005/08/24 08:53:27 ivana Exp $
+// $MpId: testPadsUpEtc.C,v 1.7 2005/09/26 16:05:25 ivana Exp $
 //
 // Test macro that starts from a given pad and prints 
 // all pads up, down, right, left from this pad
 // (up to the plane border).
 
 void testPadsUpEtc(AliMpStationType station = kStation1,
-                   AliMpPlaneType  planeType = kBendingPlane)
+                   AliMpPlaneType  plane = kBendingPlane,
+                  Bool_t rootInput = false)
 {
-  AliMpSectorReader reader(station, planeType);  
-  AliMpSector* sector = reader.BuildSector();
+  AliMpSector *sector = 0;
+  if (!rootInput) {
+    AliMpSectorReader r(station, plane);
+    sector=r.BuildSector();
+  }
+  else  {
+    TString filePath = AliMpFiles::Instance()->SectorFilePath(station,plane);
+    filePath.ReplaceAll("zones.dat", "sector.root"); 
+
+    TFile f(filePath.Data(), "READ");
+    sector = (AliMpSector*)f.Get("Sector");
+  }  
+    
   AliMpSectorSegmentation segmentation(sector);
   
   AliMpIntPair indices(85, 101);
-  if( planeType == kNonBendingPlane) indices = AliMpIntPair(129, 10);
+  if( plane == kNonBendingPlane) indices = AliMpIntPair(129, 10);
  
   AliMpPad pad;
   if (segmentation.HasPad(indices)) {
index 50f88f201fab5f26070bde022c73aed7afe87f28..f14163d228095be2c5df874366fd193639bce434 100644 (file)
@@ -1,15 +1,26 @@
 // $Id$
-// $MpId: testPrintLimits.C,v 1.7 2005/08/24 08:53:27 ivana Exp $
+// $MpId: testPrintLimits.C,v 1.8 2005/09/26 16:05:25 ivana Exp $
 //
 // Test macro for making an output file, where all mapping elements
 // indices & positions are written.
 
 void testPrintLimits(AliMpStationType station = kStation1,
-                    AliMpPlaneType plane = kBendingPlane, ostream& out=cout)
+                    AliMpPlaneType plane = kBendingPlane, 
+                   Bool_t rootInput = false,
+                   ostream& out=cout)
 {
-  AliMpSectorReader r(station, plane);
+  AliMpSector *sector = 0;
+  if (!rootInput) {
+    AliMpSectorReader r(station, plane);
+    sector=r.BuildSector();
+  }
+  else  {
+    TString filePath = AliMpFiles::Instance()->SectorFilePath(station,plane);
+    filePath.ReplaceAll("zones.dat", "sector.root"); 
 
-  AliMpSector *sector=r.BuildSector();
+    TFile f(filePath.Data(), "READ");
+    sector = (AliMpSector*)f.Get("Sector");
+  }  
 
   AliMpVPainter* painter = AliMpVPainter::CreatePainter(sector);
   painter->Draw("ZSSMP");
index b8e9261b602f13925463204bf39516e8d1ef9626..bdf743f28e72e10f512911f4971f2ca228a5d0fa 100644 (file)
@@ -1,10 +1,11 @@
 // $Id$
-// $MpId: testReadMotifType.C,v 1.11 2005/08/24 08:53:27 ivana Exp $
+// $MpId: testReadMotifType.C,v 1.12 2005/09/26 16:05:25 ivana Exp $
 //
 // Test macro for reading motif type data.
 
 void testReadMotifType(AliMpStationType station = kStation1,
-                       AliMpPlaneType plane = kBendingPlane)
+                       AliMpPlaneType plane = kBendingPlane,
+                      Bool_t rootInput = false)
 {
   AliMpMotifReader r(station, plane);
   //r.SetVerboseLevel(2);
index c186d84359b7d5e58a429c237e397fc060096551..35059c7c3bc4db3f13da611715b171611b1fb855 100644 (file)
@@ -1,19 +1,36 @@
 // $Id$
-// $MpId: testReadSector.C,v 1.13 2005/09/02 10:57:47 ivana Exp $
+// $MpId: testReadSector.C,v 1.14 2005/09/26 16:04:07 ivana Exp $
 //
 // Test macro for reading sector data.
 
 #include <iomanip>
 
 void testReadSector(AliMpStationType station = kStation1,
-                    AliMpPlaneType plane = kBendingPlane) 
+                    AliMpPlaneType plane = kBendingPlane, 
+                   Bool_t rootInput = false)
 {
-  AliMpSectorReader reader(station, plane);  
-  //reader.SetVerboseLevel(1);
-  
-  // Read data 
-  AliMpSector* sector = reader.BuildSector();
+  AliMpSector *sector = 0;
+  if (!rootInput) {
+    AliMpSectorReader r(station, plane);
+    //reader.SetVerboseLevel(1);
+    sector=r.BuildSector();
+
+    // Write sector on Root file
+    TString filePath = AliMpFiles::Instance()->SectorFilePath(station,plane);
+    filePath.ReplaceAll("zones.dat", "sector.root"); 
   
+    TFile f(filePath.Data(), "RECREATE");
+    sector->Write();
+    f.Close();
+  }
+  else  {
+    TString filePath = AliMpFiles::Instance()->SectorFilePath(station,plane);
+    filePath.ReplaceAll("zones.dat", "sector.root"); 
+
+    TFile f(filePath.Data(), "READ");
+    sector = (AliMpSector*)f.Get("Sector");
+  }  
+
   cout << endl;
 
   // Sector geometry
index f7f593b764d7cd269000deccb75be4c555f3d39a..e2176ff56d73e1149ee5ab4d9c41f8671714f068 100644 (file)
@@ -1,5 +1,5 @@
 // $Id$
-// $MpId: testSectorAreaIterator.C,v 1.3 2005/08/24 08:53:27 ivana Exp $
+// $MpId: testSectorAreaIterator.C,v 1.4 2005/09/26 16:05:25 ivana Exp $
 //
 // Test macro for iterating over the whole sector
 
@@ -48,10 +48,22 @@ void MarkPads(AliMpVPadIterator& it, Double_t xmax, Double_t ymax,
 }
 
 void testSectorAreaIterator(AliMpStationType station = kStation1,
-                            AliMpPlaneType planeType = kBendingPlane)
+                            AliMpPlaneType plane = kBendingPlane,
+                           Bool_t rootInput = false)
 {
-  AliMpSectorReader reader(station, planeType);  
-  AliMpSector* sector = reader.BuildSector();
+  AliMpSector *sector = 0;
+  if (!rootInput) {
+    AliMpSectorReader r(station, plane);
+    sector=r.BuildSector();
+  }
+  else  {
+    TString filePath = AliMpFiles::Instance()->SectorFilePath(station,plane);
+    filePath.ReplaceAll("zones.dat", "sector.root"); 
+
+    TFile f(filePath.Data(), "READ");
+    sector = (AliMpSector*)f.Get("Sector");
+  }  
+
   AliMpSectorSegmentation segmentation(sector);
 
   AliMpArea area;
index 85717627441c1cf0317513d64bc2b271b81de1ac..785685a5a3e1e957065d03d0ce01172b74daf436 100644 (file)
@@ -1,15 +1,26 @@
 // $Id$
-// $MpId: testSectorFind.C,v 1.8 2005/08/24 08:53:27 ivana Exp $
+// $MpId: testSectorFind.C,v 1.9 2005/09/26 16:05:25 ivana Exp $
 //
 // Test macro for which verify that all FindPosition, FindIndices
 // and FindLocation methods are consistents between them.
 
 void testSectorFind(AliMpStationType station = kStation1,
-                    AliMpPlaneType plane = kBendingPlane)
+                    AliMpPlaneType plane = kBendingPlane,
+                   Bool_t rootInput = false) 
 {
-  AliMpSectorReader r(station, plane);
+  AliMpSector *sector = 0;
+  if (!rootInput) {
+    AliMpSectorReader r(station, plane);
+    sector=r.BuildSector();
+  }
+  else  {
+    TString filePath = AliMpFiles::Instance()->SectorFilePath(station,plane);
+    filePath.ReplaceAll("zones.dat", "sector.root"); 
+
+    TFile f(filePath.Data(), "READ");
+    sector = (AliMpSector*)f.Get("Sector");
+  }  
 
-  AliMpSector *sector=r.BuildSector();
   AliMpSectorSegmentation segmentation(sector);
   AliMpVPainter* painter = AliMpVPainter::CreatePainter(sector);
   painter->Draw("ZSSMP");
index a0422754d5572618a265455c28687cb3036dbe15..feea8b63554e5a39c3d4aa2f6aad2fd42da2a89c 100644 (file)
@@ -1,13 +1,24 @@
 // $Id$
-// $MpId: testSectorPadIterators.C,v 1.7 2005/08/24 08:53:27 ivana Exp $
+// $MpId: testSectorPadIterators.C,v 1.8 2005/09/26 16:05:25 ivana Exp $
 //
 // Test macro for reading  sector, and iterate over it
 
 void testSectorPadIterators(AliMpStationType station = kStation1,
-                            AliMpPlaneType plane = kBendingPlane)
+                            AliMpPlaneType plane = kBendingPlane,
+                           Bool_t rootInput = false)
 {
-  AliMpSectorReader r(station, plane);
-  AliMpSector* sect = r.BuildSector();
+  AliMpSector *sector = 0;
+  if (!rootInput) {
+    AliMpSectorReader r(station, plane);
+    sector=r.BuildSector();
+  }
+  else  {
+    TString filePath = AliMpFiles::Instance()->SectorFilePath(station,plane);
+    filePath.ReplaceAll("zones.dat", "sector.root"); 
+
+    TFile f(filePath.Data(), "READ");
+    sector = (AliMpSector*)f.Get("Sector");
+  }  
   
   Int_t num=0;
   
@@ -16,7 +27,7 @@ void testSectorPadIterators(AliMpStationType station = kStation1,
   const Double_t xmax=150;
   const Double_t ymax=250;
 
-  AliMpSectorPadIterator it = AliMpSectorPadIterator(sect);
+  AliMpSectorPadIterator it = AliMpSectorPadIterator(sector);
 
   for (it.First(); ! it.IsDone(); it.Next()) {
     AliMpIntPair indices = it.CurrentItem().GetIndices();
@@ -27,5 +38,5 @@ void testSectorPadIterators(AliMpStationType station = kStation1,
     marker->Draw();
   }
   
-  delete sect;
+  delete sector;
 }