[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(),
fTracker(0),
fITSonly(kFALSE),
fVertexer(0),
fClusterFinder(0),
fPID(0),
fVtxr3DZCutWide(0.),
fVtxr3DRCutWide(0.),
fVtxr3DZCutNarrow(0.),
fVtxr3DRCutNarrow(0.),
fVtxr3DPhiCutLoose(0.),
fVtxr3DPhiCutTight(0.),
fVtxr3DDCACut(0.),
fMaxSnp(1.),
fNSigmaYLayerForRoadY(0),
fNSigmaRoadY(0),
fNSigmaZLayerForRoadZ(0),
fNSigmaRoadZ(0),
fNSigma2RoadZC(0),
fNSigma2RoadYC(0),
fNSigma2RoadZNonC(0),
fNSigma2RoadYNonC(0),
fRoadMisal(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),
fIPlanePlaneEff(0),
fReadPlaneEffFromOCDB(kFALSE),
fMinPtPlaneEff(0),
fMaxMissingClustersPlaneEff(0),
fRequireClusterInOuterLayerPlaneEff(kFALSE),
fRequireClusterInInnerLayerPlaneEff(kFALSE),
fOnlyConstraintPlaneEff(kFALSE),
fExtendedEtaAcceptance(kFALSE),
fUseBadZonesFromOCDB(kTRUE),
fUseSingleBadChannelsFromOCDB(kFALSE),
fMinFractionOfBadInRoad(0),
fAllowProlongationWithEmptyRoad(kFALSE),
fOuterStartLayerSA(0),
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),
fTrackleterOnlyOneTrackletPerC2(kTRUE),
fTrackleterPhiWindow(0.08),
fTrackleterZetaWindow(1.00),
fTrackleterRemoveClustersFromOverlaps(kFALSE),
fTrackleterPhiOverlapCut(0.005),
fTrackleterZetaOverlapCut(0.05),
fUseCosmicRunShiftsSSD(kFALSE)
{
  //
  // 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);
  SetClusterMisalError(0.);
  SetVertexer3DDefaults();
}
//_____________________________________________________________________________
AliITSRecoParam::~AliITSRecoParam() 
{
  //
  // destructor
  //  
}
//_____________________________________________________________________________
AliITSRecoParam *AliITSRecoParam::GetHighFluxParam() 
{
  //
  // make default reconstruction  parameters for hig  flux env.
  //
  AliITSRecoParam *param = new AliITSRecoParam();
  param->SetVertexerZ();

  // use of bads from OCDB
  param->SetUseBadZonesFromOCDB(kTRUE);
  param->SetUseSingleBadChannelsFromOCDB(kFALSE);
  // extended eta acceptance
  param->SetExtendedEtaAcceptance(kFALSE);
  // allow to skip layer if no cluster and no bad
  param->SetAllowProlongationWithEmptyRoad(kFALSE);

  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->SetVertexer3DDefaults();

  // full use of bads from OCDB
  param->SetUseBadZonesFromOCDB(kTRUE);
  param->SetUseSingleBadChannelsFromOCDB(kTRUE);
  // extended eta acceptance
  param->SetExtendedEtaAcceptance(kTRUE);
  // allow to skip layer if no cluster and no bad
  param->SetAllowProlongationWithEmptyRoad(kTRUE);

  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();

  // vertexer for cosmics
  param->SetVertexer(2);

  // find independently ITS SA tracks 
  param->SetSAUseAllClusters();
  param->SetOuterStartLayerSA(AliITSgeomTGeo::GetNLayers()-2);

  // 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
  //
  if (i<0 || i>=AliITSgeomTGeo::kNLayers) {
    printf("AliITSRecoParam::GetPlaneEffParam: index of ITS Plane not in the range [0,5]\n");
    printf("returning null pointer");
    return NULL;
  }
  AliITSRecoParam *param;
  param = GetHighFluxParam();
  param->SetComputePlaneEff();
  param->SetLayerToSkip(i);
  param->SetIPlanePlaneEff(i);
  // optimized setting for SPD0 (i==0)
  if (i==0 || i==1) {
    param->fMinPtPlaneEff = 0.200; // high pt particles
    param->fMaxMissingClustersPlaneEff = 1; // at most 1 layer out of 5 without cluster
    param->fRequireClusterInOuterLayerPlaneEff = kTRUE; // cluster on SPD1 must be
    //param->fOnlyConstraintPlaneEff = kTRUE;
  }
  if (i==2 || i==3) {
    param->fMinPtPlaneEff = 0.200; // high pt particles
    param->fMaxMissingClustersPlaneEff = 1; // at most 1 layer out of 5 without cluster
    param->fRequireClusterInOuterLayerPlaneEff = kTRUE;
    //param->fOnlyConstraintPlaneEff = kTRUE;
  }
  if (i==4) {
    param->fMinPtPlaneEff = 0.200; // high pt particles
    param->fMaxMissingClustersPlaneEff = 0; // at most 1 layer out of 5 without cluster
    param->fRequireClusterInOuterLayerPlaneEff = kTRUE;
    //param->fOnlyConstraintPlaneEff = kTRUE;
  }
  if (i==5) {
    param->fMinPtPlaneEff = 0.200; // high pt particles
  }
  //
  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, sigmaZ2 %f\n  sigmaMisalY %f, sigmaMisalZ %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],fClusterMisalErrorY[i],fClusterMisalErrorZ[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(),
f9119eb9	65	fTracker(0),
	66	fITSonly(kFALSE),
	67	fVertexer(0),
876026b6	68	fClusterFinder(0),
876026b6	69	fPID(0),
7203e11a	70	fVtxr3DZCutWide(0.),
	71	fVtxr3DRCutWide(0.),
	72	fVtxr3DZCutNarrow(0.),
	73	fVtxr3DRCutNarrow(0.),
	74	fVtxr3DPhiCutLoose(0.),
	75	fVtxr3DPhiCutTight(0.),
	76	fVtxr3DDCACut(0.),
e50912db	77	fMaxSnp(1.),
572f41f9	78	fNSigmaYLayerForRoadY(0),
	79	fNSigmaRoadY(0),
	80	fNSigmaZLayerForRoadZ(0),
	81	fNSigmaRoadZ(0),
	82	fNSigma2RoadZC(0),
	83	fNSigma2RoadYC(0),
	84	fNSigma2RoadZNonC(0),
	85	fNSigma2RoadYNonC(0),
1c97ce2f	86	fRoadMisal(0),
572f41f9	87	fMaxNormChi2NonCForHypothesis(0),
	88	fMaxChi2(0),
	89	fMaxRoad(0),
	90	fMaxChi2In(0),
	91	fChi2PerCluster(0),
	92	fXV(0),
	93	fYV(0),
	94	fZV(0),
	95	fSigmaXV(0),
	96	fSigmaYV(0),
	97	fSigmaZV(0),
	98	fVertexCut(0),
	99	fMaxDZforPrimTrk(0),
	100	fMaxDZToUseConstraint(0),
	101	fMaxDforV0dghtrForProlongation(0),
	102	fMaxDForProlongation(0),
	103	fMaxDZForProlongation(0),
	104	fMinPtForProlongation(0),
6518a6c5	105	fAddVirtualClustersInDeadZone(kFALSE),
572f41f9	106	fZWindowDeadZone(0),
	107	fSigmaXDeadZoneHit2(0),
	108	fSigmaZDeadZoneHit2(0),
	109	fXPassDeadZoneHits(0),
6518a6c5	110	fUseTGeoInTracker(3),
572f41f9	111	fAllowSharedClusters(kTRUE),
e50912db	112	fClusterErrorsParam(1),
4a66240a	113	fComputePlaneEff(kFALSE),
5fbd4fd6	114	fHistoPlaneEff(kFALSE),
0ed58a47	115	fIPlanePlaneEff(0),
275a301c	116	fReadPlaneEffFromOCDB(kFALSE),
0ed58a47	117	fMinPtPlaneEff(0),
	118	fMaxMissingClustersPlaneEff(0),
	119	fRequireClusterInOuterLayerPlaneEff(kFALSE),
	120	fRequireClusterInInnerLayerPlaneEff(kFALSE),
	121	fOnlyConstraintPlaneEff(kFALSE),
2755f080	122	fExtendedEtaAcceptance(kFALSE),
b406db2e	123	fUseBadZonesFromOCDB(kTRUE),
23197852	124	fUseSingleBadChannelsFromOCDB(kFALSE),
23197852	125	fMinFractionOfBadInRoad(0),
ae00569a	126	fAllowProlongationWithEmptyRoad(kFALSE),
9f9cae94	127	fOuterStartLayerSA(0),
2755f080	128	fFactorSAWindowSizes(1.),
c7d6d7b7	129	fNLoopsSA(33),
	130	fMinPhiSA(0.002),
	131	fMaxPhiSA(0.0145),
	132	fMinLambdaSA(0.003),
	133	fMaxLambdaSA(0.008),
5a03f353	134	fSAOnePointTracks(kFALSE),
5a03f353	135	fSAUseAllClusters(kFALSE),
6518a6c5	136	fFindV0s(kTRUE),
	137	fUseUnfoldingInClusterFinderSPD(kFALSE),
	138	fUseUnfoldingInClusterFinderSDD(kTRUE),
a86176e3	139	fUseUnfoldingInClusterFinderSSD(kTRUE),
7b116aa1	140	fUseChargeMatchingInClusterFinderSSD(kTRUE),
	141	fTrackleterOnlyOneTrackletPerC2(kTRUE),
	142	fTrackleterPhiWindow(0.08),
	143	fTrackleterZetaWindow(1.00),
	144	fTrackleterRemoveClustersFromOverlaps(kFALSE),
	145	fTrackleterPhiOverlapCut(0.005),
42ed6062	146	fTrackleterZetaOverlapCut(0.05),
42ed6062	147	fUseCosmicRunShiftsSSD(kFALSE)
44347160	148	{
	149	//
	150	// constructor
	151	//
6518a6c5	152	SetName("ITS");
	153	SetTitle("ITS");
	154
44347160	155	SetLayersParameters();
2755f080	156	for(Int_t i=0; i<AliITSgeomTGeo::kNLayers; i++) fLayersToSkip[i]=0;
6518a6c5	157	SetUseTGeoInTracker(3);
afd25725	158	SetAllowSharedClusters(kTRUE);
afd25725	159	SetFindV0s(kTRUE);
6518a6c5	160	SetAddVirtualClustersInDeadZone(kFALSE);
572f41f9	161	SetUseAmplitudeInfo(kTRUE);
e50912db	162	SetClusterErrorsParam(1);
f9119eb9	163	SetClusterMisalError(0.);
7203e11a	164	SetVertexer3DDefaults();
44347160	165	}
	166	//_____________________________________________________________________________
	167	AliITSRecoParam::~AliITSRecoParam()
	168	{
	169	//
	170	// destructor
	171	//
	172	}
	173	//_____________________________________________________________________________
	174	AliITSRecoParam *AliITSRecoParam::GetHighFluxParam()
	175	{
	176	//
	177	// make default reconstruction parameters for hig flux env.
	178	//
	179	AliITSRecoParam *param = new AliITSRecoParam();
14bcd755	180	param->SetVertexerZ();
44347160	181
1aa655ec	182	// use of bads from OCDB
	183	param->SetUseBadZonesFromOCDB(kTRUE);
	184	param->SetUseSingleBadChannelsFromOCDB(kFALSE);
	185	// extended eta acceptance
	186	param->SetExtendedEtaAcceptance(kFALSE);
	187	// allow to skip layer if no cluster and no bad
	188	param->SetAllowProlongationWithEmptyRoad(kFALSE);
	189
44347160	190	param->fMaxSnp = 0.95;
	191
	192	param->fNSigmaYLayerForRoadY = 4.;
	193	param->fNSigmaRoadY = 7.5;
	194	param->fNSigmaZLayerForRoadZ = 4.;
	195	param->fNSigmaRoadZ = 7.5;
	196
	197	param->fNSigma2RoadZC = 60.; //7.75^2
	198	param->fNSigma2RoadYC = 60.; //7.75^2
	199	param->fNSigma2RoadZNonC = 50.; //7.07^2
	200	param->fNSigma2RoadYNonC = 50.; //7.07^2
	201
	202	param->fMaxChi2PerCluster[0] = 11.; //7
	203	param->fMaxChi2PerCluster[1] = 12.; //5
	204	param->fMaxChi2PerCluster[2] = 12.; //8
	205	param->fMaxChi2PerCluster[3] = 5.; //8
	206	param->fMaxChi2PerCluster[4] = 12.; //6.5
	207
	208	param->fMaxNormChi2NonC[0] = 7.;
	209	param->fMaxNormChi2NonC[1] = 8.;
	210	param->fMaxNormChi2NonC[2] = 8.;
	211	param->fMaxNormChi2NonC[3] = 11.;
	212	param->fMaxNormChi2NonC[4] = 14.;
	213	param->fMaxNormChi2NonC[5] = 25.;
	214
	215	param->fMaxNormChi2C[0] = 11.;
	216	param->fMaxNormChi2C[1] = 13.;
	217	param->fMaxNormChi2C[2] = 15.;
	218	param->fMaxNormChi2C[3] = 18.;
	219	param->fMaxNormChi2C[4] = 30.;
	220	param->fMaxNormChi2C[5] = 35.;
afd25725	221
afd25725	222	param->fMaxNormChi2NonCForHypothesis = 7.;
44347160	223
	224	param->fMaxChi2 = 35.;
	225
	226	param->fMaxChi2s[0] = 25.; //40
	227	param->fMaxChi2s[1] = 25.; //40
	228	param->fMaxChi2s[2] = 25.; //40
	229	param->fMaxChi2s[3] = 25.; //40
	230	param->fMaxChi2s[4] = 40.; //40
	231	param->fMaxChi2s[5] = 50.; //40
	232
	233	param->fMaxRoad = 6.;
	234
	235	// not used
	236	param->fMaxChi2In = 16.;
	237
	238	param->fMaxChi2sR[0] = 10.;
	239	param->fMaxChi2sR[1] = 10.;
	240	param->fMaxChi2sR[2] = 10.;
	241	param->fMaxChi2sR[3] = 10.;
	242	param->fMaxChi2sR[4] = 30.;
	243	param->fMaxChi2sR[5] = 40.;
	244
	245	param->fChi2PerCluster = 9.;
	246	// not used
	247
	248	param->fXV = 0.;
	249	param->fYV = 0.;
	250	param->fZV = 0.;
	251	param->fSigmaXV = 0.0050;
	252	param->fSigmaYV = 0.0050;
	253	param->fSigmaZV = 0.0100;
	254
	255	param->fVertexCut = 25.;
afd25725	256
	257	param->fMaxDZforPrimTrk = 0.4;
	258	param->fMaxDZToUseConstraint = 3.;
	259
	260	param->fMaxDforV0dghtrForProlongation = 30.;
6518a6c5	261	param->fMaxDForProlongation = 40.;//10.;
6518a6c5	262	param->fMaxDZForProlongation = 60.;//20.;
afd25725	263	param->fMinPtForProlongation = 0.120;
	264
	265	param->fZWindowDeadZone = 2.0;
	266	param->fSigmaXDeadZoneHit2 = 0.004/12.;
	267	param->fSigmaZDeadZoneHit2 = 0.001/12.;
	268	param->fXPassDeadZoneHits = 0.018;
44347160	269
	270	return param;
	271	}
	272	//_____________________________________________________________________________
	273	AliITSRecoParam *AliITSRecoParam::GetLowFluxParam()
	274	{
	275	//
	276	// make default reconstruction parameters for low flux env.
	277	//
c7d6d7b7	278	AliITSRecoParam *param = new AliITSRecoParam();
7203e11a	279	param->SetVertexer3DDefaults();
c7d6d7b7	280
1aa655ec	281	// full use of bads from OCDB
	282	param->SetUseBadZonesFromOCDB(kTRUE);
	283	param->SetUseSingleBadChannelsFromOCDB(kTRUE);
	284	// extended eta acceptance
	285	param->SetExtendedEtaAcceptance(kTRUE);
	286	// allow to skip layer if no cluster and no bad
	287	param->SetAllowProlongationWithEmptyRoad(kTRUE);
	288
c7d6d7b7	289	param->fMaxSnp = 0.95;
	290
	291	param->fNSigmaYLayerForRoadY = 4.;
	292	param->fNSigmaRoadY = 7.5;
	293	param->fNSigmaZLayerForRoadZ = 4.;
	294	param->fNSigmaRoadZ = 7.5;
	295
	296	param->fNSigma2RoadZC = 60.; //7.75^2
	297	param->fNSigma2RoadYC = 60.; //7.75^2
	298	param->fNSigma2RoadZNonC = 50.; //7.07^2
	299	param->fNSigma2RoadYNonC = 50.; //7.07^2
	300
	301	param->fMaxChi2PerCluster[0] = 11.; //7
	302	param->fMaxChi2PerCluster[1] = 12.; //5
	303	param->fMaxChi2PerCluster[2] = 12.; //8
	304	param->fMaxChi2PerCluster[3] = 5.; //8
	305	param->fMaxChi2PerCluster[4] = 12.; //6.5
	306
	307	param->fMaxNormChi2NonC[0] = 7.;
	308	param->fMaxNormChi2NonC[1] = 8.;
	309	param->fMaxNormChi2NonC[2] = 8.;
	310	param->fMaxNormChi2NonC[3] = 11.;
	311	param->fMaxNormChi2NonC[4] = 14.;
	312	param->fMaxNormChi2NonC[5] = 25.;
	313
	314	param->fMaxNormChi2C[0] = 11.;
	315	param->fMaxNormChi2C[1] = 13.;
	316	param->fMaxNormChi2C[2] = 15.;
	317	param->fMaxNormChi2C[3] = 18.;
	318	param->fMaxNormChi2C[4] = 30.;
	319	param->fMaxNormChi2C[5] = 35.;
	320
	321	param->fMaxNormChi2NonCForHypothesis = 7.;
	322
	323	param->fMaxChi2 = 35.;
	324
	325	param->fMaxChi2s[0] = 25.; //40
	326	param->fMaxChi2s[1] = 25.; //40
	327	param->fMaxChi2s[2] = 25.; //40
	328	param->fMaxChi2s[3] = 25.; //40
	329	param->fMaxChi2s[4] = 40.; //40
	330	param->fMaxChi2s[5] = 50.; //40
	331
	332	param->fMaxRoad = 6.;
	333
	334	// not used
	335	param->fMaxChi2In = 16.;
	336
	337	param->fMaxChi2sR[0] = 10.;
	338	param->fMaxChi2sR[1] = 10.;
	339	param->fMaxChi2sR[2] = 10.;
	340	param->fMaxChi2sR[3] = 10.;
	341	param->fMaxChi2sR[4] = 30.;
	342	param->fMaxChi2sR[5] = 40.;
	343
	344	param->fChi2PerCluster = 9.;
	345	// not used
	346
	347	param->fXV = 0.;
	348	param->fYV = 0.;
	349	param->fZV = 0.;
	350	param->fSigmaXV = 0.0050;
	351	param->fSigmaYV = 0.0050;
	352	param->fSigmaZV = 0.0100;
353
354	param->fVertexCut = 25.;
355
356	param->fMaxDZforPrimTrk = 0.4;
357	param->fMaxDZToUseConstraint = 3.;
358
359	param->fMaxDforV0dghtrForProlongation = 30.;
360	param->fMaxDForProlongation = 40.;//10.;
361	param->fMaxDZForProlongation = 60.;//20.;
362	param->fMinPtForProlongation = 0.120;
363
364	param->fZWindowDeadZone = 2.0;
365	param->fSigmaXDeadZoneHit2 = 0.004/12.;
366	param->fSigmaZDeadZoneHit2 = 0.001/12.;
367	param->fXPassDeadZoneHits = 0.018;
368	param->SetNLoopsSA(10);
7203e11a	369
c7d6d7b7	370	return param;
44347160	371	}
	372	//_____________________________________________________________________________
	373	AliITSRecoParam *AliITSRecoParam::GetCosmicTestParam()
	374	{
	375	//
	376	// make default reconstruction parameters for cosmics
	377	//
e340bb86	378	AliITSRecoParam *param = new AliITSRecoParam();
e340bb86	379
f9119eb9	380	// vertexer for cosmics
	381	param->SetVertexer(2);
	382
9f9cae94	383	// find independently ITS SA tracks
	384	param->SetSAUseAllClusters();
	385	param->SetOuterStartLayerSA(AliITSgeomTGeo::GetNLayers()-2);
f9119eb9	386
2755f080	387	// larger seach windows for SA (in case of large misalignments)
2755f080	388	param->SetFactorSAWindowSizes(3.);
e340bb86	389
	390	param->fMaxSnp = 0.95;
	391
	392	param->fNSigmaYLayerForRoadY = 4.;
	393	param->fNSigmaRoadY = 7.5;
	394	param->fNSigmaZLayerForRoadZ = 4.;
	395	param->fNSigmaRoadZ = 7.5;
	396
	397	param->fNSigma2RoadZC = 60.; //7.75^2
	398	param->fNSigma2RoadYC = 60.; //7.75^2
	399	param->fNSigma2RoadZNonC = 50.; //7.07^2
	400	param->fNSigma2RoadYNonC = 50.; //7.07^2
	401
	402	param->fMaxChi2PerCluster[0] = 11.; //7
	403	param->fMaxChi2PerCluster[1] = 12.; //5
	404	param->fMaxChi2PerCluster[2] = 12.; //8
	405	param->fMaxChi2PerCluster[3] = 5.; //8
	406	param->fMaxChi2PerCluster[4] = 12.; //6.5
	407
	408	param->fMaxNormChi2NonC[0] = 7.;
	409	param->fMaxNormChi2NonC[1] = 8.;
	410	param->fMaxNormChi2NonC[2] = 8.;
	411	param->fMaxNormChi2NonC[3] = 11.;
	412	param->fMaxNormChi2NonC[4] = 14.;
	413	param->fMaxNormChi2NonC[5] = 25.;
	414
	415	param->fMaxNormChi2C[0] = 11.;
	416	param->fMaxNormChi2C[1] = 13.;
	417	param->fMaxNormChi2C[2] = 15.;
	418	param->fMaxNormChi2C[3] = 18.;
	419	param->fMaxNormChi2C[4] = 30.;
	420	param->fMaxNormChi2C[5] = 35.;
	421
	422	param->fMaxNormChi2NonCForHypothesis = 7.;
	423
	424	param->fMaxChi2 = 35.;
	425
	426	param->fMaxChi2s[0] = 25.; //40
	427	param->fMaxChi2s[1] = 25.; //40
	428	param->fMaxChi2s[2] = 25.; //40
	429	param->fMaxChi2s[3] = 25.; //40
	430	param->fMaxChi2s[4] = 40.; //40
	431	param->fMaxChi2s[5] = 50.; //40
	432
	433	param->fMaxRoad = 6.;
	434
	435	// not used
	436	param->fMaxChi2In = 16.;
	437
	438	param->fMaxChi2sR[0] = 10.;
	439	param->fMaxChi2sR[1] = 10.;
	440	param->fMaxChi2sR[2] = 10.;
	441	param->fMaxChi2sR[3] = 10.;
	442	param->fMaxChi2sR[4] = 30.;
	443	param->fMaxChi2sR[5] = 40.;
	444
	445	param->fChi2PerCluster = 9.;
	446	// not used
	447
	448	param->fXV = 0.;
	449	param->fYV = 0.;
	450	param->fZV = 0.;
	451	param->fSigmaXV = 0.0050;
	452	param->fSigmaYV = 0.0050;
453	param->fSigmaZV = 0.0100;
454
455	param->fVertexCut = 25.;
456
457	param->fMaxDZforPrimTrk = 0.4;
458	param->fMaxDZToUseConstraint = 3.;
459
460	param->fMaxDforV0dghtrForProlongation = 30.;
461	param->fMaxDForProlongation = 10.;
462	param->fMaxDZForProlongation = 20.;
463	param->fMinPtForProlongation = 0.120;
464
465	param->fZWindowDeadZone = 2.0;
466	param->fSigmaXDeadZoneHit2 = 0.004/12.;
467	param->fSigmaZDeadZoneHit2 = 0.001/12.;
468	param->fXPassDeadZoneHits = 0.018;
469
470
471	return param;
44347160	472	}
44347160	473	//_____________________________________________________________________________
ae00569a	474	AliITSRecoParam *AliITSRecoParam::GetPlaneEffParam(Int_t i)
	475	{
	476	//
	477	// make special reconstruction parameters for Plane Efficiency study on layer i
	478	//
0ed58a47	479	if (i<0 \|\| i>=AliITSgeomTGeo::kNLayers) {
	480	printf("AliITSRecoParam::GetPlaneEffParam: index of ITS Plane not in the range [0,5]\n");
	481	printf("returning null pointer");
	482	return NULL;
	483	}
ae00569a	484	AliITSRecoParam *param;
	485	param = GetHighFluxParam();
	486	param->SetComputePlaneEff();
	487	param->SetLayerToSkip(i);
0ed58a47	488	param->SetIPlanePlaneEff(i);
	489	// optimized setting for SPD0 (i==0)
	490	if (i==0 \|\| i==1) {
	491	param->fMinPtPlaneEff = 0.200; // high pt particles
	492	param->fMaxMissingClustersPlaneEff = 1; // at most 1 layer out of 5 without cluster
	493	param->fRequireClusterInOuterLayerPlaneEff = kTRUE; // cluster on SPD1 must be
	494	//param->fOnlyConstraintPlaneEff = kTRUE;
	495	}
	496	if (i==2 \|\| i==3) {
	497	param->fMinPtPlaneEff = 0.200; // high pt particles
	498	param->fMaxMissingClustersPlaneEff = 1; // at most 1 layer out of 5 without cluster
	499	param->fRequireClusterInOuterLayerPlaneEff = kTRUE;
	500	//param->fOnlyConstraintPlaneEff = kTRUE;
	501	}
	502	if (i==4) {
	503	param->fMinPtPlaneEff = 0.200; // high pt particles
	504	param->fMaxMissingClustersPlaneEff = 0; // at most 1 layer out of 5 without cluster
	505	param->fRequireClusterInOuterLayerPlaneEff = kTRUE;
	506	//param->fOnlyConstraintPlaneEff = kTRUE;
	507	}
	508	if (i==5) {
	509	param->fMinPtPlaneEff = 0.200; // high pt particles
	510	}
	511	//
ae00569a	512	return param;
	513	}
	514	//_____________________________________________________________________________
44347160	515	void AliITSRecoParam::SetLayersParameters()
	516	{
	517	//
	518	// number of layers and layers spatial resolutions
	519	//
	520
	521	// spatial resolutions of the detectors
	522	// y: 12 12 38 38 20 20 micron
	523	fSigmaY2[0]=1.44e-6;
	524	fSigmaY2[1]=1.44e-6;
	525	fSigmaY2[2]=1.444e-5;
	526	fSigmaY2[3]=1.444e-5;
	527	fSigmaY2[4]=4.0e-6;
	528	fSigmaY2[5]=4.0e-6;
	529	// z: 120 120 28 28 830 830 micron
	530	fSigmaZ2[0]=1.44e-4;
	531	fSigmaZ2[1]=1.44e-4;
	532	fSigmaZ2[2]=7.84e-6;
	533	fSigmaZ2[3]=7.84e-6;
	534	fSigmaZ2[4]=6.889e-3;
	535	fSigmaZ2[5]=6.889e-3;
	536
	537	return;
	538	}
ed446fa3	539	//_____________________________________________________________________________
	540	void AliITSRecoParam::PrintParameters() const
	541	{
	542	//
	543	// print parameters
	544	//
	545
	546	printf("============================= AliITSRecoParam::PrintParameters ");
f9119eb9	547	printf("============================= \n\n");
ed446fa3	548	for(Int_t i=0; i<AliITSgeomTGeo::kNLayers; i++) {
	549	if(!fLayersToSkip[i]) {
	550	printf("ITS Traking: using layer %d\n",i);
	551	} else {
	552	printf("ITS Traking: skipping layer %d\n",i);
	553	}
	554	}
	555	for(Int_t i=0; i<AliITSgeomTGeo::kNLayers; i++) {
	556	if(fUseAmplitudeInfo[i]) {
	557	printf("ITS Traking: use amplitude info for layer %d\n",i);
	558	} else {
	559	printf("ITS Traking: don't use amplitude info for layer %d\n",i);
	560	}
	561	}
	562	for(Int_t i=0; i<AliITSgeomTGeo::kNLayers; i++)
f9119eb9	563	printf("Layer %d:\n sigmaY2 %f, sigmaZ2 %f\n sigmaMisalY %f, sigmaMisalZ %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],fClusterMisalErrorY[i],fClusterMisalErrorZ[i],fMaxNormChi2NonC[i],fMaxNormChi2C[i],fMaxChi2s[i]);
ed446fa3	564
	565
	566	Dump();
	567
	568	return;
	569	}