Geant4_10
G4GIDI_Misc.cc
Go to the documentation of this file.
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 # <<BEGIN-copyright>>
28 # Copyright (c) 2010, Lawrence Livermore National Security, LLC.
29 # Produced at the Lawrence Livermore National Laboratory
30 # Written by Bret R. Beck, beck6@llnl.gov.
31 # CODE-461393
32 # All rights reserved.
33 #
34 # This file is part of GIDI. For details, see nuclear.llnl.gov.
35 # Please also read the "Additional BSD Notice" at nuclear.llnl.gov.
36 #
37 # Redistribution and use in source and binary forms, with or without modification,
38 # are permitted provided that the following conditions are met:
39 #
40 # 1) Redistributions of source code must retain the above copyright notice,
41 # this list of conditions and the disclaimer below.
42 # 2) Redistributions in binary form must reproduce the above copyright notice,
43 # this list of conditions and the disclaimer (as noted below) in the
44 # documentation and/or other materials provided with the distribution.
45 # 3) Neither the name of the LLNS/LLNL nor the names of its contributors may be
46 # used to endorse or promote products derived from this software without
47 # specific prior written permission.
48 #
49 # THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND ANY
50 # EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES
51 # OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT
52 # SHALL LAWRENCE LIVERMORE NATIONAL SECURITY, LLC, THE U.S. DEPARTMENT OF ENERGY OR
53 # CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
54 # CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
55 # OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED
56 # AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
57 # (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE,
58 # EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
59 # <<END-copyright>>
60 */
61 #include <iostream>
62 #include <string>
63 #include <vector>
64 #include <xData.h>
65 #include <tpia_target.h>
66 #include <tpia_misc.h>
67 #include <string.h>
68 #include "G4GIDI_Misc.hh"
69 
70 using namespace std;
71 using namespace GIDI;
72 
73 /*
74 ***************************************************************
75 */
76 char *G4GIDI_Misc_Z_A_m_ToName( int iZ, int iA, int im ) {
77 
78  const char *Z = tpia_misc_ZToSymbol( iZ );
79  char S[128], mS[32], *name;
80 
81  if( Z == NULL ) return( NULL );
82  if( iA == 0 ) {
83  if( im != 0 ) return( NULL );
84  sprintf( S, "%s_natural", Z ); }
85  else {
86  sprintf( S, "%s_%d", Z, iA );
87  if( im != 0 ) {
88  sprintf( mS, "_m%d", im );
89  strcat( S, mS );
90  }
91  }
92  name = (char *) xData_malloc2( NULL, strlen( S ) + 1, 0, "name" );
93  if( name != NULL ) strcpy( name, S );
94  return( name );
95 }
96 /*
97 ***************************************************************
98 */
99 char *G4GIDI_Misc_channelCompound( char *particle1, char *particle2 ) {
100 
101  int Z1, A1, m1, Z2, A2, m2;
102 
103  if( tpia_miscNameToZAm( NULL, particle1, &Z1, &A1, &m1 ) ) return( NULL );
104  if( tpia_miscNameToZAm( NULL, particle2, &Z2, &A2, &m2 ) ) return( NULL );
105  if( A1 == 0 ) A2 = 0;
106  if( A2 == 0 ) A1 = 0;
107  return( G4GIDI_Misc_Z_A_m_ToName( Z1 + Z2, A1 + A2, 0 ) );
108 }
109 /*
110 ***************************************************************
111 */
112 int G4GIDI_Misc_channelProductsCompare( tpia_channel *channel, int nProducts, char **productNames ) {
113 
114  int i;
115  tpia_product *product;
116 
117  if( channel->decayChannel.numberOfProducts != nProducts ) return( 0 );
118  for( product = tpia_channel_getFirstProduct( channel ), i = 0; product != NULL; product = tpia_decayChannel_getNextProduct( product ), i++ ) {
119  if( strcmp( product->productID->name, productNames[i] ) ) return( 0 );
120  }
121  return( 1 );
122 }
123 /*
124 ***************************************************************
125 */
126 //int getNamesOfAvailableTargets_walker( tpia_mapEntry *entry, int level, void *userData ) {
128 
129  vector<string> *listOfTargets = (vector<string> *) userData;
130  vector<string>::iterator iter;
131 
132  if( entry->type != tpia_mapEntry_type_target ) return( 0 );
133  for( iter = listOfTargets->begin( ); iter != listOfTargets->end( ); iter++ ) {
134  if( entry->targetName == iter->c_str( ) ) return( 0 );
135  }
136  listOfTargets->push_back( entry->targetName );
137  return( 0 );
138 }
enum tpia_mapEntry_type type
Definition: tpia_map.h:64
Double_t Z2
Definition: plot.C:266
tpia_decayChannel decayChannel
Definition: tpia_target.h:264
const XML_Char * name
Definition: expat.h:151
#define xData_malloc2(smr, size, zero, forItem)
Definition: xData.h:313
#define userData
Definition: xmlparse.cc:555
tpia_product * tpia_channel_getFirstProduct(tpia_channel *channel)
int tpia_miscNameToZAm(statusMessageReporting *smr, const char *name, int *Z, int *A, int *m)
Definition: tpia_misc.cc:105
char * targetName
Definition: tpia_map.h:69
int getNamesOfAvailableTargets_walker(GIDI::tpia_mapEntry *entry, int level, void *userData)
char * G4GIDI_Misc_Z_A_m_ToName(int iZ, int iA, int im=0)
Definition: G4GIDI_Misc.cc:76
const char * tpia_misc_ZToSymbol(int iZ)
Definition: tpia_misc.cc:83
Float_t Z
Definition: plot.C:39
tpia_particle * productID
Definition: tpia_target.h:236
Definition: tpia_map.h:62
int G4GIDI_Misc_channelProductsCompare(GIDI::tpia_channel *channel, int nProducts, char **productNames)
Double_t Z1
Definition: plot.C:266
tpia_product * tpia_decayChannel_getNextProduct(tpia_product *product)
char * G4GIDI_Misc_channelCompound(char *particle1, char *particle2)
Definition: G4GIDI_Misc.cc:99