[u/mrichter/AliRoot.git] / TRD / AliTRDrecoParam.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$ */

///////////////////////////////////////////////////////////////////////////////
//                                                                           //
//  Parameter class for the TRD reconstruction                               //
//                                                                           //
//  Authors:                                                                 //
//    Alex Bercuci <A.Bercuci@gsi.de>                                        //
//    Markus Fasel <M.Fasel@gsi.de>                                          //
//                                                                           //
///////////////////////////////////////////////////////////////////////////////

#include "AliLog.h"

#include "AliTRDrecoParam.h"

ClassImp(AliTRDrecoParam)


//______________________________________________________________
AliTRDrecoParam::AliTRDrecoParam()
  :AliDetectorRecoParam()
  ,fkdNchdy(12.)
  ,fkMaxTheta(1.0)	
  ,fkMaxPhi(2.0) 
  ,fkRoad0y(6.0)
  ,fkRoad0z(8.5) 
  ,fkRoad1y(2.0)
  ,fkRoad1z(20.0)	
  ,fkRoad2y(3.0)
  ,fkRoad2z(20.0)
  ,fkPtThreshold(2.0) 
  ,fkPlaneQualityThreshold(5.0)// 4.2? under Investigation
  ,fkRoadzMultiplicator(1.5)
  ,fkFindable(.333)
  ,fkChi2Z(30./*14.*//*12.5*/)
  ,fkChi2Y(.25)
  ,fkChi2YSlope(7.73)
  ,fkChi2ZSlope(0.069)
  ,fChi2Cut(25)
  ,fkChi2YCut(0.5)
  ,fkPhiSlope(10.6)
  ,fkNMeanClusters(20.)
  ,fkNSigmaClusters(2.)
  ,fkNClusterNoise(0.)
  ,fkNMeanTracklets(5.5)
  ,fkTrackLikelihood(-15.)
  ,fNumberOfConfigs(3)
  ,fFlags(0)
  ,fRawStreamVersion("DEFAULT")
  ,fdzdxXcrossFactor(0.)
  ,fMinMaxCutSigma(4.)
  ,fMinLeftRightCutSigma(8.)
  ,fClusMaxThresh(4.5)
  ,fClusSigThresh(3.5)
  ,fTCnexp(1)
  ,fRecEveryNTB(1)
  ,fNumberOfPresamples(0)
  ,fNumberOfPostsamples(0)
{
  //
  // Default constructor
  //
  fSysCovMatrix[0] = 0.; // y direction (1 cm)
  fSysCovMatrix[1] = 0.; // z direction (1 cm)
  fSysCovMatrix[2] = 0.; // snp
  fSysCovMatrix[3] = 0.; // tgl
  fSysCovMatrix[4] = 0.; // 1/pt

  // Xe tail cancellation parameters
  fTCParams[0] = 1.156; // r1
  fTCParams[1] = 0.130; // r2
  fTCParams[2] = 0.114; // c1
  fTCParams[3] = 0.624; // c2
  // Ar tail cancellation parameters
  fTCParams[4] = 6.;    // r1
  fTCParams[5] = 0.62;  // r2
  fTCParams[6] = 0.0087;// c1
  fTCParams[7] = 0.07;  // c2

  memset(fPIDThreshold, 0, AliTRDCalPID::kNMom*sizeof(Double_t));
  memset(fStreamLevel, 0, kTRDreconstructionTasks * sizeof(Int_t));

  SetPIDmethod(AliTRDPIDResponse::kLQ1D);
  SetEightSlices();
  SetImproveTracklets();
  SetLUT();
  SetTailCancelation();
  SetTrackletParams();
}

//______________________________________________________________
AliTRDrecoParam::AliTRDrecoParam(const AliTRDrecoParam &ref)
  :AliDetectorRecoParam(ref)
  ,fkdNchdy(ref.fkdNchdy)
  ,fkMaxTheta(ref.fkMaxTheta)
  ,fkMaxPhi(ref.fkMaxPhi)
  ,fkRoad0y(ref.fkRoad0y)
  ,fkRoad0z(ref.fkRoad0z)
  ,fkRoad1y(ref.fkRoad1y)
  ,fkRoad1z(ref.fkRoad1z)
  ,fkRoad2y(ref.fkRoad2y)
  ,fkRoad2z(ref.fkRoad2z)
  ,fkPtThreshold(ref.fkPtThreshold)
  ,fkPlaneQualityThreshold(ref.fkPlaneQualityThreshold)
  ,fkRoadzMultiplicator(ref.fkRoadzMultiplicator)
  ,fkFindable(ref.fkFindable)
  ,fkChi2Z(ref.fkChi2Z)
  ,fkChi2Y(ref.fkChi2Y)
  ,fkChi2YSlope(ref.fkChi2YSlope)
  ,fkChi2ZSlope(ref.fkChi2ZSlope)
  ,fChi2Cut(ref.fChi2Cut)
  ,fkChi2YCut(ref.fkChi2YCut)
  ,fkPhiSlope(ref.fkPhiSlope)
  ,fkNMeanClusters(ref.fkNMeanClusters)
  ,fkNSigmaClusters(ref.fkNSigmaClusters)
  ,fkNClusterNoise(ref.fkNClusterNoise)
  ,fkNMeanTracklets(ref.fkNMeanTracklets)
  ,fkTrackLikelihood(ref.fkTrackLikelihood)
  ,fNumberOfConfigs(ref.fNumberOfConfigs)
  ,fFlags(ref.fFlags)
  ,fRawStreamVersion(ref.fRawStreamVersion)
  ,fdzdxXcrossFactor(ref.fdzdxXcrossFactor)
  ,fMinMaxCutSigma(ref.fMinMaxCutSigma)
  ,fMinLeftRightCutSigma(ref.fMinLeftRightCutSigma)
  ,fClusMaxThresh(ref.fClusMaxThresh)
  ,fClusSigThresh(ref.fClusSigThresh)
  ,fTCnexp(ref.fTCnexp)
  ,fRecEveryNTB(ref.fRecEveryNTB)
  ,fNumberOfPresamples(ref.fNumberOfPresamples)
  ,fNumberOfPostsamples(ref.fNumberOfPostsamples)
{
  //
  // Copy constructor
  //
  memcpy(fSysCovMatrix, ref.fSysCovMatrix, 5*sizeof(Double_t));
  memcpy(fTCParams, ref.fTCParams, 8*sizeof(Double_t));
  memcpy(fPIDThreshold, ref.fPIDThreshold, AliTRDCalPID::kNMom*sizeof(Double_t));
  memcpy(fStreamLevel, ref.fStreamLevel, kTRDreconstructionTasks * sizeof(Int_t));

  // tracklet params
  memcpy(fdzdxCorrFactor, ref.fdzdxCorrFactor, 2*sizeof(Double_t));
  memcpy(fdzdxCorrRCbias, ref.fdzdxCorrRCbias, 2*sizeof(Double_t));
  memcpy(fYcorrTailCancel, ref.fdzdxCorrRCbias, 6*sizeof(Double_t));
  memcpy(fS2Ycorr, ref.fS2Ycorr, 2*sizeof(Double_t));
}

//______________________________________________________________
AliTRDrecoParam& AliTRDrecoParam::operator=(const AliTRDrecoParam &ref)
{
  //
  // assignment operator
  //

  if(this == &ref) return *this;
  AliDetectorRecoParam::operator=(ref);
  fkdNchdy              = ref.fkdNchdy;
  fkMaxTheta            = ref.fkMaxTheta;
  fkMaxPhi              = ref.fkMaxPhi;
  fkRoad0y              = ref.fkRoad0y;
  fkRoad0z              = ref.fkRoad0z;
  fkRoad1y              = ref.fkRoad1y;
  fkRoad1z              = ref.fkRoad1z;
  fkRoad2y              = ref.fkRoad2y;
  fkRoad2z              = ref.fkRoad2z;
  fkPtThreshold         = ref.fkPtThreshold;
  fkPlaneQualityThreshold= ref.fkPlaneQualityThreshold;
  fkRoadzMultiplicator  = ref.fkRoadzMultiplicator;
  fkFindable            = ref.fkFindable;
  fkChi2Z               = ref.fkChi2Z;
  fkChi2Y               = ref.fkChi2Y;
  fkChi2YSlope          = ref.fkChi2YSlope;
  fkChi2ZSlope          = ref.fkChi2ZSlope;
  fChi2Cut            = ref.fChi2Cut;
  fkChi2YCut            = ref.fkChi2YCut;
  fkPhiSlope            = ref.fkPhiSlope;
  fkNMeanClusters       = ref.fkNMeanClusters;
  fkNSigmaClusters      = ref.fkNSigmaClusters;
  fkNClusterNoise       = ref.fkNClusterNoise;
  fkNMeanTracklets      = ref.fkNMeanTracklets;
  fkTrackLikelihood     = ref.fkTrackLikelihood;
  fNumberOfConfigs      = ref.fNumberOfConfigs;
  fFlags                = ref.fFlags;
  fRawStreamVersion     = ref.fRawStreamVersion;
  fdzdxXcrossFactor     = ref.fdzdxXcrossFactor;
  fMinMaxCutSigma       = ref.fMinMaxCutSigma;
  fMinLeftRightCutSigma = ref.fMinLeftRightCutSigma;
  fClusMaxThresh        = ref.fClusMaxThresh;
  fClusSigThresh        = ref.fClusSigThresh;
  fTCnexp               = ref.fTCnexp;
  fRecEveryNTB          = ref.fRecEveryNTB;
  fNumberOfPresamples   = ref.fNumberOfPresamples;
  fNumberOfPostsamples  = ref.fNumberOfPostsamples;

  memcpy(fSysCovMatrix, ref.fSysCovMatrix, 5*sizeof(Double_t));
  memcpy(fTCParams, ref.fTCParams, 8*sizeof(Double_t));
  memcpy(fPIDThreshold, ref.fPIDThreshold, AliTRDCalPID::kNMom*sizeof(Double_t));
  memcpy(fStreamLevel, ref.fStreamLevel, kTRDreconstructionTasks * sizeof(Int_t));

  // tracklet params
  memcpy(fdzdxCorrFactor, ref.fdzdxCorrFactor, 2*sizeof(Double_t));
  memcpy(fdzdxCorrRCbias, ref.fdzdxCorrRCbias, 2*sizeof(Double_t));
  memcpy(fYcorrTailCancel, ref.fdzdxCorrRCbias, 6*sizeof(Double_t));
  memcpy(fS2Ycorr, ref.fS2Ycorr, 2*sizeof(Double_t));
  return *this;
}

//______________________________________________________________
AliTRDrecoParam *AliTRDrecoParam::GetLowFluxParam()
{
  //
  // Parameters for the low flux environment
  //

  AliTRDrecoParam *rec = new AliTRDrecoParam();
  rec->fkdNchdy = 12.; // pp in TRD
  rec->SetVertexConstrained();
  rec->SetCheckTimeConsistency();
  return rec;

}

//______________________________________________________________
AliTRDrecoParam *AliTRDrecoParam::GetLowFluxHLTParam()
{
  //
  // Parameters for the high flux environment in HLT
  //

  AliTRDrecoParam *rec = GetLowFluxParam();
  rec->fNumberOfConfigs = 2;
  return rec;

}

//______________________________________________________________
AliTRDrecoParam *AliTRDrecoParam::GetHighFluxParam()
{
  //
  // Parameters for the high flux environment
  //

  AliTRDrecoParam *rec = new AliTRDrecoParam();
  rec->fkdNchdy = 4000.; // PbPb in TRD
  rec->SetVertexConstrained();
  rec->SetCheckTimeConsistency();
  return rec;

}

//______________________________________________________________
AliTRDrecoParam *AliTRDrecoParam::GetHighFluxHLTParam()
{
  //
  // Parameters for the high flux environment in HLT
  //

  AliTRDrecoParam *rec = GetHighFluxParam();
  rec->fNumberOfConfigs = 1;
  return rec;

}

//______________________________________________________________
AliTRDrecoParam *AliTRDrecoParam::GetCosmicTestParam()
{
  //
  // Parameters for the cosmics data
  //

  AliTRDrecoParam *par = new AliTRDrecoParam();
  par->fSysCovMatrix[0] = 2.; // y direction (1 cm)
  par->fSysCovMatrix[1] = 2.; // z direction (1 cm)
  par->fkChi2YSlope     = 0.11853;
  par->fkChi2ZSlope     = 0.04527;
  par->fkChi2YCut       = 25.;
  par->fkChi2YCut       = 1.;
  par->fkPhiSlope       = 10.; //3.17954;
  par->fkMaxTheta       = 2.1445;
  par->fkMaxPhi         = 2.7475;
  par->fkNMeanClusters  = 12.89;
  par->fkNSigmaClusters = 2.095;
  par->fkRoadzMultiplicator = 3.;
  par->fStreamLevel[kTracker] = 1;
  par->SetCheckTimeConsistency();
  return par;

}


//______________________________________________________________
Float_t AliTRDrecoParam::GetNClusters() const
{
  // Estimate the number of clusters in the TRD detector
  
  Float_t nclusters = (fkNMeanClusters + 2*fkNSigmaClusters)*fkNMeanTracklets*fkdNchdy;
  nclusters *= 1.+fkNClusterNoise;
  return nclusters;
}

//______________________________________________________________
void AliTRDrecoParam::SetPIDLQslices(Int_t s)
{
// Setting number of slices used by the PID LQ method s={1, 2}
// If PID NN is set this function will change to PID LQ.
 
  if(IsPIDNeuralNetwork()){
    AliWarning("PID set to NN. Changing to LQ.");
    SetPIDNeuralNetwork(kFALSE);
  } 

  switch(s){
  case 1: 
    if(TESTBIT(fFlags, kLQ2D)) CLRBIT(fFlags, kLQ2D);
    break;
  case 2:
    SETBIT(fFlags, kLQ2D);
    break;
  default:
    AliWarning(Form("N[%d] PID LQ slices not implemented. Using default 2.", s));
    SETBIT(fFlags, kLQ2D);
    break;
  }
}

//___________________________________________________
void  AliTRDrecoParam::SetTrackletParams(Double_t *par)
{
  // Load tracklet reconstruction parameters. If none are set use defaults
  if(par){
    // correct dzdx for the bias in z
    fdzdxCorrFactor[0] = par[0];  // !RC 
    fdzdxCorrFactor[1] = par[1];  // RC
    // correct dzdx in RC tracklets for the bias in cluster attachment
    fdzdxCorrRCbias[0] = par[2];   // dz/dx > 0  
    fdzdxCorrRCbias[1] = par[3];   // dz/dx < 0
    /// correct x_cross for the bias in dzdx
    fdzdxXcrossFactor  = par[4];
    // y linear q/pt correction due to wrong tail cancellation. 
    fYcorrTailCancel[0][0] = par[5]; fYcorrTailCancel[0][1] = par[6];  // opposite sign !RC
    fYcorrTailCancel[1][0] = par[7]; fYcorrTailCancel[1][1] = par[8];  // same sign !RC
    fYcorrTailCancel[2][0] = par[9]; fYcorrTailCancel[2][1] = par[10]; // RC
    // inflation factor of error parameterization in r-phi due to wrong estimation of residuals. 
    fS2Ycorr[0] = par[11];  // opposite sign, 
    fS2Ycorr[1] = par[12];  // same sign
    
  } else {
    // correct dzdx for the bias in z
    fdzdxCorrFactor[0] = 1.09;  // !RC 
    fdzdxCorrFactor[1] = 1.05;  // RC
    // correct dzdx in RC tracklets for the bias in cluster attachment
    fdzdxCorrRCbias[0] = 0.;     // dz/dx > 0  
    fdzdxCorrRCbias[1] = -0.012; // dz/dx < 0
    /// correct x_cross for the bias in dzdx
    fdzdxXcrossFactor  = 0.14;
    // y linear q/pt correction due to wrong tail cancellation. 
    fYcorrTailCancel[0][0] = 0.;   fYcorrTailCancel[0][1] = 0.027;  // opposite sign !RC
    fYcorrTailCancel[1][0] = 0.04; fYcorrTailCancel[1][1] = 0.027;  // same sign !RC
    fYcorrTailCancel[2][0] = 0.013;fYcorrTailCancel[2][1] = 0.018;  // RC
    // inflation factor of error parameterization in r-phi due to wrong estimation of residuals. 
    fS2Ycorr[0] = 5.;  // opposite sign, 
    fS2Ycorr[1] = 10;  // same sign
  }
}
Commit	Line	Data
e4f2f73d	1	/**************************************************************************
	2	* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	3	* *
	4	* Author: The ALICE Off-line Project. *
	5	* Contributors are mentioned in the code where appropriate. *
	6	* *
	7	* Permission to use, copy, modify and distribute this software and its *
	8	* documentation strictly for non-commercial purposes is hereby granted *
	9	* without fee, provided that the above copyright notice appears in all *
	10	* copies and that both the copyright notice and this permission notice *
	11	* appear in the supporting documentation. The authors make no claims *
	12	* about the suitability of this software for any purpose. It is *
	13	* provided "as is" without express or implied warranty. *
	14	**************************************************************************/
	15
	16	/* $Id$ */
	17
	18	///////////////////////////////////////////////////////////////////////////////
	19	// //
	20	// Parameter class for the TRD reconstruction //
	21	// //
	22	// Authors: //
	23	// Alex Bercuci <A.Bercuci@gsi.de> //
	24	// Markus Fasel <M.Fasel@gsi.de> //
	25	// //
	26	///////////////////////////////////////////////////////////////////////////////
	27
8df1f8f5	28	#include "AliLog.h"
8df1f8f5	29
e4f2f73d	30	#include "AliTRDrecoParam.h"
	31
	32	ClassImp(AliTRDrecoParam)
	33
a2fbb6ec	34
e4f2f73d	35	//______________________________________________________________
	36	AliTRDrecoParam::AliTRDrecoParam()
	37	:AliDetectorRecoParam()
8ae98148	38	,fkdNchdy(12.)
91834b8d	39	,fkMaxTheta(1.0)
91834b8d	40	,fkMaxPhi(2.0)
e4f2f73d	41	,fkRoad0y(6.0)
	42	,fkRoad0z(8.5)
	43	,fkRoad1y(2.0)
	44	,fkRoad1z(20.0)
	45	,fkRoad2y(3.0)
	46	,fkRoad2z(20.0)
b06a50a5	47	,fkPtThreshold(2.0)
e4f2f73d	48	,fkPlaneQualityThreshold(5.0)// 4.2? under Investigation
566bf887	49	,fkRoadzMultiplicator(1.5)
e4f2f73d	50	,fkFindable(.333)
eb38ed55	51	,fkChi2Z(30./14.//12.5/)
e4f2f73d	52	,fkChi2Y(.25)
5a2e200c	53	,fkChi2YSlope(7.73)
5a2e200c	54	,fkChi2ZSlope(0.069)
39d9c03c	55	,fChi2Cut(25)
8c499dbf	56	,fkChi2YCut(0.5)
5a2e200c	57	,fkPhiSlope(10.6)
8ae98148	58	,fkNMeanClusters(20.)
	59	,fkNSigmaClusters(2.)
	60	,fkNClusterNoise(0.)
	61	,fkNMeanTracklets(5.5)
e4f2f73d	62	,fkTrackLikelihood(-15.)
8c499dbf	63	,fNumberOfConfigs(3)
a2fbb6ec	64	,fFlags(0)
9dcc64cc	65	,fRawStreamVersion("DEFAULT")
86a261ed	66	,fdzdxXcrossFactor(0.)
df83a620	67	,fMinMaxCutSigma(4.)
df83a620	68	,fMinLeftRightCutSigma(8.)
fc546d21	69	,fClusMaxThresh(4.5)
fc546d21	70	,fClusSigThresh(3.5)
56b32fd7	71	,fTCnexp(1)
8c499dbf	72	,fRecEveryNTB(1)
4e459a9d	73	,fNumberOfPresamples(0)
4e459a9d	74	,fNumberOfPostsamples(0)
e4f2f73d	75	{
	76	//
	77	// Default constructor
	78	//
d937ad7a	79	fSysCovMatrix[0] = 0.; // y direction (1 cm)
d937ad7a	80	fSysCovMatrix[1] = 0.; // z direction (1 cm)
3afdab72	81	fSysCovMatrix[2] = 0.; // snp
	82	fSysCovMatrix[3] = 0.; // tgl
	83	fSysCovMatrix[4] = 0.; // 1/pt
22a4ab0c	84
a2fbb6ec	85	// Xe tail cancellation parameters
	86	fTCParams[0] = 1.156; // r1
	87	fTCParams[1] = 0.130; // r2
	88	fTCParams[2] = 0.114; // c1
	89	fTCParams[3] = 0.624; // c2
	90	// Ar tail cancellation parameters
	91	fTCParams[4] = 6.; // r1
	92	fTCParams[5] = 0.62; // r2
	93	fTCParams[6] = 0.0087;// c1
	94	fTCParams[7] = 0.07; // c2
	95
22a4ab0c	96	memset(fPIDThreshold, 0, AliTRDCalPID::kNMom*sizeof(Double_t));
a2fbb6ec	97	memset(fStreamLevel, 0, kTRDreconstructionTasks * sizeof(Int_t));
a2fbb6ec	98
9dcc64cc	99	SetPIDmethod(AliTRDPIDResponse::kLQ1D);
a2fbb6ec	100	SetEightSlices();
	101	SetImproveTracklets();
	102	SetLUT();
	103	SetTailCancelation();
86a261ed	104	SetTrackletParams();
e4f2f73d	105	}
e4f2f73d	106
a7ac01d2	107	//______________________________________________________________
	108	AliTRDrecoParam::AliTRDrecoParam(const AliTRDrecoParam &ref)
	109	:AliDetectorRecoParam(ref)
8ae98148	110	,fkdNchdy(ref.fkdNchdy)
a7ac01d2	111	,fkMaxTheta(ref.fkMaxTheta)
	112	,fkMaxPhi(ref.fkMaxPhi)
	113	,fkRoad0y(ref.fkRoad0y)
d6171f43	114	,fkRoad0z(ref.fkRoad0z)
a7ac01d2	115	,fkRoad1y(ref.fkRoad1y)
d6171f43	116	,fkRoad1z(ref.fkRoad1z)
a7ac01d2	117	,fkRoad2y(ref.fkRoad2y)
a7ac01d2	118	,fkRoad2z(ref.fkRoad2z)
d6171f43	119	,fkPtThreshold(ref.fkPtThreshold)
a7ac01d2	120	,fkPlaneQualityThreshold(ref.fkPlaneQualityThreshold)
566bf887	121	,fkRoadzMultiplicator(ref.fkRoadzMultiplicator)
a7ac01d2	122	,fkFindable(ref.fkFindable)
	123	,fkChi2Z(ref.fkChi2Z)
	124	,fkChi2Y(ref.fkChi2Y)
5a2e200c	125	,fkChi2YSlope(ref.fkChi2YSlope)
5a2e200c	126	,fkChi2ZSlope(ref.fkChi2ZSlope)
39d9c03c	127	,fChi2Cut(ref.fChi2Cut)
8c499dbf	128	,fkChi2YCut(ref.fkChi2YCut)
5a2e200c	129	,fkPhiSlope(ref.fkPhiSlope)
8ae98148	130	,fkNMeanClusters(ref.fkNMeanClusters)
	131	,fkNSigmaClusters(ref.fkNSigmaClusters)
	132	,fkNClusterNoise(ref.fkNClusterNoise)
	133	,fkNMeanTracklets(ref.fkNMeanTracklets)
a7ac01d2	134	,fkTrackLikelihood(ref.fkTrackLikelihood)
8c499dbf	135	,fNumberOfConfigs(ref.fNumberOfConfigs)
a2fbb6ec	136	,fFlags(ref.fFlags)
a2fbb6ec	137	,fRawStreamVersion(ref.fRawStreamVersion)
86a261ed	138	,fdzdxXcrossFactor(ref.fdzdxXcrossFactor)
a7ac01d2	139	,fMinMaxCutSigma(ref.fMinMaxCutSigma)
	140	,fMinLeftRightCutSigma(ref.fMinLeftRightCutSigma)
	141	,fClusMaxThresh(ref.fClusMaxThresh)
	142	,fClusSigThresh(ref.fClusSigThresh)
	143	,fTCnexp(ref.fTCnexp)
8c499dbf	144	,fRecEveryNTB(ref.fRecEveryNTB)
4e459a9d	145	,fNumberOfPresamples(ref.fNumberOfPresamples)
4e459a9d	146	,fNumberOfPostsamples(ref.fNumberOfPostsamples)
a7ac01d2	147	{
	148	//
	149	// Copy constructor
	150	//
a7ac01d2	151	memcpy(fSysCovMatrix, ref.fSysCovMatrix, 5*sizeof(Double_t));
a2fbb6ec	152	memcpy(fTCParams, ref.fTCParams, 8*sizeof(Double_t));
22a4ab0c	153	memcpy(fPIDThreshold, ref.fPIDThreshold, AliTRDCalPID::kNMom*sizeof(Double_t));
a2fbb6ec	154	memcpy(fStreamLevel, ref.fStreamLevel, kTRDreconstructionTasks * sizeof(Int_t));
86a261ed	155
	156	// tracklet params
	157	memcpy(fdzdxCorrFactor, ref.fdzdxCorrFactor, 2*sizeof(Double_t));
	158	memcpy(fdzdxCorrRCbias, ref.fdzdxCorrRCbias, 2*sizeof(Double_t));
	159	memcpy(fYcorrTailCancel, ref.fdzdxCorrRCbias, 6*sizeof(Double_t));
	160	memcpy(fS2Ycorr, ref.fS2Ycorr, 2*sizeof(Double_t));
a7ac01d2	161	}
a7ac01d2	162
d6171f43	163	//______________________________________________________________
f09f3167	164	AliTRDrecoParam& AliTRDrecoParam::operator=(const AliTRDrecoParam &ref)
d6171f43	165	{
	166	//
	167	// assignment operator
	168	//
	169
	170	if(this == &ref) return *this;
	171	AliDetectorRecoParam::operator=(ref);
	172	fkdNchdy = ref.fkdNchdy;
	173	fkMaxTheta = ref.fkMaxTheta;
	174	fkMaxPhi = ref.fkMaxPhi;
	175	fkRoad0y = ref.fkRoad0y;
	176	fkRoad0z = ref.fkRoad0z;
	177	fkRoad1y = ref.fkRoad1y;
	178	fkRoad1z = ref.fkRoad1z;
	179	fkRoad2y = ref.fkRoad2y;
	180	fkRoad2z = ref.fkRoad2z;
	181	fkPtThreshold = ref.fkPtThreshold;
	182	fkPlaneQualityThreshold= ref.fkPlaneQualityThreshold;
	183	fkRoadzMultiplicator = ref.fkRoadzMultiplicator;
	184	fkFindable = ref.fkFindable;
	185	fkChi2Z = ref.fkChi2Z;
	186	fkChi2Y = ref.fkChi2Y;
	187	fkChi2YSlope = ref.fkChi2YSlope;
	188	fkChi2ZSlope = ref.fkChi2ZSlope;
39d9c03c	189	fChi2Cut = ref.fChi2Cut;
d6171f43	190	fkChi2YCut = ref.fkChi2YCut;
	191	fkPhiSlope = ref.fkPhiSlope;
	192	fkNMeanClusters = ref.fkNMeanClusters;
	193	fkNSigmaClusters = ref.fkNSigmaClusters;
	194	fkNClusterNoise = ref.fkNClusterNoise;
	195	fkNMeanTracklets = ref.fkNMeanTracklets;
	196	fkTrackLikelihood = ref.fkTrackLikelihood;
	197	fNumberOfConfigs = ref.fNumberOfConfigs;
	198	fFlags = ref.fFlags;
	199	fRawStreamVersion = ref.fRawStreamVersion;
86a261ed	200	fdzdxXcrossFactor = ref.fdzdxXcrossFactor;
d6171f43	201	fMinMaxCutSigma = ref.fMinMaxCutSigma;
	202	fMinLeftRightCutSigma = ref.fMinLeftRightCutSigma;
	203	fClusMaxThresh = ref.fClusMaxThresh;
	204	fClusSigThresh = ref.fClusSigThresh;
	205	fTCnexp = ref.fTCnexp;
	206	fRecEveryNTB = ref.fRecEveryNTB;
	207	fNumberOfPresamples = ref.fNumberOfPresamples;
	208	fNumberOfPostsamples = ref.fNumberOfPostsamples;
	209
	210	memcpy(fSysCovMatrix, ref.fSysCovMatrix, 5*sizeof(Double_t));
	211	memcpy(fTCParams, ref.fTCParams, 8*sizeof(Double_t));
	212	memcpy(fPIDThreshold, ref.fPIDThreshold, AliTRDCalPID::kNMom*sizeof(Double_t));
	213	memcpy(fStreamLevel, ref.fStreamLevel, kTRDreconstructionTasks * sizeof(Int_t));
86a261ed	214
	215	// tracklet params
	216	memcpy(fdzdxCorrFactor, ref.fdzdxCorrFactor, 2*sizeof(Double_t));
	217	memcpy(fdzdxCorrRCbias, ref.fdzdxCorrRCbias, 2*sizeof(Double_t));
	218	memcpy(fYcorrTailCancel, ref.fdzdxCorrRCbias, 6*sizeof(Double_t));
	219	memcpy(fS2Ycorr, ref.fS2Ycorr, 2*sizeof(Double_t));
d6171f43	220	return *this;
	221	}
	222
e4f2f73d	223	//______________________________________________________________
	224	AliTRDrecoParam *AliTRDrecoParam::GetLowFluxParam()
	225	{
	226	//
	227	// Parameters for the low flux environment
	228	//
	229
8ae98148	230	AliTRDrecoParam *rec = new AliTRDrecoParam();
8ae98148	231	rec->fkdNchdy = 12.; // pp in TRD
a2fbb6ec	232	rec->SetVertexConstrained();
828c6f80	233	rec->SetCheckTimeConsistency();
8ae98148	234	return rec;
e4f2f73d	235
	236	}
	237
8c499dbf	238	//______________________________________________________________
	239	AliTRDrecoParam *AliTRDrecoParam::GetLowFluxHLTParam()
	240	{
	241	//
	242	// Parameters for the high flux environment in HLT
	243	//
	244
	245	AliTRDrecoParam *rec = GetLowFluxParam();
82ddb093	246	rec->fNumberOfConfigs = 2;
8c499dbf	247	return rec;
	248
	249	}
	250
e4f2f73d	251	//______________________________________________________________
	252	AliTRDrecoParam *AliTRDrecoParam::GetHighFluxParam()
	253	{
	254	//
	255	// Parameters for the high flux environment
	256	//
	257
fb872574	258	AliTRDrecoParam *rec = new AliTRDrecoParam();
8ae98148	259	rec->fkdNchdy = 4000.; // PbPb in TRD
a2fbb6ec	260	rec->SetVertexConstrained();
828c6f80	261	rec->SetCheckTimeConsistency();
fb872574	262	return rec;
e4f2f73d	263
e4f2f73d	264	}
9716329b	265
8c499dbf	266	//______________________________________________________________
	267	AliTRDrecoParam *AliTRDrecoParam::GetHighFluxHLTParam()
	268	{
	269	//
	270	// Parameters for the high flux environment in HLT
	271	//
	272
	273	AliTRDrecoParam *rec = GetHighFluxParam();
	274	rec->fNumberOfConfigs = 1;
8c499dbf	275	return rec;
	276
	277	}
	278
9716329b	279	//______________________________________________________________
	280	AliTRDrecoParam *AliTRDrecoParam::GetCosmicTestParam()
	281	{
	282	//
	283	// Parameters for the cosmics data
	284	//
	285
9716329b	286	AliTRDrecoParam *par = new AliTRDrecoParam();
6e39bde4	287	par->fSysCovMatrix[0] = 2.; // y direction (1 cm)
6e39bde4	288	par->fSysCovMatrix[1] = 2.; // z direction (1 cm)
b72f4eaf	289	par->fkChi2YSlope = 0.11853;
b72f4eaf	290	par->fkChi2ZSlope = 0.04527;
39d9c03c	291	par->fkChi2YCut = 25.;
b72f4eaf	292	par->fkChi2YCut = 1.;
6e39bde4	293	par->fkPhiSlope = 10.; //3.17954;
8ae98148	294	par->fkMaxTheta = 2.1445;
8ae98148	295	par->fkMaxPhi = 2.7475;
5a2e200c	296	par->fkNMeanClusters = 12.89;
5a2e200c	297	par->fkNSigmaClusters = 2.095;
566bf887	298	par->fkRoadzMultiplicator = 3.;
a2fbb6ec	299	par->fStreamLevel[kTracker] = 1;
828c6f80	300	par->SetCheckTimeConsistency();
9716329b	301	return par;
	302
	303	}
8ae98148	304
	305
	306	//______________________________________________________________
	307	Float_t AliTRDrecoParam::GetNClusters() const
	308	{
	309	// Estimate the number of clusters in the TRD detector
	310
	311	Float_t nclusters = (fkNMeanClusters + 2fkNSigmaClusters)fkNMeanTracklets*fkdNchdy;
	312	nclusters *= 1.+fkNClusterNoise;
	313	return nclusters;
	314	}
11d80e40	315
	316	//______________________________________________________________
	317	void AliTRDrecoParam::SetPIDLQslices(Int_t s)
	318	{
	319	// Setting number of slices used by the PID LQ method s={1, 2}
	320	// If PID NN is set this function will change to PID LQ.
	321
	322	if(IsPIDNeuralNetwork()){
	323	AliWarning("PID set to NN. Changing to LQ.");
	324	SetPIDNeuralNetwork(kFALSE);
	325	}
	326
	327	switch(s){
	328	case 1:
	329	if(TESTBIT(fFlags, kLQ2D)) CLRBIT(fFlags, kLQ2D);
	330	break;
	331	case 2:
	332	SETBIT(fFlags, kLQ2D);
	333	break;
	334	default:
	335	AliWarning(Form("N[%d] PID LQ slices not implemented. Using default 2.", s));
	336	SETBIT(fFlags, kLQ2D);
	337	break;
	338	}
	339	}
	340
86a261ed	341	//___________________________________________________
	342	void AliTRDrecoParam::SetTrackletParams(Double_t *par)
	343	{
	344	// Load tracklet reconstruction parameters. If none are set use defaults
	345	if(par){
	346	// correct dzdx for the bias in z
	347	fdzdxCorrFactor[0] = par[0]; // !RC
	348	fdzdxCorrFactor[1] = par[1]; // RC
	349	// correct dzdx in RC tracklets for the bias in cluster attachment
	350	fdzdxCorrRCbias[0] = par[2]; // dz/dx > 0
	351	fdzdxCorrRCbias[1] = par[3]; // dz/dx < 0
	352	/// correct x_cross for the bias in dzdx
	353	fdzdxXcrossFactor = par[4];
	354	// y linear q/pt correction due to wrong tail cancellation.
	355	fYcorrTailCancel[0][0] = par[5]; fYcorrTailCancel[0][1] = par[6]; // opposite sign !RC
	356	fYcorrTailCancel[1][0] = par[7]; fYcorrTailCancel[1][1] = par[8]; // same sign !RC
	357	fYcorrTailCancel[2][0] = par[9]; fYcorrTailCancel[2][1] = par[10]; // RC
	358	// inflation factor of error parameterization in r-phi due to wrong estimation of residuals.
	359	fS2Ycorr[0] = par[11]; // opposite sign,
	360	fS2Ycorr[1] = par[12]; // same sign
	361
	362	} else {
	363	// correct dzdx for the bias in z
	364	fdzdxCorrFactor[0] = 1.09; // !RC
	365	fdzdxCorrFactor[1] = 1.05; // RC
	366	// correct dzdx in RC tracklets for the bias in cluster attachment
	367	fdzdxCorrRCbias[0] = 0.; // dz/dx > 0
	368	fdzdxCorrRCbias[1] = -0.012; // dz/dx < 0
	369	/// correct x_cross for the bias in dzdx
	370	fdzdxXcrossFactor = 0.14;
	371	// y linear q/pt correction due to wrong tail cancellation.
	372	fYcorrTailCancel[0][0] = 0.; fYcorrTailCancel[0][1] = 0.027; // opposite sign !RC
	373	fYcorrTailCancel[1][0] = 0.04; fYcorrTailCancel[1][1] = 0.027; // same sign !RC
	374	fYcorrTailCancel[2][0] = 0.013;fYcorrTailCancel[2][1] = 0.018; // RC
	375	// inflation factor of error parameterization in r-phi due to wrong estimation of residuals.
	376	fS2Ycorr[0] = 5.; // opposite sign,
	377	fS2Ycorr[1] = 10; // same sign
	378	}
	379	}