]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - ITS/AliITSPreprocessorSDD.cxx
documentation fixes
[u/mrichter/AliRoot.git] / ITS / AliITSPreprocessorSDD.cxx
index e12c19c41daccfe88f133ae64b7a98ea78ceeb57..8cdefee4258223f53286eff276d1cb5631b44d79 100644 (file)
@@ -4,12 +4,17 @@
 //                                     //
 ////////////////////////////////////////
 
+/* $Id$ */
+
 #include "AliITSPreprocessorSDD.h"
 #include "AliITSCalibrationSDD.h"
 #include "AliShuttleInterface.h"
 #include "AliCDBMetaData.h"
 #include "TObjArray.h"
 #include "AliLog.h"
+#include <TObjString.h>
+#include <TSystem.h>
+#include <TList.h>
 
 const Int_t AliITSPreprocessorSDD::fgkNumberOfSDD = 260;
 const Int_t AliITSPreprocessorSDD::fgkNumberOfChannels = 512;
@@ -22,38 +27,74 @@ UInt_t AliITSPreprocessorSDD::Process(TMap*/* dcsAliasMap*/){
 
   //preprocessing. 
 
-
-  AliCDBMetaData *md1= new AliCDBMetaData(); // metaData describing the object
-  md1->SetObjectClassName("AliITSCalibration");
-  md1->SetResponsible("Elisabetta Crescio, Francesco Prino");
-  md1->SetBeamPeriod(0);
-  md1->SetAliRootVersion("head 5 April 2007"); //root version
-  md1->SetComment("This is a test");
-
   TObjArray respSDD(fgkNumberOfSDD);
   respSDD.SetOwner(kFALSE);
-  
-  Char_t filid[20];
-  Float_t baseline,rawnoise,cmn,gain;
+  Float_t baseline,rawnoise,cmn,corn,gain;
   Int_t isgoodan,i,im,is,isgoodmod;
   Int_t numOfBadChannels[fgkNumberOfSDD];
+  //TString pwd = gSystem->pwd();
+//  const Char_t* tempDir=AliShuttleInterface::GetShuttleTempDir();
+  
+  TList* sourceList = GetFileSources(kDAQ, "SDD_Calib");
+  if (!sourceList){ 
+    Log("Error: no sources found for SDD_Calib");
+    return 2;
+  }
+
+  Int_t ind = 0;
+  Char_t command[100];
+  while (sourceList->At(ind)!=NULL) {
+    TObjString* tarId = (TObjString*) sourceList->At(ind);
+    TString tarName = GetFile(kDAQ, "SDD_Calib", tarId->GetString().Data());
+//    gSystem->cd(tempDir);
+    if(tarName.Length()==0){
+      Log(Form("Baseline tar file from source %d not found.",ind));
+      return 2;
+    }
+    sprintf(command,"tar -xf %s",tarName.Data());
+    gSystem->Exec(command);
+    ind++;
+  }
+   
+  sourceList = GetFileSources(kDAQ, "SDD_Injec");
+  if (!sourceList){ 
+    Log("Error: no sources found for SDD_Injec");
+    return 2;
+  }
+
+  ind = 0;
+  while (sourceList->At(ind)!=NULL) {
+   TObjString* tarId = (TObjString*) sourceList->At(ind);
+   TString tarName = GetFile(kDAQ, "SDD_Injec", tarId->GetString().Data());
+//    gSystem->cd(tempDir);
+    if(tarName.Length()==0){
+      Log(Form("Injector tar file from source %d not found.",ind));
+      return 2;
+    }
+    sprintf(command,"tar -xf %s",tarName.Data());
+    gSystem->Exec(command);
+    ind++;
+  }
+  
+  delete sourceList;
+
   for(Int_t imod=0;imod<fgkNumberOfSDD;imod++){
     AliITSCalibrationSDD *cal = new AliITSCalibrationSDD("simulated");
     numOfBadChannels[imod]=0;
     Int_t badch[fgkNumberOfChannels];
     for(Int_t isid=0;isid<=1;isid++){
-      sprintf(filid,"DAQDAm%03ds%d",imod,isid);
-      const char* filenamed= GetFile(kDAQ,filid,"GDC");
-      FILE* filed = fopen(filenamed,"read");
-      if (filed == 0) {
-       AliWarning("File not found"); 
-       continue;
-      }
-      
-      fscanf(filed,"%d %d %d\n",&im,&is,&isgoodmod);
+      Char_t basFileName[100];
+//      sprintf(basFileName,"%s/SDDbase_mod%03d_sid%d.data",tempDir,imod,isid);
+      sprintf(basFileName,"./SDDbase_mod%03d_sid%d.data",imod,isid);
+      FILE* basFil = fopen(basFileName,"read");
+      if (basFil == 0) {
+       Log(Form("File %s not found.",basFileName));
+       return 2;
+      }      
+      fscanf(basFil,"%d %d %d\n",&im,&is,&isgoodmod);
       if(!isgoodmod) cal->SetDead();
       for(Int_t ian=0;ian<(fgkNumberOfChannels/2);ian++){
-       fscanf(filed,"%d %d %f %f %f %f\n",&i,&isgoodan,&baseline,&rawnoise,&cmn,&gain);
+       fscanf(basFil,"%d %d %f %f %f %f %f\n",&i,&isgoodan,&baseline,&rawnoise,&cmn,&corn,&gain);
        Int_t ich=ian;
        if(isid==1) ich+=256;
        if(!isgoodan){ 
@@ -62,20 +103,45 @@ UInt_t AliITSPreprocessorSDD::Process(TMap*/* dcsAliasMap*/){
          badch[ibad]=ich;
        }
        cal->SetBaseline(ich,baseline);
-       cal->SetNoiseAfterElectronics(ich,rawnoise);    
-       //      cal->SetGain(gain,isid,ian/4,ian%4);
+       cal->SetNoiseAfterElectronics(ich,rawnoise);
+       Int_t iChip=cal->GetChip(ich);
+       Int_t iChInChip=cal->GetChipChannel(ich);
+       cal->SetGain(gain,isid,iChip,iChInChip);
       }
       cal->SetDeadChannels(numOfBadChannels[imod]);
       for(Int_t ibad=0;ibad<numOfBadChannels[imod];ibad++){
        cal->SetBadChannel(ibad,badch[ibad]);
       }
-      fclose(filed);
+      fclose(basFil);
+
+      Char_t injFileName[100];
+      Int_t evNumb; 
+      UInt_t timeStamp;
+      Float_t param[4];
+      sprintf(injFileName,"./SDDinj_mod%03d_sid%d.data",imod,isid);
+      FILE* injFil = fopen(injFileName,"read");
+      if (injFil == 0) {
+       Log(Form("File %s not found.",basFileName));
+       return 2;
+      }      
+      while (!feof(injFil)){
+       fscanf(injFil,"%d %d",&evNumb,&timeStamp);
+       if(feof(injFil)) break;
+       for(Int_t ic=0;ic<4;ic++) fscanf(injFil,"%f",&param[ic]);
+       cal->SetDriftSpeedParam(isid,param);
+      }
     }
     respSDD.Add(cal);
   }
 
-  Bool_t retcode = Store("Calib","Data",&respSDD,md1);
+  AliCDBMetaData *md1= new AliCDBMetaData(); // metaData describing the object
+  md1->SetObjectClassName("AliITSCalibration");
+  md1->SetResponsible("Elisabetta Crescio, Francesco Prino");
+  md1->SetBeamPeriod(0);
+  md1->SetAliRootVersion("head 5 April 2007"); //root version
+  md1->SetComment("This is a test");
 
+  Bool_t retcode = Store("Calib","CalibSDD",&respSDD,md1, 0, kTRUE);
 
   if(retcode) return 0;
   else return 1;