Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/extended/electromagnetic/TestEm2/egs4/GammaFit.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/TestEm2/egs4/GammaFit.C (Version 11.3.0) and /examples/extended/electromagnetic/TestEm2/egs4/GammaFit.C (Version 9.1.p1)


  1 //                                                  1 
  2 // Fit of longitudinal profile with Gamma func    
  3 // See Review of Particles Physics - Electroma    
  4 //                                                
  5 //....oooOO0OOooo........oooOO0OOooo........oo    
  6                                                   
  7 Double_t GammaFunction(Double_t* t, Double_t*     
  8   if (ROOT::TMath::Gamma(params[0]) != 0 ) {      
  9    return  (100 * pow(params[1],params[0]) / R    
 10            * pow(t[0],(params[0]-1)) * exp (-p    
 11   } else {                                        
 12    return 0;                                      
 13   }                                               
 14 }                                                 
 15                                                   
 16 //....oooOO0OOooo........oooOO0OOooo........oo    
 17                                                   
 18 void GammaFit() {                                 
 19  cout                                             
 20   << "Try fit [100*pow(b,a)/ROOT::TMath::Gamma    
 21   << endl;                                        
 22 TFile* f = new  TFile("93ref0.root");             
 23 TH1D* h1 = (TH1D*) f->Get("4");                   
 24 //put error in h1                                 
 25 TH1D* h2 = (TH1D*) f->Get("5");                   
 26 int nmax = h1->GetNbinsX();                       
 27 for (int n=0; n<nmax; n++) {                      
 28   Double_t er = h2->GetBinContent(n);             
 29   h1->SetBinError(n,er);                          
 30 }                                                 
 31                                                   
 32 TF1* func = new TF1("fit", GammaFunction,-0.5,    
 33 func->SetParameter(0,1.);                         
 34 func->SetParameter(1,1.);                         
 35 func->SetParNames("a", "b");                      
 36                                                   
 37 h1->Fit("fit","r","");                            
 38 h1->SetMarkerColor(kRed);                         
 39 h1->SetMarkerStyle(3);                            
 40 h1->Draw("epl");                                  
 41                                                   
 42 Double_t a = func->GetParameter(0);               
 43 Double_t b = func->GetParameter(1);               
 44 Double_t tmax = (a-1)/b;                          
 45                                                   
 46 // Chi-square distribution                        
 47 double chisq=func->GetChisquare();                
 48 double ndf=func->GetNDF();                        
 49 double chisqdf=chisq/ndf;                         
 50                                                   
 51 gStyle->SetOptFit(1111);                          
 52                                                   
 53 cout << "-------------------------------------    
 54 cout << "\n Chisquare: " << chisq << " / " <<     
 55 cout << "\n-----------------------------------    
 56 cout << "\n a = " << a;                           
 57 cout << "\n b = " << b;                           
 58 cout << "\n tmax = " << tmax;                     
 59 cout << "\n-----------------------------------    
 60 }                                                 
 61                                                   
 62 //....oooOO0OOooo........oooOO0OOooo........oo    
 63