[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(),
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),
fReadPlaneEffFromOCDB(kFALSE),
fExtendedEtaAcceptance(kFALSE),
fUseBadZonesFromOCDB(kFALSE),
fUseSingleBadChannelsFromOCDB(kFALSE),
fMinFractionOfBadInRoad(0),
fAllowProlongationWithEmptyRoad(kFALSE),
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)
{
  //
  // 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);
}
//_____________________________________________________________________________
AliITSRecoParam::~AliITSRecoParam() 
{
  //
  // destructor
  //  
}
//_____________________________________________________________________________
AliITSRecoParam *AliITSRecoParam::GetHighFluxParam() 
{
  //
  // make default reconstruction  parameters for hig  flux env.
  //
  AliITSRecoParam *param = new AliITSRecoParam();

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

  // 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
  //
  AliITSRecoParam *param;
  param = GetHighFluxParam();
  param->SetComputePlaneEff();
  param->SetLayerToSkip(i);
  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, sigma Z2 %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],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(),
e50912db	65	fMaxSnp(1.),
572f41f9	66	fNSigmaYLayerForRoadY(0),
	67	fNSigmaRoadY(0),
	68	fNSigmaZLayerForRoadZ(0),
	69	fNSigmaRoadZ(0),
	70	fNSigma2RoadZC(0),
	71	fNSigma2RoadYC(0),
	72	fNSigma2RoadZNonC(0),
	73	fNSigma2RoadYNonC(0),
1c97ce2f	74	fRoadMisal(0),
572f41f9	75	fMaxNormChi2NonCForHypothesis(0),
	76	fMaxChi2(0),
	77	fMaxRoad(0),
	78	fMaxChi2In(0),
	79	fChi2PerCluster(0),
	80	fXV(0),
	81	fYV(0),
	82	fZV(0),
	83	fSigmaXV(0),
	84	fSigmaYV(0),
	85	fSigmaZV(0),
	86	fVertexCut(0),
	87	fMaxDZforPrimTrk(0),
	88	fMaxDZToUseConstraint(0),
	89	fMaxDforV0dghtrForProlongation(0),
	90	fMaxDForProlongation(0),
	91	fMaxDZForProlongation(0),
	92	fMinPtForProlongation(0),
6518a6c5	93	fAddVirtualClustersInDeadZone(kFALSE),
572f41f9	94	fZWindowDeadZone(0),
	95	fSigmaXDeadZoneHit2(0),
	96	fSigmaZDeadZoneHit2(0),
	97	fXPassDeadZoneHits(0),
6518a6c5	98	fUseTGeoInTracker(3),
572f41f9	99	fAllowSharedClusters(kTRUE),
e50912db	100	fClusterErrorsParam(1),
4a66240a	101	fComputePlaneEff(kFALSE),
5fbd4fd6	102	fHistoPlaneEff(kFALSE),
275a301c	103	fReadPlaneEffFromOCDB(kFALSE),
2755f080	104	fExtendedEtaAcceptance(kFALSE),
23197852	105	fUseBadZonesFromOCDB(kFALSE),
	106	fUseSingleBadChannelsFromOCDB(kFALSE),
	107	fMinFractionOfBadInRoad(0),
ae00569a	108	fAllowProlongationWithEmptyRoad(kFALSE),
2755f080	109	fFactorSAWindowSizes(1.),
c7d6d7b7	110	fNLoopsSA(33),
	111	fMinPhiSA(0.002),
	112	fMaxPhiSA(0.0145),
	113	fMinLambdaSA(0.003),
	114	fMaxLambdaSA(0.008),
5a03f353	115	fSAOnePointTracks(kFALSE),
5a03f353	116	fSAUseAllClusters(kFALSE),
6518a6c5	117	fFindV0s(kTRUE),
	118	fUseUnfoldingInClusterFinderSPD(kFALSE),
	119	fUseUnfoldingInClusterFinderSDD(kTRUE),
a86176e3	120	fUseUnfoldingInClusterFinderSSD(kTRUE),
a86176e3	121	fUseChargeMatchingInClusterFinderSSD(kTRUE)
44347160	122	{
	123	//
	124	// constructor
	125	//
6518a6c5	126	SetName("ITS");
	127	SetTitle("ITS");
	128
44347160	129	SetLayersParameters();
2755f080	130	for(Int_t i=0; i<AliITSgeomTGeo::kNLayers; i++) fLayersToSkip[i]=0;
6518a6c5	131	SetUseTGeoInTracker(3);
afd25725	132	SetAllowSharedClusters(kTRUE);
afd25725	133	SetFindV0s(kTRUE);
6518a6c5	134	SetAddVirtualClustersInDeadZone(kFALSE);
572f41f9	135	SetUseAmplitudeInfo(kTRUE);
e50912db	136	SetClusterErrorsParam(1);
44347160	137	}
	138	//_____________________________________________________________________________
	139	AliITSRecoParam::~AliITSRecoParam()
	140	{
	141	//
	142	// destructor
	143	//
	144	}
	145	//_____________________________________________________________________________
	146	AliITSRecoParam *AliITSRecoParam::GetHighFluxParam()
	147	{
	148	//
	149	// make default reconstruction parameters for hig flux env.
	150	//
	151	AliITSRecoParam *param = new AliITSRecoParam();
	152
	153	param->fMaxSnp = 0.95;
	154
	155	param->fNSigmaYLayerForRoadY = 4.;
	156	param->fNSigmaRoadY = 7.5;
	157	param->fNSigmaZLayerForRoadZ = 4.;
	158	param->fNSigmaRoadZ = 7.5;
	159
	160	param->fNSigma2RoadZC = 60.; //7.75^2
	161	param->fNSigma2RoadYC = 60.; //7.75^2
	162	param->fNSigma2RoadZNonC = 50.; //7.07^2
	163	param->fNSigma2RoadYNonC = 50.; //7.07^2
	164
	165	param->fMaxChi2PerCluster[0] = 11.; //7
	166	param->fMaxChi2PerCluster[1] = 12.; //5
	167	param->fMaxChi2PerCluster[2] = 12.; //8
	168	param->fMaxChi2PerCluster[3] = 5.; //8
	169	param->fMaxChi2PerCluster[4] = 12.; //6.5
	170
	171	param->fMaxNormChi2NonC[0] = 7.;
	172	param->fMaxNormChi2NonC[1] = 8.;
	173	param->fMaxNormChi2NonC[2] = 8.;
	174	param->fMaxNormChi2NonC[3] = 11.;
	175	param->fMaxNormChi2NonC[4] = 14.;
	176	param->fMaxNormChi2NonC[5] = 25.;
	177
	178	param->fMaxNormChi2C[0] = 11.;
	179	param->fMaxNormChi2C[1] = 13.;
	180	param->fMaxNormChi2C[2] = 15.;
	181	param->fMaxNormChi2C[3] = 18.;
	182	param->fMaxNormChi2C[4] = 30.;
	183	param->fMaxNormChi2C[5] = 35.;
afd25725	184
afd25725	185	param->fMaxNormChi2NonCForHypothesis = 7.;
44347160	186
	187	param->fMaxChi2 = 35.;
	188
	189	param->fMaxChi2s[0] = 25.; //40
	190	param->fMaxChi2s[1] = 25.; //40
	191	param->fMaxChi2s[2] = 25.; //40
	192	param->fMaxChi2s[3] = 25.; //40
	193	param->fMaxChi2s[4] = 40.; //40
	194	param->fMaxChi2s[5] = 50.; //40
	195
	196	param->fMaxRoad = 6.;
	197
	198	// not used
	199	param->fMaxChi2In = 16.;
	200
	201	param->fMaxChi2sR[0] = 10.;
	202	param->fMaxChi2sR[1] = 10.;
	203	param->fMaxChi2sR[2] = 10.;
	204	param->fMaxChi2sR[3] = 10.;
	205	param->fMaxChi2sR[4] = 30.;
	206	param->fMaxChi2sR[5] = 40.;
	207
	208	param->fChi2PerCluster = 9.;
	209	// not used
	210
	211	param->fXV = 0.;
	212	param->fYV = 0.;
	213	param->fZV = 0.;
	214	param->fSigmaXV = 0.0050;
	215	param->fSigmaYV = 0.0050;
	216	param->fSigmaZV = 0.0100;
	217
	218	param->fVertexCut = 25.;
afd25725	219
	220	param->fMaxDZforPrimTrk = 0.4;
	221	param->fMaxDZToUseConstraint = 3.;
	222
	223	param->fMaxDforV0dghtrForProlongation = 30.;
6518a6c5	224	param->fMaxDForProlongation = 40.;//10.;
6518a6c5	225	param->fMaxDZForProlongation = 60.;//20.;
afd25725	226	param->fMinPtForProlongation = 0.120;
	227
	228	param->fZWindowDeadZone = 2.0;
	229	param->fSigmaXDeadZoneHit2 = 0.004/12.;
	230	param->fSigmaZDeadZoneHit2 = 0.001/12.;
	231	param->fXPassDeadZoneHits = 0.018;
44347160	232
	233	return param;
	234	}
	235	//_____________________________________________________________________________
	236	AliITSRecoParam *AliITSRecoParam::GetLowFluxParam()
	237	{
	238	//
	239	// make default reconstruction parameters for low flux env.
	240	//
c7d6d7b7	241	AliITSRecoParam *param = new AliITSRecoParam();
	242
	243	param->fMaxSnp = 0.95;
	244
	245	param->fNSigmaYLayerForRoadY = 4.;
	246	param->fNSigmaRoadY = 7.5;
	247	param->fNSigmaZLayerForRoadZ = 4.;
	248	param->fNSigmaRoadZ = 7.5;
	249
	250	param->fNSigma2RoadZC = 60.; //7.75^2
	251	param->fNSigma2RoadYC = 60.; //7.75^2
	252	param->fNSigma2RoadZNonC = 50.; //7.07^2
	253	param->fNSigma2RoadYNonC = 50.; //7.07^2
	254
	255	param->fMaxChi2PerCluster[0] = 11.; //7
	256	param->fMaxChi2PerCluster[1] = 12.; //5
	257	param->fMaxChi2PerCluster[2] = 12.; //8
	258	param->fMaxChi2PerCluster[3] = 5.; //8
	259	param->fMaxChi2PerCluster[4] = 12.; //6.5
	260
	261	param->fMaxNormChi2NonC[0] = 7.;
	262	param->fMaxNormChi2NonC[1] = 8.;
	263	param->fMaxNormChi2NonC[2] = 8.;
	264	param->fMaxNormChi2NonC[3] = 11.;
	265	param->fMaxNormChi2NonC[4] = 14.;
	266	param->fMaxNormChi2NonC[5] = 25.;
	267
	268	param->fMaxNormChi2C[0] = 11.;
	269	param->fMaxNormChi2C[1] = 13.;
	270	param->fMaxNormChi2C[2] = 15.;
	271	param->fMaxNormChi2C[3] = 18.;
	272	param->fMaxNormChi2C[4] = 30.;
	273	param->fMaxNormChi2C[5] = 35.;
	274
	275	param->fMaxNormChi2NonCForHypothesis = 7.;
	276
	277	param->fMaxChi2 = 35.;
	278
	279	param->fMaxChi2s[0] = 25.; //40
	280	param->fMaxChi2s[1] = 25.; //40
	281	param->fMaxChi2s[2] = 25.; //40
	282	param->fMaxChi2s[3] = 25.; //40
	283	param->fMaxChi2s[4] = 40.; //40
	284	param->fMaxChi2s[5] = 50.; //40
	285
	286	param->fMaxRoad = 6.;
	287
	288	// not used
	289	param->fMaxChi2In = 16.;
	290
	291	param->fMaxChi2sR[0] = 10.;
	292	param->fMaxChi2sR[1] = 10.;
	293	param->fMaxChi2sR[2] = 10.;
	294	param->fMaxChi2sR[3] = 10.;
	295	param->fMaxChi2sR[4] = 30.;
	296	param->fMaxChi2sR[5] = 40.;
	297
	298	param->fChi2PerCluster = 9.;
	299	// not used
	300
	301	param->fXV = 0.;
	302	param->fYV = 0.;
	303	param->fZV = 0.;
	304	param->fSigmaXV = 0.0050;
305	param->fSigmaYV = 0.0050;
306	param->fSigmaZV = 0.0100;
307
308	param->fVertexCut = 25.;
309
310	param->fMaxDZforPrimTrk = 0.4;
311	param->fMaxDZToUseConstraint = 3.;
312
313	param->fMaxDforV0dghtrForProlongation = 30.;
314	param->fMaxDForProlongation = 40.;//10.;
315	param->fMaxDZForProlongation = 60.;//20.;
316	param->fMinPtForProlongation = 0.120;
317
318	param->fZWindowDeadZone = 2.0;
319	param->fSigmaXDeadZoneHit2 = 0.004/12.;
320	param->fSigmaZDeadZoneHit2 = 0.001/12.;
321	param->fXPassDeadZoneHits = 0.018;
322	param->SetNLoopsSA(10);
323
324	return param;
44347160	325	}
	326	//_____________________________________________________________________________
	327	AliITSRecoParam *AliITSRecoParam::GetCosmicTestParam()
	328	{
	329	//
	330	// make default reconstruction parameters for cosmics
	331	//
e340bb86	332	AliITSRecoParam *param = new AliITSRecoParam();
e340bb86	333
2755f080	334	// larger seach windows for SA (in case of large misalignments)
2755f080	335	param->SetFactorSAWindowSizes(3.);
e340bb86	336
	337	param->fMaxSnp = 0.95;
	338
	339	param->fNSigmaYLayerForRoadY = 4.;
	340	param->fNSigmaRoadY = 7.5;
	341	param->fNSigmaZLayerForRoadZ = 4.;
	342	param->fNSigmaRoadZ = 7.5;
	343
	344	param->fNSigma2RoadZC = 60.; //7.75^2
	345	param->fNSigma2RoadYC = 60.; //7.75^2
	346	param->fNSigma2RoadZNonC = 50.; //7.07^2
	347	param->fNSigma2RoadYNonC = 50.; //7.07^2
	348
	349	param->fMaxChi2PerCluster[0] = 11.; //7
	350	param->fMaxChi2PerCluster[1] = 12.; //5
	351	param->fMaxChi2PerCluster[2] = 12.; //8
	352	param->fMaxChi2PerCluster[3] = 5.; //8
	353	param->fMaxChi2PerCluster[4] = 12.; //6.5
	354
	355	param->fMaxNormChi2NonC[0] = 7.;
	356	param->fMaxNormChi2NonC[1] = 8.;
	357	param->fMaxNormChi2NonC[2] = 8.;
	358	param->fMaxNormChi2NonC[3] = 11.;
	359	param->fMaxNormChi2NonC[4] = 14.;
	360	param->fMaxNormChi2NonC[5] = 25.;
	361
	362	param->fMaxNormChi2C[0] = 11.;
	363	param->fMaxNormChi2C[1] = 13.;
	364	param->fMaxNormChi2C[2] = 15.;
	365	param->fMaxNormChi2C[3] = 18.;
	366	param->fMaxNormChi2C[4] = 30.;
	367	param->fMaxNormChi2C[5] = 35.;
	368
	369	param->fMaxNormChi2NonCForHypothesis = 7.;
	370
	371	param->fMaxChi2 = 35.;
	372
	373	param->fMaxChi2s[0] = 25.; //40
	374	param->fMaxChi2s[1] = 25.; //40
	375	param->fMaxChi2s[2] = 25.; //40
	376	param->fMaxChi2s[3] = 25.; //40
	377	param->fMaxChi2s[4] = 40.; //40
	378	param->fMaxChi2s[5] = 50.; //40
	379
	380	param->fMaxRoad = 6.;
	381
	382	// not used
	383	param->fMaxChi2In = 16.;
	384
	385	param->fMaxChi2sR[0] = 10.;
	386	param->fMaxChi2sR[1] = 10.;
	387	param->fMaxChi2sR[2] = 10.;
	388	param->fMaxChi2sR[3] = 10.;
	389	param->fMaxChi2sR[4] = 30.;
	390	param->fMaxChi2sR[5] = 40.;
	391
	392	param->fChi2PerCluster = 9.;
	393	// not used
	394
	395	param->fXV = 0.;
	396	param->fYV = 0.;
	397	param->fZV = 0.;
	398	param->fSigmaXV = 0.0050;
	399	param->fSigmaYV = 0.0050;
400	param->fSigmaZV = 0.0100;
401
402	param->fVertexCut = 25.;
403
404	param->fMaxDZforPrimTrk = 0.4;
405	param->fMaxDZToUseConstraint = 3.;
406
407	param->fMaxDforV0dghtrForProlongation = 30.;
408	param->fMaxDForProlongation = 10.;
409	param->fMaxDZForProlongation = 20.;
410	param->fMinPtForProlongation = 0.120;
411
412	param->fZWindowDeadZone = 2.0;
413	param->fSigmaXDeadZoneHit2 = 0.004/12.;
414	param->fSigmaZDeadZoneHit2 = 0.001/12.;
415	param->fXPassDeadZoneHits = 0.018;
416
417
418	return param;
44347160	419	}
44347160	420	//_____________________________________________________________________________
ae00569a	421	AliITSRecoParam *AliITSRecoParam::GetPlaneEffParam(Int_t i)
	422	{
	423	//
	424	// make special reconstruction parameters for Plane Efficiency study on layer i
	425	//
	426	AliITSRecoParam *param;
	427	param = GetHighFluxParam();
	428	param->SetComputePlaneEff();
	429	param->SetLayerToSkip(i);
	430	return param;
	431	}
	432	//_____________________________________________________________________________
44347160	433	void AliITSRecoParam::SetLayersParameters()
	434	{
	435	//
	436	// number of layers and layers spatial resolutions
	437	//
	438
	439	// spatial resolutions of the detectors
	440	// y: 12 12 38 38 20 20 micron
	441	fSigmaY2[0]=1.44e-6;
	442	fSigmaY2[1]=1.44e-6;
	443	fSigmaY2[2]=1.444e-5;
	444	fSigmaY2[3]=1.444e-5;
	445	fSigmaY2[4]=4.0e-6;
	446	fSigmaY2[5]=4.0e-6;
	447	// z: 120 120 28 28 830 830 micron
	448	fSigmaZ2[0]=1.44e-4;
	449	fSigmaZ2[1]=1.44e-4;
	450	fSigmaZ2[2]=7.84e-6;
	451	fSigmaZ2[3]=7.84e-6;
	452	fSigmaZ2[4]=6.889e-3;
	453	fSigmaZ2[5]=6.889e-3;
	454
	455	return;
	456	}
ed446fa3	457	//_____________________________________________________________________________
	458	void AliITSRecoParam::PrintParameters() const
	459	{
	460	//
	461	// print parameters
	462	//
	463
	464	printf("============================= AliITSRecoParam::PrintParameters ");
1c97ce2f	465	printf("============================= \n \n");
ed446fa3	466	for(Int_t i=0; i<AliITSgeomTGeo::kNLayers; i++) {
	467	if(!fLayersToSkip[i]) {
	468	printf("ITS Traking: using layer %d\n",i);
	469	} else {
	470	printf("ITS Traking: skipping layer %d\n",i);
	471	}
	472	}
	473	for(Int_t i=0; i<AliITSgeomTGeo::kNLayers; i++) {
	474	if(fUseAmplitudeInfo[i]) {
	475	printf("ITS Traking: use amplitude info for layer %d\n",i);
	476	} else {
	477	printf("ITS Traking: don't use amplitude info for layer %d\n",i);
	478	}
	479	}
	480	for(Int_t i=0; i<AliITSgeomTGeo::kNLayers; i++)
	481	printf("Layer %d:\n sigmaY2 %f, sigma Z2 %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],fMaxNormChi2NonC[i],fMaxNormChi2C[i],fMaxChi2s[i]);
	482
	483
	484	Dump();
	485
	486	return;
	487	}