Message: Re: alphaInelastic products Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None Re: alphaInelastic products 

Forum: Hadronic Processes
Re: Question alphaInelastic products (Coste Benoit)
Date: 23 Feb, 2011
From: Coste Benoit <Coste Benoit>

Edit: I made a mistake. Actually the alphaInelastic process doesn't come from G4IonBinaryCascadePhysics::ConstructProcess() but from a G4HadronInelasticProcess with G4BinaryLightIonReaction registered.

Also, I forgot to mention that I made my test 3 TeV (total energy)

Here is the important piece of code:

          G4int nParticles = 5;
          G4String processName[]={"dInelastic","tInelastic","He3Inelastic","alphaInelastic","ionInelastic"};
          G4String particleName[]={"deuteron","triton","He3","alpha","GenericIon"};
          G4double Emin[]={1.*TeV,1.*TeV,20.*GeV,1.*TeV,20.*GeV};      

          G4ProcessTable* processTable=G4ProcessTable::GetProcessTable();

          for (G4int i=0; i< nParticles; i++) {
            G4BinaryLightIonReaction *theGenIonCascade = new G4BinaryLightIonReaction;
            G4HadronInelasticProcess* hadi = (G4HadronInelasticProcess*)processTable->FindProcess(processName[i],particleName[i]);
            if (hadi != 0) {
              G4cout << "Found " << hadi->GetProcessName()  << " for " << particleName[i] << G4endl;
              hadi->AddDataSet(theSihverCrossSection);
              theGenIonCascade->SetMinEnergy(Emin[i] - MeV);
              theGenIonCascade->SetMaxEnergy(100*TeV);
              hadi->RegisterMe(theGenIonCascade);

            }

 Add Message Add Message
to: "Re: alphaInelastic products"

 Subscribe Subscribe

This site runs SLAC HyperNews version 1.11-slac-98, derived from the original HyperNews


[ Geant 4 Home | Geant 4 HyperNews | Search | Request New Forum | Feedback ]