1 /////////////////////////////////////////
2 // Class for SDD digits preprocessing //
5 ////////////////////////////////////////
7 #include "AliITSPreprocessorSDD.h"
8 #include "AliITSCalibrationSDD.h"
9 #include "AliShuttleInterface.h"
10 #include "AliCDBMetaData.h"
11 #include "TObjArray.h"
13 #include <TObjString.h>
17 const Int_t AliITSPreprocessorSDD::fgkNumberOfSDD = 260;
18 const Int_t AliITSPreprocessorSDD::fgkNumberOfChannels = 512;
19 const Int_t AliITSPreprocessorSDD::fgkNumberOfChannelsPerChip = 64;
20 const TString AliITSPreprocessorSDD::fgkNameHistoPedestals = "hpedestal";
21 const TString AliITSPreprocessorSDD::fgkNameHistoNoise = "hnoise";
22 ClassImp(AliITSPreprocessorSDD)
25 UInt_t AliITSPreprocessorSDD::Process(TMap*/* dcsAliasMap*/){
29 TObjArray respSDD(fgkNumberOfSDD);
30 respSDD.SetOwner(kFALSE);
31 Float_t baseline,rawnoise,cmn,gain;
32 Int_t isgoodan,i,im,is,isgoodmod;
33 Int_t numOfBadChannels[fgkNumberOfSDD];
34 //TString pwd = gSystem->pwd();
35 const Char_t* tempDir=AliShuttleInterface::GetShuttleTempDir();
37 TList* sourceList = GetFileSources(kDAQ, "SDD_Calib");
39 Log("Error: no sources found for SDD_Calib");
44 while (sourceList->At(ind)!=NULL) {
45 TObjString* tarId = (TObjString*) sourceList->At(ind);
46 TString tarName = GetFile(kDAQ, "SDD_Calib", tarId->GetString().Data());
49 sprintf(command,"tar -xf %s",tarName.Data());
50 gSystem->Exec(command);
56 for(Int_t imod=0;imod<fgkNumberOfSDD;imod++){
57 AliITSCalibrationSDD *cal = new AliITSCalibrationSDD("simulated");
58 numOfBadChannels[imod]=0;
59 Int_t badch[fgkNumberOfChannels];
60 for(Int_t isid=0;isid<=1;isid++){
61 Char_t basFileName[100];
62 sprintf(basFileName,"%s/SDDbase_mod%03d_sid%d.data",tempDir,imod,isid);
63 FILE* basFil = fopen(basFileName,"read");
65 Log(Form("File %s not found.",basFileName));
68 fscanf(basFil,"%d %d %d\n",&im,&is,&isgoodmod);
69 if(!isgoodmod) cal->SetDead();
70 for(Int_t ian=0;ian<(fgkNumberOfChannels/2);ian++){
71 fscanf(basFil,"%d %d %f %f %f %f\n",&i,&isgoodan,&baseline,&rawnoise,&cmn,&gain);
75 Int_t ibad=numOfBadChannels[imod];
76 numOfBadChannels[imod]++;
79 cal->SetBaseline(ich,baseline);
80 cal->SetNoiseAfterElectronics(ich,rawnoise);
81 Int_t iChip=ian/fgkNumberOfChannelsPerChip;
82 Int_t iChInChip=ian%fgkNumberOfChannelsPerChip;
83 cal->SetGain(gain,isid,iChip,iChInChip);
85 cal->SetDeadChannels(numOfBadChannels[imod]);
86 for(Int_t ibad=0;ibad<numOfBadChannels[imod];ibad++){
87 cal->SetBadChannel(ibad,badch[ibad]);
94 AliCDBMetaData *md1= new AliCDBMetaData(); // metaData describing the object
95 md1->SetObjectClassName("AliITSCalibration");
96 md1->SetResponsible("Elisabetta Crescio, Francesco Prino");
97 md1->SetBeamPeriod(0);
98 md1->SetAliRootVersion("head 5 April 2007"); //root version
99 md1->SetComment("This is a test");
101 Bool_t retcode = Store("Calib","CalibSDD",&respSDD,md1, 0, kTRUE);
103 if(retcode) return 0;