[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
44347160	1	/**************************************************************************
572f41f9	2	* Copyright(c) 2007-2009, ALICE Experiment at CERN, All rights reserved. *
44347160	3	* *
	4	* Author: The ALICE Off-line Project. *
	5	* Contributors are mentioned in the code where appropriate. *
	6	* *
	7	* Permission to use, copy, modify and distribute this software and its *
	8	* documentation strictly for non-commercial purposes is hereby granted *
	9	* without fee, provided that the above copyright notice appears in all *
	10	* copies and that both the copyright notice and this permission notice *
	11	* appear in the supporting documentation. The authors make no claims *
	12	* about the suitability of this software for any purpose. It is *
	13	* provided "as is" without express or implied warranty. *
	14	**************************************************************************/
	15
572f41f9	16	#include "AliITSRecoParam.h"
	17
	18	/* $Id$ */
44347160	19
	20	///////////////////////////////////////////////////////////////////////////////
	21	// //
	22	// Class with ITS reconstruction parameters //
	23	// Origin: andrea.dainese@lnl.infn.it //
	24	// //
	25	///////////////////////////////////////////////////////////////////////////////
	26
	27
44347160	28
	29	ClassImp(AliITSRecoParam)
	30
e50912db	31	const Int_t AliITSRecoParam::fgkLayersNotToSkip[AliITSgeomTGeo::kNLayers]={0,0,0,0,0,0};
	32	const Int_t AliITSRecoParam::fgkLastLayerToTrackTo=0;
	33	const Int_t AliITSRecoParam::fgkMaxDetectorPerLayer=1000;
	34	const Double_t AliITSRecoParam::fgkriw=80.0;
	35	const Double_t AliITSRecoParam::fgkdiw=0.0053;
	36	const Double_t AliITSRecoParam::fgkX0iw=30.0;
	37	const Double_t AliITSRecoParam::fgkrcd=61.0;
	38	const Double_t AliITSRecoParam::fgkdcd=0.0053;
	39	const Double_t AliITSRecoParam::fgkX0cd=30.0;
	40	const Double_t AliITSRecoParam::fgkyr=12.8;
	41	const Double_t AliITSRecoParam::fgkdr=0.03;
	42	const Double_t AliITSRecoParam::fgkzm=0.2;
	43	const Double_t AliITSRecoParam::fgkdm=0.40;
	44	const Double_t AliITSRecoParam::fgkrs=50.0;
	45	const Double_t AliITSRecoParam::fgkds=0.001;
	46	const Double_t AliITSRecoParam::fgkrInsideITSscreen=49.0;
	47	const Double_t AliITSRecoParam::fgkrInsideSPD1=3.7;
	48	const Double_t AliITSRecoParam::fgkrPipe=3.;
	49	const Double_t AliITSRecoParam::fgkrInsidePipe=2.7;
	50	const Double_t AliITSRecoParam::fgkrOutsidePipe=3.3;
	51	const Double_t AliITSRecoParam::fgkdPipe=0.0028;
	52	const Double_t AliITSRecoParam::fgkrInsideShield[2]={7.5,25.0};
	53	const Double_t AliITSRecoParam::fgkrOutsideShield[2]={10.5,30.0};
	54	const Double_t AliITSRecoParam::fgkdshield[2]={0.0097,0.0034};
	55	const Double_t AliITSRecoParam::fgkX0shield[2]={38.6,42.0};
	56	const Double_t AliITSRecoParam::fgkX0Air=21.82;
	57	const Double_t AliITSRecoParam::fgkX0Be=65.19;
	58	const Double_t AliITSRecoParam::fgkBoundaryWidth=0.2;
	59	const Double_t AliITSRecoParam::fgkDeltaXNeighbDets=0.5;
	60	const Double_t AliITSRecoParam::fgkSPDdetzlength=6.960; // 7.072-2*0.056
	61	const Double_t AliITSRecoParam::fgkSPDdetxlength=1.298; // 1.410-2*0.056
44347160	62
44347160	63	//_____________________________________________________________________________
6518a6c5	64	AliITSRecoParam::AliITSRecoParam() : AliDetectorRecoParam(),
f9119eb9	65	fTracker(0),
	66	fITSonly(kFALSE),
	67	fVertexer(0),
876026b6	68	fClusterFinder(0),
876026b6	69	fPID(0),
7203e11a	70	fVtxr3DZCutWide(0.),
	71	fVtxr3DRCutWide(0.),
	72	fVtxr3DZCutNarrow(0.),
	73	fVtxr3DRCutNarrow(0.),
	74	fVtxr3DPhiCutLoose(0.),
	75	fVtxr3DPhiCutTight(0.),
	76	fVtxr3DDCACut(0.),
e50912db	77	fMaxSnp(1.),
572f41f9	78	fNSigmaYLayerForRoadY(0),
	79	fNSigmaRoadY(0),
	80	fNSigmaZLayerForRoadZ(0),
	81	fNSigmaRoadZ(0),
	82	fNSigma2RoadZC(0),
	83	fNSigma2RoadYC(0),
	84	fNSigma2RoadZNonC(0),
	85	fNSigma2RoadYNonC(0),
1c97ce2f	86	fRoadMisal(0),
572f41f9	87	fMaxNormChi2NonCForHypothesis(0),
	88	fMaxChi2(0),
	89	fMaxRoad(0),
	90	fMaxChi2In(0),
	91	fChi2PerCluster(0),
	92	fXV(0),
	93	fYV(0),
	94	fZV(0),
	95	fSigmaXV(0),
	96	fSigmaYV(0),
	97	fSigmaZV(0),
	98	fVertexCut(0),
	99	fMaxDZforPrimTrk(0),
	100	fMaxDZToUseConstraint(0),
	101	fMaxDforV0dghtrForProlongation(0),
	102	fMaxDForProlongation(0),
	103	fMaxDZForProlongation(0),
	104	fMinPtForProlongation(0),
6518a6c5	105	fAddVirtualClustersInDeadZone(kFALSE),
572f41f9	106	fZWindowDeadZone(0),
	107	fSigmaXDeadZoneHit2(0),
	108	fSigmaZDeadZoneHit2(0),
	109	fXPassDeadZoneHits(0),
6518a6c5	110	fUseTGeoInTracker(3),
572f41f9	111	fAllowSharedClusters(kTRUE),
e50912db	112	fClusterErrorsParam(1),
4a66240a	113	fComputePlaneEff(kFALSE),
5fbd4fd6	114	fHistoPlaneEff(kFALSE),
0ed58a47	115	fIPlanePlaneEff(0),
275a301c	116	fReadPlaneEffFromOCDB(kFALSE),
0ed58a47	117	fMinPtPlaneEff(0),
	118	fMaxMissingClustersPlaneEff(0),
	119	fRequireClusterInOuterLayerPlaneEff(kFALSE),
	120	fRequireClusterInInnerLayerPlaneEff(kFALSE),
	121	fOnlyConstraintPlaneEff(kFALSE),
2755f080	122	fExtendedEtaAcceptance(kFALSE),
b406db2e	123	fUseBadZonesFromOCDB(kTRUE),
23197852	124	fUseSingleBadChannelsFromOCDB(kFALSE),
23197852	125	fMinFractionOfBadInRoad(0),
ae00569a	126	fAllowProlongationWithEmptyRoad(kFALSE),
9f9cae94	127	fOuterStartLayerSA(0),
2755f080	128	fFactorSAWindowSizes(1.),
c7d6d7b7	129	fNLoopsSA(33),
	130	fMinPhiSA(0.002),
	131	fMaxPhiSA(0.0145),
	132	fMinLambdaSA(0.003),
	133	fMaxLambdaSA(0.008),
5a03f353	134	fSAOnePointTracks(kFALSE),
5a03f353	135	fSAUseAllClusters(kFALSE),
6518a6c5	136	fFindV0s(kTRUE),
	137	fUseUnfoldingInClusterFinderSPD(kFALSE),
	138	fUseUnfoldingInClusterFinderSDD(kTRUE),
a86176e3	139	fUseUnfoldingInClusterFinderSSD(kTRUE),
7b116aa1	140	fUseChargeMatchingInClusterFinderSSD(kTRUE),
	141	fTrackleterOnlyOneTrackletPerC2(kTRUE),
	142	fTrackleterPhiWindow(0.08),
	143	fTrackleterZetaWindow(1.00),
	144	fTrackleterRemoveClustersFromOverlaps(kFALSE),
	145	fTrackleterPhiOverlapCut(0.005),
42ed6062	146	fTrackleterZetaOverlapCut(0.05),
42ed6062	147	fUseCosmicRunShiftsSSD(kFALSE)
44347160	148	{
	149	//
	150	// constructor
	151	//
6518a6c5	152	SetName("ITS");
	153	SetTitle("ITS");
	154
44347160	155	SetLayersParameters();
2755f080	156	for(Int_t i=0; i<AliITSgeomTGeo::kNLayers; i++) fLayersToSkip[i]=0;
6518a6c5	157	SetUseTGeoInTracker(3);
afd25725	158	SetAllowSharedClusters(kTRUE);
afd25725	159	SetFindV0s(kTRUE);
6518a6c5	160	SetAddVirtualClustersInDeadZone(kFALSE);
572f41f9	161	SetUseAmplitudeInfo(kTRUE);
e50912db	162	SetClusterErrorsParam(1);
f9119eb9	163	SetClusterMisalError(0.);
7203e11a	164	SetVertexer3DDefaults();
44347160	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();
14bcd755	180	param->SetVertexerZ();
44347160	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.;
afd25725	213
afd25725	214	param->fMaxNormChi2NonCForHypothesis = 7.;
44347160	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.;
afd25725	248
	249	param->fMaxDZforPrimTrk = 0.4;
	250	param->fMaxDZToUseConstraint = 3.;
	251
	252	param->fMaxDforV0dghtrForProlongation = 30.;
6518a6c5	253	param->fMaxDForProlongation = 40.;//10.;
6518a6c5	254	param->fMaxDZForProlongation = 60.;//20.;
afd25725	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;
44347160	261
	262	return param;
	263	}
	264	//_____________________________________________________________________________
	265	AliITSRecoParam *AliITSRecoParam::GetLowFluxParam()
	266	{
	267	//
	268	// make default reconstruction parameters for low flux env.
	269	//
c7d6d7b7	270	AliITSRecoParam *param = new AliITSRecoParam();
7203e11a	271	param->SetVertexer3DDefaults();
c7d6d7b7	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);
7203e11a	353
c7d6d7b7	354	return param;
44347160	355	}
	356	//_____________________________________________________________________________
	357	AliITSRecoParam *AliITSRecoParam::GetCosmicTestParam()
	358	{
	359	//
	360	// make default reconstruction parameters for cosmics
	361	//
e340bb86	362	AliITSRecoParam *param = new AliITSRecoParam();
e340bb86	363
f9119eb9	364	// vertexer for cosmics
	365	param->SetVertexer(2);
	366
9f9cae94	367	// find independently ITS SA tracks
	368	param->SetSAUseAllClusters();
	369	param->SetOuterStartLayerSA(AliITSgeomTGeo::GetNLayers()-2);
f9119eb9	370
2755f080	371	// larger seach windows for SA (in case of large misalignments)
2755f080	372	param->SetFactorSAWindowSizes(3.);
e340bb86	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;
44347160	456	}
44347160	457	//_____________________________________________________________________________
ae00569a	458	AliITSRecoParam *AliITSRecoParam::GetPlaneEffParam(Int_t i)
	459	{
	460	//
	461	// make special reconstruction parameters for Plane Efficiency study on layer i
	462	//
0ed58a47	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	}
ae00569a	468	AliITSRecoParam *param;
	469	param = GetHighFluxParam();
	470	param->SetComputePlaneEff();
	471	param->SetLayerToSkip(i);
0ed58a47	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	//
ae00569a	496	return param;
	497	}
	498	//_____________________________________________________________________________
44347160	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	}
ed446fa3	523	//_____________________________________________________________________________
	524	void AliITSRecoParam::PrintParameters() const
	525	{
	526	//
	527	// print parameters
	528	//
	529
	530	printf("============================= AliITSRecoParam::PrintParameters ");
f9119eb9	531	printf("============================= \n\n");
ed446fa3	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++)
f9119eb9	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]);
ed446fa3	548
	549
	550	Dump();
	551
	552	return;
	553	}