[u/mrichter/AliRoot.git] / ITS / AliITSRecoParam.cxx

/**************************************************************************
 * Copyright(c) 2007-2009, 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.                  *
 **************************************************************************/

#include "AliITSRecoParam.h"

/* $Id$ */

///////////////////////////////////////////////////////////////////////////////
//                                                                           //
// Class with ITS reconstruction parameters                                  //
// Origin: andrea.dainese@lnl.infn.it                                        //
//                                                                           //
///////////////////////////////////////////////////////////////////////////////


ClassImp(AliITSRecoParam)

const Int_t AliITSRecoParam::fgkLayersNotToSkip[AliITSgeomTGeo::kNLayers]={0,0,0,0,0,0};
const Int_t AliITSRecoParam::fgkLastLayerToTrackTo=0;
const Int_t AliITSRecoParam::fgkMaxDetectorPerLayer=1000;
const Double_t AliITSRecoParam::fgkriw=80.0;
const Double_t AliITSRecoParam::fgkdiw=0.0053;
const Double_t AliITSRecoParam::fgkX0iw=30.0;
const Double_t AliITSRecoParam::fgkrcd=61.0;
const Double_t AliITSRecoParam::fgkdcd=0.0053;
const Double_t AliITSRecoParam::fgkX0cd=30.0;
const Double_t AliITSRecoParam::fgkyr=12.8;
const Double_t AliITSRecoParam::fgkdr=0.03;
const Double_t AliITSRecoParam::fgkzm=0.2;
const Double_t AliITSRecoParam::fgkdm=0.40;
const Double_t AliITSRecoParam::fgkrs=50.0;
const Double_t AliITSRecoParam::fgkds=0.001;
const Double_t AliITSRecoParam::fgkrInsideITSscreen=49.0;
const Double_t AliITSRecoParam::fgkrInsideSPD1=3.7;
const Double_t AliITSRecoParam::fgkrPipe=3.;
const Double_t AliITSRecoParam::fgkrInsidePipe=2.7;
const Double_t AliITSRecoParam::fgkrOutsidePipe=3.3;
const Double_t AliITSRecoParam::fgkdPipe=0.0028;
const Double_t AliITSRecoParam::fgkrInsideShield[2]={7.5,25.0};
const Double_t AliITSRecoParam::fgkrOutsideShield[2]={10.5,30.0};
const Double_t AliITSRecoParam::fgkdshield[2]={0.0097,0.0034};
const Double_t AliITSRecoParam::fgkX0shield[2]={38.6,42.0};
const Double_t AliITSRecoParam::fgkX0Air=21.82;
const Double_t AliITSRecoParam::fgkX0Be=65.19;
const Double_t AliITSRecoParam::fgkBoundaryWidth=0.2;
const Double_t AliITSRecoParam::fgkDeltaXNeighbDets=0.5;
const Double_t AliITSRecoParam::fgkSPDdetzlength=6.960; // 7.072-2*0.056
const Double_t AliITSRecoParam::fgkSPDdetxlength=1.298; // 1.410-2*0.056

//_____________________________________________________________________________
AliITSRecoParam::AliITSRecoParam() : AliDetectorRecoParam(),
fMaxSnp(1.),
fNSigmaYLayerForRoadY(0),
fNSigmaRoadY(0),
fNSigmaZLayerForRoadZ(0),
fNSigmaRoadZ(0),
fNSigma2RoadZC(0),
fNSigma2RoadYC(0),
fNSigma2RoadZNonC(0),
fNSigma2RoadYNonC(0),
fMaxNormChi2NonCForHypothesis(0),
fMaxChi2(0),
fMaxRoad(0),
fMaxChi2In(0),
fChi2PerCluster(0),
fXV(0), 
fYV(0),
fZV(0),
fSigmaXV(0),
fSigmaYV(0),
fSigmaZV(0),
fVertexCut(0),
fMaxDZforPrimTrk(0),
fMaxDZToUseConstraint(0), 
fMaxDforV0dghtrForProlongation(0),
fMaxDForProlongation(0),
fMaxDZForProlongation(0),
fMinPtForProlongation(0),
fAddVirtualClustersInDeadZone(kFALSE),
fZWindowDeadZone(0),
fSigmaXDeadZoneHit2(0),
fSigmaZDeadZoneHit2(0),
fXPassDeadZoneHits(0),
fUseTGeoInTracker(3),
fAllowSharedClusters(kTRUE),
fClusterErrorsParam(1),
fComputePlaneEff(kFALSE),
fExtendedEtaAcceptance(kFALSE),
fUseDeadZonesFromOCDB(kFALSE),
fAllowProlongationWithEmptyRoad(kFALSE),
fFactorSAWindowSizes(1.),
fFindV0s(kTRUE),
fUseUnfoldingInClusterFinderSPD(kFALSE),
fUseUnfoldingInClusterFinderSDD(kTRUE),
fUseUnfoldingInClusterFinderSSD(kTRUE)
{
  //
  // constructor
  //
  SetName("ITS");
  SetTitle("ITS");

  SetLayersParameters();
  for(Int_t i=0; i<AliITSgeomTGeo::kNLayers; i++) fLayersToSkip[i]=0;
  SetUseTGeoInTracker(3);
  SetAllowSharedClusters(kTRUE);
  SetFindV0s(kTRUE);
  SetAddVirtualClustersInDeadZone(kFALSE);
  SetUseAmplitudeInfo(kTRUE);
  SetClusterErrorsParam(1);
}
//_____________________________________________________________________________
AliITSRecoParam::~AliITSRecoParam() 
{
  //
  // destructor
  //  
}
//_____________________________________________________________________________
AliITSRecoParam *AliITSRecoParam::GetHighFluxParam() 
{
  //
  // make default reconstruction  parameters for hig  flux env.
  //
  AliITSRecoParam *param = new AliITSRecoParam();

  param->fMaxSnp = 0.95;

  param->fNSigmaYLayerForRoadY = 4.;
  param->fNSigmaRoadY = 7.5;
  param->fNSigmaZLayerForRoadZ = 4.;
  param->fNSigmaRoadZ = 7.5;

  param->fNSigma2RoadZC = 60.; //7.75^2
  param->fNSigma2RoadYC = 60.; //7.75^2
  param->fNSigma2RoadZNonC = 50.; //7.07^2
  param->fNSigma2RoadYNonC = 50.; //7.07^2

  param->fMaxChi2PerCluster[0] = 11.; //7
  param->fMaxChi2PerCluster[1] = 12.; //5
  param->fMaxChi2PerCluster[2] = 12.; //8
  param->fMaxChi2PerCluster[3] = 5.;  //8
  param->fMaxChi2PerCluster[4] = 12.; //6.5

  param->fMaxNormChi2NonC[0] = 7.;
  param->fMaxNormChi2NonC[1] = 8.;
  param->fMaxNormChi2NonC[2] = 8.;
  param->fMaxNormChi2NonC[3] = 11.;
  param->fMaxNormChi2NonC[4] = 14.;
  param->fMaxNormChi2NonC[5] = 25.;

  param->fMaxNormChi2C[0] = 11.;
  param->fMaxNormChi2C[1] = 13.;
  param->fMaxNormChi2C[2] = 15.;
  param->fMaxNormChi2C[3] = 18.;
  param->fMaxNormChi2C[4] = 30.;
  param->fMaxNormChi2C[5] = 35.;

  param->fMaxNormChi2NonCForHypothesis = 7.;
  
  param->fMaxChi2 = 35.;

  param->fMaxChi2s[0] = 25.; //40   
  param->fMaxChi2s[1] = 25.; //40   
  param->fMaxChi2s[2] = 25.; //40   
  param->fMaxChi2s[3] = 25.; //40   
  param->fMaxChi2s[4] = 40.; //40   
  param->fMaxChi2s[5] = 50.; //40

  param->fMaxRoad = 6.;

  // not used
  param->fMaxChi2In = 16.;
   
  param->fMaxChi2sR[0] = 10.;   
  param->fMaxChi2sR[1] = 10.;   
  param->fMaxChi2sR[2] = 10.;   
  param->fMaxChi2sR[3] = 10.;   
  param->fMaxChi2sR[4] = 30.;   
  param->fMaxChi2sR[5] = 40.;   

  param->fChi2PerCluster = 9.;
  // not used

  param->fXV = 0.;
  param->fYV = 0.;
  param->fZV = 0.;
  param->fSigmaXV = 0.0050;
  param->fSigmaYV = 0.0050;
  param->fSigmaZV = 0.0100;

  param->fVertexCut = 25.;

  param->fMaxDZforPrimTrk = 0.4;
  param->fMaxDZToUseConstraint = 3.;

  param->fMaxDforV0dghtrForProlongation = 30.;
  param->fMaxDForProlongation = 40.;//10.;
  param->fMaxDZForProlongation = 60.;//20.;
  param->fMinPtForProlongation = 0.120;

  param->fZWindowDeadZone = 2.0;
  param->fSigmaXDeadZoneHit2 = 0.004/12.;
  param->fSigmaZDeadZoneHit2 = 0.001/12.;
  param->fXPassDeadZoneHits = 0.018;

  
  return param;
}
//_____________________________________________________________________________
AliITSRecoParam *AliITSRecoParam::GetLowFluxParam() 
{
  //
  // make default reconstruction  parameters for low  flux env.
  //
  return GetHighFluxParam();
}
//_____________________________________________________________________________
AliITSRecoParam *AliITSRecoParam::GetCosmicTestParam() 
{
  //
  // make default reconstruction  parameters for cosmics
  //
  AliITSRecoParam *param = new AliITSRecoParam();

  // larger seach windows for SA (in case of large misalignments)
  param->SetFactorSAWindowSizes(3.); 

  param->fMaxSnp = 0.95;

  param->fNSigmaYLayerForRoadY = 4.;
  param->fNSigmaRoadY = 7.5;
  param->fNSigmaZLayerForRoadZ = 4.;
  param->fNSigmaRoadZ = 7.5;

  param->fNSigma2RoadZC = 60.; //7.75^2
  param->fNSigma2RoadYC = 60.; //7.75^2
  param->fNSigma2RoadZNonC = 50.; //7.07^2
  param->fNSigma2RoadYNonC = 50.; //7.07^2

  param->fMaxChi2PerCluster[0] = 11.; //7
  param->fMaxChi2PerCluster[1] = 12.; //5
  param->fMaxChi2PerCluster[2] = 12.; //8
  param->fMaxChi2PerCluster[3] = 5.;  //8
  param->fMaxChi2PerCluster[4] = 12.; //6.5

  param->fMaxNormChi2NonC[0] = 7.;
  param->fMaxNormChi2NonC[1] = 8.;
  param->fMaxNormChi2NonC[2] = 8.;
  param->fMaxNormChi2NonC[3] = 11.;
  param->fMaxNormChi2NonC[4] = 14.;
  param->fMaxNormChi2NonC[5] = 25.;

  param->fMaxNormChi2C[0] = 11.;
  param->fMaxNormChi2C[1] = 13.;
  param->fMaxNormChi2C[2] = 15.;
  param->fMaxNormChi2C[3] = 18.;
  param->fMaxNormChi2C[4] = 30.;
  param->fMaxNormChi2C[5] = 35.;

  param->fMaxNormChi2NonCForHypothesis = 7.;
  
  param->fMaxChi2 = 35.;

  param->fMaxChi2s[0] = 25.; //40   
  param->fMaxChi2s[1] = 25.; //40   
  param->fMaxChi2s[2] = 25.; //40   
  param->fMaxChi2s[3] = 25.; //40   
  param->fMaxChi2s[4] = 40.; //40   
  param->fMaxChi2s[5] = 50.; //40

  param->fMaxRoad = 6.;

  // not used
  param->fMaxChi2In = 16.;
   
  param->fMaxChi2sR[0] = 10.;   
  param->fMaxChi2sR[1] = 10.;   
  param->fMaxChi2sR[2] = 10.;   
  param->fMaxChi2sR[3] = 10.;   
  param->fMaxChi2sR[4] = 30.;   
  param->fMaxChi2sR[5] = 40.;   

  param->fChi2PerCluster = 9.;
  // not used

  param->fXV = 0.;
  param->fYV = 0.;
  param->fZV = 0.;
  param->fSigmaXV = 0.0050;
  param->fSigmaYV = 0.0050;
  param->fSigmaZV = 0.0100;

  param->fVertexCut = 25.;

  param->fMaxDZforPrimTrk = 0.4;
  param->fMaxDZToUseConstraint = 3.;

  param->fMaxDforV0dghtrForProlongation = 30.;
  param->fMaxDForProlongation = 10.;
  param->fMaxDZForProlongation = 20.;
  param->fMinPtForProlongation = 0.120;

  param->fZWindowDeadZone = 2.0;
  param->fSigmaXDeadZoneHit2 = 0.004/12.;
  param->fSigmaZDeadZoneHit2 = 0.001/12.;
  param->fXPassDeadZoneHits = 0.018;

  
  return param;
}
//_____________________________________________________________________________
AliITSRecoParam *AliITSRecoParam::GetPlaneEffParam(Int_t i)
{
  //
  // make special reconstruction parameters for Plane Efficiency study on layer i
  //
  AliITSRecoParam *param;
  param = GetHighFluxParam();
  param->SetComputePlaneEff();
  param->SetLayerToSkip(i);
  return param;
}
//_____________________________________________________________________________
void AliITSRecoParam::SetLayersParameters() 
{
  //
  // number of layers and layers spatial resolutions
  //

  // spatial resolutions of the detectors
  // y: 12 12 38 38 20 20 micron
  fSigmaY2[0]=1.44e-6;
  fSigmaY2[1]=1.44e-6;
  fSigmaY2[2]=1.444e-5;
  fSigmaY2[3]=1.444e-5;
  fSigmaY2[4]=4.0e-6;
  fSigmaY2[5]=4.0e-6;
  // z: 120 120 28 28 830 830 micron
  fSigmaZ2[0]=1.44e-4;
  fSigmaZ2[1]=1.44e-4;
  fSigmaZ2[2]=7.84e-6;
  fSigmaZ2[3]=7.84e-6;
  fSigmaZ2[4]=6.889e-3;
  fSigmaZ2[5]=6.889e-3;

  return;
}
Commit	Line	Data
44347160	1	/**************************************************************************
572f41f9	2	* Copyright(c) 2007-2009, ALICE Experiment at CERN, All rights reserved. *
44347160	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
572f41f9	16	#include "AliITSRecoParam.h"
	17
	18	/* $Id$ */
44347160	19
	20	///////////////////////////////////////////////////////////////////////////////
	21	// //
	22	// Class with ITS reconstruction parameters //
	23	// Origin: andrea.dainese@lnl.infn.it //
	24	// //
	25	///////////////////////////////////////////////////////////////////////////////
	26
	27
44347160	28
	29	ClassImp(AliITSRecoParam)
	30
e50912db	31	const Int_t AliITSRecoParam::fgkLayersNotToSkip[AliITSgeomTGeo::kNLayers]={0,0,0,0,0,0};
	32	const Int_t AliITSRecoParam::fgkLastLayerToTrackTo=0;
	33	const Int_t AliITSRecoParam::fgkMaxDetectorPerLayer=1000;
	34	const Double_t AliITSRecoParam::fgkriw=80.0;
	35	const Double_t AliITSRecoParam::fgkdiw=0.0053;
	36	const Double_t AliITSRecoParam::fgkX0iw=30.0;
	37	const Double_t AliITSRecoParam::fgkrcd=61.0;
	38	const Double_t AliITSRecoParam::fgkdcd=0.0053;
	39	const Double_t AliITSRecoParam::fgkX0cd=30.0;
	40	const Double_t AliITSRecoParam::fgkyr=12.8;
	41	const Double_t AliITSRecoParam::fgkdr=0.03;
	42	const Double_t AliITSRecoParam::fgkzm=0.2;
	43	const Double_t AliITSRecoParam::fgkdm=0.40;
	44	const Double_t AliITSRecoParam::fgkrs=50.0;
	45	const Double_t AliITSRecoParam::fgkds=0.001;
	46	const Double_t AliITSRecoParam::fgkrInsideITSscreen=49.0;
	47	const Double_t AliITSRecoParam::fgkrInsideSPD1=3.7;
	48	const Double_t AliITSRecoParam::fgkrPipe=3.;
	49	const Double_t AliITSRecoParam::fgkrInsidePipe=2.7;
	50	const Double_t AliITSRecoParam::fgkrOutsidePipe=3.3;
	51	const Double_t AliITSRecoParam::fgkdPipe=0.0028;
	52	const Double_t AliITSRecoParam::fgkrInsideShield[2]={7.5,25.0};
	53	const Double_t AliITSRecoParam::fgkrOutsideShield[2]={10.5,30.0};
	54	const Double_t AliITSRecoParam::fgkdshield[2]={0.0097,0.0034};
	55	const Double_t AliITSRecoParam::fgkX0shield[2]={38.6,42.0};
	56	const Double_t AliITSRecoParam::fgkX0Air=21.82;
	57	const Double_t AliITSRecoParam::fgkX0Be=65.19;
	58	const Double_t AliITSRecoParam::fgkBoundaryWidth=0.2;
	59	const Double_t AliITSRecoParam::fgkDeltaXNeighbDets=0.5;
	60	const Double_t AliITSRecoParam::fgkSPDdetzlength=6.960; // 7.072-2*0.056
	61	const Double_t AliITSRecoParam::fgkSPDdetxlength=1.298; // 1.410-2*0.056
44347160	62
44347160	63	//_____________________________________________________________________________
6518a6c5	64	AliITSRecoParam::AliITSRecoParam() : AliDetectorRecoParam(),
e50912db	65	fMaxSnp(1.),
572f41f9	66	fNSigmaYLayerForRoadY(0),
	67	fNSigmaRoadY(0),
	68	fNSigmaZLayerForRoadZ(0),
	69	fNSigmaRoadZ(0),
	70	fNSigma2RoadZC(0),
	71	fNSigma2RoadYC(0),
	72	fNSigma2RoadZNonC(0),
	73	fNSigma2RoadYNonC(0),
	74	fMaxNormChi2NonCForHypothesis(0),
	75	fMaxChi2(0),
	76	fMaxRoad(0),
	77	fMaxChi2In(0),
	78	fChi2PerCluster(0),
	79	fXV(0),
	80	fYV(0),
	81	fZV(0),
	82	fSigmaXV(0),
	83	fSigmaYV(0),
	84	fSigmaZV(0),
	85	fVertexCut(0),
	86	fMaxDZforPrimTrk(0),
	87	fMaxDZToUseConstraint(0),
	88	fMaxDforV0dghtrForProlongation(0),
	89	fMaxDForProlongation(0),
	90	fMaxDZForProlongation(0),
	91	fMinPtForProlongation(0),
6518a6c5	92	fAddVirtualClustersInDeadZone(kFALSE),
572f41f9	93	fZWindowDeadZone(0),
	94	fSigmaXDeadZoneHit2(0),
	95	fSigmaZDeadZoneHit2(0),
	96	fXPassDeadZoneHits(0),
6518a6c5	97	fUseTGeoInTracker(3),
572f41f9	98	fAllowSharedClusters(kTRUE),
e50912db	99	fClusterErrorsParam(1),
4a66240a	100	fComputePlaneEff(kFALSE),
2755f080	101	fExtendedEtaAcceptance(kFALSE),
6518a6c5	102	fUseDeadZonesFromOCDB(kFALSE),
ae00569a	103	fAllowProlongationWithEmptyRoad(kFALSE),
2755f080	104	fFactorSAWindowSizes(1.),
6518a6c5	105	fFindV0s(kTRUE),
	106	fUseUnfoldingInClusterFinderSPD(kFALSE),
	107	fUseUnfoldingInClusterFinderSDD(kTRUE),
	108	fUseUnfoldingInClusterFinderSSD(kTRUE)
44347160	109	{
	110	//
	111	// constructor
	112	//
6518a6c5	113	SetName("ITS");
	114	SetTitle("ITS");
	115
44347160	116	SetLayersParameters();
2755f080	117	for(Int_t i=0; i<AliITSgeomTGeo::kNLayers; i++) fLayersToSkip[i]=0;
6518a6c5	118	SetUseTGeoInTracker(3);
afd25725	119	SetAllowSharedClusters(kTRUE);
afd25725	120	SetFindV0s(kTRUE);
6518a6c5	121	SetAddVirtualClustersInDeadZone(kFALSE);
572f41f9	122	SetUseAmplitudeInfo(kTRUE);
e50912db	123	SetClusterErrorsParam(1);
44347160	124	}
	125	//_____________________________________________________________________________
	126	AliITSRecoParam::~AliITSRecoParam()
	127	{
	128	//
	129	// destructor
	130	//
	131	}
	132	//_____________________________________________________________________________
	133	AliITSRecoParam *AliITSRecoParam::GetHighFluxParam()
	134	{
	135	//
	136	// make default reconstruction parameters for hig flux env.
	137	//
	138	AliITSRecoParam *param = new AliITSRecoParam();
	139
	140	param->fMaxSnp = 0.95;
	141
	142	param->fNSigmaYLayerForRoadY = 4.;
	143	param->fNSigmaRoadY = 7.5;
	144	param->fNSigmaZLayerForRoadZ = 4.;
	145	param->fNSigmaRoadZ = 7.5;
	146
	147	param->fNSigma2RoadZC = 60.; //7.75^2
	148	param->fNSigma2RoadYC = 60.; //7.75^2
	149	param->fNSigma2RoadZNonC = 50.; //7.07^2
	150	param->fNSigma2RoadYNonC = 50.; //7.07^2
	151
	152	param->fMaxChi2PerCluster[0] = 11.; //7
	153	param->fMaxChi2PerCluster[1] = 12.; //5
	154	param->fMaxChi2PerCluster[2] = 12.; //8
	155	param->fMaxChi2PerCluster[3] = 5.; //8
	156	param->fMaxChi2PerCluster[4] = 12.; //6.5
	157
	158	param->fMaxNormChi2NonC[0] = 7.;
	159	param->fMaxNormChi2NonC[1] = 8.;
	160	param->fMaxNormChi2NonC[2] = 8.;
	161	param->fMaxNormChi2NonC[3] = 11.;
	162	param->fMaxNormChi2NonC[4] = 14.;
	163	param->fMaxNormChi2NonC[5] = 25.;
	164
	165	param->fMaxNormChi2C[0] = 11.;
	166	param->fMaxNormChi2C[1] = 13.;
	167	param->fMaxNormChi2C[2] = 15.;
	168	param->fMaxNormChi2C[3] = 18.;
	169	param->fMaxNormChi2C[4] = 30.;
	170	param->fMaxNormChi2C[5] = 35.;
afd25725	171
afd25725	172	param->fMaxNormChi2NonCForHypothesis = 7.;
44347160	173
	174	param->fMaxChi2 = 35.;
	175
	176	param->fMaxChi2s[0] = 25.; //40
	177	param->fMaxChi2s[1] = 25.; //40
	178	param->fMaxChi2s[2] = 25.; //40
	179	param->fMaxChi2s[3] = 25.; //40
	180	param->fMaxChi2s[4] = 40.; //40
	181	param->fMaxChi2s[5] = 50.; //40
	182
	183	param->fMaxRoad = 6.;
	184
	185	// not used
	186	param->fMaxChi2In = 16.;
	187
	188	param->fMaxChi2sR[0] = 10.;
	189	param->fMaxChi2sR[1] = 10.;
	190	param->fMaxChi2sR[2] = 10.;
	191	param->fMaxChi2sR[3] = 10.;
	192	param->fMaxChi2sR[4] = 30.;
	193	param->fMaxChi2sR[5] = 40.;
	194
	195	param->fChi2PerCluster = 9.;
	196	// not used
	197
	198	param->fXV = 0.;
	199	param->fYV = 0.;
	200	param->fZV = 0.;
	201	param->fSigmaXV = 0.0050;
	202	param->fSigmaYV = 0.0050;
	203	param->fSigmaZV = 0.0100;
	204
	205	param->fVertexCut = 25.;
afd25725	206
	207	param->fMaxDZforPrimTrk = 0.4;
	208	param->fMaxDZToUseConstraint = 3.;
	209
	210	param->fMaxDforV0dghtrForProlongation = 30.;
6518a6c5	211	param->fMaxDForProlongation = 40.;//10.;
6518a6c5	212	param->fMaxDZForProlongation = 60.;//20.;
afd25725	213	param->fMinPtForProlongation = 0.120;
	214
	215	param->fZWindowDeadZone = 2.0;
	216	param->fSigmaXDeadZoneHit2 = 0.004/12.;
	217	param->fSigmaZDeadZoneHit2 = 0.001/12.;
	218	param->fXPassDeadZoneHits = 0.018;
	219
44347160	220
	221	return param;
	222	}
	223	//_____________________________________________________________________________
	224	AliITSRecoParam *AliITSRecoParam::GetLowFluxParam()
	225	{
	226	//
	227	// make default reconstruction parameters for low flux env.
	228	//
	229	return GetHighFluxParam();
	230	}
	231	//_____________________________________________________________________________
	232	AliITSRecoParam *AliITSRecoParam::GetCosmicTestParam()
	233	{
	234	//
	235	// make default reconstruction parameters for cosmics
	236	//
e340bb86	237	AliITSRecoParam *param = new AliITSRecoParam();
e340bb86	238
2755f080	239	// larger seach windows for SA (in case of large misalignments)
2755f080	240	param->SetFactorSAWindowSizes(3.);
e340bb86	241
	242	param->fMaxSnp = 0.95;
	243
	244	param->fNSigmaYLayerForRoadY = 4.;
	245	param->fNSigmaRoadY = 7.5;
	246	param->fNSigmaZLayerForRoadZ = 4.;
	247	param->fNSigmaRoadZ = 7.5;
	248
	249	param->fNSigma2RoadZC = 60.; //7.75^2
	250	param->fNSigma2RoadYC = 60.; //7.75^2
	251	param->fNSigma2RoadZNonC = 50.; //7.07^2
	252	param->fNSigma2RoadYNonC = 50.; //7.07^2
	253
	254	param->fMaxChi2PerCluster[0] = 11.; //7
	255	param->fMaxChi2PerCluster[1] = 12.; //5
	256	param->fMaxChi2PerCluster[2] = 12.; //8
	257	param->fMaxChi2PerCluster[3] = 5.; //8
	258	param->fMaxChi2PerCluster[4] = 12.; //6.5
	259
	260	param->fMaxNormChi2NonC[0] = 7.;
	261	param->fMaxNormChi2NonC[1] = 8.;
	262	param->fMaxNormChi2NonC[2] = 8.;
	263	param->fMaxNormChi2NonC[3] = 11.;
	264	param->fMaxNormChi2NonC[4] = 14.;
	265	param->fMaxNormChi2NonC[5] = 25.;
	266
	267	param->fMaxNormChi2C[0] = 11.;
	268	param->fMaxNormChi2C[1] = 13.;
	269	param->fMaxNormChi2C[2] = 15.;
	270	param->fMaxNormChi2C[3] = 18.;
	271	param->fMaxNormChi2C[4] = 30.;
	272	param->fMaxNormChi2C[5] = 35.;
	273
	274	param->fMaxNormChi2NonCForHypothesis = 7.;
	275
	276	param->fMaxChi2 = 35.;
	277
	278	param->fMaxChi2s[0] = 25.; //40
	279	param->fMaxChi2s[1] = 25.; //40
	280	param->fMaxChi2s[2] = 25.; //40
	281	param->fMaxChi2s[3] = 25.; //40
	282	param->fMaxChi2s[4] = 40.; //40
	283	param->fMaxChi2s[5] = 50.; //40
	284
	285	param->fMaxRoad = 6.;
	286
	287	// not used
	288	param->fMaxChi2In = 16.;
	289
	290	param->fMaxChi2sR[0] = 10.;
	291	param->fMaxChi2sR[1] = 10.;
	292	param->fMaxChi2sR[2] = 10.;
	293	param->fMaxChi2sR[3] = 10.;
	294	param->fMaxChi2sR[4] = 30.;
	295	param->fMaxChi2sR[5] = 40.;
	296
	297	param->fChi2PerCluster = 9.;
	298	// not used
	299
	300	param->fXV = 0.;
	301	param->fYV = 0.;
	302	param->fZV = 0.;
	303	param->fSigmaXV = 0.0050;
	304	param->fSigmaYV = 0.0050;
305	param->fSigmaZV = 0.0100;
306
307	param->fVertexCut = 25.;
308
309	param->fMaxDZforPrimTrk = 0.4;
310	param->fMaxDZToUseConstraint = 3.;
311
312	param->fMaxDforV0dghtrForProlongation = 30.;
313	param->fMaxDForProlongation = 10.;
314	param->fMaxDZForProlongation = 20.;
315	param->fMinPtForProlongation = 0.120;
316
317	param->fZWindowDeadZone = 2.0;
318	param->fSigmaXDeadZoneHit2 = 0.004/12.;
319	param->fSigmaZDeadZoneHit2 = 0.001/12.;
320	param->fXPassDeadZoneHits = 0.018;
321
322
323	return param;
44347160	324	}
44347160	325	//_____________________________________________________________________________
ae00569a	326	AliITSRecoParam *AliITSRecoParam::GetPlaneEffParam(Int_t i)
	327	{
	328	//
	329	// make special reconstruction parameters for Plane Efficiency study on layer i
	330	//
	331	AliITSRecoParam *param;
	332	param = GetHighFluxParam();
	333	param->SetComputePlaneEff();
	334	param->SetLayerToSkip(i);
	335	return param;
	336	}
	337	//_____________________________________________________________________________
44347160	338	void AliITSRecoParam::SetLayersParameters()
	339	{
	340	//
	341	// number of layers and layers spatial resolutions
	342	//
	343
	344	// spatial resolutions of the detectors
	345	// y: 12 12 38 38 20 20 micron
	346	fSigmaY2[0]=1.44e-6;
	347	fSigmaY2[1]=1.44e-6;
	348	fSigmaY2[2]=1.444e-5;
	349	fSigmaY2[3]=1.444e-5;
	350	fSigmaY2[4]=4.0e-6;
	351	fSigmaY2[5]=4.0e-6;
	352	// z: 120 120 28 28 830 830 micron
	353	fSigmaZ2[0]=1.44e-4;
	354	fSigmaZ2[1]=1.44e-4;
	355	fSigmaZ2[2]=7.84e-6;
	356	fSigmaZ2[3]=7.84e-6;
	357	fSigmaZ2[4]=6.889e-3;
	358	fSigmaZ2[5]=6.889e-3;
	359
	360	return;
	361	}
ae00569a	362