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[24]={0.1,
  31                           0.1,
  32                           0.5,
  33                           //4, // minimum 4 clusters (default was 5)
  34                           3, // minimum 3 clusters (was 4)
  35                           1,
  36                           3.,
  37                           100.,
  38                           1000.,
  39                           3.,
  40                           30.,
  41                           6, // 6: MultiVertexer (was 1)
  42                           4,
  43                           // multivertexer settings
  44                           7.,
  45                           1e3.,
  46                           5.0,
  47                           0.05,
  48                           10e-4,
  49                           2.,
  50                           10.,
  51                           1.,
  52                           50.,
  53                           50.,
  54                           0.,
  55                           999999.,
  56                           3.
  57     };
  58     param->SetVertexerTracksCutsITS(24,cutsITS);
  59     param->SetVertexerTracksConstraintTPC(kTRUE);
  60     recoParamArray->AddLast(param);
  61   }
  62   {
  63     AliGRPRecoParam * param = AliGRPRecoParam::GetHighFluxParam();
  64     param->SetEventSpecie(AliRecoParam::kHighMult);
  65     param->SetVertexerTracksConstraintITS(kTRUE);
  66     Double_t cutsITS[24]={0.1,
  67                           0.1,
  68                           0.5,
  69                           4, // minimum 4 clusters (default was 5)
  70                           1,
  71                           3.,
  72                           100.,
  73                           1000.,
  74                           3.,
  75                           30.,
  76                           1,
  77                           1,
  78                           // multivertexer settings
  79                           7.,
  80                           1e3.,
  81                           5.0,
  82                           0.05,
  83                           10e-4,
  84                           2.,
  85                           10.,
  86                           1.,
  87                           50.,
  88                           0.,
  89                           999999.,
  90                           3.
  91     }; // faster finder algo for Iteration 0
  92     param->SetVertexerTracksCutsITS(24,cutsITS);
  93     param->SetVertexerTracksConstraintTPC(kTRUE);
  94     recoParamArray->AddLast(param);
  95   }
  96
  97   // Set the default
  98   Bool_t defaultIsSet = kFALSE;
  99   for(Int_t i =0; i < recoParamArray->GetEntriesFast(); i++) {
 100     AliDetectorRecoParam *par = (AliDetectorRecoParam *)recoParamArray->UncheckedAt(i);
 101     if (!par) continue;
 102     if (default & par->GetEventSpecie()) {
 103       par->SetAsDefault();
 104       defaultIsSet = kTRUE;
 105     }
 106   }
 107
 108   if (!defaultIsSet) {
 109     Error(macroname,"The default reconstruction parameters are not set! Exiting...");
 110     return;
 111   }
 112
 113   // save in CDB storage
 114   AliCDBMetaData *md= new AliCDBMetaData();
 115   md->SetResponsible("Cvetan Cheshkov");
 116   md->SetComment("GRP reconstruction parameters");
 117   md->SetAliRootVersion(gSystem->Getenv("ARVERSION"));
 118   md->SetBeamPeriod(0);
 119   AliCDBId id("GRP/Calib/RecoParam",0,AliCDBRunRange::Infinity());
 120   cdb->GetDefaultStorage()->Put(recoParamArray,id, md);
 121
 122   return;
 123 }