Load some PWG4 and PWG3muon libraries needed to read the AOD from the Analysis Train

[u/mrichter/AliRoot.git] / TOF / AliTOFRecoParam.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.                  *
 **************************************************************************/


///////////////////////////////////////////////////////////////////////////////
//                                                                           //
// Class with TOF reconstruction parameters                                  //
//                                                                           //
//                                                                           //
///////////////////////////////////////////////////////////////////////////////


#include "AliLog.h"

#include "AliTOFRecoParam.h"

ClassImp(AliTOFRecoParam)

//_____________________________________________________________________________
AliTOFRecoParam::AliTOFRecoParam():
  AliDetectorRecoParam(),       
  fTimeZero(kFALSE),       
  fTimeZerofromT0(kFALSE),       
  fTimeZerofromTOF(kFALSE),       
  fTimeWalkCorr(kFALSE),       
  fApplyPbPbCuts(kFALSE),       
  fWindowSizeMaxY(50.),
  fWindowSizeMaxZ(35.),
  fWindowScaleFact(5.),
  fDistanceCut(3.),
  fSensRadius(379.5),
  fStepSize(0.1),
  fMaxChi2(10.),
  fMaxChi2TRD(150.),
  fTimeResolution(80.),
  fTimeNSigma(5.)
{
  //
  // constructor
  //
  SetNameTitle("TOF","TOF");
}
//_____________________________________________________________________________
AliTOFRecoParam::AliTOFRecoParam(const AliTOFRecoParam &p):
  AliDetectorRecoParam(p),       
  fTimeZero(p.fTimeZero),       
  fTimeZerofromT0(p.fTimeZerofromT0),       
  fTimeZerofromTOF(p.fTimeZerofromTOF),       
  fTimeWalkCorr(p.fTimeWalkCorr),       
  fApplyPbPbCuts(p.fApplyPbPbCuts),       
  fWindowSizeMaxY(p.fWindowSizeMaxY),
  fWindowSizeMaxZ(p.fWindowSizeMaxZ),
  fWindowScaleFact(p.fWindowScaleFact),
  fDistanceCut(p.fDistanceCut),
  fSensRadius(p.fSensRadius),
  fStepSize(p.fStepSize),
  fMaxChi2(p.fMaxChi2),
  fMaxChi2TRD(p.fMaxChi2TRD),
  fTimeResolution(p.fTimeResolution),
  fTimeNSigma(p.fTimeNSigma)
 { 
  //copy Ctor


}
//_____________________________________________________________________________
AliTOFRecoParam& AliTOFRecoParam::operator=(const AliTOFRecoParam &p)
{
  //
  // assign. operator
  //

  if (this == &p)
    return *this;
  
  AliDetectorRecoParam::operator=(p);
  fTimeZero=p.fTimeZero;       
   fTimeZerofromT0=p.fTimeZerofromT0;
   fTimeZerofromTOF=p.fTimeZerofromTOF;       
   fTimeWalkCorr=p.fTimeWalkCorr;       
   fApplyPbPbCuts=p.fApplyPbPbCuts;       
   fWindowSizeMaxY=p.fWindowSizeMaxY;
   fWindowSizeMaxZ=p.fWindowSizeMaxZ;
   fDistanceCut=p.fDistanceCut;
   fWindowScaleFact=p.fWindowScaleFact;
   fStepSize=p.fStepSize;
   fSensRadius=p.fSensRadius;
   fMaxChi2=p.fMaxChi2;
   fMaxChi2TRD=p.fMaxChi2TRD;
   fTimeResolution=p.fTimeResolution;
   fTimeNSigma=p.fTimeNSigma;   
   return *this;
}
//_____________________________________________________________________________
AliTOFRecoParam::~AliTOFRecoParam() 
{
  //
  // destructor
  //  
}

//_____________________________________________________________________________
AliTOFRecoParam *AliTOFRecoParam::GetPbPbparam(){
  //
  // set default reconstruction parameters for PbPb.
  //
  AliTOFRecoParam *param = new AliTOFRecoParam();
  param->fApplyPbPbCuts = kTRUE;
  param->fWindowScaleFact = 5.;
  param->fDistanceCut = 3.;
  return param;
}

//_____________________________________________________________________________
AliTOFRecoParam *AliTOFRecoParam::GetPPparam(){
  //
  // set default reconstruction parameters for PP.
  //
  AliTOFRecoParam *param = new AliTOFRecoParam();
  param->fApplyPbPbCuts = kFALSE;
  param->fWindowScaleFact = 5.;
  param->fDistanceCut = 10.;
  return param;
}

//_____________________________________________________________________________
AliTOFRecoParam *AliTOFRecoParam::GetCosmicMuonParam(){
  //
  // set default reconstruction parameters for cosmic muon run
  //
  AliTOFRecoParam *param = new AliTOFRecoParam();
  param->fApplyPbPbCuts = kFALSE;
  param->fWindowScaleFact = 5.;
  param->fDistanceCut = 10.;
  return param;
}

//_____________________________________________________________________________
void AliTOFRecoParam::PrintParameters() const
{
  //
  // Printing of the used TOF reconstruction parameters
  //

  AliInfo(Form(" Use Time Zero info in Pid: %i", fTimeZero));
  AliInfo(Form(" Use Time Zero as determined by T0: %i", fTimeZerofromT0));
  AliInfo(Form(" Use Time Zero as determined from TOF: %i",
	       fTimeZerofromTOF));
  AliInfo(Form(" Correct for signal time walk in z: %i", fTimeWalkCorr));
  AliInfo(Form(" Apply high flux cuts: %i", fApplyPbPbCuts));

  AliInfo(Form(" Cluster search window - size, Y: %f cm", fWindowSizeMaxY));
  AliInfo(Form(" Cluster search window - size, Z: %f cm", fWindowSizeMaxZ));
  AliInfo(Form(" Cluster search window - scale factor: %f",
	       fWindowScaleFact));

  AliInfo(Form(" Cut on the closest approach distance: %f", fDistanceCut));
  AliInfo(Form(" Average radius of sensitive volumes: %f cm", fSensRadius));
  AliInfo(Form(" Propagation step size: %f cm", fStepSize));
  AliInfo(Form(" Maximum X2 track-tof clusters: %f", fMaxChi2));

  AliInfo(Form(" Maximum X2 track-tof clusters (TRD): %f", fMaxChi2TRD));
  AliInfo(Form(" Time resolution for responce function in PID: %f ps",
	       fTimeResolution));
  AliInfo(Form("  N-Sigma Range used for responce function in PID: %f",
	       fTimeNSigma));

}