BesPrimaryGeneratorAction Class Reference

#include <BesPrimaryGeneratorAction.hh>

Inheritance diagram for BesPrimaryGeneratorAction:

Public Member Functions
	BesPrimaryGeneratorAction ()
	~BesPrimaryGeneratorAction ()
void	GeneratePrimaries (G4Event *anEvent)
G4String	GetGeneratorName ()
void	SetGeneratorName (G4String genNew)
G4String	GetGenbesName ()
void	SetGenbesName (G4String genbesNew)
void	SetNParticle (G4int np)
void	SetParticleName (G4String pname)
void	SetMinCos (G4double min1)
void	SetMaxCos (G4double max1)
void	SetPhiStart (G4double ps)
void	SetPhiEnd (G4double pe)
void	SetMomentum (G4double pm)
void	SetDeltaP (G4double dp)
void	SetPosX (G4double x)
void	SetPosY (G4double y)
void	SetPosZ (G4double z)
	BesPrimaryGeneratorAction ()
	~BesPrimaryGeneratorAction ()
void	GeneratePrimaries (G4Event *anEvent)
G4String	GetGeneratorName ()
void	SetGeneratorName (G4String genNew)
G4String	GetGenbesName ()
void	SetGenbesName (G4String genbesNew)
void	SetNParticle (G4int np)
void	SetParticleName (G4String pname)
void	SetMinCos (G4double min1)
void	SetMaxCos (G4double max1)
void	SetPhiStart (G4double ps)
void	SetPhiEnd (G4double pe)
void	SetMomentum (G4double pm)
void	SetDeltaP (G4double dp)
void	SetPosX (G4double x)
void	SetPosY (G4double y)
void	SetPosZ (G4double z)

Detailed Description

Definition at line 26 of file InstallArea/include/GenSim/GenSim/BesPrimaryGeneratorAction.hh.

Constructor & Destructor Documentation

BesPrimaryGeneratorAction() [1/2]

BesPrimaryGeneratorAction::BesPrimaryGeneratorAction

(

)

Definition at line 24 of file BesPrimaryGeneratorAction.cc.

{
  nParticle = 1;
  particleName = "pi-";
  minCos = -0.8;
  maxCos = 0.8;
  phiStart = 0.;
  phiEnd = 360.;
  pMomentum = 1.; 
  deltaP = 0.;
  posX=0;
  posY=0;
  posZ=0;
 
  particleGun = new G4ParticleGun(1);
  isGenbes = true;
  HEPEvt = 0;
  generatorName = "tester";  
  //TESTER parameters passed by this messenger  
  messenger = new BesPrimaryGeneratorMessenger(this);
  
  if(generatorName == "cosmic")
  {
    std::string path = getenv("GENSIMROOT");
    G4cout<<"path: "<<path<<G4endl;
 
    path += "/root/";
    std::string pFile = path + "ppdc.root";
    std::string thetaFile = path + "theta.root";
    std::string phiFile = path +"phi.root";
  
    TFile* f1 = new TFile(pFile.c_str());
    h1 = (TH1F*)f1->Get("htemp");
  
    TFile* f2 = new TFile(thetaFile.c_str());
    h2 = (TH1F*)f2->Get("htemp");
  
    TFile* f3 = new TFile(phiFile.c_str());
    h3 = (TH1F*)f3->Get("htemp");
  
    //ftest = new TFile("ftest.root","recreate");
    //tuple = new TNtuple("test","test","p:theta:phi"); 
    //counter = 0;
  }
}

~BesPrimaryGeneratorAction() [1/2]

BesPrimaryGeneratorAction::~BesPrimaryGeneratorAction

(

)

Definition at line 70 of file BesPrimaryGeneratorAction.cc.

{
  delete particleGun;
  if(messenger) delete messenger; 
  if(generatorName=="genbes")delete HEPEvt; 
}

BesPrimaryGeneratorAction() [2/2]

BesPrimaryGeneratorAction::BesPrimaryGeneratorAction

(

)

~BesPrimaryGeneratorAction() [2/2]

BesPrimaryGeneratorAction::~BesPrimaryGeneratorAction

(

)

Member Function Documentation

GeneratePrimaries() [1/2]

void BesPrimaryGeneratorAction::GeneratePrimaries

(

G4Event *

anEvent

)

Definition at line 77 of file BesPrimaryGeneratorAction.cc.

{
  if(generatorName=="tester")
  {
    G4ParticleTable* particleTable = G4ParticleTable::GetParticleTable();
    G4ParticleDefinition* particle
      = particleTable->FindParticle(particleName);
    particleGun->SetParticleDefinition(particle);
    particleGun->SetParticlePosition(G4ThreeVector(posX,posY,posZ));
        
    //set TESTER parameters
    for(G4int i=0; i<nParticle; i++)
    {   
      G4double pMag = pMomentum;
      //randomize p
      if(deltaP>0.)pMag = pMomentum - deltaP*(1.0 - 2.0*G4UniformRand());
      pMag = pMag*GeV;
      //randomize cos(theta)
      G4double costheta = minCos +(maxCos - minCos)*G4UniformRand();
      //randomize phi
      G4double phi = phiStart+(phiEnd-phiStart)*G4UniformRand();
      phi = phi*degree;
      G4double sintheta = sqrt(1.-costheta*costheta);
      //computer 3-vector momentum
      G4ParticleMomentum aMomentum; 
      aMomentum[0] = pMag*sintheta*cos(phi);
      aMomentum[1] = pMag*sintheta*sin(phi);
      aMomentum[2] = pMag*costheta;  
      //use ParticleGun to generate event 
      particleGun->SetParticleMomentum(aMomentum);
      particleGun->GeneratePrimaryVertex(anEvent);
    }
  }
  else if(generatorName=="cosmic")
  { 
    G4cout<<"generatorName: "<<generatorName<<G4endl;
    G4ParticleTable* particleTable = G4ParticleTable::GetParticleTable();
    G4ParticleDefinition* particle
      = particleTable->FindParticle(particleName);
    G4cout<<"particleName: "<<particleName<<G4endl;
 
    particleGun->SetParticleDefinition(particle);
    particleGun->SetParticlePosition(G4ThreeVector(posX,posY,posZ));
    
    /*std::string path = getenv("GENSIMROOT");
    G4cout<<"path: "<<path<<G4endl;
 
    path += "/root/";
    std::string pFile = path + "ppdc.root";
    std::string thetaFile = path + "theta.root";
    std::string phiFile = path +"phi.root";*/
 
    //TFile* f1 = new TFile(pFile.c_str());
    //H1F* h1 = (TH1F*)f1->Get("htemp");
    G4double pMag = h1->GetRandom()*GeV;
    G4cout<<"pMag: "<<pMag<<G4endl;
    //f1->Close();
 
    //TFile* f2 = new TFile(thetaFile.c_str());
    //TH1F* h2 = (TH1F*)f2->Get("htemp");
    //randomize cos(theta)
    G4double theta =(Double_t)h2->GetRandom();
    G4cout<<"theta: "<<theta<<G4endl;
    G4double costheta = cos(theta);
    //f2->Close();
 
    //TFile* f3 = new TFile(phiFile.c_str());
    //TH1F* h3 = (TH1F*)f3->Get("htemp");
    //randomize phi
    G4double phi = (Double_t)h3->GetRandom();
    G4cout<<"phi: "<<phi<<G4endl;
    //f3->Close();
 
    //f1->Close();
    //f2->Close();
    //f3->Close();
    //ftest->ReOpen("update");
    //tuple->Fill(pMag,theta,phi);
    //counter++;
    //if(counter==2000)
    //  tuple->Write();
 
    G4double sintheta = sqrt(1.-costheta*costheta);
    //computer 3-vector momentum
    G4ParticleMomentum aMomentum; 
    aMomentum[0] = pMag*sintheta*cos(phi);
    aMomentum[1] = pMag*sintheta*sin(phi);
    aMomentum[2] = pMag*costheta;  
    //use ParticleGun to generate event 
    particleGun->SetParticleMomentum(aMomentum);
    particleGun->GeneratePrimaryVertex(anEvent);
  }
 
  else if(generatorName=="genbes")
  {
    G4cout<<"genbes called"<<G4endl;
    if(isGenbes)
    {
      isGenbes=false;
      HEPEvt=new G4HEPEvtInterface(genbesName);
    }
    HEPEvt->GeneratePrimaryVertex(anEvent);
  } 
}

GeneratePrimaries() [2/2]

void BesPrimaryGeneratorAction::GeneratePrimaries

(

G4Event *

anEvent

)

GetGenbesName() [1/2]

G4String BesPrimaryGeneratorAction::GetGenbesName ( )

inline

Definition at line 39 of file InstallArea/include/GenSim/GenSim/BesPrimaryGeneratorAction.hh.

{return genbesName;}

Referenced by BesPrimaryGeneratorMessenger::GetCurrentValue().

GetGenbesName() [2/2]

G4String BesPrimaryGeneratorAction::GetGenbesName ( )

inline

Definition at line 39 of file Simulation/BOOST/GenSim/GenSim-00-00-07/GenSim/BesPrimaryGeneratorAction.hh.

{return genbesName;}

GetGeneratorName() [1/2]

G4String BesPrimaryGeneratorAction::GetGeneratorName ( )

inline

Definition at line 37 of file InstallArea/include/GenSim/GenSim/BesPrimaryGeneratorAction.hh.

{return generatorName;}

Referenced by BesPrimaryGeneratorMessenger::GetCurrentValue().

GetGeneratorName() [2/2]

G4String BesPrimaryGeneratorAction::GetGeneratorName ( )

inline

Definition at line 37 of file Simulation/BOOST/GenSim/GenSim-00-00-07/GenSim/BesPrimaryGeneratorAction.hh.

{return generatorName;}

SetDeltaP() [1/2]

void BesPrimaryGeneratorAction::SetDeltaP ( G4double  dp )

inline

Definition at line 48 of file InstallArea/include/GenSim/GenSim/BesPrimaryGeneratorAction.hh.

{deltaP = dp;}

Referenced by BesPrimaryGeneratorMessenger::SetNewValue().

SetDeltaP() [2/2]

void BesPrimaryGeneratorAction::SetDeltaP ( G4double  dp )

inline

Definition at line 48 of file Simulation/BOOST/GenSim/GenSim-00-00-07/GenSim/BesPrimaryGeneratorAction.hh.

{deltaP = dp;}

SetGenbesName() [1/2]

void BesPrimaryGeneratorAction::SetGenbesName ( G4String  genbesNew )

inline

Definition at line 40 of file InstallArea/include/GenSim/GenSim/BesPrimaryGeneratorAction.hh.

{genbesName=genbesNew;} 

Referenced by BesPrimaryGeneratorMessenger::SetNewValue().

SetGenbesName() [2/2]

void BesPrimaryGeneratorAction::SetGenbesName ( G4String  genbesNew )

inline

Definition at line 40 of file Simulation/BOOST/GenSim/GenSim-00-00-07/GenSim/BesPrimaryGeneratorAction.hh.

{genbesName=genbesNew;} 

SetGeneratorName() [1/2]

void BesPrimaryGeneratorAction::SetGeneratorName ( G4String  genNew )

inline

Definition at line 38 of file InstallArea/include/GenSim/GenSim/BesPrimaryGeneratorAction.hh.

{generatorName=genNew;}

Referenced by BesPrimaryGeneratorMessenger::SetNewValue().

SetGeneratorName() [2/2]

void BesPrimaryGeneratorAction::SetGeneratorName ( G4String  genNew )

inline

Definition at line 38 of file Simulation/BOOST/GenSim/GenSim-00-00-07/GenSim/BesPrimaryGeneratorAction.hh.

{generatorName=genNew;}

SetMaxCos() [1/2]

void BesPrimaryGeneratorAction::SetMaxCos ( G4double  max1 )

inline

Definition at line 44 of file InstallArea/include/GenSim/GenSim/BesPrimaryGeneratorAction.hh.

{maxCos = max1;}

Referenced by BesPrimaryGeneratorMessenger::SetNewValue().

SetMaxCos() [2/2]

void BesPrimaryGeneratorAction::SetMaxCos ( G4double  max1 )

inline

Definition at line 44 of file Simulation/BOOST/GenSim/GenSim-00-00-07/GenSim/BesPrimaryGeneratorAction.hh.

{maxCos = max1;}

SetMinCos() [1/2]

void BesPrimaryGeneratorAction::SetMinCos ( G4double  min1 )

inline

Definition at line 43 of file InstallArea/include/GenSim/GenSim/BesPrimaryGeneratorAction.hh.

{minCos = min1;}

Referenced by BesPrimaryGeneratorMessenger::SetNewValue().

SetMinCos() [2/2]

void BesPrimaryGeneratorAction::SetMinCos ( G4double  min1 )

inline

Definition at line 43 of file Simulation/BOOST/GenSim/GenSim-00-00-07/GenSim/BesPrimaryGeneratorAction.hh.

{minCos = min1;}

SetMomentum() [1/2]

void BesPrimaryGeneratorAction::SetMomentum ( G4double  pm )

inline

Definition at line 47 of file InstallArea/include/GenSim/GenSim/BesPrimaryGeneratorAction.hh.

{pMomentum = pm;}

Referenced by BesPrimaryGeneratorMessenger::SetNewValue().

SetMomentum() [2/2]

void BesPrimaryGeneratorAction::SetMomentum ( G4double  pm )

inline

Definition at line 47 of file Simulation/BOOST/GenSim/GenSim-00-00-07/GenSim/BesPrimaryGeneratorAction.hh.

{pMomentum = pm;}

SetNParticle() [1/2]

void BesPrimaryGeneratorAction::SetNParticle ( G4int  np )

inline

Definition at line 41 of file InstallArea/include/GenSim/GenSim/BesPrimaryGeneratorAction.hh.

{nParticle = np;}

Referenced by BesPrimaryGeneratorMessenger::SetNewValue().

SetNParticle() [2/2]

void BesPrimaryGeneratorAction::SetNParticle ( G4int  np )

inline

Definition at line 41 of file Simulation/BOOST/GenSim/GenSim-00-00-07/GenSim/BesPrimaryGeneratorAction.hh.

{nParticle = np;}

SetParticleName() [1/2]

void BesPrimaryGeneratorAction::SetParticleName ( G4String  pname )

inline

Definition at line 42 of file InstallArea/include/GenSim/GenSim/BesPrimaryGeneratorAction.hh.

{particleName = pname;}

Referenced by BesPrimaryGeneratorMessenger::SetNewValue().

SetParticleName() [2/2]

void BesPrimaryGeneratorAction::SetParticleName ( G4String  pname )

inline

Definition at line 42 of file Simulation/BOOST/GenSim/GenSim-00-00-07/GenSim/BesPrimaryGeneratorAction.hh.

{particleName = pname;}

SetPhiEnd() [1/2]

void BesPrimaryGeneratorAction::SetPhiEnd ( G4double  pe )

inline

Definition at line 46 of file InstallArea/include/GenSim/GenSim/BesPrimaryGeneratorAction.hh.

{phiEnd = pe;}

Referenced by BesPrimaryGeneratorMessenger::SetNewValue().

SetPhiEnd() [2/2]

void BesPrimaryGeneratorAction::SetPhiEnd ( G4double  pe )

inline

Definition at line 46 of file Simulation/BOOST/GenSim/GenSim-00-00-07/GenSim/BesPrimaryGeneratorAction.hh.

{phiEnd = pe;}

SetPhiStart() [1/2]

void BesPrimaryGeneratorAction::SetPhiStart ( G4double  ps )

inline

Definition at line 45 of file InstallArea/include/GenSim/GenSim/BesPrimaryGeneratorAction.hh.

{phiStart =ps;}

Referenced by BesPrimaryGeneratorMessenger::SetNewValue().

SetPhiStart() [2/2]

void BesPrimaryGeneratorAction::SetPhiStart ( G4double  ps )

inline

Definition at line 45 of file Simulation/BOOST/GenSim/GenSim-00-00-07/GenSim/BesPrimaryGeneratorAction.hh.

{phiStart =ps;}

SetPosX() [1/2]

void BesPrimaryGeneratorAction::SetPosX ( G4double  x )

inline

Definition at line 49 of file InstallArea/include/GenSim/GenSim/BesPrimaryGeneratorAction.hh.

{posX = x;}

Referenced by BesPrimaryGeneratorMessenger::SetNewValue().

SetPosX() [2/2]

void BesPrimaryGeneratorAction::SetPosX ( G4double  x )

inline

Definition at line 49 of file Simulation/BOOST/GenSim/GenSim-00-00-07/GenSim/BesPrimaryGeneratorAction.hh.

{posX = x;}

SetPosY() [1/2]

void BesPrimaryGeneratorAction::SetPosY ( G4double  y )

inline

Definition at line 50 of file InstallArea/include/GenSim/GenSim/BesPrimaryGeneratorAction.hh.

{posY = y;}

Referenced by BesPrimaryGeneratorMessenger::SetNewValue().

SetPosY() [2/2]

void BesPrimaryGeneratorAction::SetPosY ( G4double  y )

inline

Definition at line 50 of file Simulation/BOOST/GenSim/GenSim-00-00-07/GenSim/BesPrimaryGeneratorAction.hh.

{posY = y;}

SetPosZ() [1/2]

void BesPrimaryGeneratorAction::SetPosZ ( G4double  z )

inline

Definition at line 51 of file InstallArea/include/GenSim/GenSim/BesPrimaryGeneratorAction.hh.

{posZ = z;}

Referenced by BesPrimaryGeneratorMessenger::SetNewValue().

SetPosZ() [2/2]

void BesPrimaryGeneratorAction::SetPosZ ( G4double  z )

inline

Definition at line 51 of file Simulation/BOOST/GenSim/GenSim-00-00-07/GenSim/BesPrimaryGeneratorAction.hh.

{posZ = z;}

---

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

• InstallArea/include/GenSim/GenSim/BesPrimaryGeneratorAction.hh
• Simulation/BOOST/GenSim/GenSim-00-00-07/GenSim/BesPrimaryGeneratorAction.hh
• BesPrimaryGeneratorAction.cc