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