TPC/macros/testTPC/rec.C

   1 void rec(const char *filename="data.root")\r
   2 {\r
   3   gSystem->Load("libXrdClient.so");
   4   gSystem->Load("libNetX.so");
   5   //\r
   6   // First version of the reconstruction\r
   7   // script for the FDR'07\r
   8 \r
   9   // Set the CDB storage location\r
  10   // AliLog::SetModuleDebugLevel("STEER",2);\r
  11   AliCDBManager * man = AliCDBManager::Instance();\r
  12   //man->SetDefaultStorage("alien://folder=/alice/data/2007/LHC07w/OCDB/");\r
  13   man->SetDefaultStorage("/data/test2007/OCDB");\r
  14   man->SetSpecificStorage("GRP/GRP/Data","local://$ALICE_ROOT");\r
  15   man->SetSpecificStorage("ITS/Calib/DDLMapSDD","local://$ALICE_ROOT");\r
  16   man->SetSpecificStorage("MUON/Calib/Mapping","local://$ALICE_ROOT");\r
  17   man->SetSpecificStorage("MUON/Calib/DDLStore","local://$ALICE_ROOT");\r
  18 \r
  19   // TPC settings\r
  20   AliLog::SetClassDebugLevel("AliTPCclustererMI",2);\r
  21   AliTPCRecoParam * tpcRecoParam = AliTPCRecoParam::GetCosmicTestParam(kTRUE);\r
  22   tpcRecoParam->SetTimeInterval(60,940);\r
  23   tpcRecoParam->Dump();\r
  24   AliTPCReconstructor::SetRecoParam(tpcRecoParam);\r
  25   AliTPCReconstructor::SetStreamLevel(1);\r
  26 \r
  27   // PHOS settings\r
  28   AliPHOSRecoParam* recEmc = new AliPHOSRecoParamEmc();\r
  29   recEmc->SetSubtractPedestals(kTRUE);\r
  30   recEmc->SetMinE(0.05);\r
  31   recEmc->SetClusteringThreshold(0.10);\r
  32   AliPHOSReconstructor::SetRecoParamEmc(recEmc);\r
  33 \r
  34   // T0 settings\r
  35   AliLog::SetModuleDebugLevel("T0", 10);\r
  36 \r
  37   // MUON settings\r
  38   AliLog::SetClassDebugLevel("AliMUONRawStreamTracker",3);\r
  39  \r
  40   // Tracking settings\r
  41   //  AliMagFMaps* field = new AliMagFMaps("Maps","Maps", 2, 1., 10., 1);\r
  42   AliMagFMaps* field = new AliMagFMaps("Maps","Maps", 2, 0., 10., 2);\r
  43   AliTracker::SetFieldMap(field,1);\r
  44 \r
  45   // AliReconstruction settings\r
  46   AliReconstruction rec;\r
  47   rec.SetUniformFieldTracking(kFALSE);\r
  48   rec.SetWriteESDfriend(kTRUE);\r
  49   rec.SetWriteAlignmentData();\r
  50   rec.SetInput(filename);\r
  51   //
  52   //rec.SetRunLocalReconstruction("");\r
  53   rec.SetRunReconstruction("TPC");\r
  54   rec.SetFillESD("TPC");
  55   rec.SetRunVertexFinder(kFALSE);
  56   rec.SetRunQA(kFALSE);
  57   //
  58   //rec.SetEventRange(1,5);\r
  59 \r
  60   // In case some detectors have to be switched off...\r
  61   //  rec.SetRunLocalReconstruction("ALL");\r
  62   //  rec.SetRunTracking("ALL");\r
  63   //  rec.SetFillESD("ALL");\r
  64   // Disable vertex finder for the moment\r
  65   //  rec.SetRunVertexFinder(kFALSE);\r
  66 \r
  67   // To be enabled if some equipment IDs are not set correctly by DAQ\r
  68   //  rec.SetEquipmentIdMap("EquipmentIdMap.data");\r
  69 \r
  70   // Detector options if any\r
  71   rec.SetOption("MUON","SAVEDIGITS");\r
  72   rec.SetOption("TPC","OldRCUFormat");\r
  73   rec.SetOption("PHOS","OldRCUFormat");\r
  74 \r
  75   // To be enabled when CTP readout starts\r
  76   rec.SetFillTriggerESD(kFALSE);\r
  77 \r
  78   // all events in one single file\r
  79   rec.SetNumberOfEventsPerFile(-1);\r
  80 \r
  81   // switch off cleanESD\r
  82   rec.SetCleanESD(kFALSE);\r
  83 \r
  84   //AliLog::SetGlobalDebugLevel(2);\r
  85   rec.Run();\r
  86 \r
  87   cout << "-----------------------------------------------------------------" << endl;\r
  88   cout << "-----------------------------------------------------------------" << endl;\r
  89   cout << "--------- Reconstruction Completed. Start merging QAs -----------" << endl;\r
  90   cout << "-----------------------------------------------------------------" << endl;\r
  91   cout << "-----------------------------------------------------------------" << endl;\r
  92   AliQADataMakerSteer qas;\r
  93   qas.Merge();\r
  94 }\r