Changes in reco. parameters according to revision 27529
authorhristov <hristov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 23 Jul 2008 13:05:24 +0000 (13:05 +0000)
committerhristov <hristov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 23 Jul 2008 13:05:24 +0000 (13:05 +0000)
test/PbPbbench/rec.C
test/PbPbbench/recraw/rec.C
test/cosmic/rec.C
test/fpprod/rec.C
test/genkine/sim/rec.C
test/merge/backgr/rec.C
test/merge/signal/rec.C
test/ppbench/rec.C
test/ppbench/recraw/rec.C
test/pploadlibs/rec.C

index b4b303b..936f7ed 100644 (file)
@@ -3,8 +3,10 @@ void rec() {
 
   reco.SetWriteESDfriend();
   reco.SetWriteAlignmentData();
-  AliTPCRecoParam * tpcRecoParam = AliTPCRecoParam::GetHighFluxParam();
-  AliTPCReconstructor::SetRecoParam(tpcRecoParam);
+  reco.SetRecoParam("TPC",AliTPCRecoParam::GetHighFluxParam());
+  reco.SetRecoParam("TRD",AliTRDrecoParam::GetHighFluxParam());
+  reco.SetRecoParam("PHOS",AliPHOSRecoParam::GetDefaultParameters());
+  reco.SetRecoParam("MUON",AliMUONRecoParam::GetHighFluxParam());
   //  AliTPCReconstructor::SetStreamLevel(1);
   reco.SetRunVertexFinderTracks(kFALSE);
 
index 5cd73b9..194a5ed 100644 (file)
@@ -10,9 +10,10 @@ void rec() {
   reco.SetUniformFieldTracking(kFALSE);
   reco.SetWriteESDfriend();
   reco.SetWriteAlignmentData();
-  AliTPCRecoParam * tpcRecoParam = AliTPCRecoParam::GetHighFluxParam();
-  AliTPCReconstructor::SetRecoParam(tpcRecoParam);
-  //  AliTPCReconstructor::SetStreamLevel(1);
+  reco.SetRecoParam("TPC",AliTPCRecoParam::GetHighFluxParam());
+  reco.SetRecoParam("TRD",AliTRDrecoParam::GetHighFluxParam());
+  reco.SetRecoParam("PHOS",AliPHOSRecoParam::GetDefaultParameters());
+  reco.SetRecoParam("MUON",AliMUONRecoParam::GetHighFluxParam());
 
   reco.SetRunReconstruction("ITS TPC TRD TOF HMPID PHOS EMCAL MUON FMD PMD ZDC T0 VZERO");
   reco.SetInput("raw.root");
index c66a10f..4588ffd 100644 (file)
-void rec(const char *filename="raw.root")\r
+void rec(const char *filename="raw.root", const Int_t mfield=1)\r
 {\r
   /////////////////////////////////////////////////////////////////////////////////////////\r
   //\r
-  // First version of the reconstruction\r
-  // script for the FDR'08\r
+  // Second version of the reconstruction\r
+  // script for the 2008 cosmic data (LHC08b) \r
+  //\r
+  /////////////////////////////////////////////////////////////////////////////////////////\r
+  //AliLog::SetGlobalLogLevel(AliLog::kWarning);\r
+  AliLog::SetGlobalLogLevel(AliLog::kError);\r
+\r
+  gSystem->Load("libRAliEn.so");\r
+  gSystem->Load("libNet.so");\r
+  gSystem->Load("libMonaLisa.so");\r
+  new TMonaLisaWriter(0, "GridAliRoot-rec.C", 0, 0, "global");\r
+  gSystem->Setenv("APMON_INTERVAL", "120");\r
 \r
   // Set the CDB storage location\r
-  // AliLog::SetModuleDebugLevel("STEER",2);\r
   AliCDBManager * man = AliCDBManager::Instance();\r
-  //  man->SetDefaultStorage("local://LocalCDB");\r
-  man->SetDefaultStorage("alien://folder=/alice/data/2008/LHC08a/OCDB/");\r
+  //  man->SetDefaultStorage("local://$ALICE_ROOT");\r
+  man->SetDefaultStorage("alien://folder=/alice/data/2008/LHC08b/OCDB/");\r
+  man->SetSpecificStorage("ITS/Calib/*","local://$ALICE_ROOT");\r
   \r
-  // Files that we can not read from alien...solved\r
+  // Example in case a specific CDB storage is needed\r
   //  man->SetSpecificStorage("ITS/Calib/MapsAnodeSDD","local://$ALICE_ROOT");\r
-  //  man->SetSpecificStorage("ITS/Calib/MapsTimeSDD","local://$ALICE_ROOT");\r
-  //  man->SetSpecificStorage("TPC/Calib/ExB","local://$ALICE_ROOT");\r
 \r
-  // Objects not found if using LHC07w database...solved\r
-  //  man->SetSpecificStorage("ITS/Calib/MapsAnodeSDD","local:///afs/cern.ch/user/c/cheshkov/public/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
+  // Reconstruction settings\r
+  AliReconstruction rec;\r
 \r
   // ITS settings\r
   AliITSRecoParam * itsRecoParam = AliITSRecoParam::GetCosmicTestParam();\r
+  itsRecoParam->SetFactorSAWindowSizes(20);\r
   itsRecoParam->SetClusterErrorsParam(2);\r
   itsRecoParam->SetFindV0s(kFALSE);\r
   itsRecoParam->SetAddVirtualClustersInDeadZone(kFALSE);\r
   itsRecoParam->SetUseAmplitudeInfo(kFALSE);\r
   // In case we want to switch off a layer\r
   //  itsRecoParam->SetLayerToSkip(<N>);\r
-  itsRecoParam->SetLayerToSkip(4);\r
-  itsRecoParam->SetLayerToSkip(5);\r
+  //  itsRecoParam->SetLayerToSkip(4);\r
+  //  itsRecoParam->SetLayerToSkip(5);\r
   itsRecoParam->SetLayerToSkip(2);\r
   itsRecoParam->SetLayerToSkip(3);\r
-  AliITSReconstructor::SetRecoParam(itsRecoParam);\r
+  //itsRecoParam->SetSAOnePointTracks();\r
+  itsRecoParam->SetClusterMisalError(0.1); // [cm]\r
+  itsRecoParam->SetSAUseAllClusters();\r
+  rec.SetRecoParam("ITS",itsRecoParam);\r
 \r
   // TPC settings\r
-  AliLog::SetClassDebugLevel("AliTPCclustererMI",2);\r
-  AliTPCRecoParam * tpcRecoParam = AliTPCRecoParam::GetCosmicTestParam(kTRUE);\r
+  //AliLog::SetClassDebugLevel("AliTPCclustererMI",2);\r
+  AliTPCRecoParam * tpcRecoParam = AliTPCRecoParam::GetCosmicTestParam(kFALSE);\r
   tpcRecoParam->SetTimeInterval(60,940);\r
-  tpcRecoParam->Dump();\r
-  AliTPCReconstructor::SetRecoParam(tpcRecoParam);\r
+  Double_t sysError[5]={0.3,1, 0.3/150., 1./150.,0.3/(150*150.)};\r
+  tpcRecoParam->SetSystematicError(sysError);\r
+  tpcRecoParam->SetMinMaxCutAbs(4.);\r
+  tpcRecoParam->SetMinLeftRightCutAbs(6.);\r
+  tpcRecoParam->SetMinUpDownCutAbs(6.);\r
+  //  tpcRecoParam->Dump();\r
+  rec.SetRecoParam("TPC",tpcRecoParam);\r
   AliTPCReconstructor::SetStreamLevel(1);\r
 \r
   // TRD setting\r
+  // Settings for the TRD Raw Reader\r
   AliTRDrawStreamBase::SetRawStreamVersion("TB");\r
+  AliTRDrawStreamTB::SetNoErrorWarning();\r
+  AliTRDrawStreamTB::AllowCorruptedData();\r
+  AliTRDrawStreamTB::DisableStackNumberChecker();\r
+  AliTRDrawStreamTB::DisableStackLinkNumberChecker();\r
+  AliTRDrawStreamTB::SetSubtractBaseline(10);\r
+  \r
+  // TRD reconstruction params\r
+  AliTRDrecoParam *fTRDrecoParam = AliTRDrecoParam::GetCosmicTestParam();\r
+  rec.SetRecoParam("TRD",fTRDrecoParam);\r
+  AliTRDtrackerV1::SetNTimeBins(30);\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
+  AliPHOSRecoParam* recPHOS = new AliPHOSRecoParam();\r
+  recPHOS->SetEMCSubtractPedestals(kTRUE);\r
+  recPHOS->SetEMCMinE(0.05);\r
+  recPHOS->SetEMCClusteringThreshold(0.10);\r
+  rec.SetRecoParam("PHOS",recPHOS);\r
 \r
   // T0 settings\r
-  AliLog::SetModuleDebugLevel("T0", 10);\r
+  //AliLog::SetModuleDebugLevel("T0", 10);\r
 \r
   // MUON settings\r
-  AliLog::SetClassDebugLevel("AliMUONRawStreamTracker",3);\r
+  //AliLog::SetClassDebugLevel("AliMUONRawStreamTracker",3);\r
   AliMUONRecoParam *muonRecoParam = AliMUONRecoParam::GetLowFluxParam();\r
   muonRecoParam->CombineClusterTrackReco(kTRUE);\r
   muonRecoParam->SetCalibrationMode("NOGAIN");\r
   //muonRecoParam->SetClusteringMode("PEAKFIT");\r
   //muonRecoParam->SetClusteringMode("PEAKCOG");\r
   muonRecoParam->Print("FULL");\r
-  AliRecoParam::Instance()->RegisterRecoParam(muonRecoParam);\r
+  rec.SetRecoParam("MUON",muonRecoParam);\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
+  AliMagFMaps* field;\r
+  if (mfield)\r
+     field = new AliMagFMaps("Maps","Maps", 2, 1., 10., AliMagFMaps::k5kG);\r
+   else\r
+     field = new AliMagFMaps("Maps","Maps", 2, 0., 10., 2);\r
   AliTracker::SetFieldMap(field,1);\r
   Double_t mostProbPt=0.35;\r
   AliExternalTrackParam::SetMostProbablePt(mostProbPt);\r
 \r
   // AliReconstruction settings\r
-  AliReconstruction rec;\r
   rec.SetUniformFieldTracking(kFALSE);\r
   rec.SetWriteESDfriend(kTRUE);\r
   rec.SetWriteAlignmentData();\r
   rec.SetInput(filename);\r
-  rec.SetRunReconstruction("ALL");\r
+  //  rec.SetRunReconstruction("ALL");\r
   rec.SetUseTrackingErrorsForAlignment("ITS");\r
 \r
   // In case some detectors have to be switched off...\r
-  //  rec.SetRunLocalReconstruction("ALL");\r
-  //  rec.SetRunTracking("ALL");\r
-  //  rec.SetFillESD("ALL");\r
+  rec.SetRunReconstruction("ITS TPC TRD TOF HMPID PHOS MUON FMD PMD T0 VZERO ZDC ACORDE");\r
+\r
   // Enable vertex finder - it is needed for cosmic track reco\r
   rec.SetRunVertexFinder(kTRUE);\r
 \r
@@ -95,14 +120,12 @@ void rec(const char *filename="raw.root")
   //  rec.SetEquipmentIdMap("EquipmentIdMap.data");\r
 \r
   // Detector options if any\r
-  rec.SetOption("ITS","cosmics,onlyITS");\r
+  rec.SetOption("ITS","cosmics");\r
   rec.SetOption("MUON","SAVEDIGITS");\r
-  rec.SetOption("TPC","OldRCUFormat");\r
-  rec.SetOption("PHOS","OldRCUFormat");\r
   rec.SetOption("T0","cosmic");\r
 \r
-  // To be enabled when CTP readout starts\r
-  rec.SetFillTriggerESD(kFALSE);\r
+  // Enabled when CTP readout starts\r
+  rec.SetFillTriggerESD(kTRUE);\r
 \r
   // all events in one single file\r
   rec.SetNumberOfEventsPerFile(-1);\r
@@ -110,10 +133,10 @@ void rec(const char *filename="raw.root")
   // switch off cleanESD\r
   rec.SetCleanESD(kFALSE);\r
 \r
-  // rec.SetEventRange(0,15);\r
-  // AliLog::SetGlobalDebugLevel(2);\r
+  //rec.SetEventRange(0,15);\r
 \r
-  rec.SetRunQA(kFALSE);\r
+  rec.SetRunQA("ITS TPC:ESD RECPOINT");\r
+  rec.SetRunGlobalQA(kTRUE);\r
   AliLog::Flush();\r
   rec.Run();\r
 \r
index 8564ef8..78a419f 100644 (file)
@@ -9,20 +9,18 @@ void rec() {
 
   reco.SetWriteESDfriend();
   reco.SetWriteAlignmentData();
-  AliTPCRecoParam * tpcRecoParam = AliTPCRecoParam::GetLowFluxParam();
-  AliTPCReconstructor::SetRecoParam(tpcRecoParam);
-  AliTPCReconstructor::SetStreamLevel(1);
+  reco.SetRecoParam("TPC",AliTPCRecoParam::GetLowFluxParam());
+  reco.SetRecoParam("TRD",AliTRDrecoParam::GetLowFluxParam());
+  reco.SetRecoParam("PHOS",AliPHOSRecoParam::GetDefaultParameters());
+  reco.SetRecoParam("MUON",AliMUONRecoParam::GetLowFluxParam());
   //   reco.SetInput("raw.root");
   reco.SetRunReconstruction("ITS TPC TRD TOF HMPID PHOS EMCAL MUON VZERO T0 FMD PMD ZDC");
-  reco.SetRunQA("");
-  reco.SetRunGlobalQA(kFALSE);
+  reco.SetRunQA(":");
+  reco.SetRunGlobalQA(kTRUE);
   reco.SetMeanVertexConstraint(kFALSE);
 
 // **** The field map settings must be the same as in Config.C !
   AliMagWrapCheb* field = 0x0;
-  //field = new AliMagWrapCheb("Maps","Maps", 2, 0., 10., AliMagWrapCheb::k2kG);
-  //Bool_t uniform=kTRUE;
-  //AliTracker::SetFieldMap(field,uniform); // tracking with the uniform field
   field = new AliMagWrapCheb("Maps","Maps", 2, 1., 10., AliMagWrapCheb::k5kG);
   Bool_t uniform=kFALSE;
   AliTracker::SetFieldMap(field,uniform);  // tracking with the real map
index 127d0d3..9c7219c 100644 (file)
@@ -1,6 +1,10 @@
 void rec() {
   AliReconstruction reco;
 
+  reco.SetRecoParam("TPC",AliTPCRecoParam::GetLowFluxParam());
+  reco.SetRecoParam("TRD",AliTRDrecoParam::GetLowFluxParam());
+  reco.SetRecoParam("PHOS",AliPHOSRecoParam::GetDefaultParameters());
+  reco.SetRecoParam("MUON",AliMUONRecoParam::GetLowFluxParam());
   //  reco.SetRunReconstruction("ITS TPC TRD TOF HMPID FMD PMD VZERO START MUON ZDC");
 
 // **** The field map settings must be the same as in Config.C !
index 4f50107..6b2161f 100644 (file)
@@ -3,8 +3,10 @@ void rec() {
   AliReconstruction reco;
   reco.SetWriteESDfriend();
   reco.SetWriteAlignmentData();
-  AliTPCReconstructor::SetStreamLevel(1);
-  AliTPCReconstructor::SetRecoParam(AliTPCRecoParam::GetLowFluxParam());
+  reco.SetRecoParam("TPC",AliTPCRecoParam::GetLowFluxParam());
+  reco.SetRecoParam("TRD",AliTRDrecoParam::GetLowFluxParam());
+  reco.SetRecoParam("PHOS",AliPHOSRecoParam::GetDefaultParameters());
+  reco.SetRecoParam("MUON",AliMUONRecoParam::GetLowFluxParam());
   //  reco.SetInput("raw.root");
 
 // **** The field map settings must be the same as in Config.C !
index 5888f6b..dcb4c92 100644 (file)
@@ -2,8 +2,10 @@ void rec() {
   AliReconstruction reco;
   reco.SetWriteESDfriend();
   reco.SetWriteAlignmentData();
-  AliTPCReconstructor::SetStreamLevel(1);
-  AliTPCReconstructor::SetRecoParam(AliTPCRecoParam::GetLowFluxParam());
+  reco.SetRecoParam("TPC",AliTPCRecoParam::GetLowFluxParam());
+  reco.SetRecoParam("TRD",AliTRDrecoParam::GetLowFluxParam());
+  reco.SetRecoParam("PHOS",AliPHOSRecoParam::GetDefaultParameters());
+  reco.SetRecoParam("MUON",AliMUONRecoParam::GetLowFluxParam());
   //  reco.SetInput("raw.root");
 
 // **** The field map settings must be the same as in Config.C !
index 3e7a77f..152d794 100644 (file)
@@ -3,9 +3,10 @@ void rec() {
 
   reco.SetWriteESDfriend();
   reco.SetWriteAlignmentData();
-  AliTPCRecoParam * tpcRecoParam = AliTPCRecoParam::GetLowFluxParam();
-  AliTPCReconstructor::SetRecoParam(tpcRecoParam);
-  AliTPCReconstructor::SetStreamLevel(1);
+  reco.SetRecoParam("TPC",AliTPCRecoParam::GetLowFluxParam());
+  reco.SetRecoParam("TRD",AliTRDrecoParam::GetLowFluxParam());
+  reco.SetRecoParam("PHOS",AliPHOSRecoParam::GetDefaultParameters());
+  reco.SetRecoParam("MUON",AliMUONRecoParam::GetLowFluxParam());
   //   reco.SetInput("raw.root");
 
 // **** The field map settings must be the same as in Config.C !
index 82118ee..c97d59f 100644 (file)
@@ -10,15 +10,19 @@ void rec() {
   reco.SetUniformFieldTracking(kFALSE);
   reco.SetWriteESDfriend();
   reco.SetWriteAlignmentData();
-  AliTPCRecoParam * tpcRecoParam = AliTPCRecoParam::GetLowFluxParam();
-  AliTPCReconstructor::SetRecoParam(tpcRecoParam);
-  //  AliTPCReconstructor::SetStreamLevel(1);
+  reco.SetRecoParam("TPC",AliTPCRecoParam::GetLowFluxParam());
+  reco.SetRecoParam("TRD",AliTRDrecoParam::GetLowFluxParam());
+  reco.SetRecoParam("PHOS",AliPHOSRecoParam::GetDefaultParameters());
+  reco.SetRecoParam("MUON",AliMUONRecoParam::GetLowFluxParam());
 
   reco.SetRunReconstruction("ITS TPC TRD TOF HMPID PHOS EMCAL MUON FMD PMD ZDC T0 VZERO");
   reco.SetInput("raw.root");
 
   reco.SetNumberOfEventsPerFile(-1); // all events in one single file
 
+  reco.SetRunQA(":");
+  reco.SetRunGlobalQA(kTRUE);
+
   TStopwatch timer;
   timer.Start();
   reco.Run();
index 51fa34e..d8d8b73 100644 (file)
@@ -8,9 +8,10 @@ void rec() {
 
   reco.SetWriteESDfriend();
   reco.SetWriteAlignmentData();
-  AliTPCRecoParam * tpcRecoParam = AliTPCRecoParam::GetLowFluxParam();
-  AliTPCReconstructor::SetRecoParam(tpcRecoParam);
-  AliTPCReconstructor::SetStreamLevel(1);
+  reco.SetRecoParam("TPC",AliTPCRecoParam::GetLowFluxParam());
+  reco.SetRecoParam("TRD",AliTRDrecoParam::GetLowFluxParam());
+  reco.SetRecoParam("PHOS",AliPHOSRecoParam::GetDefaultParameters());
+  reco.SetRecoParam("MUON",AliMUONRecoParam::GetLowFluxParam());
   //   reco.SetInput("raw.root");
   reco.SetRunReconstruction("ITS TPC TRD TOF HMPID PHOS EMCAL MUON VZERO T0 FMD PMD ZDC");
 // **** The field map settings must be the same as in Config.C !