//____________________________________________________________________
//
// Event writer or the LPSC event format from the Grenoble group
// (SgTop-D3PD plus BDT input variables)
// 
//  
// Author: Oliver Maria Kind <mailto: kind@mail.desy.de>
// Update: $Id: AtlEvtWriterLPSC.cxx,v 1.2 2014/04/22 08:26:03 kind Exp $
// Copyright: 2014 (C) Oliver Maria Kind
//
#ifndef ATLAS_AtlEvtWriterLPSC
#include <AtlEvtWriterLPSC.h>
#endif
#include <iostream>

using namespace std;

#ifndef __CINT__
ClassImp(AtlEvtWriterLPSC);
#endif

//____________________________________________________________________

AtlEvtWriterLPSC::AtlEvtWriterLPSC() {
    //
    // Default constructor
    //
}

//____________________________________________________________________

AtlEvtWriterLPSC::~AtlEvtWriterLPSC() {
    //
    // Default destructor
    //
}

//____________________________________________________________________

void AtlEvtWriterLPSC::BookTree(TTree *tree, AtlEvent *evt) {
    //
    // Book output tree
    //
    AtlEvtWriterD3PDSgTop::BookTree(tree, evt);

    // Book BDT variables
    fTree->Branch("bdt_MET", &fMET, "bdt_MET/F");
    fTree->Branch("bdt_MtW", &fMtW, "bdt_MtW/F");
    fTree->Branch("bdt_Pt_Lep", &fPtLep, "bdt_Pt_Lep/F");
    fTree->Branch("bdt_DeltaPhiTopBJet2_BJet1",
    		  &fDeltaPhiTopBJet2_BJet1,
    		  "bdt_DeltaPhiTopBJet2_BJet1/F");
    fTree->Branch("bdt_DeltaPhiTopBJet1_BJet2",
    		  &fDeltaPhiTopBJet1_BJet2,
    		  "bdt_DeltaPhiTopBJet1_BJet2/F");
    fTree->Branch("bdt_SumEtLep_MET", &fSumEtLep_MET,
    		  "bdt_SumEtLep_MET/F");
    fTree->Branch("bdt_SumEtLep_MET_BJets", &fSumEtLep_MET_BJets,
    		  "bdt_SumEtLep_MET_BJets/F");
    fTree->Branch("bdt_SumPt", &fSumPt, "bdt_SumPt/F");
    fTree->Branch("bdt_DeltaEtaLep_BJet1", &fDeltaEtaLep_BJet1,
    		  "bdt_DeltaEtaLep_BJet1/F");
    fTree->Branch("bdt_DeltaEtaLep_BJet2", &fDeltaEtaLep_BJet2,
    		  "bdt_DeltaEtaLep_BJet2/F");
    fTree->Branch("bdt_DeltaEtaTopBJet2_BJet1",
    		  &fDeltaEtaTopBJet2_BJet1,
    		  "bdt_DeltaEtaTopBJet2_BJet1/F");
    fTree->Branch("bdt_DeltaPhiBJet1_BJet2", &fDeltaPhiBJet1_BJet2,
    		  "bdt_DeltaPhiBJet1_BJet2/F");
    fTree->Branch("bdt_DeltaPhiLep_MET", &fDeltaPhiLep_MET,
    		  "bdt_DeltaPhiLep_MET/F");
    fTree->Branch("bdt_CosLepW_WTopBJet1", &fCosLepW_WTopBJet1,
    		  "bdt_CosLepW_WTopBJet1/F");
    fTree->Branch("bdt_CosLepTopBJet1_TopBJet1CMS",
    		  &fCosLepTopBJet1_TopBJet1CMS,
    		  "bdt_CosLepTopBJet1_TopBJet1CMS/F");
    fTree->Branch("bdt_CosLepTopBJet2_TopBJet2CMS",
    		  &fCosLepTopBJet2_TopBJet2CMS,
    		  "bdt_CosLepTopBJet2_TopBJet2CMS/F");
    fTree->Branch("bdt_Pt_BJetTop", &fPt_BJetTop,
    		  "bdt_Pt_BJetTop/F");
    fTree->Branch("bdt_DeltaEtaNu_BJetNonTop", &fDeltaEtaNu_BJetNonTop,
    		  "bdt_DeltaEtaNu_BJetNonTop/F");
    fTree->Branch("bdt_SumPtBJet1_BJet2", &fSumPtBJet1_BJet2,
    		  "bdt_SumPtBJet1_BJet2/F");
    fTree->Branch("bdt_CosMET_BJet2", &fCosMET_BJet2,
    		  "bdt_CosMET_BJet2/F");
    fTree->Branch("bdt_MLep_BJet2", &fMLep_BJet2,
    		  "bdt_MLep_BJet2/F");
}