SND@LHC Software
Loading...
Searching...
No Matches
Public Member Functions | Private Attributes | List of all members
Floor Class Reference

#include <Floor.h>

Inheritance diagram for Floor:
Inheritance graph
[legend]
Collaboration diagram for Floor:
Collaboration graph
[legend]

Public Member Functions

 Floor (const char *name, Bool_t Active)
 
 Floor ()
 
virtual ~Floor ()
 
void SetConfPar (TString name, Float_t value)
 
void SetConfPar (TString name, Int_t value)
 
void SetConfPar (TString name, TString value)
 
Float_t GetConfParF (TString name)
 
Int_t GetConfParI (TString name)
 
TString GetConfParS (TString name)
 
virtual void Initialize ()
 
virtual Bool_t ProcessHits (FairVolume *v=0)
 
virtual void Register ()
 
virtual TClonesArray * GetCollection (Int_t iColl) const
 
virtual void Reset ()
 
void ConstructGeometry ()
 
virtual void CopyClones (TClonesArray *cl1, TClonesArray *cl2, Int_t offset)
 
virtual void SetSpecialPhysicsCuts ()
 
virtual void EndOfEvent ()
 
virtual void FinishPrimary ()
 
virtual void FinishRun ()
 
virtual void BeginPrimary ()
 
virtual void PostTrack ()
 
virtual void PreTrack ()
 
virtual void BeginEvent ()
 
vetoPoint * AddHit (Int_t trackID, Int_t detID, TVector3 pos, TVector3 mom, Double_t time, Double_t length, Double_t eLoss, Int_t pdgcode, TVector3 Lpos, TVector3 Lmom)
 
void SetEmin (float E)
 
void SetZmax (float Z)
 
void SetFastMuon ()
 
void MakeSensitive ()
 
Int_t InitMedium (const char *name)
 
TVector3 crossing (TVector3 H1, TVector3 H2, TVector3 H3, TVector3 P1, TVector3 P2, TVector3 P3)
 

Private Attributes

std::map< TString, Float_t > conf_floats
 
std::map< TString, Int_t > conf_ints
 
std::map< TString, TString > conf_strings
 
Int_t fTrackID
 
Int_t fVolumeID
 track index
 
TLorentzVector fPos
 volume id
 
TLorentzVector fMom
 position at entrance
 
Double_t fTime
 momentum at entrance
 
Double_t fLength
 time
 
Double_t fzPos
 length
 
Double_t fThick
 zPos
 
Double_t fELoss
 thickness
 
Double_t fTotalEloss
 
TString fMaterial
 
TClonesArray * fFloorPointCollection
 material
 
Double_t SND_Z
 
Bool_t fFastMuon
 
Float_t fEmin
 
Bool_t fMakeSensitive
 

Detailed Description

Definition at line 22 of file Floor.h.

Constructor & Destructor Documentation

◆ Floor() [1/2]

Floor::Floor ( const char *  name,
Bool_t  Active 
)

◆ Floor() [2/2]

Floor::Floor ( )

Definition at line 33 of file Floor.cxx.

34 : FairDetector("Floor", kTRUE, kVETO),
35 fFastMuon(kFALSE),
36 fMakeSensitive(kFALSE),
37 fEmin(0),
38 fTrackID(-1),
39 fVolumeID(-1),
40 fPos(),
41 fMom(),
42 fTime(-1.),
43 fLength(-1.),
44 fThick(-1.),
45 fzPos(3E8),
46 fTotalEloss(0),
47 fFloorPointCollection(new TClonesArray("vetoPoint"))
48{}
kVETO
@ kVETO
Definition ShipDetectorList.h:12
Floor::fTrackID
Int_t fTrackID
Definition Floor.h:85
Floor::fzPos
Double_t fzPos
length
Definition Floor.h:91
Floor::fEmin
Float_t fEmin
Definition Floor.h:100
Floor::fMom
TLorentzVector fMom
position at entrance
Definition Floor.h:88
Floor::fTotalEloss
Double_t fTotalEloss
Definition Floor.h:94
Floor::fTime
Double_t fTime
momentum at entrance
Definition Floor.h:89
Floor::fFastMuon
Bool_t fFastMuon
Definition Floor.h:99
Floor::fMakeSensitive
Bool_t fMakeSensitive
Definition Floor.h:101
Floor::fVolumeID
Int_t fVolumeID
track index
Definition Floor.h:86
Floor::fLength
Double_t fLength
time
Definition Floor.h:90
Floor::fThick
Double_t fThick
zPos
Definition Floor.h:92
Floor::fFloorPointCollection
TClonesArray * fFloorPointCollection
material
Definition Floor.h:97
Floor::fPos
TLorentzVector fPos
volume id
Definition Floor.h:87

◆ ~Floor()

Floor::~Floor ( )
virtual

Definition at line 50 of file Floor.cxx.

51{
52 if (fFloorPointCollection) {
53 fFloorPointCollection->Delete();
54 delete fFloorPointCollection;
55 }
56}

Member Function Documentation

◆ AddHit()

vetoPoint * Floor::AddHit ( Int_t  trackID,
Int_t  detID,
TVector3  pos,
TVector3  mom,
Double_t  time,
Double_t  length,
Double_t  eLoss,
Int_t  pdgcode,
TVector3  Lpos,
TVector3  Lmom 
)

Definition at line 1459 of file Floor.cxx.

1463{
1464 TClonesArray& clref = *fFloorPointCollection;
1465 Int_t size = clref.GetEntriesFast();
1466 LOG(DEBUG) << "add veto point "<<size<<" "<<trackID<<" "<<detID;
1467 return new(clref[size]) vetoPoint(trackID, detID, pos, mom,
1468 time, length, eLoss, pdgCode,Lpos,Lmom);
1469}

◆ BeginEvent()

virtual void Floor::BeginEvent ( )
inlinevirtual

Definition at line 63 of file Floor.h.

63{;}

◆ BeginPrimary()

virtual void Floor::BeginPrimary ( )
inlinevirtual

Definition at line 60 of file Floor.h.

60{;}

◆ ConstructGeometry()

void Floor::ConstructGeometry ( )

from geometer, with z pointing upwards, via localSND2CCS and then CCS2SNDPhysics Double_t M[9] = {0.999978, 0.0000821516, -0.006606, -0.00660651, -0.0124347, 0.999901, -4.69368*10^-15, 0.999923, 0.0124349}

Definition at line 165 of file Floor.cxx.

166{
167 SND_Z = conf_floats["Floor/z"];
168
169 TGeoVolume *top = gGeoManager->GetTopVolume();
170 InitMedium("Concrete");
171 TGeoMedium *concrete = gGeoManager->GetMedium("Concrete");
172 InitMedium("Rock");
173 TGeoMedium *rock = gGeoManager->GetMedium("Rock");
174 InitMedium("vacuum");
175 TGeoMedium *vacuum = gGeoManager->GetMedium("vacuum");
176
177 TGeoVolume *tunnel= new TGeoVolumeAssembly("Tunnel");
178
179 TGeoVolume *detector = new TGeoVolumeAssembly("Detector"); // assembly to place detetector elements in local coordinates
185 // with y-axis upwards, R = ((-1,0,0),(0,0,1),(0,1,0)), M = M R^-1
186 Double_t M[9] = {0.999978, -0.006606, 0.0000821516,
187 0.00660651, 0.999901, -0.0124347,
188 4.69368E-15, 0.0124349, 0.999923};
189
190 auto localSND_physCS_rot = new TGeoRotation("localSND_physCS_rot");
191 localSND_physCS_rot ->SetMatrix(M);
192 auto localSND_physCS_comb = new TGeoCombiTrans("localSND_physCS",0.,0.,0.,localSND_physCS_rot); // origin is 480m downstream of IP1
193 localSND_physCS_comb->RegisterYourself();
194
195 Double_t fTunnelDX = conf_floats["Floor/DX"];
196 Double_t fTunnelDY = conf_floats["Floor/DY"];
197 Double_t fTunnelDZ = conf_floats["Floor/DZ"];
198
199// from Fluka, geo4SND.inp, 25.12.2020
200 std::vector<double> TI18_o1 = {-221.4191473578 , 76.66172460057, 48820.152816717, 131.72892289383, 33.003800302801, -543.2846767937, 200.28612174928};
201 std::vector<double> TI18_i1 = {-221.4191473578 , 76.66172460057, 48820.152816717, 131.72892289383, 33.003800302801, -543.2846767937, 175.0};
202 std::vector<double> TI18_o2 = {-92.57914735775, 108.94172460057, 48288.782816717, 130.14367334394, 56.877589206698, -479.4033418058, 200.0};
203 std::vector<double> TI18_i2 = {-92.57914735775, 108.94172460057, 48288.782816717, 130.14367334394, 56.877589206698, -479.4033418058, 175.0};
204 std::vector<double> TI18_o3 = {14.66085264225, 155.81172460057, 47893.742816717, 1387.8506462618, 818.72395745673, -4733.240081082, 200.0};
205 std::vector<double> TI18_i3 = {14.66085264225, 155.81172460057, 47893.742816717, 1387.8506462618, 818.72395745673, -4733.240081082, 175.0};
206 std::vector<double> UJ18_o = {-471.0, 25.5, 48750.0, -201.8536100583633, 0.0, 3594.3365340638884, 495.0};
207 std::vector<double> UJ18_i = {-471.0,25.5, 48800.39, -201.8536100583633, 0.0, 3494.4938525621, 445.0};
208 std::vector<double> UJ18_x = {-92.57914735775,108.94172460057,48288.782816717,-141.138131672,-35.36121461014,582.09072513608,175.0};
209 std::vector<double> tu012_i = { -596.38194369, 30., 47654.730979, -155.10509191, 0.0,2761.9020453, 220.0};
210 std::vector<double> tu012_o = { -596.38194369, 30., 4.7654730979E+04, -155.10509191, 0.0,2761.9020453, 250.0};
211 std::vector<double> tu011_i = { -522.09173359, 30., 4.6171636594E+04, -99.7891926, 0.0,1954.9587043, 220.0};
212 std::vector<double> tu011_o = { -522.09173359, 30., 4.6171636594E+04, -99.7891926, 0.0,1954.9587043, 250.0,360.,850.};
213 std::vector<double> tu010_i = { -451.24339295, 30., 4.4607412772E+04, -89.817671994, 0.0,1955.4422017, 220.0};
214 std::vector<double> tu010_o = { -451.24339295, 30., 4.4607412772E+04, -89.817671994, 0.0,1955.4422017, 250.0,750.,1100.};
215// PLA TI18_01 and TI18_02
216 std::vector<double> TI18_1_bot = { 0.67644808637366, -74.85571453913, -4.530940169849, -221.4171473578, -18.33127539943,48818.967816717,
217 0.63027909276104 ,-74.86684394728, -4.35010639392, -221.4171473578, -93.32587539943,48818.036016717};
218// side PLA TI18_08 and PLA TI18_09
219 std::vector<double> TI18_1_sid = { -95.08832719702, 21.886725033285, -21.89021009926, -39.11837578468,3.7557255037298,48858.930618475,
220 -94.65075310463, -21.38126814052, -24.16767074894,-402.3825398293, 5.6093808446198,48770.79950338};
221// PLA TI18_04 and PLA TI18_03
222 std::vector<double> TI18_2_bot = { 1.590393176461, -74.48871190132, -8.55964069504, -92.6293750264, 13.94692460057, 48287.934643717,
223 1.12692975448, -74.65145228529, -7.089125547405, -92.57914735775, -61.04827539943, 48286.672816717};
224// side PLA TI18_81 and PLA TI18_91
225 std::vector<double> TI18_2_sid = { 9.5229226656561, -2.052234306283, 2.2588218587662, 119.6722824853, -43.26523379053, 48132.959403311,
226 0.9420958238706, 0.1874442034604, 0.2780649358884, -220.8962283082, 57.49004938073, 48051.134027534};
227// PLA TI18_05 and TI18_06
228 std::vector<double> TI18_3_bot = { 2.6719179804842, -74.62099192635, -12.16834081665, 237.27547988904, 188.51303030346, 47133.251503475,
229 2.2386457107868, -74.59695305235, -12.40093386351, 14.844992780245, -22.78799545647, 47890.126221505};
230// side PLA TI18_08x and PLA TI18_09x
231 std::vector<double> TI18_3_sid = { 9.4281652604622, -2.374007701455, 2.3396124281475, 399.46744859669, 116.81598987751, 47153.021981939,
232 9.332746336299, 1.9001431836442, 3.0478354225656, 88.79276895539223, 130.4393658272924, 47064.80229883217};
233 std::vector<double> PM18_ci = {50660.2873, -1256.8854,700.0,500.};
234 std::vector<double> PM18_co = {50660.2873, -1256.8854, 750.,500.};
235 std::map<TString,std::vector<double>> geoParameters = {
236 {"TI18_o1",TI18_o1},
237 {"TI18_i1",TI18_i1},
238 {"TI18_o2",TI18_o2},
239 {"TI18_i2",TI18_i2},
240 {"TI18_o3",TI18_o3},
241 {"TI18_i3",TI18_i3},
242 {"UJ18_o",UJ18_o},
243 {"UJ18_i",UJ18_i},
244 {"tu010_o",tu010_o},
245 {"tu010_i",tu010_i},
246 {"tu011_o",tu011_o},
247 {"tu011_i",tu011_i},
248 {"tu012_o",tu012_o},
249 {"tu012_i",tu012_i},
250 {"UJ18_x",UJ18_x},
251 {"PM18_ci",PM18_ci},
252 {"PM18_co",PM18_co},
253 {"TI18_1_bot",TI18_1_bot},
254 {"TI18_1_sid",TI18_1_sid},
255 {"TI18_2_bot",TI18_2_bot},
256 {"TI18_2_sid",TI18_2_sid},
257 {"TI18_3_bot",TI18_3_bot},
258 {"TI18_3_sid",TI18_3_sid}
259 };
260 std::vector<TString> loop1 = {"TI18_1","TI18_2","TI18_3"};
261 std::vector<TGeoCompositeShape*> shapes;
262 std::vector<TGeoCombiTrans*> trans;
263 for(auto it = std::begin(loop1); it != std::end(loop1); ++it) {
264 LOG(DEBUG)<<"Floor: DEBUG start" << *it;
265 auto name = TString(*it);
266 TString name_i = TString(*it).ReplaceAll("_","_i");
267 TString name_o = TString(*it).ReplaceAll("_","_o");
268 TString name_bot = name+"_bot";
269 TString name_sid = name+"_sid";
270 LOG(DEBUG)<<"Floor: " << geoParameters[name_i.Data()][0]<< " "<< geoParameters[name_i.Data()][1]<< " "<< geoParameters[name_i.Data()][2]<<" "<<name_i;
271 TVector3 P = TVector3(geoParameters[name_i.Data()][0],geoParameters[name_i.Data()][1],geoParameters[name_i.Data()][2] - SND_Z);
272 TVector3 H = TVector3(geoParameters[name_i.Data()][3],geoParameters[name_i.Data()][4],geoParameters[name_i.Data()][5]);
273 auto T = P + 0.5*H;
274 double rotBeta = TMath::ACos(H[2]/H.Mag());
275 double rotAlpha = TMath::ATan2(H[0]/H.Mag(), -H[1]/H.Mag());
276 double rotGamma = TMath::ACos(H[1]/H.Mag()/TMath::Sqrt(1-pow(H[2]/H.Mag(),2)));
277 auto R = new TGeoRotation("R_"+name,rotAlpha*180./TMath::Pi(),rotBeta*180./TMath::Pi(),rotGamma*180./TMath::Pi());
278 auto CombiTrans = new TGeoCombiTrans("T_"+name,T[0],T[1],T[2],R);
279 CombiTrans->RegisterYourself();
280 auto X = geoParameters[name_bot.Data()];
281 TVector3 Ptop(X[3],X[4], X[5] -SND_Z);
282 TVector3 Pbot(X[9],X[10], X[11]-SND_Z);
283 TVector3 Htop(X[0], X[1],X[2]);
284 TVector3 Hbot(X[6], X[7],X[8]);
285 X = geoParameters[name_sid.Data()];
286 TVector3 Pleft(X[9],X[10], X[11] -SND_Z);
287 TVector3 Prigh(X[3],X[4], X[5]-SND_Z);
288 TVector3 Hleft(X[6], X[7],X[8]);
289 TVector3 Hrigh(X[0], X[1],X[2]);
290 std::map<TString,TVector3> face = {
291 {"front",P},{"back",P+H}
292 };
293 std::vector<TVector3> points;
294 std::vector<TVector3> GlobalPoints;
295 Double_t master[3], P1[3], P2[3], P3[3], P4[3];
296#
297 rotBeta = TMath::ACos(Htop[2]/Htop.Mag());
298 rotAlpha = TMath::ATan2(Htop[0]/H.Mag(), -Htop[1]/H.Mag()) ;
299 rotGamma = TMath::ACos(Htop[1]/Htop.Mag()/TMath::Sqrt(1-pow(Htop[2]/Htop.Mag(),2)));
300 auto RF = new TGeoRotation( "R_tube_"+name,rotAlpha*180./TMath::Pi(),rotBeta*180./TMath::Pi(),rotGamma*180./TMath::Pi());
301 Double_t dY = (Htop.Dot(Ptop)-Htop.Dot(Pbot))/Htop.Mag();
302 std::vector<TString> loopF = {"front","back"};
303 for(auto il = std::begin(loopF); il != std::end(loopF); ++il) {
304 auto X0 = crossing(Hleft,Htop,H,Pleft,Ptop,face[*il]);
305 for( int i = 0;i<3;i+=1) { master[i]=X0[i];}
306 GlobalPoints .push_back(X0);
307 RF->MasterToLocal(master,P1);
308 points.push_back(TVector3(P1[0],P1[1],P1[2]));
309 LOG(DEBUG) << "top left "<<X0[0]<<","<<X0[1]<<","<<X0[2]<<"; "<<P1[0]<<","<<P1[1]<<","<<P1[2];
310 X0 = crossing(Hrigh,Htop,H,Prigh,Ptop,face[*il]);
311 for( int i = 0;i<3;i+=1) { master[i]=X0[i];}
312 GlobalPoints .push_back(X0);
313 RF->MasterToLocal(master,P2);
314 points.push_back(TVector3(P2[0],P2[1],P2[2]));
315 LOG(DEBUG) << "top right "<<X0[0]<<","<<X0[1]<<","<<X0[2]<<"; "<<P2[0]<<","<<P2[1]<<","<<P2[2];
316 X0 = crossing(Hrigh,Htop,H,Prigh,Pbot,face[*il]);
317 for( int i = 0;i<3;i+=1) { master[i]=X0[i];}
318 GlobalPoints .push_back(X0);
319 RF->MasterToLocal(master,P3);
320 points.push_back(TVector3(P3[0],P3[1],P3[2]));
321 LOG(DEBUG) << "bottom right "<<X0[0]<<","<<X0[1]<<","<<X0[2]<<"; "<<P3[0]<<","<<P3[1]<<","<<P3[2];
322 X0 = crossing(Hleft,Htop,H,Pleft,Pbot,face[*il]);
323 for( int i = 0;i<3;i+=1) { master[i]=X0[i];}
324 GlobalPoints .push_back(X0);
325 RF->MasterToLocal(master,P4);
326 points.push_back(TVector3(P4[0],P4[1],P4[2]));
327 LOG(DEBUG) << "bottom left "<<X0[0]<<","<<X0[1]<<","<<X0[2]<<"; "<<P4[0]<<","<<P4[1]<<","<<P4[2];
328 }
329
330 Double_t dZ = (P3[2]-P1[2])/2.;
331 auto arb = new TGeoArb8(name+"_arb8",dZ);
332 arb->SetVertex(0, points[1][0], points[1][1]);
333 arb->SetVertex(1, points[5][0], points[5][1]);
334 arb->SetVertex(2, points[4][0], points[4][1]);
335 arb->SetVertex(3, points[0][0], points[0][1]);
336 arb->SetVertex(4, points[2][0], points[2][1]);
337 arb->SetVertex(5, points[6][0], points[6][1]);
338 arb->SetVertex(6, points[7][0], points[7][1]);
339 arb->SetVertex(7, points[3][0], points[3][1]);
340// to get the transformation
341 Double_t local[3] = {points[1][0],points[1][1],-dZ};
342 RF->LocalToMaster(local,master);
343 TVector3 Tf(GlobalPoints[1][0]-master[0],GlobalPoints[1][1]-master[1],GlobalPoints[1][2]-master[2]);
344 auto CombiTransF = new TGeoCombiTrans("T_floor_"+name,Tf[0],Tf[1],Tf[2],RF);
345 CombiTransF->RegisterYourself();
346 auto tube = gGeoManager->MakeTube(name+"_tube", concrete, geoParameters[name_i.Data()][6], geoParameters[name_o.Data()][6],H.Mag()/2.);
347 auto Fulltube = gGeoManager->MakeTube(name+"_fulltube", concrete, 0., geoParameters[name_o.Data()][6],H.Mag()/2.);
348 LOG(DEBUG)<<name<<" "<<Tf[0]<<" "<<Tf[1]<<" "<<Tf[2];
349 auto CombiTransF1 = new TGeoCombiTrans("T_floorS1_"+name,Tf[0],Tf[1]-2*dZ,Tf[2],RF);
350 auto CombiTransF2 = new TGeoCombiTrans("T_floorS2_"+name,Tf[0],Tf[1]-2*dZ,Tf[2]-20,RF);
351 CombiTransF1->RegisterYourself();
352 CombiTransF2->RegisterYourself();
353 auto sunion = new TGeoCompositeShape(name+"_union","("+name+"_tube:T_"+name+"+"+name+"_arb8:T_floor_"+name+"-"+
354 name+"_arb8:T_floorS1_"+name+"-"+name+"_arb8:T_floorS2_"+name+")");
355 shapes.push_back(sunion);
356 auto Fsunion = new TGeoCompositeShape(name+"_Funion","("+name+"_fulltube:T_"+name+"+"+name+"_arb8:T_floor_"+name+"-"+
357 name+"_arb8:T_floorS1_"+name+"-"+name+"_arb8:T_floorS2_"+name+")");
358 shapes.push_back(Fsunion);
359
360 }
361 LOG(DEBUG) << "shapes "<<shapes[0]->GetName()<<" "<<shapes[1]->GetName()<<" "<<shapes[2]->GetName();
362
363 // Since 2024 there is a pit on the TI18 floor hosting the veto
364 // Unless defined in the geo config file, the pit dims are 0.
365 auto vetoPit = new TGeoBBox("vetoPit",
366 conf_floats["Floor/VetoPitXdim"]/2.,
367 conf_floats["Floor/VetoPitZdim"]/2.,
368 conf_floats["Floor/VetoPitYdim"]/2.);
369 auto VetoPit_transl = new TGeoTranslation("VetoPit_transl",
370 -conf_floats["Floor/VetoPitX"]-conf_floats["Floor/VetoPitXdim"]/2.,
371 conf_floats["Floor/VetoPitZ"]-conf_floats["Floor/VetoPitZdim"]/2.,
372 conf_floats["Floor/VetoPitY"]-conf_floats["Floor/VetoPitYdim"]/2.);
373 VetoPit_transl->RegisterYourself();
374
375 //UJ18 tunnel
376 // UJ18 translation and angle rotation
377 Double_t dx_0_UJ18 = -2300.0;
378 Double_t dy_0_UJ18 = -273.0;
379 Double_t dz_0_UJ18 = 7346.0;
380 Double_t rot_X = 1.3;
381 Double_t rot_Y = -13.63;
382 Double_t rot_Z = 0;
383
384 TGeoRotation *rotA = new TGeoRotation("rotA");
385 rotA->RotateX(rot_X);
386 rotA->RotateY(rot_Y);
387 rotA->RotateZ(rot_Z);
388 TGeoCombiTrans *matA = new TGeoCombiTrans("matA",dx_0_UJ18, dy_0_UJ18, dz_0_UJ18, rotA);
389
390 Double_t dx_UJ18, dy_UJ18, dz_UJ18;
391 Double_t dx_UJ181, dx_UJ182, dy_UJ181, dy_UJ182;
392 Double_t vert[20], par[20];
393 Double_t theta, phi, h1, bl1, tl1, alpha1, h2, bl2, tl2, alpha2;
394 Double_t twist;
395 Double_t origin[3];
396 Double_t rmin, rmax, rmin1, rmax1, rmin2, rmax2;
397 Double_t r, rlo, rhi;
398 Double_t a, b;
399 Double_t point[3], norm[3];
400 Double_t rin, stin, rout, stout;
401 Double_t thx, phx, thy, phy, thz, phz;
402 Double_t alpha, theta1, theta2, phi1, phi2, dphi;
403 Double_t tr[3], rot[9];
404 Double_t z, density, radl, absl, w;
405 Double_t lx, ly, lz, tx, ty, tz;
406 Double_t xvert[50], yvert[50];
407 Double_t zsect, x0, y0, scale0;
408 Int_t nel, numed, nz, nedges, nvert;
409 TGeoBoolNode *Boolean = nullptr;
410
411 // UJ18 construction: Sequential naming represents step-by-step Boolean operations (e.g., Cut_001, Add_001, trans_1).
412 rmin = 0;
413 rmax = 506.18;
414 dz_UJ18 = 30;
415 TGeoShape *GDMLTube_001 = new TGeoTube("GDMLTube_001",rmin,rmax,dz_UJ18);
416
417 dx_UJ18 = 600;
418 dy_UJ18 = 200;
419 dz_UJ18 = 40;
420 TGeoShape *GDMLBox_001 = new TGeoBBox("GDMLBox_001", dx_UJ18,dy_UJ18,dz_UJ18);
421
422 dx_UJ18 = 0;
423 dy_UJ18 = -524.94;
424 dz_UJ18 = 0;
425 auto trans_1 = new TGeoCombiTrans("trans_1"); // trans = translation
426 trans_1->SetTranslation(dx_UJ18, dy_UJ18, dz_UJ18);
427 Boolean = new TGeoSubtraction(GDMLTube_001,GDMLBox_001,0,trans_1);
428
429 TGeoShape *Cut_001 = new TGeoCompositeShape("Cut_001", Boolean); //Cut = sustraction in boolean operation
430
431 rmin = 0;
432 rmax = 506.18;
433 dz_UJ18 = 1799.5;
434 TGeoShape *GDMLTube_002 = new TGeoTube("GDMLTube_002",rmin,rmax,dz_UJ18);
435
436 rmin = 0;
437 rmax = 450;
438 dz_UJ18 = 1950;
439 TGeoShape *GDMLTube_003 = new TGeoTube("GDMLTube_003",rmin,rmax,dz_UJ18);
440
441 dx_UJ18 = 500;
442 dy_UJ18 = 200;
443 dz_UJ18 = 2000;
444 TGeoShape *GDMLBox_002 = new TGeoBBox("GDMLBox_002", dx_UJ18,dy_UJ18,dz_UJ18);
445
446 dx_UJ18 = 0;
447 dy_UJ18 = -257.94;
448 dz_UJ18 = 0;
449 auto trans_2 = new TGeoCombiTrans("trans_2");
450 trans_2->SetTranslation(dx_UJ18, dy_UJ18, dz_UJ18);
451 Boolean = new TGeoSubtraction(GDMLTube_003,GDMLBox_002,0,trans_2);
452
453 TGeoShape *Cut_002 = new TGeoCompositeShape("Cut_002", Boolean);
454 Boolean = new TGeoSubtraction(GDMLTube_002,Cut_002,0,0);
455
456 TGeoShape *Cut_003 = new TGeoCompositeShape("Cut_003", Boolean);
457
458 dx_UJ18 = 600;
459 dy_UJ18 = 200;
460 dz_UJ18 = 2000;
461 TGeoShape *GDMLBox_003 = new TGeoBBox("GDMLBox_003", dx_UJ18,dy_UJ18,dz_UJ18);
462
463 dx_UJ18 = 0;
464 dy_UJ18 = -524.94;
465 dz_UJ18 = 0;
466 auto trans_3 = new TGeoCombiTrans("trans_3");
467 trans_3->SetTranslation(dx_UJ18, dy_UJ18, dz_UJ18);
468 Boolean = new TGeoSubtraction(Cut_003,GDMLBox_003,0,trans_3);
469
470 TGeoShape *Cut_004 = new TGeoCompositeShape("Cut_004", Boolean);
471
472 rmin = 0;
473 rmax = 506.18;
474 dz_UJ18 = 30;
475 TGeoShape *GDMLTube_004 = new TGeoTube("GDMLTube_004",rmin,rmax,dz_UJ18);
476
477 dx_UJ18 = 600;
478 dy_UJ18 = 200;
479 dz_UJ18 = 40;
480 TGeoShape *GDMLBox_004 = new TGeoBBox("GDMLBox_004", dx_UJ18,dy_UJ18,dz_UJ18);
481
482 dx_UJ18 = 0;
483 dy_UJ18 = -524.94;
484 dz_UJ18 = 0;
485 auto trans_4 = new TGeoCombiTrans("trans_4");
486 trans_4->SetTranslation(dx_UJ18, dy_UJ18, dz_UJ18);
487 Boolean = new TGeoSubtraction(GDMLTube_004,GDMLBox_004,0,trans_4);
488
489 TGeoShape *Cut_005 = new TGeoCompositeShape("Cut_005", Boolean);
490
491 dx_UJ18 = 0;
492 dy_UJ18 = 0;
493 dz_UJ18 = -1829.5;
494 auto trans_5 = new TGeoCombiTrans("trans_5");
495 trans_5->SetTranslation(dx_UJ18, dy_UJ18, dz_UJ18);
496 Boolean = new TGeoUnion(Cut_004,Cut_005,0,trans_5);
497
498 TGeoShape *Add_001 = new TGeoCompositeShape("Add_001", Boolean); // Add = union in boolean operation
499
500 dx_UJ18 = 0;
501 dy_UJ18 = 0;
502 dz_UJ18 = -1829.5;
503 auto trans_6 = new TGeoCombiTrans("trans_6");
504 trans_6->SetTranslation(dx_UJ18, dy_UJ18, dz_UJ18);
505 Boolean = new TGeoUnion(Cut_001,Add_001,0,trans_6);
506
507 TGeoShape *Add_002 = new TGeoCompositeShape("Add_002", Boolean);
508
509 rmin = 0;
510 rmax = 175.65;
511 dz_UJ18 = 50;
512 TGeoShape *GDMLTube_005 = new TGeoTube("GDMLTube_005",rmin,rmax,dz_UJ18);
513
514 dx_UJ18 = 250;
515 dy_UJ18 = 100;
516 dz_UJ18 = 60;
517 TGeoShape *GDMLBox_005 = new TGeoBBox("GDMLBox_005", dx_UJ18,dy_UJ18,dz_UJ18);
518
519 dx_UJ18 = 0;
520 dy_UJ18 = -196.2;
521 dz_UJ18 = 0;
522 auto trans_7 = new TGeoCombiTrans("trans_7");
523 trans_7->SetTranslation(dx_UJ18, dy_UJ18, dz_UJ18);
524 Boolean = new TGeoSubtraction(GDMLTube_005,GDMLBox_005,0,trans_7);
525
526 TGeoShape *Cut_006 = new TGeoCompositeShape("Cut_006", Boolean);
527
528 dx_UJ18 = 245.426;
529 dy_UJ18 = 106.189;
530 dz_UJ18 = -3658.834;
531 auto trans_8 = new TGeoCombiTrans("trans_8");
532 trans_8->SetTranslation(dx_UJ18, dy_UJ18, dz_UJ18);
533 Boolean = new TGeoSubtraction(Add_002,Cut_006,0,trans_8);
534
535 TGeoShape *Cut_007 = new TGeoCompositeShape("Cut_007", Boolean);
536
537 rmin = 0;
538 rmax = 220.39;
539 dz_UJ18 = 1900;
540 TGeoShape *GDMLTube_006 = new TGeoTube("GDMLTube_006",rmin,rmax,dz_UJ18);
541
542 dx_UJ18 = 250;
543 dy_UJ18 = 100;
544 dz_UJ18 = 1925;
545 TGeoShape *GDMLBox_006 = new TGeoBBox("GDMLBox_006", dx_UJ18,dy_UJ18,dz_UJ18);
546
547 dx_UJ18 = 0;
548 dy_UJ18 = -224.647;
549 dz_UJ18 = 0;
550 auto trans_9 = new TGeoCombiTrans("trans_9");
551 trans_9->SetTranslation(dx_UJ18, dy_UJ18, dz_UJ18);
552 Boolean = new TGeoSubtraction(GDMLTube_006,GDMLBox_006,0,trans_9);
553
554 TGeoShape *Cut_008 = new TGeoCompositeShape("Cut_008", Boolean);
555
556 dx_UJ18 = -193.584;
557 dy_UJ18 = 55.7;
558 dz_UJ18 = -1829.5;
559 auto trans_10 = new TGeoCombiTrans("trans_10");
560 trans_10->SetTranslation(dx_UJ18, dy_UJ18, dz_UJ18);
561 Boolean = new TGeoSubtraction(Cut_007,Cut_008,0,trans_10);
562
563 TGeoShape *Cut_009 = new TGeoCompositeShape("Cut_009", Boolean);
564
565 dz_UJ18 = 473.8;
566 theta = 0;
567 phi = 0;
568 h1 = 37.375;
569 bl1 = 204.985;
570 tl1 = 204.985;
571 alpha1 = 0;
572 h2 = 0;
573 bl2 = 204.985;
574 tl2 = 204.985;
575 alpha2 = 0;
576 TGeoShape *GDMLTrap_1 = new TGeoTrap("GDMLTrap_1", dz_UJ18,theta,phi,h1,bl1,tl1,alpha1,h2,bl2,tl2,alpha2);
577
578 dx_UJ18 = 236.991;
579 dy_UJ18 = -45.068;
580 dz_UJ18 = -3156.224;
581
582 thx = 90; phx = 0;
583 thy = 87.92; phy = 90;
584 thz = 2.08; phz = 270;
585 TGeoRotation *rot_1 = new TGeoRotation("rot_1",thx,phx,thy,phy,thz,phz);
586 auto trans_11 = new TGeoCombiTrans("trans_11", dx_UJ18, dy_UJ18, dz_UJ18, rot_1);
587 Boolean = new TGeoUnion(Cut_009,GDMLTrap_1,0,trans_11);
588 TGeoShape *pUJ18_shape_25 = new TGeoCompositeShape("UJ18_shape", Boolean);
589
590 //get the new rotation and translation
591 TGeoHMatrix * Matrix_setup = new TGeoHMatrix("Matrix_setup");
592 Matrix_setup->Multiply(matA);
593 Matrix_setup->Multiply(trans_11);
594
595 //create the new translation for union(with TI18) subtraction(with rock)
596 const Double_t *tr1 = Matrix_setup->GetTranslation();
597 TGeoRotation *rot2 = new TGeoRotation();
598 rot2->SetMatrix(Matrix_setup->GetRotationMatrix());
599 //for registyourself
600 auto new_translation_for_UJ18 = new TGeoCombiTrans("new_translation_for_UJ18", tr1[0], tr1[1], tr1[2], rot2);
601 new_translation_for_UJ18->RegisterYourself();
602
603 // UJ18_solid (to be subtracted from Bigbox)
604 rmin = 0;
605 rmax = 506.18;
606 dz_UJ18 = 30;
607 TGeoShape *GDMLTube_007 = new TGeoTube("GDMLTube_007",rmin,rmax,dz_UJ18);
608
609 dx_UJ18 = 600;
610 dy_UJ18 = 200;
611 dz_UJ18 = 40;
612 TGeoShape *GDMLBox_007 = new TGeoBBox("GDMLBox_007", dx_UJ18,dy_UJ18,dz_UJ18);
613
614 dx_UJ18 = 0;
615 dy_UJ18 = -524.94;
616 dz_UJ18 = 0;
617 auto trans_12 = new TGeoCombiTrans("trans_12");
618 trans_12->SetTranslation(dx_UJ18, dy_UJ18, dz_UJ18);
619 Boolean = new TGeoSubtraction(GDMLTube_007,GDMLBox_007,0,trans_12);
620
621 TGeoShape *Cut_010 = new TGeoCompositeShape("Cut_010", Boolean);
622
623 rmin = 0;
624 rmax = 506.18;
625 dz_UJ18 = 1799.5;
626 TGeoShape *GDMLTube_008 = new TGeoTube("GDMLTube_8",rmin,rmax,dz_UJ18);
627
628 dx_UJ18 = 600;
629 dy_UJ18 = 200;
630 dz_UJ18 = 2000;
631 TGeoShape *GDMLBox_008 = new TGeoBBox("GDML_Box008", dx_UJ18,dy_UJ18,dz_UJ18);
632
633 dx_UJ18 = 0;
634 dy_UJ18 = -524.94;
635 dz_UJ18 = 0;
636 auto trans_13 = new TGeoCombiTrans("trans_13");
637 trans_13->SetTranslation(dx_UJ18, dy_UJ18, dz_UJ18);
638 Boolean = new TGeoSubtraction(GDMLTube_008,GDMLBox_008,0,trans_13);
639
640 TGeoShape *Cut_011 = new TGeoCompositeShape("Cut_011", Boolean);
641
642 rmin = 0;
643 rmax = 506.18;
644 dz_UJ18 = 30;
645 TGeoShape *GDMLTube_009 = new TGeoTube("GDMLTube_009",rmin,rmax,dz_UJ18);
646
647 dx_UJ18 = 600;
648 dy_UJ18 = 200;
649 dz_UJ18 = 40;
650 TGeoShape *GDMLBox_009 = new TGeoBBox("GDMLBox_009", dx_UJ18,dy_UJ18,dz_UJ18);
651
652 dx_UJ18 = 0;
653 dy_UJ18 = -524.94;
654 dz_UJ18 = 0;
655 auto trans_14 = new TGeoCombiTrans("trans_14");
656 trans_14->SetTranslation(dx_UJ18, dy_UJ18, dz_UJ18);
657 Boolean = new TGeoSubtraction(GDMLTube_009,GDMLBox_009,0,trans_14);
658
659 TGeoShape *Cut_012 = new TGeoCompositeShape("Cut_012", Boolean);
660
661 dx_UJ18 = 0;
662 dy_UJ18 = 0;
663 dz_UJ18 = -1829.5;
664 auto trans_15 = new TGeoCombiTrans("trans_15");
665 trans_15->SetTranslation(dx_UJ18, dy_UJ18, dz_UJ18);
666 Boolean = new TGeoUnion(Cut_011,Cut_012,0,trans_15);
667
668 TGeoShape *Add_003 = new TGeoCompositeShape("Add_003", Boolean);
669
670 dx_UJ18 = 0;
671 dy_UJ18 = 0;
672 dz_UJ18 = -1829.5;
673 auto trans_16 = new TGeoCombiTrans("trans_16");
674 trans_16->SetTranslation(dx_UJ18, dy_UJ18, dz_UJ18);
675 Boolean = new TGeoUnion(Cut_010,Add_003,0,trans_16);
676 TGeoShape *Solid_UJ18_shape = new TGeoCompositeShape("Solid_UJ18_shape", Boolean);
677
678 std::string tunnel_shape_name = "(UJ18_shape:new_translation_for_UJ18+TI18_1_union+TI18_2_union+TI18_3_union)";
679 if (conf_floats["Floor/VetoPitXdim"]!=0) {
680 tunnel_shape_name+="-vetoPit:VetoPit_transl";
681 }
682 auto total = new TGeoCompositeShape("Stotal",tunnel_shape_name.c_str());
683 auto volT = new TGeoVolume("VTI18",total,concrete);
684 volT->SetTransparency(0);
685 volT->SetLineColor(kGray);
686// make tunnel sensitive for debugging
687 if (fMakeSensitive) {AddSensitiveVolume(volT);}
688 tunnel->AddNode(volT, 1);
689
690 std::vector<TString> loop2 = {"UJ18_","tu010_","tu011_","tu012_"};
691 for(auto it = std::begin(loop2); it != std::end(loop2); ++it) {
692 TString name_i = TString(*it).ReplaceAll("_","_i");
693 TString name_o= TString(*it).ReplaceAll("_","_o");
694 TString name_F = "F"+TString(*it);
695 TVector3 P = TVector3(geoParameters[name_i.Data()][0],geoParameters[name_i.Data()][1],geoParameters[name_i.Data()][2] - SND_Z);
696 TVector3 H = TVector3(geoParameters[name_i.Data()][3],geoParameters[name_i.Data()][4],geoParameters[name_i.Data()][5]);
697 auto vol = gGeoManager->MakeTube(*it, concrete, geoParameters[name_i.Data()][6], geoParameters[name_o.Data()][6],H.Mag()/2.-5.); // subtract 5cm to avoid overlap
698 auto Fvol = gGeoManager->MakeTube(name_F, concrete, 0., geoParameters[name_o.Data()][6],H.Mag()/2.-5.);
699 vol->SetTransparency(0);
700 vol->SetLineColor(kGray);
701 auto T = P + 0.5*H;
702 double rotBeta = TMath::ACos(H[2]/H.Mag());
703 double rotAlpha = TMath::ATan2(H[0]/H.Mag(), -H[1]/H.Mag());
704 double rotGamma = TMath::ACos(H[1]/H.Mag()/TMath::Sqrt(1-pow(H[2]/H.Mag(),2)));
705 TString name = "R_"+TString(*it);
706 auto R = new TGeoRotation(name,rotAlpha*180./TMath::Pi(),rotBeta*180./TMath::Pi(),rotGamma*180./TMath::Pi());
707 auto CombiTrans = new TGeoCombiTrans("TS_"+TString(*it),T[0],T[1],T[2],R);
708 CombiTrans->RegisterYourself();
709 trans.push_back(CombiTrans);
710 }
711 auto total2 = new TGeoCompositeShape("Stotal2","tu010_:TS_tu010_+tu011_:TS_tu011_");
712 auto Ftotal2 = new TGeoCompositeShape("Ftotal2","Ftu010_:TS_tu010_+Ftu011_:TS_tu011_");
713 auto volT2 = new TGeoVolume("VUJ",total2,concrete);
714 volT2->SetTransparency(0);
715 volT2->SetLineColor(kGray);
716 tunnel->AddNode(volT2, 1);
717
718/*
719 std::vector<TString> loopr = {"tu010_","tu011_"};
720 for(auto it = std::begin(loopr); it != std::end(loopr); ++it) {
721 TString name_o= TString(*it).ReplaceAll("_","_o");
722 TString name_b = TString(*it).ReplaceAll("_","_cone");
723 TString name_c = TString(*it).ReplaceAll("_","_rock");
724 TString name = TString(*it);
725 TString name_x = TString(*it).ReplaceAll("_","_x");
726
727 TVector3 H = TVector3(geoParameters[name_o.Data()][3],geoParameters[name_o.Data()][4],geoParameters[name_o.Data()][5]);
728 auto tube = new TGeoTube(name_x,0., geoParameters[name_o.Data()][6]+10.,H.Mag()/2.+0.1);
729 auto cone = gGeoManager->MakeCone(name_b,rock,H.Mag()/2.,0.,geoParameters[name_o.Data()][8],0.,geoParameters[name_o.Data()][7]);
730 auto Srock = new TGeoCompositeShape(name_c, name_b+":TS_"+name+"-"+name_x+":TS_"+name);
731 }
732 auto total3 = new TGeoCompositeShape("Stotal3","tu010_rock+tu011_rock");
733 auto volT3 = new TGeoVolume("Vrock",total3,rock);
734 volT3->SetTransparency(75);
735 volT3->SetLineColor(kRed);
736 tunnel->AddNode(volT3, 1);
737
738// add rock up to Fluka scoring plane, 409 m from IP1.
739 double zs = 40900.;
740 double dz = (geoParameters["tu010_o"][2] - zs)/2.;
741 auto fluka = gGeoManager->MakeBox("FlukaScoringPlane",rock, geoParameters["tu010_o"][8],geoParameters["tu010_o"][8],dz);
742 fluka->SetLineColor(kRed);
743 tunnel->AddNode(fluka,1, new TGeoTranslation(-350.,0,dz+zs- SND_Z-50.)); // move 50cm upstream to avoid overlap
744*/
745
746 double zs = 40900.; // scoring plane
747 double dz = (geoParameters["TI18_o1"][2] - zs)/2.;
748 double extended_z = 1860; //extend the rock to cover the new added UJ18
749
750 auto bigBox = new TGeoBBox("BigBox", 2000.,1500. , dz+extended_z);
751 auto TR_1 = new TGeoTranslation("TR_1",0.,0.,-dz+extended_z+geoParameters["TI18_o1"][2]-SND_Z - 50.); // move a bit more upstream to have free view from the back
752 TR_1->RegisterYourself();
753
754 auto cutOut = new TGeoCompositeShape("cutOut", "BigBox:TR_1-Ftotal2-(TI18_1_Funion+TI18_2_Funion+TI18_3_Funion+Solid_UJ18_shape:new_translation_for_UJ18)");
755 auto volT3 = new TGeoVolume("Vrock",cutOut,rock);
756 volT3->SetTransparency(75);
757 volT3->SetLineColor(kRed);
758 tunnel->AddNode(volT3, 1);
759
760 double z_dim_exitPlane = 1.;
761auto exitPlane = gGeoManager->MakeBox("exitScoringPlane",vacuum,1000.,1000. , z_dim_exitPlane);
762exitPlane->SetLineColor(kGreen);
763exitPlane->SetVisibility(kFALSE);
764if (fMakeSensitive) {AddSensitiveVolume(exitPlane);}
765tunnel->AddNode(exitPlane,1, new TGeoTranslation(-1500.0,0,((2*extended_z)+geoParameters["TI18_o1"][2]-SND_Z - 50 + z_dim_exitPlane)));
766
767 // COLDBOX GEOM
768 InitMedium("Borated30polyethylene");
769 TGeoMedium *Bor30Poly = gGeoManager->GetMedium("Borated30polyethylene");
770 InitMedium("PMMA");
771 TGeoMedium *Acrylic = gGeoManager->GetMedium("PMMA");
772 TVector3 displacement;
773
774 Double_t fAcrylicWidth = conf_floats["Floor/Acrylic_width"];
775 Double_t fBPolyWidth = conf_floats["Floor/BPoly_width"];
776 Double_t fCBFrontWallXDim = conf_floats["Floor/CBFrontWall_xdim"];
777 Double_t fCBFrontWallYDim = conf_floats["Floor/CBFrontWall_ydim"];
778 Double_t fCBLatWallZDim = conf_floats["Floor/CBLatWall_zdim"];
779 Double_t fCBTinyZDim = conf_floats["Floor/CBTiny_zdim"];
780 Double_t fCBExtraZDim = conf_floats["Floor/CBExtra_zdim"];
781 Double_t fCBExtraXDim = conf_floats["Floor/CBExtra_xdim"];
782 Double_t fSlopedWallZProj = conf_floats["Floor/SlopedWall_zproj"];
783 Double_t fCBRearWallXDim =
784 fCBFrontWallXDim - fSlopedWallZProj * TMath::Tan(TMath::DegToRad() * 15.) - fCBExtraXDim + fAcrylicWidth;
785 Double_t fFeBlockX = conf_floats["Floor/MFeBlockX"];
786 Double_t fFeBlockY = conf_floats["Floor/MFeBlockY"];
787 Double_t fFeBlockZ = conf_floats["Floor/MFeBlockZ"];
788
789 // ************************ ACRYLIC
790 // Shapes definition
791 TGeoBBox *FeBlock_cb = new TGeoBBox("FeBlock_cb", fFeBlockX / 2, fFeBlockY / 2, fFeBlockZ / 2);
792 TGeoBBox *CBFrontWall_a =
793 new TGeoBBox("CBFrontWall_a", fCBFrontWallXDim / 2., fCBFrontWallYDim / 2., fAcrylicWidth / 2.);
794 TGeoBBox *CBLateral_a =
795 new TGeoBBox("CBLateral_a", fAcrylicWidth / 2., fCBFrontWallYDim / 2., (fCBLatWallZDim - 2 * fAcrylicWidth) / 2.);
796 TGeoBBox *CBExtraFront_a =
797 new TGeoBBox("CBExtraFront_a", fCBExtraXDim / 2., fCBFrontWallYDim / 2., fAcrylicWidth / 2.);
798 TGeoBBox *CBExtraLat_a =
799 new TGeoBBox("CBExtraLat_a", fAcrylicWidth / 2., fCBFrontWallYDim / 2., (fCBExtraZDim - 2 * fAcrylicWidth) / 2.);
800 TGeoBBox *CBTiny1_a = new TGeoBBox("CBTiny1_a", fAcrylicWidth / 2., fCBFrontWallYDim / 2.,
801 (fCBTinyZDim - fAcrylicWidth - fBPolyWidth) / 2.);
802 TGeoBBox *CBTiny2_a =
803 new TGeoBBox("CBTiny2_a", fAcrylicWidth / 2., fCBFrontWallYDim / 2., (fCBTinyZDim - fAcrylicWidth) / 2.);
804 TGeoBBox *CBRearWall_a =
805 new TGeoBBox("CBRearWall_a", fCBRearWallXDim / 2., fCBFrontWallYDim / 2., fAcrylicWidth / 2.);
806
807 TGeoTranslation *CBWallpos = new TGeoTranslation("CBWallpos", (fCBRearWallXDim - fFeBlockX) / 2. - 28.5,
808 (fCBFrontWallYDim - fFeBlockY) / 2. - 118., 0);
809 CBWallpos->RegisterYourself();
810
811 TGeoCompositeShape *CBWallDownstream =
812 new TGeoCompositeShape("CBWallDownstream", "CBRearWall_a-(FeBlock_cb:CBWallpos)");
813 TGeoPara *CBWallSlope_a =
814 new TGeoPara("CBWallSlope_a", fAcrylicWidth / 2., fCBFrontWallYDim / 2., fSlopedWallZProj / 2., 0, -15, 0);
815
816 // Coverings
817 Double_t LatCover1_averts[8][2];
818 LatCover1_averts[0][0] = 0.;
819 LatCover1_averts[0][1] = 0.;
820 LatCover1_averts[1][0] = fAcrylicWidth;
821 LatCover1_averts[1][1] = 0;
822 LatCover1_averts[2][0] = fAcrylicWidth;
823 LatCover1_averts[2][1] = -(176.71 - 170.72);
824 LatCover1_averts[3][0] = 0.;
825 LatCover1_averts[3][1] = -(176.71 - 170.72);
826
827 LatCover1_averts[4][0] = 0;
828 LatCover1_averts[4][1] = 0;
829 LatCover1_averts[5][0] = fAcrylicWidth;
830 LatCover1_averts[5][1] = 0;
831 LatCover1_averts[6][0] = fAcrylicWidth;
832 LatCover1_averts[6][1] = -(187.03 - 170.72);
833 LatCover1_averts[7][0] = 0;
834 LatCover1_averts[7][1] = -(187.03 - 170.72);
835 TGeoArb8 *LatCov1_a = new TGeoArb8("LatCov1_a", 176. / 2., (Double_t *)LatCover1_averts);
836
837 Double_t tanalpha = (183.26 - 170.72) / 144.;
838 Double_t LatCover21_averts[8][2];
839 LatCover21_averts[0][0] = 0.;
840 LatCover21_averts[0][1] = 0.;
841 LatCover21_averts[1][0] = fAcrylicWidth;
842 LatCover21_averts[1][1] = 0;
843 LatCover21_averts[2][0] = fAcrylicWidth;
844 LatCover21_averts[2][1] = 0;
845 LatCover21_averts[3][0] = 0.;
846 LatCover21_averts[3][1] = 0;
847
848 LatCover21_averts[4][0] = 0;
849 LatCover21_averts[4][1] = 0;
850 LatCover21_averts[5][0] = fAcrylicWidth;
851 LatCover21_averts[5][1] = 0;
852 LatCover21_averts[6][0] = fAcrylicWidth;
853 LatCover21_averts[6][1] = -fCBTinyZDim * tanalpha;
854 LatCover21_averts[7][0] = 0;
855 LatCover21_averts[7][1] = -fCBTinyZDim * tanalpha;
856 TGeoArb8 *LatCov21_a = new TGeoArb8("LatCov21_a", fCBTinyZDim / 2., (Double_t *)LatCover21_averts);
857
858 Double_t LatCover22_averts[8][2];
859 LatCover22_averts[0][0] = 0.;
860 LatCover22_averts[0][1] = 0.;
861 LatCover22_averts[1][0] = fAcrylicWidth;
862 LatCover22_averts[1][1] = 0;
863 LatCover22_averts[2][0] = fAcrylicWidth;
864 LatCover22_averts[2][1] = -fCBTinyZDim * tanalpha;
865 LatCover22_averts[3][0] = 0.;
866 LatCover22_averts[3][1] = -fCBTinyZDim * tanalpha;
867
868 LatCover22_averts[4][0] = -fSlopedWallZProj * TMath::Tan(TMath::DegToRad() * 15.);
869 LatCover22_averts[4][1] = 0;
870 LatCover22_averts[5][0] = -fSlopedWallZProj * TMath::Tan(TMath::DegToRad() * 15.) + fAcrylicWidth;
871 LatCover22_averts[5][1] = 0;
872 LatCover22_averts[6][0] = -fSlopedWallZProj * TMath::Tan(TMath::DegToRad() * 15.) + fAcrylicWidth;
873 LatCover22_averts[6][1] = -(fCBTinyZDim + fSlopedWallZProj) * tanalpha;
874 LatCover22_averts[7][0] = -fSlopedWallZProj * TMath::Tan(TMath::DegToRad() * 15.);
875 LatCover22_averts[7][1] = -(fCBTinyZDim + fSlopedWallZProj) * tanalpha;
876 TGeoArb8 *LatCov22_a = new TGeoArb8("LatCov22_a", fSlopedWallZProj / 2., (Double_t *)LatCover22_averts);
877
878 Double_t LatCover23_averts[8][2];
879 LatCover23_averts[0][0] = 0.;
880 LatCover23_averts[0][1] = 0.;
881 LatCover23_averts[1][0] = fAcrylicWidth;
882 LatCover23_averts[1][1] = 0;
883 LatCover23_averts[2][0] = fAcrylicWidth;
884 LatCover23_averts[2][1] = -(fCBTinyZDim + fSlopedWallZProj) * tanalpha;
885 LatCover23_averts[3][0] = 0.;
886 LatCover23_averts[3][1] = -(fCBTinyZDim + fSlopedWallZProj) * tanalpha;
887
888 LatCover23_averts[4][0] = 0;
889 LatCover23_averts[4][1] = 0;
890 LatCover23_averts[5][0] = fAcrylicWidth;
891 LatCover23_averts[5][1] = 0;
892 LatCover23_averts[6][0] = fAcrylicWidth;
893 LatCover23_averts[6][1] = -(183.26 - 170.72);
894 LatCover23_averts[7][0] = 0;
895 LatCover23_averts[7][1] = -(183.26 - 170.72);
896 TGeoArb8 *LatCov23_a = new TGeoArb8("LatCov23_a", (fCBTinyZDim - fBPolyWidth) / 2., (Double_t *)LatCover23_averts);
897
898 Double_t tanbeta = (176.74 - 170.72) / 219.0;
899 Double_t FrontCover1_averts[8][2];
900 FrontCover1_averts[0][0] = (fCBFrontWallXDim-2*fAcrylicWidth)/2.;
901 FrontCover1_averts[0][1] = 0.;
902 FrontCover1_averts[1][0] = (fCBFrontWallXDim-2*fAcrylicWidth)/2.;
903 FrontCover1_averts[1][1] = -(fAcrylicWidth)*tanbeta;
904 FrontCover1_averts[2][0] = -(fCBFrontWallXDim-2*fAcrylicWidth)/2.;
905 FrontCover1_averts[2][1] = -(fCBFrontWallXDim-fAcrylicWidth)*tanbeta;
906 FrontCover1_averts[3][0] = -(fCBFrontWallXDim-2*fAcrylicWidth)/2.;
907 FrontCover1_averts[3][1] = 0;
908
909 FrontCover1_averts[4][0] = (fCBFrontWallXDim-2*fAcrylicWidth)/2.;
910 FrontCover1_averts[4][1] = 0;
911 FrontCover1_averts[5][0] = (fCBFrontWallXDim-2*fAcrylicWidth)/2.;
912 FrontCover1_averts[5][1] = -(fAcrylicWidth)*tanbeta;
913 FrontCover1_averts[6][0] = -(fCBFrontWallXDim-2*fAcrylicWidth)/2.;
914 FrontCover1_averts[6][1] = -(fCBFrontWallXDim-fAcrylicWidth)*tanbeta;
915 FrontCover1_averts[7][0] = -(fCBFrontWallXDim-2*fAcrylicWidth)/2.;
916 FrontCover1_averts[7][1] = 0;
917 TGeoArb8 *FrontCov1_a = new TGeoArb8("FrontCov1_a", fAcrylicWidth / 2., (Double_t *)FrontCover1_averts);
918
919 Double_t tangamma = (187.03 - 183.26) / (fCBRearWallXDim - fAcrylicWidth + fCBExtraXDim);
920 Double_t xdim_projected = (187.03 - 170.72) / tangamma;
921 Double_t extra = xdim_projected - (fCBRearWallXDim - fAcrylicWidth + fCBExtraXDim);
922 Double_t RearCover_11_averts[8][2];
923 RearCover_11_averts[0][0] = (fCBRearWallXDim-2*fAcrylicWidth)/2.;
924 RearCover_11_averts[0][1] = 0.;
925 RearCover_11_averts[1][0] = (fCBRearWallXDim-2*fAcrylicWidth)/2.;
926 RearCover_11_averts[1][1] = -(183.26-170.72); // approx
927 RearCover_11_averts[2][0] = -(fCBRearWallXDim-2*fAcrylicWidth)/2.;
928 RearCover_11_averts[2][1] = -(fCBRearWallXDim-fAcrylicWidth+extra)*tangamma;
929 RearCover_11_averts[3][0] = -(fCBRearWallXDim-2*fAcrylicWidth)/2.;
930 RearCover_11_averts[3][1] = 0;
931
932 RearCover_11_averts[4][0] = (fCBRearWallXDim-2*fAcrylicWidth)/2.;
933 RearCover_11_averts[4][1] = 0;
934 RearCover_11_averts[5][0] = (fCBRearWallXDim-2*fAcrylicWidth)/2.;
935 RearCover_11_averts[5][1] = -(183.26-170.72); // approx
936 RearCover_11_averts[6][0] = -(fCBRearWallXDim-2*fAcrylicWidth)/2.;
937 RearCover_11_averts[6][1] = -(fCBRearWallXDim-fAcrylicWidth+extra)*tangamma;
938 RearCover_11_averts[7][0] = -(fCBRearWallXDim-2*fAcrylicWidth)/2.;
939 RearCover_11_averts[7][1] = 0;
940 TGeoArb8 *RearCov11_atot = new TGeoArb8("RearCov11_atot", fAcrylicWidth / 2., (Double_t *)RearCover_11_averts);
941 TGeoTranslation *CBWallpos2 =
942 new TGeoTranslation("CBWallpos2", (fCBRearWallXDim - fFeBlockX) / 2. - 28.5, 0, -fAcrylicWidth / 2.);
943 CBWallpos2->RegisterYourself();
944 TGeoCompositeShape *RearCov11_a = new TGeoCompositeShape("RearCov11_a", "RearCov11_atot-(FeBlock_cb:CBWallpos2)");
945
946 Double_t RearCover_12_averts[8][2];
947 RearCover_12_averts[0][0] = fAcrylicWidth/2.;
948 RearCover_12_averts[0][1] = 0.;
949 RearCover_12_averts[1][0] = fAcrylicWidth/2.;
950 RearCover_12_averts[1][1] = -(fCBRearWallXDim-fAcrylicWidth+extra)*tangamma;
951 RearCover_12_averts[2][0] = -fAcrylicWidth/2.;
952 RearCover_12_averts[2][1] = -(fCBRearWallXDim+extra)*tangamma;
953 RearCover_12_averts[3][0] = -fAcrylicWidth/2.;
954 RearCover_12_averts[3][1] = 0;
955
956 RearCover_12_averts[4][0] = fAcrylicWidth/2.;
957 RearCover_12_averts[4][1] = 0;
958 RearCover_12_averts[5][0] = fAcrylicWidth/2.;
959 RearCover_12_averts[5][1] = -(fCBRearWallXDim-fAcrylicWidth+extra)*tangamma;
960 RearCover_12_averts[6][0] = -fAcrylicWidth/2.;
961 RearCover_12_averts[6][1] = -(fCBRearWallXDim+extra)*tangamma;
962 RearCover_12_averts[7][0] = -fAcrylicWidth/2.;
963 RearCover_12_averts[7][1] = 0;
964 TGeoArb8 *RearCov12_a = new TGeoArb8("RearCov12_a", fCBExtraZDim / 2., (Double_t *)RearCover_12_averts);
965
966 Double_t RearCover_13_averts[8][2];
967 RearCover_13_averts[0][0] = (fCBExtraXDim-2*fAcrylicWidth)/2.;
968 RearCover_13_averts[0][1] = 0.;
969 RearCover_13_averts[1][0] = (fCBExtraXDim-2*fAcrylicWidth)/2.;
970 RearCover_13_averts[1][1] = -(fCBRearWallXDim+extra)*tangamma;
971 RearCover_13_averts[2][0] = -(fCBExtraXDim-2*fAcrylicWidth)/2.;
972 RearCover_13_averts[2][1] = -(187.03-170.72);
973 RearCover_13_averts[3][0] = -(fCBExtraXDim-2*fAcrylicWidth)/2.;
974 RearCover_13_averts[3][1] = 0;
975
976 RearCover_13_averts[4][0] = (fCBExtraXDim-2*fAcrylicWidth)/2.;
977 RearCover_13_averts[4][1] = 0;
978 RearCover_13_averts[5][0] = (fCBExtraXDim-2*fAcrylicWidth)/2.;
979 RearCover_13_averts[5][1] = -(fCBRearWallXDim+extra)*tangamma;
980 RearCover_13_averts[6][0] = -(fCBExtraXDim-2*fAcrylicWidth)/2.;
981 RearCover_13_averts[6][1] = -(187.03-170.72);
982 RearCover_13_averts[7][0] = -(fCBExtraXDim-2*fAcrylicWidth)/2.;
983 RearCover_13_averts[7][1] = 0;
984 TGeoArb8 *RearCov13_a = new TGeoArb8("RearCov13_a", fAcrylicWidth / 2., (Double_t *)RearCover_13_averts);
985
986 // Acrylic mother shape definition
987 TGeoTranslation *FrontWallpos =
988 new TGeoTranslation("FrontWallpos", -fCBRearWallXDim / 2. - fCBExtraXDim + fAcrylicWidth + fCBFrontWallXDim / 2.,
989 0, -(fSlopedWallZProj + 2 * (fCBTinyZDim - fAcrylicWidth) + fAcrylicWidth) + fBPolyWidth);
990 FrontWallpos->RegisterYourself();
991 TGeoTranslation *Tiny1pos =
992 new TGeoTranslation("Tiny1pos", (fCBRearWallXDim - fAcrylicWidth) / 2., 0, -fCBTinyZDim / 2. + fBPolyWidth / 2.);
993 Tiny1pos->RegisterYourself();
994 TGeoTranslation *SlopeWallpos =
995 new TGeoTranslation("SlopeWallpos",
996 (fCBRearWallXDim + fAcrylicWidth) / 2. + fAcrylicWidth +
997 fSlopedWallZProj / (2 * TMath::Tan(TMath::DegToRad() * 85.)),
998 0, -fCBTinyZDim - fSlopedWallZProj / 2. + fAcrylicWidth / 2. + fBPolyWidth);
999 SlopeWallpos->RegisterYourself();
1000 TGeoTranslation *Tiny2pos = new TGeoTranslation(
1001 "Tiny2pos",
1002 3 * fAcrylicWidth + fCBRearWallXDim / 2. + fAcrylicWidth / 2. +
1003 fSlopedWallZProj / (TMath::Tan(TMath::DegToRad() * 85.)),
1004 0, -(fSlopedWallZProj + 2 * (fCBTinyZDim - fAcrylicWidth) + fAcrylicWidth) + fCBTinyZDim / 2. + fBPolyWidth);
1005 Tiny2pos->RegisterYourself();
1006 TGeoTranslation *CBExtraLatpos = new TGeoTranslation("CBExtraLatpos", -fCBRearWallXDim / 2. + fAcrylicWidth / 2., 0,
1007 fAcrylicWidth / 2. + (fCBExtraZDim - 2 * fAcrylicWidth) / 2.);
1008 CBExtraLatpos->RegisterYourself();
1009 TGeoTranslation *CBExtraFrontpos = new TGeoTranslation(
1010 "CBExtraFrontpos", -fCBRearWallXDim / 2. + fAcrylicWidth - fCBExtraXDim / 2., 0, fCBExtraZDim - fAcrylicWidth);
1011 CBExtraFrontpos->RegisterYourself();
1012 TGeoTranslation *CBLateralpos =
1013 new TGeoTranslation("CBLateralpos", -fCBRearWallXDim / 2. - fCBExtraXDim + fAcrylicWidth + fAcrylicWidth / 2., 0,
1014 fCBExtraZDim - fCBLatWallZDim / 2. - fAcrylicWidth / 2.);
1015 CBLateralpos->RegisterYourself();
1016 TGeoTranslation *CBLatCov1pos =
1017 new TGeoTranslation("CBLatCov1pos", -fCBRearWallXDim / 2. - fCBExtraXDim + fAcrylicWidth, -fCBFrontWallYDim / 2.,
1018 fCBExtraZDim - fCBLatWallZDim / 2. - fAcrylicWidth / 2.);
1019 CBLatCov1pos->RegisterYourself();
1020 TGeoTranslation *CBLatCov21pos = new TGeoTranslation(
1021 "CBLatCov21pos",
1022 3 * fAcrylicWidth + fCBRearWallXDim / 2. + fSlopedWallZProj / (TMath::Tan(TMath::DegToRad() * 85.)),
1023 -fCBFrontWallYDim / 2.,
1024 -(fSlopedWallZProj + 2 * (fCBTinyZDim - fAcrylicWidth) + fAcrylicWidth) + fCBTinyZDim / 2. + fBPolyWidth -
1025 fAcrylicWidth / 2.);
1026 CBLatCov21pos->RegisterYourself();
1027 TGeoTranslation *CBLatCov22pos = new TGeoTranslation(
1028 "CBLatCov22pos",
1029 3 * fAcrylicWidth + fCBRearWallXDim / 2. + fSlopedWallZProj / (TMath::Tan(TMath::DegToRad() * 85.)),
1030 -fCBFrontWallYDim / 2.,
1031 -(fSlopedWallZProj / 2 + 2 * (fCBTinyZDim - fAcrylicWidth) + fAcrylicWidth) + fCBTinyZDim / 2. + fBPolyWidth -
1032 fAcrylicWidth / 2. + fCBTinyZDim / 2.);
1033 CBLatCov22pos->RegisterYourself();
1034 TGeoTranslation *CBLatCov23pos =
1035 new TGeoTranslation("CBLatCov23pos", (fCBRearWallXDim - fAcrylicWidth) / 2. - fAcrylicWidth / 2.,
1036 -fCBFrontWallYDim / 2., -fCBTinyZDim / 2. + fBPolyWidth + 0.5);
1037 CBLatCov23pos->RegisterYourself();
1038 TGeoTranslation *CBFrontCov1pos = new TGeoTranslation(
1039 "CBFrontCov1pos", -fCBRearWallXDim / 2. - fCBExtraXDim + fAcrylicWidth + fCBFrontWallXDim / 2.,
1040 -fCBFrontWallYDim / 2., -(fSlopedWallZProj + 2 * (fCBTinyZDim - fAcrylicWidth) + fAcrylicWidth) + fBPolyWidth);
1041 CBFrontCov1pos->RegisterYourself();
1042 TGeoTranslation *CBRearCov11pos = new TGeoTranslation("CBRearCov11pos", 0, -fCBFrontWallYDim / 2., 0);
1043 CBRearCov11pos->RegisterYourself();
1044 TGeoTranslation *CBRearCov12pos =
1045 new TGeoTranslation("CBRearCov12pos", -fCBRearWallXDim / 2. + fAcrylicWidth / 2., -fCBFrontWallYDim / 2.,
1046 fAcrylicWidth / 2. + (fCBExtraZDim - 2 * fAcrylicWidth) / 2.);
1047 CBRearCov12pos->RegisterYourself();
1048 TGeoTranslation *CBRearCov13pos =
1049 new TGeoTranslation("CBRearCov13pos", -fCBRearWallXDim / 2. + fAcrylicWidth - fCBExtraXDim / 2.,
1050 -fCBFrontWallYDim / 2., fCBExtraZDim - fAcrylicWidth);
1051 CBRearCov13pos->RegisterYourself();
1052
1053 // Acrylic mother volume definition
1054 TGeoCompositeShape *COLDBOXA = new TGeoCompositeShape(
1055 "COLDBOXA",
1056 "CBWallDownstream+(CBFrontWall_a:FrontWallpos)+(CBTiny1_a:Tiny1pos)+(CBWallSlope_a:SlopeWallpos)+(CBTiny2_a:"
1057 "Tiny2pos)+(CBExtraLat_a:CBExtraLatpos)+(CBExtraFront_a:CBExtraFrontpos)+(CBLateral_a:CBLateralpos)+(LatCov1_a:"
1058 "CBLatCov1pos)+(LatCov21_a:CBLatCov21pos)+(LatCov22_a:CBLatCov22pos)+(LatCov23_a:CBLatCov23pos)+(FrontCov1_a:"
1059 "CBFrontCov1pos)+(RearCov11_a:CBRearCov11pos)+(RearCov12_a:CBRearCov12pos)+(RearCov13_a:CBRearCov13pos)");
1060 TGeoVolume *volCOLDBOXA = new TGeoVolume("volCOLDBOXA", COLDBOXA, Acrylic);
1061
1062 // ************************ BORATED POLYETHYLENE
1063 Double_t fCBFrontWallXDim_b = fCBFrontWallXDim - 2 * fAcrylicWidth - fBPolyWidth; // cm
1064 Double_t fCBFrontWallYDim_b = fCBFrontWallYDim - fBPolyWidth;
1065 Double_t fCBLatWallZDim_b = fCBLatWallZDim - 2 * fAcrylicWidth; // cm
1066 Double_t fCBExtraXDim_b = fCBExtraXDim - 2 * fAcrylicWidth; // cm
1067 Double_t fCBRearWallXDim_b = fCBRearWallXDim - fAcrylicWidth;
1068 // Shapes definition
1069 TGeoBBox *CBFrontWall_b = new TGeoBBox(
1070 "CBFrontWall_b", fCBFrontWallXDim_b / 2. + (fAcrylicWidth - fBPolyWidth) / 10., fCBFrontWallYDim_b / 2.,
1071 fBPolyWidth / 2.); // (fAcrylicWidth-fBPolyWidth)/10. is due to approximations, I guess
1072 TGeoBBox *CBLateral_b = new TGeoBBox("CBLateral_b", fBPolyWidth / 2., fCBFrontWallYDim_b / 2.,
1073 (fCBLatWallZDim - 2 * fBPolyWidth) / 2. - fAcrylicWidth);
1074 TGeoBBox *CBExtraFront_b =
1075 new TGeoBBox("CBExtraFront_b", fCBExtraXDim_b / 2., fCBFrontWallYDim_b / 2., fBPolyWidth / 2.);
1076 TGeoBBox *CBExtraLat_b = new TGeoBBox("CBExtraLat_b", fBPolyWidth / 2., fCBFrontWallYDim_b / 2.,
1077 (fCBExtraZDim - fAcrylicWidth - fBPolyWidth) / 2.);
1078 TGeoBBox *CBTiny1_b = new TGeoBBox("CBTiny1_b", fBPolyWidth / 2., fCBFrontWallYDim_b / 2.,
1079 (fCBTinyZDim - fAcrylicWidth - fBPolyWidth) / 2.);
1080 TGeoBBox *CBTiny2_b =
1081 new TGeoBBox("CBTiny2_b", fBPolyWidth / 2., fCBFrontWallYDim_b / 2., (fCBTinyZDim - fAcrylicWidth) / 2.);
1082 TGeoBBox *CBRearWall_b =
1083 new TGeoBBox("CBRearWall_b", fCBRearWallXDim_b / 2., fCBFrontWallYDim_b / 2., fBPolyWidth / 2.);
1084 TGeoPara *CBWallSlope_b =
1085 new TGeoPara("CBWallSlope_b", fBPolyWidth / 2., fCBFrontWallYDim_b / 2., fSlopedWallZProj / 2., 0, -15, 0);
1086
1087 // Coverings
1088 Double_t tanomega = (187.03 - 176.71) / 176.;
1089 Double_t LatCover1_bverts[8][2];
1090 LatCover1_bverts[0][0] = 0.;
1091 LatCover1_bverts[0][1] = 0.;
1092 LatCover1_bverts[1][0] = fBPolyWidth;
1093 LatCover1_bverts[1][1] = 0;
1094 LatCover1_bverts[2][0] = fBPolyWidth;
1095 LatCover1_bverts[2][1] = -(176.71 - 170.72) - fAcrylicWidth * tanomega;
1096 LatCover1_bverts[3][0] = 0.;
1097 LatCover1_bverts[3][1] = -(176.71 - 170.72) - fAcrylicWidth * tanomega;
1098
1099 LatCover1_bverts[4][0] = 0;
1100 LatCover1_bverts[4][1] = 0;
1101 LatCover1_bverts[5][0] = fBPolyWidth;
1102 LatCover1_bverts[5][1] = 0;
1103 LatCover1_bverts[6][0] = fBPolyWidth;
1104 LatCover1_bverts[6][1] = -(187.03 - 170.72) + fAcrylicWidth * tanomega;
1105 LatCover1_bverts[7][0] = 0;
1106 LatCover1_bverts[7][1] = -(187.03 - 170.72) + fAcrylicWidth * tanomega;
1107 TGeoArb8 *LatCov1_b = new TGeoArb8("LatCov1_b", 176. / 2. - fAcrylicWidth, (Double_t *)LatCover1_bverts);
1108
1109 Double_t LatCover21_bverts[8][2];
1110 LatCover21_bverts[0][0] = 0.;
1111 LatCover21_bverts[0][1] = 0.;
1112 LatCover21_bverts[1][0] = fBPolyWidth;
1113 LatCover21_bverts[1][1] = 0;
1114 LatCover21_bverts[2][0] = fBPolyWidth;
1115 LatCover21_bverts[2][1] = 0.; //-fAcrylicWidth*tanalpha;
1116 LatCover21_bverts[3][0] = 0.;
1117 LatCover21_bverts[3][1] = 0.; //-fAcrylicWidth*tanalpha;
1118
1119 LatCover21_bverts[4][0] = 0;
1120 LatCover21_bverts[4][1] = 0;
1121 LatCover21_bverts[5][0] = fBPolyWidth;
1122 LatCover21_bverts[5][1] = 0;
1123 LatCover21_bverts[6][0] = fBPolyWidth;
1124 LatCover21_bverts[6][1] = -fCBTinyZDim * tanalpha;
1125 LatCover21_bverts[7][0] = 0;
1126 LatCover21_bverts[7][1] = -fCBTinyZDim * tanalpha;
1127 TGeoArb8 *LatCov21_b = new TGeoArb8("LatCov21_b", (fCBTinyZDim - fAcrylicWidth) / 2., (Double_t *)LatCover21_bverts);
1128
1129 Double_t LatCover22_bverts[8][2];
1130 LatCover22_bverts[0][0] = 0.;
1131 LatCover22_bverts[0][1] = 0.;
1132 LatCover22_bverts[1][0] = fBPolyWidth;
1133 LatCover22_bverts[1][1] = 0;
1134 LatCover22_bverts[2][0] = fBPolyWidth;
1135 LatCover22_bverts[2][1] = -fCBTinyZDim * tanalpha;
1136 LatCover22_bverts[3][0] = 0.;
1137 LatCover22_bverts[3][1] = -fCBTinyZDim * tanalpha;
1138
1139 LatCover22_bverts[4][0] = -fSlopedWallZProj * TMath::Tan(TMath::DegToRad() * 15.);
1140 LatCover22_bverts[4][1] = 0;
1141 LatCover22_bverts[5][0] = -fSlopedWallZProj * TMath::Tan(TMath::DegToRad() * 15.) + fBPolyWidth;
1142 LatCover22_bverts[5][1] = 0;
1143 LatCover22_bverts[6][0] = -fSlopedWallZProj * TMath::Tan(TMath::DegToRad() * 15.) + fBPolyWidth;
1144 LatCover22_bverts[6][1] = -(fCBTinyZDim + fSlopedWallZProj) * tanalpha;
1145 LatCover22_bverts[7][0] = -fSlopedWallZProj * TMath::Tan(TMath::DegToRad() * 15.);
1146 LatCover22_bverts[7][1] = -(fCBTinyZDim + fSlopedWallZProj) * tanalpha;
1147 TGeoArb8 *LatCov22_b = new TGeoArb8("LatCov22_b", fSlopedWallZProj / 2., (Double_t *)LatCover22_bverts);
1148
1149 Double_t LatCover23_bverts[8][2];
1150 LatCover23_bverts[0][0] = 0.;
1151 LatCover23_bverts[0][1] = 0.;
1152 LatCover23_bverts[1][0] = fBPolyWidth;
1153 LatCover23_bverts[1][1] = 0;
1154 LatCover23_bverts[2][0] = fBPolyWidth;
1155 LatCover23_bverts[2][1] = -(fCBTinyZDim + fSlopedWallZProj) * tanalpha;
1156 LatCover23_bverts[3][0] = 0.;
1157 LatCover23_bverts[3][1] = -(fCBTinyZDim + fSlopedWallZProj) * tanalpha;
1158
1159 LatCover23_bverts[4][0] = 0;
1160 LatCover23_bverts[4][1] = 0;
1161 LatCover23_bverts[5][0] = fBPolyWidth;
1162 LatCover23_bverts[5][1] = 0;
1163 LatCover23_bverts[6][0] = fBPolyWidth;
1164 LatCover23_bverts[6][1] = -(183.26 - 170.72) + fAcrylicWidth * tanalpha;
1165 LatCover23_bverts[7][0] = 0;
1166 LatCover23_bverts[7][1] = -(183.26 - 170.72) + fAcrylicWidth * tanalpha;
1167 TGeoArb8 *LatCov23_b =
1168 new TGeoArb8("LatCov23_b", (fCBTinyZDim - fBPolyWidth - fAcrylicWidth) / 2., (Double_t *)LatCover23_bverts);
1169
1170 Double_t FrontCover1_bverts[8][2];
1171 FrontCover1_bverts[0][0] = (fCBFrontWallXDim_b-2*fBPolyWidth)/2.;
1172 FrontCover1_bverts[0][1] = 0.;
1173 FrontCover1_bverts[1][0] = (fCBFrontWallXDim_b-2*fBPolyWidth)/2.;
1174 FrontCover1_bverts[1][1] = -(fBPolyWidth)*tanbeta;
1175 FrontCover1_bverts[2][0] = -(fCBFrontWallXDim_b-2*fBPolyWidth)/2.;
1176 FrontCover1_bverts[2][1] = -(176.71-170.72)-fAcrylicWidth*tanomega;
1177 FrontCover1_bverts[3][0] = -(fCBFrontWallXDim_b-2*fBPolyWidth)/2.;
1178 FrontCover1_bverts[3][1] = 0;
1179
1180 FrontCover1_bverts[4][0] = (fCBFrontWallXDim_b-2*fBPolyWidth)/2.;
1181 FrontCover1_bverts[4][1] = 0;
1182 FrontCover1_bverts[5][0] = (fCBFrontWallXDim_b-2*fBPolyWidth)/2.;
1183 FrontCover1_bverts[5][1] = -(fBPolyWidth)*tanbeta;
1184 FrontCover1_bverts[6][0] = -(fCBFrontWallXDim_b-2*fBPolyWidth)/2.;
1185 FrontCover1_bverts[6][1] = -(176.71-170.72)-fAcrylicWidth*tanomega;
1186 FrontCover1_bverts[7][0] = -(fCBFrontWallXDim_b-2*fBPolyWidth)/2.;
1187 FrontCover1_bverts[7][1] = 0;
1188 TGeoArb8 *FrontCov1_b = new TGeoArb8("FrontCov1_b", fBPolyWidth / 2., (Double_t *)FrontCover1_bverts);
1189
1190 Double_t RearCover_11_bverts[8][2];
1191 RearCover_11_bverts[0][0] = (fCBRearWallXDim_b-2*fBPolyWidth)/2.+fBPolyWidth/2.;
1192 RearCover_11_bverts[0][1] = 0.;
1193 RearCover_11_bverts[1][0] = (fCBRearWallXDim_b-2*fBPolyWidth)/2.+fBPolyWidth/2.;
1194 RearCover_11_bverts[1][1] = -(183.26-170.72)+fAcrylicWidth*tanalpha; // approx
1195 RearCover_11_bverts[2][0] = -(fCBRearWallXDim_b-2*fBPolyWidth)/2.-fBPolyWidth/2.;
1196 RearCover_11_bverts[2][1] = -(fCBRearWallXDim_b-fBPolyWidth+extra)*tangamma;
1197 RearCover_11_bverts[3][0] = -(fCBRearWallXDim_b-2*fBPolyWidth)/2.-fBPolyWidth/2.;
1198 RearCover_11_bverts[3][1] = 0;
1199
1200 RearCover_11_bverts[4][0] = (fCBRearWallXDim_b-2*fBPolyWidth)/2.+fBPolyWidth/2.;
1201 RearCover_11_bverts[4][1] = 0;
1202 RearCover_11_bverts[5][0] = (fCBRearWallXDim_b-2*fBPolyWidth)/2.+fBPolyWidth/2.;
1203 RearCover_11_bverts[5][1] = -(183.26-170.72)+fAcrylicWidth*tanalpha; // approx
1204 RearCover_11_bverts[6][0] = -(fCBRearWallXDim_b-2*fBPolyWidth)/2.-fBPolyWidth/2.;
1205 RearCover_11_bverts[6][1] = -(fCBRearWallXDim_b-fBPolyWidth+extra)*tangamma;
1206 RearCover_11_bverts[7][0] = -(fCBRearWallXDim_b-2*fBPolyWidth)/2.-fBPolyWidth/2.;
1207 RearCover_11_bverts[7][1] = 0;
1208 TGeoArb8 *RearCov11_b = new TGeoArb8("RearCov11_b", fBPolyWidth / 2., (Double_t *)RearCover_11_bverts);
1209
1210 Double_t RearCover_12_bverts[8][2];
1211 RearCover_12_bverts[0][0] = fBPolyWidth/2.;
1212 RearCover_12_bverts[0][1] = 0.;
1213 RearCover_12_bverts[1][0] = fBPolyWidth/2.;
1214 RearCover_12_bverts[1][1] = -(fCBRearWallXDim_b-fBPolyWidth+extra)*tangamma;
1215 RearCover_12_bverts[2][0] = -fBPolyWidth/2.;
1216 RearCover_12_bverts[2][1] = -(fCBRearWallXDim_b+extra)*tangamma;
1217 RearCover_12_bverts[3][0] = -fBPolyWidth/2.;
1218 RearCover_12_bverts[3][1] = 0;
1219
1220 RearCover_12_bverts[4][0] = fBPolyWidth/2.;
1221 RearCover_12_bverts[4][1] = 0;
1222 RearCover_12_bverts[5][0] = fBPolyWidth/2.;
1223 RearCover_12_bverts[5][1] = -(fCBRearWallXDim_b-fAcrylicWidth+extra)*tangamma;
1224 RearCover_12_bverts[6][0] = -fBPolyWidth/2.;
1225 RearCover_12_bverts[6][1] = -(fCBRearWallXDim_b+extra)*tangamma;
1226 RearCover_12_bverts[7][0] = -fBPolyWidth/2.;
1227 RearCover_12_bverts[7][1] = 0;
1228 TGeoArb8 *RearCov12_b =
1229 new TGeoArb8("RearCov12_b", (fCBExtraZDim - fAcrylicWidth + fBPolyWidth) / 2., (Double_t *)RearCover_12_bverts);
1230
1231 Double_t RearCover_13_bverts[8][2];
1232 RearCover_13_bverts[0][0] = (fCBExtraXDim_b-2*fBPolyWidth)/2.;
1233 RearCover_13_bverts[0][1] = 0.;
1234 RearCover_13_bverts[1][0] = (fCBExtraXDim_b-2*fBPolyWidth)/2.;
1235 RearCover_13_bverts[1][1] = -(fCBRearWallXDim_b+extra)*tangamma;
1236 RearCover_13_bverts[2][0] = -(fCBExtraXDim_b-2*fBPolyWidth)/2.;
1237 RearCover_13_bverts[2][1] = -(187.03-170.72)+fAcrylicWidth*tanomega;
1238 RearCover_13_bverts[3][0] = -(fCBExtraXDim_b-2*fBPolyWidth)/2.;
1239 RearCover_13_bverts[3][1] = 0;
1240
1241 RearCover_13_bverts[4][0] = (fCBExtraXDim_b-2*fBPolyWidth)/2.;
1242 RearCover_13_bverts[4][1] = 0;
1243 RearCover_13_bverts[5][0] = (fCBExtraXDim_b-2*fBPolyWidth)/2.;
1244 RearCover_13_bverts[5][1] = -(fCBRearWallXDim_b+extra)*tangamma;
1245 RearCover_13_bverts[6][0] = -(fCBExtraXDim_b-2*fBPolyWidth)/2.;
1246 RearCover_13_bverts[6][1] = -(187.03-170.72)+fAcrylicWidth*tanomega;
1247 RearCover_13_bverts[7][0] = -(fCBExtraXDim_b-2*fBPolyWidth)/2.;
1248 RearCover_13_bverts[7][1] = 0;
1249 TGeoArb8 *RearCov13_b = new TGeoArb8("RearCov13_b", fBPolyWidth / 2., (Double_t *)RearCover_13_bverts);
1250
1251 // Borated Polyethylene mother shape definition
1252 TGeoTranslation *FrontWallpos_b = new TGeoTranslation(
1253 "FrontWallpos_b", -fCBRearWallXDim_b / 2. - fCBExtraXDim_b + fBPolyWidth + fCBFrontWallXDim_b / 2. + 0.1, 0,
1254 -fSlopedWallZProj - (fCBTinyZDim - fAcrylicWidth + fBPolyWidth)); // +0.1 is due to approximations, I guess
1255 FrontWallpos_b->RegisterYourself();
1256 TGeoTranslation *Tiny1pos_b =
1257 new TGeoTranslation("Tiny1pos_b", fCBRearWallXDim_b / 2. + fBPolyWidth / 2., 0, -fBPolyWidth / 2.);
1258 Tiny1pos_b->RegisterYourself();
1259 TGeoTranslation *SlopeWallpos_b = new TGeoTranslation(
1260 "SlopeWallpos_b",
1261 fSlopedWallZProj / (2 * TMath::Tan(TMath::DegToRad() * 85.)) + fCBRearWallXDim_b / 2. + 3 * fBPolyWidth, 0,
1262 -fBPolyWidth / 2. - fSlopedWallZProj / 2. - (fCBTinyZDim - fAcrylicWidth - fBPolyWidth) / 2.);
1263 SlopeWallpos_b->RegisterYourself();
1264 TGeoTranslation *Tiny2pos_b = new TGeoTranslation("Tiny2pos_b",
1265 5 * fBPolyWidth + fCBRearWallXDim / 2. +
1266 fSlopedWallZProj / (TMath::Tan(TMath::DegToRad() * 85.)) -
1267 (fAcrylicWidth - fBPolyWidth) / 2.,
1268 0, -fSlopedWallZProj - (fCBTinyZDim - fAcrylicWidth));
1269 Tiny2pos_b->RegisterYourself();
1270 TGeoTranslation *CBExtraLatpos_b = new TGeoTranslation("CBExtraLatpos_b", -fCBRearWallXDim_b / 2. + fBPolyWidth / 2.,
1271 0, (fCBExtraZDim - fAcrylicWidth) / 2.);
1272 CBExtraLatpos_b->RegisterYourself();
1273 TGeoTranslation *CBExtraFrontpos_b =
1274 new TGeoTranslation("CBExtraFrontpos_b", -fCBRearWallXDim_b / 2. + fBPolyWidth - fCBExtraXDim_b / 2., 0,
1275 fCBExtraZDim - fAcrylicWidth);
1276 CBExtraFrontpos_b->RegisterYourself();
1277 TGeoTranslation *CBLateralpos_b =
1278 new TGeoTranslation("CBLateralpos_b", -fCBRearWallXDim_b / 2. - fCBExtraXDim_b + fBPolyWidth + fBPolyWidth / 2.,
1279 0, fCBExtraZDim - fCBLatWallZDim_b / 2. + fBPolyWidth / 2. - fAcrylicWidth);
1280 CBLateralpos_b->RegisterYourself();
1281 TGeoTranslation *CBLatCov1pos_b = new TGeoTranslation(
1282 "CBLatCov1pos_b", -fCBRearWallXDim / 2. - fCBExtraXDim + 2 * fAcrylicWidth + fBPolyWidth + fAcrylicWidth / 2.,
1283 -fCBFrontWallYDim_b / 2., fCBExtraZDim - fCBLatWallZDim / 2. + fAcrylicWidth / 2 - 0.5);
1284 CBLatCov1pos_b->RegisterYourself();
1285 TGeoTranslation *CBLatCov21pos_b =
1286 new TGeoTranslation("CBLatCov21pos_b",
1287 3 * fAcrylicWidth + fCBRearWallXDim / 2. +
1288 fSlopedWallZProj / (TMath::Tan(TMath::DegToRad() * 85.)) + fBPolyWidth / 2. + 0.4,
1289 -fCBFrontWallYDim_b / 2.,
1290 -(fSlopedWallZProj + 2 * (fCBTinyZDim - fAcrylicWidth) + fAcrylicWidth) + fCBTinyZDim / 2. +
1291 fBPolyWidth + fAcrylicWidth - 0.5);
1292 CBLatCov21pos_b->RegisterYourself();
1293 TGeoTranslation *CBLatCov22pos_b = new TGeoTranslation(
1294 "CBLatCov22pos_b",
1295 3 * fAcrylicWidth + fCBRearWallXDim / 2. + fSlopedWallZProj / (TMath::Tan(TMath::DegToRad() * 85.)) +
1296 fBPolyWidth / 2. + 0.4,
1297 -fCBFrontWallYDim_b / 2.,
1298 -(fSlopedWallZProj / 2 + 2 * (fCBTinyZDim - fAcrylicWidth) + fAcrylicWidth) + fCBTinyZDim / 2. + fBPolyWidth -
1299 fAcrylicWidth / 2. + fCBTinyZDim / 2. + fAcrylicWidth - 0.5);
1300 CBLatCov22pos_b->RegisterYourself();
1301 TGeoTranslation *CBLatCov23pos_b = new TGeoTranslation(
1302 "CBLatCov23pos_b", (fCBRearWallXDim - fAcrylicWidth) / 2. - fAcrylicWidth / 2. + fBPolyWidth / 2. + 0.5,
1303 -fCBFrontWallYDim_b / 2., -fCBTinyZDim / 2. + fBPolyWidth + fBPolyWidth / 2. + 0.5);
1304 CBLatCov23pos_b->RegisterYourself();
1305 TGeoTranslation *CBFrontCov1pos_b = new TGeoTranslation(
1306 "CBFrontCov1pos_b", -fCBRearWallXDim_b / 2. - fCBExtraXDim_b + fBPolyWidth + fCBFrontWallXDim_b / 2.,
1307 -fCBFrontWallYDim_b / 2., -fSlopedWallZProj - (fCBTinyZDim - fAcrylicWidth + fBPolyWidth));
1308 CBFrontCov1pos_b->RegisterYourself();
1309 TGeoTranslation *CBRearCov11pos_b =
1310 new TGeoTranslation("CBRearCov11pos_b", +fBPolyWidth / 2., -fCBFrontWallYDim_b / 2., 0);
1311 CBRearCov11pos_b->RegisterYourself();
1312 TGeoTranslation *CBRearCov12pos_b =
1313 new TGeoTranslation("CBRearCov12pos_b", -fCBRearWallXDim_b / 2. + fAcrylicWidth / 2. - 0.5,
1314 -fCBFrontWallYDim_b / 2., fAcrylicWidth / 2. + (fCBExtraZDim - 2 * fAcrylicWidth) / 2.);
1315 CBRearCov12pos_b->RegisterYourself();
1316 TGeoTranslation *CBRearCov13pos_b =
1317 new TGeoTranslation("CBRearCov13pos_b", -fCBRearWallXDim_b / 2. + fAcrylicWidth - fCBExtraXDim / 2 + fBPolyWidth,
1318 -fCBFrontWallYDim_b / 2., fCBExtraZDim - fAcrylicWidth);
1319 CBRearCov13pos_b->RegisterYourself();
1320
1321 // Borated Polyethylene mother volume definition
1322 TGeoCompositeShape *COLDBOXB = new TGeoCompositeShape(
1323 "COLDBOXB", "CBRearWall_b+(CBTiny1_b:Tiny1pos_b)+(CBExtraLat_b:CBExtraLatpos_b)+(CBWallSlope_b:SlopeWallpos_b)+("
1324 "CBTiny2_b:Tiny2pos_b)+(CBExtraFront_b:CBExtraFrontpos_b)+(CBLateral_b:CBLateralpos_b)+(CBFrontWall_"
1325 "b:FrontWallpos_b)+(LatCov1_b:CBLatCov1pos_b)+(LatCov21_b:CBLatCov21pos_b)+(LatCov22_b:CBLatCov22pos_"
1326 "b)+(LatCov23_b:CBLatCov23pos_b)+(FrontCov1_b:CBFrontCov1pos_b)+(RearCov11_b:CBRearCov11pos_b)+("
1327 "RearCov12_b:CBRearCov12pos_b)+(RearCov13_b:CBRearCov13pos_b)");
1328 TGeoVolume *volCOLDBOXB = new TGeoVolume("volCOLDBOXB", COLDBOXB, Bor30Poly);
1329
1330 // Acrylic Roof shape definition
1331 Double_t Roof4_averts[8][2];
1332 Roof4_averts[0][0] = 0.;
1333 Roof4_averts[0][1] = 0.;
1334 Roof4_averts[1][0] = 0.;
1335 Roof4_averts[1][1] = fAcrylicWidth;
1336 Roof4_averts[2][0] = fSlopedWallZProj * (TMath::Tan(TMath::DegToRad() * 15.));
1337 Roof4_averts[2][1] = fAcrylicWidth;
1338 Roof4_averts[3][0] = fSlopedWallZProj * (TMath::Tan(TMath::DegToRad() * 15.));
1339 Roof4_averts[3][1] = 0;
1340 Roof4_averts[4][0] = 0;
1341 Roof4_averts[4][1] = 0;
1342 Roof4_averts[5][0] = 0;
1343 Roof4_averts[5][1] = fAcrylicWidth;
1344 Roof4_averts[6][0] = 0;
1345 Roof4_averts[6][1] = fAcrylicWidth;
1346 Roof4_averts[7][0] = 0;
1347 Roof4_averts[7][1] = 0;
1348
1349 TGeoBBox *CBRoof1_a = new TGeoBBox("CBRoof1_a", fCBExtraXDim / 2., fAcrylicWidth / 2., fCBLatWallZDim / 2.);
1350 TGeoBBox *CBRoof2_a = new TGeoBBox("CBRoof2_a", (fCBRearWallXDim - fAcrylicWidth) / 2., fAcrylicWidth / 2.,
1351 (fCBLatWallZDim - fCBExtraZDim + fAcrylicWidth) / 2.);
1352 TGeoBBox *CBRoof3_a = new TGeoBBox("CBRoof3_a", (fSlopedWallZProj * (TMath::Tan(TMath::DegToRad() * 15.))) / 2.,
1353 fAcrylicWidth / 2., fCBTinyZDim / 2.);
1354 TGeoArb8 *CBRoof4_a = new TGeoArb8("CBRoof4_a", fSlopedWallZProj / 2., (Double_t *)Roof4_averts);
1355
1356 TGeoTranslation *Roof1_apos =
1357 new TGeoTranslation("Roof1_apos", -(fCBRearWallXDim - fAcrylicWidth) / 2. - fCBExtraXDim / 2., 0,
1358 fCBExtraZDim / 2. - fAcrylicWidth / 2.);
1359 Roof1_apos->RegisterYourself();
1360 TGeoTranslation *Roof3_apos = new TGeoTranslation(
1361 "Roof3_apos",
1362 (fCBRearWallXDim - fAcrylicWidth) / 2. + (fSlopedWallZProj * (TMath::Tan(TMath::DegToRad() * 15.))) / 2., 0,
1363 -(fCBLatWallZDim - fCBExtraZDim + fAcrylicWidth) / 2. + fCBTinyZDim / 2.);
1364 Roof3_apos->RegisterYourself();
1365 TGeoTranslation *Roof4_apos =
1366 new TGeoTranslation("Roof4_apos", (fCBRearWallXDim - fAcrylicWidth) / 2., -fAcrylicWidth / 2.,
1367 -(fCBLatWallZDim - fCBExtraZDim + fAcrylicWidth) / 2. + fCBTinyZDim + fSlopedWallZProj / 2.);
1368 Roof4_apos->RegisterYourself();
1369
1370 // Acrylic roof volume definition
1371 TGeoCompositeShape *CBRoof_a = new TGeoCompositeShape(
1372 "CBRoof_a", "CBRoof2_a+(CBRoof1_a:Roof1_apos)+(CBRoof3_a:Roof3_apos)+(CBRoof4_a:Roof4_apos)");
1373 TGeoVolume *volCBRoof_a = new TGeoVolume("volCBRoof_a", CBRoof_a, Acrylic);
1374
1375 // Borated Polythylene Roof shape definition
1376 Double_t Roof4_bverts[8][2];
1377 Roof4_bverts[0][0] = 0.;
1378 Roof4_bverts[0][1] = 0.;
1379 Roof4_bverts[1][0] = 0.;
1380 Roof4_bverts[1][1] = fBPolyWidth;
1381 Roof4_bverts[2][0] = fSlopedWallZProj * (TMath::Tan(TMath::DegToRad() * 15.));
1382 Roof4_bverts[2][1] = fBPolyWidth;
1383 Roof4_bverts[3][0] = fSlopedWallZProj * (TMath::Tan(TMath::DegToRad() * 15.));
1384 Roof4_bverts[3][1] = 0;
1385 Roof4_bverts[4][0] = 0;
1386 Roof4_bverts[4][1] = 0;
1387 Roof4_bverts[5][0] = 0;
1388 Roof4_bverts[5][1] = fBPolyWidth;
1389 Roof4_bverts[6][0] = 0;
1390 Roof4_bverts[6][1] = fBPolyWidth;
1391 Roof4_bverts[7][0] = 0;
1392 Roof4_bverts[7][1] = 0;
1393
1394 TGeoBBox *CBRoof1_b = new TGeoBBox("CBRoof1_b", fCBExtraXDim_b / 2., fBPolyWidth / 2., fCBLatWallZDim_b / 2.);
1395 TGeoBBox *CBRoof2_b = new TGeoBBox("CBRoof2_b", (fCBRearWallXDim_b - fBPolyWidth) / 2. + fBPolyWidth / 2.,
1396 fBPolyWidth / 2., (fCBLatWallZDim_b - fCBExtraZDim + fAcrylicWidth) / 2.);
1397 TGeoBBox *CBRoof3_b = new TGeoBBox("CBRoof3_b", (fSlopedWallZProj * (TMath::Tan(TMath::DegToRad() * 15.))) / 2.,
1398 fBPolyWidth / 2., (fCBTinyZDim - fAcrylicWidth) / 2.);
1399 TGeoArb8 *CBRoof4_b = new TGeoArb8("CBRoof4_b", fSlopedWallZProj / 2., (Double_t *)Roof4_bverts);
1400
1401 TGeoTranslation *Roof1_bpos =
1402 new TGeoTranslation("Roof1_bpos", -(fCBRearWallXDim - fAcrylicWidth) / 2. - fCBExtraXDim_b / 2., 0,
1403 fCBExtraZDim / 2. - fAcrylicWidth / 2.);
1404 Roof1_bpos->RegisterYourself();
1405 TGeoTranslation *Roof3_bpos = new TGeoTranslation(
1406 "Roof3_bpos",
1407 (fCBRearWallXDim - fAcrylicWidth) / 2. + (fSlopedWallZProj * (TMath::Tan(TMath::DegToRad() * 15.))) / 2., 0,
1408 -(fCBLatWallZDim_b - fCBExtraZDim + fAcrylicWidth) / 2. + (fCBTinyZDim - fAcrylicWidth) / 2.);
1409 Roof3_bpos->RegisterYourself();
1410 TGeoTranslation *Roof4_bpos = new TGeoTranslation(
1411 "Roof4_bpos", (fCBRearWallXDim - fAcrylicWidth) / 2., -fBPolyWidth / 2.,
1412 -(fCBLatWallZDim_b - fCBExtraZDim + fAcrylicWidth) / 2. + fCBTinyZDim - fAcrylicWidth + fSlopedWallZProj / 2.);
1413 Roof4_bpos->RegisterYourself();
1414
1415 // Borated Polyethylene volume definition
1416 TGeoCompositeShape *CBRoof_b = new TGeoCompositeShape(
1417 "CBRoof_b", "CBRoof2_b+(CBRoof1_b:Roof1_bpos)+(CBRoof3_b:Roof3_bpos)+(CBRoof4_b:Roof4_bpos)");
1418 TGeoVolume *volCBRoof_b = new TGeoVolume("volCBRoof_b", CBRoof_b, Bor30Poly);
1419
1420 // Volumes positioning
1421 TGeoVolumeAssembly *volColdBox = new TGeoVolumeAssembly("volColdBox");
1422 volCOLDBOXA->SetLineColor(kGray - 1);
1423 volCOLDBOXA->SetTransparency(60);
1424 volCOLDBOXB->SetLineColor(kGray - 1);
1425 volCOLDBOXB->SetTransparency(60);
1426 volCBRoof_a->SetLineColor(kGray - 1);
1427 volCBRoof_a->SetTransparency(60);
1428 volCBRoof_b->SetLineColor(kGray - 1);
1429 volCBRoof_b->SetTransparency(60);
1430
1431 volColdBox->AddNode(volCOLDBOXA, 0, 0);
1432 volColdBox->AddNode(volCOLDBOXB, 0,
1433 new TGeoTranslation(-fBPolyWidth - fAcrylicWidth / 2., -fBPolyWidth / 2.,
1434 -fAcrylicWidth / 2. - fBPolyWidth / 2.));
1435 volColdBox->AddNode(volCBRoof_a, 0,
1436 new TGeoTranslation(fAcrylicWidth / 2., fCBFrontWallYDim / 2. + fAcrylicWidth / 2.,
1437 -(fCBLatWallZDim - fCBExtraZDim + fAcrylicWidth) / 2. + fAcrylicWidth / 2.));
1438 volColdBox->AddNode(volCBRoof_b, 0,
1439 new TGeoTranslation(-fAcrylicWidth / 2., fCBFrontWallYDim / 2. - fBPolyWidth / 2.,
1440 -(fCBLatWallZDim - fCBExtraZDim + fAcrylicWidth) / 2. + fAcrylicWidth / 2.));
1441
1442 displacement =
1443 TVector3(-37.79 - 1.40082, 44.66,
1444 368.11); // taken from MuFilter.cxx "edge_Iron[1]-TVector3(FeX, FeY, FeZ)" EDIT: 1.40082 for overlap fix
1445 tunnel->AddNode(volColdBox, 0,
1446 new TGeoTranslation(displacement.X() - (fCBRearWallXDim - fFeBlockX) / 2. + 28.5,
1447 displacement.Y() - (fCBFrontWallYDim - fFeBlockY) / 2. + 121,
1448 displacement.Z() + fAcrylicWidth - fFeBlockZ / 2. - fBPolyWidth + 1.));
1449
1450if (SND_Z<0.1){ // for H6 and H8 testbeam setup
1451 top->AddNode(detector, 0);
1452 }else{
1453 top->AddNode(detector, 0,localSND_physCS_comb);
1454 top->AddNode(tunnel , 1,new TGeoTranslation(fTunnelDX,fTunnelDY,fTunnelDZ));
1455 }
1456
1457}
Floor::conf_floats
std::map< TString, Float_t > conf_floats
Definition Floor.h:81
Floor::crossing
TVector3 crossing(TVector3 H1, TVector3 H2, TVector3 H3, TVector3 P1, TVector3 P2, TVector3 P3)
Definition Floor.cxx:143
Floor::SND_Z
Double_t SND_Z
Definition Floor.h:98
Floor::InitMedium
Int_t InitMedium(const char *name)
Definition Floor.cxx:123
ShipAna.i
int i
Definition ShipAna.py:86
ana_ShipMuon.origin
origin(sTree, it)
Definition ana_ShipMuon.py:377
charm-geometry_config.H
H
Definition charm-geometry_config.py:186
charm-geometry_config.z
z
Definition charm-geometry_config.py:7
conditionsDatabase.getconditions_condDB.detector
str detector
Definition getconditions_condDB.py:38
g4Ex_gap_mergeFiles.norm
norm
Definition g4Ex_gap_mergeFiles.py:199
geometry_config.X
X
Definition geometry_config.py:517
neutrinoFilterGoldenSample_stage2.b
b
Definition neutrinoFilterGoldenSample_stage2.py:493
neutrinoFilterGoldenSample_stage2.a
a
Definition neutrinoFilterGoldenSample_stage2.py:492
loop1
subroutine loop1
First data loop (get global labels).
Definition pede.f90:4822
loop2
subroutine loop2
Second data loop (number of derivatives, global label pairs).
Definition pede.f90:5216

◆ CopyClones()

virtual void Floor::CopyClones ( TClonesArray *  cl1,
TClonesArray *  cl2,
Int_t  offset 
)
inlinevirtual

Definition at line 54 of file Floor.h.

55 {;}

◆ crossing()

TVector3 Floor::crossing ( TVector3  H1,
TVector3  H2,
TVector3  H3,
TVector3  P1,
TVector3  P2,
TVector3  P3 
)

Definition at line 143 of file Floor.cxx.

144{
145 TMatrixD M(3,3);
146 TMatrixD P(3,1);
147 TMatrixD X(3,1);
148 std::map<int,TVector3> Row = {
149 {0,H1},{1,H2},{2,H3}
150 };
151 std::map<int,TVector3> P123 = {
152 {0,P1},{1,P2},{2,P3}
153 };
154 for( int j = 0;j<3;j+=1) {
155 P[j][0] = Row[j].Dot(P123[j]);
156 for( int k = 0;k<3;k+=1) {
157 M[j][k] = Row[j][k];
158 }
159 }
160 auto Minv = M.Invert();
161 X.Mult(Minv,P);
162 TVector3 Crossing(X[0][0],X[1][0],X[2][0]);
163 return Crossing;
164}
conditionsDatabase.snd_conDB.j
int j
Definition snd_conDB.py:197

◆ EndOfEvent()

void Floor::EndOfEvent ( )
virtual

Definition at line 98 of file Floor.cxx.

99{
100 fFloorPointCollection->Clear();
101 fTotalEloss=0;
102}

◆ FinishPrimary()

virtual void Floor::FinishPrimary ( )
inlinevirtual

Definition at line 58 of file Floor.h.

58{;}

◆ FinishRun()

virtual void Floor::FinishRun ( )
inlinevirtual

Definition at line 59 of file Floor.h.

59{;}

◆ GetCollection()

TClonesArray * Floor::GetCollection ( Int_t  iColl) const
virtual

Gets the produced collections

Definition at line 1478 of file Floor.cxx.

1479{
1480 if (iColl == 0) { return fFloorPointCollection; }
1481 else { return NULL; }
1482}

◆ GetConfParF()

Float_t Floor::GetConfParF ( TString  name)
inline

Definition at line 31 of file Floor.h.

31{return conf_floats[name];}

◆ GetConfParI()

Int_t Floor::GetConfParI ( TString  name)
inline

Definition at line 32 of file Floor.h.

32{return conf_ints[name];}
Floor::conf_ints
std::map< TString, Int_t > conf_ints
Definition Floor.h:82

◆ GetConfParS()

TString Floor::GetConfParS ( TString  name)
inline

Definition at line 33 of file Floor.h.

33{return conf_strings[name];}
Floor::conf_strings
std::map< TString, TString > conf_strings
Definition Floor.h:83

◆ Initialize()

void Floor::Initialize ( )
virtual

Initialization of the detector is done here

Definition at line 118 of file Floor.cxx.

119{
120 FairDetector::Initialize();
121}

◆ InitMedium()

Int_t Floor::InitMedium ( const char *  name)

Definition at line 123 of file Floor.cxx.

124{
125 static FairGeoLoader *geoLoad=FairGeoLoader::Instance();
126 static FairGeoInterface *geoFace=geoLoad->getGeoInterface();
127 static FairGeoMedia *media=geoFace->getMedia();
128 static FairGeoBuilder *geoBuild=geoLoad->getGeoBuilder();
129
130 FairGeoMedium *ShipMedium=media->getMedium(name);
131
132 if (!ShipMedium)
133 {
134 Fatal("InitMedium","Material %s not defined in media file.", name);
135 return -1111;
136 }
137 TGeoMedium* medium=gGeoManager->GetMedium(name);
138 if (medium!=NULL)
139 return ShipMedium->getMediumIndex();
140 return geoBuild->createMedium(ShipMedium);
141}

◆ MakeSensitive()

void Floor::MakeSensitive ( )
inline

Definition at line 73 of file Floor.h.

73{fMakeSensitive=true;} // make Tunnel sensitive

◆ PostTrack()

virtual void Floor::PostTrack ( )
inlinevirtual

Definition at line 61 of file Floor.h.

61{;}

◆ PreTrack()

void Floor::PreTrack ( )
virtual

Definition at line 104 of file Floor.cxx.

104 {
105 TLorentzVector pos;
106 gMC->TrackPosition(pos);
107 if (fFastMuon && pos.Z()<fzPos && TMath::Abs(gMC->TrackPid())!=13){
108 gMC->StopTrack();
109 }
110 TLorentzVector mom;
111 gMC->TrackMomentum(mom);
112 if ( (mom.E()-mom.M() )<fEmin && pos.Z()<fzPos){
113 gMC->StopTrack();
114 return;
115 }
116}
run_simPgun.mom
mom
Definition run_simPgun.py:84
run_simPgun.pos
pos
Definition run_simPgun.py:85
run_simScript.gMC
gMC
Definition run_simScript.py:590

◆ ProcessHits()

Bool_t Floor::ProcessHits ( FairVolume *  v = 0)
virtual

this method is called for each step during simulation (see FairMCApplication::Stepping())

This method is called from the MC stepping

Definition at line 58 of file Floor.cxx.

59{
61 //Set parameters at entrance of volume. Reset ELoss.
62 if ( gMC->IsTrackEntering() ) {
63 fELoss = 0.;
64 fTime = gMC->TrackTime() * 1.0e09;
65 fLength = gMC->TrackLength();
66 gMC->TrackPosition(fPos);
67 gMC->TrackMomentum(fMom);
68 }
69 // Sum energy loss for all steps in the active volume
70 fELoss += gMC->Edep();
71
72 // Create vetoPoint at exit of active volume
73 if ( gMC->IsTrackExiting() ||
74 gMC->IsTrackStop() ||
75 gMC->IsTrackDisappeared() ) {
76 fTrackID = gMC->GetStack()->GetCurrentTrackNumber();
77
78 Int_t veto_uniqueId = 47.;
79
80 TParticle* p=gMC->GetStack()->GetCurrentTrack();
81 Int_t pdgCode = p->GetPdgCode();
82 TLorentzVector Pos;
83 gMC->TrackPosition(Pos);
84 TLorentzVector Mom;
85 gMC->TrackMomentum(Mom);
86 Double_t xmean = (fPos.X()+Pos.X())/2. ;
87 Double_t ymean = (fPos.Y()+Pos.Y())/2. ;
88 Double_t zmean = (fPos.Z()+Pos.Z())/2. ;
89 LOG(DEBUG)<< veto_uniqueId << " :(" << xmean << ", " << ymean << ", " << zmean << "): " << fELoss;
90 AddHit(fTrackID, veto_uniqueId, TVector3(xmean, ymean, zmean),
91 TVector3(fMom.Px(), fMom.Py(), fMom.Pz()), fTime, fLength,
92 fELoss,pdgCode,TVector3(Pos.X(), Pos.Y(), Pos.Z()),TVector3(Mom.Px(), Mom.Py(), Mom.Pz()) );
93 }
94
95 return kTRUE;
96}
Floor::AddHit
vetoPoint * AddHit(Int_t trackID, Int_t detID, TVector3 pos, TVector3 mom, Double_t time, Double_t length, Double_t eLoss, Int_t pdgcode, TVector3 Lpos, TVector3 Lmom)
Definition Floor.cxx:1459
Floor::fELoss
Double_t fELoss
thickness
Definition Floor.h:93

◆ Register()

void Floor::Register ( )
virtual

Registers the produced collections in FAIRRootManager.

Definition at line 1471 of file Floor.cxx.

1472{
1473
1474 FairRootManager::Instance()->Register("vetoPoint", "veto",
1475 fFloorPointCollection, kTRUE);
1476}

◆ Reset()

void Floor::Reset ( )
virtual

has to be called after each event to reset the containers

Definition at line 1483 of file Floor.cxx.

1484{
1485 fFloorPointCollection->Clear();
1486}

◆ SetConfPar() [1/3]

void Floor::SetConfPar ( TString  name,
Float_t  value 
)
inline

Definition at line 28 of file Floor.h.

28{conf_floats[name]=value;}

◆ SetConfPar() [2/3]

void Floor::SetConfPar ( TString  name,
Int_t  value 
)
inline

Definition at line 29 of file Floor.h.

29{conf_ints[name]=value;}

◆ SetConfPar() [3/3]

void Floor::SetConfPar ( TString  name,
TString  value 
)
inline

Definition at line 30 of file Floor.h.

30{conf_strings[name]=value;}

◆ SetEmin()

void Floor::SetEmin ( float  E)
inline

Definition at line 70 of file Floor.h.

70{fEmin = E;} // set min kin energy for tracking

◆ SetFastMuon()

void Floor::SetFastMuon ( )
inline

Definition at line 72 of file Floor.h.

72{fFastMuon=true;} // kill all tracks except of muons

◆ SetSpecialPhysicsCuts()

virtual void Floor::SetSpecialPhysicsCuts ( )
inlinevirtual

Definition at line 56 of file Floor.h.

56{;}

◆ SetZmax()

void Floor::SetZmax ( float  Z)
inline

Definition at line 71 of file Floor.h.

71{fzPos = Z;} // set max z position for which energy cut is applied.
geometry_config.Z
Z
Definition geometry_config.py:519

Member Data Documentation

◆ conf_floats

std::map<TString,Float_t> Floor::conf_floats
private

configuration parameters

Definition at line 81 of file Floor.h.

◆ conf_ints

std::map<TString,Int_t> Floor::conf_ints
private

Definition at line 82 of file Floor.h.

◆ conf_strings

std::map<TString,TString> Floor::conf_strings
private

Definition at line 83 of file Floor.h.

◆ fELoss

Double_t Floor::fELoss
private

thickness

Definition at line 93 of file Floor.h.

◆ fEmin

Float_t Floor::fEmin
private

Definition at line 100 of file Floor.h.

◆ fFastMuon

Bool_t Floor::fFastMuon
private

Definition at line 99 of file Floor.h.

◆ fFloorPointCollection

TClonesArray* Floor::fFloorPointCollection
private

material

container for data points

Definition at line 97 of file Floor.h.

◆ fLength

Double_t Floor::fLength
private

time

Definition at line 90 of file Floor.h.

◆ fMakeSensitive

Bool_t Floor::fMakeSensitive
private

Definition at line 101 of file Floor.h.

◆ fMaterial

TString Floor::fMaterial
private

Definition at line 95 of file Floor.h.

◆ fMom

TLorentzVector Floor::fMom
private

position at entrance

Definition at line 88 of file Floor.h.

◆ fPos

TLorentzVector Floor::fPos
private

volume id

Definition at line 87 of file Floor.h.

◆ fThick

Double_t Floor::fThick
private

zPos

Definition at line 92 of file Floor.h.

◆ fTime

Double_t Floor::fTime
private

momentum at entrance

Definition at line 89 of file Floor.h.

◆ fTotalEloss

Double_t Floor::fTotalEloss
private

Definition at line 94 of file Floor.h.

◆ fTrackID

Int_t Floor::fTrackID
private

Definition at line 85 of file Floor.h.

◆ fVolumeID

Int_t Floor::fVolumeID
private

track index

Definition at line 86 of file Floor.h.

◆ fzPos

Double_t Floor::fzPos
private

length

Definition at line 91 of file Floor.h.

◆ SND_Z

Double_t Floor::SND_Z
private

Definition at line 98 of file Floor.h.

The documentation for this class was generated from the following files: