CosmicMuonParameters.h

Go to the documentation of this file.

#ifndef CosmicMuonParameters_h
#define CosmicMuonParameters_h
//
// Parameters for CosmicMuonGenerator by droll (05/DEC/2005)
//
//
// added plug and clay(moraine) specific constants, sonne (15/Jan/2009)
//
#include "TMath.h"

// flags
const bool Debug = false; // debugging printout
const bool EventDisplay = true; // display single events (if ROOT_INTERACTIVE is defined as true)

// algorithmic constants
const double MinStepSize = 10.; // minimal propagation step size [mm] must be small compared to target size
// mathematical constants
const double Pi = acos(-1.); // [rad]
const double TwoPi = 2.0*Pi; // [rad]
const double Deg2Rad = Pi/180.; // [deg] -> [rad]
const double Rad2Deg = 180./Pi; // [rad] -> [deg]
// physical constants
const double SpeedOfLight = 299.792458; // [mm/ns]
const double MuonMass = 0.105658357; // [GeV/c^2]
//const double ChargeFrac = 0.545454545; // n(mu+)/n(mu-) ~ 1.2 defined in CMSCGEN
// geometry
const double SurfaceOfEarth = 88874.; // Y-distance to surface of earth [mm]
const double Z_PX56 = -14000.; // [mm] Z position of PX56 centre [mm]
/*
// densities of materials
const double RhoAir  = 0.001214; // [g cm^-3]
const double RhoWall = 2.5; // [g cm^-3]
const double RhoRock = 2.5; // [g cm^-3]
const double RhoClay = 2.3; // [g cm^-3]
const double RhoPlug = 2.5; // [g cm^-3]
*/
// width of clay layer between surface and rock
const double DefaultClayWidth = 50000.; // [mm]

//plug constants
const double PlugWidth = 2250.; // [mm]
const double PlugXlength = 20600.; // [mm]
const double PlugZlength = 16000.; // [mm]
const double PlugNoseXlength = 6400.; // [mm]
const double PlugNoseZlength = 1800.; // [mm]
const double PlugOnShaftVx = 0.; // [mm]
const double PlugOnShaftVz = Z_PX56; // [mm]

// cylinder around CMS (with R, +-Z)
// WARNING: These values will be set to tracker-only setup if "TrackerOnly=true" in .cfg-file. 
// This means R=1200 and Z=2800, no material or B-field outside is considered
const double RadiusCMS =  8000.; // [mm]
const double Z_DistCMS = 15000.; // [mm]
const double RadiusTracker =  1200.; // [mm]
const double Z_DistTracker = 2800.; // [mm]
// cylinder actually used in the code
//const double RadiusTarget = RadiusCMS; // [mm]  //now controlled by cfg-file!!!
//const double Z_DistTarget = Z_DistCMS; // [mm]  //now controlled by cfg-file!!!



//define different materials
enum {Unknown=0, Plug, Wall, Air, Clay, Rock};


//Parameters for upward muons from neutrinos
const double N_A = 6.022e23; //mol^-1
const double alpha = 2.; //MeV/(g/cm^2)
const double beta_const = 3.9e-6; //(g/cm^2)^-1
const double epsilon = alpha/beta_const;
const double Rearth = 6370.e6; //mm

//Multi Muon relevant parameters
const double NorthCMSzDeltaPhi = 3./8.*Pi; //rad (Pi/2 if CMS -x = North)
const int max_Trials = 200000;


#endif