Trigger board name according to PRR
[u/mrichter/AliRoot.git] / Flugg / WrapReg.cxx
1
2 // Flugg tag 
3
4 ///////////////////////////////////////////////////////////////////
5 //
6 // WrapReg.hh - Sara Vanini 
7 //
8 // Wrapper for scoring hits: previous step end-point is taken from 
9 // history (and compared with fluka region index, flukaReg),
10 // then the wrapper returns all the information regarding the 
11 // volume tree, i.e. returns indMother[] array with all the 
12 // mother volumes index and repMother[] array with all the 
13 // mother volumes repetition number.   
14 //
15 // modified: 16/III/99
16 // modified: 14/IV/00 ptrLttc included 
17 // modified: 24.10.00: by I. Hrivnacova
18 //   functions declarations separated from implementation
19 //   (moved to Wrappers.hh);
20 // modified: 17/06/02 by I. Gonzalez. STL migration.
21 //
22 ///////////////////////////////////////////////////////////////////
23
24 #include "Wrappers.hh"
25 #include "FGeometryInit.hh"
26 #include "NavHistWithCount.hh"
27 #include "G4VPhysicalVolume.hh"
28 #include "G4ThreeVector.hh"
29 #include "G4PhysicalVolumeStore.hh"
30 #include "globals.hh"
31
32
33 void rgrpwr(const G4int& flukaReg, const G4int& ptrLttc, G4int& g4Reg,
34             G4int* indMother, G4int* repMother, G4int& depthFluka)
35 {
36   //flag
37 #ifdef G4GEOMETRY_DEBUG
38   G4cout << "============= RGRPWR ==============" << G4endl;    
39   G4cout << "ptrLttc=" << ptrLttc << G4endl;
40 #endif 
41   
42   //Geoinit, Navigator, VolStore pointers
43   static FGeometryInit * ptrGeoInit = FGeometryInit::GetInstance();
44   G4PhysicalVolumeStore * pVolStore = G4PhysicalVolumeStore::GetInstance();
45   
46   //get jrLtGeant array and flag
47   G4int * jrLtGeant = ptrGeoInit->GetJrLtGeantArray();
48   G4int LttcFlagGeant = ptrGeoInit->GetLttcFlagGeant();
49   
50   G4bool foundHistory = false;
51   G4int i = LttcFlagGeant;
52   
53   while(!foundHistory && i>=0) {
54     if(jrLtGeant[i]==ptrLttc) foundHistory = true;
55     i -= 1;
56   }
57   
58   if(!foundHistory) {
59     G4cout << "* ERROR! History not in jrLtGeant!" << G4endl;
60     //only in debugging version....
61     assert(foundHistory);
62   }
63   else {
64     //get history pointer from ptrLttc
65     NavHistWithCount* ptrNavHistCount = 
66       reinterpret_cast<NavHistWithCount*>(ptrLttc);
67     G4NavigationHistory* ptrNavHist = ptrNavHistCount->GetNavHistPtr();
68     
69     G4VPhysicalVolume* ptrVolHistory = 0;
70     G4int volHistIndex = 0;
71     G4int depth = ptrNavHist->GetDepth();
72     for(G4int h=0; h<=depth; h++) {
73       ptrVolHistory = ptrNavHist->GetVolume(h);
74       //
75       volHistIndex = ~0;
76       for (unsigned int i=0; i<pVolStore->size(); i++)
77         if ((*pVolStore)[i] == ptrVolHistory) 
78           volHistIndex = i; 
79       if (volHistIndex==(~0)) {
80         G4cerr << "FLUGG: Problem in routine WrapReg tryingto find volume after step" << G4endl;
81         exit(-999);
82       }
83       //volHistIndex = G4int(pVolStore->index(ptrVolHistory));
84       //
85       indMother[h] = volHistIndex+1;
86       if(ptrVolHistory->IsReplicated()) {
87         //true if volume is replica or parameterized,
88         //false for placements; repetition numbers 
89         //are set: 1,2,3,etc.; 0 for placed volumes.
90         repMother[h] = 1+G4int(ptrNavHist->GetReplicaNo(h));
91       }
92       else {
93         repMother[h] = 0;
94       }
95       
96 #ifdef  G4GEOMETRY_DEBUG
97       //  G4cout<<"Level="<<h<<"   : ";
98       //  G4cout<<"Region="<<indMother[h]<<"   (repetition="<<
99       //    repMother[h]<<")"<<G4endl;
100 #endif
101     }
102     
103     //compute new region index
104     G4int volIndex = ~0;
105     for (unsigned int i=0; i<pVolStore->size(); i++)
106       if ((*pVolStore)[i] == ptrVolHistory) 
107         volIndex = i;
108     //G4int volIndex=G4int(pVolStore->index(ptrVolHistory));
109     if (volIndex==(~0)) {
110       G4cerr << "FLUGG: Problem in routine WrapReg tryingto find volume after step" << G4endl;
111       exit(-999);
112     }
113     
114     g4Reg=volIndex+1;
115     
116     depthFluka = depth;
117   }
118   
119 }
120
121
122
123
124