#include <CMSHadronPhysicsFTFP_BERT.h>

Inheritance diagram for CMSHadronPhysicsFTFP_BERT:

Public Member Functions
	CMSHadronPhysicsFTFP_BERT (G4int verbose=1)

void	ConstructParticle () override

void	ConstructProcess () override

void	TerminateWorker () override

	~CMSHadronPhysicsFTFP_BERT () override

Protected Member Functions
virtual void	CreateModels ()

virtual void	DumpBanner ()

virtual void	ExtraConfiguration ()

virtual void	Kaon ()

virtual void	Neutron ()

virtual void	Others ()

virtual void	Pion ()

virtual void	Proton ()

Protected Attributes
G4double	maxBERT_kaon

G4double	maxBERT_neutron

G4double	maxBERT_pion

G4double	maxBERT_proton

G4double	minFTFP_kaon

G4double	minFTFP_neutron

G4double	minFTFP_pion

G4double	minFTFP_proton

G4bool	QuasiElastic

G4VectorCache< G4VCrossSectionDataSet * >	xs_ds

G4Cache< G4ComponentGGHadronNucleusXsc * >	xs_k

Detailed Description

Definition at line 24 of file CMSHadronPhysicsFTFP_BERT.h.

Constructor & Destructor Documentation

CMSHadronPhysicsFTFP_BERT::CMSHadronPhysicsFTFP_BERT ( G4int verbose = 1 )

explicit

Definition at line 49 of file CMSHadronPhysicsFTFP_BERT.cc.

References GeV, maxBERT_kaon, maxBERT_neutron, maxBERT_pion, maxBERT_proton, minFTFP_kaon, minFTFP_neutron, minFTFP_pion, and minFTFP_proton.

     : G4VPhysicsConstructor("hInelastic CMS FTFP_BERT"), QuasiElastic(false) {
   minFTFP_pion = 4.0 * GeV;
   maxBERT_pion = 5.0 * GeV;
   minFTFP_kaon = 4.0 * GeV;
   maxBERT_kaon = 5.0 * GeV;
   minFTFP_proton = 4.0 * GeV;
   maxBERT_proton = 5.0 * GeV;
   minFTFP_neutron = 4.0 * GeV;
   maxBERT_neutron = 5.0 * GeV;
 }

CMSHadronPhysicsFTFP_BERT::~CMSHadronPhysicsFTFP_BERT ( )

override

Definition at line 61 of file CMSHadronPhysicsFTFP_BERT.cc.

References xs_ds, and xs_k.

                                                       {
   //Detele master-owned stuff
   delete xs_k.Get();
   std::for_each(xs_ds.Begin(), xs_ds.End(), [](G4VCrossSectionDataSet* el) { delete el; });
 }

Member Function Documentation

void CMSHadronPhysicsFTFP_BERT::ConstructParticle ( )

override

Definition at line 67 of file CMSHadronPhysicsFTFP_BERT.cc.

                                                   {
   G4MesonConstructor pMesonConstructor;
   pMesonConstructor.ConstructParticle();
 
   G4BaryonConstructor pBaryonConstructor;
   pBaryonConstructor.ConstructParticle();
 
   G4ShortLivedConstructor pShortLivedConstructor;
   pShortLivedConstructor.ConstructParticle();
 }

void CMSHadronPhysicsFTFP_BERT::ConstructProcess ( )

override

Definition at line 178 of file CMSHadronPhysicsFTFP_BERT.cc.

References CreateModels(), DumpBanner(), and ExtraConfiguration().

                                                  {
   if (G4Threading::IsMasterThread()) {
     DumpBanner();
   }
   CreateModels();
   ExtraConfiguration();
 }

void CMSHadronPhysicsFTFP_BERT::CreateModels ( )

protectedvirtual

Definition at line 93 of file CMSHadronPhysicsFTFP_BERT.cc.

References Kaon(), Neutron(), Others(), Pion(), and Proton().

Referenced by ConstructProcess().

                                              {
   Neutron();
   Proton();
   Pion();
   Kaon();
   Others();
 }

void CMSHadronPhysicsFTFP_BERT::DumpBanner ( )

protectedvirtual

Definition at line 85 of file CMSHadronPhysicsFTFP_BERT.cc.

References ecalTB2006H4_GenSimDigiReco_cfg::G4cout, GeV, maxBERT_kaon, maxBERT_neutron, maxBERT_pion, maxBERT_proton, minFTFP_kaon, minFTFP_neutron, minFTFP_pion, and minFTFP_proton.

Referenced by ConstructProcess().

                                            {
   G4cout << G4endl << " FTFP_BERT : new threshold between BERT and FTFP is over the interval " << G4endl
          << " for pions :   " << minFTFP_pion / GeV << " to " << maxBERT_pion / GeV << " GeV" << G4endl
          << " for kaons :   " << minFTFP_kaon / GeV << " to " << maxBERT_kaon / GeV << " GeV" << G4endl
          << " for proton :  " << minFTFP_proton / GeV << " to " << maxBERT_proton / GeV << " GeV" << G4endl
          << " for neutron : " << minFTFP_neutron / GeV << " to " << maxBERT_neutron / GeV << " GeV" << G4endl << G4endl;
 }

void CMSHadronPhysicsFTFP_BERT::ExtraConfiguration ( )

protectedvirtual

Definition at line 186 of file CMSHadronPhysicsFTFP_BERT.cc.

References ztee::capture(), mps_fire::i, MetAnalyzer::pv(), xs_ds, and xs_k.

Referenced by ConstructProcess().

                                                    {
   //Modify XS for kaons
   auto xsk = new G4ComponentGGHadronNucleusXsc();
   xs_k.Put(xsk);
   G4VCrossSectionDataSet* kaonxs = new G4CrossSectionInelastic(xsk);
   xs_ds.Push_back(kaonxs);
   G4PhysListUtil::FindInelasticProcess(G4KaonMinus::KaonMinus())->AddDataSet(kaonxs);
   G4PhysListUtil::FindInelasticProcess(G4KaonPlus::KaonPlus())->AddDataSet(kaonxs);
   G4PhysListUtil::FindInelasticProcess(G4KaonZeroShort::KaonZeroShort())->AddDataSet(kaonxs);
   G4PhysListUtil::FindInelasticProcess(G4KaonZeroLong::KaonZeroLong())->AddDataSet(kaonxs);
 
   //Modify Neutrons
   auto xs_n_in = (G4NeutronInelasticXS*)G4CrossSectionDataSetRegistry::Instance()->GetCrossSectionDataSet(
       G4NeutronInelasticXS::Default_Name());
   xs_ds.Push_back(xs_n_in);  //TODO: Is this needed? Who owns the pointer?
   G4PhysListUtil::FindInelasticProcess(G4Neutron::Neutron())->AddDataSet(xs_n_in);
   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);
   }
   auto xs_n_c = (G4NeutronCaptureXS*)G4CrossSectionDataSetRegistry::Instance()->GetCrossSectionDataSet(
       G4NeutronCaptureXS::Default_Name());
   xs_ds.Push_back(xs_n_c);
   capture->AddDataSet(xs_n_c);
   capture->RegisterMe(new G4NeutronRadCapture());
 }

void CMSHadronPhysicsFTFP_BERT::Kaon ( )

protectedvirtual

Definition at line 149 of file CMSHadronPhysicsFTFP_BERT.cc.

References gen::k, maxBERT_kaon, minFTFP_kaon, and QuasiElastic.

Referenced by CreateModels().

                                      {
   auto k = new G4KaonBuilder;
   AddBuilder(k);
   auto ftfpk = new G4FTFPKaonBuilder(QuasiElastic);
   AddBuilder(ftfpk);
   k->RegisterMe(ftfpk);
   ftfpk->SetMinEnergy(minFTFP_kaon);
   auto bertk = new G4BertiniKaonBuilder;
   AddBuilder(bertk);
   k->RegisterMe(bertk);
   bertk->SetMaxEnergy(maxBERT_kaon);
   k->Build();
 }

void CMSHadronPhysicsFTFP_BERT::Neutron ( )

protectedvirtual

Definition at line 101 of file CMSHadronPhysicsFTFP_BERT.cc.

References GeV, maxBERT_neutron, minFTFP_neutron, and QuasiElastic.

Referenced by CreateModels().

                                         {
   //General schema:
   // 1) Create a builder
   // 2) Call AddBuilder
   // 3) Configure the builder, possibly with sub-builders
   // 4) Call builder->Build()
   auto neu = new G4NeutronBuilder;
   AddBuilder(neu);
   auto ftfpn = new G4FTFPNeutronBuilder(QuasiElastic);
   AddBuilder(ftfpn);
   neu->RegisterMe(ftfpn);
   ftfpn->SetMinEnergy(minFTFP_neutron);
   auto bertn = new G4BertiniNeutronBuilder;
   AddBuilder(bertn);
   neu->RegisterMe(bertn);
   bertn->SetMinEnergy(0. * GeV);
   bertn->SetMaxEnergy(maxBERT_neutron);
   neu->Build();
 }

void CMSHadronPhysicsFTFP_BERT::Others ( )

protectedvirtual

===== Anti-barions==== //

Definition at line 163 of file CMSHadronPhysicsFTFP_BERT.cc.

References QuasiElastic.

Referenced by CreateModels().

                                        {
   //===== Hyperons ====== //
   auto hyp = new G4HyperonFTFPBuilder;
   AddBuilder(hyp);
   hyp->Build();
 
   auto abar = new G4AntiBarionBuilder;
   AddBuilder(abar);
   auto ftfpabar = new G4FTFPAntiBarionBuilder(QuasiElastic);
   AddBuilder(ftfpabar);
   abar->RegisterMe(ftfpabar);
   abar->Build();
 }

void CMSHadronPhysicsFTFP_BERT::Pion ( )

protectedvirtual

Definition at line 135 of file CMSHadronPhysicsFTFP_BERT.cc.

References maxBERT_pion, minFTFP_pion, pi, and QuasiElastic.

Referenced by CreateModels().

                                      {
   auto pi = new G4PionBuilder;
   AddBuilder(pi);
   auto ftfppi = new G4FTFPPionBuilder(QuasiElastic);
   AddBuilder(ftfppi);
   pi->RegisterMe(ftfppi);
   ftfppi->SetMinEnergy(minFTFP_pion);
   auto bertpi = new G4BertiniPionBuilder;
   AddBuilder(bertpi);
   pi->RegisterMe(bertpi);
   bertpi->SetMaxEnergy(maxBERT_pion);
   pi->Build();
 }

void CMSHadronPhysicsFTFP_BERT::Proton ( )

protectedvirtual

Definition at line 121 of file CMSHadronPhysicsFTFP_BERT.cc.

References maxBERT_proton, minFTFP_proton, and QuasiElastic.

Referenced by CreateModels().

                                        {
   auto pro = new G4ProtonBuilder;
   AddBuilder(pro);
   auto ftfpp = new G4FTFPProtonBuilder(QuasiElastic);
   AddBuilder(ftfpp);
   pro->RegisterMe(ftfpp);
   ftfpp->SetMinEnergy(minFTFP_proton);
   auto bertp = new G4BertiniProtonBuilder;
   AddBuilder(bertp);
   pro->RegisterMe(bertp);
   bertp->SetMaxEnergy(maxBERT_proton);
   pro->Build();
 }

void CMSHadronPhysicsFTFP_BERT::TerminateWorker ( )

override

Definition at line 78 of file CMSHadronPhysicsFTFP_BERT.cc.

References xs_ds, and xs_k.

                                                 {
   delete xs_k.Get();
   std::for_each(xs_ds.Begin(), xs_ds.End(), [](G4VCrossSectionDataSet* el) { delete el; });
   xs_ds.Clear();
   G4VPhysicsConstructor::TerminateWorker();
 }