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