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