Added the possibility to do diagnostics histograms.

[u/mrichter/AliRoot.git] / MUON / mapping / macros / testMotifTypeIterators.C

diff --git a/MUON/mapping/macros/testMotifTypeIterators.C b/MUON/mapping/macros/testMotifTypeIterators.C
index 7ff7167f1d40fd47a570c7775d96a3baaccbbcb2..2c2d998d1e54882fd5f5bd6212730f78343f6ae0 100644 (file)
--- a/MUON/mapping/macros/testMotifTypeIterators.C
+++ b/MUON/mapping/macros/testMotifTypeIterators.C
@@ -1,17 +1,38 @@
 // $Id$
+// $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="ABCDEFGHI";
-  //TString names="FEG";
+  TString names;
+  TString names2;
+  Int_t nv =0;
+  if ( station == kStation1 )
+    if ( plane == kBendingPlane ) 
+      names ="ABCDEFGHI";
+    else
+      names = "ABCDEFGHIJKLMN";
+  else if ( station == kStation2 ) 
+    if ( plane == kBendingPlane ) {
+      names ="ABCDEFGHIJKLMNOPQRSTUVWXY";
+      names2 ="abcdefghimnptuvvvvv";
+      nv = 5;
+    }  
+    else {
+      names = "ABCEFGHIJKLMN";
+      names2 ="abcdefgijklmnopqrstuwvvvvv";
+      nv = 5;
+    }  
+  Int_t nofMotifs = names.Length() + names2.Length(); 
+  // cout << " nofMotifs: " << nofMotifs << endl;   
     
-  TH2C* histos[] = new TH2C* [names.Length()];
-  TCanvas* canv[] = new TCanvas* [1+(names.Length()-1)/4];
+  TH2C* histos[] = new TH2C* [nofMotifs];
+  TCanvas* canv[] = new TCanvas* [1+(nofMotifs-1)/4];
   Int_t i;
-  for (i=0;i<1+(names.Length()-1)/4;++i){
+  for (i=0;i<1+(nofMotifs-1)/4;++i){
     // canv[i] = new TCanvas(Form("canv%d",i),"Iterator viewing...");
                // CINT limitation on DEC
               
@@ -21,11 +42,27 @@ void testMotifTypeIterators(AliMpStationType station = kStation1,
     canv[i]->Divide(2,2); 
   }
     
-  AliMpReader r(station, plane);
+  AliMpMotifReader r(station, plane);
   //r.SetVerboseLevel(2);
 
-  for (i=0;i<names.Length();++i){
-    AliMpMotifType *mt = r.BuildMotifType(names[i]);
+  for (i=0;i<nofMotifs;++i){
+
+    // Get motif name
+    TString mname;
+    if (i<names.Length())
+      mname = names(i, 1);
+    else {
+      mname = names2(i-names.Length(), 1); 
+      if (mname == "v")
+        mname += i - names.Length() - (names2.Length()-nv-1);
+      else        
+        mname += "1";
+    }  
+    //if (i==36) continue;  
+        // break for these motifs (St2, BP) - to be investigated
+   
+    AliMpMotifType *mt = r.BuildMotifType(mname);
+
     canv[i/4]->cd(1+ (i%4));
     //histos[i] = new TH2C(Form("h%d",i),Form("Motif type %c",names[i]),
     //                     mt->GetNofPadsX(),-0.5,mt->GetNofPadsX()-0.5,
@@ -33,7 +70,6 @@ void testMotifTypeIterators(AliMpStationType station = kStation1,
                // CINT limitation on DEC
 
     TString hname("h"); hname += i;
-    TString mname = names(i,1);               
 
     histos[i] = new TH2C(hname.Data(), mname.Data(),
                          mt->GetNofPadsX(),-0.5,mt->GetNofPadsX()-0.5,
@@ -52,10 +88,8 @@ void testMotifTypeIterators(AliMpStationType station = kStation1,
                       it.CurrentItem().GetIndices().GetSecond(),num);
     }
 
-    delete mt;
-
+    //delete mt;
     histos[i]->Draw("text");
     canv[i/4]->Update();
-    
   }
 }