G4XTankFastCerenkov.h

Go to the documentation of this file.

//
// ********************************************************************
// * DISCLAIMER                                                       *
// *                                                                  *
// * The following disclaimer summarizes all the specific disclaimers *
// * of contributors to this software. The specific disclaimers,which *
// * govern, are listed with their locations in:                      *
// *   http://cern.ch/geant4/license                                  *
// *                                                                  *
// * Neither the authors of this software system, nor their employing *
// * institutes,nor the agencies providing financial support for this *
// * work  make  any representation or  warranty, express or implied, *
// * regarding  this  software system or assume any liability for its *
// * use.                                                             *
// *                                                                  *
// * This  code  implementation is the  intellectual property  of the *
// * GEANT4 collaboration.                                            *
// * By copying,  distributing  or modifying the Program (or any work *
// * based  on  the Program)  you indicate  your  acceptance of  this *
// * statement, and all its terms.                                    *
// ********************************************************************

#ifndef _G4XTankSimulatorAG_G4XTankFastCerenkov_h_
#define _G4XTankSimulatorAG_G4XTankFastCerenkov_h_

#include <sdet/PMT.h>
using sdet::PMT;
#include <utl/Point.h>
using utl::Point;
#include <sevt/Station.h>

#include <globals.hh>
#include <templates.hh>
#include <Randomize.hh>
#include <G4ThreeVector.hh>
#include <G4ParticleMomentum.hh>
#include <G4Step.hh>
#include <G4VContinuousProcess.hh>
#include <G4OpticalPhoton.hh>
#include <G4DynamicParticle.hh>
#include <G4Material.hh>
#include <G4PhysicsTable.hh>
#include <G4MaterialPropertiesTable.hh>
#include <G4PhysicsOrderedFreeVector.hh>
#include <G4ios.hh>
#include <G4Poisson.hh>

//#include <G4FastTankSimulatorPS/G4FastTankConstruction.h>

#define IN_TANK 0
#define DOWN_IN_TANK 1
#define IN_DOME_1 11
#define IN_DOME_2 12
#define IN_DOME_3 13
#define IN_INTERFACE_1 21
#define IN_INTERFACE_2 22
#define IN_INTERFACE_3 23
#define INWARD 31
#define OUTWARD 32
#define TARGET_PMT1 0
#define TARGET_PMT2 1
#define TARGET_PMT3 2
#define TARGET_WALL 3
#define TARGET_FLOOR 4
#define BOUNCED 10


namespace G4XTankSimulatorAG {

  class G4XTankFastCerenkov : public G4VContinuousProcess {

  public:

    G4XTankFastCerenkov(const G4String& processName = "Cerenkov");
    ~G4XTankFastCerenkov();     
    
    static void GetDataFromConstruction();
    G4bool IsApplicable(const G4ParticleDefinition &);
    void SetTrackSecondariesFirst(const G4bool);
    void SetMaxNumPhotonsPerStep(const G4int);
    G4PhysicsTable* GetPhysicsTable() const;
    void DumpPhysicsTable() const; 
    G4VParticleChange* AlongStepDoIt(const G4Track &, const G4Step &);
    G4double GetContinuousStepLimit(const G4Track &, G4double, G4double,
                                    G4double &);
   
  private:
    
    void BuildThePhysicsTable();
    G4double GetAverageNumberOfPhotons(const G4DynamicParticle *,
                                       const G4Material *,
                                       const G4MaterialPropertyVector *) const;
    void TrackPhotonInTank();
    void PropogateInTank(G4int);

    template<int pmtId>
    G4bool TransitionToDome(G4int);
    template<int pmtId>
    G4bool PropogateInDome(G4int &);
    template<int pmtId>
    G4bool TransitionToInterface(G4int &);
    template<int pmtId>
    G4bool PropogateInInterface(G4int &);

    G4double GetSphereIntersect(const G4ThreeVector &, const G4double) const;
    G4double GetEllipsoidIntersect(const G4ThreeVector &, 
                                   const G4double, const G4double) const;
    void DoRayleighScatter();
    G4bool ScatterOffRoof();
    G4bool ScatterOffWall();
    G4bool ScatterOffFloor();
    void BackScatter();
    void DiffuseScatterHorizontal(const G4ThreeVector &);
    void DiffuseScatterVertical(const G4ThreeVector &);
    void SpikeScatter(const G4ThreeVector &);
    void LobeScatterHorizontal(const G4ThreeVector &);
    void LobeScatterVertical(const G4ThreeVector &);
    
    G4bool fTrackSecondariesFirst;
    G4int fMaxPhotons;
    G4ThreeVector fPhotonMomentum, fPhotonPolarization, fPhotonPosition;
    G4double fPhotonTime, fSampledMomentum;

    G4double  fInvMFP, fRayleighFrac, fAbsorbFrac;

    static double fTankHeight, fTankHalfHeight, fTankThickness, fTankRadius, fTankRadiusSq;
    static double fDomeRadius, fInterfaceRadius, fPMTRadius, fHeightz;
    static double fDomeRadiusz, fInterfaceRadiusz, fPMTRadiusz;
    static double fDomeRadiusSq, fInterfaceRadiusSq, fPMTRadiusSq;
    static double fDomeRadiuszSq, fInterfaceRadiuszSq, fPMTRadiuszSq;
    
    static double  fSigmaAlpha;

    static utl::TabulatedFunction fQE;
  
    static utl::TabulatedFunction fWaterAbsLength;
    static utl::TabulatedFunction fInterfaceRIndex, fDomeRIndex;
  
    static utl::TabulatedFunction fInterfaceAbsLength, fDomeAbsLength;
    static utl::TabulatedFunction fLinerReflectivity, fLinerSpecularLobe;
    static utl::TabulatedFunction fLinerSpecularSpike, fLinerBackscatter;
      
    
    // PMT positions (the zeroth position is dummy)
    static G4ThreeVector fPMTPos[4];
    // The height of the roof, referred to each PMT
    static G4double fRoofPos[4];

    G4double fSampledRI;
    G4double fRayleighConst;

  protected:
    
    G4PhysicsTable* fThePhysicsTable;
    
  };

}

#endif