LCOpticalPhysics.cc

Go to the documentation of this file.

//
//                                                            //
//  Title:  Optical  physics for a Linear Collider Detector      //
//  Date:                                         //
//  Author:                                 //
//                                                            //
//

#include "LCOpticalPhysics.hh"

// #include "globals.hh"
// #include "G4ios.hh"
// #include <iomanip>   

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

LCOpticalPhysics::~LCOpticalPhysics()
{
}



void LCOpticalPhysics::ConstructParticle()
{

//  G4cout << "Inside LCOpticalPhysics::ConstructParticle()" << G4endl;

  // optical photon
  G4OpticalPhoton::OpticalPhotonDefinition();
}

void LCOpticalPhysics::ConstructProcess()
{
  ConstructOp();

}

void LCOpticalPhysics::ConstructOp()
{
  theCerenkovProcess           = new G4Cerenkov("Cerenkov");
  theScintillationProcess      = new G4Scintillation("Scintillation");
  theAbsorptionProcess         = new G4OpAbsorption();
  theRayleighScatteringProcess = new G4OpRayleigh();
  theBoundaryProcess           = new G4OpBoundaryProcess();

//  theCerenkovProcess->DumpPhysicsTable();
//  theScintillationProcess->DumpPhysicsTable();
//  theAbsorptionProcess->DumpPhysicsTable();
//  theRayleighScatteringProcess->DumpPhysicsTable();



//#ifdef debug  
//  G4cout << "VERBOSE = " << theCerenkovProcess->GetVerboseLevel() << G4endl;
//  G4cout << "Inside LCOpticalPhysics::ConstructOp()" << G4endl;
  //#endif
  
  
  theCerenkovProcess->SetMaxNumPhotonsPerStep(300);
  theCerenkovProcess->SetTrackSecondariesFirst(true);
  
  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 (theCerenkovProcess->IsApplicable(*particle)) {
      //#ifdef debug
      //G4cout << "applicable to Cerenkov" << G4endl;
      //#endif
      pmanager->AddContinuousProcess(theCerenkovProcess);
    }
    if (theScintillationProcess->IsApplicable(*particle)) {
      //#ifdef debug
      //G4cout << "applicable to Scintillator" << G4endl;
      //#endif
      pmanager->AddProcess(theScintillationProcess);
      pmanager->SetProcessOrderingToLast(theScintillationProcess, idxAtRest);
      pmanager->SetProcessOrderingToLast(theScintillationProcess, idxPostStep);
    }
    if (particleName == "opticalphoton") {
      //#ifdef debug
      //G4cout << " AddDiscreteProcess to OpticalPhoton " << G4endl;
      //#endif
      pmanager->AddDiscreteProcess(theAbsorptionProcess);
      pmanager->AddDiscreteProcess(theRayleighScatteringProcess);
      pmanager->AddDiscreteProcess(theBoundaryProcess);
    }
  }
}

---