PhotonFullBlueTracker.cxx

Go to the documentation of this file.

#include "PhotonFullBlueTracker.h"
#include "PhotonTransportMaker.h"
#include "ScintPhoton.h"
#include "MessageService/MsgService.h"
#include "Plex/PlexHandle.h"

CVSID( " $Id: PhotonFullBlueTracker.cxx,v 1.10 2006/12/01 20:12:11 rhatcher Exp $" );

ClassImp(PhotonFullBlueTracker)

PhotonTransportMakerProxy<PhotonFullBlueTracker>
        gPhotonFullBlueTrackerProxy("photonFullBlueTracker"); 


// Constructor
PhotonFullBlueTracker::PhotonFullBlueTracker() :
  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);
  }
}

// Configure the model.
void 
PhotonFullBlueTracker::Configure( const Registry& r )
{
  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 )
{
  // 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 )
{
  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()); 
}

---