MJWaveformPlanarDetector.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: MJWaveformPlanarDetector.cc,v 1.2 2004/11/09 13:42:39 xliu Exp $ 
//      
// CLASS IMPLEMENTATION:  MJWaveformPlanarDetector.cc
//
//---------------------------------------------------------------------------//
// 
//---------------------------------------------------------------------------//
//---------------------------------------------------------------------------//
//

#include <CLHEP/Units/SystemOfUnits.h>

#include "io/MJLogger.hh"
#include "waveform/MJWaveformChargePoint.hh"
#include "waveform/MJWaveformPlanarCrystal.hh"
#include "waveform/MJWaveformPlanarSegment.hh"

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

#include "waveform/MJWaveformPlanarDetector.hh"  

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

MJWaveformPlanarDetector::MJWaveformPlanarDetector():
  fData(0), fPlanarCrystal1(0), fPlanarCrystal2(0), fStep(0)
{;}

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

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

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

MJWaveformPlanarDetector::~MJWaveformPlanarDetector()
{
  delete fData;
}

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

void MJWaveformPlanarDetector::BeginOfEventAction()
{
  Hep3Vector   coord;


  coord.setX(5.0 * cm);
  coord.setY(5.0 * cm);
  coord.setZ((2.5 + 0.5*fStep++) * cm);
  HepInt crap[3];
  HepInt *location;
  location = crap;
  fPlanarCrystal1->ConvertCoordinateToGridpoint(&coord, location);
  HepInt sign;
  HepDouble Eloss;
  MJWaveformChargePoint *epoint =  new MJWaveformChargePoint(
                                       location,
                                       sign = -1,
                                       Eloss = 2.0 * MeV,
                                       fPlanarCrystal1);
  MJWaveformChargePoint *hpoint =  new MJWaveformChargePoint(
                                       location,
                                       sign = +1,
                                       Eloss = 2.0 * MeV,
                                       fPlanarCrystal1);
  epoint->InitializePropagation();
  hpoint->InitializePropagation();
  fPlanarCrystal1->AddChargePoint(epoint);
  fPlanarCrystal1->AddChargePoint(hpoint);

  coord.setX(5.0 * cm);
  coord.setY(5.0 * cm);
  coord.setZ(6.3 * cm);
  fPlanarCrystal1->ConvertCoordinateToGridpoint(&coord, location);
  epoint =  new MJWaveformChargePoint(
                                       location,
                                       sign = -1,
                                       Eloss = 4.0 * MeV,
                                       fPlanarCrystal1);
  hpoint =  new MJWaveformChargePoint(
                                       location,
                                       sign = +1,
                                       Eloss = 4.0 * MeV,
                                       fPlanarCrystal1);
  epoint->InitializePropagation();
  hpoint->InitializePropagation();
  fPlanarCrystal1->AddChargePoint(epoint);
  fPlanarCrystal1->AddChargePoint(hpoint);

  fPlanarCrystal1->GetFirstSegment()->ResetCounters();
  fPlanarCrystal1->GetFirstSegment()->GetNextSegment()->ResetCounters();
}

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

void MJWaveformPlanarDetector::Construct()
{
  // Construct two planar crystals
  fPlanarCrystal1 = new MJWaveformPlanarCrystal;
  fPlanarCrystal2 = new MJWaveformPlanarCrystal;
  fPlanarCrystal1->Construct();
  fPlanarCrystal2->Construct();
  AddCrystal(fPlanarCrystal1);
  AddCrystal(fPlanarCrystal2);
  
  // Create Root output files and trees.
  fData = new TreeStructure;
  fRootFile = new TFile("crap.root", "RECREATE");
  fRootTree = new TTree("fTree", "A Tree, A Beautiful Tree");
  fRootTree->Branch("Time", fData->Time, "Time[500]/F");
  fRootTree->Branch("Current", fData->Current, "Current[2][500]/F");
  fRootTree->Branch("Charge", fData->Charge, "Charge[2][500]/F");
  fRootTree->Branch("InitialPosition", fData->InitialPosition, 
                    "InitialPosition[10][3]/F");
  fRootTree->Branch("InitialCharge", fData->InitialCharge, 
                    "InitialCharge[10]/F");
  fRootTree->Branch("NumberOfPoints", &fData->NumberOfPoints,
                    "NumberOfPoints/I");
}

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

void MJWaveformPlanarDetector::EndOfEventAction()
{
  GenerateSignals();
  MJVWaveformSegment *seg1 = fPlanarCrystal1->GetFirstSegment();
  MJVWaveformSegment *seg2 = seg1->GetNextSegment();
  HepDouble timestep = fPlanarCrystal1->GetTimeStep();
  HepDouble *curr1 = seg1->GetCurrentSignal();
  HepDouble *chrg1 = seg1->GetChargeSignal();
  HepDouble *curr2 = seg2->GetCurrentSignal();
  HepDouble *chrg2 = seg2->GetChargeSignal();
  HepInt np = fData->NumberOfPoints = seg1->GetChargeStep();
  for(HepInt i = 0; i < np; i++) {
    fData->Time[i] = CLHEPtoROOT(timestep * i / ns);
    fData->Current[0][i] = CLHEPtoROOT(curr1[i] / microampere);
    fData->Current[1][i] = CLHEPtoROOT(curr2[i] / microampere);
    fData->Charge[0][i] = CLHEPtoROOT(chrg1[i] / eplus);
    fData->Charge[1][i] = CLHEPtoROOT(chrg2[i] / eplus);
  }
  fRootTree->Fill();
}

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

void MJWaveformPlanarDetector::EndOfRunAction()
{
  fRootTree->Write();
  fRootFile->Write();
  fRootFile->Close();
  fRootFile->Delete();
}

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

---