CMSSW Reference Manual

00001 #ifndef Point5MaterialMap_cc
00002 #define Point5MaterialMap_cc
00003 
00004 #include <iostream>
00005 
00006 #include "GeneratorInterface/CosmicMuonGenerator/interface/CosmicMuonParameters.h"
00007 
00008 
00009 
00010 inline int inPlug(double vx, double vy, double vz,
00011                   double PlugVx = PlugOnShaftVx, double PlugVz = PlugOnShaftVz) {
00012   if (vy > SurfaceOfEarth && vy < SurfaceOfEarth + PlugWidth) {
00013     if (vx > PlugVx - PlugXlength/2. && vx < PlugVx + PlugXlength/2. &&
00014         vz > PlugVz - PlugZlength/2. && vz < PlugVz + PlugZlength/2.) return Plug;
00015     if (vz >= PlugVz - PlugZlength/2. - PlugNoseZlength && vz < PlugVz - PlugZlength/2. &&
00016         vx > PlugVx - PlugNoseXlength/2. && vx < PlugVx + PlugNoseXlength/2.) return Plug;
00017   }
00018   return Unknown;
00019 }
00020 
00021 
00022 inline int inAirAfterPlug(double vx, double vy, double vz) {
00023   // particles above surface of earth
00024   if (vy >= SurfaceOfEarth) return Air;
00025   
00026   // CMS cavern (UXC 55)
00027   if (fabs(vz) < 26548. && sqrt((vx*1.1576)*(vx*1.1576) + vy*vy) < 15460. &&
00028       vy > -8762) return Air;
00029   
00030   // access shaft (PX 56)
00031   if (vy > 0. && vy < (SurfaceOfEarth-2250.) && 
00032       sqrt(vx*vx + (vz-Z_PX56)*(vz-Z_PX56)) < 10250.) return Air;
00033   
00034   //surface hall ground floor
00035   if (vy >= SurfaceOfEarth-2250. && vy < SurfaceOfEarth) {
00036     if (sqrt(vx*vx + (vz-Z_PX56)*(vz-Z_PX56)) < 10250.
00037         && vz-Z_PX56 > -7000. && vz-Z_PX56 < 7000.) return Air;
00038     if (vx > -2400. && vx < 2400. && vz-Z_PX56 >= -9800. && vz-Z_PX56 < -7000. )
00039       return Air;
00040   }
00041   
00042     // Shaft (PM 54)
00043   if (vy > 3233. && vy < (SurfaceOfEarth) && 
00044       sqrt((vx-26600.)*(vx-26600.) + (vz-30100.-Z_PX56)*(vz-30100.-Z_PX56)) < 6050.)
00045     //sqrt((vx-5000.)*(vx-5000.) + (vz-30100.-Z_PX56)*(vz-30100.-Z_PX56)) < 10050.)//@@@@@
00046     return Air;
00047   
00048   
00049   // Shaft (PM 56)
00050   if (vy > -3700. && vy < (SurfaceOfEarth) &&
00051       sqrt((vx-18220.)*(vx-18220.) + (vz+24227.-Z_PX56)*(vz+24227.-Z_PX56)) < 3550.)
00052       return Air;
00053   
00054   //Service cavern USC 55
00055   if (vz > -22050. && vz < 62050. &&
00056       sqrt((vx-29550.)*(vx-29550.) + (vy-3233.)*(vy-3233.)) < 9050. && 
00057       vy > -3650.) return Air;
00058   
00059   //UXC55 cavern endcap beam openings
00060   if (((vz >= -29278. && vz < -26548.) || (vz >= 26548 && vz <= 29278.)) &&
00061       sqrt(vy*vy + vx*vx) < 800.) return Air; 
00062   
00063   //Pillar access between CMS collision cavern and service cavern
00064   if (vx > 10000. && vx < 20500.+3050. && //TX54 galerie
00065       vz > 14460. && vz < 16260. &&//estimated wall thickness 1m
00066       vy > -8680. && vy < -1507.) return Air;
00067   
00068   if (vx > 10000. && vx < 16500. && //TX54 galerie 30.9 degree opening to UXC55
00069       vy > -8680. && vy < -1507. && //@vx=13300 delta(vz)=1865., 2vx=1650. delta(vz)=0
00070       vz > 14460.-1865.*(16500.-vx)/(3200.) && vz <= 14460.) return Air;
00071   
00072   if (vx > 13300. && vx < 20500.+3050. && //TX54 galerie
00073       vz > 14460. && vz < 16260. && //estimated wall thickness 1m
00074       vy > -8680. && vy < -1507.) return Air;
00075   
00076   if (vx > 26600-6050. && vx < 20500.+3050. && //TX54 going up in sewrvice cavern
00077       vz > 14460. && vz < 16260. && //estimated wall thickness 1m
00078       vy >= -1507. && vy < 1000.) return Air;
00079   
00080   //R56, LHC beam tunnel East
00081   if (vz > -85000. && vz <= -29278. && //UJ57 junction to UXC55 cavern
00082       sqrt(vy*vy + (vx-350.)*(vx-350.)) < 1900. &&
00083       vy > -1000.) return Air;
00084   
00085   //R54, LHC beam tunnel West
00086   if (vz >= 29278. && vz < 63000. && //UJ57 junction to UXC55 cavern
00087       sqrt(vy*vy + (vx-350.)*(vx-350.)) < 1900. &&
00088       vy > -1000.) return Air;
00089   
00090   //UJ56 cavern
00091   if (vz > -58875. && vz < -33927. &&
00092       sqrt(vy*vy + (vx-4450.)*(vx-4450.)) < 6750. && vy > -1000. && //and beam shielding
00093       !(vx > 2250. && vx < 4250. && (vz > -(33927.+18000.) || vz < -(33927.+19500.))) && 
00094       !(vx >= 4250. && vx < 6650. && vz > -(33927.+18000.) && vz < -(33927.+16000.)))
00095     return Air;
00096   
00097   //connection between PM56 shaft and UJ56 cavern
00098   if (vx > 9000. && vx < 18220. &&
00099       sqrt((vy-50.)*(vy-50.)+(vz+24227.-Z_PX56)*(vz+24227.-Z_PX56)) < 3550.
00100       && vy > -1000.) return Air;
00101   
00102   return Unknown;
00103   
00104 }
00105 
00106 
00107 inline int inWallAfterAir(double vx, double vy, double vz) {
00108   // phase II surface building
00109   if (vy < SurfaceOfEarth && vy >= (SurfaceOfEarth-2250.)) {
00110     if (fabs(vz-Z_PX56) < 30000. && fabs(vx) < 10950) return Wall;
00111     // foundation of crane
00112     if (fabs(vz-Z_PX56) < 9000. && fabs(vx) >= 10950 && fabs(vx) < 16950) 
00113       return Wall;
00114   }
00115   
00116   // CMS cavern (UXC 55)
00117   if (fabs(vz) < 29278. && sqrt((vx*1.1576)*(vx*1.1576) + vy*vy) < 16830. &&
00118       vy > -11762.) return Wall;
00119   
00120   // access shaft (PX 56)
00121   if (vy > 0. && vy < (SurfaceOfEarth-2250.) && //t(shaft wall)=2150.
00122       sqrt(vx*vx + (vz-Z_PX56)*(vz-Z_PX56)) < 12400.) return Wall;
00123   
00124   // Shaft (PM 54)
00125   if (vy > 3233. && vy < (SurfaceOfEarth-1000.) && //t~=t(PX56)/R(PX56)*R(PM54)
00126       sqrt((vx-26600.)*(vx-26600.) + (vz-30100.-Z_PX56)*(vz-30100.-Z_PX56)) 
00127       < 6050.+2150./10250.*6050.) return Wall;
00128   //sqrt((vx-5000.)*(vx-5000.) + (vz-30100.-Z_PX56)*(vz-30100.-Z_PX56))//@@@@@ 
00129   //< 10050.+2150./10250.*6050.) return Wall;//@@@@@
00130   else if (vy >= SurfaceOfEarth-1000. && vy < SurfaceOfEarth && //t~=t(PX56)/R(PX56)*R(PM54)
00131            sqrt((vx-26600.)*(vx-26600.) + (vz-30100.-Z_PX56)*(vz-30100.-Z_PX56)) 
00132            < 6050.+2150./10250.*6050. +1800.) return Wall;
00133   //sqrt((vx-5000.)*(vx-5000.) + (vz-30100.-Z_PX56)*(vz-30100.-Z_PX56))//@@@@@ 
00134   //< 10050.+2150./10250.*10050. +1800.) return Wall; //@@@@@@
00135   
00136   // Shaft (PM 56)
00137   if (vy > -5450. && vy < (SurfaceOfEarth-1000.) && //t~=t(PX56)/R(PX56)*R(PM56)
00138       sqrt((vx-18220.)*(vx-18220.) + (vz+24227.-Z_PX56)*(vz+24227.-Z_PX56)) 
00139       < 3550.+2150./10250.*3550.) return Wall;
00140   else if (vy > SurfaceOfEarth-1000. && vy < SurfaceOfEarth && //t~=t(PX56)/R(PX56)*R(PM56)
00141            sqrt((vx-18220.)*(vx-18220.) + (vz+24227.-Z_PX56)*(vz+24227.-Z_PX56)) 
00142            < 3550.+2150./10250.*3550. +1800.) return Wall;
00143   
00144   //Service cavern USC 55
00145   if (vz > -(22050.+1150.) && vz < (62050.+1150.) &&
00146       sqrt((vx-29550.)*(vx-29550.) + (vy-3233.)*(vy-3233.)) < 9050.+950. &&
00147       vy > -3650.-2000.) return Wall; //8762=estimate, to be checked 
00148   
00149   //Pillar between CMS collision cavern and service cavern
00150   if (vz > -29278.+1000. && vz < 29278.+1000.) {
00151     if (vy > -17985. && vy < 10410. && vx > 13300. && vx < 20500.) 
00152       return Wall;
00153     
00154     if (vy > 0. && vy < 10410. && vx > 10000. && vx <= 13300.)
00155       return Wall;
00156     
00157     if (vy > -3650.-2000. && vy < -3233. && // bottom edge between pillar and service cavern
00158         vx > 20000. && vx < 24000.) return Wall;
00159     if (vy > -11762. && vy < -5000. && // bottom edge between pillar and UXC55 cavern
00160         vx > 10500. && vx < 14000.) return Wall;
00161   }
00162   
00163   if (vy > -14000. && vy < -1450. && //TX54 galerie surrounding
00164       vz > 13460. && vz < 17260. && vx >= 20500. && vx < 24550.) 
00165     return Wall;
00166   
00167   //R56, LHC beam tunnel East
00168   if (vz > -85000. && vz < -28510.) { //UJ57 junction to UXC55 cavern
00169     if (sqrt(vy*vy + (vx-350.)*(vx-350.)) < 2250.) return Wall;
00170   }
00171   
00172   //R54, LHC beam tunnel West
00173   if (vz > 26550. && vz < 63000. && //UJ57 junction to UXC55 cavern
00174       sqrt(vy*vy + (vx-350.)*(vx-350.)) < 2250.) return Wall;
00175   
00176   //UJ56 cavern
00177   if (vz > -(58875.+500.) && vz < -(33927.-500.) &&
00178       sqrt(vy*vy + (vx-4450.)*(vx-4450.)) < (6750.+500.) && vy > -3650.) 
00179     return Wall;
00180   
00181   //connection between PM56 shaft and UJ56 cavern
00182   if (vx > 9000. && vx < 18220. &&
00183       sqrt((vy-50.)*(vy-50.)+(vz+24227.-Z_PX56)*(vz+24227.-Z_PX56)) < 3550.+500. 
00184       && vy > -3650.) return Wall;
00185   
00186   return Unknown;
00187 }
00188 
00189 
00190 
00191 inline int inClayOrRockAfterWall(double vx, double vy, double vz) {
00192   
00193   //So, it is not plug, air and wall, Check for clay
00194   if (vy >= SurfaceOfEarth - ClayWidth && vy < SurfaceOfEarth)
00195     return Clay;
00196   
00197   //So, it is not plug, air, wall and clay, Check for rock
00198   if (vy < SurfaceOfEarth - ClayWidth)
00199     return Rock;
00200   
00201   return Unknown;
00202   
00203 }
00204 
00205 
00206 
00207 inline int inClayAfterWall(double vx, double vy, double vz) {
00208   
00209   //So, it is not plug, air and wall, Check for clay
00210   if (vy >= SurfaceOfEarth - ClayWidth && vy < SurfaceOfEarth)
00211     return Clay;
00212   
00213   return Unknown;
00214   
00215 }
00216 
00217 
00218 
00219 inline int inRockAfterClay(double vx, double vy, double vz) {
00220   
00221   //So, it is not plug, air, wall and clay, Check for rock
00222   if (vy < SurfaceOfEarth - ClayWidth)
00223     return Rock;
00224   
00225   return Unknown;
00226   
00227 }
00228 
00229 
00230 
00231 
00232 
00233 inline int inMat(double vx, double vy, double vz,
00234                  double PlugVx = PlugOnShaftVx, double PlugVz = PlugOnShaftVz) {
00235   
00236   //check for Plug
00237   if (inPlug(vx, vy, vz, PlugVx, PlugVz)) return Plug;
00238   
00239   //So, it is not plug, Check for air
00240   if (inAirAfterPlug(vx, vy, vz)) return Air;
00241   
00242   //So, it is not plug and air, Check for wall
00243   if (inWallAfterAir(vx, vy, vz)) return Wall;
00244   
00245   //So, it is not plug, air and wall, Check for clay
00246   if (vy >= SurfaceOfEarth - ClayWidth && vy < SurfaceOfEarth)
00247     return Clay;
00248   
00249   //So, it is not plug, air, wall and clay, Check for rock
00250   if (vy < SurfaceOfEarth - ClayWidth)
00251     return Rock;
00252   
00253   
00254   std::cout << "Point5MaterialMap.h: Warning! No Material recognised for point: vx="
00255             << vx << " vy=" << vy << " vz=" << vz << std::endl;
00256   //Something went wrong
00257   return Unknown;
00258   
00259 }
00260 
00261 
00262 
00263 
00264 #endif
Point5MaterialMap.cc
