update

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

/**************************************************************************
 * Copyright(c) 1998-1999, 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$ */
/***********************************************************************
 *      this class doing calibration during reconstruction 
 *      2 steps:
 *      - equalizing channels
 *      - applying walk corrections
 *
 * Alla.Maevskaya@cern.ch
 *
 **********************************************************************/      
 
//#include <Riostream.h>

#include "AliLog.h"
#include "AliT0Parameters.h"
#include "AliT0Calibrator.h"
#include <TGraph.h>
#include <TH1F.h>
#include <TMath.h>
#include <Riostream.h>

ClassImp(AliT0Calibrator)

//____________________________________________________________________
  AliT0Calibrator::AliT0Calibrator():TNamed(),
				     fChannelWidth(0),  
				     fWalk(0)
				     
{
  //constructor

   AliT0Parameters* param = AliT0Parameters::Instance();
   param->Init();
  
  fChannelWidth = param->GetChannelWidth() ;  
  Double_t *grY ; //= new grY[2500] ;
  // Double_t *grX ;
  Int_t index[25000];
  Bool_t down=true;
  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);
 	fWalk.AddAtAndExpand(fu,i);
	
  if(fu) {
      Int_t np=fu->GetN();
      if(np>0) {
	grY = fu->GetY();
	//	grX[i] = fu->GetX();
	TMath::Sort(np, grY, index,down);
	fMaxValue[i]=Int_t(grY[index[0]]);

      }
    }
	
  }

  //  delete [] grY;
  //  delete [] grX;
  
  //
}
//_____________________________________________________________________________

AliT0Calibrator::AliT0Calibrator(const AliT0Calibrator &r): TNamed(),
  fChannelWidth(0),  
  fWalk(0)

{
  //
  // AliT0calibartor copy constructor
  //

  ((AliT0Calibrator &) r).Copy(*this);

}

//_____________________________________________________________________________
AliT0Calibrator &AliT0Calibrator::operator=(const AliT0Calibrator &r)
{
  //
  // Assignment operator
  //

  if (this != &r) ((AliT0Calibrator &) r).Copy(*this);
  return *this;

}



//____________________________________________________________________
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 timeEq=0, timeWalk=0;  
  //  TGraph* fu1 = param ->GetWalk(ipmt);
  //   TGraph* fu1  = param ->GetAmpLEDRec(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) ;
    //    timeEq= timeWalk - fTimeDelayCFD[ipmt];
    timeEq= timeWalk - fTimeDelayCFD[ipmt]; //for the same as cosmic
     AliDebug(10,Form(" ipmt %i time before %i timeWalk %i ,  qt %i timeEq %i \n ",
		 ipmt, time,timeWalk, qt, timeEq ));
  }
  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;
}