musrIBackwardHit.cc

Go to the documentation of this file.

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

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

G4Allocator<musrIBackwardHit> musrIBackwardHitAllocator;

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

musrIBackwardHit::musrIBackwardHit() {}

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

musrIBackwardHit::~musrIBackwardHit() {}

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
G4int musrIBackwardHit::eventIDoldIB = -1;
G4int musrIBackwardHit::NIB = 0;
G4int musrIBackwardHit::BMultihit = 0;
G4int musrIBackwardHit::runIDoldIB = -1;
G4int musrIBackwardHit::IBchamberNbold = -1;

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

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

}

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

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

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

  return *this;
}

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

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

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

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

  G4RunManager* fRunManager = G4RunManager::GetRunManager();
  G4int eventID = fRunManager->GetCurrentEvent()->GetEventID();
  G4int runID   = fRunManager->GetCurrentRun()->GetRunID();
  BMultihitSwitch=0;
        if (runID != runIDoldIB)
          {
            NIB=0;
            BMultihit = 0;
            eventIDoldIB = -1;
            IBchamberNbold = -1;
          }

    if (particle_name== "e+" and (eventIDoldIB != eventID or IBchamberNbold != IBchamberNb))
      {



         if (eventIDoldIB == eventID)
           {
             BMultihit++;
             BMultihitSwitch=1;

           }
         NIB++;

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

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

  ofstream posfile1;
  posfile1.open ("ibackward.dat", ios::out | ios::app);
  posfile1 << fabs(B[2]) <<" "<< pos.x()/cm<<" "<<pos.y()/cm <<" "<< pos.z()/cm  
           << " "  << globaltime/s << " "<< IBchamberNb 
           << " " << BMultihitSwitch 
           << G4endl;
  posfile1.close();

      }

  eventIDoldIB=eventID;
  runIDoldIB = runID;
  IBchamberNbold = IBchamberNb;  


}

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

---