//-----------------------------------------------------------------------------
//  TopSummaryBlock
//  Author:    Phillip Koehn
//  Date:      Sept. 10 2001
// 
//-----------------------------------------------------------------------------
#include <string>
#include <iostream>
#include <iomanip>

#include <Stntuple/obj/TStnTopSummaryBlock.hh>

ClassImp(TStnTopSummaryBlock)
//______________________________________________________________________________
void TStnTopSummaryBlock::Streamer(TBuffer &R__b)
{
   // Stream an object of class TStnTopSummaryBlock.

   if (R__b.IsReading()) {
     Version_t R__v = R__b.ReadVersion(); 

       R__b >> fTopEventClass;
       R__b >> fnTightEle;
       R__b >> fnTightMuo;
       R__b >> fnTightJet;
       R__b >> fnTightLepton;
       R__b >> fnLooseEle;
       R__b >> fnLooseMuo;
       R__b >> fnLooseJet;
       R__b >> fnLooseLepton;
       R__b >> frawMissingEt;
       R__b >> frawMissingEtPhi;
       R__b >> fvertexCorrMissingEt;
       R__b >> fvertexCorrMissingEtPhi;
       R__b >> fmuonCorrMissingEt;
       R__b >> fmuonCorrMissingEtPhi;
       R__b >> fleptonZVertex;
       R__b >> fjetZVertex;

       //new for version 2
       if(R__v == 2){
	 
	 R__b >> fhadSumEt;
	 R__b >> fhadSumEt3;
	 R__b >> fhadCentrality;
	 R__b >> fhadAplanarity;
	 R__b >> fhadSphericity;
       }

   } else { 
     R__b.WriteVersion(TStnTopSummaryBlock::IsA());
     R__b << fTopEventClass;
     R__b << fnTightEle;
     R__b << fnTightMuo;
     R__b << fnTightJet;
     R__b << fnTightLepton;
     R__b << fnLooseEle;
     R__b << fnLooseMuo;
     R__b << fnLooseJet;
     R__b << fnLooseLepton;
     R__b << frawMissingEt;
     R__b << frawMissingEtPhi;
     R__b << fvertexCorrMissingEt;
     R__b << fvertexCorrMissingEtPhi;
     R__b << fmuonCorrMissingEt;
     R__b << fmuonCorrMissingEtPhi;
     R__b << fleptonZVertex;
     R__b << fjetZVertex;

     //new for version 2
     R__b << fhadSumEt;
     R__b << fhadSumEt3;
     R__b << fhadCentrality;
     R__b << fhadAplanarity;
     R__b << fhadSphericity;


     
   }
}

//_____________________________________________________________________________
TStnTopSummaryBlock::TStnTopSummaryBlock() {

   fTopEventClass = 0;
   fnTightEle = 0;
   fnTightMuo = 0;
   fnTightJet = 0;
   fnTightLepton = 0;

   fnLooseEle = 0;
   fnLooseMuo = 0;
   fnLooseJet = 0;
   fnLooseLepton = 0;

   frawMissingEt= 0.0;
   frawMissingEtPhi = 0.0;
   fvertexCorrMissingEt= 0.0;
   fvertexCorrMissingEtPhi = 0.0;
   fmuonCorrMissingEt= 0.0;
   fmuonCorrMissingEtPhi = 0.0;

   fleptonZVertex = 0.0;
   fjetZVertex = 0.0;

   //new for version 2
   fhadSumEt = 0.0;
   fhadSumEt3 = 0.0;
   fhadCentrality = 0.0;
   fhadAplanarity = 0.0;
   fhadSphericity = 0.0;


   for(int i = 0; i < 5; i++){
     fTightEleIndices[i] = -1;
     fLooseEleIndices[i] = -1;
     fTightMuoIndices[i] = -1;
     fLooseMuoIndices[i] = -1;
     fTightJetIndices[i] = -1;
     fLooseJetIndices[i] = -1;
   }

}

//_____________________________________________________________________________
TStnTopSummaryBlock::~TStnTopSummaryBlock() {
}


//_____________________________________________________________________________
void TStnTopSummaryBlock::Clear(Option_t* opt) {

   fTopEventClass = 0;
   fnTightEle = 0;
   fnTightMuo = 0;
   fnTightJet = 0;
   fnTightLepton = 0;

   fnLooseEle = 0;
   fnLooseMuo = 0;
   fnLooseJet = 0;
   fnLooseLepton = 0;

   frawMissingEt= 0.0;
   frawMissingEtPhi = 0.0;
   fvertexCorrMissingEt= 0.0;
   fvertexCorrMissingEtPhi = 0.0;
   fmuonCorrMissingEt= 0.0;
   fmuonCorrMissingEtPhi = 0.0;

   fleptonZVertex = 0.0;
   fjetZVertex = 0.0;

   //new for version 2
   fhadSumEt = 0.0;
   fhadSumEt3 = 0.0;
   fhadCentrality = 0.0;
   fhadAplanarity = 0.0;
   fhadSphericity = 0.0;


   for(int i = 0; i < 5; i++){
     fTightEleIndices[i] = -1;
     fLooseEleIndices[i] = -1;
     fTightMuoIndices[i] = -1;
     fLooseMuoIndices[i] = -1;
     fTightJetIndices[i] = -1;
     fLooseJetIndices[i] = -1;
   }
}

//_____________________________________________________________________________
void TStnTopSummaryBlock::Print(Option_t* opt) const {
  std::cout << " ************** print_TopSummary ************ " << std::endl << std::endl;
  std::cout.precision(5);
  std::cout.setf(std::ios::right,std::ios::adjustfield);
  std::cout << "" << std::endl;
  std::cout << " TopEventClass  = " << fTopEventClass  << std::endl;
  std::cout << " -------------------------------------------- " << std::endl;

  std::cout << " nTightEle = " << fnTightEle << " nLooseEle = " << fnLooseEle << std::endl;
  std::cout << " nTightMuo = " << fnTightMuo << " nLooseMuo = " << fnLooseMuo << std::endl;
  std::cout << " nTightJet = " << fnTightJet << " nLooseJet = " << fnLooseJet << std::endl;
  std::cout << " nTightLep = " << fnTightLepton << " nLooseLep = " << fnLooseLepton << std::endl;

  std::cout << " Raw Met        = " << frawMissingEt        << " Phi = " << frawMissingEtPhi <<  std::endl;
  std::cout << " VertexCorr Met = " << fvertexCorrMissingEt << " Phi = " << fvertexCorrMissingEtPhi <<  std::endl;
  std::cout << " MuonCorr Met   = " << fmuonCorrMissingEt   << " Phi = " << fmuonCorrMissingEtPhi <<  std::endl;

  std::cout << " LeptonZVertex  = " << fleptonZVertex  << std::endl;
  std::cout << " JetZVertex  = " << fjetZVertex  << std::endl;

  //new for version 2
  std::cout << " Had Sum Et  = " << fhadSumEt  << std::endl;
  std::cout << " Had Sum Et3  = " << fhadSumEt3  << std::endl;
  std::cout << " Had Centrality = " << fhadCentrality  << std::endl;
  std::cout << " Had Aplanarity  = " << fhadAplanarity  << std::endl;
  std::cout << " Had Sphericity  = " << fhadSphericity  << std::endl;


  /*  std::cout << " -----------Not filled anymore --------------------------------- " << std::endl;

  for(int i = 0; i < fnTightEle; i++){
    std::cout << "   fTightEleIndices["<< i << "] = " << fTightEleIndices[i] << std::endl;
  }  
  for(int i = 0; i < fnLooseEle; i++){
    std::cout << "   fLooseEleIndices["<< i << "] = " << fLooseEleIndices[i] << std::endl;
  }
  for(int i = 0; i < fnTightMuo; i++){
    std::cout << "   fTightMuoIndices["<< i << "] = " << fTightMuoIndices[i] << std::endl;
  }  
  for(int i = 0; i < fnLooseMuo; i++){
    std::cout << "   fLooseMuoIndices["<< i << "] = " << fLooseMuoIndices[i] << std::endl;
  }  
  for(int i = 0; i < fnTightJet; i++){
    std::cout << "   fTightJetIndices["<< i << "] = " << fTightJetIndices[i] << std::endl;
  }  
  for(int i = 0; i < fnLooseJet; i++){
    std::cout << "   fLooseJetIndices[i] = " << fLooseJetIndices[i] << std::endl;
  }
  */
  
  
  std::cout << "" << std::endl;
  std::cout << " ************ end print_TopSummary ********** " << std::endl << std::endl;
}