test/cosmic/rec.C

   1 void rec(Int_t runNumber = 0, const char* year = "08", const char *localFileName = NULL)\r
   2 {\r
   3   // Offline shifter reconstruction macro\r
   4 \r
   5   cout << "Going to run the reconstruction for run: " << runNumber << endl;\r
   6 \r
   7   // connect to the grid \r
   8   TGrid * grid = 0x0 ; \r
   9   grid = TGrid::Connect("alien://") ; \r
  10                 \r
  11   // make the file name pattern year and run number\r
  12   TString pattern;\r
  13   pattern.Form("%9d",runNumber);\r
  14   pattern.ReplaceAll(" ", "0") ; \r
  15   pattern.Prepend(year);\r
  16   pattern.Append("*0.root");\r
  17 \r
  18   // find the files associated to this run\r
  19   // get the list of files from AliEn directly \r
  20   TString baseDir; \r
  21   baseDir.Form("/alice/data/20%s/",year);\r
  22 \r
  23   cout << "Looking for raw-data files with pattern " << pattern << " in folder " << baseDir << endl;\r
  24 \r
  25   TGridResult *result = grid->Query(baseDir, pattern);\r
  26 \r
  27   TList *fileList = result->GetFileInfoList();\r
  28 \r
  29   cout << fileList->GetEntries() << " raw-data files found" << endl;\r
  30   if ( fileList->GetEntries() == 0) {\r
  31     cout << "Exiting..." << endl;\r
  32     return;\r
  33   }\r
  34 \r
  35   // Take the first (or last?) file...\r
  36   TFileInfo *fi =  (TFileInfo *)fileList->At(0); \r
  37   //  TFileInfo *fi =  (TFileInfo *)fileList->At(fileList->GetEntries()-1); \r
  38 \r
  39   cout << "Getting the file:" << fi->GetCurrentUrl()->GetUrl() << endl;\r
  40   fi->Dump();\r
  41 \r
  42   TString filename(fi->GetCurrentUrl()->GetUrl());\r
  43 \r
  44   AliLog::Flush();\r
  45 \r
  46   /////////////////////////////////////////////////////////////////////////////////////////\r
  47   //\r
  48   // First version of the reconstruction\r
  49   // script for the FDR'08\r
  50 \r
  51   // Set the CDB storage location\r
  52   // AliLog::SetModuleDebugLevel("STEER",2);\r
  53   AliCDBManager * man = AliCDBManager::Instance();\r
  54   //  man->SetDefaultStorage("local://LocalCDB");\r
  55   man->SetDefaultStorage("alien://folder=/alice/data/2008/LHC08a/OCDB/");\r
  56   \r
  57   // Files that we can not read from alien...\r
  58   man->SetSpecificStorage("ITS/Calib/MapsAnodeSDD","local://$ALICE_ROOT");\r
  59   man->SetSpecificStorage("ITS/Calib/MapsTimeSDD","local://$ALICE_ROOT");\r
  60   man->SetSpecificStorage("TPC/Calib/ExB","local://$ALICE_ROOT");\r
  61 \r
  62   // Objects not found if using LHC07w database\r
  63   //  man->SetSpecificStorage("ITS/Calib/MapsAnodeSDD","local:///afs/cern.ch/user/c/cheshkov/public/OCDB");\r
  64   // man->SetSpecificStorage("GRP/GRP/Data","local://$ALICE_ROOT");\r
  65   // man->SetSpecificStorage("ITS/Calib/DDLMapSDD","local://$ALICE_ROOT");\r
  66   // man->SetSpecificStorage("MUON/Calib/Mapping","local://$ALICE_ROOT");\r
  67   // man->SetSpecificStorage("MUON/Calib/DDLStore","local://$ALICE_ROOT");\r
  68 \r
  69   // ITS settings\r
  70   AliITSRecoParam * itsRecoParam = AliITSRecoParam::GetCosmicTestParam();\r
  71   itsRecoParam->SetClusterErrorsParam(2);\r
  72   itsRecoParam->SetFindV0s(kFALSE);\r
  73   itsRecoParam->SetAddVirtualClustersInDeadZone(kFALSE);\r
  74   itsRecoParam->SetUseAmplitudeInfo(kFALSE);\r
  75   // In case we want to switch off a layer\r
  76   //  itsRecoParam->SetLayerToSkip(<N>);\r
  77   itsRecoParam->SetLayerToSkip(4);\r
  78   itsRecoParam->SetLayerToSkip(5);\r
  79   itsRecoParam->SetLayerToSkip(2);\r
  80   itsRecoParam->SetLayerToSkip(3);\r
  81   AliITSReconstructor::SetRecoParam(itsRecoParam);\r
  82 \r
  83   // TPC settings\r
  84   AliLog::SetClassDebugLevel("AliTPCclustererMI",2);\r
  85   AliTPCRecoParam * tpcRecoParam = AliTPCRecoParam::GetCosmicTestParam(kTRUE);\r
  86   tpcRecoParam->SetTimeInterval(60,940);\r
  87   tpcRecoParam->Dump();\r
  88   AliTPCReconstructor::SetRecoParam(tpcRecoParam);\r
  89   AliTPCReconstructor::SetStreamLevel(1);\r
  90 \r
  91   // PHOS settings\r
  92   AliPHOSRecoParam* recEmc = new AliPHOSRecoParamEmc();\r
  93   recEmc->SetSubtractPedestals(kTRUE);\r
  94   recEmc->SetMinE(0.05);\r
  95   recEmc->SetClusteringThreshold(0.10);\r
  96   AliPHOSReconstructor::SetRecoParamEmc(recEmc);\r
  97 \r
  98   // T0 settings\r
  99   AliLog::SetModuleDebugLevel("T0", 10);\r
 100 \r
 101   // MUON settings\r
 102   AliLog::SetClassDebugLevel("AliMUONRawStreamTracker",3);\r
 103   AliMUONRecoParam *muonRecoParam = AliMUONRecoParam::GetLowFluxParam();\r
 104   muonRecoParam->CombineClusterTrackReco(kTRUE);\r
 105   muonRecoParam->SetCalibrationMode("NOGAIN");\r
 106   //muonRecoParam->SetClusteringMode("PEAKFIT");\r
 107   //muonRecoParam->SetClusteringMode("PEAKCOG");\r
 108   muonRecoParam->Print("FULL");\r
 109   AliRecoParam::Instance()->RegisterRecoParam(muonRecoParam);\r
 110  \r
 111   // Tracking settings\r
 112   //  AliMagFMaps* field = new AliMagFMaps("Maps","Maps", 2, 1., 10., 1);\r
 113   AliMagFMaps* field = new AliMagFMaps("Maps","Maps", 2, 0., 10., 2);\r
 114   AliTracker::SetFieldMap(field,1);\r
 115   Double_t mostProbPt=0.35;\r
 116   AliExternalTrackParam::SetMostProbablePt(mostProbPt);\r
 117 \r
 118   // AliReconstruction settings\r
 119   AliReconstruction rec;\r
 120   rec.SetUniformFieldTracking(kFALSE);\r
 121   rec.SetWriteESDfriend(kTRUE);\r
 122   rec.SetWriteAlignmentData();\r
 123   rec.SetInput(filename.Data());\r
 124   rec.SetRunReconstruction("ALL");\r
 125   rec.SetUseTrackingErrorsForAlignment("ITS");\r
 126 \r
 127   // In case some detectors have to be switched off...\r
 128   //  rec.SetRunLocalReconstruction("ALL");\r
 129   //  rec.SetRunTracking("ALL");\r
 130   //  rec.SetFillESD("ALL");\r
 131   // Disable vertex finder for the moment\r
 132   rec.SetRunVertexFinder(kFALSE);\r
 133 \r
 134   // To be enabled if some equipment IDs are not set correctly by DAQ\r
 135   //  rec.SetEquipmentIdMap("EquipmentIdMap.data");\r
 136 \r
 137   // Detector options if any\r
 138   rec.SetOption("ITS","cosmics,onlyITS");\r
 139   rec.SetOption("MUON","SAVEDIGITS");\r
 140   rec.SetOption("TPC","OldRCUFormat");\r
 141   rec.SetOption("PHOS","OldRCUFormat");\r
 142 \r
 143   // To be enabled when CTP readout starts\r
 144   rec.SetFillTriggerESD(kFALSE);\r
 145 \r
 146   // all events in one single file\r
 147   rec.SetNumberOfEventsPerFile(-1);\r
 148 \r
 149   // switch off cleanESD\r
 150   rec.SetCleanESD(kFALSE);\r
 151 \r
 152   // rec.SetEventRange(0,15);\r
 153   // AliLog::SetGlobalDebugLevel(2);\r
 154 \r
 155   rec.SetRunQA(kFALSE);\r
 156   AliLog::Flush();\r
 157   rec.Run();\r
 158 \r
 159   //cout << "-----------------------------------------------------------------" << endl;\r
 160   //cout << "-----------------------------------------------------------------" << endl;\r
 161   //cout << "--------- Reconstruction Completed. Start merging QAs -----------" << endl;\r
 162   //cout << "-----------------------------------------------------------------" << endl;\r
 163   //cout << "-----------------------------------------------------------------" << endl;\r
 164   //AliQADataMakerSteer qas;\r
 165   //qas.Merge();\r
 166 }\r