ChargeDrifterFP420 Class Reference

#include <ChargeDrifterFP420.h>

Inheritance diagram for ChargeDrifterFP420:

Public Member Functions
	ChargeDrifterFP420 (double, double, double, double, double, double, double, double, double, int)
CDrifterFP420::collection_type	drift (const CDrifterFP420::ionization_type &, const G4ThreeVector &, const int &)
Public Member Functions inherited from CDrifterFP420
virtual	~CDrifterFP420 ()

Private Member Functions
AmplitudeSegmentFP420	drift (const EnergySegmentFP420 &, const G4ThreeVector &, const int &)

Private Attributes
double	appV
double	constDe
double	constTe
double	depV
double	ldriftcurrX
double	ldriftcurrY
double	modulePath
double	startT0
double	temperature
int	verbo

Additional Inherited Members
Public Types inherited from CDrifterFP420
typedef std::vector < AmplitudeSegmentFP420 >	collection_type
typedef std::vector < EnergySegmentFP420 >	ionization_type

Detailed Description

Definition at line 7 of file ChargeDrifterFP420.h.

Constructor & Destructor Documentation

ChargeDrifterFP420::ChargeDrifterFP420	(	double	mt,
		double	tn,
		double	dc,
		double	tm,
		double	cdr,
		double	dv,
		double	av,
		double	ptx,
		double	pty,
		int	verbosity
	)

Definition at line 10 of file ChargeDrifterFP420.cc.

References gather_cfg::cout, and PFRecoTauDiscriminationAgainstMuon2_cfi::verbosity.

                                                                {
  //
  //
  verbo=verbosity; 
  modulePath = mt;
  constTe = tn;
  constDe = dc;
  temperature = tm;// keep just in case
  startT0 = cdr;
  depV = dv;
  appV = av;
  
  ldriftcurrX = ptx;
  ldriftcurrY = pty;
  //
  //
  
  
  
  //      edm::LogInfo("ChargeDrifterFP420") << "call constructor";
  if(verbo>0) {
    std::cout << "ChargeDrifterFP420: call constructor" << std::endl;
    std::cout << "ldriftcurrX= " << ldriftcurrX << "ldriftcurrY= " << ldriftcurrY << std::endl;
    std::cout << "modulePath= " << modulePath << "constTe= " << constTe << std::endl;
    std::cout << "ChargeDrifterFP420:----------------------" << std::endl;
  }
}

Member Function Documentation

CDrifterFP420::collection_type ChargeDrifterFP420::drift ( const CDrifterFP420::ionization_type &  ion, const G4ThreeVector &  driftDir, const int &  xytype  )

virtual

Implements CDrifterFP420.

Definition at line 48 of file ChargeDrifterFP420.cc.

References gather_cfg::cout, shallow::drift(), and i.

                                                                                                                                              {
  //
  //
  if(verbo>0) {
    std::cout << "ChargeDrifterFP420: collection_type call drift" << std::endl;
  }
  //
  //
  collection_type _temp;
  _temp.resize(ion.size());
  
  for (unsigned int i=0; i<ion.size(); i++){
    _temp[i] = drift(ion[i], driftDir, xytype);
  }
  return _temp;
}

AmplitudeSegmentFP420 ChargeDrifterFP420::drift ( const EnergySegmentFP420 &  edu, const G4ThreeVector &  drift, const int &  xytype  )

private

Definition at line 66 of file ChargeDrifterFP420.cc.

References gather_cfg::cout, relval_parameters_module::energy, log, mathSSE::sqrt(), x, detailsBasic3DVector::y, and detailsBasic3DVector::z.

                                                                              {
  //
  //
  //    sliX sliY sliZ are LOCALl(!) coordinates coming from ...hit.getEntryLocalP()...  
  //
  //
  //                Exchange xytype 1 <-> 2
  double sliX = (edu).x();
  double sliY = (edu).y();
  double sliZ = (edu).z();
  
  double currX = sliX;
  double currY = sliY;
  double currZ = sliZ;
  
  double pathFraction=0., pathValue=0.;
  double tanShiftAngleX=0.,tanShiftAngleY=0.,tanShiftAngleZ=0.;
  double xDriftDueToField=0.,yDriftDueToField=0.,zDriftDueToField=0.;
  if(verbo>0) {
    std::cout << "=================================================================== " << std::endl;
    std::cout << "ChargeDrifterFP420: xytype= " << xytype << std::endl;
    std::cout << "constTe= " <<  constTe << "ldriftcurrX= " <<  ldriftcurrX << "ldriftcurrY= " <<  ldriftcurrY << std::endl;
    std::cout << "1currX = " << currX  << "  1currY = " << currY  << "  1currZ = " << currZ  << std::endl;
    std::cout << "drift.x() = " << drift.x()  << "  drift.y() = " << drift.y()  << "  drift.z() = " << drift.z()  << std::endl;
  }
  
  // Yglobal Xlocal:
  //
  //
  // will define drift time of electrons along Xlocal or Ygobal with ldriftcurrY from Yglobal
  // change of ldriftcurrX to ldriftcurrY is done in intialization of ChargeDrifterFP420, so use old names for ldriftcurres
  if(xytype == 2) {
    tanShiftAngleY = drift.y()/drift.x();
    tanShiftAngleZ = drift.z()/drift.x();
    
    //    pathValue = fabs(sliX-ldriftcurrX*int(sliX/ldriftcurrX));
    pathValue = fabs(sliX)-int(fabs(sliX/(2*ldriftcurrX))+0.5)*(2*ldriftcurrX);
    pathValue = fabs(pathValue);
    pathFraction = pathValue/ldriftcurrX; 
    if(verbo>0) {
      std::cout << "==================================" << std::endl;
      std::cout << "fabs(sliX)= " << fabs(sliX) << "ldriftcurrX= " << ldriftcurrX << std::endl;
      std::cout << "fabs(sliX/(2*ldriftcurrX))+0.5= " << fabs(sliX/(2*ldriftcurrX))+0.5 << std::endl;
      std::cout << "int(fabs(sliX/(2*ldriftcurrX))+0.5)*(2*ldriftcurrX)= " << int(fabs(sliX/(2*ldriftcurrX))+0.5)*(2*ldriftcurrX) << std::endl;
      std::cout << "pathValue= " << pathValue << std::endl;
      std::cout << "pathFraction= " << pathFraction << std::endl;
      std::cout << "==================================" << std::endl;
    }
    
    pathFraction = pathFraction>0. ? pathFraction : 0. ;
    pathFraction = pathFraction<1. ? pathFraction : 1. ;
    
    yDriftDueToField // Drift along Y due to BField
      = pathValue*tanShiftAngleY;
    zDriftDueToField // Drift along Z due to BField
      = pathValue*tanShiftAngleZ;
    // will define Y and Z ccordinates (E along X)
    currY = sliY + yDriftDueToField;
    currZ = sliZ + zDriftDueToField;
  }
  
  // Xglobal Ylocal:
  // will define drift time of electrons along Ylocal
  else if(xytype == 1) {
    tanShiftAngleX = drift.x()/drift.y();
    tanShiftAngleZ = drift.z()/drift.y();
    
    //pathValue = fabs(sliY-ldriftcurrY*int(sliY/ldriftcurrY));
    pathValue = fabs(sliY)-int(fabs(sliY/(2*ldriftcurrY))+0.5)*(2*ldriftcurrY);
    pathValue = fabs(pathValue);
    pathFraction = pathValue/ldriftcurrY; 
    //
    //
    
    if(verbo>0) {
      std::cout << "==================================" << std::endl;
      std::cout << "fabs(sliY)= " << fabs(sliY) << "ldriftcurrY= " << ldriftcurrY << std::endl;
      std::cout << "fabs(sliY/(2*ldriftcurrY))+0.5= " << fabs(sliY/(2*ldriftcurrY))+0.5 << std::endl;
      std::cout << "int(fabs(sliY/(2*ldriftcurrY))+0.5)*(2*ldriftcurrY)= " << int(fabs(sliY/(2*ldriftcurrY))+0.5)*(2*ldriftcurrY) << std::endl;
      std::cout << "pathValue= " << pathValue << std::endl;
      std::cout << "pathFraction= " << pathFraction << std::endl;
      std::cout << "==================================" << std::endl;
    }
    
    if(pathFraction<0. || pathFraction>1.) std::cout << "ChargeDrifterFP420: ERROR:pathFraction=" << pathFraction << std::endl;
    pathFraction = pathFraction>0. ? pathFraction : 0. ;
    pathFraction = pathFraction<1. ? pathFraction : 1. ;
    xDriftDueToField // Drift along X due to BField
      = pathValue*tanShiftAngleX;
    //      = (ldriftcurrY-sliY)*tanShiftAngleX;
    zDriftDueToField // Drift along Z due to BField
      = pathValue*tanShiftAngleZ;
    // will define X and Z ccordinates (E along Y)
    currX = sliX + xDriftDueToField;
    currZ = sliZ + zDriftDueToField;
  }
  //  double tanShiftAngleX = drift.x()/drift.z();
  //  double tanShiftAngleY = drift.y()/drift.z();
  // double pathFraction = (modulePath/2.-sliZ)/modulePath ; 
  // pathFraction = pathFraction>0. ? pathFraction : 0. ;
  // pathFraction = pathFraction<1. ? pathFraction : 1. ;
  
  //uble xDriftDueToField // Drift along X due to BField
  //= (modulePath/2. - sliZ)*tanShiftAngleX;
  //uble yDriftDueToField // Drift along Y due to BField
  //= (modulePath/2. - sliZ)*tanShiftAngleY;
  //uble currX = sliX + xDriftDueToField;
  //uble currY = sliY + yDriftDueToField;
  //  
  //  
  // log is a ln 
  //  std::cout << "ChargeDrifterFP420: depV=" <<depV  << "  appV=" << appV << " startT0 =" <<startT0  << " 1.-2*depV*pathFraction/(depV+appV) =" <<1.-2*depV*pathFraction/(depV+appV)  << " log(1.-2*depV*pathFraction/(depV+appV)) =" <<log(1.-2*depV*pathFraction/(depV+appV))  << std::endl;
  double bbb = 1.-2*depV*pathFraction/(depV+appV);
  if(bbb<0.) std::cout << "ChargeDrifterFP420:ERROR: check your Voltage for log(bbb) bbb=" << bbb << std::endl;
  double driftTime = -constTe*log(bbb) + startT0;  
  //    log(1.-2*depV*pathFraction/(depV+appV)) + startT0;  
  // since no magnetic field the Sigma_x, Sigma_y are the same =  sigma  !!!!!!!!!!!!!!!!! 
  double sigma = sqrt(2.*constDe*driftTime*100.);  //  * 100.  - since constDe is [cm2/sec], but i want [mm2/sec]
  
  //  std::cout << "ChargeDrifterFP420: driftTime=  " << driftTime << "  pathFraction=  " << pathFraction << "  constTe=  " << constTe << "  sigma=  " << sigma << std::endl;
  if(verbo>0) {
    std::cout << "ChargeDrifterFP420: drift: xytype=" << xytype << "pathFraction=" << pathFraction << std::endl;
    std::cout << "  constTe= " << constTe  << "  driftTime = " << driftTime << " startT0  = " << startT0 <<  std::endl;
    std::cout << " log = " << log(1.-2*depV*pathFraction/(depV+appV))  << std::endl;
    std::cout << " negativ inside log = " << -2*depV*pathFraction/(depV+appV)  << std::endl;
    std::cout << " constDe = " << constDe  << "  sigma = " << sigma << std::endl;
    
    std::cout << "ChargeDrifterFP420: drift: xytype=" << xytype << "pathValue=" << pathValue << std::endl;
    std::cout << " tanShiftAngleX = " << tanShiftAngleX  << "  tanShiftAngleY = " << tanShiftAngleY  << "  tanShiftAngleZ = " << tanShiftAngleZ  << std::endl;
    std::cout << "sliX = " << sliX  << "  sliY = " << sliY  << "  sliZ = " << sliZ  << std::endl;
    std::cout << "pathFraction = " << pathFraction  << "  driftTime = " << driftTime  << std::endl;
    std::cout << "sigma = " << sigma  << std::endl;
    std::cout << "xDriftDueToField = " << xDriftDueToField  << "  yDriftDueToField = " << yDriftDueToField  << "  zDriftDueToField = " << zDriftDueToField  << std::endl;
    std::cout << "2currX = " << currX  << "  2currY = " << currY  << "  2currZ = " << currZ  << std::endl;
    std::cout << "ChargeDrifterFP420: drift; finally, rETURN AmplitudeSlimentFP420" << std::endl;
    std::cout << "===================================================================" << std::endl;
    std::cout << " (edu).energy()= " << (edu).energy() << std::endl;
    std::cout << "==" << std::endl;
  }
  //  std::cout << "ChargeDrifterFP420: (edu).energy()= " << (edu).energy() << std::endl;
  return AmplitudeSegmentFP420(currX,currY,currZ,sigma,
                   (edu).energy());  
}

Member Data Documentation

double ChargeDrifterFP420::appV

private

Definition at line 22 of file ChargeDrifterFP420.h.

double ChargeDrifterFP420::constDe

private

Definition at line 18 of file ChargeDrifterFP420.h.

double ChargeDrifterFP420::constTe

private

Definition at line 17 of file ChargeDrifterFP420.h.

double ChargeDrifterFP420::depV

private

Definition at line 21 of file ChargeDrifterFP420.h.

double ChargeDrifterFP420::ldriftcurrX

private

Definition at line 23 of file ChargeDrifterFP420.h.

double ChargeDrifterFP420::ldriftcurrY

private

Definition at line 24 of file ChargeDrifterFP420.h.

double ChargeDrifterFP420::modulePath

private

Definition at line 16 of file ChargeDrifterFP420.h.

double ChargeDrifterFP420::startT0

private

Definition at line 20 of file ChargeDrifterFP420.h.

double ChargeDrifterFP420::temperature

private

Definition at line 19 of file ChargeDrifterFP420.h.

int ChargeDrifterFP420::verbo

private

Definition at line 25 of file ChargeDrifterFP420.h.