SDD DA updated: baseline equalized to 20 for each hybrid. Added Carlos thresholds...

[u/mrichter/AliRoot.git] / ITS / AliITSPreprocessorSDD.cxx

diff --git a/ITS/AliITSPreprocessorSDD.cxx b/ITS/AliITSPreprocessorSDD.cxx
index 56ea4d3bcb0b95068e2e9820af089bfa8f9f8f8a..aa15a56c4aa101f1ba198405fffe91e0fe69f5c1 100644 (file)
--- a/ITS/AliITSPreprocessorSDD.cxx
+++ b/ITS/AliITSPreprocessorSDD.cxx
@@ -1,16 +1,35 @@
-/////////////////////////////////////////
-// Class for SDD digits preprocessing  //
-//                                     //
-//                                     //
-////////////////////////////////////////
+/**************************************************************************
+ * Copyright(c) 2007-2009, ALICE Experiment at CERN, All rights reserved. *
+ *                                                                        *
+ * Author: The ALICE Off-line Project.                                    *
+ * Contributors are mentioned in the code where appropriate.              *
+ *                                                                        *
+ * Permission to use, copy, modify and distribute this software and its   *
+ * documentation strictly for non-commercial purposes is hereby granted   *
+ * without fee, provided that the above copyright notice appears in all   *
+ * copies and that both the copyright notice and this permission notice   *
+ * appear in the supporting documentation. The authors make no claims     *
+ * about the suitability of this software for any purpose. It is          *
+ * provided "as is" without express or implied warranty.                  *
+ **************************************************************************/

 
 

-/* $Id$ */
+/* $Id:$ */
+
+///////////////////////////////////////////////////////////////////
+//                                                               //
+// Implementation of the class for SDD preprocessing             //
+// Origin: E.Crescio, Torino, crescio@to.infn.it                 //
+//         F.Prino, Torino, prino@to.infn.it                     //
+//                                                               //
+///////////////////////////////////////////////////////////////////

 
 #include "AliITSPreprocessorSDD.h"
 #include "AliITSCalibrationSDD.h"
 #include "AliITSDriftSpeedSDD.h"
 #include "AliITSDriftSpeedArraySDD.h"
 
 #include "AliITSPreprocessorSDD.h"
 #include "AliITSCalibrationSDD.h"
 #include "AliITSDriftSpeedSDD.h"
 #include "AliITSDriftSpeedArraySDD.h"

+#include "AliITSDCSAnalyzerSDD.h"

 #include "AliShuttleInterface.h"
 #include "AliShuttleInterface.h"

+#include "AliCDBEntry.h"

 #include "AliCDBMetaData.h"
 #include "TObjArray.h"
 #include "AliLog.h"
 #include "AliCDBMetaData.h"
 #include "TObjArray.h"
 #include "AliLog.h"


@@ -18,127 +37,170 @@
 #include <TSystem.h>
 #include <TList.h>
 
 #include <TSystem.h>
 #include <TList.h>
 

-const Int_t AliITSPreprocessorSDD::fgkNumberOfSDD = 260;
-const Int_t AliITSPreprocessorSDD::fgkNumberOfChannels = 512;

 const TString AliITSPreprocessorSDD::fgkNameHistoPedestals = "hpedestal";
 const TString AliITSPreprocessorSDD::fgkNameHistoNoise = "hnoise";
 ClassImp(AliITSPreprocessorSDD)
 
 
 const TString AliITSPreprocessorSDD::fgkNameHistoPedestals = "hpedestal";
 const TString AliITSPreprocessorSDD::fgkNameHistoNoise = "hnoise";
 ClassImp(AliITSPreprocessorSDD)
 
 

-UInt_t AliITSPreprocessorSDD::Process(TMap*/* dcsAliasMap*/){
+//______________________________________________________________________
+UInt_t AliITSPreprocessorSDD::Process(TMap* dcsAliasMap){

 
 

-  //preprocessing. 
+  // Get DDL map from OCDB
+  AliCDBEntry* entry = GetFromOCDB("Calib", "DDLMapSDD");
+  if(!entry){
+    Log("DDL map file not found in OCDB.");  
+    return 2;
+  }
+  AliITSDDLModuleMapSDD* ddlmap = (AliITSDDLModuleMapSDD*)entry->GetObject();
+  if(!ddlmap){ 
+    Log("AliITSDDLModuleMapSDD object not in file.");
+    return 2;
+  }

 
 

+  //preprocessing

   TString runType = GetRunType();
   Bool_t retcode;
   TString runType = GetRunType();
   Bool_t retcode;

+
+  if (runType == "PULSER"){
+    retcode=ProcessPulser(ddlmap);
+  }else if(runType== "INJECTOR"){
+    retcode=ProcessInjector(ddlmap);
+  }else{
+    // do nothing for other run types
+    retcode=1;
+  }
+  if(retcode){
+    Bool_t retcodedcs =ProcessDCSDataPoints(dcsAliasMap);
+    if(retcodedcs) return 0;
+  }
+  return 1;
+}
+//______________________________________________________________________
+Bool_t AliITSPreprocessorSDD::ProcessPulser(AliITSDDLModuleMapSDD* ddlmap){
+  // Process FXS files from PULSER run (baseline, noise, gain)
+  TObjArray calSDD(kNumberOfSDD);
+  calSDD.SetOwner(kFALSE);

   Char_t command[100];
   Char_t inpFileName[100];
   Char_t command[100];
   Char_t inpFileName[100];

-  AliCDBMetaData *md1= new AliCDBMetaData();
-  md1->SetResponsible("Elisabetta Crescio, Francesco Prino");
-  md1->SetBeamPeriod(0);
-  md1->SetAliRootVersion("head 30 November 2007"); //root version
-  md1->SetComment("This is a test");
-
-  if (runType == "PULSER_RUN"){
-    TObjArray calSDD(fgkNumberOfSDD);
-    calSDD.SetOwner(kFALSE);
-    Float_t baseline,rawnoise,cmn,corn,gain;
-    Int_t isgoodan,i,im,is,isgoodmod;
-    Int_t numOfBadChannels[fgkNumberOfSDD];
+  Float_t baseline,rawnoise,cmn,corn,gain;
+  Int_t isgoodan,i,im,is,isgoodmod,basmin,basoff;
+  Int_t th,tl;
+  Int_t numOfBadChannels[kNumberOfSDD];

   
   

-    TList* sourceList = GetFileSources(kDAQ, "SDD_Calib");
-    if (!sourceList){ 
-      Log("Error: no sources found for SDD_Calib");
-      return 2;
-    }
+  TList* sourceList = GetFileSources(kDAQ, "SDD_Calib");
+  if (!sourceList){ 
+    Log("Error: no sources found for SDD_Calib");
+    return 2;
+  }

 
 

-    Int_t ind = 0;
-    while (sourceList->At(ind)!=NULL) {
-      TObjString* tarId = (TObjString*) sourceList->At(ind);
-      TString tarName = GetFile(kDAQ, "SDD_Calib", tarId->GetString().Data());
-      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++;
+  Int_t ind = 0;
+  while (sourceList->At(ind)!=NULL) {
+    TObjString* tarId = (TObjString*) sourceList->At(ind);
+    TString tarName = GetFile(kDAQ, "SDD_Calib", tarId->GetString().Data());
+    if(tarName.Length()==0){
+      Log(Form("Baseline tar file from source %d not found.",ind));
+      return 2;

     }
     }

-    delete sourceList;
-
-    for(Int_t imod=0;imod<fgkNumberOfSDD;imod++){
+    sprintf(command,"tar -xf %s",tarName.Data());
+    gSystem->Exec(command);
+    ind++;
+  }
+  delete sourceList;
+  
+  for(Int_t iddl=0;iddl<kNumberOfDDL;iddl++){
+    for(Int_t imod=0;imod<kModulesPerDDL;imod++){
+      Int_t modID=ddlmap->GetModuleNumber(iddl,imod);
+      if(modID==-1) continue;
+      modID-=240; // to have SDD modules numbering from 0 to 260

       AliITSCalibrationSDD *cal = new AliITSCalibrationSDD("simulated");
       AliITSCalibrationSDD *cal = new AliITSCalibrationSDD("simulated");

-      numOfBadChannels[imod]=0;
-      Int_t badch[fgkNumberOfChannels];
+      numOfBadChannels[modID]=0;
+      Int_t badch[kNumberOfChannels];

       for(Int_t isid=0;isid<=1;isid++){
       for(Int_t isid=0;isid<=1;isid++){

-       sprintf(inpFileName,"./SDDbase_mod%03d_sid%d.data",imod,isid);
-       FILE* basFil = fopen(inpFileName,"read");
-       if (basFil == 0) {
-        Log(Form("File %s not found.",inpFileName));
-        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(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){ 
-          Int_t ibad=numOfBadChannels[imod];
-          numOfBadChannels[imod]++;
-          badch[ibad]=ich;
-        }
-        cal->SetBaseline(ich,baseline);
-        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(basFil);
+       sprintf(inpFileName,"./SDDbase_ddl%02dc%02d_sid%d.data",iddl,imod,isid);
+       FILE* basFil = fopen(inpFileName,"read");
+       if (basFil == 0) {
+         Log(Form("File %s not found.",inpFileName));
+         cal->SetBad();
+         continue;
+       }
+       fscanf(basFil,"%d %d %d\n",&im,&is,&isgoodmod);
+       fscanf(basFil,"%d\n",&th);
+       fscanf(basFil,"%d\n",&tl);
+       if(!isgoodmod) cal->SetBad();
+       for(Int_t ian=0;ian<(kNumberOfChannels/2);ian++){
+         fscanf(basFil,"%d %d %f %d %d %f %f %f %f\n",&i,&isgoodan,&baseline,&basmin,&basoff,&rawnoise,&cmn,&corn,&gain);
+         Int_t ich=ian;
+         if(isid==1) ich+=256;
+         if(!isgoodan){ 
+           Int_t ibad=numOfBadChannels[modID];
+           badch[ibad]=ich;
+           numOfBadChannels[modID]++;
+         }
+         cal->SetBaseline(ich,baseline);
+         cal->SetNoiseAfterElectronics(ich,rawnoise);
+         cal->SetGain(ich,gain);
+       }
+       cal->SetDeadChannels(numOfBadChannels[modID]);
+       for(Int_t ibad=0;ibad<numOfBadChannels[modID];ibad++){
+         cal->SetBadChannel(ibad,badch[ibad]);
+       }
+       fclose(basFil);

       }
       }

-      calSDD.Add(cal);
+      calSDD.AddAt(cal,modID);

     }
     }

-    md1->SetObjectClassName("AliITSCalibration");
-    retcode = Store("Calib","CalibSDD",&calSDD,md1, 0, kTRUE);
-  }else if(runType == "PHYSICS"){
-
-    TObjArray vdrift(2*fgkNumberOfSDD);
-    vdrift.SetOwner(kFALSE);
-    Int_t evNumb,polDeg; 
-    UInt_t timeStamp;
-    Float_t param[4];
+  }
+  AliCDBMetaData *md= new AliCDBMetaData();
+  md->SetResponsible("Francesco Prino");
+  md->SetBeamPeriod(0);
+  md->SetComment("AliITSCalibrationSDD from PEDESTAL+PULSER runs");
+  Bool_t retCode = Store("Calib","CalibSDD",&calSDD,md, 0, kTRUE);
+  return retCode;
+}
+//______________________________________________________________________
+Bool_t AliITSPreprocessorSDD::ProcessInjector(AliITSDDLModuleMapSDD* ddlmap){
+  // Process FXS files from injector events (INJECTOR or PHYSICS runs)
+  TObjArray vdrift(2*kNumberOfSDD);
+  vdrift.SetOwner(kFALSE);
+  Char_t command[100];
+  Char_t inpFileName[100];
+  Int_t evNumb,polDeg; 
+  UInt_t timeStamp;
+  Float_t param[4];

 
 

-    TList* sourceList = GetFileSources(kDAQ, "SDD_Injec");
-    if (!sourceList){ 
-      Log("Error: no sources found for SDD_Injec");
+  TList* sourceList = GetFileSources(kDAQ, "SDD_Injec");
+  if (!sourceList){ 
+    Log("Error: no sources found for SDD_Injec");
+    return 2;
+  }
+  Int_t ind = 0;
+  while (sourceList->At(ind)!=NULL) {
+    TObjString* tarId = (TObjString*) sourceList->At(ind);
+    TString tarName = GetFile(kDAQ, "SDD_Injec", tarId->GetString().Data());
+    if(tarName.Length()==0){
+      Log(Form("Injector tar file from source %d not found.",ind));

       return 2;
     }
       return 2;
     }

-    Int_t ind = 0;
-    while (sourceList->At(ind)!=NULL) {
-      TObjString* tarId = (TObjString*) sourceList->At(ind);
-      TString tarName = GetFile(kDAQ, "SDD_Injec", tarId->GetString().Data());
-      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++){
+    sprintf(command,"tar -xf %s",tarName.Data());
+    gSystem->Exec(command);
+    ind++;
+  }
+  delete sourceList;
+  
+  for(Int_t iddl=0;iddl<kNumberOfDDL;iddl++){
+    for(Int_t imod=0;imod<kModulesPerDDL;imod++){
+      Int_t modID=ddlmap->GetModuleNumber(iddl,imod);
+      if(modID==-1) continue;
+      modID-=240; // to have SDD modules numbering from 0 to 260

       for(Int_t isid=0;isid<=1;isid++){
        AliITSDriftSpeedArraySDD *arr=new AliITSDriftSpeedArraySDD();
       for(Int_t isid=0;isid<=1;isid++){
        AliITSDriftSpeedArraySDD *arr=new AliITSDriftSpeedArraySDD();

-       sprintf(inpFileName,"./SDDinj_mod%03d_sid%d.data",imod,isid);
+       sprintf(inpFileName,"./SDDinj_ddl%02dc%02d_sid%d.data",iddl,imod,isid);

        FILE* injFil = fopen(inpFileName,"read");
        if (injFil == 0) {
          Log(Form("File %s not found.",inpFileName));
        FILE* injFil = fopen(inpFileName,"read");
        if (injFil == 0) {
          Log(Form("File %s not found.",inpFileName));

-         return 2;
+         AliITSDriftSpeedSDD *dsp=new AliITSDriftSpeedSDD();
+         arr->AddDriftSpeed(dsp);
+         vdrift.AddAt(arr,2*modID+isid);
+         continue;

        }
        fscanf(injFil,"%d",&polDeg);
        while (!feof(injFil)){
        }
        fscanf(injFil,"%d",&polDeg);
        while (!feof(injFil)){


@@ -148,15 +210,32 @@ UInt_t AliITSPreprocessorSDD::Process(TMap*/* dcsAliasMap*/){
          AliITSDriftSpeedSDD *dsp=new AliITSDriftSpeedSDD(evNumb,timeStamp,polDeg,param);
          arr->AddDriftSpeed(dsp);
        }
          AliITSDriftSpeedSDD *dsp=new AliITSDriftSpeedSDD(evNumb,timeStamp,polDeg,param);
          arr->AddDriftSpeed(dsp);
        }

-       vdrift.Add(arr);
+       vdrift.AddAt(arr,2*modID+isid);

       }
     }
       }
     }

-    md1->SetObjectClassName("AliITSDriftSpeedArraySDD");
-    retcode = Store("Calib","DriftSpeedSDD",&vdrift,md1,0, kTRUE);    
-  }else{
-    // do nothing for other run types
-    retcode=1;

   }
   }

-  if(retcode) return 0;
-  else return 1;
+  AliCDBMetaData *md= new AliCDBMetaData();
+  md->SetResponsible("Francesco Prino");
+  md->SetBeamPeriod(0);
+  md->SetComment("AliITSDriftSpeedSDD from injector events");
+  Bool_t retCode = Store("Calib","DriftSpeedSDD",&vdrift,md,0, kTRUE);
+  return retCode;
+}
+//______________________________________________________________________
+Bool_t AliITSPreprocessorSDD::ProcessDCSDataPoints(TMap* dcsAliasMap){
+  // Process DCS data
+  AliITSDCSAnalyzerSDD *dcs=new AliITSDCSAnalyzerSDD();
+  dcs->AnalyzeData(dcsAliasMap);
+  TObjArray refDCS(kNumberOfSDD);
+  refDCS.SetOwner(kFALSE);
+  for(Int_t imod=0;imod<kNumberOfSDD;imod++){
+    AliITSDCSDataSDD *dcsdata=dcs->GetDCSData(imod);
+    refDCS.Add(dcsdata);
+  }    
+  AliCDBMetaData *md= new AliCDBMetaData();
+  md->SetResponsible("Francesco Prino");
+  md->SetBeamPeriod(0);
+  md->SetComment("AliITSDCSDataSDD objects from DCS DB");
+  Bool_t retCode = StoreReferenceData("DCS","DataSDD",&refDCS,md);
+  return retCode;

 }
 }