[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),
fReadPlaneEffFromOCDB(kFALSE),
fExtendedEtaAcceptance(kFALSE),
fUseBadZonesFromOCDB(kFALSE),
fUseSingleBadChannelsFromOCDB(kFALSE),
fMinFractionOfBadInRoad(0),
fAllowProlongationWithEmptyRoad(kFALSE),
fFactorSAWindowSizes(1.),
fNLoopsSA(33),
fMinPhiSA(0.002),
fMaxPhiSA(0.0145),
fMinLambdaSA(0.003),
fMaxLambdaSA(0.008),
fSAOnePointTracks(kFALSE),
fSAUseAllClusters(kFALSE),
fFindV0s(kTRUE),
fUseUnfoldingInClusterFinderSPD(kFALSE),
fUseUnfoldingInClusterFinderSDD(kTRUE),
fUseUnfoldingInClusterFinderSSD(kTRUE),
fUseChargeMatchingInClusterFinderSSD(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.
  //
  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;
  param->SetNLoopsSA(10);
  
  return param;
}
//_____________________________________________________________________________
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;
}
//_____________________________________________________________________________
void AliITSRecoParam::PrintParameters() const 
{
  //
  // print parameters
  //

  printf("=============================  AliITSRecoParam::PrintParameters ");
  printf("============================= \n \n ");
  for(Int_t i=0; i<AliITSgeomTGeo::kNLayers; i++) {
    if(!fLayersToSkip[i]) {
      printf("ITS Traking: using layer %d\n",i);
    } else {
      printf("ITS Traking: skipping layer %d\n",i);
    }
  }
  for(Int_t i=0; i<AliITSgeomTGeo::kNLayers; i++) {
    if(fUseAmplitudeInfo[i]) {
      printf("ITS Traking: use amplitude info for layer %d\n",i);
    } else {
      printf("ITS Traking: don't use amplitude info for layer %d\n",i);
    }
  }
  for(Int_t i=0; i<AliITSgeomTGeo::kNLayers; i++)
    printf("Layer %d:\n  sigmaY2 %f, sigma Z2 %f\n  max norm chi2 for non constrained tracks %f\n  max norm chi2 for constrained tracks %f\n  max predicted chi2 (cluster & track prol.) %f\n",i,fSigmaY2[i],fSigmaZ2[i],fMaxNormChi2NonC[i],fMaxNormChi2C[i],fMaxChi2s[i]);


  Dump();

  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),
275a301c	102	fReadPlaneEffFromOCDB(kFALSE),
2755f080	103	fExtendedEtaAcceptance(kFALSE),
23197852	104	fUseBadZonesFromOCDB(kFALSE),
	105	fUseSingleBadChannelsFromOCDB(kFALSE),
	106	fMinFractionOfBadInRoad(0),
ae00569a	107	fAllowProlongationWithEmptyRoad(kFALSE),
2755f080	108	fFactorSAWindowSizes(1.),
c7d6d7b7	109	fNLoopsSA(33),
	110	fMinPhiSA(0.002),
	111	fMaxPhiSA(0.0145),
	112	fMinLambdaSA(0.003),
	113	fMaxLambdaSA(0.008),
5a03f353	114	fSAOnePointTracks(kFALSE),
5a03f353	115	fSAUseAllClusters(kFALSE),
6518a6c5	116	fFindV0s(kTRUE),
	117	fUseUnfoldingInClusterFinderSPD(kFALSE),
	118	fUseUnfoldingInClusterFinderSDD(kTRUE),
a86176e3	119	fUseUnfoldingInClusterFinderSSD(kTRUE),
a86176e3	120	fUseChargeMatchingInClusterFinderSSD(kTRUE)
44347160	121	{
	122	//
	123	// constructor
	124	//
6518a6c5	125	SetName("ITS");
	126	SetTitle("ITS");
	127
44347160	128	SetLayersParameters();
2755f080	129	for(Int_t i=0; i<AliITSgeomTGeo::kNLayers; i++) fLayersToSkip[i]=0;
6518a6c5	130	SetUseTGeoInTracker(3);
afd25725	131	SetAllowSharedClusters(kTRUE);
afd25725	132	SetFindV0s(kTRUE);
6518a6c5	133	SetAddVirtualClustersInDeadZone(kFALSE);
572f41f9	134	SetUseAmplitudeInfo(kTRUE);
e50912db	135	SetClusterErrorsParam(1);
44347160	136	}
	137	//_____________________________________________________________________________
	138	AliITSRecoParam::~AliITSRecoParam()
	139	{
	140	//
	141	// destructor
	142	//
	143	}
	144	//_____________________________________________________________________________
	145	AliITSRecoParam *AliITSRecoParam::GetHighFluxParam()
	146	{
	147	//
	148	// make default reconstruction parameters for hig flux env.
	149	//
	150	AliITSRecoParam *param = new AliITSRecoParam();
	151
	152	param->fMaxSnp = 0.95;
	153
	154	param->fNSigmaYLayerForRoadY = 4.;
	155	param->fNSigmaRoadY = 7.5;
	156	param->fNSigmaZLayerForRoadZ = 4.;
	157	param->fNSigmaRoadZ = 7.5;
	158
	159	param->fNSigma2RoadZC = 60.; //7.75^2
	160	param->fNSigma2RoadYC = 60.; //7.75^2
	161	param->fNSigma2RoadZNonC = 50.; //7.07^2
	162	param->fNSigma2RoadYNonC = 50.; //7.07^2
	163
	164	param->fMaxChi2PerCluster[0] = 11.; //7
	165	param->fMaxChi2PerCluster[1] = 12.; //5
	166	param->fMaxChi2PerCluster[2] = 12.; //8
	167	param->fMaxChi2PerCluster[3] = 5.; //8
	168	param->fMaxChi2PerCluster[4] = 12.; //6.5
	169
	170	param->fMaxNormChi2NonC[0] = 7.;
	171	param->fMaxNormChi2NonC[1] = 8.;
	172	param->fMaxNormChi2NonC[2] = 8.;
	173	param->fMaxNormChi2NonC[3] = 11.;
	174	param->fMaxNormChi2NonC[4] = 14.;
	175	param->fMaxNormChi2NonC[5] = 25.;
	176
	177	param->fMaxNormChi2C[0] = 11.;
	178	param->fMaxNormChi2C[1] = 13.;
	179	param->fMaxNormChi2C[2] = 15.;
	180	param->fMaxNormChi2C[3] = 18.;
	181	param->fMaxNormChi2C[4] = 30.;
	182	param->fMaxNormChi2C[5] = 35.;
afd25725	183
afd25725	184	param->fMaxNormChi2NonCForHypothesis = 7.;
44347160	185
	186	param->fMaxChi2 = 35.;
	187
	188	param->fMaxChi2s[0] = 25.; //40
	189	param->fMaxChi2s[1] = 25.; //40
	190	param->fMaxChi2s[2] = 25.; //40
	191	param->fMaxChi2s[3] = 25.; //40
	192	param->fMaxChi2s[4] = 40.; //40
	193	param->fMaxChi2s[5] = 50.; //40
	194
	195	param->fMaxRoad = 6.;
	196
	197	// not used
	198	param->fMaxChi2In = 16.;
	199
	200	param->fMaxChi2sR[0] = 10.;
	201	param->fMaxChi2sR[1] = 10.;
	202	param->fMaxChi2sR[2] = 10.;
	203	param->fMaxChi2sR[3] = 10.;
	204	param->fMaxChi2sR[4] = 30.;
	205	param->fMaxChi2sR[5] = 40.;
	206
	207	param->fChi2PerCluster = 9.;
	208	// not used
	209
	210	param->fXV = 0.;
	211	param->fYV = 0.;
	212	param->fZV = 0.;
	213	param->fSigmaXV = 0.0050;
	214	param->fSigmaYV = 0.0050;
	215	param->fSigmaZV = 0.0100;
	216
	217	param->fVertexCut = 25.;
afd25725	218
	219	param->fMaxDZforPrimTrk = 0.4;
	220	param->fMaxDZToUseConstraint = 3.;
	221
	222	param->fMaxDforV0dghtrForProlongation = 30.;
6518a6c5	223	param->fMaxDForProlongation = 40.;//10.;
6518a6c5	224	param->fMaxDZForProlongation = 60.;//20.;
afd25725	225	param->fMinPtForProlongation = 0.120;
	226
	227	param->fZWindowDeadZone = 2.0;
	228	param->fSigmaXDeadZoneHit2 = 0.004/12.;
	229	param->fSigmaZDeadZoneHit2 = 0.001/12.;
	230	param->fXPassDeadZoneHits = 0.018;
44347160	231
	232	return param;
	233	}
	234	//_____________________________________________________________________________
	235	AliITSRecoParam *AliITSRecoParam::GetLowFluxParam()
	236	{
	237	//
	238	// make default reconstruction parameters for low flux env.
	239	//
c7d6d7b7	240	AliITSRecoParam *param = new AliITSRecoParam();
	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 = 40.;//10.;
314	param->fMaxDZForProlongation = 60.;//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	param->SetNLoopsSA(10);
322
323	return param;
44347160	324	}
	325	//_____________________________________________________________________________
	326	AliITSRecoParam *AliITSRecoParam::GetCosmicTestParam()
	327	{
	328	//
	329	// make default reconstruction parameters for cosmics
	330	//
e340bb86	331	AliITSRecoParam *param = new AliITSRecoParam();
e340bb86	332
2755f080	333	// larger seach windows for SA (in case of large misalignments)
2755f080	334	param->SetFactorSAWindowSizes(3.);
e340bb86	335
	336	param->fMaxSnp = 0.95;
	337
	338	param->fNSigmaYLayerForRoadY = 4.;
	339	param->fNSigmaRoadY = 7.5;
	340	param->fNSigmaZLayerForRoadZ = 4.;
	341	param->fNSigmaRoadZ = 7.5;
	342
	343	param->fNSigma2RoadZC = 60.; //7.75^2
	344	param->fNSigma2RoadYC = 60.; //7.75^2
	345	param->fNSigma2RoadZNonC = 50.; //7.07^2
	346	param->fNSigma2RoadYNonC = 50.; //7.07^2
	347
	348	param->fMaxChi2PerCluster[0] = 11.; //7
	349	param->fMaxChi2PerCluster[1] = 12.; //5
	350	param->fMaxChi2PerCluster[2] = 12.; //8
	351	param->fMaxChi2PerCluster[3] = 5.; //8
	352	param->fMaxChi2PerCluster[4] = 12.; //6.5
	353
	354	param->fMaxNormChi2NonC[0] = 7.;
	355	param->fMaxNormChi2NonC[1] = 8.;
	356	param->fMaxNormChi2NonC[2] = 8.;
	357	param->fMaxNormChi2NonC[3] = 11.;
	358	param->fMaxNormChi2NonC[4] = 14.;
	359	param->fMaxNormChi2NonC[5] = 25.;
	360
	361	param->fMaxNormChi2C[0] = 11.;
	362	param->fMaxNormChi2C[1] = 13.;
	363	param->fMaxNormChi2C[2] = 15.;
	364	param->fMaxNormChi2C[3] = 18.;
	365	param->fMaxNormChi2C[4] = 30.;
	366	param->fMaxNormChi2C[5] = 35.;
	367
	368	param->fMaxNormChi2NonCForHypothesis = 7.;
	369
	370	param->fMaxChi2 = 35.;
	371
	372	param->fMaxChi2s[0] = 25.; //40
	373	param->fMaxChi2s[1] = 25.; //40
	374	param->fMaxChi2s[2] = 25.; //40
	375	param->fMaxChi2s[3] = 25.; //40
	376	param->fMaxChi2s[4] = 40.; //40
	377	param->fMaxChi2s[5] = 50.; //40
	378
	379	param->fMaxRoad = 6.;
	380
	381	// not used
	382	param->fMaxChi2In = 16.;
	383
	384	param->fMaxChi2sR[0] = 10.;
	385	param->fMaxChi2sR[1] = 10.;
	386	param->fMaxChi2sR[2] = 10.;
	387	param->fMaxChi2sR[3] = 10.;
	388	param->fMaxChi2sR[4] = 30.;
	389	param->fMaxChi2sR[5] = 40.;
	390
	391	param->fChi2PerCluster = 9.;
	392	// not used
	393
	394	param->fXV = 0.;
	395	param->fYV = 0.;
	396	param->fZV = 0.;
	397	param->fSigmaXV = 0.0050;
	398	param->fSigmaYV = 0.0050;
399	param->fSigmaZV = 0.0100;
400
401	param->fVertexCut = 25.;
402
403	param->fMaxDZforPrimTrk = 0.4;
404	param->fMaxDZToUseConstraint = 3.;
405
406	param->fMaxDforV0dghtrForProlongation = 30.;
407	param->fMaxDForProlongation = 10.;
408	param->fMaxDZForProlongation = 20.;
409	param->fMinPtForProlongation = 0.120;
410
411	param->fZWindowDeadZone = 2.0;
412	param->fSigmaXDeadZoneHit2 = 0.004/12.;
413	param->fSigmaZDeadZoneHit2 = 0.001/12.;
414	param->fXPassDeadZoneHits = 0.018;
415
416
417	return param;
44347160	418	}
44347160	419	//_____________________________________________________________________________
ae00569a	420	AliITSRecoParam *AliITSRecoParam::GetPlaneEffParam(Int_t i)
	421	{
	422	//
	423	// make special reconstruction parameters for Plane Efficiency study on layer i
	424	//
	425	AliITSRecoParam *param;
	426	param = GetHighFluxParam();
	427	param->SetComputePlaneEff();
	428	param->SetLayerToSkip(i);
	429	return param;
	430	}
	431	//_____________________________________________________________________________
44347160	432	void AliITSRecoParam::SetLayersParameters()
	433	{
	434	//
	435	// number of layers and layers spatial resolutions
	436	//
	437
	438	// spatial resolutions of the detectors
	439	// y: 12 12 38 38 20 20 micron
	440	fSigmaY2[0]=1.44e-6;
	441	fSigmaY2[1]=1.44e-6;
	442	fSigmaY2[2]=1.444e-5;
	443	fSigmaY2[3]=1.444e-5;
	444	fSigmaY2[4]=4.0e-6;
	445	fSigmaY2[5]=4.0e-6;
	446	// z: 120 120 28 28 830 830 micron
	447	fSigmaZ2[0]=1.44e-4;
	448	fSigmaZ2[1]=1.44e-4;
	449	fSigmaZ2[2]=7.84e-6;
	450	fSigmaZ2[3]=7.84e-6;
	451	fSigmaZ2[4]=6.889e-3;
	452	fSigmaZ2[5]=6.889e-3;
	453
	454	return;
	455	}
ed446fa3	456	//_____________________________________________________________________________
	457	void AliITSRecoParam::PrintParameters() const
	458	{
	459	//
	460	// print parameters
	461	//
	462
	463	printf("============================= AliITSRecoParam::PrintParameters ");
	464	printf("============================= \n \n ");
	465	for(Int_t i=0; i<AliITSgeomTGeo::kNLayers; i++) {
	466	if(!fLayersToSkip[i]) {
	467	printf("ITS Traking: using layer %d\n",i);
	468	} else {
	469	printf("ITS Traking: skipping layer %d\n",i);
	470	}
	471	}
	472	for(Int_t i=0; i<AliITSgeomTGeo::kNLayers; i++) {
	473	if(fUseAmplitudeInfo[i]) {
	474	printf("ITS Traking: use amplitude info for layer %d\n",i);
	475	} else {
	476	printf("ITS Traking: don't use amplitude info for layer %d\n",i);
	477	}
	478	}
	479	for(Int_t i=0; i<AliITSgeomTGeo::kNLayers; i++)
	480	printf("Layer %d:\n sigmaY2 %f, sigma Z2 %f\n max norm chi2 for non constrained tracks %f\n max norm chi2 for constrained tracks %f\n max predicted chi2 (cluster & track prol.) %f\n",i,fSigmaY2[i],fSigmaZ2[i],fMaxNormChi2NonC[i],fMaxNormChi2C[i],fMaxChi2s[i]);
	481
	482
	483	Dump();
	484
	485	return;
	486	}