//-----------------------------------------------------------------------------
//  Dec 28 2000 P.Murat: Run II STNTUPLE input module
//-----------------------------------------------------------------------------

#include "TROOT.h"
#include "TFile.h"
#include "TChain.h"

#include "Stntuple/obj/TStnNode.hh"

#include "Stntuple/loop/TStnAna.hh"

#include "Stntuple/run1/StnRun1InitDataBlocks.hh"
#include "Stntuple/run1/TStnRun1Event.hh"
#include "Stntuple/run1/TStnRun1InputModule.hh"

ClassImp(TStnRun1InputModule)
//_____________________________________________________________________________
TStnRun1InputModule::TStnRun1InputModule(): TStnInputModule() {
  fTreeName = "/STNTUPLE/h1";
}

//_____________________________________________________________________________
TStnRun1InputModule::TStnRun1InputModule(const char* FileName, 
					 const char* TreeName):
  TStnInputModule("Run1InputModule","Run I Input Module")
{
  // run1 STNTUPLE resides in /STNTUPLE subdirectory and has its name `h1'

  InitChain(FileName,TreeName);
  fCurrent = -1;
  //  SetBranches();
}


//_____________________________________________________________________________
TStnRun1InputModule::TStnRun1InputModule(TChain* Chain):
  TStnInputModule("Run1InputModule","Run I Input Module")
{
  // input module owns chain and deletes it on exit

  fChain = Chain;
  fCurrent = -1;
  //  SetBranches();
}


//_____________________________________________________________________________
TStnRun1InputModule::~TStnRun1InputModule() {
  // destructor: module owns its histograms and a chain

  fHistogramList->Delete();
  delete fHistogramList;
  delete fChain;
}


//_____________________________________________________________________________
TStnRun1Event* TStnRun1InputModule::GetRun1Event() {
  return (TStnRun1Event*) fAna->GetEvent();
}

//_____________________________________________________________________________
Double_t TStnRun1InputModule::GetEntries() {
  return fChain->GetEntries();
}

//_____________________________________________________________________________
Int_t TStnRun1InputModule::LoadEntry(Int_t entry) {
  // deferred input: set the environment to read one entry

  if (!fChain) return -5;
  Int_t centry = fChain->LoadTree(entry);
  if (centry < 0) return centry;
  Int_t itree = fChain->GetTreeNumber();
  if (itree != fCurrent) {
    fCurrent = itree;
    SetBranches();
  }
  return centry;
}

//_____________________________________________________________________________
Int_t TStnRun1InputModule::NextEvent(Int_t Ientry) {
  LoadEntry(Ientry);
  return fChain->GetEntry(Ientry);
}

//_____________________________________________________________________________
TStnNode* TStnRun1InputModule::GetNode(const char* BranchName,
				       const char* ClassName)
{
  char message[100];
				// check if a node with the name exists
  TStnNode* node;

  TStnEvent*  ev = fAna->GetEvent();

  node = ev->FindNode(BranchName);
  if (node) return node;
				// node doesn't exist,
				// check if class exists

  TClass* cl = gROOT->GetClass(ClassName);
  if (! cl) return NULL;
				// make sure branch exists; run 1 STNTUPLE has 
				// fixed number of blocks!, however, a block
				// has several branches corresponding to it

  Int_t (*func)(TStnDataBlock *, TStnEvent *, Int_t);

  if (strcmp(BranchName,"ElectronBlock") == 0) {
    func = StnRun1InitElectronBlock;
  }
  else if (strcmp(BranchName,"GenpBlock") == 0) {
    func = StnRun1InitGenpBlock;
  }
  else if (strcmp(BranchName,"HeaderBlock") == 0) {
    func = StnRun1InitHeaderBlock;
  }
  else if (strcmp(BranchName,"JetBlock") == 0) {
    func = StnRun1InitJetBlock;
  }
  else if (strcmp(BranchName,"MetBlock") == 0) {
    func = StnRun1InitMetBlock;
  }
  else if (strcmp(BranchName,"FatJetBlock") == 0) {
    func = StnRun1InitJetBlock;
  }
  else if (strcmp(BranchName,"MuonBlock") == 0) {
    func = StnRun1InitMuonBlock;
  }
  else if (strcmp(BranchName,"PhotonBlock") == 0) {
    func = StnRun1InitPhotonBlock;
  }
  else if (strcmp(BranchName,"TagBlock") == 0) {
    func = StnRun1InitTagBlock;
  }
  else if (strcmp(BranchName,"TauBlock") == 0) {
    func = StnRun1InitTauBlock;
  }
  else if (strcmp(BranchName,"TrackBlock") == 0) {
    func = StnRun1InitTrackBlock;
  }
  else if (strcmp(BranchName,"VertexBlock") == 0) {
    func = StnRun1InitVertexBlock;
  }
  else {
    func = 0;
  }

  if (func) {
    node = new TStnNode(BranchName,cl,func);
    ev->GetListOfNodes()->Add(node);
  }

  return (TStnNode*) node;
}

//_____________________________________________________________________________
int TStnRun1InputModule::BeginJob() {
  // I'd expect that from now on it the event is owned by the BigLoop

  TStnRun1Event* ev = new TStnRun1Event();
  fAna->SetEvent(ev);
  TStnNode::SetEvent(ev);

  fChain->SetBranchStatus("*",1);
  return 0;
}

//_____________________________________________________________________________
int TStnRun1InputModule::BeginRun() {
  return 0;
}

//_____________________________________________________________________________
int TStnRun1InputModule::Event(Int_t i) {
  return 0;
}

//_____________________________________________________________________________
int TStnRun1InputModule::EndRun() {
  return 0;
}

//_____________________________________________________________________________
int TStnRun1InputModule::EndJob() {
  return 0;
}


//_____________________________________________________________________________
int TStnRun1InputModule::SetBranches() 
{
  // set addresses of the branches for Run I STNTUPLE

  TStnRun1Event* run1ev = GetRun1Event();
  StnGeneralBlock_t* gen_block = run1ev->GeneralBlock();

  fChain->GetBranch("Version" )->SetAddress(&gen_block->Version);
  fChain->GetBranch("Run"     )->SetAddress(&gen_block->Run);
  fChain->GetBranch("Event"   )->SetAddress(&gen_block->Event);
  fChain->GetBranch("Ilum"    )->SetAddress(&gen_block->Ilum);
  fChain->GetBranch("Trigbits")->SetAddress(&gen_block->Trigbits);
  fChain->GetBranch("Goodrun" )->SetAddress(&gen_block->Goodrun);
  fChain->GetBranch("Brcode"  )->SetAddress(&gen_block->Brcode);
  fChain->GetBranch("Goodtrig")->SetAddress(&gen_block->Goodtrig);
  fChain->GetBranch("Trigword")->SetAddress(&gen_block->Trigword);
  fChain->GetBranch("Njet10"  )->SetAddress(&gen_block->Njet10);
  fChain->GetBranch("Njet15"  )->SetAddress(&gen_block->Njet15);
  fChain->GetBranch("Njet20"  )->SetAddress(&gen_block->Njet20);
  fChain->GetBranch("Njet25"  )->SetAddress(&gen_block->Njet25);
  fChain->GetBranch("Nhpte"   )->SetAddress(&gen_block->Nhpte);
  fChain->GetBranch("Nlcent"  )->SetAddress(&gen_block->Nlcent);
  fChain->GetBranch("Nlplug"  )->SetAddress(&gen_block->Nlplug);
  fChain->GetBranch("Nhptm"   )->SetAddress(&gen_block->Nhptm);
  fChain->GetBranch("Ncmio"   )->SetAddress(&gen_block->Ncmio);
  fChain->GetBranch("Ngtau"   )->SetAddress(&gen_block->Ngtau);
  fChain->GetBranch("Nlepton" )->SetAddress(&gen_block->Nlepton);
  fChain->GetBranch("Etout"   )->SetAddress(&gen_block->Etout);

  StnJetsBlock_t* jets_block = run1ev->JetsBlock();

  fChain->GetBranch("Njt")->SetAddress(&jets_block->Njt);
  fChain->GetBranch("Jtet")->SetAddress(jets_block->Jtet);
  fChain->GetBranch("Jtm")->SetAddress(jets_block->Jtm);
  fChain->GetBranch("Jtp4")->SetAddress(jets_block->Jtp4);
  fChain->GetBranch("Jteveta")->SetAddress(jets_block->Jteveta);
  fChain->GetBranch("Jtdteta")->SetAddress(jets_block->Jtdteta);
  fChain->GetBranch("Jtphi")->SetAddress(jets_block->Jtphi);
  fChain->GetBranch("Jtemf")->SetAddress(jets_block->Jtemf);
  fChain->GetBranch("Jtword")->SetAddress(jets_block->Jtword);
  fChain->GetBranch("Jtcorfm")->SetAddress(jets_block->Jtcorfm);
  fChain->GetBranch("Jtcorfd")->SetAddress(jets_block->Jtcorfd);
  fChain->GetBranch("Jtcorfa")->SetAddress(jets_block->Jtcorfa);
  fChain->GetBranch("Jtntr")->SetAddress(jets_block->Jtntr);
  fChain->GetBranch("Jtntrp")->SetAddress(jets_block->Jtntrp);
  fChain->GetBranch("Jtep")->SetAddress(jets_block->Jtep);
  fChain->GetBranch("Jtnjtrk0")->SetAddress(jets_block->Jtnjtrk0);
  fChain->GetBranch("Jttag")->SetAddress(jets_block->Jttag);
  fChain->GetBranch("Jtegrd")->SetAddress(jets_block->Jtegrd);
  fChain->GetBranch("Jvbpb")->SetAddress(jets_block->Jvbpb);
  fChain->GetBranch("Jpbpb")->SetAddress(jets_block->Jpbpb);
  fChain->GetBranch("Scbpb")->SetAddress(jets_block->Scbpb);
  fChain->GetBranch("Scnpb")->SetAddress(jets_block->Scnpb);
  fChain->GetBranch("Jpppb")->SetAddress(jets_block->Jpppb);
  fChain->GetBranch("Jpnpb")->SetAddress(jets_block->Jpnpb);
  fChain->GetBranch("Jtau1")->SetAddress(jets_block->Jtau1);
  fChain->GetBranch("Jtau2")->SetAddress(jets_block->Jtau2);
  fChain->GetBranch("Jtzvtx")->SetAddress(jets_block->Jtzvtx);
  fChain->GetBranch("Jtdzvtx")->SetAddress(jets_block->Jtdzvtx);
  fChain->GetBranch("Jt_ptout")->SetAddress(jets_block->Jt_ptout);
  fChain->GetBranch("Jtvpt")->SetAddress(jets_block->Jtvpt);
  fChain->GetBranch("Jtnovpt")->SetAddress(jets_block->Jtnovpt);

				// ****** fat jets block

  StnFatJetsBlock_t* fat_block = run1ev->FatJetsBlock();

  fChain->GetBranch("Nj1")->SetAddress(&fat_block->Nj1);
  fChain->GetBranch("J1et")->SetAddress(fat_block->J1et);
  fChain->GetBranch("J1m")->SetAddress(fat_block->J1m);
  fChain->GetBranch("J1p4")->SetAddress(fat_block->J1p4);
  fChain->GetBranch("J1eveta")->SetAddress(fat_block->J1eveta);
  fChain->GetBranch("J1dteta")->SetAddress(fat_block->J1dteta);
  fChain->GetBranch("J1phi")->SetAddress(fat_block->J1phi);
  fChain->GetBranch("J1word")->SetAddress(fat_block->J1word);
  fChain->GetBranch("J1jts")->SetAddress(fat_block->J1jts);
  fChain->GetBranch("J1corfm")->SetAddress(fat_block->J1corfm);
  fChain->GetBranch("J1corfd")->SetAddress(fat_block->J1corfd);
  fChain->GetBranch("j1egrd")->SetAddress(fat_block->j1egrd);

				// ****** MET block

  StnMetBlock_t* met_block = run1ev->MetBlock();

  fChain->GetBranch("Met")->SetAddress(met_block->Met);
  fChain->GetBranch("Metphi")->SetAddress(met_block->Metphi);
  fChain->GetBranch("Metjt")->SetAddress(&met_block->Metjt);
  fChain->GetBranch("Metjtphi")->SetAddress(&met_block->Metjtphi);
  fChain->GetBranch("Metj1")->SetAddress(&met_block->Metj1);
  fChain->GetBranch("Metj1phi")->SetAddress(&met_block->Metj1phi);
  fChain->GetBranch("metc")->SetAddress(&met_block->metc);
  fChain->GetBranch("metcphi")->SetAddress(&met_block->metcphi);
  fChain->GetBranch("htc")->SetAddress(&met_block->htc);
  fChain->GetBranch("metcf")->SetAddress(&met_block->metcf);
  fChain->GetBranch("metcfphi")->SetAddress(&met_block->metcfphi);
  fChain->GetBranch("htcf")->SetAddress(&met_block->htcf);
  fChain->GetBranch("Etatrks")->SetAddress(&met_block->Etatrks);
  fChain->GetBranch("Phitrks")->SetAddress(&met_block->Phitrks);
  fChain->GetBranch("Pttrks")->SetAddress(&met_block->Pttrks);
  fChain->GetBranch("Seedpt")->SetAddress(&met_block->Seedpt);
  fChain->GetBranch("Sumet")->SetAddress(&met_block->Sumet);
  fChain->GetBranch("Sumetjet")->SetAddress(&met_block->Sumetjet);
  fChain->GetBranch("Metsig")->SetAddress(&met_block->Metsig);

				// ****** electron block

  StnElectronBlock_t* ele_block = run1ev->ElectronBlock();

  fChain->GetBranch("Nel")->SetAddress(&ele_block->Nel);
  fChain->GetBranch("Ep4")->SetAddress(ele_block->Ep4);
  fChain->GetBranch("Eetcor")->SetAddress(ele_block->Eetcor);
  fChain->GetBranch("Eet")->SetAddress(ele_block->Eet);
  fChain->GetBranch("Echarge")->SetAddress(ele_block->Echarge);
  fChain->GetBranch("Edteta")->SetAddress(ele_block->Edteta);
  fChain->GetBranch("Eeveta")->SetAddress(ele_block->Eeveta);
  fChain->GetBranch("Ephi")->SetAddress(ele_block->Ephi);
  fChain->GetBranch("Ept")->SetAddress(ele_block->Ept);
  fChain->GetBranch("Estat")->SetAddress(ele_block->Estat);
  fChain->GetBranch("Ezv")->SetAddress(ele_block->Ezv);
  fChain->GetBranch("Eep")->SetAddress(ele_block->Eep);
  fChain->GetBranch("Ehadem")->SetAddress(ele_block->Ehadem);
  fChain->GetBranch("Exces")->SetAddress(ele_block->Exces);
  fChain->GetBranch("Edelx")->SetAddress(ele_block->Edelx);
  fChain->GetBranch("Ezces")->SetAddress(ele_block->Ezces);
  fChain->GetBranch("Edelz")->SetAddress(ele_block->Edelz);
  fChain->GetBranch("Echis")->SetAddress(ele_block->Echis);
  fChain->GetBranch("Echiw")->SetAddress(ele_block->Echiw);
  fChain->GetBranch("Eiso")->SetAddress(ele_block->Eiso);
  fChain->GetBranch("Etiso")->SetAddress(ele_block->Etiso);
  fChain->GetBranch("Eidwrd")->SetAddress(ele_block->Eidwrd);
  fChain->GetBranch("Eiswrd")->SetAddress(ele_block->Eiswrd);
  fChain->GetBranch("Econwrd")->SetAddress(ele_block->Econwrd);
  fChain->GetBranch("Eztrk")->SetAddress(ele_block->Eztrk);
  fChain->GetBranch("Efid")->SetAddress(ele_block->Efid);
  fChain->GetBranch("Edet")->SetAddress(ele_block->Edet);
  fChain->GetBranch("Ebte")->SetAddress(ele_block->Ebte);
  fChain->GetBranch("Enasl")->SetAddress(ele_block->Enasl);
  fChain->GetBranch("Enssl")->SetAddress(ele_block->Enssl);
  fChain->GetBranch("eptrk")->SetAddress(ele_block->eptrk);
  fChain->GetBranch("elshr2")->SetAddress(ele_block->elshr2);
  fChain->GetBranch("echi3")->SetAddress(ele_block->echi3);
  fChain->GetBranch("echid")->SetAddress(ele_block->echid);
  fChain->GetBranch("evtxocc")->SetAddress(ele_block->evtxocc);
  fChain->GetBranch("etime")->SetAddress(ele_block->etime);
  fChain->GetBranch("etrind")->SetAddress(ele_block->etrind);

				// ****** muons block

  StnMuonsBlock_t* muon_block = run1ev->MuonsBlock();

  fChain->GetBranch("Nmu")->SetAddress(&muon_block->Nmu);
  fChain->GetBranch("Murawpt")->SetAddress(muon_block->Murawpt);
  fChain->GetBranch("Mup4")->SetAddress(muon_block->Mup4);
  fChain->GetBranch("Mupt")->SetAddress(muon_block->Mupt);
  fChain->GetBranch("Muq")->SetAddress(muon_block->Muq);
  fChain->GetBranch("Muhad")->SetAddress(muon_block->Muhad);
  fChain->GetBranch("Muem")->SetAddress(muon_block->Muem);
  fChain->GetBranch("Mueveta")->SetAddress(muon_block->Mueveta);
  fChain->GetBranch("Muphi")->SetAddress(muon_block->Muphi);
  fChain->GetBranch("Mustat")->SetAddress(muon_block->Mustat);
  fChain->GetBranch("Muzv")->SetAddress(muon_block->Muzv);
  fChain->GetBranch("Mud0")->SetAddress(muon_block->Mud0);
  fChain->GetBranch("Mudelz")->SetAddress(muon_block->Mudelz);
  fChain->GetBranch("Mudelx")->SetAddress(muon_block->Mudelx);
  fChain->GetBranch("Muiso")->SetAddress(muon_block->Muiso);
  fChain->GetBranch("Mutiso")->SetAddress(muon_block->Mutiso);
  fChain->GetBranch("Mucray")->SetAddress(muon_block->Mucray);
  fChain->GetBranch("Muidwrd")->SetAddress(muon_block->Muidwrd);
  fChain->GetBranch("Muiswrd")->SetAddress(muon_block->Muiswrd);
  fChain->GetBranch("Muztrk")->SetAddress(muon_block->Muztrk);
  fChain->GetBranch("Mufid")->SetAddress(muon_block->Mufid);
  fChain->GetBranch("Mudet")->SetAddress(muon_block->Mudet);
  fChain->GetBranch("Mubte")->SetAddress(muon_block->Mubte);
  fChain->GetBranch("Munasl")->SetAddress(muon_block->Munasl);
  fChain->GetBranch("Munssl")->SetAddress(muon_block->Munssl);
  fChain->GetBranch("muchz0")->SetAddress(muon_block->muchz0);
  fChain->GetBranch("muchd0")->SetAddress(muon_block->muchd0);
  fChain->GetBranch("mutime")->SetAddress(muon_block->mutime);
  fChain->GetBranch("muvtxocc")->SetAddress(muon_block->muvtxocc);
  fChain->GetBranch("mudteta")->SetAddress(muon_block->mudteta);
  fChain->GetBranch("mutrind")->SetAddress(muon_block->mutrind);

				// ****** taus block

  StnTausBlock_t* tau_block = run1ev->TausBlock();

  fChain->GetBranch("Ntau")->SetAddress(&tau_block->Ntau);
  fChain->GetBranch("Tp4")->SetAddress(tau_block->Tp4);
  fChain->GetBranch("Tet")->SetAddress(tau_block->Tet);
  fChain->GetBranch("Tq")->SetAddress(tau_block->Tq);
  fChain->GetBranch("Tzv")->SetAddress(tau_block->Tzv);
  fChain->GetBranch("Tdteta")->SetAddress(tau_block->Tdteta);
  fChain->GetBranch("Teveta")->SetAddress(tau_block->Teveta);
  fChain->GetBranch("Tphi")->SetAddress(tau_block->Tphi);
  fChain->GetBranch("Tpt")->SetAddress(tau_block->Tpt);
  fChain->GetBranch("Tptvis")->SetAddress(tau_block->Tptvis);
  fChain->GetBranch("Tstat")->SetAddress(tau_block->Tstat);
  fChain->GetBranch("Tdelr")->SetAddress(tau_block->Tdelr);
  fChain->GetBranch("Temf")->SetAddress(tau_block->Temf);
  fChain->GetBranch("Tmass")->SetAddress(tau_block->Tmass);
  fChain->GetBranch("tpttr")->SetAddress(tau_block->tpttr);
  fChain->GetBranch("tntr10")->SetAddress(tau_block->tntr10);
  fChain->GetBranch("tnpi0")->SetAddress(tau_block->tnpi0);
  fChain->GetBranch("tmustub")->SetAddress(tau_block->tmustub);
  fChain->GetBranch("ttrind")->SetAddress(tau_block->ttrind);

				// ****** dilepton block

  StnDileptonBlock_t* dil_block = run1ev->DileptonBlock();

  fChain->GetBranch("Ndil")->SetAddress(&dil_block->Ndil);
  fChain->GetBranch("Tzflag")->SetAddress(&dil_block->Tzflag);
  fChain->GetBranch("Dilwrd")->SetAddress(&dil_block->Dilwrd);
  fChain->GetBranch("Metdilep")->SetAddress(&dil_block->Metdilep);
  fChain->GetBranch("Njtdilep")->SetAddress(&dil_block->Njtdilep);
  fChain->GetBranch("Dphiminj")->SetAddress(&dil_block->Dphiminj);
  fChain->GetBranch("Dphiminl")->SetAddress(&dil_block->Dphiminl);
  fChain->GetBranch("Dilmode")->SetAddress(&dil_block->Dilmode);
  fChain->GetBranch("Diltype1")->SetAddress(&dil_block->Diltype1);

				// ****** mass block

  StnMassBlock_t* mass_block = run1ev->MassBlock();

  fChain->GetBranch("Eemass")->SetAddress(&mass_block->Eemass);
  fChain->GetBranch("Mumumass")->SetAddress(&mass_block->Mumumass);
  fChain->GetBranch("Taumass")->SetAddress(&mass_block->Taumass);
  fChain->GetBranch("Phomass")->SetAddress(&mass_block->Phomass);

				// ****** photon block

  StnPhotonsBlock_t* pho_block = run1ev->PhotonsBlock();

  fChain->GetBranch("Npho")->SetAddress(&pho_block->Npho);
  fChain->GetBranch("Ncencl")->SetAddress(&pho_block->Ncencl);
  fChain->GetBranch("Nplugcl")->SetAddress(&pho_block->Nplugcl);
  fChain->GetBranch("Nphoidl")->SetAddress(&pho_block->Nphoidl);
  fChain->GetBranch("Nphoisol")->SetAddress(&pho_block->Nphoisol);
  fChain->GetBranch("Nphoid")->SetAddress(&pho_block->Nphoid);
  fChain->GetBranch("Nphoiso")->SetAddress(&pho_block->Nphoiso);
  fChain->GetBranch("Ncmus")->SetAddress(&pho_block->Ncmus);
  fChain->GetBranch("Costub")->SetAddress(&pho_block->Costub);
  fChain->GetBranch("Pdetect")->SetAddress(pho_block->Pdetect);
  fChain->GetBranch("Pzv")->SetAddress(pho_block->Pzv);
  fChain->GetBranch("Pet")->SetAddress(pho_block->Pet);
  fChain->GetBranch("Petc")->SetAddress(pho_block->Petc);
  fChain->GetBranch("Pp4")->SetAddress(pho_block->Pp4);
  fChain->GetBranch("Pco4")->SetAddress(pho_block->Pco4);
  fChain->GetBranch("Pcesx")->SetAddress(pho_block->Pcesx);
  fChain->GetBranch("Pcesz")->SetAddress(pho_block->Pcesz);
  fChain->GetBranch("Pcpr5ph")->SetAddress(pho_block->Pcpr5ph);
  fChain->GetBranch("Pcpr5ps")->SetAddress(pho_block->Pcpr5ps);
  fChain->GetBranch("Pdteta")->SetAddress(pho_block->Pdteta);
  fChain->GetBranch("Peveta")->SetAddress(pho_block->Peveta);
  fChain->GetBranch("Pcese")->SetAddress(pho_block->Pcese);
  fChain->GetBranch("Pphi")->SetAddress(pho_block->Pphi);
  fChain->GetBranch("Pco7")->SetAddress(pho_block->Pco7);
  fChain->GetBranch("Phadem")->SetAddress(pho_block->Phadem);
  fChain->GetBranch("Plshr")->SetAddress(pho_block->Plshr);
  fChain->GetBranch("Pstr2")->SetAddress(pho_block->Pstr2);
  fChain->GetBranch("Pstre2")->SetAddress(pho_block->Pstre2);
  fChain->GetBranch("Pwir2")->SetAddress(pho_block->Pwir2);
  fChain->GetBranch("Pwire2")->SetAddress(pho_block->Pwire2);
  fChain->GetBranch("Psumpt4")->SetAddress(pho_block->Psumpt4);
  fChain->GetBranch("Pvtxocc")->SetAddress(pho_block->Pvtxocc);
  fChain->GetBranch("Ppt")->SetAddress(pho_block->Ppt);
  fChain->GetBranch("Pn3d")->SetAddress(pho_block->Pn3d);
  fChain->GetBranch("Pchi")->SetAddress(pho_block->Pchi);
  fChain->GetBranch("Pchi3x3")->SetAddress(pho_block->Pchi3x3);
  fChain->GetBranch("Pchieta")->SetAddress(pho_block->Pchieta);
  fChain->GetBranch("Pchiphi")->SetAddress(pho_block->Pchiphi);
  fChain->GetBranch("Pstat")->SetAddress(pho_block->Pstat);
  fChain->GetBranch("pboxiso")->SetAddress(pho_block->pboxiso);
  fChain->GetBranch("pcesflg")->SetAddress(pho_block->pcesflg);
  fChain->GetBranch("pe")->SetAddress(pho_block->pe);
  fChain->GetBranch("pceswht")->SetAddress(pho_block->pceswht);
  fChain->GetBranch("pcprwht")->SetAddress(pho_block->pcprwht);
  fChain->GetBranch("pcesslide")->SetAddress(pho_block->pcesslide);
  fChain->GetBranch("pwrd")->SetAddress(pho_block->pwrd);
  fChain->GetBranch("psth")->SetAddress(pho_block->psth);
  fChain->GetBranch("pelind")->SetAddress(pho_block->pelind);
  fChain->GetBranch("pchistr")->SetAddress(pho_block->pchistr);
  fChain->GetBranch("pchiwir")->SetAddress(pho_block->pchiwir);
  fChain->GetBranch("pcessig")->SetAddress(pho_block->pcessig);
  fChain->GetBranch("pcespg")->SetAddress(pho_block->pcespg);
  fChain->GetBranch("pcespb")->SetAddress(pho_block->pcespb);
  fChain->GetBranch("pcprpg")->SetAddress(pho_block->pcprpg);
  fChain->GetBranch("pcprpb")->SetAddress(pho_block->pcprpb);
  fChain->GetBranch("ppt2")->SetAddress(pho_block->ppt2);
  fChain->GetBranch("pntrwd3")->SetAddress(pho_block->pntrwd3);
  fChain->GetBranch("pcescprx")->SetAddress(pho_block->pcescprx);
  fChain->GetBranch("ptkcprx")->SetAddress(pho_block->ptkcprx);
  fChain->GetBranch("ptime")->SetAddress(pho_block->ptime);
  fChain->GetBranch("pvcor")->SetAddress(pho_block->pvcor);
  fChain->GetBranch("plcor")->SetAddress(pho_block->plcor);

				// ****** tagging block

  StnTaggingBlock_t* tag_block = run1ev->TaggingBlock();

  fChain->GetBranch("Ntag")->SetAddress(&tag_block->Ntag);
  fChain->GetBranch("Btgwrd")->SetAddress(&tag_block->Btgwrd);
  fChain->GetBranch("Lxy")->SetAddress(tag_block->Lxy);
  fChain->GetBranch("Lxyerr")->SetAddress(tag_block->Lxyerr);
  fChain->GetBranch("Ctau")->SetAddress(tag_block->Ctau);
  fChain->GetBranch("Tagchi")->SetAddress(tag_block->Tagchi);

				// ****** vertex block

  StnVertexBlock_t* vtx_block = run1ev->VertexBlock();

  fChain->GetBranch("Npvert")->SetAddress(&vtx_block->Npvert);
  fChain->GetBranch("Vxprim")->SetAddress(vtx_block->Vxprim);
  fChain->GetBranch("Vclass")->SetAddress(vtx_block->Vclass);
  fChain->GetBranch("Vz")->SetAddress(vtx_block->Vz);
  fChain->GetBranch("Vq")->SetAddress(vtx_block->Vq);
  fChain->GetBranch("ptnov")->SetAddress(&vtx_block->ptnov);
  fChain->GetBranch("vzev")->SetAddress(&vtx_block->vzev);
  fChain->GetBranch("vclassev")->SetAddress(&vtx_block->vclassev);

				// ****** GENP block

  StnGenpBlock_t* genp_block = run1ev->GenpBlock();

  fChain->GetBranch("Ngenp")->SetAddress(&genp_block->Ngenp);
  fChain->GetBranch("Ngid")->SetAddress(genp_block->Ngid);
  fChain->GetBranch("Ngn")->SetAddress(genp_block->Ngn);
  fChain->GetBranch("Ngnp")->SetAddress(genp_block->Ngnp);
  fChain->GetBranch("Ngpid")->SetAddress(genp_block->Ngpid);
  fChain->GetBranch("Gnp4")->SetAddress(genp_block->Gnp4);

				// ****** track block

  StnTrackBlock_t* trk_block = run1ev->TrackBlock();

  fChain->GetBranch("ntr")->SetAddress(&trk_block->ntr);
  fChain->GetBranch("trvind")->SetAddress(trk_block->trvind);
  fChain->GetBranch("trnsvx")->SetAddress(trk_block->trnsvx);
  fChain->GetBranch("trnasl")->SetAddress(trk_block->trnasl);
  fChain->GetBranch("trnssl")->SetAddress(trk_block->trnssl);
  fChain->GetBranch("trcot")->SetAddress(trk_block->trcot);
  fChain->GetBranch("trcurv")->SetAddress(trk_block->trcurv);
  fChain->GetBranch("trz0")->SetAddress(trk_block->trz0);
  fChain->GetBranch("trd0")->SetAddress(trk_block->trd0);
  fChain->GetBranch("trq")->SetAddress(trk_block->trq);
  fChain->GetBranch("trphi0")->SetAddress(trk_block->trphi0);
  fChain->GetBranch("trp4")->SetAddress(trk_block->trp4);
  fChain->GetBranch("trpt")->SetAddress(trk_block->trpt);
  fChain->GetBranch("treta")->SetAddress(trk_block->treta);
  fChain->GetBranch("trexeta")->SetAddress(trk_block->trexeta);
  fChain->GetBranch("trexphi")->SetAddress(trk_block->trexphi);
  fChain->GetBranch("trbcvind")->SetAddress(trk_block->trbcvind);
  fChain->GetBranch("trbcz0")->SetAddress(trk_block->trbcz0);
  fChain->GetBranch("trbcd0")->SetAddress(trk_block->trbcd0);
  fChain->GetBranch("trbcq")->SetAddress(trk_block->trbcq);
  fChain->GetBranch("trbcp4")->SetAddress(trk_block->trbcp4);
  fChain->GetBranch("trbcpt")->SetAddress(trk_block->trbcpt);
  fChain->GetBranch("trbceta")->SetAddress(trk_block->trbceta);
  fChain->GetBranch("triso4")->SetAddress(trk_block->triso4);

//    fChain->GetBranch("triggers")->SetAddress(&fEvent->MoreTrigBlock()->triggers);

//    fChain->GetBranch("Nm")->SetAddress(&fEvent->MoreMuonBlock()->Nm);
//    fChain->GetBranch("Mpt")->SetAddress(fEvent->MoreMuonBlock()->Mpt);
//    fChain->GetBranch("Mx0")->SetAddress(fEvent->MoreMuonBlock()->Mx0);
//    fChain->GetBranch("My0")->SetAddress(fEvent->MoreMuonBlock()->My0);
//    fChain->GetBranch("Mz0")->SetAddress(fEvent->MoreMuonBlock()->Mz0);
//    fChain->GetBranch("Macsm")->SetAddress(fEvent->MoreMuonBlock()->Macsm);
//    fChain->GetBranch("Mtcsm")->SetAddress(fEvent->MoreMuonBlock()->Mtcsm);
//    fChain->GetBranch("Mtcsp")->SetAddress(fEvent->MoreMuonBlock()->Mtcsp);
//    fChain->GetBranch("Mtcsx")->SetAddress(fEvent->MoreMuonBlock()->Mtcsx);
//    fChain->GetBranch("Maftm")->SetAddress(fEvent->MoreMuonBlock()->Maftm);
//    fChain->GetBranch("Mtftm")->SetAddress(fEvent->MoreMuonBlock()->Mtftm);
//    fChain->GetBranch("Mtftp")->SetAddress(fEvent->MoreMuonBlock()->Mtftp);
//    fChain->GetBranch("Maftx")->SetAddress(fEvent->MoreMuonBlock()->Maftx);
//    fChain->GetBranch("Mtftx")->SetAddress(fEvent->MoreMuonBlock()->Mtftx);
//    fChain->GetBranch("Mzftm")->SetAddress(fEvent->MoreMuonBlock()->Mzftm);
//    fChain->GetBranch("Mxftm")->SetAddress(fEvent->MoreMuonBlock()->Mxftm);
//    fChain->GetBranch("Mxftp")->SetAddress(fEvent->MoreMuonBlock()->Mxftp);
//    fChain->GetBranch("Mzftx")->SetAddress(fEvent->MoreMuonBlock()->Mzftx);
//    fChain->GetBranch("Mxftx")->SetAddress(fEvent->MoreMuonBlock()->Mxftx);
//    fChain->GetBranch("Mfidm")->SetAddress(fEvent->MoreMuonBlock()->Mfidm);
//    fChain->GetBranch("Mfidp")->SetAddress(fEvent->MoreMuonBlock()->Mfidp);
//    fChain->GetBranch("Mfidx")->SetAddress(fEvent->MoreMuonBlock()->Mfidx);
//    fChain->GetBranch("Maphm")->SetAddress(fEvent->MoreMuonBlock()->Maphm);
//    fChain->GetBranch("Mrphm")->SetAddress(fEvent->MoreMuonBlock()->Mrphm);
//    fChain->GetBranch("Mathm")->SetAddress(fEvent->MoreMuonBlock()->Mathm);
//    fChain->GetBranch("Mrthm")->SetAddress(fEvent->MoreMuonBlock()->Mrthm);
//    fChain->GetBranch("Maphp")->SetAddress(fEvent->MoreMuonBlock()->Maphp);
//    fChain->GetBranch("Mrphp")->SetAddress(fEvent->MoreMuonBlock()->Mrphp);
//    fChain->GetBranch("Maphx")->SetAddress(fEvent->MoreMuonBlock()->Maphx);
//    fChain->GetBranch("Mrphx")->SetAddress(fEvent->MoreMuonBlock()->Mrphx);
//    fChain->GetBranch("Mathx")->SetAddress(fEvent->MoreMuonBlock()->Mathx);
//    fChain->GetBranch("Mrthx")->SetAddress(fEvent->MoreMuonBlock()->Mrthx);

  return 0;
}