Adding rec.C for test TPC data
authormarian <marian@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 8 Feb 2008 17:29:32 +0000 (17:29 +0000)
committermarian <marian@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 8 Feb 2008 17:29:32 +0000 (17:29 +0000)
TPC/macros/testTPC/rec.C [new file with mode: 0644]

diff --git a/TPC/macros/testTPC/rec.C b/TPC/macros/testTPC/rec.C
new file mode 100644 (file)
index 0000000..5784784
--- /dev/null
@@ -0,0 +1,94 @@
+void rec(const char *filename="data.root")\r
+{\r
+  gSystem->Load("libXrdClient.so");
+  gSystem->Load("libNetX.so");
+  //\r
+  // First version of the reconstruction\r
+  // script for the FDR'07\r
+\r
+  // Set the CDB storage location\r
+  // AliLog::SetModuleDebugLevel("STEER",2);\r
+  AliCDBManager * man = AliCDBManager::Instance();\r
+  //man->SetDefaultStorage("alien://folder=/alice/data/2007/LHC07w/OCDB/");\r
+  man->SetDefaultStorage("/data/test2007/OCDB");\r
+  man->SetSpecificStorage("GRP/GRP/Data","local://$ALICE_ROOT");\r
+  man->SetSpecificStorage("ITS/Calib/DDLMapSDD","local://$ALICE_ROOT");\r
+  man->SetSpecificStorage("MUON/Calib/Mapping","local://$ALICE_ROOT");\r
+  man->SetSpecificStorage("MUON/Calib/DDLStore","local://$ALICE_ROOT");\r
+\r
+  // TPC settings\r
+  AliLog::SetClassDebugLevel("AliTPCclustererMI",2);\r
+  AliTPCRecoParam * tpcRecoParam = AliTPCRecoParam::GetCosmicTestParam(kTRUE);\r
+  tpcRecoParam->SetTimeInterval(60,940);\r
+  tpcRecoParam->Dump();\r
+  AliTPCReconstructor::SetRecoParam(tpcRecoParam);\r
+  AliTPCReconstructor::SetStreamLevel(1);\r
+\r
+  // PHOS settings\r
+  AliPHOSRecoParam* recEmc = new AliPHOSRecoParamEmc();\r
+  recEmc->SetSubtractPedestals(kTRUE);\r
+  recEmc->SetMinE(0.05);\r
+  recEmc->SetClusteringThreshold(0.10);\r
+  AliPHOSReconstructor::SetRecoParamEmc(recEmc);\r
+\r
+  // T0 settings\r
+  AliLog::SetModuleDebugLevel("T0", 10);\r
+\r
+  // MUON settings\r
+  AliLog::SetClassDebugLevel("AliMUONRawStreamTracker",3);\r
\r
+  // Tracking settings\r
+  //  AliMagFMaps* field = new AliMagFMaps("Maps","Maps", 2, 1., 10., 1);\r
+  AliMagFMaps* field = new AliMagFMaps("Maps","Maps", 2, 0., 10., 2);\r
+  AliTracker::SetFieldMap(field,1);\r
+\r
+  // AliReconstruction settings\r
+  AliReconstruction rec;\r
+  rec.SetUniformFieldTracking(kFALSE);\r  
+  rec.SetWriteESDfriend(kTRUE);\r
+  rec.SetWriteAlignmentData();\r
+  rec.SetInput(filename);\r
+  //
+  //rec.SetRunLocalReconstruction("");\r
+  rec.SetRunReconstruction("TPC");\r
+  rec.SetFillESD("TPC");
+  rec.SetRunVertexFinder(kFALSE);
+  rec.SetRunQA(kFALSE);
+  //
+  //rec.SetEventRange(1,5);\r
+\r
+  // In case some detectors have to be switched off...\r
+  //  rec.SetRunLocalReconstruction("ALL");\r
+  //  rec.SetRunTracking("ALL");\r
+  //  rec.SetFillESD("ALL");\r
+  // Disable vertex finder for the moment\r
+  //  rec.SetRunVertexFinder(kFALSE);\r
+\r
+  // To be enabled if some equipment IDs are not set correctly by DAQ\r
+  //  rec.SetEquipmentIdMap("EquipmentIdMap.data");\r
+\r
+  // Detector options if any\r
+  rec.SetOption("MUON","SAVEDIGITS");\r
+  rec.SetOption("TPC","OldRCUFormat");\r
+  rec.SetOption("PHOS","OldRCUFormat");\r
+\r
+  // To be enabled when CTP readout starts\r
+  rec.SetFillTriggerESD(kFALSE);\r
+\r
+  // all events in one single file\r
+  rec.SetNumberOfEventsPerFile(-1);\r
+\r
+  // switch off cleanESD\r
+  rec.SetCleanESD(kFALSE);\r
+\r
+  //AliLog::SetGlobalDebugLevel(2);\r
+  rec.Run();\r
+\r
+  cout << "-----------------------------------------------------------------" << endl;\r
+  cout << "-----------------------------------------------------------------" << endl;\r
+  cout << "--------- Reconstruction Completed. Start merging QAs -----------" << endl;\r
+  cout << "-----------------------------------------------------------------" << endl;\r
+  cout << "-----------------------------------------------------------------" << endl;\r
+  AliQADataMakerSteer qas;\r
+  qas.Merge();\r
+}\r