//______________________________________________________________________________
void ChangeITSRecoParam(Int_t startRun=157560,const char* outputOCDB="alien://folder=/alice/cern.ch/user/l/laphecet/OCDB")
{
AliCDBManager::Instance()->SetDefaultStorage("raw://");
AliCDBManager::Instance()->SetRun(startRun);
AliCDBEntry* entry = AliCDBManager::Instance()->Get("ITS/Calib/RecoParam");
if (!entry) return;
TObject* o = entry->GetObject();
if (!o)
{
cout << "Could not get recoparams ? Oups" << endl;
return;
}
if ( o->IsA() != TObjArray::Class() )
{
cout << "This code only works with TObjArray recoparams. Sorry" << endl;
return;
}
TObjArray* array = static_cast<TObjArray*>(o);
for ( Int_t i = 0; i <= array->GetLast(); ++i )
{
AliDetectorRecoParam* p = static_cast<AliDetectorRecoParam*>(array->At(i));
if (AliRecoParam::Convert(p->GetEventSpecie())==AliRecoParam::kLowMult)
{
cout << Form("array[%d]=%s %s %s",i,
p ? p->ClassName() : "",
p ? AliRecoParam::GetEventSpecieName(AliRecoParam::Convert(p->GetEventSpecie())) :"",
p ? ( p->IsDefault() ? "default" : "") : "" ) << endl;
AliITSRecoParam* rp = dynamic_cast<AliITSRecoParam*>(p);
if (!rp)
{
cout << "OUPS. OUPS" << endl;
return;
}
rp->SetVertexer(4); // should be fast (i.e. MC) vertexer (an enum would be better for sure)
}
}
AliCDBManager::Instance()->SetDefaultStorage(outputOCDB);
AliMUONCDB::WriteToCDB(array, "ITS/Calib/RecoParam", startRun, AliCDBRunRange::Infinity(),
"reconstruction parameters for ITS, patched for fake vertex finder", "L. Aphecetche");
}
//______________________________________________________________________________
void ChangeRecoParam(Int_t startRun=195767,const char* outputOCDB="alien://folder=/alice/cern.ch/user/l/laphecet/OCDB")
{
AliCDBManager::Instance()->SetDefaultStorage("raw://");
AliCDBManager::Instance()->SetRun(99999999);
AliCDBEntry* entry = AliCDBManager::Instance()->Get("MUON/Calib/RecoParam");
if (!entry) return;
TObject* o = entry->GetObject();
if (!o)
{
cout << "Could not get recoparams ? Oups" << endl;
return;
}
if ( o->IsA() != TObjArray::Class() )
{
cout << "This code only works with TObjArray recoparams. Sorry" << endl;
return;
}
TObjArray* array = static_cast<TObjArray*>(o);
for ( Int_t i = 0; i <= array->GetLast(); ++i )
{
AliDetectorRecoParam* p = static_cast<AliDetectorRecoParam*>(array->At(i));
if (AliRecoParam::Convert(p->GetEventSpecie())==AliRecoParam::kLowMult)
{
cout << Form("array[%d]=%s %s %s",i,
p ? p->ClassName() : "",
p ? AliRecoParam::GetEventSpecieName(AliRecoParam::Convert(p->GetEventSpecie())) :"",
p ? ( p->IsDefault() ? "default" : "") : "" ) << endl;
p->Print("");
AliMUONRecoParam* rp = dynamic_cast<AliMUONRecoParam*>(p);
if (!rp)
{
cout << "OUPS. OUPS" << endl;
return;
}
rp->SetHVLimit(1,-99999);
rp->SetPadGoodnessMask(0x400be93);
rp->Print("");
}
}
AliCDBManager::Instance()->SetDefaultStorage(outputOCDB);
AliMUONCDB::WriteToCDB(array, "MUON/Calib/RecoParam", startRun, AliCDBRunRange::Infinity(),
"reconstruction parameters for MUON, patched for event size limits (for DQM)", "L. Aphecetche");
}
void MakeADRecoParamEntry(AliRecoParam::EventSpecie_t defaultEventSpecie=AliRecoParam::kLowMult, const char *outputCDB = "local://$ALICE_ROOT/OCDB") {
//========================================================================
//
// Steering macro for AD reconstruction parameters
//
// Author: Michal Broz
//
//========================================================================
const char* macroname = "MakeADRecoParam.C";
// Activate CDB storage and load geometry from CDB
AliCDBManager* cdb = AliCDBManager::Instance();
cdb->SetDefaultStorage(outputCDB);
cdb->SetRun(0);
TObjArray *recoParamArray = new TObjArray();
{
AliADRecoParam * ADRecoParam = new AliADRecoParam;
ADRecoParam->SetEventSpecie(AliRecoParam::kCosmic);
ADRecoParam->SetStartClock(0);
ADRecoParam->SetEndClock(20);
ADRecoParam->SetNPreClocks(1);
ADRecoParam->SetNPostClocks(10);
ADRecoParam->SetTimeWindowBBALow(-2.5);
ADRecoParam->SetTimeWindowBBAUp(2.5);
ADRecoParam->SetTimeWindowBGALow(-4.0);
ADRecoParam->SetTimeWindowBGAUp(4.0);
ADRecoParam->SetTimeWindowBBCLow(-1.5);
ADRecoParam->SetTimeWindowBBCUp(1.5);
ADRecoParam->SetTimeWindowBGCLow(-2.0);
ADRecoParam->SetTimeWindowBGCUp(2.0);
ADRecoParam->SetAdcThresHold(5);
ADRecoParam->SetMaxResid(1.5);
ADRecoParam->SetResidRise(0.02);
recoParamArray->AddLast(ADRecoParam);
}
{
AliADRecoParam * ADRecoParam = new AliADRecoParam;
ADRecoParam->SetStartClock(0);
ADRecoParam->SetEndClock(20);
ADRecoParam->SetNPreClocks(1);
ADRecoParam->SetNPostClocks(10);
ADRecoParam->SetTimeWindowBBALow(-2.5);
ADRecoParam->SetTimeWindowBBAUp(2.5);
ADRecoParam->SetTimeWindowBGALow(-4.0);
ADRecoParam->SetTimeWindowBGAUp(4.0);
ADRecoParam->SetTimeWindowBBCLow(-1.5);
ADRecoParam->SetTimeWindowBBCUp(1.5);
ADRecoParam->SetTimeWindowBGCLow(-2.0);
ADRecoParam->SetTimeWindowBGCUp(2.0);
ADRecoParam->SetAdcThresHold(5);
ADRecoParam->SetMaxResid(1.5);
ADRecoParam->SetResidRise(0.02);
ADRecoParam->SetEventSpecie(AliRecoParam::kLowMult);
recoParamArray->AddLast(ADRecoParam);
}
{
AliADRecoParam * ADRecoParam = new AliADRecoParam;
ADRecoParam->SetStartClock(9);
ADRecoParam->SetEndClock(11);
ADRecoParam->SetNPostClocks(6);
ADRecoParam->SetTimeWindowBBALow(-2.5);
ADRecoParam->SetTimeWindowBBAUp(2.5);
ADRecoParam->SetTimeWindowBGALow(-4.0);
ADRecoParam->SetTimeWindowBGAUp(4.0);
ADRecoParam->SetTimeWindowBBCLow(-1.5);
ADRecoParam->SetTimeWindowBBCUp(1.5);
ADRecoParam->SetTimeWindowBGCLow(-2.0);
ADRecoParam->SetTimeWindowBGCUp(2.0);
ADRecoParam->SetAdcThresHold(5);
ADRecoParam->SetMaxResid(1.5);
ADRecoParam->SetResidRise(0.02);
ADRecoParam->SetEventSpecie(AliRecoParam::kHighMult);
recoParamArray->AddLast(ADRecoParam);
}
// Set the defaultEventSpecie
Bool_t defaultIsSet = kFALSE;
for(Int_t i =0; i < recoParamArray->GetEntriesFast(); i++) {
AliDetectorRecoParam *param = (AliDetectorRecoParam *)recoParamArray->UncheckedAt(i);
if (!param) continue;
if (defaultEventSpecie & param->GetEventSpecie()) {
param->SetAsDefault();
defaultIsSet = kTRUE;
}
}
if (!defaultIsSet) {
Error(macroname,"The default reconstruction parameters are not set! Exiting...");
return;
}
// save in CDB storage
AliCDBMetaData *md= new AliCDBMetaData();
md->SetResponsible("Michal Broz");
md->SetComment("Reconstruction parameters for AD");
md->SetAliRootVersion(gSystem->Getenv("ARVERSION"));
md->SetBeamPeriod(0);
AliCDBId id("AD/Calib/RecoParam", 0, AliCDBRunRange::Infinity());
cdb->Put(recoParamArray, id, md);
return;
}
