Rev 84 | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed
| Rev 84 | Rev 85 | ||
|---|---|---|---|
| Line 543... | Line 543... | ||
| 543 | 543 | ||
| 544 | // ---------------------------------------------------------------------------- |
544 | // ---------------------------------------------------------------------------- |
| 545 | // --------------- reflection at "reflection" probability --------------------- |
545 | // --------------- reflection at "reflection" probability --------------------- |
| 546 | // ---------------------------------------------------------------------------- |
546 | // ---------------------------------------------------------------------------- |
| 547 | case SURF_REFLE: |
547 | case SURF_REFLE: |
| 548 | p_ref |
548 | //p_ref = rand.Uniform(0.0, 1.0); |
| 549 | if(p_ref < reflection) { // se odbije |
549 | if(p_ref < reflection) { // se odbije |
| 550 | cosTi = in.GetK() * n; |
550 | cosTi = in.GetK() * n; |
| 551 | transmit = in.GetK() - 2*cosTi*n; |
551 | transmit = in.GetK() - 2*cosTi*n; |
| 552 | out.Set(intersect, transmit); |
552 | out.Set(intersect, transmit); |
| 553 | return REFLECTION; //sdhfvjhsdbfjhsdbcvjhsb |
553 | return REFLECTION; //sdhfvjhsdbfjhsdbcvjhsb |
| Line 558... | Line 558... | ||
| 558 | } |
558 | } |
| 559 | break; |
559 | break; |
| 560 | 560 | ||
| 561 | // total reflection from n1 to n2 with R probbability |
561 | // total reflection from n1 to n2 with R probbability |
| 562 | case SURF_IMPER: |
562 | case SURF_IMPER: |
| 563 | p_ref |
563 | //p_ref = rand.Uniform(0.0, 1.0); |
| 564 | if(p_ref < reflection) { // se odbije |
564 | if(p_ref < reflection) { // se odbije |
| 565 | cosTi = in.GetK() * n; |
565 | cosTi = in.GetK() * n; |
| 566 | if(TMath::Abs(cosTi) < cosTtotal) { // totalni odboj |
566 | if(TMath::Abs(cosTi) < cosTtotal) { // totalni odboj |
| 567 | transmit = in.GetK() - 2*cosTi*n; |
567 | transmit = in.GetK() - 2*cosTi*n; |
| 568 | out.Set(intersect, transmit); |
568 | out.Set(intersect, transmit); |
| Line 641... | Line 641... | ||
| 641 | noCoupling.FlipN(); |
641 | noCoupling.FlipN(); |
| 642 | // grease = specific pattern area of coupling |
642 | // grease = specific pattern area of coupling |
| 643 | TVector3 activePosition(center); |
643 | TVector3 activePosition(center); |
| 644 | activePosition += TVector3(_d, 0, 0); |
644 | activePosition += TVector3(_d, 0, 0); |
| 645 | TVector3 normal(1,0,0); |
645 | TVector3 normal(1,0,0); |
| 646 | grease.Set(activePosition, normal, 0. |
646 | grease.Set(activePosition, normal, 0.80*a/2.0); |
| 647 | 647 | ||
| 648 | if(fresnel) for(int i=0; i<6; i++) s_side[i].SetFresnel(1); |
648 | if(fresnel) for(int i=0; i<6; i++) s_side[i].SetFresnel(1); |
| 649 | 649 | ||
| 650 | // statistics histograms |
650 | // statistics histograms |
| 651 | //hfate = (TH1F)gROOT->FindObject("hfate"); //if(hfate) delete hfate; |
651 | //hfate = (TH1F)gROOT->FindObject("hfate"); //if(hfate) delete hfate; |