GRP/MakeGRPRecoParam.C

   1 void MakeGRPRecoParam(AliRecoParam::EventSpecie_t default=AliRecoParam::kLowMult) {
   2 //========================================================================
   3 //
   4 // Steering macro for GRP reconstruction parameters
   5 //
   6 //
   7 //========================================================================
   8
   9
  10   const char* macroname = "MakeGRPRecoParam.C";
  11
  12   // Activate CDB storage and load geometry from CDB
  13   AliCDBManager* cdb = AliCDBManager::Instance();
  14   if(!cdb->IsDefaultStorageSet()) cdb->SetDefaultStorage("local://OCDB");
  15
  16   TObjArray *recoParamArray = new TObjArray();
  17
  18   {
  19     AliGRPRecoParam * param = AliGRPRecoParam::GetCosmicTestParam();
  20     param->SetEventSpecie(AliRecoParam::kCosmic);
  21     param->SetVertexerTracksConstraintITS(kFALSE);
  22     param->SetVertexerTracksConstraintTPC(kFALSE);
  23     recoParamArray->AddLast(param);
  24   }
  25   {
  26     // new settings for pass 2reco of Dec09 pp data
  27     AliGRPRecoParam * param = AliGRPRecoParam::GetLowFluxParam();
  28     param->SetEventSpecie(AliRecoParam::kLowMult);
  29     param->SetVertexerTracksConstraintITS(kTRUE);
  30     Double_t cutsITS[12]={0.1,
  31                           0.1,
  32                           0.5,
  33                           4, // minimum 4 clusters (default was 5)
  34                           1,
  35                           3.,
  36                           100.,
  37                           1000.,
  38                           3.,
  39                           30.,
  40                           1,
  41                           4};
  42     param->SetVertexerTracksCutsITS(12,cutsITS);
  43     param->SetVertexerTracksConstraintTPC(kTRUE);
  44     recoParamArray->AddLast(param);
  45   }
  46   {
  47     AliGRPRecoParam * param = AliGRPRecoParam::GetHighFluxParam();
  48     param->SetEventSpecie(AliRecoParam::kHighMult);
  49     param->SetVertexerTracksConstraintITS(kTRUE);
  50     Double_t cutsITS[12]={0.1,
  51                           0.1,
  52                           0.5,
  53                           4, // minimum 4 clusters (default was 5)
  54                           1,
  55                           3.,
  56                           100.,
  57                           1000.,
  58                           3.,
  59                           30.,
  60                           1,
  61                           1}; // faster finder algo for Iteration 0
  62     param->SetVertexerTracksCutsITS(12,cutsITS);
  63     param->SetVertexerTracksConstraintTPC(kTRUE);
  64     recoParamArray->AddLast(param);
  65   }
  66
  67   // Set the default
  68   Bool_t defaultIsSet = kFALSE;
  69   for(Int_t i =0; i < recoParamArray->GetEntriesFast(); i++) {
  70     AliDetectorRecoParam *par = (AliDetectorRecoParam *)recoParamArray->UncheckedAt(i);
  71     if (!par) continue;
  72     if (default & par->GetEventSpecie()) {
  73       par->SetAsDefault();
  74       defaultIsSet = kTRUE;
  75     }
  76   }
  77
  78   if (!defaultIsSet) {
  79     Error(macroname,"The default reconstruction parameters are not set! Exiting...");
  80     return;
  81   }
  82
  83   // save in CDB storage
  84   AliCDBMetaData *md= new AliCDBMetaData();
  85   md->SetResponsible("Cvetan Cheshkov");
  86   md->SetComment("GRP reconstruction parameters");
  87   md->SetAliRootVersion(gSystem->Getenv("ARVERSION"));
  88   md->SetBeamPeriod(0);
  89   AliCDBId id("GRP/Calib/RecoParam",0,AliCDBRunRange::Infinity());
  90   cdb->GetDefaultStorage()->Put(recoParamArray,id, md);
  91
  92   return;
  93 }