60 likes | 168 Views
Cosmic-ray generator for End-to-End simulation. Tsunefumi Mizuno Apr 09, 2002. Structure of Cosmic-ray generator. CrProtonPrimary CrProtonReentrant CrProtonSplash CrElectronPrimary …. user. CrGenerator. each compoment is derived from CrSpectrum class.
E N D
Cosmic-ray generator for End-to-End simulation Tsunefumi Mizuno Apr 09, 2002
Structure of Cosmic-ray generator CrProtonPrimary CrProtonReentrant CrProtonSplash CrElectronPrimary … user CrGenerator each compoment is derived from CrSpectrum class • CrGenerator: interfaces with cosmic-ray components • main methods and members(p3) • CrSpectrum: base class of cosmic-ray components • calculate COR/phi from position/time(p4 and 5)
methods and members private: std::vector<CrSpectrum*> m_subComponents; CrSpectrum* m_component; double m_energy; std::pair<double,double> m_dir; CrGenerator src("CrProtonPrimary");src.addcomponent("CrProtonReentrant");src.addcomponent("CrProtonSplash"); HepRandomEngine* engine = new HepJamesRandom; // 31.78deg, -95.73deg (GLAST Balloon Experiment)src.setPosition(0.5547, -1.671, 0.0); src.ShootParticle(engine);G4double energy = src.GetEnergy();std::pair<G4double,G4double> dir = src.GetDir();G4double theta = acos(dir.first); // [0, pi]G4double phi = dir.second; // * rad; will be corrected to meet the naming convention
Calculate COR from geographic coordinate (1) Cutoff rigidity is calculated in CrSpectrum class Geographic latitude/longitude Geomagnetic latitude/longitude http://hpamsmi2.mi.infn.it/~wwwams/geo.html
Calculate COR from geographic coordinate (2) geograhic latitude and longitude of the geomagnetic north pole in 2000http://swdcdb.kugi.kyoto-u.ac.jp/trans/index.html
ToDo List • calculate solar potential from time • interface with XML