[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),
fHistoPlaneEff(kFALSE),
fExtendedEtaAcceptance(kFALSE),
fUseDeadZonesFromOCDB(kFALSE),
fAllowProlongationWithEmptyRoad(kFALSE),
fFactorSAWindowSizes(1.),
fSAOnePointTracks(kFALSE),
fSAUseAllClusters(kFALSE),
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),
5fbd4fd6	101	fHistoPlaneEff(kFALSE),
2755f080	102	fExtendedEtaAcceptance(kFALSE),
6518a6c5	103	fUseDeadZonesFromOCDB(kFALSE),
ae00569a	104	fAllowProlongationWithEmptyRoad(kFALSE),
2755f080	105	fFactorSAWindowSizes(1.),
5a03f353	106	fSAOnePointTracks(kFALSE),
5a03f353	107	fSAUseAllClusters(kFALSE),
6518a6c5	108	fFindV0s(kTRUE),
	109	fUseUnfoldingInClusterFinderSPD(kFALSE),
	110	fUseUnfoldingInClusterFinderSDD(kTRUE),
	111	fUseUnfoldingInClusterFinderSSD(kTRUE)
44347160	112	{
	113	//
	114	// constructor
	115	//
6518a6c5	116	SetName("ITS");
	117	SetTitle("ITS");
	118
44347160	119	SetLayersParameters();
2755f080	120	for(Int_t i=0; i<AliITSgeomTGeo::kNLayers; i++) fLayersToSkip[i]=0;
6518a6c5	121	SetUseTGeoInTracker(3);
afd25725	122	SetAllowSharedClusters(kTRUE);
afd25725	123	SetFindV0s(kTRUE);
6518a6c5	124	SetAddVirtualClustersInDeadZone(kFALSE);
572f41f9	125	SetUseAmplitudeInfo(kTRUE);
e50912db	126	SetClusterErrorsParam(1);
44347160	127	}
	128	//_____________________________________________________________________________
	129	AliITSRecoParam::~AliITSRecoParam()
	130	{
	131	//
	132	// destructor
	133	//
	134	}
	135	//_____________________________________________________________________________
	136	AliITSRecoParam *AliITSRecoParam::GetHighFluxParam()
	137	{
	138	//
	139	// make default reconstruction parameters for hig flux env.
	140	//
	141	AliITSRecoParam *param = new AliITSRecoParam();
	142
	143	param->fMaxSnp = 0.95;
	144
	145	param->fNSigmaYLayerForRoadY = 4.;
	146	param->fNSigmaRoadY = 7.5;
	147	param->fNSigmaZLayerForRoadZ = 4.;
	148	param->fNSigmaRoadZ = 7.5;
	149
	150	param->fNSigma2RoadZC = 60.; //7.75^2
	151	param->fNSigma2RoadYC = 60.; //7.75^2
	152	param->fNSigma2RoadZNonC = 50.; //7.07^2
	153	param->fNSigma2RoadYNonC = 50.; //7.07^2
	154
	155	param->fMaxChi2PerCluster[0] = 11.; //7
	156	param->fMaxChi2PerCluster[1] = 12.; //5
	157	param->fMaxChi2PerCluster[2] = 12.; //8
	158	param->fMaxChi2PerCluster[3] = 5.; //8
	159	param->fMaxChi2PerCluster[4] = 12.; //6.5
	160
	161	param->fMaxNormChi2NonC[0] = 7.;
	162	param->fMaxNormChi2NonC[1] = 8.;
	163	param->fMaxNormChi2NonC[2] = 8.;
	164	param->fMaxNormChi2NonC[3] = 11.;
	165	param->fMaxNormChi2NonC[4] = 14.;
	166	param->fMaxNormChi2NonC[5] = 25.;
	167
	168	param->fMaxNormChi2C[0] = 11.;
	169	param->fMaxNormChi2C[1] = 13.;
	170	param->fMaxNormChi2C[2] = 15.;
	171	param->fMaxNormChi2C[3] = 18.;
	172	param->fMaxNormChi2C[4] = 30.;
	173	param->fMaxNormChi2C[5] = 35.;
afd25725	174
afd25725	175	param->fMaxNormChi2NonCForHypothesis = 7.;
44347160	176
	177	param->fMaxChi2 = 35.;
	178
	179	param->fMaxChi2s[0] = 25.; //40
	180	param->fMaxChi2s[1] = 25.; //40
	181	param->fMaxChi2s[2] = 25.; //40
	182	param->fMaxChi2s[3] = 25.; //40
	183	param->fMaxChi2s[4] = 40.; //40
	184	param->fMaxChi2s[5] = 50.; //40
	185
	186	param->fMaxRoad = 6.;
	187
	188	// not used
	189	param->fMaxChi2In = 16.;
	190
	191	param->fMaxChi2sR[0] = 10.;
	192	param->fMaxChi2sR[1] = 10.;
	193	param->fMaxChi2sR[2] = 10.;
	194	param->fMaxChi2sR[3] = 10.;
	195	param->fMaxChi2sR[4] = 30.;
	196	param->fMaxChi2sR[5] = 40.;
	197
	198	param->fChi2PerCluster = 9.;
	199	// not used
	200
	201	param->fXV = 0.;
	202	param->fYV = 0.;
	203	param->fZV = 0.;
	204	param->fSigmaXV = 0.0050;
	205	param->fSigmaYV = 0.0050;
	206	param->fSigmaZV = 0.0100;
	207
	208	param->fVertexCut = 25.;
afd25725	209
	210	param->fMaxDZforPrimTrk = 0.4;
	211	param->fMaxDZToUseConstraint = 3.;
	212
	213	param->fMaxDforV0dghtrForProlongation = 30.;
6518a6c5	214	param->fMaxDForProlongation = 40.;//10.;
6518a6c5	215	param->fMaxDZForProlongation = 60.;//20.;
afd25725	216	param->fMinPtForProlongation = 0.120;
	217
	218	param->fZWindowDeadZone = 2.0;
	219	param->fSigmaXDeadZoneHit2 = 0.004/12.;
	220	param->fSigmaZDeadZoneHit2 = 0.001/12.;
	221	param->fXPassDeadZoneHits = 0.018;
	222
44347160	223
	224	return param;
	225	}
	226	//_____________________________________________________________________________
	227	AliITSRecoParam *AliITSRecoParam::GetLowFluxParam()
	228	{
	229	//
	230	// make default reconstruction parameters for low flux env.
	231	//
	232	return GetHighFluxParam();
	233	}
	234	//_____________________________________________________________________________
	235	AliITSRecoParam *AliITSRecoParam::GetCosmicTestParam()
	236	{
	237	//
	238	// make default reconstruction parameters for cosmics
	239	//
e340bb86	240	AliITSRecoParam *param = new AliITSRecoParam();
e340bb86	241
2755f080	242	// larger seach windows for SA (in case of large misalignments)
2755f080	243	param->SetFactorSAWindowSizes(3.);
e340bb86	244
	245	param->fMaxSnp = 0.95;
	246
	247	param->fNSigmaYLayerForRoadY = 4.;
	248	param->fNSigmaRoadY = 7.5;
	249	param->fNSigmaZLayerForRoadZ = 4.;
	250	param->fNSigmaRoadZ = 7.5;
	251
	252	param->fNSigma2RoadZC = 60.; //7.75^2
	253	param->fNSigma2RoadYC = 60.; //7.75^2
	254	param->fNSigma2RoadZNonC = 50.; //7.07^2
	255	param->fNSigma2RoadYNonC = 50.; //7.07^2
	256
	257	param->fMaxChi2PerCluster[0] = 11.; //7
	258	param->fMaxChi2PerCluster[1] = 12.; //5
	259	param->fMaxChi2PerCluster[2] = 12.; //8
	260	param->fMaxChi2PerCluster[3] = 5.; //8
	261	param->fMaxChi2PerCluster[4] = 12.; //6.5
	262
	263	param->fMaxNormChi2NonC[0] = 7.;
	264	param->fMaxNormChi2NonC[1] = 8.;
	265	param->fMaxNormChi2NonC[2] = 8.;
	266	param->fMaxNormChi2NonC[3] = 11.;
	267	param->fMaxNormChi2NonC[4] = 14.;
	268	param->fMaxNormChi2NonC[5] = 25.;
	269
	270	param->fMaxNormChi2C[0] = 11.;
	271	param->fMaxNormChi2C[1] = 13.;
	272	param->fMaxNormChi2C[2] = 15.;
	273	param->fMaxNormChi2C[3] = 18.;
	274	param->fMaxNormChi2C[4] = 30.;
	275	param->fMaxNormChi2C[5] = 35.;
	276
	277	param->fMaxNormChi2NonCForHypothesis = 7.;
	278
	279	param->fMaxChi2 = 35.;
	280
	281	param->fMaxChi2s[0] = 25.; //40
	282	param->fMaxChi2s[1] = 25.; //40
	283	param->fMaxChi2s[2] = 25.; //40
	284	param->fMaxChi2s[3] = 25.; //40
	285	param->fMaxChi2s[4] = 40.; //40
	286	param->fMaxChi2s[5] = 50.; //40
	287
	288	param->fMaxRoad = 6.;
	289
	290	// not used
	291	param->fMaxChi2In = 16.;
	292
	293	param->fMaxChi2sR[0] = 10.;
	294	param->fMaxChi2sR[1] = 10.;
	295	param->fMaxChi2sR[2] = 10.;
	296	param->fMaxChi2sR[3] = 10.;
	297	param->fMaxChi2sR[4] = 30.;
	298	param->fMaxChi2sR[5] = 40.;
	299
	300	param->fChi2PerCluster = 9.;
	301	// not used
	302
	303	param->fXV = 0.;
	304	param->fYV = 0.;
	305	param->fZV = 0.;
	306	param->fSigmaXV = 0.0050;
	307	param->fSigmaYV = 0.0050;
308	param->fSigmaZV = 0.0100;
309
310	param->fVertexCut = 25.;
311
312	param->fMaxDZforPrimTrk = 0.4;
313	param->fMaxDZToUseConstraint = 3.;
314
315	param->fMaxDforV0dghtrForProlongation = 30.;
316	param->fMaxDForProlongation = 10.;
317	param->fMaxDZForProlongation = 20.;
318	param->fMinPtForProlongation = 0.120;
319
320	param->fZWindowDeadZone = 2.0;
321	param->fSigmaXDeadZoneHit2 = 0.004/12.;
322	param->fSigmaZDeadZoneHit2 = 0.001/12.;
323	param->fXPassDeadZoneHits = 0.018;
324
325
326	return param;
44347160	327	}
44347160	328	//_____________________________________________________________________________
ae00569a	329	AliITSRecoParam *AliITSRecoParam::GetPlaneEffParam(Int_t i)
	330	{
	331	//
	332	// make special reconstruction parameters for Plane Efficiency study on layer i
	333	//
	334	AliITSRecoParam *param;
	335	param = GetHighFluxParam();
	336	param->SetComputePlaneEff();
	337	param->SetLayerToSkip(i);
	338	return param;
	339	}
	340	//_____________________________________________________________________________
44347160	341	void AliITSRecoParam::SetLayersParameters()
	342	{
	343	//
	344	// number of layers and layers spatial resolutions
	345	//
	346
	347	// spatial resolutions of the detectors
	348	// y: 12 12 38 38 20 20 micron
	349	fSigmaY2[0]=1.44e-6;
	350	fSigmaY2[1]=1.44e-6;
	351	fSigmaY2[2]=1.444e-5;
	352	fSigmaY2[3]=1.444e-5;
	353	fSigmaY2[4]=4.0e-6;
	354	fSigmaY2[5]=4.0e-6;
	355	// z: 120 120 28 28 830 830 micron
	356	fSigmaZ2[0]=1.44e-4;
	357	fSigmaZ2[1]=1.44e-4;
	358	fSigmaZ2[2]=7.84e-6;
	359	fSigmaZ2[3]=7.84e-6;
	360	fSigmaZ2[4]=6.889e-3;
	361	fSigmaZ2[5]=6.889e-3;
	362
	363	return;
	364	}
ae00569a	365