@@ -65,19 +65,27 @@ void Geometry::setCsideModules()
6565 for (Int_t i = 0 ; i < 6 ; i++) {
6666 gridpoints[i] = crad * TMath::Sin ((1 - 1 / (2 * TMath::Abs (grdin[i]))) * grdin[i] * btta);
6767 }
68+ Double_t xi[NCellsC] = {-15.038271418735729 , 15.038271418735729 ,
69+ -15.003757581112167 , 15.003757581112167 , -9.02690018974363 ,
70+ 9.02690018974363 , -9.026897413747076 , 9.026897413747076 ,
71+ -9.026896531935773 , 9.026896531935773 , -3.0004568618531313 ,
72+ 3.0004568618531313 , -3.0270795197907225 , 3.0270795197907225 ,
73+ 3.0003978432927543 , -3.0003978432927543 , 3.0270569670429572 ,
74+ -3.0270569670429572 , 9.026750365564254 , -9.026750365564254 ,
75+ 9.026837450695885 , -9.026837450695885 , 9.026849243816981 ,
76+ -9.026849243816981 , 15.038129472387304 , -15.038129472387304 ,
77+ 15.003621961057961 , -15.003621961057961 };
78+ Double_t yi[NCellsC] = {3.1599494336464455 , -3.1599494336464455 ,
79+ 9.165191680982874 , -9.165191680982874 , 3.1383331772537426 ,
80+ -3.1383331772537426 , 9.165226363918643 , -9.165226363918643 ,
81+ 15.141616002932361 , -15.141616002932361 , 9.16517861649866 ,
82+ -9.16517861649866 , 15.188854859073416 , -15.188854859073416 ,
83+ 9.165053319552113 , -9.165053319552113 , 15.188703787345304 ,
84+ -15.188703787345304 , 3.138263189805292 , -3.138263189805292 ,
85+ 9.165104089644917 , -9.165104089644917 , 15.141494417823818 ,
86+ -15.141494417823818 , 3.1599158563428644 , -3.1599158563428644 ,
87+ 9.165116302773846 , -9.165116302773846 };
6888
69- Double_t xi[NCellsC] = {gridpoints[1 ], gridpoints[2 ], gridpoints[3 ], gridpoints[4 ], gridpoints[0 ],
70- gridpoints[1 ], gridpoints[2 ], gridpoints[3 ], gridpoints[4 ], gridpoints[5 ],
71- gridpoints[0 ], gridpoints[1 ], gridpoints[4 ], gridpoints[5 ], gridpoints[0 ],
72- gridpoints[1 ], gridpoints[4 ], gridpoints[5 ], gridpoints[0 ], gridpoints[1 ],
73- gridpoints[2 ], gridpoints[3 ], gridpoints[4 ], gridpoints[5 ], gridpoints[1 ],
74- gridpoints[2 ], gridpoints[3 ], gridpoints[4 ]};
75- Double_t yi[NCellsC] = {gridpoints[5 ], gridpoints[5 ], gridpoints[5 ], gridpoints[5 ], gridpoints[4 ],
76- gridpoints[4 ], gridpoints[4 ], gridpoints[4 ], gridpoints[4 ], gridpoints[4 ],
77- gridpoints[3 ], gridpoints[3 ], gridpoints[3 ], gridpoints[3 ], gridpoints[2 ],
78- gridpoints[2 ], gridpoints[2 ], gridpoints[2 ], gridpoints[1 ], gridpoints[1 ],
79- gridpoints[1 ], gridpoints[1 ], gridpoints[1 ], gridpoints[1 ], gridpoints[0 ],
80- gridpoints[0 ], gridpoints[0 ], gridpoints[0 ]};
8189 Double_t zi[NCellsC];
8290 for (Int_t i = 0 ; i < NCellsC; i++) {
8391 zi[i] = TMath::Sqrt (TMath::Power (crad, 2 ) - TMath::Power (xi[i], 2 ) - TMath::Power (yi[i], 2 ));
0 commit comments