loon: PhotonFullFibreModel Class Reference

PhotonFullFibreModel::PhotonFullFibreModel ( )

Definition at line 17 of file PhotonFullFibreModel.cxx.

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

00017                                            :
00018   fFibreAbsorbTable("PHOTONFIBREABSORB"),
00019   fGreenSpectrumTable("PHOTONGREENSPECTRUM")
00020 {
00021   Configure(PhotonConfiguration());
00022     
00023   if (fFibreAbsorbTable.GetNumRows()<1) {
00024     MSG("Photon",Msg::kError) << "Can't find data for table PHOTONFIBREABSORB. Aborting!\n";
00025     assert(0);
00026   }
00027   if (fGreenSpectrumTable.GetNumRows()<1) {
00028     MSG("Photon",Msg::kError) << "Can't find data for table PHOTONFIBREABSORB. Aborting!\n";
00029     assert(0);
00030   }
00031 }

PhotonFullFibreModel::ClassDef	(	PhotonFullFibreModel	,
		1
	)			`[private]`

void PhotonFullFibreModel::Configure ( const Registry & r ) [virtual]

Reimplemented from PhotonTransportModule.

Definition at line 36 of file PhotonFullFibreModel.cxx.

References fFibreCoreRadius, fFibreFlurorDecayTime, fFibreIndex, fFibreRadius, fInnerCladdingIndex, fInnerCladdingRadius, fOuterCladdingIndex, fScintIndex, Registry::Get(), Msg::kDebug, MSG, and PhotonTransportModule::SetWlsPrescaleFactor().

00037 {
00038   MSG("Photon",Msg::kDebug) << "Configuring FullFibreModel." << std::endl;
00039   double tmpd;
00040   if (r.Get("FibreRadius",tmpd))          fFibreRadius  = tmpd;
00041   if (r.Get("FibreCoreRadius",tmpd))      fFibreCoreRadius  = tmpd;
00042   if (r.Get("InnerCladdingRadius",tmpd))  fInnerCladdingRadius  = tmpd;
00043 
00044   if (r.Get("FibreFlurorDecayTime",tmpd)) fFibreFlurorDecayTime  = tmpd;
00045 
00046   if (r.Get("FibreIndex",tmpd))           fFibreIndex  = tmpd;
00047   if (r.Get("ScintIndex",tmpd))           fScintIndex  = tmpd;
00048   if (r.Get("OuterCladdingIndex",tmpd))   fOuterCladdingIndex  = tmpd;
00049   if (r.Get("InnerCladdingIndex",tmpd))   fInnerCladdingIndex  = tmpd;
00050 
00051   SetWlsPrescaleFactor(1.0); // No prescaling.
00052 }

Bool_t PhotonFullFibreModel::FibreHitToGreenPhoton	(	const UgliStripHandle &	inStrip,
		StripEnd::StripEnd_t	inDirection,
		DigiPhoton *	inBluePhoton,
		DigiPhoton *&	outGreenPhoton
	)			`[virtual]`

Reimplemented from PhotonTransportModule.

Definition at line 229 of file PhotonFullFibreModel.cxx.

References FindBlueAbsorptionPosition(), fLastResult, IsGreenTrapped(), and MakeGreenPhoton().

00234 {
00235   // Return a contained green photon.
00236   // This contains the full photon simulation.
00237   // Blue photons are tracked through the cladding and into the core
00238   // where they are absorbed by the correct attenuation.
00239   fLastResult = -1;
00240 
00241   outGreenPhoton = NULL;
00242   if(!inBluePhoton) { return false; };
00243 
00244   TVector3 flashPos;
00245   fLastResult = FindBlueAbsorptionPosition( inStrip,
00246                                             inBluePhoton,
00247                                             flashPos );
00248 
00249   if(fLastResult !=0 ) return false;
00250 
00251   // Loop until we have a valid solution.
00252 
00253   // Find the new direction of the green photon.
00254   // Go in the forward or backward direction, depending on the 
00255   // given parameter.
00256   
00257   if(outGreenPhoton) delete outGreenPhoton;
00258   outGreenPhoton = MakeGreenPhoton( inBluePhoton, 
00259                                     flashPos,
00260                                     inDirection);
00261 
00262   if(!outGreenPhoton) return false; // Huh?
00263 
00264   if(IsGreenTrapped(outGreenPhoton)) return true;
00265   
00266   // else:
00267   return false; // No trapped photon!
00268 }

Int_t PhotonFullFibreModel::FindBlueAbsorptionPosition	(	const UgliStripHandle &	inStrip,
		DigiPhoton *	inBluePhoton,
		TVector3 &	outPos
	)

Definition at line 64 of file PhotonFullFibreModel.cxx.

References fBlue_InnerHit, fBlue_OuterHit, fBlue_Start, fBlue_Stop, fFibreAbsorbTable, fFibreCoreRadius, fFibreIndex, fFibreRadius, fInnerCladdingIndex, fInnerCladdingRadius, fOuterCladdingIndex, PhotonTransportModule::fRandom, fScintIndex, DigiPhoton::GetDir(), UgliStripHandle::GetHalfThickness(), DigiPhoton::GetPos(), PhotonUtil::InteractWithCylinder(), PhotonLookupTable::Interpolate(), PhotonUtil::IntersectCylinder(), Msg::kWarning, MSG, Munits::nm, and DigiPhoton::WaveLength().

Referenced by FibreHitToGreenPhoton().

00067 {
00068   // Track the blue photon until it gets absorbed.
00069   // Return -ve number for a condition that shouldn't happen 
00070   // Return 0 if photon stopped
00071   // Return +ve number if photon escaped
00072 
00073   // Reset the blue photon coordinates so they're centered on our fibre.
00074   TVector3 x = inBluePhoton->GetPos();
00075   
00076   x -= TVector3(0,0,
00077                 inStrip.GetHalfThickness() - fFibreRadius);
00078   fBlue_Start = fBlue_Stop = fBlue_OuterHit = fBlue_InnerHit = x;
00079   
00080   TVector3 v = inBluePhoton->GetDir();
00081   
00082   // Interact with the outer cladding.
00083   int interact  = PhotonUtil::InteractWithCylinder( x,v,
00084                                                     fFibreRadius,
00085                                                     fScintIndex,
00086                                                     fOuterCladdingIndex);
00087   if(interact!=0) { 
00088     // It reflected, not refracted. How did this happen?
00089     fBlue_OuterHit = x + v*0.0001;
00090     MSG("Photon",Msg::kWarning) << "Didn't refract into fibre. Huh? x=(" 
00091                                 << x.x() << " " << x.y() << " " << x.z() << ") "
00092                                 << "v=(" << v.x() << " " << v.y() << " "  << v.z() <<")\n";
00093                                 
00094     return -2;
00095   } 
00096   // Interact with the inner cladding
00097   int inclad1 = PhotonUtil::InteractWithCylinder( x,v,
00098                                                    fInnerCladdingRadius,
00099                                                    fOuterCladdingIndex,                 
00100                                                    fInnerCladdingIndex );
00101 
00102   fBlue_OuterHit = x;
00103   if(inclad1!=0) { return 1; } // It got out.
00104 
00105   // Intersect with the core:
00106   int incore = PhotonUtil::InteractWithCylinder( x, v,
00107                                                  fFibreCoreRadius,
00108                                                  fInnerCladdingIndex,
00109                                                  fFibreIndex );
00110   fBlue_InnerHit = x;
00111   // Missed the core or reflected off it. Shucks.
00112   if(incore!=0) { return 2; };
00113   
00114   // Photon hit the core. See how far it could go.
00115   double d1, d2;
00116   incore = PhotonUtil::IntersectCylinder( x, v, fFibreCoreRadius, d1, d2 );
00117   
00118   // This shouldn't happen.. we've already said it hit the core.
00119   if(!incore) { return -3;}
00120 
00121   // Find wavelength of blue photon in nm.
00122   double wavelength = inBluePhoton->WaveLength()/Munits::nm;
00123   if(wavelength<=0) wavelength = 420; // Something reasonable.
00124   
00125   // Find absorption length of this photon in nm.
00126   double atten = fFibreAbsorbTable.Interpolate(wavelength);
00127   
00128   // find where the photon got absorbed.
00129   double d = d1 + fRandom->Exp(atten);
00130 
00131   // Compute the final postion.
00132   outPos = x + v*d; 
00133   fBlue_Stop = outPos;;
00134   if(d>d2) { return 3; } // the blue photon managed to escape through the other side.
00135   
00136   return 0;
00137 }

Int_t PhotonFullFibreModel::GetLastResult ( ) [inline]

Definition at line 43 of file PhotonFullFibreModel.h.

References fLastResult.

00043 { return fLastResult; }; // For debugging, plot making.

Bool_t PhotonFullFibreModel::IsGreenTrapped ( DigiPhoton * inGreen )

Definition at line 185 of file PhotonFullFibreModel.cxx.

References fFibreCoreRadius, fFibreIndex, fGreen_InnerHit, fGreen_OuterHit, fGreen_Start, fInnerCladdingIndex, fInnerCladdingRadius, fOuterCladdingIndex, DigiPhoton::GetDir(), DigiPhoton::GetPos(), PhotonUtil::InteractWithCylinder(), Msg::kError, and MSG.

Referenced by FibreHitToGreenPhoton().

00186 {
00187   TVector3 x   = inGreen->GetPos();
00188   TVector3 dir = inGreen->GetDir();
00189 
00190   // Special cases:
00191   if(dir.x()==0) return false; // Going sideways.
00192   if(fabs(dir.x())>=1.) return true; // Going straight down the pipe.
00193 
00194   fGreen_Start = x;
00195 
00196   int interact;
00197   // Interact with the inner/outer cladding interface.
00198   interact = PhotonUtil::InteractWithCylinder( x,dir,
00199                                                fFibreCoreRadius,
00200                                                fFibreIndex,
00201                                                fInnerCladdingIndex);
00202   if(interact==-1) {
00203     MSG("Photon",Msg::kError) << "Didn't hit inner wall IsGreenTrapped(). Huh?!?" << endl;
00204   };
00205   fGreen_InnerHit = x;
00206   if(interact==1) {
00207     // It reflected off the inner surface. Bonus.
00208     fGreen_OuterHit = x;
00209     return true;
00210   }
00211 
00212 
00213   // Interact with the outer cladding/scint interface. 
00214   interact = PhotonUtil::InteractWithCylinder( x,dir,
00215                                                fInnerCladdingRadius,
00216                                                fInnerCladdingIndex,
00217                                                fOuterCladdingIndex);  
00218   if(interact==-1) {
00219     MSG("Photon",Msg::kError) << "Didn't hit outer wall IsGreenTrapped(). Huh?!?" << endl;
00220   };
00221   fGreen_OuterHit = x;
00222 
00223   if(interact==0) return false; // Refracted out.
00224   return true; // Reflected in.
00225 }

DigiPhoton * PhotonFullFibreModel::MakeGreenPhoton	(	DigiPhoton *	inBluePhoton,
		const TVector3 &	inPos,
		StripEnd::StripEnd_t	inDirection
	)

Definition at line 141 of file PhotonFullFibreModel.cxx.

References fFibreFlurorDecayTime, fGreenSpectrumTable, PhotonTransportModule::fRandom, StripEnd::kNegative, StripEnd::kPositive, Munits::nm, DigiPhoton::ParentHit(), PhotonLookupTable::Pick(), and DigiPhoton::T().

Referenced by FibreHitToGreenPhoton().

00145 {
00146   // Choose a direction.
00147   double cosTheta  = 2.*fRandom->Rndm()-1.0;
00148   double sinTheta = TMath::Sin(TMath::ACos(cosTheta));
00149   double phi = 2.*TMath::Pi()*fRandom->Rndm();
00150   TVector3 dir( cosTheta, 
00151                 sinTheta*TMath::Cos(phi),
00152                 sinTheta*TMath::Sin(phi) );
00153   
00154   // Let's ensure it's going in the correct direction:
00155   switch(inDirection) {
00156   case StripEnd::kPositive:
00157     dir.SetX(fabs(dir.X()));
00158     break;
00159     
00160   case StripEnd::kNegative:
00161     dir.SetX(-fabs(dir.X()));
00162     break;
00163     
00164   default:
00165     // Do nothing.
00166     break;    
00167   }
00168 
00169   // Choose wavelength.
00170   double wavelength = fGreenSpectrumTable.Pick(fRandom->Rndm())*Munits::nm;
00171   
00172   double time = inBluePhoton->T() + fRandom->Exp(fFibreFlurorDecayTime);
00173 
00174   DigiPhoton* green =  new DigiPhoton( inBluePhoton->ParentHit(),
00175                                        wavelength,
00176                                        time,
00177                                        inPos,
00178                                        dir
00179                                        );
00180   return green;
00181 }

void PhotonFullFibreModel::Reset ( const VldContext & cx ) [protected, virtual]

Reimplemented from PhotonTransportModule.

Definition at line 55 of file PhotonFullFibreModel.cxx.

References fFibreAbsorbTable, fGreenSpectrumTable, and PhotonLookupTable::Reset().

00056 {
00057   // Ensure tables are up-to-date.
00058   fFibreAbsorbTable.Reset(cx);
00059   fGreenSpectrumTable.Reset(cx);
00060 }

TVector3 PhotonFullFibreModel::fBlue_InnerHit

Definition at line 47 of file PhotonFullFibreModel.h.

Referenced by FindBlueAbsorptionPosition().

TVector3 PhotonFullFibreModel::fBlue_OuterHit

Definition at line 46 of file PhotonFullFibreModel.h.

Referenced by FindBlueAbsorptionPosition().

TVector3 PhotonFullFibreModel::fBlue_Start

Definition at line 43 of file PhotonFullFibreModel.h.

Referenced by FindBlueAbsorptionPosition().

TVector3 PhotonFullFibreModel::fBlue_Stop

Definition at line 48 of file PhotonFullFibreModel.h.

Referenced by FindBlueAbsorptionPosition().

PhotonLookupTable PhotonFullFibreModel::fFibreAbsorbTable [private]

Definition at line 59 of file PhotonFullFibreModel.h.

Referenced by FindBlueAbsorptionPosition(), and Reset().

Double_t PhotonFullFibreModel::fFibreCoreRadius [private]

Definition at line 67 of file PhotonFullFibreModel.h.

Referenced by Configure(), FindBlueAbsorptionPosition(), and IsGreenTrapped().

Double_t PhotonFullFibreModel::fFibreFlurorDecayTime [private]

Definition at line 63 of file PhotonFullFibreModel.h.

Referenced by Configure(), and MakeGreenPhoton().

Double_t PhotonFullFibreModel::fFibreIndex [private]

Definition at line 66 of file PhotonFullFibreModel.h.

Referenced by Configure(), FindBlueAbsorptionPosition(), and IsGreenTrapped().

Double_t PhotonFullFibreModel::fFibreRadius [private]

Definition at line 62 of file PhotonFullFibreModel.h.

Referenced by Configure(), and FindBlueAbsorptionPosition().

TVector3 PhotonFullFibreModel::fGreen_InnerHit

Definition at line 50 of file PhotonFullFibreModel.h.

Referenced by IsGreenTrapped().

TVector3 PhotonFullFibreModel::fGreen_OuterHit

Definition at line 51 of file PhotonFullFibreModel.h.

Referenced by IsGreenTrapped().

TVector3 PhotonFullFibreModel::fGreen_Start

Definition at line 49 of file PhotonFullFibreModel.h.

Referenced by IsGreenTrapped().

PhotonLookupTable PhotonFullFibreModel::fGreenSpectrumTable [private]

Definition at line 60 of file PhotonFullFibreModel.h.

Referenced by MakeGreenPhoton(), and Reset().

Double_t PhotonFullFibreModel::fInnerCladdingIndex [private]

Definition at line 65 of file PhotonFullFibreModel.h.

Referenced by Configure(), FindBlueAbsorptionPosition(), and IsGreenTrapped().

Double_t PhotonFullFibreModel::fInnerCladdingRadius [private]

Definition at line 68 of file PhotonFullFibreModel.h.

Referenced by Configure(), FindBlueAbsorptionPosition(), and IsGreenTrapped().

Int_t PhotonFullFibreModel::fLastResult [private]

Definition at line 57 of file PhotonFullFibreModel.h.

Referenced by FibreHitToGreenPhoton(), and GetLastResult().

Double_t PhotonFullFibreModel::fOuterCladdingIndex [private]

Definition at line 64 of file PhotonFullFibreModel.h.

Referenced by Configure(), FindBlueAbsorptionPosition(), and IsGreenTrapped().

Double_t PhotonFullFibreModel::fScintIndex [private]

Definition at line 69 of file PhotonFullFibreModel.h.

Referenced by Configure(), and FindBlueAbsorptionPosition().

PhotonFullFibreModel Class Reference

Public Member Functions

Public Attributes

Protected Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation


Public Member Functions
	PhotonFullFibreModel ()
virtual void	Configure (const Registry &r)
Int_t	FindBlueAbsorptionPosition (const UgliStripHandle &inStrip, DigiPhoton *inBluePhoton, TVector3 &outPos)
DigiPhoton *	MakeGreenPhoton (DigiPhoton *inBluePhoton, const TVector3 &inPos, StripEnd::StripEnd_t inDirection)
Bool_t	IsGreenTrapped (DigiPhoton *inGreen)
virtual Bool_t	FibreHitToGreenPhoton (const UgliStripHandle &inStrip, StripEnd::StripEnd_t inDirection, DigiPhoton inBluePhoton, DigiPhoton &outGreenPhoton)
Int_t	GetLastResult ()
Public Attributes
TVector3	fBlue_Start
TVector3	fBlue_OuterHit
TVector3	fBlue_InnerHit
TVector3	fBlue_Stop
TVector3	fGreen_Start
TVector3	fGreen_InnerHit
TVector3	fGreen_OuterHit
Protected Member Functions
virtual void	Reset (const VldContext &cx)
Private Member Functions
	ClassDef (PhotonFullFibreModel, 1)
Private Attributes
Int_t	fLastResult
PhotonLookupTable	fFibreAbsorbTable
PhotonLookupTable	fGreenSpectrumTable
Double_t	fFibreRadius
Double_t	fFibreFlurorDecayTime
Double_t	fOuterCladdingIndex
Double_t	fInnerCladdingIndex
Double_t	fFibreIndex
Double_t	fFibreCoreRadius
Double_t	fInnerCladdingRadius
Double_t	fScintIndex