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