[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 "AliT0Reconstructor.h"
#include "AliT0RecoParam.h"
#include "AliCDBEntry.h"
#include "AliCDBManager.h"

#include <TGraph.h>
#include <TH1F.h>
#include <TMath.h>
#include <Riostream.h>

ClassImp(AliT0Calibrator)

//____________________________________________________________________
  AliT0Calibrator::AliT0Calibrator():TNamed(),
				     fChannelWidth(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 *grX ; 
  for (Int_t i=0; i<24; i++){
    fMaxValue[i]=0;
    fTimeDelayCFD[i] = Int_t (param->GetTimeDelayCFD(i));
     TGraph* fu = param ->GetWalk(i);
    // fWalk.AddAtAndExpand(fu,i);
    //TGraph* fu = param ->GetAmpLEDRec(i);
    fWalk.AddAtAndExpand(fu,i);
  }
  
}
//_____________________________________________________________________________

AliT0Calibrator::AliT0Calibrator(const AliT0Calibrator &r): TNamed(),
							    fChannelWidth(0),  
							    fWalk(0), 
							    fEqualized(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 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=(TGraph*) fWalk.At(ipmt);
  if(fu1 && fu1->GetN()>0) {
    walk = Int_t(fu1->Eval(Double_t(qt)));
  }
  
  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;
}
Commit	Line	Data
	1
	2	/**************************************************************************
	3	* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	4	* *
	5	* Author: The ALICE Off-line Project. *
	6	* Contributors are mentioned in the code where appropriate. *
	7	* *
	8	* Permission to use, copy, modify and distribute this software and its *
	9	* documentation strictly for non-commercial purposes is hereby granted *
	10	* without fee, provided that the above copyright notice appears in all *
	11	* copies and that both the copyright notice and this permission notice *
	12	* appear in the supporting documentation. The authors make no claims *
	13	* about the suitability of this software for any purpose. It is *
	14	* provided "as is" without express or implied warranty. *
	15	**************************************************************************/
	16
	17	/* $Id$ */
	18	/***********************************************************************
	19	* this class doing calibration during reconstruction
	20	* 2 steps:
	21	* - equalizing channels
	22	* - applying walk corrections
	23	*
	24	* Alla.Maevskaya@cern.ch
	25	*
	26	**********************************************************************/
	27
	28	//#include <Riostream.h>
	29
	30	#include "AliLog.h"
	31	#include "AliT0Parameters.h"
	32	#include "AliT0Calibrator.h"
	33	#include "AliT0Reconstructor.h"
	34	#include "AliT0RecoParam.h"
	35	#include "AliCDBEntry.h"
	36	#include "AliCDBManager.h"
	37
	38	#include <TGraph.h>
	39	#include <TH1F.h>
	40	#include <TMath.h>
	41	#include <Riostream.h>
	42
	43	ClassImp(AliT0Calibrator)
	44
	45	//____________________________________________________________________
	46	AliT0Calibrator::AliT0Calibrator():TNamed(),
	47	fChannelWidth(0),
	48	fWalk(0),
	49	fEqualized(0)
	50
	51	{
	52	//constructor
	53	printf(" AliT0Calibrator ::: AliT0RecoParam GetEq() %i\n", fEqualized);
	54	AliT0Parameters* param = AliT0Parameters::Instance();
	55	param->Init();
	56	//slewing correcion and equalizing channels
	57
	58	fChannelWidth = param->GetChannelWidth() ;
	59	// Double_t *grX ;
	60	for (Int_t i=0; i<24; i++){
	61	fMaxValue[i]=0;
	62	fTimeDelayCFD[i] = Int_t (param->GetTimeDelayCFD(i));
	63	TGraph* fu = param ->GetWalk(i);
	64	// fWalk.AddAtAndExpand(fu,i);
	65	//TGraph* fu = param ->GetAmpLEDRec(i);
	66	fWalk.AddAtAndExpand(fu,i);
	67	}
	68
	69	}
	70	//_____________________________________________________________________________
	71
	72	AliT0Calibrator::AliT0Calibrator(const AliT0Calibrator &r): TNamed(),
	73	fChannelWidth(0),
	74	fWalk(0),
	75	fEqualized(0)
	76
	77	{
	78	//
	79	// AliT0calibartor copy constructor
	80	//
	81
	82	((AliT0Calibrator &) r).Copy(*this);
	83
	84	}
	85
	86	//_____________________________________________________________________________
	87	AliT0Calibrator &AliT0Calibrator::operator=(const AliT0Calibrator &r)
	88	{
	89	//
	90	// Assignment operator
	91	//
	92
	93	if (this != &r) ((AliT0Calibrator &) r).Copy(*this);
	94	return *this;
	95
	96	}
	97
	98
	99	//____________________________________________________________________
	100
	101	Int_t AliT0Calibrator::WalkCorrection(Int_t refAmp, Int_t ipmt, Int_t qt, Int_t time)
	102
	103	{
	104	//
	105	// referemce amplitude for walk correction now read from RecoParam
	106
	107	Int_t walk=0;
	108	Int_t timeEq=0, timeWalk=0;
	109	TGraph fu1=(TGraph) fWalk.At(ipmt);
	110	if(fu1 && fu1->GetN()>0) {
	111	walk = Int_t(fu1->Eval(Double_t(qt)));
	112	}
	113
	114	if (fEqualized == 0)
	115	timeEq= time - fTimeDelayCFD[ipmt]-walk;
	116	else
	117	timeEq = time - walk - refAmp;
	118
	119	// printf(" ipmt %i time before %i timeWalk %i , walk %i qt %i fTimeDelayCFD[ipmt] %i timeEq %i \n ",
	120	// ipmt, time,timeWalk, walk, qt,fTimeDelayCFD[ipmt], timeEq );
	121	AliDebug(2,Form(" fEqualized %i ipmt %i refAmp %i time before %i timeWalk %i , walk %i qt %i timeEq %i, diff %i \n ",
	122	fEqualized, ipmt, refAmp, time,timeWalk, walk, qt, timeEq , fTimeDelayCFD[ipmt]));
	123
	124	return timeEq;
	125	}
	126
	127
	128