MJWaveformChargePoint.cc

Go to the documentation of this file.

//---------------------------------------------------------------------------//
//bb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nu//
//                                                                           //
//                         MAJORANA Simulation                               //
//                                                                           //
//      This code implementation is the intellectual property of the         //
//      MAJORANA Collaboration. It is based on Geant4, an intellectual       //
//      property of the RD44 GEANT4 collaboration.                           //
//                                                                           //
//                        *********************                              //
//                                                                           //
//    Neither the authors of this software system, nor their employing       //
//    institutes, nor the agencies providing financial support for this      //
//    work  make  any representation or  warranty, express or implied,       //
//    regarding this software system or assume any liability for its use.    //
//    By copying, distributing or modifying the Program (or any work based   //
//    on on the Program) you indicate your acceptance of this statement,     //
//    and all its terms.                                                     //
//                                                                           //
//bb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nu//
//---------------------------------------------------------------------------//
//                                                          
// $Id: MJWaveformChargePoint.cc,v 1.2 2004/11/09 13:42:39 xliu Exp $ 
//      
// CLASS IMPLEMENTATION:  MJWaveformChargePoint.cc
//
//---------------------------------------------------------------------------//
// 
//---------------------------------------------------------------------------//
//---------------------------------------------------------------------------//
//

#include <iostream> 
#include <CLHEP/Vector/ThreeVector.h>

#include "waveform/MJVWaveformCrystal.hh"

//---------------------------------------------------------------------------//

#include "waveform/MJWaveformChargePoint.hh"

using namespace std;

//---------------------------------------------------------------------------//

MJWaveformChargePoint::MJWaveformChargePoint(HepInt *location,
                                             HepInt sign, 
                                             HepDouble eloss,
                                             MJVWaveformCrystal *Xtal):
fChargeSign(sign), fEnergyLoss(eloss),
fNextActiveChargePoint(0), fNextChargePoint(0),
fPreviousActiveChargePoint(0), fPreviousChargePoint(0), 
fXtal(Xtal)
{
  fGridLocation[0] = location[0];
  fGridLocation[1] = location[1];
  fGridLocation[2] = location[2];
}

//---------------------------------------------------------------------------//

MJWaveformChargePoint::MJWaveformChargePoint(const MJWaveformChargePoint&other)
{;}

//---------------------------------------------------------------------------//

MJWaveformChargePoint::~MJWaveformChargePoint()
{;}

//---------------------------------------------------------------------------//
 
void MJWaveformChargePoint::DoDiffusion(void)
{
  if(!(fNumberOfSteps % fDiffusionStepInterval)) {
    // Do diffusion. Not implemented yet.
  }
}

//---------------------------------------------------------------------------//
 
void MJWaveformChargePoint::DoOneStep(void)
{
  if(fActive) {
    PropagateOneStep();
    if(fTrappingOn) DoTrapping();
    if(fDiffusionOn) DoDiffusion();
  }
}

//---------------------------------------------------------------------------//

void MJWaveformChargePoint::DoTrapping(void)
{
  if(!(fNumberOfSteps % fTrappingStepInterval)) {
    // Do Trapping. Not implemented yet.
  }
}

//---------------------------------------------------------------------------//

void MJWaveformChargePoint::Dump(void)
{
  cout << "\n\nActive: " << fActive;
  cout << "\nEnergy loss (MeV): " << fEnergyLoss / MeV;
  cout << "\nCharge (C): " << fCharge / coulomb;
  cout << "\nCharge sign : " << fChargeSign;
  cout << "\nCurrent pos (cm): " << fCurrentPosition / cm;
  cout << "\nPrevious pos (cm): " << fPreviousPosition / cm;
  cout << "\nOriginal pos (cm): " << fOriginalPosition / cm;
  cout << "\nGridlocation : " << fGridLocation[0] << fGridLocation[1];
  cout << fGridLocation[2];
  cout << "\nGrid Point Size (cm) : " << fGridPointSize / cm;
  cout << "\nNumber Of Steps : " << fNumberOfSteps;
  cout << "\nSigma : " << fSigma;
  cout << "\nDiffusion Step Interval : " << fDiffusionStepInterval;
  cout << "\nTrapping Step Interval : " << fTrappingStepInterval;
  cout << "\nTime Step (ns): " << fTimeStep / ns;
  cout << "\nTotal Time (ns) : " << fTotalTime / ns;
  Hep3Vector vel;
  GetVelocity(&vel);
  cout << "\nVelocity (cm/us) : " << vel / cm * microsecond;
  cout << "\n\n";
}

//---------------------------------------------------------------------------//

void MJWaveformChargePoint::InitializePropagation(void)
{
  fActive = !(fXtal->OutsideCrystal(&fCurrentPosition));
  fCharge = fXtal->ConvertEnergyLossToCharge(fEnergyLoss);
  fOriginalPosition = fXtal->ConvertGridpointToCoordinate(fGridLocation);
  fPreviousPosition = fCurrentPosition = fOriginalPosition;
  fGridPointSize = fXtal->GetGridPointSize();
  fNumberOfSteps = 0;
  fSigma = fXtal->GetInitialGridPointSigma();
  if(fDiffusionOn = fXtal->DiffusionOn()) 
    fDiffusionStepInterval = fXtal->GetDiffusionStepInterval();
  if(fTrappingOn = fXtal->TrappingOn())
    fTrappingStepInterval = fXtal->GetTrappingStepInterval();
  fTimeStep = fXtal->GetTimeStep();
  
}

//---------------------------------------------------------------------------//

void MJWaveformChargePoint::PropagateOneStep(void)
{
  // Note: Calling method has to determine whether this charge point is still
  // active.
    static Hep3Vector vel;

    GetVelocity(&vel);
    fPreviousPosition = fCurrentPosition;
    fCurrentPosition +=  -fTimeStep * fChargeSign * vel; 
    fNumberOfSteps++;
    fTotalTime += fTimeStep;
    if(fXtal->OutsideCrystal(&fCurrentPosition)) {
      fActive = false;
      RemoveFromActiveList();
    }
}

//---------------------------------------------------------------------------//

void MJWaveformChargePoint::RemoveFromActiveList(void)
{ 
  if(fPreviousActiveChargePoint && fNextActiveChargePoint) {
    fPreviousActiveChargePoint->
      SetNextActiveChargePoint(fNextActiveChargePoint);
    fNextActiveChargePoint->
      SetPreviousActiveChargePoint(fPreviousActiveChargePoint);
  } else {
    if(fNextActiveChargePoint) 
      fNextActiveChargePoint->SetPreviousActiveChargePoint(NULL);
    else if (fPreviousActiveChargePoint) 
      fPreviousActiveChargePoint->SetNextActiveChargePoint(NULL);
  }
}

//---------------------------------------------------------------------------//
//---------------------------------------------------------------------------//

---