]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - ITS/AliITSOnlineSDDTP.cxx
New class for common mode noise correction
[u/mrichter/AliRoot.git] / ITS / AliITSOnlineSDDTP.cxx
index 4b66f52b0b5691f8125d9f5013d01c22cfe697a7..f23f3c30c27175f8308eb00272a9e477c3c7696a 100644 (file)
@@ -13,6 +13,7 @@
  * provided "as is" without express or implied warranty.                  *
  **************************************************************************/
 #include "AliITSOnlineSDDTP.h"
+#include "AliLog.h"
 #include <TH2F.h>
 #include <TMath.h>
 
@@ -81,10 +82,10 @@ void AliITSOnlineSDDTP::AddEvent(TH2F* hrawd){
 void AliITSOnlineSDDTP::ReadBaselines(){
   // assume baselines and good anodes are taken from previous run
   Char_t basfilnam[100];
-  sprintf(basfilnam,"SDDbase_step1_mod%03d_sid%d.data",fModuleId,fSide);
+  sprintf(basfilnam,"SDDbase_step2_mod%03d_sid%d.data",fModuleId,fSide);
   FILE* basf=fopen(basfilnam,"r");
   if(basf==0){
-    printf("Baselinefile not present, Set all baselines to 50\n");
+    AliWarning("Baselinefile not present, Set all baselines to 50\n");
     for(Int_t ian=0;ian<fgkNAnodes;ian++){ 
       fBaseline[ian]=50.;
       fGoodAnode[ian]=1;
@@ -105,7 +106,6 @@ void AliITSOnlineSDDTP::ReadBaselines(){
 void AliITSOnlineSDDTP::ValidateAnodes(){
   Float_t meang,rmsg;
   StatGain(meang,rmsg);
-  printf("<gain>=%f,rms=%f\n",meang,rmsg);
   Float_t lowlim=meang-fNSigmaGain*rmsg;
   Float_t hilim=meang+fNSigmaGain*rmsg;
 
@@ -140,8 +140,12 @@ void AliITSOnlineSDDTP::StatGain(Float_t &mean, Float_t  &rms){
 void AliITSOnlineSDDTP::WriteToFXS(){
   //
   Char_t basfilnam[100];
-  sprintf(basfilnam,"SDDbase_step1_mod%03d_sid%d.data",fModuleId,fSide);
+  sprintf(basfilnam,"SDDbase_step2_mod%03d_sid%d.data",fModuleId,fSide);
   FILE* basf=fopen(basfilnam,"r");
+  if(basf==0){
+    AliWarning("Baseline file not present, launch baseline analysis first\n");
+    return;
+  }
   Int_t n,ok;
   Float_t base,rms,cmn,corrnoi;
   Float_t noise[fgkNAnodes],cmncoef[fgkNAnodes],corrnoise[fgkNAnodes];
@@ -152,7 +156,6 @@ void AliITSOnlineSDDTP::WriteToFXS(){
     corrnoise[ian]=corrnoi;
   }
   fclose(basf);
-  printf("Read All******************\n");
   Char_t outfilnam[100];
   sprintf(outfilnam,"SDDbase_mod%03d_sid%d.data",fModuleId,fSide);
   FILE* outf=fopen(outfilnam,"w");