HadronPhysicsQGSPCMS_FTFP_BERT Class Reference

#include <HadronPhysicsQGSPCMS_FTFP_BERT.h>

Inheritance diagram for HadronPhysicsQGSPCMS_FTFP_BERT:

Classes
struct	ThreadPrivate

Public Member Functions
void	ConstructParticle () override
void	ConstructProcess () override
	HadronPhysicsQGSPCMS_FTFP_BERT (G4int verbose)
	~HadronPhysicsQGSPCMS_FTFP_BERT () override

Private Member Functions
void	CreateModels ()

Static Private Attributes
static G4ThreadLocal ThreadPrivate *	tpdata = nullptr

Detailed Description

Definition at line 28 of file HadronPhysicsQGSPCMS_FTFP_BERT.h.

Constructor & Destructor Documentation

HadronPhysicsQGSPCMS_FTFP_BERT::HadronPhysicsQGSPCMS_FTFP_BERT

(

G4int

verbose

)

Definition at line 25 of file HadronPhysicsQGSPCMS_FTFP_BERT.cc.

  :  G4VPhysicsConstructor("hInelasticQGSPCMS_FTFP_BERT")
{}

HadronPhysicsQGSPCMS_FTFP_BERT::~HadronPhysicsQGSPCMS_FTFP_BERT ( )

override

Definition at line 95 of file HadronPhysicsQGSPCMS_FTFP_BERT.cc.

References HadronPhysicsQGSPCMS_FTFP_BERT::ThreadPrivate::theAntiBaryon, HadronPhysicsQGSPCMS_FTFP_BERT::ThreadPrivate::theBertiniNeutron, HadronPhysicsQGSPCMS_FTFP_BERT::ThreadPrivate::theBertiniPiK, HadronPhysicsQGSPCMS_FTFP_BERT::ThreadPrivate::theBertiniPro, HadronPhysicsQGSPCMS_FTFP_BERT::ThreadPrivate::theFTFPAntiBaryon, HadronPhysicsQGSPCMS_FTFP_BERT::ThreadPrivate::theFTFPNeutron, HadronPhysicsQGSPCMS_FTFP_BERT::ThreadPrivate::theFTFPPiK, HadronPhysicsQGSPCMS_FTFP_BERT::ThreadPrivate::theFTFPPro, HadronPhysicsQGSPCMS_FTFP_BERT::ThreadPrivate::theHyperon, HadronPhysicsQGSPCMS_FTFP_BERT::ThreadPrivate::theNeutrons, HadronPhysicsQGSPCMS_FTFP_BERT::ThreadPrivate::thePiK, HadronPhysicsQGSPCMS_FTFP_BERT::ThreadPrivate::thePro, HadronPhysicsQGSPCMS_FTFP_BERT::ThreadPrivate::theQGSPNeutron, HadronPhysicsQGSPCMS_FTFP_BERT::ThreadPrivate::theQGSPPiK, HadronPhysicsQGSPCMS_FTFP_BERT::ThreadPrivate::theQGSPPro, and tpdata.

{
  if(nullptr != tpdata) {
   delete tpdata->theQGSPNeutron;
   delete tpdata->theFTFPNeutron;
   delete tpdata->theBertiniNeutron;
   delete tpdata->theNeutrons;

   delete tpdata->theQGSPPro;
   delete tpdata->theFTFPPro;
   delete tpdata->thePro;
   delete tpdata->theBertiniPro;

   delete tpdata->theQGSPPiK;
   delete tpdata->theFTFPPiK;
   delete tpdata->theBertiniPiK;
   delete tpdata->thePiK;

   delete tpdata->theHyperon;
   delete tpdata->theAntiBaryon;
   delete tpdata->theFTFPAntiBaryon;
   
   delete tpdata;
   tpdata = nullptr;
  }
}

Member Function Documentation

void HadronPhysicsQGSPCMS_FTFP_BERT::ConstructParticle ( )

override

Definition at line 122 of file HadronPhysicsQGSPCMS_FTFP_BERT.cc.

{
  G4MesonConstructor pMesonConstructor;
  pMesonConstructor.ConstructParticle();

  G4BaryonConstructor pBaryonConstructor;
  pBaryonConstructor.ConstructParticle();

  G4ShortLivedConstructor pShortLivedConstructor;
  pShortLivedConstructor.ConstructParticle();  

  G4IonConstructor pIonConstructor;
  pIonConstructor.ConstructParticle();
}

void HadronPhysicsQGSPCMS_FTFP_BERT::ConstructProcess ( )

override

Definition at line 138 of file HadronPhysicsQGSPCMS_FTFP_BERT.cc.

References ztee::capture(), CreateModels(), mps_fire::i, MetAnalyzer::pv(), HadronPhysicsQGSPCMS_FTFP_BERT::ThreadPrivate::theAntiBaryon, HadronPhysicsQGSPCMS_FTFP_BERT::ThreadPrivate::theHyperon, HadronPhysicsQGSPCMS_FTFP_BERT::ThreadPrivate::theNeutrons, HadronPhysicsQGSPCMS_FTFP_BERT::ThreadPrivate::thePiK, HadronPhysicsQGSPCMS_FTFP_BERT::ThreadPrivate::thePro, and tpdata.

{
  if ( tpdata == nullptr ) { tpdata = new ThreadPrivate; }
  CreateModels();
  tpdata->theNeutrons->Build();
  tpdata->thePro->Build();
  tpdata->thePiK->Build();
  tpdata->theHyperon->Build(); 
  tpdata->theAntiBaryon->Build(); 

  // --- Neutrons ---
  G4PhysListUtil::FindInelasticProcess(G4Neutron::Neutron())
    ->AddDataSet(new G4NeutronInelasticXS());

  G4HadronicProcess* capture = nullptr;
  G4ProcessManager* pmanager = G4Neutron::Neutron()->GetProcessManager();
  G4ProcessVector*  pv = pmanager->GetProcessList();
  for ( size_t i=0; i < static_cast<size_t>(pv->size()); ++i ) {
    if ( fCapture == ((*pv)[i])->GetProcessSubType() ) {
      capture = static_cast<G4HadronicProcess*>((*pv)[i]);
    }
  }
  if ( ! capture ) {
    capture = new G4HadronCaptureProcess("nCapture");
    pmanager->AddDiscreteProcess(capture);
  }
  capture->AddDataSet(new G4NeutronCaptureXS());
  capture->RegisterMe(new G4NeutronRadCapture());
}

void HadronPhysicsQGSPCMS_FTFP_BERT::CreateModels ( )

private

Definition at line 29 of file HadronPhysicsQGSPCMS_FTFP_BERT.cc.

References ecalTB2006H4_GenSimDigiReco_cfg::G4cout, GeV, HadronPhysicsQGSPCMS_FTFP_BERT::ThreadPrivate::theAntiBaryon, HadronPhysicsQGSPCMS_FTFP_BERT::ThreadPrivate::theBertiniNeutron, HadronPhysicsQGSPCMS_FTFP_BERT::ThreadPrivate::theBertiniPiK, HadronPhysicsQGSPCMS_FTFP_BERT::ThreadPrivate::theBertiniPro, HadronPhysicsQGSPCMS_FTFP_BERT::ThreadPrivate::theFTFPAntiBaryon, HadronPhysicsQGSPCMS_FTFP_BERT::ThreadPrivate::theFTFPNeutron, HadronPhysicsQGSPCMS_FTFP_BERT::ThreadPrivate::theFTFPPiK, HadronPhysicsQGSPCMS_FTFP_BERT::ThreadPrivate::theFTFPPro, HadronPhysicsQGSPCMS_FTFP_BERT::ThreadPrivate::theHyperon, HadronPhysicsQGSPCMS_FTFP_BERT::ThreadPrivate::theNeutrons, HadronPhysicsQGSPCMS_FTFP_BERT::ThreadPrivate::thePiK, HadronPhysicsQGSPCMS_FTFP_BERT::ThreadPrivate::thePro, HadronPhysicsQGSPCMS_FTFP_BERT::ThreadPrivate::theQGSPNeutron, HadronPhysicsQGSPCMS_FTFP_BERT::ThreadPrivate::theQGSPPiK, HadronPhysicsQGSPCMS_FTFP_BERT::ThreadPrivate::theQGSPPro, and tpdata.

Referenced by ConstructProcess().

{
  // First transition, between BERT and FTF/P
  G4double minFTFP= 6.0 * GeV;  
  G4double maxBERT= 8.0 * GeV;  
  // Second transition, between FTF/P and QGS/P
  G4double minQGSP= 12.0 * GeV;
  G4double maxFTFP= 25.0 * GeV; 

  G4bool   quasiElasFTF= false;   // Use built-in quasi-elastic (not add-on)
  G4bool   quasiElasQGS= true;    // For QGS, it must use it.

  G4cout << " New QGSPCMS_FTFP_BERT hadronic inealstic physics" << G4endl;
  G4cout << "   Thresholds: " << G4endl;
  G4cout << "     1) between BERT  and FTFP over the interval " 
     << minFTFP/GeV << " to " << maxBERT/GeV << " GeV. " << G4endl;
  G4cout << "     2) between FTFP and QGSP over the interval " 
     << minQGSP/GeV << " to " << maxFTFP/GeV << " GeV. " << G4endl;
  G4cout << "   QuasiElastic: " << quasiElasQGS << " for QGS "
     << " and " << quasiElasFTF << " for FTF " << G4endl;

  tpdata->theNeutrons=new G4NeutronBuilder;
  tpdata->theNeutrons->RegisterMe(
    tpdata->theQGSPNeutron=new G4QGSPNeutronBuilder(quasiElasQGS));
  tpdata->theQGSPNeutron->SetMinEnergy(minQGSP);   
  tpdata->theNeutrons->RegisterMe(
    tpdata->theFTFPNeutron=new G4FTFPNeutronBuilder(quasiElasFTF));
  tpdata->theFTFPNeutron->SetMinEnergy(minFTFP);   
  tpdata->theFTFPNeutron->SetMaxEnergy(maxFTFP);   

  tpdata->theNeutrons->RegisterMe(
    tpdata->theBertiniNeutron=new G4BertiniNeutronBuilder);
  tpdata->theBertiniNeutron->SetMinEnergy(0.0*GeV);
  tpdata->theBertiniNeutron->SetMaxEnergy(maxBERT);  

  tpdata->thePro=new G4ProtonBuilder;
  tpdata->thePro->RegisterMe(
    tpdata->theQGSPPro=new G4QGSPProtonBuilder(quasiElasQGS));
  tpdata->theQGSPPro->SetMinEnergy(minQGSP);   
  tpdata->thePro->RegisterMe(
    tpdata->theFTFPPro=new G4FTFPProtonBuilder(quasiElasFTF));
  tpdata->theFTFPPro->SetMinEnergy(minFTFP);   
  tpdata->theFTFPPro->SetMaxEnergy(maxFTFP);   
  tpdata->thePro->RegisterMe(
    tpdata->theBertiniPro=new G4BertiniProtonBuilder);
  tpdata->theBertiniPro->SetMaxEnergy(maxBERT); 
  
  tpdata->thePiK=new G4PiKBuilder;
  tpdata->thePiK->RegisterMe(
    tpdata->theQGSPPiK=new G4QGSPPiKBuilder(quasiElasQGS));
  tpdata->theQGSPPiK->SetMinEnergy(minQGSP);   
  tpdata->thePiK->RegisterMe(
    tpdata->theFTFPPiK=new G4FTFPPiKBuilder(quasiElasFTF));
  tpdata->theFTFPPiK->SetMaxEnergy(maxFTFP);  
  tpdata->theFTFPPiK->SetMinEnergy(minFTFP); 
  tpdata->thePiK->RegisterMe(tpdata->theBertiniPiK=new G4BertiniPiKBuilder);
  tpdata->theBertiniPiK->SetMaxEnergy(maxBERT); 
  
  // Hyperons use FTF
  tpdata->theHyperon=new G4HyperonFTFPBuilder;

  tpdata->theAntiBaryon=new G4AntiBarionBuilder;
  tpdata->theAntiBaryon->RegisterMe(
    tpdata->theFTFPAntiBaryon=new G4FTFPAntiBarionBuilder(quasiElasFTF));
}

Member Data Documentation

G4ThreadLocal HadronPhysicsQGSPCMS_FTFP_BERT::ThreadPrivate * HadronPhysicsQGSPCMS_FTFP_BERT::tpdata = nullptr

staticprivate

Definition at line 61 of file HadronPhysicsQGSPCMS_FTFP_BERT.h.

Referenced by ConstructProcess(), CreateModels(), and ~HadronPhysicsQGSPCMS_FTFP_BERT().