[u/mrichter/AliRoot.git] / ITS / AliITSSimuParam.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.                  *
 **************************************************************************/

/* $Id:$ */

///////////////////////////////////////////////////////////////////
//                                                               //
// Implementation of the class to store the parameters used in   //
// the simulation of SPD, SDD and SSD detectors                  //
// Origin: F.Prino, Torino, prino@to.infn.it                     //
//                                                               //
///////////////////////////////////////////////////////////////////

#include "AliITSSimuParam.h"
#include <TMath.h>

const Float_t  AliITSSimuParam::fgkSPDBiasVoltageDefault = 18.182;
const Double_t AliITSSimuParam::fgkSPDThreshDefault = 3000.;
const Double_t AliITSSimuParam::fgkSPDSigmaDefault = 250.;
const TString  AliITSSimuParam::fgkSPDCouplingOptDefault = "old";
const Double_t AliITSSimuParam::fgkSPDCouplColDefault = 0.;
const Double_t AliITSSimuParam::fgkSPDCouplRowDefault = 0.055;
const Float_t  AliITSSimuParam::fgkSPDEccDiffDefault = 0.85;
const Float_t  AliITSSimuParam::fgkSDDDiffCoeffDefault = 3.23;
const Float_t  AliITSSimuParam::fgkSDDDiffCoeff1Default = 30.;
const Float_t  AliITSSimuParam::fgkSDDJitterErrorDefault = 20.; // 20 um from beam test 2001
const Float_t  AliITSSimuParam::fgkSDDDynamicRangeDefault = 132.;
const Int_t    AliITSSimuParam::fgkSDDMaxAdcDefault = 1024;
const Float_t  AliITSSimuParam::fgkSDDChargeLossDefault = 0.;
const Double_t AliITSSimuParam::fgkSSDCouplingPRDefault = 0.01;
const Double_t AliITSSimuParam::fgkSSDCouplingPLDefault = 0.01;
const Double_t AliITSSimuParam::fgkSSDCouplingNRDefault = 0.01;
const Double_t AliITSSimuParam::fgkSSDCouplingNLDefault = 0.01;
const Int_t    AliITSSimuParam::fgkSSDZSThresholdDefault = 3;

const Float_t AliITSSimuParam::fgkNsigmasDefault = 3.;
const Int_t AliITSSimuParam::fgkNcompsDefault = 121;

ClassImp(AliITSSimuParam)

//______________________________________________________________________
AliITSSimuParam::AliITSSimuParam():
  TObject(),
fGeVcharge(0.),
fDOverV(0.),
fSPDBiasVoltage(fgkSPDBiasVoltageDefault),
fSPDThresh(fgkSPDThreshDefault),
fSPDSigma(fgkSPDSigmaDefault),
fSPDCouplOpt(0),
fSPDCouplCol(fgkSPDCouplColDefault),
fSPDCouplRow(fgkSPDCouplRowDefault),
fSPDEccDiff(0.),
fSDDElectronics(0),
fSDDDiffCoeff(0.),
fSDDDiffCoeff1(0.),
fSDDJitterError(fgkSDDJitterErrorDefault),
fSDDDynamicRange(fgkSDDDynamicRangeDefault),
fSDDMaxAdc(0.),
fSDDChargeLoss(fgkSDDChargeLossDefault),
fSSDADCpereV(0.),
fSSDCouplingPR(0),
fSSDCouplingPL(0),
fSSDCouplingNR(0),
fSSDCouplingNL(0),
fSSDZSThreshold(fgkSSDZSThresholdDefault),
fNsigmas(fgkNsigmasDefault),
fNcomps(fgkNcompsDefault),
fGaus()
{  
  // default constructor
  SetGeVToCharge();
  SetDistanceOverVoltage();
  SetSPDCouplingOption(fgkSPDCouplingOptDefault);
  SetSPDSigmaDiffusionAsymmetry(fgkSPDEccDiffDefault);
  SetSDDElectronics();
  SetSDDDiffCoeff(fgkSDDDiffCoeffDefault,fgkSDDDiffCoeff1Default);
  SetSDDMaxAdc((Double_t)fgkSDDMaxAdcDefault);
  SetSSDADCpereV();
  SetSSDCouplings(fgkSSDCouplingPRDefault,fgkSSDCouplingPLDefault,fgkSSDCouplingNRDefault,fgkSSDCouplingNLDefault);
}
//______________________________________________________________________
AliITSSimuParam::AliITSSimuParam(const AliITSSimuParam &simpar):
TObject(),
fGeVcharge(simpar.fGeVcharge),
fDOverV(simpar.fDOverV),
fSPDBiasVoltage(simpar.fSPDBiasVoltage),
fSPDThresh(simpar.fSPDThresh),
fSPDSigma(simpar.fSPDSigma),
fSPDCouplOpt(simpar.fSPDCouplOpt),
fSPDCouplCol(simpar.fSPDCouplCol),
fSPDCouplRow(simpar.fSPDCouplRow),
fSPDEccDiff(simpar.fSPDEccDiff),
fSDDElectronics(simpar.fSDDElectronics),
fSDDDiffCoeff(simpar.fSDDDiffCoeff),
fSDDDiffCoeff1(simpar.fSDDDiffCoeff1),
fSDDJitterError(simpar.fSDDJitterError),
fSDDDynamicRange(simpar.fSDDDynamicRange),
fSDDMaxAdc(simpar.fSDDMaxAdc),
fSDDChargeLoss(simpar.fSDDChargeLoss),
fSSDADCpereV(simpar.fSSDADCpereV),
fSSDCouplingPR(simpar.fSSDCouplingPR),
fSSDCouplingPL(simpar.fSSDCouplingPL),
fSSDCouplingNR(simpar.fSSDCouplingNR),
fSSDCouplingNL(simpar.fSSDCouplingNL),
fSSDZSThreshold(simpar.fSSDZSThreshold),
fNsigmas(simpar.fNsigmas),
fNcomps(simpar.fNcomps),
fGaus(){
  // copy constructor
}

//______________________________________________________________________
AliITSSimuParam& AliITSSimuParam::operator=(const AliITSSimuParam& source){
    // Assignment operator. 
    this->~AliITSSimuParam();
    new(this) AliITSSimuParam(source);
    return *this;

}


//______________________________________________________________________
AliITSSimuParam::~AliITSSimuParam() {
  // destructor
  if(fGaus) delete fGaus;
}
//________________________________________________________________________
void AliITSSimuParam::SetNLookUp(Int_t p1){
  // Set number of sigmas over which cluster disintegration is performed
  fNcomps=p1;
  if (fGaus) delete fGaus;
  fGaus = new TArrayF(fNcomps+1);
  for(Int_t i=0; i<=fNcomps; i++) {
    Float_t x = -fNsigmas + (2.*i*fNsigmas)/(fNcomps-1);
    (*fGaus)[i] = exp(-((x*x)/2));
  }
}
//________________________________________________________________________
void AliITSSimuParam::PrintParameters() const{
  printf("GeVToCharge               = %G\n",fGeVcharge);
  printf("DistanveOverVoltage       = %f \n",fDOverV);
  printf("\n");
  printf("=====  SPD parameters  =====\n");
  printf("Bias Voltage              = %f \n",fSPDBiasVoltage);
  printf("Threshold and sigma       = %f %f\n",fSPDThresh,fSPDSigma);
  printf("Coupling Option           = %s\n",fSPDCouplOpt.Data());
  printf("Coupling value (column)   = %f\n",fSPDCouplCol);
  printf("Coupling value (row)      = %f\n",fSPDCouplRow);
  printf("Eccentricity in diffusion = %f\n",fSPDEccDiff);
  printf("\n");
  printf("=====  SDD parameters  =====\n");
  printf("Electronic chips          = %d\n",fSDDElectronics);
  printf("Diffusion Coefficients    = %f %f\n",fSDDDiffCoeff,fSDDDiffCoeff1);
  printf("Jitter Error              = %f um\n",fSDDJitterError);
  printf("Dynamic Range             = %f\n",fSDDDynamicRange);
  printf("Max. ADC                  = %f\n",fSDDMaxAdc);
  printf("Charge Loss               = %f\n",fSDDChargeLoss);  
  printf("\n");
  printf("=====  SSD parameters  =====\n");
  printf("ADC per eV                = %f\n",fSSDADCpereV);
  printf("Coupling PR               = %f\n",fSSDCouplingPR);
  printf("Coupling PL               = %f\n",fSSDCouplingPL);
  printf("Coupling NR               = %f\n",fSSDCouplingNR);
  printf("Coupling NL               = %f\n",fSSDCouplingNL);
  printf("Zero Supp threshold       = %d\n",fSSDZSThreshold);
}
Commit	Line	Data
cd2a0045	1	/**************************************************************************
	2	* Copyright(c) 2007-2009, ALICE Experiment at CERN, All rights reserved. *
	3	* *
	4	* Author: The ALICE Off-line Project. *
	5	* Contributors are mentioned in the code where appropriate. *
	6	* *
	7	* Permission to use, copy, modify and distribute this software and its *
	8	* documentation strictly for non-commercial purposes is hereby granted *
	9	* without fee, provided that the above copyright notice appears in all *
	10	* copies and that both the copyright notice and this permission notice *
	11	* appear in the supporting documentation. The authors make no claims *
	12	* about the suitability of this software for any purpose. It is *
	13	* provided "as is" without express or implied warranty. *
	14	**************************************************************************/
	15
	16	/* $Id:$ */
	17
	18	///////////////////////////////////////////////////////////////////
	19	// //
	20	// Implementation of the class to store the parameters used in //
	21	// the simulation of SPD, SDD and SSD detectors //
	22	// Origin: F.Prino, Torino, prino@to.infn.it //
	23	// //
	24	///////////////////////////////////////////////////////////////////
	25
	26	#include "AliITSSimuParam.h"
	27	#include <TMath.h>
	28
	29	const Float_t AliITSSimuParam::fgkSPDBiasVoltageDefault = 18.182;
	30	const Double_t AliITSSimuParam::fgkSPDThreshDefault = 3000.;
	31	const Double_t AliITSSimuParam::fgkSPDSigmaDefault = 250.;
	32	const TString AliITSSimuParam::fgkSPDCouplingOptDefault = "old";
	33	const Double_t AliITSSimuParam::fgkSPDCouplColDefault = 0.;
	34	const Double_t AliITSSimuParam::fgkSPDCouplRowDefault = 0.055;
	35	const Float_t AliITSSimuParam::fgkSPDEccDiffDefault = 0.85;
	36	const Float_t AliITSSimuParam::fgkSDDDiffCoeffDefault = 3.23;
	37	const Float_t AliITSSimuParam::fgkSDDDiffCoeff1Default = 30.;
	38	const Float_t AliITSSimuParam::fgkSDDJitterErrorDefault = 20.; // 20 um from beam test 2001
	39	const Float_t AliITSSimuParam::fgkSDDDynamicRangeDefault = 132.;
	40	const Int_t AliITSSimuParam::fgkSDDMaxAdcDefault = 1024;
	41	const Float_t AliITSSimuParam::fgkSDDChargeLossDefault = 0.;
	42	const Double_t AliITSSimuParam::fgkSSDCouplingPRDefault = 0.01;
	43	const Double_t AliITSSimuParam::fgkSSDCouplingPLDefault = 0.01;
	44	const Double_t AliITSSimuParam::fgkSSDCouplingNRDefault = 0.01;
	45	const Double_t AliITSSimuParam::fgkSSDCouplingNLDefault = 0.01;
	46	const Int_t AliITSSimuParam::fgkSSDZSThresholdDefault = 3;
	47
	48	const Float_t AliITSSimuParam::fgkNsigmasDefault = 3.;
	49	const Int_t AliITSSimuParam::fgkNcompsDefault = 121;
	50
	51	ClassImp(AliITSSimuParam)
	52
	53	//______________________________________________________________________
	54	AliITSSimuParam::AliITSSimuParam():
	55	TObject(),
	56	fGeVcharge(0.),
	57	fDOverV(0.),
	58	fSPDBiasVoltage(fgkSPDBiasVoltageDefault),
	59	fSPDThresh(fgkSPDThreshDefault),
	60	fSPDSigma(fgkSPDSigmaDefault),
	61	fSPDCouplOpt(0),
	62	fSPDCouplCol(fgkSPDCouplColDefault),
	63	fSPDCouplRow(fgkSPDCouplRowDefault),
	64	fSPDEccDiff(0.),
65	fSDDElectronics(0),
66	fSDDDiffCoeff(0.),
67	fSDDDiffCoeff1(0.),
68	fSDDJitterError(fgkSDDJitterErrorDefault),
69	fSDDDynamicRange(fgkSDDDynamicRangeDefault),
70	fSDDMaxAdc(0.),
71	fSDDChargeLoss(fgkSDDChargeLossDefault),
72	fSSDADCpereV(0.),
73	fSSDCouplingPR(0),
74	fSSDCouplingPL(0),
75	fSSDCouplingNR(0),
76	fSSDCouplingNL(0),
77	fSSDZSThreshold(fgkSSDZSThresholdDefault),
78	fNsigmas(fgkNsigmasDefault),
79	fNcomps(fgkNcompsDefault),
80	fGaus()
81	{
82	// default constructor
83	SetGeVToCharge();
84	SetDistanceOverVoltage();
85	SetSPDCouplingOption(fgkSPDCouplingOptDefault);
86	SetSPDSigmaDiffusionAsymmetry(fgkSPDEccDiffDefault);
87	SetSDDElectronics();
88	SetSDDDiffCoeff(fgkSDDDiffCoeffDefault,fgkSDDDiffCoeff1Default);
89	SetSDDMaxAdc((Double_t)fgkSDDMaxAdcDefault);
90	SetSSDADCpereV();
91	SetSSDCouplings(fgkSSDCouplingPRDefault,fgkSSDCouplingPLDefault,fgkSSDCouplingNRDefault,fgkSSDCouplingNLDefault);
92	}
93	//______________________________________________________________________
94	AliITSSimuParam::AliITSSimuParam(const AliITSSimuParam &simpar):
95	TObject(),
96	fGeVcharge(simpar.fGeVcharge),
97	fDOverV(simpar.fDOverV),
98	fSPDBiasVoltage(simpar.fSPDBiasVoltage),
99	fSPDThresh(simpar.fSPDThresh),
100	fSPDSigma(simpar.fSPDSigma),
101	fSPDCouplOpt(simpar.fSPDCouplOpt),
102	fSPDCouplCol(simpar.fSPDCouplCol),
103	fSPDCouplRow(simpar.fSPDCouplRow),
104	fSPDEccDiff(simpar.fSPDEccDiff),
105	fSDDElectronics(simpar.fSDDElectronics),
106	fSDDDiffCoeff(simpar.fSDDDiffCoeff),
107	fSDDDiffCoeff1(simpar.fSDDDiffCoeff1),
108	fSDDJitterError(simpar.fSDDJitterError),
109	fSDDDynamicRange(simpar.fSDDDynamicRange),
110	fSDDMaxAdc(simpar.fSDDMaxAdc),
111	fSDDChargeLoss(simpar.fSDDChargeLoss),
112	fSSDADCpereV(simpar.fSSDADCpereV),
113	fSSDCouplingPR(simpar.fSSDCouplingPR),
114	fSSDCouplingPL(simpar.fSSDCouplingPL),
115	fSSDCouplingNR(simpar.fSSDCouplingNR),
116	fSSDCouplingNL(simpar.fSSDCouplingNL),
117	fSSDZSThreshold(simpar.fSSDZSThreshold),
118	fNsigmas(simpar.fNsigmas),
119	fNcomps(simpar.fNcomps),
120	fGaus(){
121	// copy constructor
122	}
123
124	//______________________________________________________________________
125	AliITSSimuParam& AliITSSimuParam::operator=(const AliITSSimuParam& source){
126	// Assignment operator.
127	this->~AliITSSimuParam();
128	new(this) AliITSSimuParam(source);
129	return *this;
130
131	}
132
133
134	//______________________________________________________________________
135	AliITSSimuParam::~AliITSSimuParam() {
136	// destructor
137	if(fGaus) delete fGaus;
138	}
139	//________________________________________________________________________
140	void AliITSSimuParam::SetNLookUp(Int_t p1){
141	// Set number of sigmas over which cluster disintegration is performed
142	fNcomps=p1;
143	if (fGaus) delete fGaus;
144	fGaus = new TArrayF(fNcomps+1);
145	for(Int_t i=0; i<=fNcomps; i++) {
146	Float_t x = -fNsigmas + (2.ifNsigmas)/(fNcomps-1);
147	(fGaus)[i] = exp(-((xx)/2));
148	}
149	}
150	//________________________________________________________________________
151	void AliITSSimuParam::PrintParameters() const{
152	printf("GeVToCharge = %G\n",fGeVcharge);
153	printf("DistanveOverVoltage = %f \n",fDOverV);
154	printf("\n");
155	printf("===== SPD parameters =====\n");
156	printf("Bias Voltage = %f \n",fSPDBiasVoltage);
157	printf("Threshold and sigma = %f %f\n",fSPDThresh,fSPDSigma);
158	printf("Coupling Option = %s\n",fSPDCouplOpt.Data());
159	printf("Coupling value (column) = %f\n",fSPDCouplCol);
160	printf("Coupling value (row) = %f\n",fSPDCouplRow);
161	printf("Eccentricity in diffusion = %f\n",fSPDEccDiff);
162	printf("\n");
163	printf("===== SDD parameters =====\n");
164	printf("Electronic chips = %d\n",fSDDElectronics);
165	printf("Diffusion Coefficients = %f %f\n",fSDDDiffCoeff,fSDDDiffCoeff1);
166	printf("Jitter Error = %f um\n",fSDDJitterError);
167	printf("Dynamic Range = %f\n",fSDDDynamicRange);
168	printf("Max. ADC = %f\n",fSDDMaxAdc);
169	printf("Charge Loss = %f\n",fSDDChargeLoss);
170	printf("\n");
171	printf("===== SSD parameters =====\n");
172	printf("ADC per eV = %f\n",fSSDADCpereV);
173	printf("Coupling PR = %f\n",fSSDCouplingPR);
174	printf("Coupling PL = %f\n",fSSDCouplingPL);
175	printf("Coupling NR = %f\n",fSSDCouplingNR);
176	printf("Coupling NL = %f\n",fSSDCouplingNL);
177	printf("Zero Supp threshold = %d\n",fSSDZSThreshold);
178	}