ATO-18 - a.) Removed upper limit on number of electrons. b.)Use qMax correction if...

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

diff --git a/TPC/recTPC.C b/TPC/recTPC.C
index 4503d58a827976c20de2b41ea7e86a19790ca059..de411864d36b67f3f92824721fdfd9af78bddb27 100644 (file)
--- a/TPC/recTPC.C
+++ b/TPC/recTPC.C
@@ -6,58 +6,66 @@
 //
 //
 
-void recTPC(const char *filename="data.root")
+void recTPC(Int_t type, const char *filename="data.root")
 {
   //
   // Set path to calibration data
   //
+  // type variable = 0 - cosmic test
+  //               = 1 - laser test   
   AliCDBManager * man = AliCDBManager::Instance();
-  man->SetDefaultStorage("local://$ALICE_ROOT");
-  man->SetRun(0);
-  man->SetSpecificStorage("TPC","local:///afs/cern.ch/user/m/mivanov/public/Calib");
+  man->SetDefaultStorage("local://$ALICE_ROOT/OCDB");
+  //man->SetRun(0);
+  //man->SetSpecificStorage("TPC/*/*","local:///afs/cern.ch/user/m/mivanov/public/Calib");
   //
   // Set reconstruction parameters
   //
-  AliLog::SetClassDebugLevel("AliTPCclustererMI",2);
-  AliTPCRecoParam * tpcRecoParam = AliTPCRecoParam::GetLaserTestParam(kTRUE);
+  AliLog::SetClassDebugLevel("AliTPCclusterer",2);
+  AliTPCRecoParam * tpcRecoParam = 0;
+  if (type==0)  tpcRecoParam = AliTPCRecoParam::GetCosmicTestParam(kTRUE);
+  if (type>0)  tpcRecoParam = AliTPCRecoParam::GetLaserTestParam(kTRUE);
+  tpcRecoParam->Dump();
   AliTPCReconstructor::SetRecoParam(tpcRecoParam);
   AliTPCReconstructor::SetStreamLevel(1);
   //
   //
   //
   AliReconstruction rec;  
-  rec.SetSpecificStorage("TPC","local:///afs/cern.ch/user/m/mivanov/public/Calib");
+  rec.SetDefaultStorage("local://$ALICE_ROOT/OCDB");
+  rec.SetSpecificStorage("TPC/*/*","local:///afs/cern.ch/user/m/mivanov/public/Calib");
   rec.SetLoadAlignData("");
   rec.SetWriteESDfriend(kTRUE);
   rec.SetInput(filename);
   rec.SetEquipmentIdMap("EquipmentIdMap.data");
   rec.SetRunReconstruction("TPC");
-  rec.SetOption("TPC","PedestalSubtraction OldRCUFormat");
+  rec.SetOption("TPC","PedestalSubtraction");
   //  rec.SetRunLocalReconstruction("");
   //  rec.SetRunTracking("TPC");
   rec.SetFillESD("TPC");
   rec.SetFillTriggerESD(kFALSE);
   rec.SetRunVertexFinder(kFALSE);
-  AliMagFMaps* field = new AliMagFMaps("Maps","Maps", 2, 1., 10., 1);
-  AliTracker::SetFieldMap(field,1);
   rec.SetWriteAlignmentData(kTRUE);
   rec.Run();
 }
 
-void recTracking(const char *filename="data.root", Int_t nevents=1)
+void recTracking(Int_t type, const char *filename="data.root")
 {
   //
   // Set path to calibration data
   //
   AliCDBManager * man = AliCDBManager::Instance();
-  man->SetDefaultStorage("local://$ALICE_ROOT");
+  man->SetDefaultStorage("local://$ALICE_ROOT/OCDB");
   man->SetRun(0);
-  man->SetSpecificStorage("TPC","local:///afs/cern.ch/user/m/mivanov/public/Calib");
+  man->SetSpecificStorage("TPC/*/*","local:///afs/cern.ch/user/m/mivanov/public/Calib");
   //
   // Set reconstruction parameters
   //
-  AliLog::SetClassDebugLevel("AliTPCclustererMI",2);
-  AliTPCRecoParam * tpcRecoParam = AliTPCRecoParam::GetLaserTestParam(kTRUE);
+  AliLog::SetClassDebugLevel("AliTPCclusterer",2);
+
+  AliTPCRecoParam * tpcRecoParam = 0;
+  if (type==0)  tpcRecoParam = AliTPCRecoParam::GetCosmicTestParam(kTRUE);
+  if (type>0)  tpcRecoParam = AliTPCRecoParam::GetLaserTestParam(kTRUE);
+
   AliTPCReconstructor::SetRecoParam(tpcRecoParam);
   AliTPCReconstructor::SetStreamLevel(1);
 
@@ -65,22 +73,20 @@ void recTracking(const char *filename="data.root", Int_t nevents=1)
   //
   //
   AliReconstruction rec;
-  //rec.SetSpecificStorage("TPC","local:///afs/cern.ch/user/m/mivanov/public/Calib");
+  rec.SetSpecificStorage("TPC/*/*","local:///afs/cern.ch/user/m/mivanov/public/Calib");
   rec.SetLoadAlignData("");
   rec.SetWriteESDfriend(kTRUE);
   rec.SetInput(filename);
   rec.SetEquipmentIdMap("EquipmentIdMap.data");
   //rec.SetRunReconstruction("TPC");
-  rec.SetOption("TPC","PedestalSubtraction OldRCUFormat");
+  rec.SetOption("TPC","PedestalSubtraction");
   rec.SetRunLocalReconstruction("");
   rec.SetRunTracking("TPC");
   rec.SetFillESD("TPC");
   rec.SetFillTriggerESD(kFALSE);
   rec.SetRunVertexFinder(kFALSE);
-  AliMagFMaps* field = new AliMagFMaps("Maps","Maps", 2, 1., 10., 1);
-  AliTracker::SetFieldMap(field,1);
   rec.SetWriteAlignmentData(kTRUE);
-  rec.Run(0,nevents);
+  rec.Run(0);
 }