bool isIsomorphic(string s, string t) {
vector<int> mapS(256, -1), mapT(256, -1);
int n = s.size();
for (int i = 0; i < n; ++i) {
if (mapS[s[i]] != mapT[t[i]])
return false;
mapS[s[i]] = mapT[t[i]] = i;
}
return true;
}
void StoreMapsSDD(Int_t firstRun=0,Int_t lastRun=AliCDBRunRange::Infinity(), Bool_t optSmear=kFALSE){
///////////////////////////////////////////////////////////////////////
// Macro to generate and store the correction maps for SDD //
// Generates: //
// 1 file with 520 AliITSCorrMapSDD drift maps (MapsTimeSDD) //
///////////////////////////////////////////////////////////////////////
if(!AliCDBManager::Instance()->IsDefaultStorageSet()) {
AliCDBManager::Instance()->SetDefaultStorage("local://OCDB");
}
AliCDBMetaData *md = new AliCDBMetaData();
md->SetObjectClassName("TObjArray");
md->SetResponsible("Francesco Prino");
md->SetBeamPeriod(0);
md->SetComment("Simulated data");
AliCDBId mapT("ITS/Calib/MapsTimeSDD",firstRun,lastRun);
TObjArray tmap(520);
tmap.SetOwner(kFALSE);
TRandom3 *gran = new TRandom3();
AliITSCorrMapSDD* mapTime0;
AliITSCorrMapSDD* mapTime1;
for(Int_t mod=0;mod<260;mod++){
// maps
Char_t name[20];
sprintf(name,"DriftTimeMap_%d_%d\n",mod,0);
Int_t nbinsan=1;
if(optSmear && (mod==10 || mod==240)){
nbinsan=256;
sprintf(name,"DriftTimeMap_%d_%d\n",mod,0);
mapTime0 = new AliITSCorrMap2DSDD(name,nbinsan,72);
sprintf(name,"DriftTimeMap_%d_%d\n",mod,1);
mapTime1 = new AliITSCorrMap2DSDD(name,nbinsan,72);
}else{
sprintf(name,"DriftTimeMap_%d_%d\n",mod,0);
mapTime0 = new AliITSCorrMap1DSDD(name,72);
sprintf(name,"DriftTimeMap_%d_%d\n",mod,1);
mapTime1 = new AliITSCorrMap1DSDD(name,72);
}
for(Int_t nan = 0;nan< nbinsan;nan++){
for(Int_t nt = 0;nt<36*2;nt++){
Double_t cnt0=0.;
Double_t cnt1=0.;
if(optSmear){
cnt0=gran->Gaus(0,20);
cnt1=gran->Gaus(0,20);
}
mapTime0->SetCellContent(nan,nt,cnt0);
mapTime1->SetCellContent(nan,nt,cnt1);
}
}
tmap.Add(mapTime0);
tmap.Add(mapTime1);
printf("Added module %d\n",mod);
}
AliCDBManager::Instance()->GetDefaultStorage()->Put(&tmap, mapT, md);
}
