musrTrackerHit.cc

Go to the documentation of this file.

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

#include "musrTrackerHit.hh"
#include "G4UnitsTable.hh"
#include "G4VVisManager.hh"
#include "G4Circle.hh"
#include "G4Colour.hh"
#include "G4VisAttributes.hh"
#include <fstream.h>
#include <iostream.h>
#include "G4ios.hh"
#include "G4MagneticField.hh"
#include "G4FieldManager.hh"
#include "G4TransportationManager.hh"
#include <iomanip>
#include "globals.hh"
#include "G4Transform3D.hh"
#include "G4ProcessManager.hh"
#include "G4Track.hh"
#include "G4ThreeVector.hh"

G4Allocator<musrTrackerHit> musrTrackerHitAllocator;

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

musrTrackerHit::musrTrackerHit() {}

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

musrTrackerHit::~musrTrackerHit() {}

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

musrTrackerHit::musrTrackerHit(const musrTrackerHit& right)
  : G4VHit()
{
  particle_name = right.particle_name;
  trackID   = right.trackID;
  edep      = right.edep;
  pos       = right.pos;
  pol       = right.pol;
}

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

const musrTrackerHit& musrTrackerHit::operator=(const musrTrackerHit& right)
{
  particle_name = right.particle_name;
  trackID   = right.trackID;
  edep      = right.edep;
  pos       = right.pos;
  pol       = right.pol;
  return *this;
}

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

G4int musrTrackerHit::operator==(const musrTrackerHit& right) const
{
  return (this==&right) ? 1 : 0;
}

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

void musrTrackerHit::Draw()
{
  G4VVisManager* pVVisManager = G4VVisManager::GetConcreteInstance();
  if(pVVisManager)
  {
    G4Circle circle(pos);
    circle.SetScreenSize(0.04);
    circle.SetFillStyle(G4Circle::filled);
    G4Colour colour(1.,0.,0.);
    G4VisAttributes attribs(colour);
    circle.SetVisAttributes(attribs);
    pVVisManager->Draw(circle);
  }
}

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

void musrTrackerHit::Print()
{

  if (trackID==1)
    {

        G4FieldManager *fMgr=G4TransportationManager::GetTransportationManager()->GetFieldManager();


          point[0]=0.;
          point[1]=0.;
          point[2]=0.;

          B[2]=0.0;
              

      if(!fMgr->DoesFieldChangeEnergy())//then we have a magnetic field
        {
          mfield = fMgr->GetDetectorField();
          mfield->GetFieldValue(point,B);
          B[0]=B[0]*1000;
          B[1]=B[1]*1000;
          B[2]=B[2]*1000;
          //          G4cout <<"Magnetic field at origin B ="<< B[2] << " T"<<G4endl;
}
      else{
        //          G4cout <<"Magnetic field at origin B = 0 T" << G4endl;
      }


  ofstream mufile1;
  mufile1.open ("muposition.dat", ios::out | ios::app);
  mufile1 << fabs(B[2]) <<" "<< pos.x()/cm<<" "<<pos.y()/cm <<" "<< pos.z()/cm  << G4endl;
  mufile1.close();

    ofstream mufile2;
    mufile2.open ("mupol.dat", ios::out | ios::app);
  mufile2 << pos.x()/m<<" "<<pos.y()/m <<" "<< pos.z()/m<<" " 
          << pol.x()<<" "<<pol.y()<<" "<< pol.z() <<G4endl;
  mufile2.close();
    }
}

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

---