HLT/TPCLib/test/testAliHLTTPCMapping.C

   1 // $Id$
   2
   3 /**************************************************************************
   4  * This file is property of and copyright by the ALICE HLT Project        *
   5  * ALICE Experiment at CERN, All rights reserved.                         *
   6  *                                                                        *
   7  * Primary Authors: Matthias Richter <Matthias.Richter@ift.uib.no>        *
   8  *                  for The ALICE HLT Project.                            *
   9  *                                                                        *
  10  * Permission to use, copy, modify and distribute this software and its   *
  11  * documentation strictly for non-commercial purposes is hereby granted   *
  12  * without fee, provided that the above copyright notice appears in all   *
  13  * copies and that both the copyright notice and this permission notice   *
  14  * appear in the supporting documentation. The authors make no claims     *
  15  * about the suitability of this software for any purpose. It is          *
  16  * provided "as is" without express or implied warranty.                  *
  17  **************************************************************************/
  18
  19 /** @file   testAliHLTTPCMapping.C
  20     @author Matthias Richter
  21     @date
  22     @brief  Test macro/program for the AliHLTTPCMapping class
  23  */
  24
  25 #ifndef __CINT__
  26 #include "TSystem.h"
  27 #include "AliHLTTPCMapping.h"
  28 #include "AliHLTTPCTransform.h"
  29 #include <ostream>
  30 #include <istream>
  31 #endif //__CINT__
  32
  33 /////////////////////////////////////////////////////////////////
  34 /////////////////////////////////////////////////////////////////
  35 /////////////////////////////////////////////////////////////////
  36 //
  37 // configuration of the test program
  38 //
  39
  40 // printouts or not
  41 const bool bVerbose=true;
  42
  43 /////////////////////////////////////////////////////////////////
  44 /////////////////////////////////////////////////////////////////
  45 /////////////////////////////////////////////////////////////////
  46
  47 class AliHLTTPCMapping;
  48 bool compareMapping(int patch, AliHLTTPCMapping* mapper);
  49
  50 int testAliHLTTPCMapping()
  51 {
  52   int iResult=0;
  53 #ifdef __CINT__
  54   gSystem->Load("libAliHLTUtil.so");
  55   gSystem->Load("libAliHLTRCU.so");
  56   gSystem->Load("libAliHLTTPC.so");
  57 #endif
  58   //AliHLTSystem gHLT;
  59
  60   const int nofMappers=6;
  61   AliHLTTPCMapping* mappers[nofMappers];
  62   AliHLTTPCMapping* mappers2[nofMappers];
  63   for (int i=0; i<nofMappers; i++) {
  64     mappers[i]=NULL;
  65     mappers2[i]=NULL;
  66   }
  67
  68   // create mappers
  69   for (int i=0; i<nofMappers; i++) {
  70     mappers[i]=new AliHLTTPCMapping(i);
  71   }
  72
  73   // check mappers
  74   for (int i=0; i<nofMappers; i++) {
  75     if (!compareMapping(i, mappers[i])) {
  76       iResult=-1;
  77     }
  78   }
  79   cout << "checking: 1st instance:";
  80   if (iResult<0)
  81     cout << " failed" << endl;
  82   else {
  83     cout << " ok" << endl;
  84
  85     // create 2nd instance
  86     for (int i=0; i<nofMappers; i++) {
  87       mappers2[i]=new AliHLTTPCMapping(i);
  88     }
  89
  90     // check 2nd instance
  91     for (int i=0; i<nofMappers; i++) {
  92       if (!compareMapping(i, mappers2[i])) {
  93         iResult=-1;
  94       }
  95     }
  96     cout << "checking: 2nd instance:";
  97     if (iResult<0)
  98       cout << " failed" << endl;
  99     else
 100       cout << " ok" << endl;
 101   }
 102
 103   // delete mappers
 104   for (int i=0; i<nofMappers; i++) {
 105     if (mappers[i]) delete mappers[i];
 106     if (mappers2[i]) delete mappers2[i];
 107   }
 108   return iResult;
 109 }
 110
 111 bool compareMapping(int patch, AliHLTTPCMapping* mapper)
 112 {
 113   bool result=true;
 114   if (!mapper) return false;
 115   ifstream inFile;
 116   TString filename;
 117   const char* basePath=getenv("ALICE_ROOT");
 118   if (basePath) {
 119     filename.Form("%s/TPC/mapping/Patch%d.data", basePath,patch);
 120   }
 121   inFile.open(filename.Data());
 122   if (!inFile) {
 123     cout << "Unable to open file: " << filename << endl;
 124     return false;
 125   }
 126
 127   UInt_t nHWAdd=0;
 128   UInt_t maxHWAdd=0;
 129   UInt_t hwAdd=0;
 130   UInt_t row=0;
 131   UInt_t pad=0;
 132   Int_t dummy=0;
 133   Int_t rowOffset=0;
 134   result=AliHLTTPCTransform::Slice2Sector(0, AliHLTTPCTransform::GetFirstRow(patch), dummy, rowOffset);
 135
 136   if(inFile >> nHWAdd && inFile >> maxHWAdd) {
 137     while(result && inFile>>hwAdd && inFile>>row && inFile>>pad){
 138       row-=rowOffset;
 139       if (row!=mapper->GetRow(hwAdd) || pad!=mapper->GetPad(hwAdd)) {
 140         cout << "mismatch at channel " << hwAdd << ": expected " << row << "/" << pad << "  got " << mapper->GetRow(hwAdd) << "/" << mapper->GetPad(hwAdd) << endl;
 141         result=false;
 142         break;
 143       }
 144     }
 145   }
 146   inFile.close();
 147   return result;
 148 }
 149
 150 int main(int /*argc*/, const char** /*argv*/)
 151 {
 152   return testAliHLTTPCMapping();
 153 }