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