54{
55 if (
m_verboseLevel>5) { G4cout <<
"######### Method TRTParametersForBarrelHits::DefineParameters" << G4endl; }
56
57
58 int numberOfStrawLayersA =
60 int numberOfStrawLayersB =
62 int numberOfStrawLayersC =
64
65
66 int* numberOfStrawsInLayersA = new int[numberOfStrawLayersA];
67 m_pParameters->GetIntegerArray(
"NumberOfStrawsInLayersA", numberOfStrawLayersA,
68 numberOfStrawsInLayersA);
69
70 int* integralDistributionOfStrawsA = new int[numberOfStrawLayersA];
71 integralDistributionOfStrawsA[0] = numberOfStrawsInLayersA[0];
73 for (i = 1;
i < numberOfStrawLayersA; ++
i)
74 integralDistributionOfStrawsA[i] = integralDistributionOfStrawsA[i - 1] +
75 numberOfStrawsInLayersA[i];
76
77 for (i = 0;
i < numberOfStrawLayersA; ++
i)
79 integralDistributionOfStrawsA[i];
80
81 delete [] numberOfStrawsInLayersA;
82
83
84 int* numberOfStrawsInLayersB = new int[numberOfStrawLayersB];
85 m_pParameters->GetIntegerArray(
"NumberOfStrawsInLayersB", numberOfStrawLayersB,
86 numberOfStrawsInLayersB);
87
88 int* integralDistributionOfStrawsB = new int[numberOfStrawLayersB];
89 integralDistributionOfStrawsB[0] = numberOfStrawsInLayersB[0];
90 for (i = 1;
i < numberOfStrawLayersB; ++
i)
91 integralDistributionOfStrawsB[i] = integralDistributionOfStrawsB[i - 1] +
92 numberOfStrawsInLayersB[i];
93
94 for (i = 0;
i < numberOfStrawLayersB; ++
i)
96 integralDistributionOfStrawsB[i];
97
98 delete [] numberOfStrawsInLayersB;
99
100
101 int* numberOfStrawsInLayersC = new int[numberOfStrawLayersC];
102 m_pParameters->GetIntegerArray(
"NumberOfStrawsInLayersC", numberOfStrawLayersC,
103 numberOfStrawsInLayersC);
104
105 int* integralDistributionOfStrawsC = new int[numberOfStrawLayersC];
106 integralDistributionOfStrawsC[0] = numberOfStrawsInLayersC[0];
107 for (i = 1;
i < numberOfStrawLayersC; ++
i)
108 integralDistributionOfStrawsC[i] = integralDistributionOfStrawsC[i - 1] +
109 numberOfStrawsInLayersC[i];
110
111 for (i = 0;
i < numberOfStrawLayersC; ++
i)
113 integralDistributionOfStrawsC[i];
114
115 delete [] numberOfStrawsInLayersC;
116
118 {
119
120 int layerID = 0;
122 {
123 if (i + 1 > integralDistributionOfStrawsA[layerID])
124 ++layerID;
126 }
127
128 layerID = 0;
130 {
131 if (i + 1 > integralDistributionOfStrawsB[layerID])
132 ++layerID;
134 }
135
136 layerID = 0;
138 {
139 if (i + 1 > integralDistributionOfStrawsC[layerID])
140 ++layerID;
142 }
143 }
144
145 delete [] integralDistributionOfStrawsA;
146 delete [] integralDistributionOfStrawsB;
147 delete [] integralDistributionOfStrawsC;
148
149
150 if (
m_verboseLevel>5) { G4cout <<
"######### Method TRTParametersForBarrelHits::DefineParameters done" << G4endl; }
151
152}