00001 #ifndef CosmicMuonParameters_h 00002 #define CosmicMuonParameters_h 00003 // 00004 // Parameters for CosmicMuonGenerator by droll (05/DEC/2005) 00005 // 00006 // 00007 // added plug and clay(moraine) specific constants, sonne (15/Jan/2009) 00008 // 00009 #include "TMath.h" 00010 00011 // flags 00012 const bool Debug = false; // debugging printout 00013 const bool EventDisplay = true; // display single events (if ROOT_INTERACTIVE is defined as true) 00014 00015 // algorithmic constants 00016 const double MinStepSize = 10.; // minimal propagation step size [mm] must be small compared to target size 00017 // mathematical constants 00018 const double Pi = acos(-1.); // [rad] 00019 const double TwoPi = 2.0*Pi; // [rad] 00020 const double Deg2Rad = Pi/180.; // [deg] -> [rad] 00021 const double Rad2Deg = 180./Pi; // [rad] -> [deg] 00022 // physical constants 00023 const double SpeedOfLight = 299.792458; // [mm/ns] 00024 const double MuonMass = 0.105658357; // [GeV/c^2] 00025 //const double ChargeFrac = 0.545454545; // n(mu+)/n(mu-) ~ 1.2 defined in CMSCGEN 00026 // geometry 00027 const double SurfaceOfEarth = 88874.; // Y-distance to surface of earth [mm] 00028 const double Z_PX56 = -14000.; // [mm] Z position of PX56 centre [mm] 00029 /* 00030 // densities of materials 00031 const double RhoAir = 0.001214; // [g cm^-3] 00032 const double RhoWall = 2.5; // [g cm^-3] 00033 const double RhoRock = 2.5; // [g cm^-3] 00034 const double RhoClay = 2.3; // [g cm^-3] 00035 const double RhoPlug = 2.5; // [g cm^-3] 00036 */ 00037 // width of clay layer between surface and rock 00038 const double DefaultClayWidth = 50000.; // [mm] 00039 00040 //plug constants 00041 const double PlugWidth = 2250.; // [mm] 00042 const double PlugXlength = 20600.; // [mm] 00043 const double PlugZlength = 16000.; // [mm] 00044 const double PlugNoseXlength = 6400.; // [mm] 00045 const double PlugNoseZlength = 1800.; // [mm] 00046 const double PlugOnShaftVx = 0.; // [mm] 00047 const double PlugOnShaftVz = Z_PX56; // [mm] 00048 00049 // cylinder around CMS (with R, +-Z) 00050 // WARNING: These values will be set to tracker-only setup if "TrackerOnly=true" in .cfg-file. 00051 // This means R=1200 and Z=2800, no material or B-field outside is considered 00052 const double RadiusCMS = 8000.; // [mm] 00053 const double Z_DistCMS = 15000.; // [mm] 00054 const double RadiusTracker = 1200.; // [mm] 00055 const double Z_DistTracker = 2800.; // [mm] 00056 // cylinder actually used in the code 00057 //const double RadiusTarget = RadiusCMS; // [mm] //now controlled by cfg-file!!! 00058 //const double Z_DistTarget = Z_DistCMS; // [mm] //now controlled by cfg-file!!! 00059 00060 00061 00062 //define different materials 00063 enum {Unknown=0, Plug, Wall, Air, Clay, Rock}; 00064 00065 00066 //Parameters for upward muons from neutrinos 00067 const double N_A = 6.022e23; //mol^-1 00068 const double alpha = 2.; //MeV/(g/cm^2) 00069 const double beta = 3.9e-6; //(g/cm^2)^-1 00070 const double epsilon = alpha/beta; 00071 const double Rearth = 6370.e6; //mm 00072 00073 //Multi Muon relevant parameters 00074 const double NorthCMSzDeltaPhi = 3./8.*Pi; //rad (Pi/2 if CMS -x = North) 00075 const int max_Trials = 200000; 00076 00077 00078 #endif