void CAngularSpectralProperties::calculateAngularProperties( std::shared_ptr<CSpectralSample> const & t_SpectralSample, MaterialType const t_Type) { assert(t_SpectralSample != nullptr); auto aMeasuredData = t_SpectralSample->getMeasuredData(); if(m_Angle != 0) { auto aSourceData = t_SpectralSample->getSourceData(); auto aWavelengths = aMeasuredData->getWavelengths(); auto aT = aMeasuredData->properties(SampleData::T); assert(aT->size() == aWavelengths.size()); auto aRf = aMeasuredData->properties(SampleData::Rf); assert(aRf->size() == aWavelengths.size()); auto aRb = aMeasuredData->properties(SampleData::Rb); assert(aRb->size() == aWavelengths.size()); auto lowLambda = 0.3; auto highLambda = 2.5; // TODO: Only one side is measured and it is considered that front properties are equal // to back properties auto aTSolNorm = t_SpectralSample->getProperty(lowLambda, highLambda, Property::T, Side::Front); for(size_t i = 0; i < aWavelengths.size(); ++i) { auto ww = aWavelengths[i] * 1e-6; auto T = (*aT)[i].value(); auto Rf = (*aRf)[i].value(); auto Rb = (*aRb)[i].value(); auto aSurfaceType = coatingType.at(t_Type); auto aFrontFactory = CAngularPropertiesFactory(T, Rf, m_Thickness, aTSolNorm); auto aBackFactory = CAngularPropertiesFactory(T, Rb, m_Thickness, aTSolNorm); auto aFrontProperties = aFrontFactory.getAngularProperties(aSurfaceType); auto aBackProperties = aBackFactory.getAngularProperties(aSurfaceType); auto Tangle = aFrontProperties->transmittance(m_Angle, ww); auto Rfangle = aFrontProperties->reflectance(m_Angle, ww); auto Rbangle = aBackProperties->reflectance(m_Angle, ww); m_AngularData->addRecord(ww * 1e6, Tangle, Rfangle, Rbangle); } } else { m_AngularData = aMeasuredData; } }
//#TPT-Directive ElementHeader Element_FILT static int interactWavelengths(Particle* cpart, int origWl) // Returns the wavelengths in a particle after FILT interacts with it (e.g. a photon) // cpart is the FILT particle, origWl the original wavelengths in the interacting particle int Element_FILT::interactWavelengths(Particle* cpart, int origWl) { const int mask = 0x3FFFFFFF; int filtWl = getWavelengths(cpart); switch (cpart->tmp) { case 0: return filtWl; //Assign Colour case 1: return origWl & filtWl; //Filter Colour case 2: return origWl | filtWl; //Add Colour case 3: return origWl & (~filtWl); //Subtract colour of filt from colour of photon case 4: { int shift = int((cpart->temp-273.0f)*0.025f); if (shift<=0) shift = 1; return (origWl << shift) & mask; // red shift } case 5: { int shift = int((cpart->temp-273.0f)*0.025f); if (shift<=0) shift = 1; return (origWl >> shift) & mask; // blue shift } case 6: return origWl; // No change case 7: return origWl ^ filtWl; // XOR colours case 8: return (~origWl) & mask; // Invert colours case 9: { int t1 = (origWl & 0x0000FF)+(rand()%5)-2; int t2 = ((origWl & 0x00FF00)>>8)+(rand()%5)-2; int t3 = ((origWl & 0xFF0000)>>16)+(rand()%5)-2; return (origWl & 0xFF000000) | (t3<<16) | (t2<<8) | t1; } default: return filtWl; } }