Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/extended/electromagnetic/TestEm5/tramu/tramu.C

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

Diff markup

Differences between /examples/extended/electromagnetic/TestEm5/tramu/tramu.C (Version 11.3.0) and /examples/extended/electromagnetic/TestEm5/tramu/tramu.C (Version 3.0)


  1 {                                                   1 
  2    gROOT->Reset();                                
  3                                                   
  4    // Draw histograms fill by Geant4 TestEm5 s    
  5    TFile f1("./opt3.root");                       
  6    TH1D* h1 = (TH1D*) f1.Get("10");               
  7    h1->SetTitle("1 TeV muon in 3 m iron : kine    
  8    h1->GetXaxis()->SetTitle("Ekine (GeV)");       
  9    h1->GetYaxis()->SetTitle("nb/GeV");            
 10    h1->SetStats(kFALSE);  // Eliminate statist    
 11    h1->SetLineColor(kBlue);                       
 12    h1->Draw("HIST");                              
 13 /*                                                
 14    TFile f2("./local.root");                      
 15    TH1D* h2 = (TH1D*) f2.Get("10");               
 16    h2->SetStats(kFALSE);  // Eliminate statist    
 17    h2->SetLineColor(kRed);                        
 18    h2->Draw("SAME HIST");                         
 19 */                                                
 20 /*                                                
 21 * muon 1 TeV/c in 3 m Iron                        
 22 * Particle Data Group. Physics Letters B 592 (    
 23 * distribution of the muon kinetic energy         
 24 * (from 950 GeV to 1000 GeV by bin of 0.5 GeV     
 25 */                                                
 26                                                   
 27    ifstream in;                                   
 28    in.open("mars14.ascii");                       
 29                                                   
 30    // Create a new histogram with mars14.acsii    
 31    int nb_bins = 100;                             
 32    float x_min = 950.;                            
 33    float x_max = 1000.;                           
 34    TH1F* h1f = new TH1F("h1f","",nb_bins,x_min    
 35                                                   
 36    Float_t x, y;                                  
 37    while (1) {                                    
 38       in >> x >> y ;                              
 39       if (!in.good()) break;                      
 40       h1f->Fill(x,y);                             
 41    }                                              
 42    in.close();                                    
 43                                                   
 44    // Draw histogram fill by mars14.acsii valu    
 45    h1f->SetLineColor(kRed);                       
 46    h1f->Draw("SAME");                             
 47                                                   
 48    // Print the histograms legend                 
 49    TLegend* legend = new TLegend(0.2,0.55,0.45    
 50    legend->AddEntry(h1,"local (Urban90)","l");    
 51    //legend->AddEntry(h2,"local (msc90)","l");    
 52    legend->AddEntry(h1f,"Mars14 simul ","L");     
 53    legend->Draw();                                
 54                                                   
 55                                                   
 56 }                                                 
 57