[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() : TObject(),
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(kTRUE),
fZWindowDeadZone(0),
fSigmaXDeadZoneHit2(0),
fSigmaZDeadZoneHit2(0),
fXPassDeadZoneHits(0),
fUseTGeoInTracker(0),
fAllowSharedClusters(kTRUE),
fClusterErrorsParam(1),
fFindV0s(kTRUE)
{
  //
  // constructor
  //
  SetLayersParameters();
  SetUseTGeoInTracker(0);
  SetAllowSharedClusters(kTRUE);
  SetFindV0s(kTRUE);
  SetAddVirtualClustersInDeadZone(kTRUE);
  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 = 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::GetLowFluxParam() 
{
  //
  // make default reconstruction  parameters for low  flux env.
  //
  return GetHighFluxParam();
}
//_____________________________________________________________________________
AliITSRecoParam *AliITSRecoParam::GetCosmicTestParam() 
{
  //
  // make default reconstruction  parameters for cosmics
  //
  return GetHighFluxParam();
}
//_____________________________________________________________________________
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	//_____________________________________________________________________________
572f41f9	64	AliITSRecoParam::AliITSRecoParam() : TObject(),
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),
	92	fAddVirtualClustersInDeadZone(kTRUE),
	93	fZWindowDeadZone(0),
	94	fSigmaXDeadZoneHit2(0),
	95	fSigmaZDeadZoneHit2(0),
	96	fXPassDeadZoneHits(0),
e50912db	97	fUseTGeoInTracker(0),
572f41f9	98	fAllowSharedClusters(kTRUE),
e50912db	99	fClusterErrorsParam(1),
572f41f9	100	fFindV0s(kTRUE)
44347160	101	{
	102	//
	103	// constructor
	104	//
	105	SetLayersParameters();
e50912db	106	SetUseTGeoInTracker(0);
afd25725	107	SetAllowSharedClusters(kTRUE);
	108	SetFindV0s(kTRUE);
	109	SetAddVirtualClustersInDeadZone(kTRUE);
572f41f9	110	SetUseAmplitudeInfo(kTRUE);
e50912db	111	SetClusterErrorsParam(1);
44347160	112	}
	113	//_____________________________________________________________________________
	114	AliITSRecoParam::~AliITSRecoParam()
	115	{
	116	//
	117	// destructor
	118	//
	119	}
	120	//_____________________________________________________________________________
	121	AliITSRecoParam *AliITSRecoParam::GetHighFluxParam()
	122	{
	123	//
	124	// make default reconstruction parameters for hig flux env.
	125	//
	126	AliITSRecoParam *param = new AliITSRecoParam();
	127
	128	param->fMaxSnp = 0.95;
	129
	130	param->fNSigmaYLayerForRoadY = 4.;
	131	param->fNSigmaRoadY = 7.5;
	132	param->fNSigmaZLayerForRoadZ = 4.;
	133	param->fNSigmaRoadZ = 7.5;
	134
	135	param->fNSigma2RoadZC = 60.; //7.75^2
	136	param->fNSigma2RoadYC = 60.; //7.75^2
	137	param->fNSigma2RoadZNonC = 50.; //7.07^2
	138	param->fNSigma2RoadYNonC = 50.; //7.07^2
	139
	140	param->fMaxChi2PerCluster[0] = 11.; //7
	141	param->fMaxChi2PerCluster[1] = 12.; //5
	142	param->fMaxChi2PerCluster[2] = 12.; //8
	143	param->fMaxChi2PerCluster[3] = 5.; //8
	144	param->fMaxChi2PerCluster[4] = 12.; //6.5
	145
	146	param->fMaxNormChi2NonC[0] = 7.;
	147	param->fMaxNormChi2NonC[1] = 8.;
	148	param->fMaxNormChi2NonC[2] = 8.;
	149	param->fMaxNormChi2NonC[3] = 11.;
	150	param->fMaxNormChi2NonC[4] = 14.;
	151	param->fMaxNormChi2NonC[5] = 25.;
	152
	153	param->fMaxNormChi2C[0] = 11.;
	154	param->fMaxNormChi2C[1] = 13.;
	155	param->fMaxNormChi2C[2] = 15.;
	156	param->fMaxNormChi2C[3] = 18.;
	157	param->fMaxNormChi2C[4] = 30.;
	158	param->fMaxNormChi2C[5] = 35.;
afd25725	159
afd25725	160	param->fMaxNormChi2NonCForHypothesis = 7.;
44347160	161
	162	param->fMaxChi2 = 35.;
	163
	164	param->fMaxChi2s[0] = 25.; //40
	165	param->fMaxChi2s[1] = 25.; //40
	166	param->fMaxChi2s[2] = 25.; //40
	167	param->fMaxChi2s[3] = 25.; //40
	168	param->fMaxChi2s[4] = 40.; //40
	169	param->fMaxChi2s[5] = 50.; //40
	170
	171	param->fMaxRoad = 6.;
	172
	173	// not used
	174	param->fMaxChi2In = 16.;
	175
	176	param->fMaxChi2sR[0] = 10.;
	177	param->fMaxChi2sR[1] = 10.;
	178	param->fMaxChi2sR[2] = 10.;
	179	param->fMaxChi2sR[3] = 10.;
	180	param->fMaxChi2sR[4] = 30.;
	181	param->fMaxChi2sR[5] = 40.;
	182
	183	param->fChi2PerCluster = 9.;
	184	// not used
	185
	186	param->fXV = 0.;
	187	param->fYV = 0.;
	188	param->fZV = 0.;
	189	param->fSigmaXV = 0.0050;
	190	param->fSigmaYV = 0.0050;
	191	param->fSigmaZV = 0.0100;
	192
	193	param->fVertexCut = 25.;
afd25725	194
	195	param->fMaxDZforPrimTrk = 0.4;
	196	param->fMaxDZToUseConstraint = 3.;
	197
	198	param->fMaxDforV0dghtrForProlongation = 30.;
	199	param->fMaxDForProlongation = 10.;
	200	param->fMaxDZForProlongation = 20.;
	201	param->fMinPtForProlongation = 0.120;
	202
	203	param->fZWindowDeadZone = 2.0;
	204	param->fSigmaXDeadZoneHit2 = 0.004/12.;
	205	param->fSigmaZDeadZoneHit2 = 0.001/12.;
	206	param->fXPassDeadZoneHits = 0.018;
	207
44347160	208
	209	return param;
	210	}
	211	//_____________________________________________________________________________
	212	AliITSRecoParam *AliITSRecoParam::GetLowFluxParam()
	213	{
	214	//
	215	// make default reconstruction parameters for low flux env.
	216	//
	217	return GetHighFluxParam();
	218	}
	219	//_____________________________________________________________________________
	220	AliITSRecoParam *AliITSRecoParam::GetCosmicTestParam()
	221	{
	222	//
	223	// make default reconstruction parameters for cosmics
	224	//
	225	return GetHighFluxParam();
	226	}
	227	//_____________________________________________________________________________
	228	void AliITSRecoParam::SetLayersParameters()
	229	{
	230	//
	231	// number of layers and layers spatial resolutions
	232	//
	233
	234	// spatial resolutions of the detectors
	235	// y: 12 12 38 38 20 20 micron
	236	fSigmaY2[0]=1.44e-6;
	237	fSigmaY2[1]=1.44e-6;
	238	fSigmaY2[2]=1.444e-5;
	239	fSigmaY2[3]=1.444e-5;
	240	fSigmaY2[4]=4.0e-6;
	241	fSigmaY2[5]=4.0e-6;
	242	// z: 120 120 28 28 830 830 micron
	243	fSigmaZ2[0]=1.44e-4;
	244	fSigmaZ2[1]=1.44e-4;
	245	fSigmaZ2[2]=7.84e-6;
	246	fSigmaZ2[3]=7.84e-6;
	247	fSigmaZ2[4]=6.889e-3;
	248	fSigmaZ2[5]=6.889e-3;
	249
	250	return;
	251	}
	252