/* ======================================================================== *\
!
! *
! * This file is part of MARS, the MAGIC Analysis and Reconstruction
! * Software. It is distributed to you in the hope that it can be a useful
! * and timesaving tool in analysing Data of imaging Cerenkov telescopes.
! * It is distributed WITHOUT ANY WARRANTY.
! *
! * Permission to use, copy, modify and distribute this software and its
! * documentation for any purpose is hereby granted without fee,
! * provided that the above copyright notice appear in all copies and
! * that both that copyright notice and this permission notice appear
! * in supporting documentation. It is provided "as is" without express
! * or implied warranty.
! *
!
!
!   Author(s): Keiichi Mase 09/2004 <mailto:mase@mppmu.mpg.de>
!              Markus Meyer 09/2004 <mailto:meyer@astro.uni-wuerzburg.de>
!
!   Copyright: MAGIC Software Development, 2000-2004
!
!
\* ======================================================================== */

//////////////////////////////////////////////////////////////////////////////
//
// MMuonSearchParCalc
//
// Task to calculate the muon parameters
//
//  This class search for muons and store the information into the container 
// of MMuonSearchPar. The actual calculation is done in the class of 
// MMuonSearchPar. Please see the manual.
//
//  In order to use further information of muons such as the width of arcs,
// the arc length along it, the muons size. Use the infomation stored in 
// MMuonCalibPar. The information will be available by using the task of 
// MMuonCalibParCalc.
//
//
//  Input Containers:
//   [MGeomCam]
//   [MHillas]
//   [MCerPhotEvt]
//
//  Output Containers:
//   [MMuonSearchPar]
//
//////////////////////////////////////////////////////////////////////////////
#include "MMuonSearchParCalc.h"

#include <fstream>

#include "MParList.h"
#include "MGeomCam.h"
#include "MCerPhotEvt.h"
#include "MMuonSearchPar.h"
#include "MLog.h"
#include "MLogManip.h"

using namespace std;

ClassImp(MMuonSearchParCalc);

static const TString gsDefName  = "MMuonSearchParCalc";
static const TString gsDefTitle = "Calculate muon parameters";

// -------------------------------------------------------------------------
//
// Default constructor. 
//
MMuonSearchParCalc::MMuonSearchParCalc
(const char *mupar, const char *name, const char *title)
  : fHillas(NULL), fMuonPar(NULL), fCerPhotName("MCerPhotEvt")
{
    fName  = name  ? name  : gsDefName.Data();
    fTitle = title ? title : gsDefTitle.Data();

    fMuparName  = mupar;
    fHillasInput = "MHillas";
}

// -------------------------------------------------------------------------
//
Int_t MMuonSearchParCalc::PreProcess(MParList *pList)
{
    fHillas = (MHillas*)pList->FindObject(fHillasInput, "MHillas");
    if (!fHillas)
    {
        *fLog << err << fHillasInput << " [MHillas] not found... aborting." << endl;
        return kFALSE;
    }

    fCerPhotEvt = (MCerPhotEvt*)pList->FindObject(AddSerialNumber(fCerPhotName), "MCerPhotEvt");
    if (!fCerPhotEvt)
    {
        *fLog << err << fCerPhotName << " [MCerPhotEvt] not found... aborting." << endl;
        return kFALSE;
    }

    fGeomCam = (MGeomCam*)pList->FindObject("MGeomCam");
    if (!fGeomCam)
    {
        *fLog << err << "MGeomCam not found... aborting." << endl;
        return kFALSE;
    }

    fMuonPar = (MMuonSearchPar*)pList->FindCreateObj("MMuonSearchPar", fMuparName);
    if (!fMuonPar)
        return kFALSE;

    return kTRUE;
}

// -------------------------------------------------------------------------
//
Int_t MMuonSearchParCalc::Process()
{

  fMuonPar->Calc(*fGeomCam, *fCerPhotEvt, *fHillas);

  return kTRUE;
}