musrIForwardHit.cc

Go to the documentation of this file.

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

#include "musrIForwardHit.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"
#include "G4RunManager.hh"
#include "G4Run.hh"

G4Allocator<musrIForwardHit> musrIForwardHitAllocator;

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

musrIForwardHit::musrIForwardHit() {}

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

musrIForwardHit::~musrIForwardHit() {}

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
G4int musrIForwardHit::eventIDoldIF = -1;
G4int musrIForwardHit::NIF = 0;
G4int musrIForwardHit::FMultihit = 0;
G4int musrIForwardHit::runIDoldIF = -1;
G4int musrIForwardHit::IFchamberNbold = -1;

musrIForwardHit::musrIForwardHit(const musrIForwardHit& right)
  : G4VHit()
{

  trackID   = right.trackID;
  particle_name = right.particle_name;
  edep      = right.edep;
  pos       = right.pos;
  pol       = right.pol;
  logicalvolume = right.logicalvolume;
  IFchamberNb = right.IFchamberNb;
  globaltime  = right.globaltime;

}

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

const musrIForwardHit& musrIForwardHit::operator=(const musrIForwardHit& right)
{

  trackID   = right.trackID;
  particle_name = right.particle_name;
  edep      = right.edep;
  pos       = right.pos;
  pol       = right.pol;
  logicalvolume=right.logicalvolume;
  IFchamberNb = right.IFchamberNb;
  globaltime = right.globaltime;

  return *this;
}

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

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

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

void musrIForwardHit::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 musrIForwardHit::Print()
{

  G4RunManager* fRunManager = G4RunManager::GetRunManager();
  G4int eventID = fRunManager->GetCurrentEvent()->GetEventID();
  G4int runID   = fRunManager->GetCurrentRun()->GetRunID();
  FMultihitSwitch=0;



        if (runID != runIDoldIF)
          {
            NIF=0;
            FMultihit = 0;
            eventIDoldIF = -1;
            IFchamberNbold = -1;
          }

    if (particle_name== "e+" and (eventIDoldIF != eventID or IFchamberNbold != IFchamberNb))

      {



         if (eventIDoldIF == eventID)
           {
             FMultihitSwitch=1;
             FMultihit++;
             //        G4cout << "Forward Multiple Hit!, No = "<< FMultihit<<G4endl;
           }

         NIF++;
         //      G4cout << "NIF is "<< NIF<<G4endl;
         //      G4cout << "(Forward) eventID is "<< eventID<<G4endl;



        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 << "  Segment: " << IFchamberNb << G4endl;
      //      G4cout <<"Position " << pos.x()/cm<<" "<<pos.y()/cm <<" "<< pos.z()/cm <<G4endl;

      //         G4cout << "Field is "<< B[2]<<G4endl;

  ofstream posfile2;
  posfile2.open ("iforward.dat", ios::out | ios::app);
  posfile2 << fabs(B[2]) <<" "<< pos.x()/cm<<" "<<pos.y()/cm <<" "<< pos.z()/cm  
           << " "  << globaltime/s <<" "<< IFchamberNb 
           << " "<< FMultihitSwitch
           << G4endl;
  posfile2.close();

      }

  eventIDoldIF=eventID;
  runIDoldIF = runID;
  IFchamberNbold = IFchamberNb;  



}

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

---