[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)


//_____________________________________________________________________________
AliITSRecoParam::AliITSRecoParam() : TObject(),
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(kFALSE),
fAllowSharedClusters(kTRUE),
fUseNominalClusterErrors(kFALSE),
fFindV0s(kTRUE)
{
  //
  // constructor
  //
  SetLayersParameters();
  SetUseTGeoInTracker(kFALSE);
  SetAllowSharedClusters(kTRUE);
  SetFindV0s(kTRUE);
  SetAddVirtualClustersInDeadZone(kTRUE);
  SetUseAmplitudeInfo(kTRUE);
}
//_____________________________________________________________________________
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
	31
	32	//_____________________________________________________________________________
572f41f9	33	AliITSRecoParam::AliITSRecoParam() : TObject(),
	34	fNSigmaYLayerForRoadY(0),
	35	fNSigmaRoadY(0),
	36	fNSigmaZLayerForRoadZ(0),
	37	fNSigmaRoadZ(0),
	38	fNSigma2RoadZC(0),
	39	fNSigma2RoadYC(0),
	40	fNSigma2RoadZNonC(0),
	41	fNSigma2RoadYNonC(0),
	42	fMaxNormChi2NonCForHypothesis(0),
	43	fMaxChi2(0),
	44	fMaxRoad(0),
	45	fMaxChi2In(0),
	46	fChi2PerCluster(0),
	47	fXV(0),
	48	fYV(0),
	49	fZV(0),
	50	fSigmaXV(0),
	51	fSigmaYV(0),
	52	fSigmaZV(0),
	53	fVertexCut(0),
	54	fMaxDZforPrimTrk(0),
	55	fMaxDZToUseConstraint(0),
	56	fMaxDforV0dghtrForProlongation(0),
	57	fMaxDForProlongation(0),
	58	fMaxDZForProlongation(0),
	59	fMinPtForProlongation(0),
	60	fAddVirtualClustersInDeadZone(kTRUE),
	61	fZWindowDeadZone(0),
	62	fSigmaXDeadZoneHit2(0),
	63	fSigmaZDeadZoneHit2(0),
	64	fXPassDeadZoneHits(0),
	65	fUseTGeoInTracker(kFALSE),
	66	fAllowSharedClusters(kTRUE),
	67	fUseNominalClusterErrors(kFALSE),
	68	fFindV0s(kTRUE)
44347160	69	{
	70	//
	71	// constructor
	72	//
	73	SetLayersParameters();
afd25725	74	SetUseTGeoInTracker(kFALSE);
	75	SetAllowSharedClusters(kTRUE);
	76	SetFindV0s(kTRUE);
	77	SetAddVirtualClustersInDeadZone(kTRUE);
572f41f9	78	SetUseAmplitudeInfo(kTRUE);
44347160	79	}
	80	//_____________________________________________________________________________
	81	AliITSRecoParam::~AliITSRecoParam()
	82	{
	83	//
	84	// destructor
	85	//
	86	}
	87	//_____________________________________________________________________________
	88	AliITSRecoParam *AliITSRecoParam::GetHighFluxParam()
	89	{
	90	//
	91	// make default reconstruction parameters for hig flux env.
	92	//
	93	AliITSRecoParam *param = new AliITSRecoParam();
	94
	95	param->fMaxSnp = 0.95;
	96
	97	param->fNSigmaYLayerForRoadY = 4.;
	98	param->fNSigmaRoadY = 7.5;
	99	param->fNSigmaZLayerForRoadZ = 4.;
	100	param->fNSigmaRoadZ = 7.5;
	101
	102	param->fNSigma2RoadZC = 60.; //7.75^2
	103	param->fNSigma2RoadYC = 60.; //7.75^2
	104	param->fNSigma2RoadZNonC = 50.; //7.07^2
	105	param->fNSigma2RoadYNonC = 50.; //7.07^2
	106
	107	param->fMaxChi2PerCluster[0] = 11.; //7
	108	param->fMaxChi2PerCluster[1] = 12.; //5
	109	param->fMaxChi2PerCluster[2] = 12.; //8
	110	param->fMaxChi2PerCluster[3] = 5.; //8
	111	param->fMaxChi2PerCluster[4] = 12.; //6.5
	112
	113	param->fMaxNormChi2NonC[0] = 7.;
	114	param->fMaxNormChi2NonC[1] = 8.;
	115	param->fMaxNormChi2NonC[2] = 8.;
	116	param->fMaxNormChi2NonC[3] = 11.;
	117	param->fMaxNormChi2NonC[4] = 14.;
	118	param->fMaxNormChi2NonC[5] = 25.;
	119
	120	param->fMaxNormChi2C[0] = 11.;
	121	param->fMaxNormChi2C[1] = 13.;
	122	param->fMaxNormChi2C[2] = 15.;
	123	param->fMaxNormChi2C[3] = 18.;
	124	param->fMaxNormChi2C[4] = 30.;
	125	param->fMaxNormChi2C[5] = 35.;
afd25725	126
afd25725	127	param->fMaxNormChi2NonCForHypothesis = 7.;
44347160	128
	129	param->fMaxChi2 = 35.;
	130
	131	param->fMaxChi2s[0] = 25.; //40
	132	param->fMaxChi2s[1] = 25.; //40
	133	param->fMaxChi2s[2] = 25.; //40
	134	param->fMaxChi2s[3] = 25.; //40
	135	param->fMaxChi2s[4] = 40.; //40
	136	param->fMaxChi2s[5] = 50.; //40
	137
	138	param->fMaxRoad = 6.;
	139
	140	// not used
	141	param->fMaxChi2In = 16.;
	142
	143	param->fMaxChi2sR[0] = 10.;
	144	param->fMaxChi2sR[1] = 10.;
	145	param->fMaxChi2sR[2] = 10.;
	146	param->fMaxChi2sR[3] = 10.;
	147	param->fMaxChi2sR[4] = 30.;
	148	param->fMaxChi2sR[5] = 40.;
	149
	150	param->fChi2PerCluster = 9.;
	151	// not used
	152
	153	param->fXV = 0.;
	154	param->fYV = 0.;
	155	param->fZV = 0.;
	156	param->fSigmaXV = 0.0050;
	157	param->fSigmaYV = 0.0050;
	158	param->fSigmaZV = 0.0100;
	159
	160	param->fVertexCut = 25.;
afd25725	161
	162	param->fMaxDZforPrimTrk = 0.4;
	163	param->fMaxDZToUseConstraint = 3.;
	164
	165	param->fMaxDforV0dghtrForProlongation = 30.;
	166	param->fMaxDForProlongation = 10.;
	167	param->fMaxDZForProlongation = 20.;
	168	param->fMinPtForProlongation = 0.120;
	169
	170	param->fZWindowDeadZone = 2.0;
	171	param->fSigmaXDeadZoneHit2 = 0.004/12.;
	172	param->fSigmaZDeadZoneHit2 = 0.001/12.;
	173	param->fXPassDeadZoneHits = 0.018;
	174
44347160	175
	176	return param;
	177	}
	178	//_____________________________________________________________________________
	179	AliITSRecoParam *AliITSRecoParam::GetLowFluxParam()
	180	{
	181	//
	182	// make default reconstruction parameters for low flux env.
	183	//
	184	return GetHighFluxParam();
	185	}
	186	//_____________________________________________________________________________
	187	AliITSRecoParam *AliITSRecoParam::GetCosmicTestParam()
	188	{
	189	//
	190	// make default reconstruction parameters for cosmics
	191	//
	192	return GetHighFluxParam();
	193	}
	194	//_____________________________________________________________________________
	195	void AliITSRecoParam::SetLayersParameters()
	196	{
	197	//
	198	// number of layers and layers spatial resolutions
	199	//
	200
	201	// spatial resolutions of the detectors
	202	// y: 12 12 38 38 20 20 micron
	203	fSigmaY2[0]=1.44e-6;
	204	fSigmaY2[1]=1.44e-6;
	205	fSigmaY2[2]=1.444e-5;
	206	fSigmaY2[3]=1.444e-5;
	207	fSigmaY2[4]=4.0e-6;
	208	fSigmaY2[5]=4.0e-6;
	209	// z: 120 120 28 28 830 830 micron
	210	fSigmaZ2[0]=1.44e-4;
	211	fSigmaZ2[1]=1.44e-4;
	212	fSigmaZ2[2]=7.84e-6;
	213	fSigmaZ2[3]=7.84e-6;
	214	fSigmaZ2[4]=6.889e-3;
	215	fSigmaZ2[5]=6.889e-3;
	216
	217	return;
	218	}
	219