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