Updated version of raw data code (D.Favretto)

[u/mrichter/AliRoot.git] / TPC / AliTPCH5OptimizedTables.C

diff --git a/TPC/AliTPCH5OptimizedTables.C b/TPC/AliTPCH5OptimizedTables.C
index d53420caabf0c8b9011daaa381d8be976addde48..dfb6b33b821c8a64a3a1e121398bbcd9bef6b015 100644 (file)
--- a/TPC/AliTPCH5OptimizedTables.C
+++ b/TPC/AliTPCH5OptimizedTables.C
@@ -1,5 +1,5 @@
 #if !defined(__CINT__) || defined(__MAKECINT__)
-#include <fstream.h>
+#include <Riotream.h>
 #include <alles.h>
 #include "AliTPCCompression.h"
 #endif
@@ -14,24 +14,65 @@ void AliTPCH5OptimizedTables(const char* fSource="AltroFormat.dat",const char* f
   AliTPCCompression *util = new AliTPCCompression();
   TStopwatch timer;
   //verbose level can be: 0=silent 1=few messages 2=pedantic output
-  util->SetVerbose(2);
-  //Tables are created
-  util->CreateTables(fSource,NumTable);
-  //util->ReadAltroFormat("File1.txt","AltroFormat.dat");
-  //The source file is compressed 
+  util->SetVerbose(1);
+  Int_t choice;
+  do{
+    cout<<"**** Chose the tables set **** "<<endl;
+    cout<<"1==> Create tables from the input file "<<endl;
+    cout<<"2==> Use external optimized tables (txt format)"<<endl;
+    cout<<"3==> Time gap and Bunch length tables generated using formulas "<<endl;
+    cout<<"Insert the corresponding number: ";
+    cin>>choice;
+    cout<<endl;
+  }while((choice<1)||(choice>3));
+  switch(choice){
+  case 1:{
+    //Tables are created
+    util->CreateTables(fSource,NumTable);
+    cout<<"Tables have been created"<<endl;
+    break;
+  }
+  case 2:{
+    util->CreateTablesFromTxtFiles(NumTable);
+    break;
+  }
+  case 3:{
+    Double_t beta,gamma=0;
+    ULong_t  mul=0;    
+    cout<<"Multiplicity (suggested 20000) ==> ";
+    cin>>mul;
+    cout<<"Gamma (suggested 4.77) ==> ";
+    cin>>gamma;
+    cout<<"Beta (suggested 0.37) ==> ";
+    cin>>beta;
+    util->CreateTables(fSource,NumTable);
+    util->CreateTableFormula(gamma,mul,300,0);
+    util->CreateTableFormula(beta,mul,445,1);
+    break;
+  }
+  };
+
+  //BE CAREFUL, the following method must be used only for debugging and
+  //it is highly suggested to use it only for debugging files
+  //reasonably small, because otherwise the size of the txt files can reach
+  //quickly several MB wasting time and disk space.
+
+  //  util->ReadAltroFormat("File1.txt",fSource);
   
+  //The source file is compressed 
   timer.Start();
   util->CompressDataOptTables(NumTable,fSource,fDest);
   timer.Stop();
   timer.Print();
   
-  /*  
+
+  /*
   //The Compressed file is decompressed  
   timer.Start();
   util->DecompressDataOptTables(NumTable,fDest);
   timer.Stop();
   timer.Print();
-  //util->ReadAltroFormat("File2.txt","SourceDecompressed.dat");
+  util->ReadAltroFormat("File2.txt","SourceDecompressed.dat");
   */
   delete util;
 }