# The energy of the gamma photons generated is very high

Dear experts,
I am working on a radiator, it functions as follows
It absorbs the incoming electrons and generates transition radiation photons (or gamma photons). But I am facing an issue, the energy of the gamma photons generated is very high, I am firing electrons of energy 3GeV and the energy of gamma photons generated is in the MeV to GeV range which in reality should be in 0-200 keV range. I have tried different EM modular physics list still the result remains the same, I also tried to set lower energy and higher energy limit for G4eBremsstrahlung, but still, the energy of gamma photons remains very high
I have attached an image of the output.

######################################################
G4cout <<"G4Material* CH2 = "<<CH2<<G4endl;
G4cout <<"G4Material* Mylar = "<<Mylar<<G4endl;

G4double foilDensity = 0.91g/cm3; // CH2 1.39g/cm3; // Mylar // 0.534g/cm3; //Li
G4double gasDensity = 1.2928
mg/cm3; // Air // 1.977mg/cm3; // CO2 0.178mg/cm3; // He

fGasGap = 0.600
mm;
fDetGap = 0.01*mm ;

G4cout<<"foil gas ratio is :- "<<foilGasRatio<<G4endl;

fRadThick = .3*cm - fGasGap + fDetGap;
G4cout<<"The foil number is "<<fFoilNumber;

G4double totDensity = foilDensityfoilGasRatio + gasDensity(1.0-foilGasRatio) ;

G4cout<<"tot density is :- "<<totDensity<<G4endl;
G4cout<<"Rad totDensity = “<<totDensity/(g/cm3)<<” g/cm3 "<<G4endl ;
G4cout<<"use for Rad totDensity = “<<totDensity/(g/cm3)<<” g/cm3 "<<G4endl ;

G4double fractionFoil = foilDensityfoilGasRatio/totDensity ;
G4double fractionGas = gasDensity
(1.0-foilGasRatio)/totDensity ;
G4cout<<"Fraction gas "<<fractionGas<<G4endl;
G4cout<<"Fraction foil "<<fractionFoil<<G4endl;

G4double NewDensity = 0.083*(g/cm3);
G4cout<<"new Rad with totDensity = “<<NewDensity/(g/cm3)<<” g/cm3 "<<G4endl ;

// default materials of the detector and TR radiator
fFoilMat = CH2; // Kapton; // Mylar ; // Li ; // CH2 ;
fGasMat = Air; // CO2; // He; //

G4Box* solidEnv =
new G4Box(“Envelope”, //its name
7.5cm,7.5cm,7.5*cm); //its size

new G4LogicalVolume(solidEnv, //its solid
“Envelope”); //its name

new G4PVPlacement(0, //no rotation
G4ThreeVector(0.,0.,7.5*cm), //at (0,0,0)
“Envelope”, //its name
logicWorld, //its mother volume
false, //no boolean operation
1, //copy number
checkOverlaps);
##################################################
Physics list – (Partial code)
else if (particleName == “e-”)
{
// Construct processes for electron
theeminusStepCut = new VtxStepCut();
theeminusStepCut->SetMaxStep(step_cut) ;

``````  G4eIonisation* eioni = new G4eIonisation();
G4PAIModel*     pai = new G4PAIModel(particle,"PAIModel");

/*
G4eBremsstrahlung* msc = new G4eBremsstrahlung;
G4SeltzerBergerModel* msc1 = new G4SeltzerBergerModel();
G4eBremsstrahlungRelModel* msc2 = new G4eBremsstrahlungRelModel();
msc1->SetHighEnergyLimit(highEnergyLimit);
msc1->SetSecondaryThreshold(100*eV);
msc2->SetHighEnergyLimit(highEnergyLimit);
msc2->SetSecondaryThreshold(100*eV);

}
``````

If someone wants to run the complete code. The following is the GitHub repository link :-

Kindly help me out with this, and thanks in advance.

##################################################
construction.cc (18.6 KB)
physics.cc (16.6 KB)

Hello,