|
|
|
#include <HitDigitizerFP420.h>
Public Types | |
| typedef std::map< int, float, std::less< int > > | hit_map_type |
Public Member Functions | |
| HitDigitizerFP420 (float in, float ild, float ildx, float ildy, float in0, float in2, float in3, int verbosity) | |
| hit_map_type | processHit (const PSimHit &, G4ThreeVector, int, int, double, int, double, double, int) |
| void | setChargeCollectionDrifter (CDrifterFP420 *cd) |
| void | setChargeDivider (CDividerFP420 *cd) |
| void | setInduceChargeOnElectrods (IChargeFP420 *cd) |
| ~HitDigitizerFP420 () | |
Private Member Functions | |
| G4ThreeVector | DriftDirection (G4ThreeVector, int, int) |
Private Attributes | |
| double | appliedVoltage |
| double | chargeDistributionRMS |
| double | chargeMobility |
| double | depletionVoltage |
| double | gevperelectron |
| double | moduleThickness |
| bool | noDiffusion |
| float | tanLorentzAnglePerTesla |
| double | temperature |
| CDividerFP420 * | theCDividerFP420 |
| CDrifterFP420 * | theCDrifterFP420 |
| IChargeFP420 * | theIChargeFP420 |
Definition at line 17 of file HitDigitizerFP420.h.
| typedef std::map<int, float, std::less<int> > HitDigitizerFP420::hit_map_type |
Definition at line 20 of file HitDigitizerFP420.h.
| HitDigitizerFP420::HitDigitizerFP420 | ( | float | in, |
| float | ild, | ||
| float | ildx, | ||
| float | ildy, | ||
| float | in0, | ||
| float | in2, | ||
| float | in3, | ||
| int | verbosity | ||
| ) |
Definition at line 28 of file HitDigitizerFP420.cc.
References CBOLTZ, gather_cfg::cout, e_SI, recoMuon::in, and funct::pow().
{
moduleThickness =in;
double bz420 = in0;
double bzD2 = in2;
double bzD3 = in3;
// double pitch =inp;
// double pitchX =inpx;
// double pitchY =inpy;
double ldrift =ild;
double ldriftX =ildx;
double ldriftY =ildy;
//
// Construct default classes
//
//theCDividerFP420 = new ChargeDividerFP420(pitch);
theCDividerFP420 = new ChargeDividerFP420(moduleThickness, bz420, bzD2, bzD3, verbosity);
depletionVoltage=20.0; //
appliedVoltage=25.0; // a bit bigger than depletionVoltage to have positive value A for logA, A=1-2*Tfract.*Vd/(Vd+Vb)
// chargeMobility=480.0;// = 480.0 !holes mobility [cm**2/V/sec] p-side; = 1350.0 !electron mobility - n-side
chargeMobility=1350.0;// = 480.0 !holes mobility [cm**2/V/sec] p-side; = 1350.0 !electron mobility - n-side
//temperature=297.; // 24 C degree +273 = 297 ---->diffusion const for electrons= (1.38E-23/1.6E-19)*1350.0*297=34.6
// diffusion const for holes = 12.3 [cm**2/sec]
temperature=263.; // -10 C degree +273 = 263 ---->diffusion const for electrons= (1.38E-23/1.6E-19)*1350.0*263=30.6
double diffusionConstant = CBOLTZ/e_SI * chargeMobility * temperature;
// noDiffusion=true; // true if no Diffusion
noDiffusion=false; // false if Diffusion
if (noDiffusion) diffusionConstant *= 1.0e-3;
chargeDistributionRMS=6.5e-10;
// arbitrary:
tanLorentzAnglePerTesla = 0. ; // try =0.106 if B field exist
// double timeNormalisation = pow(moduleThickness,2)/(2.*depletionVoltage*chargeMobility);
// double timeNormalisation = pow(pitch,2)/(2.*depletionVoltage*chargeMobility);
double timeNormalisation = pow(ldrift,2)/(2.*depletionVoltage*chargeMobility);//
// double timeNormalisation = pow(ldrift,2)/(2.*depletionVoltage*chargeMobility);//
timeNormalisation = timeNormalisation*0.01; // because ldrift in [mm] but mu_e in [cm2/V/sec
// double timeNormalisation = pow(pitch/2.,2)/(2.*depletionVoltage*chargeMobility);// i entered pitch as a distance between 2 read-out strips, not real distance between any (p or n) electrods which will be in 2 times less. But in expression for timeNormalisation the real distance of charge collection must be, so use pitch/2. !
double clusterWidth=5.;// was = 3
gevperelectron= 3.61e-09; // double GevPerElectron = 3.61e-09
// GevPerElectron AZ:average deposited energy per e-h pair [keV]??? =0.0036
if(verbosity>0) {
std::cout << "HitDigitizerFP420: constructor ldrift= " << ldrift << std::endl;
std::cout << "ldriftY= " <<ldriftY << "ldriftX= " <<ldriftX << std::endl;
std::cout << "depletionVoltage" << depletionVoltage << "appliedVoltage" << appliedVoltage << "chargeMobility" << chargeMobility << "temperature" << temperature << "diffusionConstant" << diffusionConstant << "chargeDistributionRMS" << chargeDistributionRMS << "moduleThickness" << moduleThickness << "timeNormalisation" << timeNormalisation << "gevperelectron" << gevperelectron << std::endl;
}
//ndif
theCDrifterFP420 =
new ChargeDrifterFP420(moduleThickness,
timeNormalisation,
diffusionConstant,
temperature,
chargeDistributionRMS,
depletionVoltage,
appliedVoltage,
ldriftX,
ldriftY, verbosity);
// pitchX,
// pitchY);
theIChargeFP420 = new InduceChargeFP420(clusterWidth,gevperelectron);
}
| HitDigitizerFP420::~HitDigitizerFP420 | ( | ) |
Definition at line 102 of file HitDigitizerFP420.cc.
{
delete theCDividerFP420;
delete theCDrifterFP420;
delete theIChargeFP420;
}
| G4ThreeVector HitDigitizerFP420::DriftDirection | ( | G4ThreeVector | _bfield, |
| int | xytype, | ||
| int | verbosity | ||
| ) | [private] |
Definition at line 192 of file HitDigitizerFP420.cc.
References gather_cfg::cout.
{
// LOCAL hit: exchange xytype: 1 <-> 2
// Frame detFrame(_detp->surface().position(),_detp->surface().rotation());
// G4ThreeVector Bfield=detFrame.toLocal(_bfield);
G4ThreeVector Bfield=_bfield;
float dir_x, dir_y, dir_z;
// this lines with dir_... have to be specified(sign?) in dependence of field direction:
/*
float dir_x = tanLorentzAnglePerTesla * Bfield.y();
float dir_y = -tanLorentzAnglePerTesla * Bfield.x();
float dir_z = 1.; // if E field is in z direction
*/
// for electrons:
// E field is either in X or Y direction with change vector to opposite
// global Y or localX
// E field is in Xlocal direction with change vector to opposite
if(xytype == 2) {
dir_x = 1.; // E field in Xlocal direction
dir_y = +tanLorentzAnglePerTesla * Bfield.z();
dir_z = -tanLorentzAnglePerTesla * Bfield.y();
}
// global X
// E field is in Ylocal direction with change vector to opposite
else if(xytype == 1) {
dir_x = +tanLorentzAnglePerTesla * Bfield.z();
dir_y = 1.; // E field in Ylocal direction
dir_z = -tanLorentzAnglePerTesla * Bfield.x();
}
else{
dir_x = 0.;
dir_y = 0.;
dir_z = 0.;
std::cout << "HitDigitizerFP420: ERROR - wrong xytype=" << xytype << std::endl;
}
// G4ThreeVector theDriftDirection = LocalVector(dir_x,dir_y,dir_z);
G4ThreeVector theDriftDirection(dir_x,dir_y,dir_z);
// Local3DPoint EntryPo(aHit->getEntry().x(),aHit->getEntry().y(),aHit->getEntry().z());
if(verbosity>0) {
std::cout << "HitDigitizerFP420:DriftDirection tanLorentzAnglePerTesla= " << tanLorentzAnglePerTesla << std::endl;
std::cout << "HitDigitizerFP420:DriftDirection The drift direction in local coordinate is " << theDriftDirection << std::endl;
}
return theDriftDirection;
}
| HitDigitizerFP420::hit_map_type HitDigitizerFP420::processHit | ( | const PSimHit & | hit, |
| G4ThreeVector | bfield, | ||
| int | xytype, | ||
| int | numStrips, | ||
| double | pitch, | ||
| int | numStripsW, | ||
| double | pitchW, | ||
| double | moduleThickness, | ||
| int | verbosity | ||
| ) |
Definition at line 110 of file HitDigitizerFP420.cc.
References gather_cfg::cout, PSimHit::entryPoint(), PSimHit::exitPoint(), PFRecoTauDiscriminationAgainstMuon2_cfi::verbosity, PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().
Referenced by FP420DigiMain::run().
{
// use chargePosition just for cross-check in "induce" method
// hit center in 3D-detector r.f.
float middlex = (hit.exitPoint().x() + hit.entryPoint().x() )/2.;
float middley = (hit.exitPoint().y() + hit.entryPoint().y() )/2.;
float chargePosition= -100.;
// Y:
if(xytype == 1) {
// chargePosition = fabs(-numStrips/2. - ( int(middle.x()/pitch) +1.) );
//chargePosition = fabs(int(middle.x()/pitch+0.5*(numStrips+1)) + 1.);
// chargePosition = fabs(int(middle.y()/pitch+0.5*(numStrips+1)) + 1.);
// local and global reference frames are rotated in 90 degree, so global X and local Y are collinear
// chargePosition = int(fabs(middle.x()/pitch + 0.5*numStrips + 1.));// charge in strip coord
chargePosition = 0.5*(numStrips) + middlex/pitch ;// charge in strip coord 0 - numStrips-1
}
// X:
else if(xytype == 2) {
// chargePosition = fabs(-numStrips/2. - ( int(middle.y()/pitch) +1.) );
//chargePosition = fabs(int(middle.y()/pitch+0.5*(numStrips+1)) + 1.);
// chargePosition = fabs(int(middle.x()/pitch+0.5*(numStrips+1)) + 1.);
// local and global reference frames are rotated in 90 degree, so global X and local Y are collinear
// chargePosition = int(fabs(middle.y()/pitch + 0.5*numStrips + 1.));
chargePosition = 0.5*(numStrips) + middley/pitch ;// charge in strip coord 0 - numStrips-1
// std::cout << " chargePosition SiHitDi... = " << chargePosition << std::endl;
}
else {
std::cout <<"================================================================"<<std::endl;
std::cout << "**** HitDigitizerFP420: !!! ERROR: you have not to be here !!! xytype=" << xytype << std::endl;
// std::cout << "**** HitDigitizerFP420: !!! ERROR: you have not to be here !!! xytype=" << xytype << std::endl;
// break;
}
// if(chargePosition > numStrips || chargePosition<1) {
if(chargePosition > numStrips || chargePosition < 0) {
std::cout << "**** HitDigitizerFP420: !!! ERROR: check correspondence of XY detector dimensions in XML and here !!! chargePosition = " << chargePosition << std::endl;
// break;
}
if(verbosity>0) {
std::cout << " ====== **** HitDigitizerFP420: !!! processHit !!! : input: xytype=" << xytype << " numStrips= " << numStrips << " pitch= " << pitch << " Calculated chargePosition= " << chargePosition << std::endl;
std::cout << "The middle of hit point on input was: middlex = " << middlex << std::endl;
std::cout << "The middle of hit point on input was: middley = " << middley << std::endl;
// std::cout << "For checks: hit point Entry = " << hit.getEntry() << std::endl;
std::cout << " ====== **** HitDigitizerFP420:processHit: start CDrifterFP420 divide" << std::endl;
}
//
// Fully process one SimHit
//
// CDrifterFP420::ionization_type ion = theCDividerFP420->divide(hit,pitch);
CDrifterFP420::ionization_type ion = theCDividerFP420->divide(hit,moduleThickness);
//
// Compute the drift direction for this det
//
// G4ThreeVector driftDir = DriftDirection(&det,bfield);
G4ThreeVector driftDir = DriftDirection(bfield,xytype,verbosity);
if(verbosity>0) {
std::cout << " ====== **** HitDigitizerFP420:processHit: driftDir= " << driftDir << std::endl;
std::cout << " ====== **** HitDigitizerFP420:processHit: start induce , CDrifterFP420 drift " << std::endl;
}
// if(driftDir.z() ==0.) {
// std::cout << " pxlx: drift in z is zero " << std::endl;
// } else
//
return theIChargeFP420->induce(theCDrifterFP420->drift(ion,driftDir,xytype), numStrips, pitch, numStripsW, pitchW, xytype, verbosity);
//
}
| void HitDigitizerFP420::setChargeCollectionDrifter | ( | CDrifterFP420 * | cd | ) | [inline] |
Definition at line 33 of file HitDigitizerFP420.h.
References theCDrifterFP420.
{
if (theCDrifterFP420) delete theCDrifterFP420;
theCDrifterFP420 = cd;
}
| void HitDigitizerFP420::setChargeDivider | ( | CDividerFP420 * | cd | ) | [inline] |
Definition at line 29 of file HitDigitizerFP420.h.
References theCDividerFP420.
{
if (theCDividerFP420) delete theCDividerFP420;
theCDividerFP420 = cd;
}
| void HitDigitizerFP420::setInduceChargeOnElectrods | ( | IChargeFP420 * | cd | ) | [inline] |
Definition at line 37 of file HitDigitizerFP420.h.
References theIChargeFP420.
{
if (theIChargeFP420) delete theIChargeFP420;
theIChargeFP420 = cd;
}
double HitDigitizerFP420::appliedVoltage [private] |
Definition at line 59 of file HitDigitizerFP420.h.
double HitDigitizerFP420::chargeDistributionRMS [private] |
Definition at line 63 of file HitDigitizerFP420.h.
double HitDigitizerFP420::chargeMobility [private] |
Definition at line 60 of file HitDigitizerFP420.h.
double HitDigitizerFP420::depletionVoltage [private] |
Definition at line 58 of file HitDigitizerFP420.h.
double HitDigitizerFP420::gevperelectron [private] |
Definition at line 64 of file HitDigitizerFP420.h.
double HitDigitizerFP420::moduleThickness [private] |
Definition at line 54 of file HitDigitizerFP420.h.
bool HitDigitizerFP420::noDiffusion [private] |
Definition at line 62 of file HitDigitizerFP420.h.
float HitDigitizerFP420::tanLorentzAnglePerTesla [private] |
Definition at line 70 of file HitDigitizerFP420.h.
double HitDigitizerFP420::temperature [private] |
Definition at line 61 of file HitDigitizerFP420.h.
CDividerFP420* HitDigitizerFP420::theCDividerFP420 [private] |
Definition at line 50 of file HitDigitizerFP420.h.
Referenced by setChargeDivider().
CDrifterFP420* HitDigitizerFP420::theCDrifterFP420 [private] |
Definition at line 51 of file HitDigitizerFP420.h.
Referenced by setChargeCollectionDrifter().
IChargeFP420* HitDigitizerFP420::theIChargeFP420 [private] |
Definition at line 52 of file HitDigitizerFP420.h.
Referenced by setInduceChargeOnElectrods().
1.7.3