PhotonFullBlueTracker Class Reference

#include <PhotonFullBlueTracker.h>

Inheritance diagram for PhotonFullBlueTracker:

List of all members.

Public Member Functions
	PhotonFullBlueTracker ()
virtual void	Configure (const Registry &r)
virtual Bool_t	ScintPhotonToFibreHit (const DigiScintHit *hit, const UgliStripHandle &inStrip, DigiPhoton *&outPhoton)
Protected Member Functions
virtual void	Reset (const VldContext &cx)
Private Member Functions
	ClassDef (PhotonFullBlueTracker, 1)
Private Attributes
PhotonLookupTable	fBlueSpectrumTable
PhotonLookupTable	fScintAttenTable
PhotonLookupTable	fTiO2ReflecTable
PhotonLookupTable	fFibreAbsorbTable
Double_t	fScintDecayTime
Double_t	fScintAttenTuneFactor
Double_t	fTiO2ReflecTuneFactor

---

Detailed Description

Definition at line 20 of file PhotonFullBlueTracker.h.

---

Constructor & Destructor Documentation

PhotonFullBlueTracker::PhotonFullBlueTracker

(

)

Definition at line 17 of file PhotonFullBlueTracker.cxx.

References Msg::kError, MSG, and PhotonConfiguration().

                                             :
  fBlueSpectrumTable("PHOTONBLUESPECTRUM"),
  fScintAttenTable("PHOTONSCINTATTEN"),
  fTiO2ReflecTable("PHOTONTIO2PAINTREFLEC"),
  fFibreAbsorbTable("PHOTONFIBREABSORB"),
  fScintAttenTuneFactor(1.0),
  fTiO2ReflecTuneFactor(1.0) 
{
  Configure(PhotonConfiguration());

  if ( (fBlueSpectrumTable.GetNumRows()<1) ||
       (fScintAttenTable.GetNumRows()<1) ||
       (fTiO2ReflecTable.GetNumRows()<1) ||
       (fFibreAbsorbTable.GetNumRows()<1) ) {
    MSG("Photon",Msg::kError) << "Can't find data for one of tables! Aborting!\n";
    assert(0);
  }
}

---

Member Function Documentation

PhotonFullBlueTracker::ClassDef	(	PhotonFullBlueTracker	,
		1
	)			[private]

void PhotonFullBlueTracker::Configure

(

const Registry &

r

)

[virtual]

Reimplemented from PhotonTransportModule.

Definition at line 39 of file PhotonFullBlueTracker.cxx.

References ScintPhoton::Configure(), fScintAttenTuneFactor, fScintDecayTime, fTiO2ReflecTuneFactor, Registry::Get(), Msg::kDebug, MSG, and PhotonTransportModule::SetBluePrescaleFactor().

{
  MSG("Photon",Msg::kDebug) << "Configuring FullBlueTrackerModel." << std::endl;

  double tmpd;
  if(r.Get("ScintDecayTime",tmpd)) fScintDecayTime = tmpd;
  if(r.Get("ScintAttenTuneFactor",tmpd)) {
    fScintAttenTuneFactor = tmpd;
    MSG("Photon",Msg::kDebug) << "Setting fScintAttenTuneFactor to " << fScintAttenTuneFactor << endl;
  }
  if(r.Get("TiO2ReflecTuneFactor",tmpd))
    fTiO2ReflecTuneFactor = tmpd;

  // The other classses in the model need configuring:
  ScintPhoton::Configure(r);

  SetBluePrescaleFactor(1.0); // No prescaling.
}

void PhotonFullBlueTracker::Reset

(

const VldContext &

cx

)

[protected, virtual]

Reimplemented from PhotonTransportModule.

Definition at line 59 of file PhotonFullBlueTracker.cxx.

References fBlueSpectrumTable, fFibreAbsorbTable, fScintAttenTable, fTiO2ReflecTable, and PhotonLookupTable::Reset().

{
  // Ensure tables are up-to-date.
  fBlueSpectrumTable.Reset(cx);
  fScintAttenTable.Reset(cx);
  fTiO2ReflecTable.Reset(cx);
  fFibreAbsorbTable.Reset(cx);
}

Bool_t PhotonFullBlueTracker::ScintPhotonToFibreHit	(	const DigiScintHit *	hit,
		const UgliStripHandle &	inStrip,
		DigiPhoton *&	outPhoton
	)			[virtual]

Reimplemented from PhotonTransportModule.

Definition at line 72 of file PhotonFullBlueTracker.cxx.

References fBlueSpectrumTable, PhotonTransportModule::fRandom, fScintAttenTable, fScintAttenTuneFactor, fScintDecayTime, fTiO2ReflecTable, fTiO2ReflecTuneFactor, ScintPhoton::GeomError(), UgliStripHandle::GetHalfLength(), ScintPhoton::InFibre(), PhotonLookupTable::Interpolate(), StripEnd::kNegative, StripEnd::kPositive, Munits::nm, UgliStripHandle::PartialLength(), PhotonLookupTable::Pick(), ScintPhoton::Propagate(), ScintPhoton::Reset(), ScintPhoton::SetScintAtten(), ScintPhoton::SetTiO2Reflec(), DigiPhoton::SetWavelength(), DigiScintHit::T1(), DigiScintHit::T2(), UgliStripHandle::WlsBypass(), DigiPhoton::X(), DigiScintHit::X1(), DigiScintHit::X2(), DigiScintHit::Y1(), DigiScintHit::Y2(), DigiScintHit::Z1(), and DigiScintHit::Z2().

{
  outPhoton = NULL;
  // This default method uses the ScintPhoton class to do the dirty work.

  
  if(!hit) return false;

  // Make the photon.
  ScintPhoton*  photon = new ScintPhoton();  

  // Find the location/time for the blue photon.

  Double_t delta = fRandom->Rndm();
  TVector3 x1(hit->X1(),hit->Y1(),hit->Z1());
  TVector3 x2(hit->X2(),hit->Y2(),hit->Z2());

  TVector3 p = x1 + delta*(x2-x1);
  Double_t t = hit->T1() + delta * (hit->T2() - hit->T1()) +
    fRandom->Exp(fScintDecayTime);



  // The point of this loop is to ensure that we actually manage to 
  // simulate a photon successfully.. it doesn't just disappear, but
  // is successfully tracked to SOMEWHERE.
  bool ok;
  do {
    TVector3 dir(0,0,0); // let the scint photon do it.
    photon->Reset(fRandom,
                  hit,
                  inStrip,
                  p, t,
                  dir);

      // Find a wavelength for the photon.
    double wavelength_nm = fBlueSpectrumTable.Pick(fRandom->Rndm());
    photon->SetWavelength( wavelength_nm * Munits::nm );
    
    // Find the relevant constants.
    // Must set these AFTER the Reset() command, which re-initialized the
    // defaults.
    photon->SetScintAtten( fScintAttenTable.Interpolate(wavelength_nm) *
                           fScintAttenTuneFactor );
    photon->SetTiO2Reflec( fTiO2ReflecTable.Interpolate(wavelength_nm) *
                           fTiO2ReflecTuneFactor );
    
    photon->Propagate();
    ok = ! photon->GeomError(); // Geometry errors are bad.
  } while(!ok);

  // Successful photon.. may or may not have hit green.
  outPhoton = (DigiPhoton*) photon;


  // Simple check for intra-strip light leakage. Not completely correct; this SHOULD
  // be part of the photon tracking code.. but simple and straightforward.
  if(inStrip.WlsBypass() > 0.) {
    float x = outPhoton->X();
    // Deal with the case where the hit was near the edge of a two-part strip.
    // See if the hit is contained or not.
    if( ( (x + inStrip.GetHalfLength()) > inStrip.PartialLength(StripEnd::kNegative) ) 
        && ( (inStrip.GetHalfLength() - x) > inStrip.PartialLength(StripEnd::kPositive) )
        ) {
      
      // The photon got leaked out into the coil hole bypass
      return false;
    }    
  }

  return (photon->InFibre()); 
}

---

Member Data Documentation

PhotonLookupTable PhotonFullBlueTracker::fBlueSpectrumTable [private]

Definition at line 34 of file PhotonFullBlueTracker.h.

Referenced by Reset(), and ScintPhotonToFibreHit().

PhotonLookupTable PhotonFullBlueTracker::fFibreAbsorbTable [private]

Definition at line 37 of file PhotonFullBlueTracker.h.

Referenced by Reset().

PhotonLookupTable PhotonFullBlueTracker::fScintAttenTable [private]

Definition at line 35 of file PhotonFullBlueTracker.h.

Referenced by Reset(), and ScintPhotonToFibreHit().

Double_t PhotonFullBlueTracker::fScintAttenTuneFactor [private]

Definition at line 40 of file PhotonFullBlueTracker.h.

Referenced by Configure(), and ScintPhotonToFibreHit().

Double_t PhotonFullBlueTracker::fScintDecayTime [private]

Definition at line 39 of file PhotonFullBlueTracker.h.

Referenced by Configure(), and ScintPhotonToFibreHit().

PhotonLookupTable PhotonFullBlueTracker::fTiO2ReflecTable [private]

Definition at line 36 of file PhotonFullBlueTracker.h.

Referenced by Reset(), and ScintPhotonToFibreHit().

Double_t PhotonFullBlueTracker::fTiO2ReflecTuneFactor [private]

Definition at line 41 of file PhotonFullBlueTracker.h.

Referenced by Configure(), and ScintPhotonToFibreHit().

---

The documentation for this class was generated from the following files:

• PhotonFullBlueTracker.h
• PhotonFullBlueTracker.cxx

---