How to setting up NaI scintillator

[repan]

Hi all,
I defined a NaI scintillation detectors. It can generate photons but all these photons have wave length 200 nm. May I ask what’s the problem. My codes are:

    G4double energy[2] = {1.239841939*eV/0.2, 1.239841939*eV/0.9}; 
    G4double rindexNaI[2] = {1.78, 1.78};
    G4double fraction[2] = {1.0, 1.0};
    G4MaterialPropertiesTable *mptNaI=new G4MaterialPropertiesTable();
    mptNaI->AddProperty("RINDEX",energy,rindexNaI,2);
    mptNaI->AddProperty("SCINTILLATIONCOMPONENT1",energy,fraction,2);
    mptNaI->AddConstProperty("SCINTILLATIONYIELD",38.0/keV);
    mptNaI->AddConstProperty("RESOLUTIONSCALE",1.0);
    mptNaI->AddConstProperty("SCINTILLATIONTIMECONSTANT1",250*ns);
    mptNaI->AddConstProperty("SCINTILLATIONYIELD1",1.);
    NaI->SetMaterialPropertiesTable(mptNaI);

[repan]

setting up G4double energy[2] = {1.239841939eV/0.9, 1.239841939eV/0.2};

[phirippu]

Just a note, Saint-Gobain lists 1.85 as the refractive index.

[KACheminant]

Just a little tip:

You can also use vectors instead of arrays to define your energy, fraction and refractive index to look like that:

std::vector<G4double> energy = {1.239841939*eV/0.9, 1.239841939*eV/0.2}

One perk of doing so is that you don’t have to specify the number of entries in your AddProperty for the RINDEX and SCINTILLATIONCOMPONENT1