New raw-reader class which deals with events taken from shared memory via the DATE...

[u/mrichter/AliRoot.git] / T0 / AliT0CalibTimeEq.cxx

diff --git a/T0/AliT0CalibTimeEq.cxx b/T0/AliT0CalibTimeEq.cxx
index 128af62bd2e81d2de3ce0746542e06e54220d7c7..e37a82a6ac309534ece2b214679f3a726cd03969 100644 (file)
--- a/T0/AliT0CalibTimeEq.cxx
+++ b/T0/AliT0CalibTimeEq.cxx
@@ -17,43 +17,33 @@
 
 ///////////////////////////////////////////////////////////////////////////////
 //                                                                           //
-// class for T0 calibration                       TM-AC-AM_6-02-2006         //
+// class for T0 calibration                       TM-AC-AM_6-02-2006  
+// equalize time shift for each time CFD channel
 //                                                                           //
 ///////////////////////////////////////////////////////////////////////////////
 
 #include "AliT0CalibTimeEq.h"
-#include "AliT0LookUpValue.h"
-#include "AliLog.h"
-#include "AliRun.h"
 
-#include <TCanvas.h>
 #include <TFile.h>
 #include <TMath.h>
 #include <TF1.h>
-#include <TSystem.h>
-#include <Riostream.h>
 #include <TSpectrum.h>
-#include <TVirtualFitter.h>
 #include <TProfile.h>
-
-#include <Riostream.h>
-#include <string>
+#include <iostream>
 
 ClassImp(AliT0CalibTimeEq)
 
 //________________________________________________________________
-  AliT0CalibTimeEq::AliT0CalibTimeEq():   TNamed(),
-                                      fTimeDelayTVD(0),
-                                      fMeanT0(0) 
+  AliT0CalibTimeEq::AliT0CalibTimeEq():TNamed()
 {
   //
 }
 
 //________________________________________________________________
-AliT0CalibTimeEq::AliT0CalibTimeEq(const char* name):TNamed(),
-                                      fTimeDelayTVD(0),
-                                      fMeanT0(0)
+AliT0CalibTimeEq::AliT0CalibTimeEq(const char* name):TNamed()
 {
+  //constructor
+
   TString namst = "Calib_";
   namst += name;
   SetName(namst.Data());
@@ -61,10 +51,7 @@ AliT0CalibTimeEq::AliT0CalibTimeEq(const char* name):TNamed(),
 }
 
 //________________________________________________________________
-AliT0CalibTimeEq::AliT0CalibTimeEq(const AliT0CalibTimeEq& calibda) :
-  TNamed(calibda),             
-  fTimeDelayTVD(0),
-  fMeanT0(0)
+AliT0CalibTimeEq::AliT0CalibTimeEq(const AliT0CalibTimeEq& calibda):TNamed(calibda)            
 {
 // copy constructor
   SetName(calibda.GetName());
@@ -87,10 +74,13 @@ AliT0CalibTimeEq &AliT0CalibTimeEq::operator =(const AliT0CalibTimeEq& calibda)
 AliT0CalibTimeEq::~AliT0CalibTimeEq()
 {
   //
+  // destrictor
 }
 //________________________________________________________________
 void AliT0CalibTimeEq::Reset()
 {
+  //reset values
+
   memset(fCFDvalue,0,120*sizeof(Float_t));
   memset(fTimeEq,1,24*sizeof(Float_t));
 }
@@ -99,6 +89,7 @@ void AliT0CalibTimeEq::Reset()
 //________________________________________________________________
 void  AliT0CalibTimeEq::Print(Option_t*) const
 {
+  // print time values
 
   printf("\n   ----    PM Arrays       ----\n\n");
   printf(" Time delay CFD \n");
@@ -107,59 +98,74 @@ void  AliT0CalibTimeEq::Print(Option_t*) const
 
 
 //________________________________________________________________
-void AliT0CalibTimeEq::ComputeOnlineParams(char* name1, char* name2, char* canv, Int_t npeaks, Double_t sigma, const char* filePhys)
+void AliT0CalibTimeEq::ComputeOnlineParams(const char* filePhys)
 {
+  // compute online equalized time
+//  Int_t npeaks = 20;
+//  Double_t sigma = 4.;
+
   TFile *gFile = TFile::Open(filePhys);
-  gSystem->Load("libSpectrum");
-  npeaks = 20;
-  sigma=3.;
-  Bool_t down=false;
-  Int_t index[20];
-  TCanvas *c1 = new TCanvas(canv, canv,0,48,1280,951);
-  c1->Divide(4,3);
+//  Bool_t down=false;
+//  Int_t index[20];
   Char_t buf1[15];
-  Char_t temp[10];
-  Float_t p[12][3]={0.,0.,0.};
-  for (Int_t i=12; i<24; i++)
+//  Char_t temp[10];
+//  Float_t p[24][3]={0.,0.,0.};
+  for (Int_t i=0; i<24; i++)
+  {
+    sprintf(buf1,"CFD1-CFD%d",i+1);
+    TH1F *cfd = (TH1F*) gFile->Get(buf1);
+    //    printf(" i = %d buf1 = %s\n", i, buf1);
+    Double_t mean=cfd->GetMean();
+    SetTimeEq(i,mean);
+    delete cfd;
+  }
+    /*
+
+    TSpectrum *s = new TSpectrum(2*npeaks,1.);
+    printf(" buf1 = %s cfd = %x\n", buf1, cfd);
+    Int_t nfound = s->Search(cfd,sigma,"goff",0.2);
+    printf(" nfound = %d\n", nfound);
+    if(nfound!=0)
     {
-      c1->cd(i+1);
-      sprintf(buf1,name1);
-      sprintf(temp,"%i",i+1);
-      strcat (buf1,temp);
-      strcat (buf1,name2);
-      TH1F *cfd = (TH1F*) gFile->Get(buf1);
-      TSpectrum *s = new TSpectrum(2*npeaks,1.);
-      Int_t nfound = s->Search(cfd,sigma," ",0.2);
-      if(nfound!=0){
-        Float_t *xpeak = s->GetPositionX();
-        TMath::Sort(nfound, xpeak, index,down);
-        Float_t xp = xpeak[index[0]];
-        Float_t hmax = xp+3*sigma;
-        Float_t hmin = xp-3*sigma;
-        cfd->GetXaxis()->SetRange((Int_t)hmin-20,(Int_t)hmax+20);
-        TF1 *g1 = new TF1("g1", "gaus", hmin, hmax);
-              cfd->Fit("g1","IR");
-
-        for(Int_t j =0; j<3; j++){
-                p[i][j] = g1->GetParameter(j);
-                SetCFDvalue(i, j, p[i][j]);
-        }
-        SetCFDvalue(i, 3, hmin);
-        SetCFDvalue(i, 4, hmax);
-        SetTimeEq(i,p[i][2]);
-        // if(p[i][0]==0)
-        // cfd->Draw();
-
-        // cfd->Draw();
-      } // else
-        // cfd->Draw();
-    }
-  // TFile *onl = new TFile("onl.root","RECREATE");
-  // this->Write("Values");
-  // onl->Close();
-  // delete onl;
-  gFile->Close();
-  delete gFile;
+      Float_t *xpeak = s->GetPositionX();
+      TMath::Sort(nfound, xpeak, index,down);
+      Float_t xp = xpeak[index[0]];
+      Float_t hmax = xp+3*sigma;
+      Float_t hmin = xp-3*sigma;
+      cfd->GetXaxis()->SetRange((Int_t)hmin-20,(Int_t)hmax+20);
+      TF1 *g1 = new TF1("g1", "gaus", hmin, hmax);
+      cfd->Fit("g1","IRQN");
+      for(Int_t j =0; j<3; j++)
+      {
+        p[i][j] = g1->GetParameter(j);
+       //      SetCFDvalue(i, j, p[i][j]);
+      }
+        SetTimeEq(i,(p[i][2]));
+       cout<<" mean "<<p[i][2]<<endl;
+    */
+      /*
+      SetCFDvalue(i, 3, hmin);
+      SetCFDvalue(i, 4, hmax);
+
+      if (i<12)
+      {
+        SetTimeEq(i,(p[i][2]-p[0][2]));
+      }
+      else
+      {
+       SetTimeEq(i,(p[i][2]-p[12][2]));
+      }        
+      
+    } 
+  }
+      */  
+   gFile->Close();
+   delete gFile;
+   printf("\n\n");
+   for(int j=0;j<24;j++)
+   {
+                printf("fTimeEq[%d]=%f\n",j,fTimeEq[j]);
+   }
 }