{\rtf1\ansi\ansicpg1252\deff0\deflang2057{\fonttbl{\f0\fswiss\fcharset0 Arial;}}
{\*\generator Msftedit 5.41.15.1507;}\viewkind4\uc1\pard\f0\fs20 #include <iostream.h>\par
#include <vector.h>\par
#include <TClass.h>\par
#include <Riostream.h>\par
#include <TH1F.h>\par
#include <TChain.h>\par
#include <fstream>\par
\par
void LoadTree02(const char* filename)\{\par
\par
ofstream myfile;\par
\par
myfile.open ("LoadtreeOutput/B0sPhiGamma01.txt");\par
\par
TChain chain("100");\par
 \par
chain.Add(filename);\par
\par
int num=25000;\par
int particles;\par
int  id[num];\par
float pt[num];\par
float pseudorapidity[num];\par
float phi[num];\par
int child_of[num];\par
\par
int countk=0;\par
int countg=0;\par
int countp=0;\par
int counta=0;\par
int countb=0; \par
int count=0;\par
\par
float etaga[num];\par
float etagb[num];\par
float etap[num];\par
float etak[num];\par
float etadiffa[num];\par
float etadiffb[num];\par
\par
float phiga[num];\par
float phigb[num];\par
float phip[num];\par
float phik[num];\par
float phidiffa[num];\par
float phidiffb[num];\par
\par
\par
\par
int nentries = (Int_t) chain.GetEntries();\par
\par
chain.SetBranchAddress("particles",&particles);\par
chain.SetBranchAddress("id",id);\par
chain.SetBranchAddress("pt",pt);\par
chain.SetBranchAddress("child_of",child_of);\par
chain.SetBranchAddress("pseudorapidity",pseudorapidity);\par
chain.SetBranchAddress("phi",phi);\par
\par
TH1F* eta = new TH1F("h_ptsecond","pt of first particle in chain",100,-3.5,3.5); //creating the histogram i am going to fill  \par
\par
\par
\par
for (int i=0; i<nentries ; i++) \{\par
      chain.GetEntry(i);\par
      \par
      myfile << "\\n Entry: " << i << " number of particles " << particles << endl;\par
       for (int ipar=0; ipar<(particles) ; ++ipar) \{\par
\tab  \par
          if (id[ipar]==313&&id[child_of[ipar]-1]==511)\{\par
\tab    countk=1;\par
\tab    etak[counta]=pseudorapidity[ipar];\par
\tab    phik[counta]=phi[ipar];\par
\tab      \}\par
\tab   if (id[ipar]==333&&id[child_of[ipar]-1]==531)\{\par
\tab    countp=1;\par
\tab    etap[countb]= pseudorapidity[ipar];\par
\tab    phip[countb]=phi[ipar];\par
\tab        \}\par
\tab   if (id[ipar]==22&&id[child_of[ipar]-1]==531)\{\par
\tab    countg=1;\par
\tab    etagb[countb]=pseudorapidity[ipar];\par
\tab    phigb[countb]=phi[ipar];\par
\tab      \}\par
          if (id[ipar]==22&&id[child_of[ipar]-1]==511)\{\par
\tab    countg=1;\par
\tab    etaga[counta]=pseudorapidity[ipar];\par
\tab    phiga[counta]=phi[ipar];\par
\tab      \}\par
\tab   if(child_of[ipar]==0)\{\par
\tab     myfile << "particle " << ipar << " id : " << id[ipar] << "     child of:   " << child_of[ipar]-1 <<"phi of:   " << phi[ipar] << endl;\par
 \tab     \}\par
\tab   else\{\par
\tab    myfile << "particle " << ipar << " id : " << id[ipar] << "     child of:   " << child_of[ipar]-1 <<"        "<<  id[child_of[ipar]-1] <<"phi of:   " << phi[ipar] <<endl;\par
\tab    \}\par
\tab   \par
           \} \par
          if(countg==1&&countk==1)\{\par
\tab     ++counta;\par
\tab    \}\par
          if(countg==1&&countp==1)\{\par
\tab     ++countb;\par
           \}\par
\tab   countk=0;\par
\tab   countg=0;\par
\tab   countp=0;\par
      \}\par
\par
if (counta!=0)\{\par
\tab     count=counta;\par
\tab  \}\par
if (countb!=0) \{\par
    count=countb;\par
      \}\par
\par
for (int s=0; s<count ; s++) \{\par
       \tab etadiffa[s]= etak[s]-etaga[s];\par
       \tab etadiffb[s]= etap[s]-etagb[s];\par
\tab phidiffa[s]= phik[s]-phiga[s];\par
      \tab phidiffb[s]= phip[s]-phigb[s];\par
\tab if(phidiffa[s]>3.14159265)\{\par
\tab   phidiffa[s]=phidiffa[s]-2*3.14159265;\par
\tab     \}\par
\tab if(phidiffb[s]>3.14159265)\{\par
\tab   phidiffb[s]=phidiffb[s]-2*3.14159265;\par
\tab     \}\par
\tab if(phidiffa[s]<-3.14159265)\{\par
\tab   phidiffa[s]=phidiffa[s]+2*3.14159265;\par
\tab     \}\par
\tab if(phidiffb[s]<-3.14159265)\{\par
\tab   phidiffb[s]=phidiffb[s]+2*3.14159265;\par
\tab     \}\par
\tab myfile << etadiffa[s] << "    " << etadiffb[s] <<"    " << phidiffa[s] <<"    " << phidiffb[s] <<  endl;\par
\par
\tab if (count==counta)\{\par
\tab  eta.Fill(etadiffa[s]);\par
      \}\par
\tab if (count==countb)\{\par
\tab   eta.Fill(etadiffb[s]);\par
\tab\}\par
\par
\tab\par
for (int i=0;i<100>SetBinContent(i+1),etadiffb[i]);\par
\par
      myfile << "No of Decays (B0d -> K* Gamma):    " <<  counta <<endl;//a\par
      myfile << "No of Decays (B0s -> Phi Gamma):    " <<  countb <<endl;//b\par
 myfile.close();\par
  count=0;\par
\par
\}\par
}