Correct shell scipt and macro to be used with 2009 data.
[u/mrichter/AliRoot.git] / test / cosmic / rec.C
index 8d54b56..b2fc5c0 100644 (file)
@@ -1,14 +1,11 @@
-void rec(const char *filename="raw.root", const Int_t mfield=1)\r
+void rec(const char *filename="raw.root")\r
 {\r
   /////////////////////////////////////////////////////////////////////////////////////////\r
   //\r
-  // Second version of the reconstruction\r
-  // script for the 2008 cosmic data (LHC08b) \r
+  // First version of the reconstruction\r
+  // script for the 2009 data (LHC09b) \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
@@ -17,130 +14,32 @@ void rec(const char *filename="raw.root", const Int_t mfield=1)
 \r
   // Set the CDB storage location\r
   AliCDBManager * man = AliCDBManager::Instance();\r
-//     man->SetDefaultStorage("local://$ALICE_ROOT/OCDB");\r
-    man->SetDefaultStorage("alien://folder=/alice/data/2008/LHC08d/OCDB/");\r
-    man->SetSpecificStorage("ITS/Calib/*","local://$ALICE_ROOT/OCDB");\r
+  man->SetDefaultStorage("raw://");\r
   \r
-  // Example in case a specific CDB storage is needed\r
-  //  man->SetSpecificStorage("ITS/Calib/MapsAnodeSDD","local://$ALICE_ROOT/OCDB");\r
-\r
-  AliQA::SetQARefStorage("local://$ALICE_ROOT/QAref") ;\r
-\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(2);\r
-  itsRecoParam->SetLayerToSkip(3);\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(kFALSE);\r
-  tpcRecoParam->SetTimeInterval(60,940);\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
-  \r
-  // TRD reconstruction params\r
-  AliTRDrecoParam *fTRDrecoParam = AliTRDrecoParam::GetCosmicTestParam();\r
-  rec.SetRecoParam("TRD",fTRDrecoParam);\r
-  AliTRDtrackerV1::SetNTimeBins(30);\r
-\r
-  // PHOS settings\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
+  // QA options\r
+  rec.SetRunQA("ALL:ALL") ;\r
+  rec.SetQARefDefaultStorage("local://$ALICE_ROOT/QAref") ;\r
 \r
-  // MUON settings\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
-  rec.SetRecoParam("MUON",muonRecoParam);\r
\r
-  // Tracking settings\r
-  AliMagF* field;\r
-  if (mfield)\r
-    field = new AliMagF("Maps","Maps", 2, 1., 1., 10., AliMagF::k5kG);\r
-   else\r
-     field = new AliMagF("Maps","Maps", 2, 0., 0., 10., AliMagF::k2kG);\r
-\r
-  Double_t mostProbPt=0.35;\r
-  AliExternalTrackParam::SetMostProbablePt(mostProbPt);\r
+  AliGRPRecoParam *grpRecoParam = AliGRPRecoParam::GetLowFluxParam();\r
+  grpRecoParam->SetVertexerTracksConstraintITS(kFALSE);\r
+  grpRecoParam->SetVertexerTracksConstraintTPC(kFALSE);\r
+  grpRecoParam->SetMostProbablePt(3.0);\r
+  rec.SetRecoParam("GRP",grpRecoParam);\r
 \r
   // AliReconstruction settings\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.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
-  // 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("ITS","cosmics");\r
-  rec.SetOption("MUON","SAVEDIGITS");\r
-  rec.SetOption("T0","cosmic");\r
-\r
-  // Enabled when CTP readout starts\r
-  rec.SetFillTriggerESD(kTRUE);\r
-\r
-  // all events in one single file\r
-  rec.SetNumberOfEventsPerFile(-1);\r
-\r
   // switch off cleanESD\r
   rec.SetCleanESD(kFALSE);\r
 \r
-  //rec.SetEventRange(0,15);\r
-\r
-  rec.SetRunQA(":");\r
-  rec.SetRunGlobalQA(kTRUE);\r
   AliLog::Flush();\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