#include "Stntuple/run1/StnRun1InitDataBlocks.hh"
#include "Stntuple/run1/TStnRun1Event.hh"
#include "Stntuple/obj/TStnPhoton.hh"
#include "Stntuple/obj/TStnPhotonBlock.hh"

//_____________________________________________________________________________
Int_t StnRun1InitPhoton(TStnPhoton* Photon, StnPhotonsBlock_t* R1Phot, Int_t I) 
{
  //  init photon object from a UC standard Ntuple

  Photon->fMomentum.SetXYZT(R1Phot->Pp4[I][0],R1Phot->Pp4[I][1],
			    R1Phot->Pp4[I][2],R1Phot->Pp4[I][3]);

  Photon->fDetector      =     R1Phot->Pdetect[I];     
  Photon->fZv      	 =     R1Phot->Pzv[I];
  Photon->fEt      	 =     R1Phot->Pet[I];
  Photon->fEtc    	 =     R1Phot->Petc[I];
  Photon->fCo4     	 =     R1Phot->Pco4[I];
  Photon->fCesx    	 =     R1Phot->Pcesx[I];
  Photon->fCesz    	 =     R1Phot->Pcesz[I];      
  Photon->fCpr5ph  	 =     R1Phot->Pcpr5ph[I];    
  Photon->fCpr5ps  	 =     R1Phot->Pcpr5ps[I];    
  Photon->fDteta   	 =     R1Phot->Pdteta[I];     
  Photon->fEveta   	 =     R1Phot->Peveta[I];     
  Photon->fCese    	 =     R1Phot->Pcese[I];      
  Photon->fPhi     	 =     R1Phot->Pphi[I];       
  Photon->fCo7     	 =     R1Phot->Pco7[I];       
  Photon->fHadem   	 =     R1Phot->Phadem[I];     
  Photon->fLshr    	 =     R1Phot->Plshr[I];      
  Photon->fStr2    	 =     R1Phot->Pstr2[I];      
  Photon->fStre2   	 =     R1Phot->Pstre2[I];     
  Photon->fWwir2   	 =     R1Phot->Pwir2[I];      
  Photon->fWire2   	 =     R1Phot->Pwire2[I];     
  Photon->fSumpt4  	 =     R1Phot->Psumpt4[I];    
  Photon->fVtxocc  	 =     R1Phot->Pvtxocc[I];    
  Photon->fPt      	 =     R1Phot->Ppt[I];        
  Photon->fN3d     	 =     R1Phot->Pn3d[I];       
  Photon->fChi     	 =     R1Phot->Pchi[I];       
  Photon->fChi3x3  	 =     R1Phot->Pchi3x3[I];    
  Photon->fChieta  	 =     R1Phot->Pchieta[I];    
  Photon->fChiphi  	 =     R1Phot->Pchiphi[I];    
  Photon->fStat    	 =     R1Phot->Pstat[I];      
  Photon->fBoxiso  	 =     R1Phot->pboxiso[I];    
  Photon->fCesflg 	 =     R1Phot->pcesflg[I];    
  Photon->fE       	 =     R1Phot->pe[I];         
  Photon->fCeswht  	 =     R1Phot->pceswht[I];    
  Photon->fCprwht  	 =     R1Phot->pcprwht[I];    
  Photon->fCesslide	 =     R1Phot->pcesslide[I];  
  Photon->fWrd     	 =     R1Phot->pwrd[I];       
  Photon->fSth     	 =     R1Phot->psth[I];       
  Photon->fElind  	 =     R1Phot->pelind[I];     
  Photon->fChistr  	 =     R1Phot->pchistr[I];    
  Photon->fChiwir  	 =     R1Phot->pchiwir[I];    
  Photon->fCessig  	 =     R1Phot->pcessig[I];    
  Photon->fCespg   	 =     R1Phot->pcespg[I];     
  Photon->fCespb   	 =     R1Phot->pcespb[I];     
  Photon->fCprpg   	 =     R1Phot->pcprpg[I];     
  Photon->fCprpb   	 =     R1Phot->pcprpb[I];     
  Photon->fPt2     	 =     R1Phot->ppt2[I];       
  Photon->fNtrwd3  	 =     R1Phot->pntrwd3[I];    
  Photon->fCescprx 	 =     R1Phot->pcescprx[I];   
  Photon->fTkcprx  	 =     R1Phot->ptkcprx[I];    
  Photon->fTime    	 =     R1Phot->ptime[I];
  Photon->fVcor   	 =     R1Phot->pvcor[I];
  Photon->fLcor    	 =     R1Phot->plcor[I];

  return 0;
}		      


//_____________________________________________________________________________
int StnRun1InitPhotonBlock(TStnDataBlock* Block, TStnEvent* Event, int Mode) 
{
  // initialize photon data block starting from Run1 UC Stntuple

  TStnRun1Event* event = (TStnRun1Event*) Event;

  TStnPhotonBlock* blk = (TStnPhotonBlock*) Block;
  blk->Clear();

  StnPhotonsBlock_t* r1pho    = event->PhotonsBlock();

  blk->fNcencl     = r1pho->Ncencl;
  blk->fNphoidl    = r1pho->Nphoidl;
  blk->fNphoisol   = r1pho->Nphoisol;
  blk->fNphoid     = r1pho->Nphoid;
  blk->fNphoiso    = r1pho->Nphoiso;

  for (int i=0; i<r1pho->Npho; i++) {
    TStnPhoton* photon = blk->NewPhoton();
    StnRun1InitPhoton(photon,r1pho,i);
  }
  return 0;
}