|
version 1.1, 2003/10/16 22:17:04
|
version 1.2, 2004/02/19 12:47:30
|
|
|
|
| | |
| integer ntid | integer ntid |
| integer max_events, fit_index | integer max_events, fit_index |
| integer histid /100/ |
integer histid /200/ |
| character*3 fittyp | character*3 fittyp |
| logical retcode | logical retcode |
| write(*,*) ' in me_fitter' | write(*,*) ' in me_fitter' |
|
|
|
| c if(fittyp(1:3).eq.'The')fit_index = 2 | c if(fittyp(1:3).eq.'The')fit_index = 2 |
| c if(fittyp(1:3).eq.'Del')fit_index = 5 | c if(fittyp(1:3).eq.'Del')fit_index = 5 |
| write(*,*) 'fit_index= ',fit_index | write(*,*) 'fit_index= ',fit_index |
| write(*,*) ' calling me_fit_histos' |
write(*,*) ' calling me_fit_histos histid =',histid |
| call me_fit_histos(histid) ! define histograms | call me_fit_histos(histid) ! define histograms |
| write(*,*) ' calling me_fit_init',ntid, max_events | write(*,*) ' calling me_fit_init',ntid, max_events |
| call me_fit_init(ntid, max_events, retcode, fit_index) ! read in ntuple | call me_fit_init(ntid, max_events, retcode, fit_index) ! read in ntuple |
| if (.not. retcode) return |
if (.not. retcode) then |
| |
write(*,*) ' bad return code = ',retcode |
| |
return |
| |
endif |
| npass = 1 | npass = 1 |
| |
write(*,*) ' npass,ipass',npass,ipass |
| do while (npass.le.ipass) | do while (npass.le.ipass) |
| write(*,*) ' calling me_fit_setup' | write(*,*) ' calling me_fit_setup' |
| call me_fit_setup(fit_index) ! set up minuit for fitting | call me_fit_setup(fit_index) ! set up minuit for fitting |
| npass = npass + 1 | npass = npass + 1 |
| enddo | enddo |
| |
write(*,*) ' call me_fit_fill id = ',histid |
| call me_fit_fill(histid) ! fill histograms and ntuple | call me_fit_fill(histid) ! fill histograms and ntuple |
| write(*,*) ' ' | write(*,*) ' ' |
| write(*,*) ' DONE !!! ' | write(*,*) ' DONE !!! ' |
Return to
fitter.f
CVS log
Up to [HallC] / Optics / source_code