[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"
#include "AliLog.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(5.),
fAllowSharedClusters(kTRUE),
fClusterErrorsParam(1),
fComputePlaneEff(kFALSE),
fHistoPlaneEff(kFALSE),
fUseTrackletsPlaneEff(kFALSE),
fMCTrackletsPlaneEff(kFALSE),
fBkgTrackletsPlaneEff(kFALSE),
fTrackleterPhiWindowL1(0.10),
fTrackleterZetaWindowL1(0.6),
fUpdateOncePerEventPlaneEff(kTRUE),
fMinContVtxPlaneEff(3),
fIPlanePlaneEff(0),
fReadPlaneEffFromOCDB(kFALSE),
fMinPtPlaneEff(0),
fMaxMissingClustersPlaneEff(0),
fRequireClusterInOuterLayerPlaneEff(kFALSE),
fRequireClusterInInnerLayerPlaneEff(kFALSE),
fOnlyConstraintPlaneEff(kFALSE),
fExtendedEtaAcceptance(kFALSE),
fUseBadZonesFromOCDB(kTRUE),
fUseSingleBadChannelsFromOCDB(kFALSE),
fMinFractionOfBadInRoad(0),
fAllowProlongationWithEmptyRoad(kFALSE),
fInwardFlagSA(kFALSE),
fOuterStartLayerSA(2),
fInnerStartLayerSA(3),
fMinNPointsSA(3),
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),
fSPDRemoveNoisyFlag(kTRUE),
fSPDRemoveDeadFlag(kTRUE)
{
  //
  // constructor
  //
  SetName("ITS");
  SetTitle("ITS");

  SetLayersParameters();
  for(Int_t i=0; i<AliITSgeomTGeo::kNLayers; i++) fLayersToSkip[i]=0;
  SetUseTGeoInTracker(3);
  SetStepSizeTGeo(5.);
  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 (0,5)
  // 
  // if i=-1, then the evaluation for both pixel layers is tried with the tracklet method
  //
 if (i<-1 || i>=AliITSgeomTGeo::kNLayers) {
    printf("AliITSRecoParam::GetPlaneEffParam: index of ITS Plane nor in the range [0,5] neither =-1\n");
    printf("returning null pointer");
    return NULL;
  }
  if(i>=0) {  // Method using tracks (remove given plane from tracking)
    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;
  }
  else if (i==-1) { // Method using tracklets
    AliITSRecoParam *param;
    param = GetLowFluxParam();
    param->SetIPlanePlaneEff(i);
    param->SetComputePlaneEff(kTRUE,kFALSE);
    param->SetUseTrackletsPlaneEff(kTRUE);
    param->SetTrackleterPhiWindow(0.07);
    param->SetTrackleterZetaWindow(0.4);
    param->SetTrackleterPhiWindowL1(0.10);
    param->SetTrackleterZetaWindowL1(0.6);
    param->SetUpdateOncePerEventPlaneEff(kTRUE);
    param->SetMinContVtxPlaneEff(3);
    return param;
  }
  else {
    AliErrorGeneral("AliITSRecoParam",Form("Unrecognised value of i %d\n",i));
    return 0;
  }
}
//_____________________________________________________________________________
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"
5bd7ec3a	17	#include "AliLog.h"
572f41f9	18
572f41f9	19	/* $Id$ */
44347160	20
	21	///////////////////////////////////////////////////////////////////////////////
	22	// //
	23	// Class with ITS reconstruction parameters //
	24	// Origin: andrea.dainese@lnl.infn.it //
	25	// //
	26	///////////////////////////////////////////////////////////////////////////////
	27
	28
44347160	29
	30	ClassImp(AliITSRecoParam)
	31
e50912db	32	const Int_t AliITSRecoParam::fgkLayersNotToSkip[AliITSgeomTGeo::kNLayers]={0,0,0,0,0,0};
	33	const Int_t AliITSRecoParam::fgkLastLayerToTrackTo=0;
	34	const Int_t AliITSRecoParam::fgkMaxDetectorPerLayer=1000;
	35	const Double_t AliITSRecoParam::fgkriw=80.0;
	36	const Double_t AliITSRecoParam::fgkdiw=0.0053;
	37	const Double_t AliITSRecoParam::fgkX0iw=30.0;
	38	const Double_t AliITSRecoParam::fgkrcd=61.0;
	39	const Double_t AliITSRecoParam::fgkdcd=0.0053;
	40	const Double_t AliITSRecoParam::fgkX0cd=30.0;
	41	const Double_t AliITSRecoParam::fgkyr=12.8;
	42	const Double_t AliITSRecoParam::fgkdr=0.03;
	43	const Double_t AliITSRecoParam::fgkzm=0.2;
	44	const Double_t AliITSRecoParam::fgkdm=0.40;
	45	const Double_t AliITSRecoParam::fgkrs=50.0;
	46	const Double_t AliITSRecoParam::fgkds=0.001;
	47	const Double_t AliITSRecoParam::fgkrInsideITSscreen=49.0;
	48	const Double_t AliITSRecoParam::fgkrInsideSPD1=3.7;
	49	const Double_t AliITSRecoParam::fgkrPipe=3.;
	50	const Double_t AliITSRecoParam::fgkrInsidePipe=2.7;
	51	const Double_t AliITSRecoParam::fgkrOutsidePipe=3.3;
	52	const Double_t AliITSRecoParam::fgkdPipe=0.0028;
	53	const Double_t AliITSRecoParam::fgkrInsideShield[2]={7.5,25.0};
	54	const Double_t AliITSRecoParam::fgkrOutsideShield[2]={10.5,30.0};
	55	const Double_t AliITSRecoParam::fgkdshield[2]={0.0097,0.0034};
	56	const Double_t AliITSRecoParam::fgkX0shield[2]={38.6,42.0};
	57	const Double_t AliITSRecoParam::fgkX0Air=21.82;
	58	const Double_t AliITSRecoParam::fgkX0Be=65.19;
	59	const Double_t AliITSRecoParam::fgkBoundaryWidth=0.2;
	60	const Double_t AliITSRecoParam::fgkDeltaXNeighbDets=0.5;
	61	const Double_t AliITSRecoParam::fgkSPDdetzlength=6.960; // 7.072-2*0.056
	62	const Double_t AliITSRecoParam::fgkSPDdetxlength=1.298; // 1.410-2*0.056
44347160	63
44347160	64	//_____________________________________________________________________________
6518a6c5	65	AliITSRecoParam::AliITSRecoParam() : AliDetectorRecoParam(),
f9119eb9	66	fTracker(0),
	67	fITSonly(kFALSE),
	68	fVertexer(0),
876026b6	69	fClusterFinder(0),
876026b6	70	fPID(0),
7203e11a	71	fVtxr3DZCutWide(0.),
	72	fVtxr3DRCutWide(0.),
	73	fVtxr3DZCutNarrow(0.),
	74	fVtxr3DRCutNarrow(0.),
	75	fVtxr3DPhiCutLoose(0.),
	76	fVtxr3DPhiCutTight(0.),
	77	fVtxr3DDCACut(0.),
e50912db	78	fMaxSnp(1.),
572f41f9	79	fNSigmaYLayerForRoadY(0),
	80	fNSigmaRoadY(0),
	81	fNSigmaZLayerForRoadZ(0),
	82	fNSigmaRoadZ(0),
	83	fNSigma2RoadZC(0),
	84	fNSigma2RoadYC(0),
	85	fNSigma2RoadZNonC(0),
	86	fNSigma2RoadYNonC(0),
1c97ce2f	87	fRoadMisal(0),
572f41f9	88	fMaxNormChi2NonCForHypothesis(0),
	89	fMaxChi2(0),
	90	fMaxRoad(0),
	91	fMaxChi2In(0),
	92	fChi2PerCluster(0),
	93	fXV(0),
	94	fYV(0),
	95	fZV(0),
	96	fSigmaXV(0),
	97	fSigmaYV(0),
	98	fSigmaZV(0),
	99	fVertexCut(0),
	100	fMaxDZforPrimTrk(0),
	101	fMaxDZToUseConstraint(0),
	102	fMaxDforV0dghtrForProlongation(0),
	103	fMaxDForProlongation(0),
	104	fMaxDZForProlongation(0),
	105	fMinPtForProlongation(0),
6518a6c5	106	fAddVirtualClustersInDeadZone(kFALSE),
572f41f9	107	fZWindowDeadZone(0),
	108	fSigmaXDeadZoneHit2(0),
	109	fSigmaZDeadZoneHit2(0),
	110	fXPassDeadZoneHits(0),
6518a6c5	111	fUseTGeoInTracker(3),
908753cb	112	fStepSizeTGeo(5.),
572f41f9	113	fAllowSharedClusters(kTRUE),
e50912db	114	fClusterErrorsParam(1),
4a66240a	115	fComputePlaneEff(kFALSE),
5fbd4fd6	116	fHistoPlaneEff(kFALSE),
58e8dc31	117	fUseTrackletsPlaneEff(kFALSE),
	118	fMCTrackletsPlaneEff(kFALSE),
	119	fBkgTrackletsPlaneEff(kFALSE),
	120	fTrackleterPhiWindowL1(0.10),
	121	fTrackleterZetaWindowL1(0.6),
	122	fUpdateOncePerEventPlaneEff(kTRUE),
	123	fMinContVtxPlaneEff(3),
0ed58a47	124	fIPlanePlaneEff(0),
275a301c	125	fReadPlaneEffFromOCDB(kFALSE),
0ed58a47	126	fMinPtPlaneEff(0),
	127	fMaxMissingClustersPlaneEff(0),
	128	fRequireClusterInOuterLayerPlaneEff(kFALSE),
	129	fRequireClusterInInnerLayerPlaneEff(kFALSE),
	130	fOnlyConstraintPlaneEff(kFALSE),
2755f080	131	fExtendedEtaAcceptance(kFALSE),
b406db2e	132	fUseBadZonesFromOCDB(kTRUE),
23197852	133	fUseSingleBadChannelsFromOCDB(kFALSE),
23197852	134	fMinFractionOfBadInRoad(0),
ae00569a	135	fAllowProlongationWithEmptyRoad(kFALSE),
12b1afb7	136	fInwardFlagSA(kFALSE),
	137	fOuterStartLayerSA(2),
	138	fInnerStartLayerSA(3),
	139	fMinNPointsSA(3),
2755f080	140	fFactorSAWindowSizes(1.),
c7d6d7b7	141	fNLoopsSA(33),
	142	fMinPhiSA(0.002),
	143	fMaxPhiSA(0.0145),
	144	fMinLambdaSA(0.003),
	145	fMaxLambdaSA(0.008),
3733ccd2	146	fMinClusterChargeSA(0.),
5a03f353	147	fSAOnePointTracks(kFALSE),
5a03f353	148	fSAUseAllClusters(kFALSE),
6518a6c5	149	fFindV0s(kTRUE),
b0160ebc	150	fStoreLikeSignV0s(kFALSE),
6518a6c5	151	fUseUnfoldingInClusterFinderSPD(kFALSE),
6518a6c5	152	fUseUnfoldingInClusterFinderSDD(kTRUE),
a86176e3	153	fUseUnfoldingInClusterFinderSSD(kTRUE),
b779c7dc	154	fUseSDDClusterSizeSelection(kFALSE),
b779c7dc	155	fMinClusterChargeSDD(0.),
7b116aa1	156	fUseChargeMatchingInClusterFinderSSD(kTRUE),
	157	fTrackleterOnlyOneTrackletPerC2(kTRUE),
	158	fTrackleterPhiWindow(0.08),
	159	fTrackleterZetaWindow(1.00),
	160	fTrackleterRemoveClustersFromOverlaps(kFALSE),
	161	fTrackleterPhiOverlapCut(0.005),
42ed6062	162	fTrackleterZetaOverlapCut(0.05),
ad7f2bfa	163	fUseCosmicRunShiftsSSD(kFALSE),
	164	fSPDRemoveNoisyFlag(kTRUE),
	165	fSPDRemoveDeadFlag(kTRUE)
44347160	166	{
	167	//
	168	// constructor
	169	//
6518a6c5	170	SetName("ITS");
	171	SetTitle("ITS");
	172
44347160	173	SetLayersParameters();
2755f080	174	for(Int_t i=0; i<AliITSgeomTGeo::kNLayers; i++) fLayersToSkip[i]=0;
6518a6c5	175	SetUseTGeoInTracker(3);
908753cb	176	SetStepSizeTGeo(5.);
afd25725	177	SetAllowSharedClusters(kTRUE);
afd25725	178	SetFindV0s(kTRUE);
6518a6c5	179	SetAddVirtualClustersInDeadZone(kFALSE);
572f41f9	180	SetUseAmplitudeInfo(kTRUE);
e50912db	181	SetClusterErrorsParam(1);
f9119eb9	182	SetClusterMisalError(0.);
7203e11a	183	SetVertexer3DDefaults();
44347160	184	}
	185	//_____________________________________________________________________________
	186	AliITSRecoParam::~AliITSRecoParam()
	187	{
	188	//
	189	// destructor
	190	//
	191	}
	192	//_____________________________________________________________________________
	193	AliITSRecoParam *AliITSRecoParam::GetHighFluxParam()
	194	{
	195	//
	196	// make default reconstruction parameters for hig flux env.
	197	//
	198	AliITSRecoParam *param = new AliITSRecoParam();
14bcd755	199	param->SetVertexerZ();
44347160	200
1aa655ec	201	// use of bads from OCDB
	202	param->SetUseBadZonesFromOCDB(kTRUE);
	203	param->SetUseSingleBadChannelsFromOCDB(kFALSE);
	204	// extended eta acceptance
	205	param->SetExtendedEtaAcceptance(kFALSE);
	206	// allow to skip layer if no cluster and no bad
	207	param->SetAllowProlongationWithEmptyRoad(kFALSE);
	208
9be1d1c7	209
44347160	210	param->fMaxSnp = 0.95;
	211
	212	param->fNSigmaYLayerForRoadY = 4.;
	213	param->fNSigmaRoadY = 7.5;
	214	param->fNSigmaZLayerForRoadZ = 4.;
	215	param->fNSigmaRoadZ = 7.5;
	216
	217	param->fNSigma2RoadZC = 60.; //7.75^2
	218	param->fNSigma2RoadYC = 60.; //7.75^2
	219	param->fNSigma2RoadZNonC = 50.; //7.07^2
	220	param->fNSigma2RoadYNonC = 50.; //7.07^2
	221
	222	param->fMaxChi2PerCluster[0] = 11.; //7
	223	param->fMaxChi2PerCluster[1] = 12.; //5
	224	param->fMaxChi2PerCluster[2] = 12.; //8
	225	param->fMaxChi2PerCluster[3] = 5.; //8
	226	param->fMaxChi2PerCluster[4] = 12.; //6.5
	227
	228	param->fMaxNormChi2NonC[0] = 7.;
	229	param->fMaxNormChi2NonC[1] = 8.;
	230	param->fMaxNormChi2NonC[2] = 8.;
	231	param->fMaxNormChi2NonC[3] = 11.;
	232	param->fMaxNormChi2NonC[4] = 14.;
	233	param->fMaxNormChi2NonC[5] = 25.;
	234
	235	param->fMaxNormChi2C[0] = 11.;
	236	param->fMaxNormChi2C[1] = 13.;
	237	param->fMaxNormChi2C[2] = 15.;
	238	param->fMaxNormChi2C[3] = 18.;
	239	param->fMaxNormChi2C[4] = 30.;
	240	param->fMaxNormChi2C[5] = 35.;
afd25725	241
afd25725	242	param->fMaxNormChi2NonCForHypothesis = 7.;
44347160	243
	244	param->fMaxChi2 = 35.;
	245
	246	param->fMaxChi2s[0] = 25.; //40
	247	param->fMaxChi2s[1] = 25.; //40
	248	param->fMaxChi2s[2] = 25.; //40
	249	param->fMaxChi2s[3] = 25.; //40
	250	param->fMaxChi2s[4] = 40.; //40
	251	param->fMaxChi2s[5] = 50.; //40
	252
	253	param->fMaxRoad = 6.;
	254
	255	// not used
	256	param->fMaxChi2In = 16.;
	257
	258	param->fMaxChi2sR[0] = 10.;
	259	param->fMaxChi2sR[1] = 10.;
	260	param->fMaxChi2sR[2] = 10.;
	261	param->fMaxChi2sR[3] = 10.;
	262	param->fMaxChi2sR[4] = 30.;
	263	param->fMaxChi2sR[5] = 40.;
	264
	265	param->fChi2PerCluster = 9.;
	266	// not used
	267
	268	param->fXV = 0.;
	269	param->fYV = 0.;
	270	param->fZV = 0.;
	271	param->fSigmaXV = 0.0050;
	272	param->fSigmaYV = 0.0050;
	273	param->fSigmaZV = 0.0100;
	274
	275	param->fVertexCut = 25.;
afd25725	276
	277	param->fMaxDZforPrimTrk = 0.4;
	278	param->fMaxDZToUseConstraint = 3.;
	279
	280	param->fMaxDforV0dghtrForProlongation = 30.;
6518a6c5	281	param->fMaxDForProlongation = 40.;//10.;
6518a6c5	282	param->fMaxDZForProlongation = 60.;//20.;
afd25725	283	param->fMinPtForProlongation = 0.120;
	284
	285	param->fZWindowDeadZone = 2.0;
	286	param->fSigmaXDeadZoneHit2 = 0.004/12.;
	287	param->fSigmaZDeadZoneHit2 = 0.001/12.;
	288	param->fXPassDeadZoneHits = 0.018;
44347160	289
	290	return param;
	291	}
	292	//_____________________________________________________________________________
	293	AliITSRecoParam *AliITSRecoParam::GetLowFluxParam()
	294	{
	295	//
	296	// make default reconstruction parameters for low flux env.
	297	//
c7d6d7b7	298	AliITSRecoParam *param = new AliITSRecoParam();
7203e11a	299	param->SetVertexer3DDefaults();
c7d6d7b7	300
1aa655ec	301	// full use of bads from OCDB
	302	param->SetUseBadZonesFromOCDB(kTRUE);
	303	param->SetUseSingleBadChannelsFromOCDB(kTRUE);
	304	// extended eta acceptance
	305	param->SetExtendedEtaAcceptance(kTRUE);
	306	// allow to skip layer if no cluster and no bad
	307	param->SetAllowProlongationWithEmptyRoad(kTRUE);
	308
9be1d1c7	309
c7d6d7b7	310	param->fMaxSnp = 0.95;
	311
	312	param->fNSigmaYLayerForRoadY = 4.;
	313	param->fNSigmaRoadY = 7.5;
	314	param->fNSigmaZLayerForRoadZ = 4.;
	315	param->fNSigmaRoadZ = 7.5;
	316
	317	param->fNSigma2RoadZC = 60.; //7.75^2
	318	param->fNSigma2RoadYC = 60.; //7.75^2
	319	param->fNSigma2RoadZNonC = 50.; //7.07^2
	320	param->fNSigma2RoadYNonC = 50.; //7.07^2
	321
	322	param->fMaxChi2PerCluster[0] = 11.; //7
	323	param->fMaxChi2PerCluster[1] = 12.; //5
	324	param->fMaxChi2PerCluster[2] = 12.; //8
	325	param->fMaxChi2PerCluster[3] = 5.; //8
	326	param->fMaxChi2PerCluster[4] = 12.; //6.5
	327
	328	param->fMaxNormChi2NonC[0] = 7.;
	329	param->fMaxNormChi2NonC[1] = 8.;
	330	param->fMaxNormChi2NonC[2] = 8.;
	331	param->fMaxNormChi2NonC[3] = 11.;
	332	param->fMaxNormChi2NonC[4] = 14.;
	333	param->fMaxNormChi2NonC[5] = 25.;
	334
	335	param->fMaxNormChi2C[0] = 11.;
	336	param->fMaxNormChi2C[1] = 13.;
	337	param->fMaxNormChi2C[2] = 15.;
	338	param->fMaxNormChi2C[3] = 18.;
	339	param->fMaxNormChi2C[4] = 30.;
	340	param->fMaxNormChi2C[5] = 35.;
	341
	342	param->fMaxNormChi2NonCForHypothesis = 7.;
	343
	344	param->fMaxChi2 = 35.;
	345
	346	param->fMaxChi2s[0] = 25.; //40
	347	param->fMaxChi2s[1] = 25.; //40
	348	param->fMaxChi2s[2] = 25.; //40
	349	param->fMaxChi2s[3] = 25.; //40
	350	param->fMaxChi2s[4] = 40.; //40
	351	param->fMaxChi2s[5] = 50.; //40
	352
	353	param->fMaxRoad = 6.;
	354
	355	// not used
	356	param->fMaxChi2In = 16.;
	357
	358	param->fMaxChi2sR[0] = 10.;
	359	param->fMaxChi2sR[1] = 10.;
	360	param->fMaxChi2sR[2] = 10.;
	361	param->fMaxChi2sR[3] = 10.;
	362	param->fMaxChi2sR[4] = 30.;
	363	param->fMaxChi2sR[5] = 40.;
	364
	365	param->fChi2PerCluster = 9.;
	366	// not used
	367
	368	param->fXV = 0.;
	369	param->fYV = 0.;
	370	param->fZV = 0.;
	371	param->fSigmaXV = 0.0050;
	372	param->fSigmaYV = 0.0050;
	373	param->fSigmaZV = 0.0100;
374
375	param->fVertexCut = 25.;
376
377	param->fMaxDZforPrimTrk = 0.4;
378	param->fMaxDZToUseConstraint = 3.;
379
380	param->fMaxDforV0dghtrForProlongation = 30.;
381	param->fMaxDForProlongation = 40.;//10.;
382	param->fMaxDZForProlongation = 60.;//20.;
383	param->fMinPtForProlongation = 0.120;
384
385	param->fZWindowDeadZone = 2.0;
386	param->fSigmaXDeadZoneHit2 = 0.004/12.;
387	param->fSigmaZDeadZoneHit2 = 0.001/12.;
388	param->fXPassDeadZoneHits = 0.018;
3733ccd2	389	param->SetNLoopsSA(20);
	390	param->fMaxPhiSA = 0.07;
	391	param->fMaxLambdaSA = 0.04;
c7d6d7b7	392	return param;
44347160	393	}
	394	//_____________________________________________________________________________
	395	AliITSRecoParam *AliITSRecoParam::GetCosmicTestParam()
	396	{
	397	//
	398	// make default reconstruction parameters for cosmics
	399	//
e340bb86	400	AliITSRecoParam *param = new AliITSRecoParam();
e340bb86	401
f9119eb9	402	// vertexer for cosmics
	403	param->SetVertexer(2);
	404
9f9cae94	405	// find independently ITS SA tracks
	406	param->SetSAUseAllClusters();
	407	param->SetOuterStartLayerSA(AliITSgeomTGeo::GetNLayers()-2);
f9119eb9	408
2755f080	409	// larger seach windows for SA (in case of large misalignments)
2755f080	410	param->SetFactorSAWindowSizes(3.);
e340bb86	411
	412	param->fMaxSnp = 0.95;
	413
	414	param->fNSigmaYLayerForRoadY = 4.;
	415	param->fNSigmaRoadY = 7.5;
	416	param->fNSigmaZLayerForRoadZ = 4.;
	417	param->fNSigmaRoadZ = 7.5;
	418
	419	param->fNSigma2RoadZC = 60.; //7.75^2
	420	param->fNSigma2RoadYC = 60.; //7.75^2
	421	param->fNSigma2RoadZNonC = 50.; //7.07^2
	422	param->fNSigma2RoadYNonC = 50.; //7.07^2
	423
	424	param->fMaxChi2PerCluster[0] = 11.; //7
	425	param->fMaxChi2PerCluster[1] = 12.; //5
	426	param->fMaxChi2PerCluster[2] = 12.; //8
	427	param->fMaxChi2PerCluster[3] = 5.; //8
	428	param->fMaxChi2PerCluster[4] = 12.; //6.5
	429
	430	param->fMaxNormChi2NonC[0] = 7.;
	431	param->fMaxNormChi2NonC[1] = 8.;
	432	param->fMaxNormChi2NonC[2] = 8.;
	433	param->fMaxNormChi2NonC[3] = 11.;
	434	param->fMaxNormChi2NonC[4] = 14.;
	435	param->fMaxNormChi2NonC[5] = 25.;
	436
	437	param->fMaxNormChi2C[0] = 11.;
	438	param->fMaxNormChi2C[1] = 13.;
	439	param->fMaxNormChi2C[2] = 15.;
	440	param->fMaxNormChi2C[3] = 18.;
	441	param->fMaxNormChi2C[4] = 30.;
	442	param->fMaxNormChi2C[5] = 35.;
	443
	444	param->fMaxNormChi2NonCForHypothesis = 7.;
	445
	446	param->fMaxChi2 = 35.;
	447
	448	param->fMaxChi2s[0] = 25.; //40
	449	param->fMaxChi2s[1] = 25.; //40
	450	param->fMaxChi2s[2] = 25.; //40
	451	param->fMaxChi2s[3] = 25.; //40
	452	param->fMaxChi2s[4] = 40.; //40
	453	param->fMaxChi2s[5] = 50.; //40
	454
	455	param->fMaxRoad = 6.;
	456
	457	// not used
	458	param->fMaxChi2In = 16.;
	459
	460	param->fMaxChi2sR[0] = 10.;
	461	param->fMaxChi2sR[1] = 10.;
	462	param->fMaxChi2sR[2] = 10.;
	463	param->fMaxChi2sR[3] = 10.;
	464	param->fMaxChi2sR[4] = 30.;
	465	param->fMaxChi2sR[5] = 40.;
	466
	467	param->fChi2PerCluster = 9.;
	468	// not used
	469
	470	param->fXV = 0.;
	471	param->fYV = 0.;
	472	param->fZV = 0.;
	473	param->fSigmaXV = 0.0050;
	474	param->fSigmaYV = 0.0050;
475	param->fSigmaZV = 0.0100;
476
477	param->fVertexCut = 25.;
478
479	param->fMaxDZforPrimTrk = 0.4;
480	param->fMaxDZToUseConstraint = 3.;
481
482	param->fMaxDforV0dghtrForProlongation = 30.;
483	param->fMaxDForProlongation = 10.;
484	param->fMaxDZForProlongation = 20.;
485	param->fMinPtForProlongation = 0.120;
486
487	param->fZWindowDeadZone = 2.0;
488	param->fSigmaXDeadZoneHit2 = 0.004/12.;
489	param->fSigmaZDeadZoneHit2 = 0.001/12.;
490	param->fXPassDeadZoneHits = 0.018;
491
492
493	return param;
44347160	494	}
44347160	495	//_____________________________________________________________________________
ae00569a	496	AliITSRecoParam *AliITSRecoParam::GetPlaneEffParam(Int_t i)
	497	{
	498	//
58e8dc31	499	// make special reconstruction parameters for Plane Efficiency study on layer i (0,5)
	500	//
	501	// if i=-1, then the evaluation for both pixel layers is tried with the tracklet method
ae00569a	502	//
58e8dc31	503	if (i<-1 \|\| i>=AliITSgeomTGeo::kNLayers) {
58e8dc31	504	printf("AliITSRecoParam::GetPlaneEffParam: index of ITS Plane nor in the range [0,5] neither =-1\n");
0ed58a47	505	printf("returning null pointer");
	506	return NULL;
	507	}
58e8dc31	508	if(i>=0) { // Method using tracks (remove given plane from tracking)
	509	AliITSRecoParam *param;
	510	param = GetHighFluxParam();
	511	param->SetComputePlaneEff();
	512	param->SetLayerToSkip(i);
	513	param->SetIPlanePlaneEff(i);
	514	// optimized setting for SPD0 (i==0)
	515	if (i==0 \|\| i==1) {
	516	param->fMinPtPlaneEff = 0.200; // high pt particles
	517	param->fMaxMissingClustersPlaneEff = 1; // at most 1 layer out of 5 without cluster
	518	param->fRequireClusterInOuterLayerPlaneEff = kTRUE; // cluster on SPD1 must be
	519	//param->fOnlyConstraintPlaneEff = kTRUE;
	520	}
	521	if (i==2 \|\| i==3) {
	522	param->fMinPtPlaneEff = 0.200; // high pt particles
	523	param->fMaxMissingClustersPlaneEff = 1; // at most 1 layer out of 5 without cluster
	524	param->fRequireClusterInOuterLayerPlaneEff = kTRUE;
	525	//param->fOnlyConstraintPlaneEff = kTRUE;
	526	}
	527	if (i==4) {
	528	param->fMinPtPlaneEff = 0.200; // high pt particles
	529	param->fMaxMissingClustersPlaneEff = 0; // at most 1 layer out of 5 without cluster
	530	param->fRequireClusterInOuterLayerPlaneEff = kTRUE;
	531	//param->fOnlyConstraintPlaneEff = kTRUE;
	532	}
	533	if (i==5) {
	534	param->fMinPtPlaneEff = 0.200; // high pt particles
	535	}
	536	//
	537	return param;
0ed58a47	538	}
58e8dc31	539	else if (i==-1) { // Method using tracklets
	540	AliITSRecoParam *param;
	541	param = GetLowFluxParam();
	542	param->SetIPlanePlaneEff(i);
	543	param->SetComputePlaneEff(kTRUE,kFALSE);
	544	param->SetUseTrackletsPlaneEff(kTRUE);
	545	param->SetTrackleterPhiWindow(0.07);
	546	param->SetTrackleterZetaWindow(0.4);
	547	param->SetTrackleterPhiWindowL1(0.10);
	548	param->SetTrackleterZetaWindowL1(0.6);
	549	param->SetUpdateOncePerEventPlaneEff(kTRUE);
	550	param->SetMinContVtxPlaneEff(3);
	551	return param;
0ed58a47	552	}
5bd7ec3a	553	else {
	554	AliErrorGeneral("AliITSRecoParam",Form("Unrecognised value of i %d\n",i));
	555	return 0;
	556	}
ae00569a	557	}
ae00569a	558	//_____________________________________________________________________________
44347160	559	void AliITSRecoParam::SetLayersParameters()
	560	{
	561	//
	562	// number of layers and layers spatial resolutions
	563	//
	564
	565	// spatial resolutions of the detectors
	566	// y: 12 12 38 38 20 20 micron
	567	fSigmaY2[0]=1.44e-6;
	568	fSigmaY2[1]=1.44e-6;
	569	fSigmaY2[2]=1.444e-5;
	570	fSigmaY2[3]=1.444e-5;
	571	fSigmaY2[4]=4.0e-6;
	572	fSigmaY2[5]=4.0e-6;
	573	// z: 120 120 28 28 830 830 micron
	574	fSigmaZ2[0]=1.44e-4;
	575	fSigmaZ2[1]=1.44e-4;
	576	fSigmaZ2[2]=7.84e-6;
	577	fSigmaZ2[3]=7.84e-6;
	578	fSigmaZ2[4]=6.889e-3;
	579	fSigmaZ2[5]=6.889e-3;
	580
	581	return;
	582	}
ed446fa3	583	//_____________________________________________________________________________
	584	void AliITSRecoParam::PrintParameters() const
	585	{
	586	//
	587	// print parameters
	588	//
	589
	590	printf("============================= AliITSRecoParam::PrintParameters ");
f9119eb9	591	printf("============================= \n\n");
ed446fa3	592	for(Int_t i=0; i<AliITSgeomTGeo::kNLayers; i++) {
	593	if(!fLayersToSkip[i]) {
	594	printf("ITS Traking: using layer %d\n",i);
	595	} else {
	596	printf("ITS Traking: skipping layer %d\n",i);
	597	}
	598	}
	599	for(Int_t i=0; i<AliITSgeomTGeo::kNLayers; i++) {
	600	if(fUseAmplitudeInfo[i]) {
	601	printf("ITS Traking: use amplitude info for layer %d\n",i);
	602	} else {
	603	printf("ITS Traking: don't use amplitude info for layer %d\n",i);
	604	}
	605	}
	606	for(Int_t i=0; i<AliITSgeomTGeo::kNLayers; i++)
f9119eb9	607	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	608
	609
	610	Dump();
	611
	612	return;
	613	}