emc.h 2.0 KB

12345678910111213141516171819202122232425262728293031323334353637383940414243444546474849505152
  1. #include "TMath.h"
  2. #ifndef EMC_H
  3. #define EMC_H 1
  4. namespace emc {
  5. Double_t z = 0.0;
  6. Double_t draw_st_pos = 90.;
  7. UShort_t nModInSuperModByZ = 3;
  8. UShort_t nModInSuperModByPhi = 3;
  9. UShort_t nSuperModInRow = 4;
  10. Double_t h = 407.8; //lenght of all module (Pb+Sc)*n + plastics
  11. Double_t box_pb_h = 0.3; //thickness of Pb
  12. Double_t box_sc_h = 1.5; //thickness of Sc
  13. Double_t box_pad = 10; //thickness of base of the module
  14. Double_t base = 40.; //width of Channel
  15. Int_t NScPb = 221; //number of layers Sc+Pb
  16. Int_t NmodZ = 23; //modules in one sector
  17. Int_t NSector = 28; //number of sector in 1/2 Barrel
  18. Int_t NmodPHY = 112; //number of modules in Phi plane
  19. Double_t emc1_box1_d = 98.; //size of down side of Trap
  20. Double_t stThM = 0.5; //thickness of steel ~between modules - 0.5mm; ~between sector 1mm
  21. Double_t angleSector = 360. / NSector; //rotation angel in Phi plane for sector
  22. Double_t angleSuperModule = angleSector / nSuperModInRow; //rotation angel in Phi plane for super-module
  23. Double_t angleModule = angleSuperModule / nModInSuperModByPhi; //rotation angel in Phi plane for module
  24. Double_t emc1_box1_z = base * nModInSuperModByZ; //120 mm
  25. Double_t emc1Chamber_z = NmodZ*emc1_box1_z; //lenght "z" in 1/2 Barrel
  26. Double_t emc1Chamber_z_l = emc1Chamber_z * 2; //lenght "z" in all Barrel
  27. Double_t emc1Chamber_z_th = emc1Chamber_z / 2; //"z" start position of 1/2 Barrel
  28. Double_t lenghtM = NScPb*(box_pb_h + box_sc_h); //real lenght of module
  29. Double_t inr = 1748; //Inner rarius
  30. Double_t outr = inr + h; //Outr radius -> module + fron's plastic and end's plastic
  31. Double_t outro = inr + 464.4; //modul+base+space for electronics
  32. //Double_t angleModule; //rotation module in sector
  33. Double_t angleSteel; //rotation steel in sector
  34. Double_t angleModulePad; //cut angle in module
  35. Double_t angleTrap = TMath::ATan(((emc1_box1_z - emc1_box1_d) / 2) / lenghtM);
  36. Double_t apb = box_pb_h*TMath::Tan(angleTrap);
  37. Double_t asc = box_sc_h*TMath::Tan(angleTrap);
  38. Float_t cs, si, dx, dx1, dxM, dxSt, dAlf, Alf, dAlfM, AlfM, dAlfSt, AlfSt, dA;
  39. Int_t sch;
  40. }
  41. #endif