Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/advanced/ChargeExchangeMC/src/CexmcAngularRange.cc

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

  1 //
  2 // ********************************************************************
  3 // * License and Disclaimer                                           *
  4 // *                                                                  *
  5 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
  6 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
  7 // * conditions of the Geant4 Software License,  included in the file *
  8 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
  9 // * include a list of copyright holders.                             *
 10 // *                                                                  *
 11 // * Neither the authors of this software system, nor their employing *
 12 // * institutes,nor the agencies providing financial support for this *
 13 // * work  make  any representation or  warranty, express or implied, *
 14 // * regarding  this  software system or assume any liability for its *
 15 // * use.  Please see the license in the file  LICENSE  and URL above *
 16 // * for the full disclaimer and the limitation of liability.         *
 17 // *                                                                  *
 18 // * This  code  implementation is the result of  the  scientific and *
 19 // * technical work of the GEANT4 collaboration.                      *
 20 // * By using,  copying,  modifying or  distributing the software (or *
 21 // * any work based  on the software)  you  agree  to acknowledge its *
 22 // * use  in  resulting  scientific  publications,  and indicate your *
 23 // * acceptance of all terms of the Geant4 Software license.          *
 24 // ********************************************************************
 25 //
 26 /*
 27  * ============================================================================
 28  *
 29  *       Filename:  CexmcAngularRange.cc
 30  *
 31  *    Description:  angular range object
 32  *
 33  *        Version:  1.0
 34  *        Created:  28.12.2009 22:35:07
 35  *       Revision:  none
 36  *       Compiler:  gcc
 37  *
 38  *         Author:  Alexey Radkov (), 
 39  *        Company:  PNPI
 40  *
 41  * ============================================================================
 42  */
 43 
 44 #include <algorithm>
 45 #include <iostream>
 46 #include <iomanip>
 47 #include <cmath>
 48 #include "CexmcAngularRange.hh"
 49 
 50 
 51 void  GetNormalizedAngularRange( const CexmcAngularRangeList &  src,
 52                                  CexmcAngularRangeList &  dst )
 53 {
 54     dst = src;
 55     if ( dst.size() < 2 )
 56         return;
 57 
 58     std::sort( dst.begin(), dst.end() );
 59 
 60     const G4double  epsilon( 1E-7 );
 61 
 62     for ( CexmcAngularRangeList::iterator  k( dst.begin() + 1 );
 63                                                             k != dst.end(); )
 64     {
 65         if ( std::fabs( k->top - ( k - 1 )->top ) < epsilon ||
 66              k->bottom + epsilon >= ( k - 1 )->bottom )
 67         {
 68             k = dst.erase( k );
 69             continue;
 70         }
 71         if ( k->top + epsilon >= ( k - 1 )->bottom )
 72         {
 73             ( k - 1 )->bottom = k->bottom;
 74             k = dst.erase( k );
 75             continue;
 76         }
 77         ++k;
 78     }
 79 }
 80 
 81 
 82 void  GetAngularGaps( const CexmcAngularRangeList &  src,
 83                       CexmcAngularRangeList &  dst )
 84 {
 85     if ( src.empty() )
 86     {
 87         dst.push_back( CexmcAngularRange( 1.0, -1.0 , 0 ) );
 88         return;
 89     }
 90 
 91     CexmcAngularRangeList  normalizedAngularRanges;
 92     GetNormalizedAngularRange( src, normalizedAngularRanges );
 93 
 94     G4int  index( 0 );
 95     if ( normalizedAngularRanges[ 0 ].top < 1.0 )
 96         dst.push_back( CexmcAngularRange(
 97                         1.0, normalizedAngularRanges[ 0 ].top, index++ ) );
 98 
 99     for ( CexmcAngularRangeList::iterator
100             k( normalizedAngularRanges.begin() );
101             k != normalizedAngularRanges.end(); ++k )
102     {
103         if ( k + 1 == normalizedAngularRanges.end() )
104             break;
105         dst.push_back( CexmcAngularRange(
106                         k->bottom, ( k + 1 )->top, index++ ) );
107     }
108 
109     if ( normalizedAngularRanges.back().bottom > -1.0 )
110         dst.push_back( CexmcAngularRange(
111                         normalizedAngularRanges.back().bottom, -1.0, index ) );
112 }
113 
114 
115 std::ostream &  operator<<( std::ostream &  out,
116                             const CexmcAngularRange &  angularRange )
117 {
118     std::ostream::fmtflags  savedFlags( out.flags() );
119     std::streamsize         prec( out.precision() );
120 
121     out.precision( 4 );
122     out.flags( std::ios::fixed );
123 
124     out << std::setw( 2 ) << angularRange.index + 1 << " [" << std::setw( 7 ) <<
125            angularRange.top << ", " << std::setw( 7 ) << angularRange.bottom <<
126            ")";
127 
128     out.precision( prec );
129     out.flags( savedFlags );
130 
131     return out;
132 }
133 
134 
135 std::ostream &  operator<<( std::ostream &  out,
136                             const CexmcAngularRangeList &  angularRanges )
137 {
138     out << std::endl;
139     for ( CexmcAngularRangeList::const_iterator  k( angularRanges.begin() );
140                                                 k != angularRanges.end(); ++k )
141     {
142         out << "                 " << *k << std::endl;
143     }
144 
145     return out;
146 }
147 
148