#if !defined(__CINT__) || defined(__MAKECINT__) #include "AliCDBManager.h" #include "AliITSCalibrationSDD.h" #include "AliITSresponseSDD.h" #include "AliCDBMetaData.h" #include "AliCDBStorage.h" #include "AliCDBId.h" #include "AliCDBPath.h" #include "AliCDBEntry.h" #include #include #endif void StoreCalibSDD(Int_t firstRun=0,Int_t lastRun=9999999 ){ /////////////////////////////////////////////////////////////////////// // Macro to generate and store the calibration files for SDD // // Generates: // // 1 file with 260 AliITSCalibrationSDD objects with // // baselines, noise, gain, drift speed for each module (CalibSDD) // // 1 file with the AliITSrespionseSDD object (RespSDD) // /////////////////////////////////////////////////////////////////////// if(!AliCDBManager::Instance()->IsDefaultStorageSet()) { AliCDBManager::Instance()->SetDefaultStorage("local://OCDB"); } AliCDBMetaData *md1= new AliCDBMetaData(); // metaData describing the object md1->SetObjectClassName("AliITSCalibration"); md1->SetResponsible("Elisabetta Crescio, Francesco Prino"); md1->SetBeamPeriod(0); md1->SetAliRootVersion("Head 20 nov. 2007"); //root version md1->SetComment("This is a test"); // TObject* str; // md1->SetProperty("key1",str); AliCDBMetaData *md2 = new AliCDBMetaData(); md2->SetObjectClassName("AliITSresponse"); md2->SetResponsible("Elisabetta Crescio, Francesco Prino"); md2->SetBeamPeriod(0); md2->SetAliRootVersion("Head 20 nov. 2007"); //root version md2->SetComment("This is a test"); AliCDBId idCalSDD("ITS/Calib/CalibSDD",firstRun, lastRun); TObjArray respSDD(260); respSDD.SetOwner(kFALSE); AliCDBId idRespSDD("ITS/Calib/RespSDD",firstRun, lastRun); AliITSresponseSDD* rd = new AliITSresponseSDD(); const Int_t nData = 209; Int_t anodeUp[209] = {0,36,0,12,20,32,0,0,12,76,28,8,16,0,0,0,8,0,0,0,20,4,0,0,0,0,0,0 ,0,0,8,0,0,0,0,0,0,0,0,0,0,0,12,0,8,0,4,4,0,160,0,0,0,252,16,0,8,8 ,0,0,12,0,0,0,12,0,15,20,0,0,0,0,0,0,0,12,12,0,0,0,0,0,8,8,3,240,212 ,12,0,8,0,12,0,0,12,0,8,24,0,0,12,0,0,0,0,40,0,0,40,12,28,0,0,12,12 ,0,0,0,0,20,0,0,0,12,0,24,0,0,0,0,0,0,0,8,16,0,0,0,0,0,256,0,0,0,0,0,20 ,0,12,0,0,0,0,24,0,0,0,0,0,0,0,20,0,0,16,0,0,0,0,24,0,0,0,8,0,16,40,0 ,0,0,0,0,0,0,0,0,4,0,32,12,8,28,0,76,0,0,0,12,60,144,0,0,0,0,16,0,16,0,3 }; Int_t anodeDown[209] = {0,8,0,48,0,12,0,8,0,80,12,0,4,4,0,0,24,0,0,20,0,0,0,0,20,0,0,0,0,0,0,0 ,0,0,0,0,4,0,24,4,0,0,8,0,0,36,20,0,12,236,0,0,12,256,16,8,32,4,0,0,24 ,24,10,0,16,8,0,2,40,0,0,0,24,0,0,0,8,0,0,0,0,0,32,240,0,92,60,28,8,0,0 ,2,0,0,0,0,12,48,0,0,0,0,0,36,11,12,0,0,0,12,12,11,0,20,0,0,12,20,0,0,4 ,0,8,12,0,0,0,16,16,0,32,72,12,0,88,20,16,112,8,0,244,28,256,28,0,24,236 ,56,0,68,0,4,20,208,20,12,4,28,12,0,0,20,12,0,100,0,16,8,8,0,24,16,0,12,12 ,16,0,16,20,0,28,0,8,24,0,12,8,4,40,0,104,96,32,140,20,12,8,20,24,16,16,20 ,8,140,96,0,32,20,44}; Float_t drVelParam[4]={7.75,0.002344,-0.000009,0}; Float_t edrVelParam[4]={0,0,0,0}; Float_t drVel[4]; TRandom3 *gran = new TRandom3(); for(Int_t mod=0;mod<260;mod++){ AliITSCalibrationSDD* resd = new AliITSCalibrationSDD("simulated"); Int_t nBadUp = 0; Int_t nBadDown = 0; // gain for(Int_t iWing=0; iWing<2;iWing++){ for(Int_t iChip=0; iChip<4;iChip++){ Float_t chipgain=gran->Gaus(1.,0.1); if(chipgain<0.1) chipgain=0.1; for(Int_t iChan=0; iChan<64;iChan++){ Float_t gain=gran->Gaus(chipgain,0.01); if(gain<0.1) gain=0.1; resd->SetGain(gain,iWing,iChip,iChan); } } } // bad channels Int_t count=0; do { Int_t ranMod = Int_t(nData*gran->Uniform()); nBadUp = anodeUp[ranMod]; nBadDown = anodeDown[ranMod]; }while (nBadUp+nBadDown>25); resd->SetDeadChannels(nBadDown+nBadUp); Int_t remainingBad = nBadDown; while (remainingBad>0) { Int_t nBadChannel; if (remainingBad<5) { nBadChannel = remainingBad; } else { Int_t width = remainingBad-(5-1); if (width>4) width = 4; nBadChannel = 5 + Int_t(width*gran->Uniform()); } Int_t iChannelPos = Int_t( (4*64-nBadChannel)*gran->Uniform() ); // Int_t iChip = iChannelPos/64; // Int_t iChan = iChannelPos - iChip*64; if(resd->IsBadChannel(iChannelPos)) continue; Int_t *clus = new Int_t[nBadChannel]; Int_t ich = iChannelPos; for(Int_t i=0;iIsBadChannel(clus[i])) break; resd->SetBadChannel(count,clus[i]); count++; } remainingBad -= nBadChannel; delete [] clus; } // May happen that, due to overlapping clusters, we // have less bad channels than requested // Let's put the remaining one per one ... Int_t nSeToBad = 0; for (Int_t i=0; i<4; i++){ for(Int_t j=0;j<64;j++){ if (resd->Gain(0,i,j)<0.0001) nSeToBad++; } } while (nSeToBadUniform()); if(resd->IsBadChannel(i)==kFALSE){ resd->SetBadChannel(count,i); count++; nSeToBad++; } } remainingBad = nBadUp; while (remainingBad>0) { Int_t nBadChannel; if (remainingBad<5) { nBadChannel = remainingBad; } else { Int_t width = remainingBad-(5-1); if (width>4) width = 4; nBadChannel = 5 + Int_t(width*gran->Uniform()); } Int_t iChannelPos = Int_t( (4*64-nBadChannel)*gran->Uniform() ); // Int_t iChip = iChannelPos/64; // Int_t iChan = iChannelPos - iChip*64; if(resd->IsBadChannel(iChannelPos)) continue; Int_t *clus = new Int_t[nBadChannel]; Int_t ich = iChannelPos; for(Int_t i=0;iIsBadChannel(clus[i])) break; resd->SetBadChannel(count,clus[i]); count++; } remainingBad -= nBadChannel; delete [] clus; } nSeToBad = 0; for (Int_t i=0; i<4; i++){ for(Int_t j=0;j<64;j++){ if (resd->Gain(1,i,j)<0.0001) nSeToBad++; } } while (nSeToBadUniform()); if(resd->IsBadChannel(i+256)==kFALSE){ resd->SetBadChannel(count,i+256); count++; nSeToBad++; } } //baselines /* for(Int_t nan=0;nan<512;nan++){ Int_t baseline = resd->GetBaseline(0); Double_t noise = resd->GetNoiseAfterElectronics(0); resd->SetBaseline(nan,baseline+(Int_t)gran->Gaus(0,4)); //baseline di def. e' 20, ho messo variazione di 4 resd->SetNoiseAfterElectronics(nan,noise+(Double_t)gran->Gaus(0,0.6)); } */ if(mod==88) resd->SetDead(); respSDD.Add(resd); printf("Added module %d\n",mod); } FILE* out = fopen("deadchannels.dat","w"); for(Int_t i=0;i<260;i++){ fprintf(out,"MODULE=%d\n",i); AliITSCalibrationSDD* cl = (AliITSCalibrationSDD*)respSDD.At(i); Int_t ndead=cl->GetDeadChannels(); fprintf(out,"n %d\n",ndead); for(Int_t n=0;nGetBadChannel(n)); } } fclose(out); AliCDBManager::Instance()->GetDefaultStorage()->Put(&respSDD, idCalSDD, md1); AliCDBManager::Instance()->GetDefaultStorage()->Put(rd, idRespSDD, md2); }