flyopunta/src/src_main/fitscntrl.cpp

/***************************************************************************
                          gdevfis.cpp  -  description
                             -------------------
    begin                : Sun Nov 26 2000
    copyright            : (C) 2000 by Giuseppe m Pierazzini
    email                : pierazzini@unipi.it
 ***************************************************************************
 *                                                                         *
 *   NA48  simulation program.                                             *
 *                                                                         *
 ***************************************************************************/


//
//      Fisica    Pisa 4- 11-2000  GmP
//
#define FISICA

#define APRI    0
#define RIEMPI  1
#define CHIUDI  2


#define NEUTRI  1
#define CARICHI 2
#define ALLCHR  3
#include "flyoh.h"
#include "flyoana.h"



//-------------------------------------------------------------
//                F i t s C o n t r o l
//-------------------------------------------------------------
void FitsControl(int caso)
{
  // caso 0 inizializza
  // caso 1 fits
  // caso 2 summary


  Analisi *ana;


  //  int ht,hts;
  int ko=0;

  switch (caso)
  {
  case APRI:
    //int K4g,K6g,K2p,K3p,K3l,Ka2pi,Pi2g,E2g,E6g,K2p2g,Bea;
    //int Kenu,Kpc2g,Kpcp0,Kpcp0g,Kpc2p0,K2p0enu,K2pcenu,Kpc3g;

    if(K4g>0)     Ak4gam  = new AnaK4g;     // tipo 1
    if(K6g>0)     Ak6gam  = new AnaK6g;     // tipo 2
    if(K2p>0)     Ak2pi   = new AnaK2pi;    // tipo 3
    if(K3p>0)     Ak3pi   = new AnaK3pi;    // tipo 4 
    if(Kl2>0)     Akl2    = new AnaKl2;     // tipo 23 
    if(Ke2g>0)    Ake2g  = new AnaKe2g;     // tipo 24
    if(Kpivv>0)   Akpivv  = new AnaKpivv;   // tipo 25 
    if(Kmunu>0)   Akmunu = new AnaKmunu;   // tipo 29
    if(Kpip0>0)   Akpip0  = new AnaKpip0;  // tipo 28
    if(Kpip0ee>0) Akpip0ee  = new AnaKpip0ee;  // tipo 31
    if(Kpovv>0)   Akpovv  = new AnaKpovv;   // tipo 26
    if(Bea>0)     Beamprt = new AnaBeam;    // tipo 30; 
    if(Ausr>0)    Auser   = new AnaUser;    // tipo 0
    if(Klambda>0) Aklambda= new AnaKlambda; // tipo 27
    break;

  case RIEMPI:

    // analizza i dati se richiesto
    ana=Primofit;
    FitDone=-10;
    while (ana != 0)
    {
      //       ana->Write_titol(); Gout<<"\n FisicaContrl    "<<endl;
      
      if(FitDone<(ko=ana->Fisica())) FitDone=ko;
      //     Gout<<"FitFlag "<<FitDone<<endl;

      ana=ana->next;
    }
    // still FitDone=-100  No fit requested
    // FitDone=ko;     // >= -1  bad fit ; 1 fit ok;
    // Gout<<"\n FitCtr "<<evento_.Event<<" FitDone="<<FitDone<<" WrtNt="<<WrtNt;
    break;

  case CHIUDI:

    ana=Primofit;
  while (ana != 0){ana->print_scale(); ana=ana->next; }

    break;
  default:
    break;
  }
}