[u/mrichter/AliRoot.git] / ITS / AliITSRecoParam.cxx

/**************************************************************************
 * Copyright(c) 2007-2009, ALICE Experiment at CERN, All rights reserved. *
 *                                                                        *
 * Author: The ALICE Off-line Project.                                    *
 * Contributors are mentioned in the code where appropriate.              *
 *                                                                        *
 * Permission to use, copy, modify and distribute this software and its   *
 * documentation strictly for non-commercial purposes is hereby granted   *
 * without fee, provided that the above copyright notice appears in all   *
 * copies and that both the copyright notice and this permission notice   *
 * appear in the supporting documentation. The authors make no claims     *
 * about the suitability of this software for any purpose. It is          *
 * provided "as is" without express or implied warranty.                  *
 **************************************************************************/

#include "AliITSRecoParam.h"

/* $Id$ */

///////////////////////////////////////////////////////////////////////////////
//                                                                           //
// Class with ITS reconstruction parameters                                  //
// Origin: andrea.dainese@lnl.infn.it                                        //
//                                                                           //
///////////////////////////////////////////////////////////////////////////////


ClassImp(AliITSRecoParam)

const Int_t AliITSRecoParam::fgkLayersNotToSkip[AliITSgeomTGeo::kNLayers]={0,0,0,0,0,0};
const Int_t AliITSRecoParam::fgkLastLayerToTrackTo=0;
const Int_t AliITSRecoParam::fgkMaxDetectorPerLayer=1000;
const Double_t AliITSRecoParam::fgkriw=80.0;
const Double_t AliITSRecoParam::fgkdiw=0.0053;
const Double_t AliITSRecoParam::fgkX0iw=30.0;
const Double_t AliITSRecoParam::fgkrcd=61.0;
const Double_t AliITSRecoParam::fgkdcd=0.0053;
const Double_t AliITSRecoParam::fgkX0cd=30.0;
const Double_t AliITSRecoParam::fgkyr=12.8;
const Double_t AliITSRecoParam::fgkdr=0.03;
const Double_t AliITSRecoParam::fgkzm=0.2;
const Double_t AliITSRecoParam::fgkdm=0.40;
const Double_t AliITSRecoParam::fgkrs=50.0;
const Double_t AliITSRecoParam::fgkds=0.001;
const Double_t AliITSRecoParam::fgkrInsideITSscreen=49.0;
const Double_t AliITSRecoParam::fgkrInsideSPD1=3.7;
const Double_t AliITSRecoParam::fgkrPipe=3.;
const Double_t AliITSRecoParam::fgkrInsidePipe=2.7;
const Double_t AliITSRecoParam::fgkrOutsidePipe=3.3;
const Double_t AliITSRecoParam::fgkdPipe=0.0028;
const Double_t AliITSRecoParam::fgkrInsideShield[2]={7.5,25.0};
const Double_t AliITSRecoParam::fgkrOutsideShield[2]={10.5,30.0};
const Double_t AliITSRecoParam::fgkdshield[2]={0.0097,0.0034};
const Double_t AliITSRecoParam::fgkX0shield[2]={38.6,42.0};
const Double_t AliITSRecoParam::fgkX0Air=21.82;
const Double_t AliITSRecoParam::fgkX0Be=65.19;
const Double_t AliITSRecoParam::fgkBoundaryWidth=0.2;
const Double_t AliITSRecoParam::fgkDeltaXNeighbDets=0.5;
const Double_t AliITSRecoParam::fgkSPDdetzlength=6.960; // 7.072-2*0.056
const Double_t AliITSRecoParam::fgkSPDdetxlength=1.298; // 1.410-2*0.056

//_____________________________________________________________________________
AliITSRecoParam::AliITSRecoParam() : AliDetectorRecoParam(),
fTracker(0),
fITSonly(kFALSE),
fVertexer(0),
fClusterFinder(0),
fPID(0),
fVtxr3DZCutWide(0.),
fVtxr3DRCutWide(0.),
fVtxr3DZCutNarrow(0.),
fVtxr3DRCutNarrow(0.),
fVtxr3DPhiCutLoose(0.),
fVtxr3DPhiCutTight(0.),
fVtxr3DDCACut(0.),
fMaxSnp(1.),
fNSigmaYLayerForRoadY(0),
fNSigmaRoadY(0),
fNSigmaZLayerForRoadZ(0),
fNSigmaRoadZ(0),
fNSigma2RoadZC(0),
fNSigma2RoadYC(0),
fNSigma2RoadZNonC(0),
fNSigma2RoadYNonC(0),
fRoadMisal(0),
fMaxNormChi2NonCForHypothesis(0),
fMaxChi2(0),
fMaxRoad(0),
fMaxChi2In(0),
fChi2PerCluster(0),
fXV(0), 
fYV(0),
fZV(0),
fSigmaXV(0),
fSigmaYV(0),
fSigmaZV(0),
fVertexCut(0),
fMaxDZforPrimTrk(0),
fMaxDZToUseConstraint(0), 
fMaxDforV0dghtrForProlongation(0),
fMaxDForProlongation(0),
fMaxDZForProlongation(0),
fMinPtForProlongation(0),
fAddVirtualClustersInDeadZone(kFALSE),
fZWindowDeadZone(0),
fSigmaXDeadZoneHit2(0),
fSigmaZDeadZoneHit2(0),
fXPassDeadZoneHits(0),
fUseTGeoInTracker(3),
fAllowSharedClusters(kTRUE),
fClusterErrorsParam(1),
fComputePlaneEff(kFALSE),
fHistoPlaneEff(kFALSE),
fIPlanePlaneEff(0),
fReadPlaneEffFromOCDB(kFALSE),
fMinPtPlaneEff(0),
fMaxMissingClustersPlaneEff(0),
fRequireClusterInOuterLayerPlaneEff(kFALSE),
fRequireClusterInInnerLayerPlaneEff(kFALSE),
fOnlyConstraintPlaneEff(kFALSE),
fExtendedEtaAcceptance(kFALSE),
fUseBadZonesFromOCDB(kTRUE),
fUseSingleBadChannelsFromOCDB(kFALSE),
fMinFractionOfBadInRoad(0),
fAllowProlongationWithEmptyRoad(kFALSE),
fOuterStartLayerSA(0),
fFactorSAWindowSizes(1.),
fNLoopsSA(33),
fMinPhiSA(0.002),
fMaxPhiSA(0.0145),
fMinLambdaSA(0.003),
fMaxLambdaSA(0.008),
fSAOnePointTracks(kFALSE),
fSAUseAllClusters(kFALSE),
fFindV0s(kTRUE),
fUseUnfoldingInClusterFinderSPD(kFALSE),
fUseUnfoldingInClusterFinderSDD(kTRUE),
fUseUnfoldingInClusterFinderSSD(kTRUE),
fUseChargeMatchingInClusterFinderSSD(kTRUE),
fTrackleterOnlyOneTrackletPerC2(kTRUE),
fTrackleterPhiWindow(0.08),
fTrackleterZetaWindow(1.00),
fTrackleterRemoveClustersFromOverlaps(kFALSE),
fTrackleterPhiOverlapCut(0.005),
fTrackleterZetaOverlapCut(0.05),
fUseCosmicRunShiftsSSD(kFALSE)
{
  //
  // constructor
  //
  SetName("ITS");
  SetTitle("ITS");

  SetLayersParameters();
  for(Int_t i=0; i<AliITSgeomTGeo::kNLayers; i++) fLayersToSkip[i]=0;
  SetUseTGeoInTracker(3);
  SetAllowSharedClusters(kTRUE);
  SetFindV0s(kTRUE);
  SetAddVirtualClustersInDeadZone(kFALSE);
  SetUseAmplitudeInfo(kTRUE);
  SetClusterErrorsParam(1);
  SetClusterMisalError(0.);
  SetVertexer3DDefaults();
}
//_____________________________________________________________________________
AliITSRecoParam::~AliITSRecoParam() 
{
  //
  // destructor
  //  
}
//_____________________________________________________________________________
AliITSRecoParam *AliITSRecoParam::GetHighFluxParam() 
{
  //
  // make default reconstruction  parameters for hig  flux env.
  //
  AliITSRecoParam *param = new AliITSRecoParam();
  param->SetVertexerZ();

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

  param->fMaxSnp = 0.95;

  param->fNSigmaYLayerForRoadY = 4.;
  param->fNSigmaRoadY = 7.5;
  param->fNSigmaZLayerForRoadZ = 4.;
  param->fNSigmaRoadZ = 7.5;

  param->fNSigma2RoadZC = 60.; //7.75^2
  param->fNSigma2RoadYC = 60.; //7.75^2
  param->fNSigma2RoadZNonC = 50.; //7.07^2
  param->fNSigma2RoadYNonC = 50.; //7.07^2

  param->fMaxChi2PerCluster[0] = 11.; //7
  param->fMaxChi2PerCluster[1] = 12.; //5
  param->fMaxChi2PerCluster[2] = 12.; //8
  param->fMaxChi2PerCluster[3] = 5.;  //8
  param->fMaxChi2PerCluster[4] = 12.; //6.5

  param->fMaxNormChi2NonC[0] = 7.;
  param->fMaxNormChi2NonC[1] = 8.;
  param->fMaxNormChi2NonC[2] = 8.;
  param->fMaxNormChi2NonC[3] = 11.;
  param->fMaxNormChi2NonC[4] = 14.;
  param->fMaxNormChi2NonC[5] = 25.;

  param->fMaxNormChi2C[0] = 11.;
  param->fMaxNormChi2C[1] = 13.;
  param->fMaxNormChi2C[2] = 15.;
  param->fMaxNormChi2C[3] = 18.;
  param->fMaxNormChi2C[4] = 30.;
  param->fMaxNormChi2C[5] = 35.;

  param->fMaxNormChi2NonCForHypothesis = 7.;
  
  param->fMaxChi2 = 35.;

  param->fMaxChi2s[0] = 25.; //40   
  param->fMaxChi2s[1] = 25.; //40   
  param->fMaxChi2s[2] = 25.; //40   
  param->fMaxChi2s[3] = 25.; //40   
  param->fMaxChi2s[4] = 40.; //40   
  param->fMaxChi2s[5] = 50.; //40

  param->fMaxRoad = 6.;

  // not used
  param->fMaxChi2In = 16.;
   
  param->fMaxChi2sR[0] = 10.;   
  param->fMaxChi2sR[1] = 10.;   
  param->fMaxChi2sR[2] = 10.;   
  param->fMaxChi2sR[3] = 10.;   
  param->fMaxChi2sR[4] = 30.;   
  param->fMaxChi2sR[5] = 40.;   

  param->fChi2PerCluster = 9.;
  // not used

  param->fXV = 0.;
  param->fYV = 0.;
  param->fZV = 0.;
  param->fSigmaXV = 0.0050;
  param->fSigmaYV = 0.0050;
  param->fSigmaZV = 0.0100;

  param->fVertexCut = 25.;

  param->fMaxDZforPrimTrk = 0.4;
  param->fMaxDZToUseConstraint = 3.;

  param->fMaxDforV0dghtrForProlongation = 30.;
  param->fMaxDForProlongation = 40.;//10.;
  param->fMaxDZForProlongation = 60.;//20.;
  param->fMinPtForProlongation = 0.120;

  param->fZWindowDeadZone = 2.0;
  param->fSigmaXDeadZoneHit2 = 0.004/12.;
  param->fSigmaZDeadZoneHit2 = 0.001/12.;
  param->fXPassDeadZoneHits = 0.018;
  param->SetNLoopsSA(10);

  return param;
}
//_____________________________________________________________________________
AliITSRecoParam *AliITSRecoParam::GetCosmicTestParam() 
{
  //
  // make default reconstruction  parameters for cosmics
  //
  AliITSRecoParam *param = new AliITSRecoParam();

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

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

  // larger seach windows for SA (in case of large misalignments)
  param->SetFactorSAWindowSizes(3.); 

  param->fMaxSnp = 0.95;

  param->fNSigmaYLayerForRoadY = 4.;
  param->fNSigmaRoadY = 7.5;
  param->fNSigmaZLayerForRoadZ = 4.;
  param->fNSigmaRoadZ = 7.5;

  param->fNSigma2RoadZC = 60.; //7.75^2
  param->fNSigma2RoadYC = 60.; //7.75^2
  param->fNSigma2RoadZNonC = 50.; //7.07^2
  param->fNSigma2RoadYNonC = 50.; //7.07^2

  param->fMaxChi2PerCluster[0] = 11.; //7
  param->fMaxChi2PerCluster[1] = 12.; //5
  param->fMaxChi2PerCluster[2] = 12.; //8
  param->fMaxChi2PerCluster[3] = 5.;  //8
  param->fMaxChi2PerCluster[4] = 12.; //6.5

  param->fMaxNormChi2NonC[0] = 7.;
  param->fMaxNormChi2NonC[1] = 8.;
  param->fMaxNormChi2NonC[2] = 8.;
  param->fMaxNormChi2NonC[3] = 11.;
  param->fMaxNormChi2NonC[4] = 14.;
  param->fMaxNormChi2NonC[5] = 25.;

  param->fMaxNormChi2C[0] = 11.;
  param->fMaxNormChi2C[1] = 13.;
  param->fMaxNormChi2C[2] = 15.;
  param->fMaxNormChi2C[3] = 18.;
  param->fMaxNormChi2C[4] = 30.;
  param->fMaxNormChi2C[5] = 35.;

  param->fMaxNormChi2NonCForHypothesis = 7.;
  
  param->fMaxChi2 = 35.;

  param->fMaxChi2s[0] = 25.; //40   
  param->fMaxChi2s[1] = 25.; //40   
  param->fMaxChi2s[2] = 25.; //40   
  param->fMaxChi2s[3] = 25.; //40   
  param->fMaxChi2s[4] = 40.; //40   
  param->fMaxChi2s[5] = 50.; //40

  param->fMaxRoad = 6.;

  // not used
  param->fMaxChi2In = 16.;
   
  param->fMaxChi2sR[0] = 10.;   
  param->fMaxChi2sR[1] = 10.;   
  param->fMaxChi2sR[2] = 10.;   
  param->fMaxChi2sR[3] = 10.;   
  param->fMaxChi2sR[4] = 30.;   
  param->fMaxChi2sR[5] = 40.;   

  param->fChi2PerCluster = 9.;
  // not used

  param->fXV = 0.;
  param->fYV = 0.;
  param->fZV = 0.;
  param->fSigmaXV = 0.0050;
  param->fSigmaYV = 0.0050;
  param->fSigmaZV = 0.0100;

  param->fVertexCut = 25.;

  param->fMaxDZforPrimTrk = 0.4;
  param->fMaxDZToUseConstraint = 3.;

  param->fMaxDforV0dghtrForProlongation = 30.;
  param->fMaxDForProlongation = 10.;
  param->fMaxDZForProlongation = 20.;
  param->fMinPtForProlongation = 0.120;

  param->fZWindowDeadZone = 2.0;
  param->fSigmaXDeadZoneHit2 = 0.004/12.;
  param->fSigmaZDeadZoneHit2 = 0.001/12.;
  param->fXPassDeadZoneHits = 0.018;

  
  return param;
}
//_____________________________________________________________________________
AliITSRecoParam *AliITSRecoParam::GetPlaneEffParam(Int_t i)
{
  //
  // make special reconstruction parameters for Plane Efficiency study on layer i
  //
  if (i<0 || i>=AliITSgeomTGeo::kNLayers) {
    printf("AliITSRecoParam::GetPlaneEffParam: index of ITS Plane not in the range [0,5]\n");
    printf("returning null pointer");
    return NULL;
  }
  AliITSRecoParam *param;
  param = GetHighFluxParam();
  param->SetComputePlaneEff();
  param->SetLayerToSkip(i);
  param->SetIPlanePlaneEff(i);
  // optimized setting for SPD0 (i==0)
  if (i==0 || i==1) {
    param->fMinPtPlaneEff = 0.200; // high pt particles
    param->fMaxMissingClustersPlaneEff = 1; // at most 1 layer out of 5 without cluster
    param->fRequireClusterInOuterLayerPlaneEff = kTRUE; // cluster on SPD1 must be
    //param->fOnlyConstraintPlaneEff = kTRUE;
  }
  if (i==2 || i==3) {
    param->fMinPtPlaneEff = 0.200; // high pt particles
    param->fMaxMissingClustersPlaneEff = 1; // at most 1 layer out of 5 without cluster
    param->fRequireClusterInOuterLayerPlaneEff = kTRUE;
    //param->fOnlyConstraintPlaneEff = kTRUE;
  }
  if (i==4) {
    param->fMinPtPlaneEff = 0.200; // high pt particles
    param->fMaxMissingClustersPlaneEff = 0; // at most 1 layer out of 5 without cluster
    param->fRequireClusterInOuterLayerPlaneEff = kTRUE;
    //param->fOnlyConstraintPlaneEff = kTRUE;
  }
  if (i==5) {
    param->fMinPtPlaneEff = 0.200; // high pt particles
  }
  //
  return param;
}
//_____________________________________________________________________________
void AliITSRecoParam::SetLayersParameters() 
{
  //
  // number of layers and layers spatial resolutions
  //

  // spatial resolutions of the detectors
  // y: 12 12 38 38 20 20 micron
  fSigmaY2[0]=1.44e-6;
  fSigmaY2[1]=1.44e-6;
  fSigmaY2[2]=1.444e-5;
  fSigmaY2[3]=1.444e-5;
  fSigmaY2[4]=4.0e-6;
  fSigmaY2[5]=4.0e-6;
  // z: 120 120 28 28 830 830 micron
  fSigmaZ2[0]=1.44e-4;
  fSigmaZ2[1]=1.44e-4;
  fSigmaZ2[2]=7.84e-6;
  fSigmaZ2[3]=7.84e-6;
  fSigmaZ2[4]=6.889e-3;
  fSigmaZ2[5]=6.889e-3;

  return;
}
//_____________________________________________________________________________
void AliITSRecoParam::PrintParameters() const 
{
  //
  // print parameters
  //

  printf("=============================  AliITSRecoParam::PrintParameters ");
  printf("============================= \n\n");
  for(Int_t i=0; i<AliITSgeomTGeo::kNLayers; i++) {
    if(!fLayersToSkip[i]) {
      printf("ITS Traking: using layer %d\n",i);
    } else {
      printf("ITS Traking: skipping layer %d\n",i);
    }
  }
  for(Int_t i=0; i<AliITSgeomTGeo::kNLayers; i++) {
    if(fUseAmplitudeInfo[i]) {
      printf("ITS Traking: use amplitude info for layer %d\n",i);
    } else {
      printf("ITS Traking: don't use amplitude info for layer %d\n",i);
    }
  }
  for(Int_t i=0; i<AliITSgeomTGeo::kNLayers; i++)
    printf("Layer %d:\n  sigmaY2 %f, sigmaZ2 %f\n  sigmaMisalY %f, sigmaMisalZ %f\n  max norm chi2 for non constrained tracks %f\n  max norm chi2 for constrained tracks %f\n  max predicted chi2 (cluster & track prol.) %f\n",i,fSigmaY2[i],fSigmaZ2[i],fClusterMisalErrorY[i],fClusterMisalErrorZ[i],fMaxNormChi2NonC[i],fMaxNormChi2C[i],fMaxChi2s[i]);


  Dump();

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