#include "PndTrkProt.h"
#include "PndTrkBoundaryParStraws2.h"
#include "PndTrkChi2Fits.h"
#include "PndTrkComparisonMCtruth.h"
#include "PndTrkSttConformalFilling.h"
#include "PndTrkLegendreFits.h"
#include "PndTrkCleanup.h"
#include "PndTrkCTFindTrackInXY.h"
#include "PndTrkCTFindTrackInXY2.h"
#include "PndTrkCTGeometryCalculations.h"
#include "PndTrkMergeSort.h"
#include "PndTrkPlotMacros2.h"
#include "PndTrkPrintouts.h"
#include "PndTrkSttAdjacencies.h"
#include "PndTrkSttClusterFinder.h"

#include "PndSttHit.h"
#include "PndSciTHit.h"
#include "PndSttPoint.h"
#include "PndSttTrack.h"
#include "PndSttPoint.h"
#include "PndSttHelixHit.h"
#include "PndSttSingleStraw.h"
#include "PndSttTube.h"
#include "PndSttMapCreator.h"

#include "PndSdsHit.h"
#include "PndSdsMCPoint.h"

#include "PndTrackCand.h"
#include "PndTrackCandHit.h"
#include "PndTrack.h"

#include "FairRootManager.h"
#include "FairRunAna.h"
#include "FairRuntimeDb.h"
#include "FairTrackParP.h"
#include "FairField.h"

#include "TGeoManager.h"
#include "TClonesArray.h"
#include "TGeoVolume.h"
#include "TVector3.h"
#include "TRandom.h"
#include "TH1F.h"
#include "TMath.h"
#include "TCanvas.h"
#include "TGeoTube.h"
#include "TStopwatch.h"

#include <iostream>
#include <cmath>

using namespace std;

// -----   Default constructor   -------------------------------------------
PndTrkProt::PndTrkProt() : FairTask("Tracking") { 
  sprintf(fSttBranch,"STTHit");
}
// -------------------------------------------------------------------------

PndTrkProt::PndTrkProt(Int_t verbose) : FairTask("Tracking") { 
  sprintf(fSttBranch,"STTHit");
}

// -----   Destructor   ----------------------------------------------------
PndTrkProt::~PndTrkProt() {}

//----------------------------------------- begin PndTrkProt::Init

InitStatus PndTrkProt::Init() {

  // // Get RootManager
  FairRootManager* ioman = FairRootManager::Instance();
  if ( ! ioman ) {
    cout << "-E- PndTrkProt::Init: "
	 << "RootManager not instantiated, return!" << endl;
    return kFATAL;
  }
  //  -----   maps of STT tubes
  PndSttMapCreator *mapper = new PndSttMapCreator(fSttParameters);
  fSttTubeArray = mapper->FillTubeArray();
 
  //    get   the MCTrack  array
  fMCTrackArray = (TClonesArray*) ioman->GetObject("MCTrack");
  if ( ! fMCTrackArray)
    cout << "-E- PndTrkProt::Init: No MCTrack array!" << endl;

  // Get input array   these are the MC point of STT
  fSttPointArray = (TClonesArray*) ioman->GetObject("STTPoint");
  if ( ! fSttPointArray ) {
    cout << "-W- PndSttHelixHitProducer::Init: "
	 << "No STTPoint array, return!" << endl;
    return kERROR;
  }

  // Get input array   hit of STT after digi
  fSttHitArray = (TClonesArray*) ioman->GetObject(fSttBranch);
  if ( ! fSttHitArray ) {
    cout << "-W- PndTrkProt::Init: "
	 << "No STTHit array, return!" << endl;
    return kERROR;
  }

  cout << "-I- PndTrkProt: Initialization successfull" << endl;
  return kSUCCESS;
}

//----------------------------------------- end of PndTrkProt::Init




void PndTrkProt::SetParContainers() {
  FairRuntimeDb* rtdb = FairRunAna::Instance()->GetRuntimeDb();
  fSttParameters = (PndGeoSttPar*) rtdb->getContainer("PndGeoSttPar");
}

// -----   Public method Exec   --------------------------------------------

void PndTrkProt::Exec(Option_t* opt) {
  cout << "Sto entrando in Exec" << endl;
}
//---------------------- end of   PndTrkProt::Exec




ClassImp(PndTrkProt)