#include <Ext_track.h>

Public Member Functions
	Ext_track ()

	Ext_track (const bool msgFlag, const bool myBFieldOn, const bool GeomOptimization, const int m_tofversion, const bool UseMucKal, const int MucWindow)

	~Ext_track ()

void	Initialization (const bool aMsgFlag, const bool Bfield, const bool GeomOptimization, const bool aUseMucKal, const int aMucWindow)

bool	Set (const Hep3Vector &xv3, const Hep3Vector &pv3, const HepSymMatrix &err, const std::string &particleName, const double pathInMDC, const double tofInMdc)

void	TrackExtrapotation ()

void	Dir (const int dir)

ExtSteppingAction *	GetStepAction ()

Detailed Description

Definition at line 40 of file Ext_track.h.

Constructor & Destructor Documentation

◆ Ext_track() [1/2]

Ext_track::Ext_track ( )

Definition at line 35 of file Ext_track.cxx.

                     : myMsgFlag(true),myBFieldOn(true),myGeomOptimization(true),m_dir(0),m_tofversion(2),myUseMucKal(true),myMucWindow(6) 
{ 
//   BesSensitiveManager *besSensitiveManager = new BesSensitiveManager;
   bes3DetectorConstruction = new ExtBesDetectorConstruction(myBFieldOn,m_tofversion);
   bes3WorldVolume = bes3DetectorConstruction->Construct();
   extPhysicsList = new ExtPhysicsList;
   extTrack = new G4Track;
   
   //for geant4.8.1, move this line to Initialization, extSteppingAction = new ExtSteppingAction;
   extTrackingManager = new G4TrackingManager;
   //RunManagerKernel for geant4.9.0
   extRunManagerKernel = new G4RunManagerKernel;
}

◆ Ext_track() [2/2]

Ext_track::Ext_track	(	const bool	msgFlag,
		const bool	myBFieldOn,
		const bool	GeomOptimization,
		const int	m_tofversion,
		const bool	UseMucKal,
		const int	MucWindow
	)

Definition at line 49 of file Ext_track.cxx.

                                                                                                                                                            : myMsgFlag(msgFlag),myBFieldOn(BFieldOn),myGeomOptimization(GeomOptimization),m_dir(0),myUseMucKal(aUseMucKal),myMucWindow(aMucWindow) 
{ 
//   BesSensitiveManager *besSensitiveManager = new BesSensitiveManager;
   bes3DetectorConstruction = new ExtBesDetectorConstruction(myBFieldOn,m_tofversion);
   bes3WorldVolume = bes3DetectorConstruction->Construct();
   extPhysicsList = new ExtPhysicsList;
   extTrack = new G4Track;
   
   //for geant4.8.1, move this line to Initialization, extSteppingAction = new ExtSteppingAction;
   extTrackingManager = new G4TrackingManager;
   //RunManagerKernel for geant4.9.0
   extRunManagerKernel = new G4RunManagerKernel;
}

◆ ~Ext_track()

Ext_track::~Ext_track ( )

Definition at line 65 of file Ext_track.cxx.

{
   if(extRunManagerKernel) delete extRunManagerKernel;
   if(extTrackingManager) delete extTrackingManager;
// if(extSteppingAction) delete extSteppingAction;
   if(extTrack) delete extTrack;
   if(bes3DetectorConstruction) delete bes3DetectorConstruction;
   if(extPhysicsList) delete extPhysicsList;
   
   // open geometry for deletion
   G4GeometryManager::GetInstance()->OpenGeometry();
 
   // deletion of Geant4 kernel classes
   G4SDManager* fSDM = G4SDManager::GetSDMpointerIfExist();
   if(fSDM)   
     {
    delete fSDM;
     }
}

Member Function Documentation

◆ Dir()

void Ext_track::Dir ( const int dir )

inline

Termination of the track extrapolation. Accessors for the parameters. Modifiers for the parameters. Modifiers for the track direction. (0/1)=(outgoing/incoming).

Definition at line 154 of file Ext_track.h.

154{ m_dir = dir; }

◆ GetStepAction()

ExtSteppingAction * Ext_track::GetStepAction ( )

inline

Definition at line 90 of file Ext_track.h.

90{return extSteppingAction;};

Referenced by TrkExtAlg::execute().

◆ Initialization()

void Ext_track::Initialization	(	const bool	aMsgFlag,
		const bool	Bfield,
		const bool	GeomOptimization,
		const bool	aUseMucKal,
		const int	aMucWindow
	)

Definition at line 89 of file Ext_track.cxx.

{
    myMsgFlag=aMsgFlag;
    myGeomOptimization = GeomOptimization;
    myBFieldOn=Bfield,
        myUseMucKal=aUseMucKal;
    myMucWindow = aMucWindow;
  //add for geant4.8.1
  G4ParticleTable::GetParticleTable()->SetReadiness();
  extPhysicsList->ConstructParticle();
 
     if(myMsgFlag) cout << "Ext_track::Init will execute geant initialization." << endl;
     if(!GeometryInitialization()) cout << "Error in Ext_track::GeometryInitialization()" << endl;
     PhysicsInitialization();  
 
  extSteppingAction = new ExtSteppingAction;  
    extSteppingAction->SetMsgFlag(aMsgFlag);
    extSteppingAction->SetMucKalFlag(aUseMucKal);
    extSteppingAction->SetMucWindow(aMucWindow);
  //Set extSteppingAction
  extTrackingManager->SetUserAction(extSteppingAction);
 
   
 
}

Referenced by TrkExtAlg::initialize().

◆ Set()

bool Ext_track::Set	(	const Hep3Vector &	xv3,
		const Hep3Vector &	pv3,
		const HepSymMatrix &	err,
		const std::string &	particleName,
		const double	pathInMDC,
		const double	tofInMdc
	)

Setup the error skip list as a part of the initialization. errskip_list – Name of the error skip file. Print out the contents of the error skip list.

Definition at line 234 of file Ext_track.cxx.

{ 
  if( err.num_row() != 6 ){          // ?static const int Ndim_err=6, see Ext_errmx.h line58 
    std::cerr << "%ERROR at Ext_track::Set. Dimension of error matrix: "
    << err.num_row() << " should be 6"  << std::endl;
    exit( 0 );
  }
 
  m_vect[0] = xv3.x();// ?set starting position,private data
  m_vect[1] = xv3.y();
  m_vect[2] = xv3.z();
 
//  m_errskip_flag = 0;
//  m_errskip_level = 0;
 
 
  //  Check the starting point is inside the setup.   
  if(!CheckVertexInsideWorld(xv3)) return 0; 
 
  float p( pv3.mag() );
  m_vect[3] = pv3.x()/p; //?set direction of momentum
  m_vect[4] = pv3.y()/p;
  m_vect[5] = pv3.z()/p;
  m_vect[6] = p;
 
  // check Particlename
  if(particleName!="e+"&&particleName!="e-"&&
     particleName!="mu+"&&particleName!="mu-"&&
     particleName!="pi+"&&particleName!="pi-"&&
     particleName!="kaon+"&&particleName!="kaon-"&&
     particleName!="proton"&&particleName!="anti_proton"&&
     particleName!="gamma")
     {
    std::cerr <<"Unknown or unconstructed Particle."<< std::endl;
    return 0;
     }
     
     
  double        mass;
  double        Q; 
   
  G4ParticleDefinition* particleDefinition=G4ParticleTable::GetParticleTable()->FindParticle(particleName); 
  Q = particleDefinition->GetPDGCharge(); 
  mass = particleDefinition->GetPDGMass();
   
 
  Hep3Vector xv( m_vect[0], m_vect[1], m_vect[2] );
  Hep3Vector pv(m_vect[3]*m_vect[6], m_vect[4]*m_vect[6], m_vect[5]*m_vect[6]);
 
  m_xp_err.set_err( err, xv, pv, Q, mass ); // Set error matrix.
   
  extSteppingAction->SetXpErrPointer(&m_xp_err);
  extSteppingAction->SetInitialPath(pathInMDC);
  extSteppingAction->SetInitialTof(tofInMdc);
      
  double betaInMDC = p/sqrt(mass*mass+p*p);//velocity
  extSteppingAction->SetBetaInMDC(betaInMDC);
//  double tofInMDC = pathInMDC/(betaInMDC*299.792458);
//  if(myMsgFlag) cout<<"TOF in MDC: "<<tofInMDC<<endl;
   
//  extSteppingAction->Reset();
    extSteppingAction->MucReset();
  // extTrack Initialization.
 
/*  // comment 2008.04.07 due to memory loss       
    // Initialize a G4PrimaryParticle.   
     G4PrimaryParticle* primaryParticle = new G4PrimaryParticle(particleDefinition,pv3.x(),pv3.y(),pv3.z());
     primaryParticle->SetMass(mass);
     primaryParticle->SetCharge(Q);
       
    // Initialize a G4DynamicParticle.
//   G4DynamicParticle* DP = new G4DynamicParticle(particleDefinition,primaryParticle->GetMomentum());
//   DP->SetPrimaryParticle(primaryParticle);
*/
     G4DynamicParticle* DP = new G4DynamicParticle(particleDefinition,pv);
  
   delete extTrack; // add on 2008.04.07 to avoid memory loss
   extTrack = new G4Track(DP,0.0,xv3);
// extTrack->CopyTrackInfo(G4Track::G4Track(DP,0.0,xv3));
  
  // Reset navigator   
  Hep3Vector center(0,0,0);
  G4Navigator* navigator= G4TransportationManager::GetTransportationManager()->GetNavigatorForTracking();
  navigator->LocateGlobalPointAndSetup(center,0,false);
 
  return 1;
}

Referenced by TrkExtAlg::execute().

◆ TrackExtrapotation()

void Ext_track::TrackExtrapotation ( )

Definition at line 342 of file Ext_track.cxx.

{ 
   extTrackingManager->ProcessOneTrack(extTrack);
}

Referenced by TrkExtAlg::execute().

The documentation for this class was generated from the following files:

Public Member Functions

Detailed Description

Constructor & Destructor Documentation

◆ Ext_track() [1/2]

◆ Ext_track() [2/2]

◆ ~Ext_track()

Member Function Documentation

◆ Dir()

◆ GetStepAction()

◆ Initialization()

◆ Set()

◆ TrackExtrapotation()