#include <iostream>
#include <iomanip>

class CalculLumiEvent {

public:

    CalculLumiEvent(){
    
		LF = 0.;                // Lumi Fill
        Phase = 0.;             // Phase 
                                //   INJECTION  = 100,
                                //   ERAMPING   = 101,
                                //   LUMITUNING = 102,
                                //   LUMIRUN    = 103,
                                //   ALLHVON    = 104, 
                                //   ENDLUMIRUN = 105,
                                //   UNKNOWN    = 0
        date = 0;               // Time Stamp
        RunNumber = 0;          // Run Number

        EmeanX = 0.;            // Mean energy in X (Calibrated)
        EmeanY = 0.;            // Mean energy in Y (Calibrated)
   
        // Beam
        Xm = 0.;                // X mean position
        Ym = 0.;                // Y mean position
        SigmaX = 0.;            // X sigma
        SigmaY = 0.;
        // Errors on Beam paramaters
        Err_Xm = 0.;
        Err_Ym = 0.;
        Err_SigmaX = 0.;
        Err_SigmaY = 0.;

        // Acceptance
        IntAccX = 0.;	        // Fitted Gaussian integral (~1/Acc) in X
        IntAccY = 0.;	        // Fitted Gaussian integral (~1/Acc) in Y
        IntAcc = 0.;	        // Unused
		
        Acceptance_X = 0.;      // Acceptance in X
        Acceptance_Y = 0.;
        Acceptance = 0.;

        Err_Acceptance_X = 0.;
        Err_Acceptance_Y = 0.;

        Chi2_Acc_X  = 0.; 
        Chi2_Acc_Y  = 0.; 
		

        LumiEtot_X = 0.;
        LumiEtot_Y = 0.;
        LumiFitMC_X = 0.;
        LumiFitMC_Y = 0.;
        Err_LumiEtot_X = 0.;
        Err_LumiEtot_Y = 0.;
        Err_LumiFitMC_X = 0.;
        Err_LumiFitMC_Y = 0.;

    
        LumiRawX = 0.;
        LumiRawY = 0.;
  
        Chi2_MCFit_X = 0.;
        Chi2_MCFit_Y = 0.;   

  
        A2G_MC_X = 0.;
        A2G_MC_Y = 0.;

        Err_A2G_MC_X = 0.;
        Err_A2G_MC_Y = 0.;
 
        for(Int_t l1 = 0; l1<4; l1++)
            for(Int_t l2 = 0; l2<4; l2++) Currents[l1][l2]=0.;
    
        for(Int_t i6 = 0; i6<13; i6++)
            {
                A2G_X[i6] =0.;
                A2G_Y[i6] =0.;
                Err_A2G_X[i6] =0.;
                Err_A2G_Y[i6] =0.;
                Chi2_A2G_X[i6] =0.;
                Chi2_A2G_Y[i6] =0.;
            }
    };

    ~CalculLumiEvent(){};

    Double_t LF;
    Double_t Phase;
    Int_t date;
    Int_t RunNumber;

    Double_t Currents[4][4];
  	
    Double_t A2G_X[13];
    Double_t A2G_Y[13];

    Double_t Err_A2G_X[13];
    Double_t Err_A2G_Y[13];
	 
    Double_t Chi2_A2G_X[13];
    Double_t Chi2_A2G_Y[13];

    Double_t A2G_MC_X;
    Double_t A2G_MC_Y;

    Double_t Err_A2G_MC_X;
    Double_t Err_A2G_MC_Y;

	
    Double_t Acceptance_X;
    Double_t Acceptance_Y;
    Double_t Acceptance;

    Double_t Err_Acceptance_X;
    Double_t Err_Acceptance_Y;
  
    Double_t IntAccX;	 
    Double_t IntAccY;	 
    Double_t IntAcc;	 

	Double_t Chi2_Acc_X;
    Double_t Chi2_Acc_Y;

	
    Double_t EmeanX;
    Double_t EmeanY;

	 
    Double_t Xm;
    Double_t Ym;
    Double_t SigmaX;
    Double_t SigmaY;

    Double_t Err_Xm;
    Double_t Err_Ym;
    Double_t Err_SigmaX;
    Double_t Err_SigmaY;
	
	
    Double_t LumiEtot_X;
    Double_t LumiEtot_Y;

    Double_t LumiFitMC_X;
    Double_t LumiFitMC_Y;

    Double_t Err_LumiEtot_X;
    Double_t Err_LumiEtot_Y;

    Double_t Err_LumiFitMC_X;
    Double_t Err_LumiFitMC_Y;
	
    Double_t LumiRawX;
    Double_t LumiRawY;
  
	Double_t Chi2_MCFit_X;
	Double_t Chi2_MCFit_Y;
	
	
    void Set(CalculLumiEvent * tempo)
    {
        LF = tempo->LF;      
        Phase = tempo->Phase;   
        date = tempo->date;
        RunNumber = tempo->RunNumber;
    

        A2G_MC_X = tempo->A2G_MC_X;
        A2G_MC_Y = tempo->A2G_MC_Y;

        Err_A2G_MC_X = tempo->Err_A2G_MC_X;
        Err_A2G_MC_Y = tempo->Err_A2G_MC_Y;
    

		
        Xm = tempo->Xm;
        Ym = tempo->Ym;
        SigmaX = tempo->SigmaX;
        SigmaY = tempo->SigmaY;

        Err_Xm = tempo->Err_Xm;
        Err_Ym = tempo->Err_Ym;
        Err_SigmaX = tempo->Err_SigmaX;
        Err_SigmaY = tempo->Err_SigmaY;		
		

        EmeanX = tempo->EmeanX;
        EmeanY = tempo->EmeanY;   


        IntAccY  = tempo->IntAccY;
        IntAccX  = tempo->IntAccX;
        IntAcc  = tempo->IntAcc;
		
        Acceptance_X = tempo->Acceptance_X;
        Acceptance_Y = tempo->Acceptance_Y;
        Acceptance = tempo->Acceptance;

        Err_Acceptance_X =tempo->Err_Acceptance_X;
        Err_Acceptance_Y =tempo->Err_Acceptance_Y;		

        Chi2_Acc_X = tempo->Chi2_Acc_X;
        Chi2_Acc_X = tempo->Chi2_Acc_Y;
		
		
		
        LumiEtot_X = tempo->LumiEtot_X;
        LumiEtot_Y = tempo->LumiEtot_Y;
        LumiFitMC_X = tempo->LumiFitMC_X;
        LumiFitMC_Y = tempo->LumiFitMC_Y;

        Err_LumiEtot_X = tempo->Err_LumiEtot_X;
        Err_LumiEtot_Y = tempo->Err_LumiEtot_Y;
        Err_LumiFitMC_X = tempo->Err_LumiFitMC_X;
        Err_LumiFitMC_Y = tempo->Err_LumiFitMC_Y;

        LumiRawX = tempo->LumiRawX;
        LumiRawY = tempo->LumiRawY;

		Chi2_MCFit_X = tempo->Chi2_MCFit_X;
		Chi2_MCFit_Y = tempo->Chi2_MCFit_Y;
    
        for(Int_t i6 = 0; i6<13; i6++)
            {
                A2G_X[i6] = tempo->A2G_X[i6];
                A2G_Y[i6] = tempo->A2G_Y[i6];
                Err_A2G_X[i6] = tempo->Err_A2G_X[i6];
                Err_A2G_Y[i6] = tempo->Err_A2G_Y[i6];
                Chi2_A2G_X[i6] = tempo->Chi2_A2G_X[i6];
                Chi2_A2G_Y[i6] = tempo->Chi2_A2G_Y[i6];
            }     
  
        for(Int_t l1 = 0; l1<4; l1++)
            for(Int_t l2 = 0; l2<4; l2++) Currents[l1][l2] = tempo->Currents[l1][l2];
  
  
    };
};