Adding a naive back2back jet-like structure for the toy model

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

diff --git a/T0/AliT0Calibrator.cxx b/T0/AliT0Calibrator.cxx
index 655a4aa2b58cc245039c664aa5624ae57ee4f53d..f5b4966631a891d38cd3da3058859656bccd24ce 100644 (file)
--- a/T0/AliT0Calibrator.cxx
+++ b/T0/AliT0Calibrator.cxx
@@ -1,3 +1,4 @@
+
 /**************************************************************************
  * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
  *                                                                        *
@@ -29,56 +30,49 @@
 #include "AliLog.h"
 #include "AliT0Parameters.h"
 #include "AliT0Calibrator.h"
+#include "AliT0Reconstructor.h"
+#include "AliT0RecoParam.h"
+#include "AliCDBEntry.h"
+#include "AliCDBManager.h"
+
 #include <TGraph.h>
 #include <TH1F.h>
 #include <TMath.h>
-#include "iostream.h"
+#include <Riostream.h>
 
 ClassImp(AliT0Calibrator)
 
 //____________________________________________________________________
   AliT0Calibrator::AliT0Calibrator():TNamed(),
                                     fChannelWidth(0),  
-                                    fWalk(0)
+                                    fWalk(0),
+                                    fEqualized(0)
                                     
 {
   //constructor
-
+   printf(" AliT0Calibrator ::: AliT0RecoParam GetEq() %i\n", fEqualized);
    AliT0Parameters* param = AliT0Parameters::Instance();
    param->Init();
-  
+   //slewing correcion and equalizing channels
+
   fChannelWidth = param->GetChannelWidth() ;  
-  Double_t *grY ;
-  Double_t *grX ;
-  Int_t index[2500];
-  Bool_t down=true;
+   //   Double_t *grX ; 
   for (Int_t i=0; i<24; i++){
     fMaxValue[i]=0;
     fTimeDelayCFD[i] = Int_t (param->GetTimeDelayCFD(i));
-       
-    //  TGraph* fu = param ->GetWalk(i);
-     TGraph* fu  = param ->GetAmpLEDRec(i);
-    if(fu) {
-      Int_t np=fu->GetN();
-      if(np>0) {
-       grY = fu->GetY();
-       grX = fu->GetX();
-       TMath::Sort(np, grY, index,down);
-       fMaxValue[i]=Int_t(grY[index[0]]);
-       fWalk.AddAtAndExpand(fu,i);
-      }
-    }
+     TGraph* fu = param ->GetWalk(i);
+    // fWalk.AddAtAndExpand(fu,i);
+    //TGraph* fu = param ->GetAmpLEDRec(i);
+    fWalk.AddAtAndExpand(fu,i);
   }
-  //  delete [] grY;
-  //  delete [] grX;
   
-  //
 }
 //_____________________________________________________________________________
 
 AliT0Calibrator::AliT0Calibrator(const AliT0Calibrator &r): TNamed(),
-  fChannelWidth(0),  
-  fWalk(0)
+                                                           fChannelWidth(0),  
+                                                           fWalk(0), 
+                                                           fEqualized(0)
 
 {
   //
@@ -102,36 +96,32 @@ AliT0Calibrator &AliT0Calibrator::operator=(const AliT0Calibrator &r)
 }
 
 
-
 //____________________________________________________________________
-Int_t  AliT0Calibrator::WalkCorrection(Int_t ipmt, Int_t qt, Int_t time, TString option) 
-{
-  //slewing correcion and equalizing channels
 
-  Int_t walk=0;
+Int_t  AliT0Calibrator::WalkCorrection(Int_t refAmp,  Int_t ipmt, Int_t qt, Int_t time) 
 
+{
+  //
+  // referemce amplitude for walk correction now read from RecoParam
+
+   Int_t walk=0;
   Int_t timeEq=0, timeWalk=0;  
-  //  TGraph* fu1 = param ->GetWalk(ipmt);
-  //   TGraph* fu1  = param ->GetAmpLEDRec(ipmt);
-   TGraph *fu1=(TGraph*) fWalk.At(ipmt);
+  TGraph *fu1=(TGraph*) fWalk.At(ipmt);
   if(fu1 && fu1->GetN()>0) {
-    walk=Int_t(fu1->Eval(Double_t(qt)));
-  }
-  if (option == "pdc") {
-    timeWalk = time + Int_t((fMaxValue[ipmt]-walk)/fChannelWidth) ;
-    timeEq= timeWalk - (fTimeDelayCFD[ipmt]-fTimeDelayCFD[0]);
-     AliDebug(10,Form(" ipmt %i time before %i timeWalk %i ,  qt %i timeEq %i \n ",
-                ipmt, time,timeWalk, qt, timeEq ));
+    walk = Int_t(fu1->Eval(Double_t(qt)));
   }
-  if (option == "cosmic") {
-    timeWalk = time + Int_t((fMaxValue[ipmt]-walk)) ;
-    if(walk <1  )  timeWalk = time  ;
-     timeEq= timeWalk - fTimeDelayCFD[ipmt];
-     AliDebug(10,Form(" ipmt %i time before %i timeWalk %i ,  qt %i timeEq %i \n ",
-                ipmt, time,timeWalk, qt, timeEq ));
-  }
- 
-  return timeEq;
+  
+  if (fEqualized == 0)
+    timeEq= time - fTimeDelayCFD[ipmt]-walk;
+  else 
+    timeEq = time - walk -  refAmp;
+
+  //   printf(" ipmt %i time before %i timeWalk %i , walk %i  qt %i fTimeDelayCFD[ipmt] %i timeEq %i \n ",
+  //    ipmt, time,timeWalk, walk, qt,fTimeDelayCFD[ipmt], timeEq );
+     AliDebug(2,Form(" fEqualized %i ipmt %i refAmp %i time before %i timeWalk %i , walk %i  qt %i timeEq %i, diff %i \n ",
+                    fEqualized,   ipmt, refAmp, time,timeWalk, walk, qt, timeEq , fTimeDelayCFD[ipmt]));
+  
+   return timeEq;
 }