[u/mrichter/AliRoot.git] / TPC / TPCrec / AliTPCRecoParam.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 TPC reconstruction parameters                                  //
//                                                                           //  
//
/*
  The reconstruction parameters are used in the AliTPCclusterer and AliTPCtracker
  
  They are retrieved:
  0. User speciefied it in reconstruction macro
  1. if (not 0) from OCDB  - AliTPCcalibDB::GetRecoParam(eventtype)
  2. if (not 0 or 1) default parameter - High flux enevironment used  


  // Setting for systematic errors addition
  [0] - systematic RMSY
  [1] - systematic RMSZ
  [2] - systematic RMSSNP
  [3] - systematic RMSTheta
  [4] - systematic RMSCuravture -  systematic error in 1/cm not in 1/pt
  //
  //  How to add it example - 3 mm systematic error y, 3 cm systematic error z (drift)
  Double_t sysError[5]={0.3,3, 0.3/150., 3./150.,0.3/(150*150.)}
  param->SetSystematicError(sysError);

*/
                                                                           //
///////////////////////////////////////////////////////////////////////////////


#include "AliTPCRecoParam.h"

ClassImp(AliTPCRecoParam)


Bool_t AliTPCRecoParam::fgUseTimeCalibration=kTRUE; // flag usage the time dependent calibration
                                      // to be switched off for pass 0 reconstruction
                                      // Use static function, other option will be to use 
                                      // additional specific storage ?

//_____________________________________________________________________________
AliTPCRecoParam::AliTPCRecoParam():
  AliDetectorRecoParam(),
  fUseHLTClusters(4),  // use HLTorRAW data
  fUseHLTPreSeeding(0), // no pre-seeding for now
  fBClusterSharing(kTRUE),
  fCtgRange(1.05),       
  fMaxSnpTracker(0.95),
  fMaxSnpTrack(0.999),
  fUseOuterDetectors(kFALSE),
  fMaxChi2TPCTRD(36),     // maximal allowed chi2 between the TRD in and TPC out to be accepted for refit
  fMaxChi2TPCITS(36),     // maximal allowed chi2 between the ITS in and TPC out to be accepted for backpropagation
  fUseOulierClusterFilter(0),  // swith to use outlier cluster filter
  fDumpSignal(kFALSE),
  fFirstBin(0),
  fLastBin(-1),
  fBCalcPedestal(kFALSE),
  fBDoUnfold(kTRUE),
  fDumpAmplitudeMin(100),
  fMaxNoise(2.),
  //
  fUseOnePadCluster(kTRUE),
  fUseHLTOnePadCluster(kFALSE),
  fMinMaxCutAbs(4.),
  fMinLeftRightCutAbs(6.),
  fMinUpDownCutAbs(6.),
  //
  fMinMaxCutSigma(4.),
  fMinLeftRightCutSigma(7.),
  fMinUpDownCutSigma(8.),
  fMaxC(0.3),
  fBSpecialSeeding(kFALSE),
  fBKinkFinder(kTRUE),
  fLastSeedRowSec(120),
  fSeedGapPrim(6),
  fSeedGapSec(6),
  fUseFieldCorrection(0),      // use field correction
  fUseComposedCorrection(kFALSE),      // use field correction
  fUseRPHICorrection(0),      // use rphi correction
  fUseRadialCorrection(0),    // use radial correction
  fUseQuadrantAlignment(0),   // use quadrant alignment
  fUseSectorAlignment(0),     // use sector alignment
  fUseDriftCorrectionTime(1), // use drift correction time
  fUseDriftCorrectionGY(1),   // use drif correction global y
  fUseGainCorrectionTime(0),  // use gain correction time
  fUseExBCorrection(1),  // use ExB correction
  fUseMultiplicityCorrectionDedx(kTRUE), // use Dedx multiplicity correction
  fUseAlignmentTime(kTRUE),              // use time dependent alignment correction
  fUseIonTailCorrection(0),   // no ion tail correction for now
  fCrosstalkCorrection(0),   // crosstalk correction factor (from each signal substracted by (mean signal in wite patch)xfCrosstalkCorrection) - Effect important only after removing oc capacitors in 2012
  fCrosstalkCorrectionMissingCharge(1),   // crosstalk correction factor - missing charge (from each signal substracted by (mean signal in wite patch)xfCrosstalkCorrection) - Effect important only after removing oc capacitors in 2012
 //
  fUseTotCharge(kTRUE),          // switch use total or max charge
  fMinFraction(0.01),           // truncated mean - lower threshold
  fMaxFaction(0.7),            // truncated mean - upper threshold
  fNeighborRowsDedx(2),           // neighbour rows for below threshold dEdx calculation
  fGainCorrectionHVandPTMode(0), // switch for the usage of GainCorrectionHVandPT (see AliTPCcalibDB::GetGainCorrectionHVandPT   
  fSkipTimeBins(5),              // number of time bins to be skiiped (corrupted signal druing gating opening)
  fUseTOFCorrection(kTRUE),
  fUseSystematicCorrelation(kTRUE)
{
  //
  // constructor
  //
  SetName("TPC");
  SetTitle("TPC");
  for (Int_t i=0;i<5;i++) fSystematicErrors[i]=0;
  // systematic error parameterization at INNER wall of the TPC
  fSystematicErrorClusterInner[0]=0.5;   // 0.5 cm
  fSystematicErrorClusterInner[1]=5;     // 5 cm slope  
  //
  fSystematicErrorCluster[0]=0;   // sy cluster error
  fSystematicErrorCluster[1]=0;   // sz cluster error  

  //
  fCutSharedClusters[0]=0.5; // maximal allowed fraction of shared clusters - shorter track
  fCutSharedClusters[1]=0.25; // maximal allowed fraction of shared clusters - longer  track
  fClusterMaxRange[0]=1;     // y - pad      range
  fClusterMaxRange[1]=1;     // z - time bin range
  fKinkAngleCutChi2[0]=9;    // angular cut for kink finder - to create a kink
                             // ~ about 5 % rate  for high pt kink finder
  fKinkAngleCutChi2[1]=12;    // angular cut for kink finder - to use the partial track                             // form kink 
                             // ~ about 2 % rate  for high pt kink finder
}

//_____________________________________________________________________________
AliTPCRecoParam::~AliTPCRecoParam() 
{
  //
  // destructor
  //  
}

void AliTPCRecoParam::Print(const Option_t* /*option*/) const{
  //
  //
  //
  AliTPCRecoParam::Dump();
  printf("Systematic errors:\n");
  const char * cherrs[5]={"sy=","sz=","ssnp=","stheta=","s1pt="};
  for (Int_t i=0; i<5; i++){
    printf("%s%f\n",cherrs[i],fSystematicErrors[i]);
  }
}


AliTPCRecoParam *AliTPCRecoParam::GetLowFluxParam(){
  //
  // make default reconstruction  parameters for low  flux env.
  //
  AliTPCRecoParam *param = new AliTPCRecoParam;
  param->fCtgRange = 10;
  param->fFirstBin = 0;
  param->fLastBin  = 1000;
  param->SetName("Low Flux");
  param->SetTitle("Low Flux");
  return param;
}

AliTPCRecoParam *AliTPCRecoParam::GetHighFluxParam(){
  //
  // make reco parameters for high flux env.
  //
  AliTPCRecoParam *param = new AliTPCRecoParam;
  param->fCtgRange = 1.05;
  param->fFirstBin = 0;
  param->fLastBin  = 1000;  
  param->fUseTotCharge=kFALSE;
  param->SetName("High Flux");
  param->SetTitle("High Flux");
  return param;
}

AliTPCRecoParam *AliTPCRecoParam::GetHLTParam(){
  //
  // make reco parameters for high flux env.
  //
  AliTPCRecoParam *param = new AliTPCRecoParam;
  param->fCtgRange = 1.05;
  param->fFirstBin = 80;
  param->fLastBin  = 1000;  
  param->fMaxSnpTracker = 0.9; 
  param->fMaxC          = 0.06; 
  //
  param->SetName("Hlt Param");
  param->SetTitle("Hlt Param"); 
  param->fBKinkFinder   = kFALSE;
  return param;
}

AliTPCRecoParam *AliTPCRecoParam::GetLaserTestParam(Bool_t bPedestal){
  //
  // special setting for laser
  //
  AliTPCRecoParam *param = new AliTPCRecoParam;
  param->fDumpSignal=kTRUE;
  param->fCtgRange = 10.05;
  param->fFirstBin = 0;
  param->fLastBin  = 1000;
  param->fBCalcPedestal = bPedestal;
  param->fBDoUnfold     = kFALSE;
  param->fDumpAmplitudeMin = 150;
  param->fBKinkFinder   = kFALSE;
  param->fMaxSnpTracker = 0.98;
  param->fMaxC          = 0.02;
  param->fBSpecialSeeding = kTRUE;
  param->fUseTOFCorrection=kFALSE;
  param->fUseHLTClusters=1; // always RAW data
  //
  //
  param->SetName("Laser Flux");
  param->SetTitle("Laser Flux");
  return param;
}

AliTPCRecoParam *AliTPCRecoParam::GetCosmicTestParam(Bool_t bPedestal){
  //
  // special setting for cosmic 
  // 
  AliTPCRecoParam *param = new AliTPCRecoParam;
  param->fDumpSignal=kTRUE;
  param->fCtgRange = 10.05;    // full TPC
  param->fFirstBin = 60;
  param->fLastBin  = 1000;
  param->fBCalcPedestal = bPedestal;
  param->fBDoUnfold     = kFALSE;
  param->fBSpecialSeeding = kTRUE;
  param->fMaxC          = 0.07;
  param->fBKinkFinder   = kFALSE;
  param->fUseTOFCorrection =kFALSE;
  param->SetName("Cosmic Flux");
  param->SetTitle("Cosmic Flux");

  return param;
}


Bool_t  AliTPCRecoParam::GetUseTimeCalibration(){ 
  //
  // get
  //
  return fgUseTimeCalibration;
}
void    AliTPCRecoParam::SetUseTimeCalibration(Bool_t useTimeCalibration) {
  //
  // set 
  //
  fgUseTimeCalibration = useTimeCalibration;
}
Commit	Line	Data
669ce3c3	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
	17	///////////////////////////////////////////////////////////////////////////////
	18	// //
	19	// Class with TPC reconstruction parameters //
	20	// //
7a9e557b	21	//
7a9e557b	22	/*
829455ad	23	The reconstruction parameters are used in the AliTPCclusterer and AliTPCtracker
7a9e557b	24
	25	They are retrieved:
	26	0. User speciefied it in reconstruction macro
	27	1. if (not 0) from OCDB - AliTPCcalibDB::GetRecoParam(eventtype)
	28	2. if (not 0 or 1) default parameter - High flux enevironment used
	29
1cafd6f4	30
	31	// Setting for systematic errors addition
	32	[0] - systematic RMSY
	33	[1] - systematic RMSZ
	34	[2] - systematic RMSSNP
	35	[3] - systematic RMSTheta
	36	[4] - systematic RMSCuravture - systematic error in 1/cm not in 1/pt
	37	//
	38	// How to add it example - 3 mm systematic error y, 3 cm systematic error z (drift)
a69162c3	39	Double_t sysError[5]={0.3,3, 0.3/150., 3./150.,0.3/(150*150.)}
1cafd6f4	40	param->SetSystematicError(sysError);
1cafd6f4	41
7a9e557b	42	*/
7a9e557b	43	//
669ce3c3	44	///////////////////////////////////////////////////////////////////////////////
	45
	46
	47	#include "AliTPCRecoParam.h"
	48
	49	ClassImp(AliTPCRecoParam)
	50
	51
502d13b5	52	Bool_t AliTPCRecoParam::fgUseTimeCalibration=kTRUE; // flag usage the time dependent calibration
	53	// to be switched off for pass 0 reconstruction
	54	// Use static function, other option will be to use
	55	// additional specific storage ?
669ce3c3	56
	57	//_____________________________________________________________________________
	58	AliTPCRecoParam::AliTPCRecoParam():
b1c50d5c	59	AliDetectorRecoParam(),
1598ba75	60	fUseHLTClusters(4), // use HLTorRAW data
fa90aa2b	61	fUseHLTPreSeeding(0), // no pre-seeding for now
018e4f8a	62	fBClusterSharing(kTRUE),
669ce3c3	63	fCtgRange(1.05),
	64	fMaxSnpTracker(0.95),
	65	fMaxSnpTrack(0.999),
eba5bca1	66	fUseOuterDetectors(kFALSE),
cca0a82b	67	fMaxChi2TPCTRD(36), // maximal allowed chi2 between the TRD in and TPC out to be accepted for refit
8755834c	68	fMaxChi2TPCITS(36), // maximal allowed chi2 between the ITS in and TPC out to be accepted for backpropagation
44ca7282	69	fUseOulierClusterFilter(0), // swith to use outlier cluster filter
b127a65f	70	fDumpSignal(kFALSE),
669ce3c3	71	fFirstBin(0),
	72	fLastBin(-1),
	73	fBCalcPedestal(kFALSE),
	74	fBDoUnfold(kTRUE),
	75	fDumpAmplitudeMin(100),
837809ab	76	fMaxNoise(2.),
625fcba5	77	//
f47588e0	78	fUseOnePadCluster(kTRUE),
1598ba75	79	fUseHLTOnePadCluster(kFALSE),
625fcba5	80	fMinMaxCutAbs(4.),
	81	fMinLeftRightCutAbs(6.),
	82	fMinUpDownCutAbs(6.),
	83	//
940ed1f0	84	fMinMaxCutSigma(4.),
	85	fMinLeftRightCutSigma(7.),
	86	fMinUpDownCutSigma(8.),
b9cc16e4	87	fMaxC(0.3),
b9cc16e4	88	fBSpecialSeeding(kFALSE),
7d27c1df	89	fBKinkFinder(kTRUE),
9430b11a	90	fLastSeedRowSec(120),
5f638994	91	fSeedGapPrim(6),
5f638994	92	fSeedGapSec(6),
c2bb8829	93	fUseFieldCorrection(0), // use field correction
cf5b0aa0	94	fUseComposedCorrection(kFALSE), // use field correction
25617779	95	fUseRPHICorrection(0), // use rphi correction
	96	fUseRadialCorrection(0), // use radial correction
	97	fUseQuadrantAlignment(0), // use quadrant alignment
	98	fUseSectorAlignment(0), // use sector alignment
6be27415	99	fUseDriftCorrectionTime(1), // use drift correction time
6be27415	100	fUseDriftCorrectionGY(1), // use drif correction global y
25617779	101	fUseGainCorrectionTime(0), // use gain correction time
07c43acb	102	fUseExBCorrection(1), // use ExB correction
56b51ad9	103	fUseMultiplicityCorrectionDedx(kTRUE), // use Dedx multiplicity correction
56b51ad9	104	fUseAlignmentTime(kTRUE), // use time dependent alignment correction
fa90aa2b	105	fUseIonTailCorrection(0), // no ion tail correction for now
3449a97c	106	fCrosstalkCorrection(0), // crosstalk correction factor (from each signal substracted by (mean signal in wite patch)xfCrosstalkCorrection) - Effect important only after removing oc capacitors in 2012
	107	fCrosstalkCorrectionMissingCharge(1), // crosstalk correction factor - missing charge (from each signal substracted by (mean signal in wite patch)xfCrosstalkCorrection) - Effect important only after removing oc capacitors in 2012
	108	//
8f4188a7	109	fUseTotCharge(kTRUE), // switch use total or max charge
	110	fMinFraction(0.01), // truncated mean - lower threshold
	111	fMaxFaction(0.7), // truncated mean - upper threshold
6d64657a	112	fNeighborRowsDedx(2), // neighbour rows for below threshold dEdx calculation
00e32b66	113	fGainCorrectionHVandPTMode(0), // switch for the usage of GainCorrectionHVandPT (see AliTPCcalibDB::GetGainCorrectionHVandPT
efbb6a56	114	fSkipTimeBins(5), // number of time bins to be skiiped (corrupted signal druing gating opening)
e1dadcd0	115	fUseTOFCorrection(kTRUE),
e1dadcd0	116	fUseSystematicCorrelation(kTRUE)
669ce3c3	117	{
	118	//
	119	// constructor
	120	//
b1c50d5c	121	SetName("TPC");
b1c50d5c	122	SetTitle("TPC");
1cafd6f4	123	for (Int_t i=0;i<5;i++) fSystematicErrors[i]=0;
0b2c141b	124	// systematic error parameterization at INNER wall of the TPC
	125	fSystematicErrorClusterInner[0]=0.5; // 0.5 cm
	126	fSystematicErrorClusterInner[1]=5; // 5 cm slope
f183c81f	127	//
	128	fSystematicErrorCluster[0]=0; // sy cluster error
	129	fSystematicErrorCluster[1]=0; // sz cluster error
	130
0b2c141b	131	//
6fbe1e5c	132	fCutSharedClusters[0]=0.5; // maximal allowed fraction of shared clusters - shorter track
6fbe1e5c	133	fCutSharedClusters[1]=0.25; // maximal allowed fraction of shared clusters - longer track
eba5bca1	134	fClusterMaxRange[0]=1; // y - pad range
6fbe1e5c	135	fClusterMaxRange[1]=1; // z - time bin range
eba5bca1	136	fKinkAngleCutChi2[0]=9; // angular cut for kink finder - to create a kink
	137	// ~ about 5 % rate for high pt kink finder
	138	fKinkAngleCutChi2[1]=12; // angular cut for kink finder - to use the partial track // form kink
	139	// ~ about 2 % rate for high pt kink finder
669ce3c3	140	}
	141
	142	//_____________________________________________________________________________
	143	AliTPCRecoParam::~AliTPCRecoParam()
	144	{
	145	//
	146	// destructor
	147	//
	148	}
	149
90e46a63	150	void AliTPCRecoParam::Print(const Option_t* /option/) const{
	151	//
	152	//
	153	//
	154	AliTPCRecoParam::Dump();
	155	printf("Systematic errors:\n");
	156	const char * cherrs[5]={"sy=","sz=","ssnp=","stheta=","s1pt="};
	157	for (Int_t i=0; i<5; i++){
	158	printf("%s%f\n",cherrs[i],fSystematicErrors[i]);
	159	}
	160	}
669ce3c3	161
	162
	163	AliTPCRecoParam *AliTPCRecoParam::GetLowFluxParam(){
	164	//
	165	// make default reconstruction parameters for low flux env.
	166	//
	167	AliTPCRecoParam *param = new AliTPCRecoParam;
	168	param->fCtgRange = 10;
	169	param->fFirstBin = 0;
	170	param->fLastBin = 1000;
b447cbf9	171	param->SetName("Low Flux");
b447cbf9	172	param->SetTitle("Low Flux");
669ce3c3	173	return param;
	174	}
	175
	176	AliTPCRecoParam *AliTPCRecoParam::GetHighFluxParam(){
	177	//
	178	// make reco parameters for high flux env.
	179	//
	180	AliTPCRecoParam *param = new AliTPCRecoParam;
	181	param->fCtgRange = 1.05;
	182	param->fFirstBin = 0;
ec214961	183	param->fLastBin = 1000;
8f4188a7	184	param->fUseTotCharge=kFALSE;
b447cbf9	185	param->SetName("High Flux");
	186	param->SetTitle("High Flux");
	187	return param;
	188	}
	189
	190	AliTPCRecoParam *AliTPCRecoParam::GetHLTParam(){
	191	//
	192	// make reco parameters for high flux env.
	193	//
	194	AliTPCRecoParam *param = new AliTPCRecoParam;
	195	param->fCtgRange = 1.05;
	196	param->fFirstBin = 80;
	197	param->fLastBin = 1000;
	198	param->fMaxSnpTracker = 0.9;
1cafd6f4	199	param->fMaxC = 0.06;
b447cbf9	200	//
	201	param->SetName("Hlt Param");
	202	param->SetTitle("Hlt Param");
	203	param->fBKinkFinder = kFALSE;
669ce3c3	204	return param;
	205	}
	206
	207	AliTPCRecoParam *AliTPCRecoParam::GetLaserTestParam(Bool_t bPedestal){
	208	//
	209	// special setting for laser
	210	//
	211	AliTPCRecoParam *param = new AliTPCRecoParam;
b127a65f	212	param->fDumpSignal=kTRUE;
669ce3c3	213	param->fCtgRange = 10.05;
	214	param->fFirstBin = 0;
	215	param->fLastBin = 1000;
	216	param->fBCalcPedestal = bPedestal;
	217	param->fBDoUnfold = kFALSE;
fe5055e5	218	param->fDumpAmplitudeMin = 150;
669ce3c3	219	param->fBKinkFinder = kFALSE;
b9cc16e4	220	param->fMaxSnpTracker = 0.98;
	221	param->fMaxC = 0.02;
	222	param->fBSpecialSeeding = kTRUE;
9430b11a	223	param->fUseTOFCorrection=kFALSE;
9abc660b	224	param->fUseHLTClusters=1; // always RAW data
ec214961	225	//
ec214961	226	//
b447cbf9	227	param->SetName("Laser Flux");
b447cbf9	228	param->SetTitle("Laser Flux");
669ce3c3	229	return param;
	230	}
	231
	232	AliTPCRecoParam *AliTPCRecoParam::GetCosmicTestParam(Bool_t bPedestal){
	233	//
	234	// special setting for cosmic
	235	//
	236	AliTPCRecoParam *param = new AliTPCRecoParam;
b127a65f	237	param->fDumpSignal=kTRUE;
669ce3c3	238	param->fCtgRange = 10.05; // full TPC
	239	param->fFirstBin = 60;
	240	param->fLastBin = 1000;
	241	param->fBCalcPedestal = bPedestal;
	242	param->fBDoUnfold = kFALSE;
b9cc16e4	243	param->fBSpecialSeeding = kTRUE;
b9cc16e4	244	param->fMaxC = 0.07;
669ce3c3	245	param->fBKinkFinder = kFALSE;
9430b11a	246	param->fUseTOFCorrection =kFALSE;
b447cbf9	247	param->SetName("Cosmic Flux");
	248	param->SetTitle("Cosmic Flux");
	249
669ce3c3	250	return param;
	251	}
	252
	253
01e5d820	254	Bool_t AliTPCRecoParam::GetUseTimeCalibration(){
	255	//
	256	// get
	257	//
	258	return fgUseTimeCalibration;
	259	}
	260	void AliTPCRecoParam::SetUseTimeCalibration(Bool_t useTimeCalibration) {
	261	//
	262	// set
	263	//
	264	fgUseTimeCalibration = useTimeCalibration;
	265	}
669ce3c3	266