Message: Re: No Optical Photon is produced Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None Re: No Optical Photon is produced 

Forum: Processes Involving Optical Photons
Re: None No Optical Photon is produced (Davinder Siwal)
Re: Feedback Re: No Optical Photon is produced (Gumplinger Peter)
Date: 15 Sep, 2010
From: Davinder Siwal <Davinder Siwal>

Hiiii

I am Using the following physics list for the optical process

NdetOpticalPhysics::NdetOpticalPhysics(const G4String& name)
  : G4VPhysicsConstructor(name)
{
}

NdetOpticalPhysics::~NdetOpticalPhysics() {}

void NdetOpticalPhysics::ConstructProcess()
{
  ConstructOpticalProcess();
}
void NdetOpticalPhysics::ConstructOpticalProcess()
{
  theScintillationProcess = new G4Scintillation("Scintillation");
  theAbsorptionProcess = new G4OpAbsorption();
  theRayleighScatteringProcess = new G4OpRayleigh();
  theBoundaryProcess  = new G4OpBoundaryProcess();

  //-------- Set Verbosity 
  theScintillationProcess->SetVerboseLevel(1);
  theAbsorptionProcess->SetVerboseLevel(1);
  theRayleighScatteringProcess->SetVerboseLevel(1);
  theBoundaryProcess->SetVerboseLevel(1);
  //-----------------------

  G4ProcessManager* pManager = 0;
  pManager = G4OpticalPhoton::OpticalPhoton()->GetProcessManager();
  pManager->AddDiscreteProcess(theAbsorptionProcess);
  pManager->AddDiscreteProcess(theRayleighScatteringProcess);
  pManager->AddDiscreteProcess(theBoundaryProcess);

  theScintillationProcess->SetScintillationYieldFactor(1);
  theScintillationProcess->SetTrackSecondariesFirst(true);

   G4OpticalSurfaceModel themodel = unified;
  theBoundaryProcess->SetModel(themodel);

  theParticleIterator->reset();
  while( (*theParticleIterator)())
    {
   G4ParticleDefinition* particle = theParticleIterator->value();
   G4ProcessManager* pmanager = particle->GetProcessManager();
   //  G4String particleName = particle->GetParticleName();

   if (theScintillationProcess->IsApplicable(*particle)) 
     {
       pmanager->AddProcess(theScintillationProcess);
       pmanager->SetProcessOrderingToLast(theScintillationProcess, idxAtRest);
       pmanager->SetProcessOrderingToLast(theScintillationProcess, idxPostStep);
     }/*
   if (particleName == "opticalphoton") 
     {
       G4cout << " AddDiscreteProcess to OpticalPhoton " << G4endl;
       pmanager->AddDiscreteProcess(theAbsorptionProcess);
       pmanager->AddDiscreteProcess(theRayleighScatteringProcess);
       pmanager->AddDiscreteProcess(theBoundaryProcess);
     }
      */
    }
}

and using the following scintillation properties for the detector

 const G4int nEntries = 2;
 double PhotonEnergy[nEntries] = {2.47*eV,3.09*eV};
 double RefractiveIndex[nEntries] = {1.538,1.538};
 double AbsLength[nEntries] = {4*m,4*m};

 G4MaterialPropertiesTable* myDet = new G4MaterialPropertiesTable();
 myDet->AddProperty("RINDEX", PhotonEnergy, RefractiveIndex, nEntries);
 //myDet->AddProperty("ABSLENGTH", PhotonEnergy, AbsLength, nEntries);

 myDet->AddConstProperty("SCINTILLATIONYIELD",16000./MeV);
 myDet->AddConstProperty("DECAYTIME",3.3*ns);

 Scinti->SetMaterialPropertiesTable(myDet);

apart from this i don't have any other information about the property

i am not able to figure out where is the fault

thanks in advance

Inline Depth:
 1 1
 All All
Outline Depth:
 1 1
 2 2
 All All
Add message: (add)

1 Feedback: Re: No Optical Photon is produced   (Gumplinger Peter - 15 Sep, 2010)
(_ None: Re: No Optical Photon is produced   (Davinder Siwal - 15 Sep, 2010)
(_ None: Re: No Optical Photon is produced(Solved)   (Davinder Siwal - 20 Sep, 2010)
 Add Message Add Message
to: "Re: No Optical Photon is produced"

 Subscribe Subscribe

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