[Clas12_verystrange] REC_ForwardTagger

[Michael Dugger]

Hi,

I now have the correct way to get the forward tagger information using REC 
banks.

Here is what I now do:

int nFT         = myTree.REC_ForwardTagger_time->size();
int nREC        = myTree.REC_Particle_pid->size();

for (int j=0;j<nREC;j++){
   for (int ftIndex=0;ftIndex<nFT;ftIndex++){
     int recParticleIndex = myTree.REC_ForwardTagger_pindex->at(ftIndex);
     int ftDet            = myTree.REC_ForwardTagger_detector->at(ftIndex);
     if (recParticleIndex == j && ftDet == 10){
       ftTime     = myTree.REC_ForwardTagger_time->at(ftIndex);
       ftEnergy   = myTree.REC_ForwardTagger_energy->at(ftIndex);
       double pxTmp = myTree.REC_Particle_px->at(recParticleIndex);
       double pyTmp = myTree.REC_Particle_py->at(recParticleIndex);
       double pzTmp = myTree.REC_Particle_pz->at(recParticleIndex);
       double pTmp  = sqrt(pow(pxTmp,2)+pow(pyTmp,2)+pow(pzTmp,2));
      ftCx         = pxTmp/pTmp;
      ftCy         = pyTmp/pTmp;
      ftCz         = pzTmp/pTmp;
     }
   }
}

Note: This ends up setting the forward tagger variables of interest only 
for the last ft event. However, the vast majority of the time when there 
is more thant a single ft hit with detector = 10, the other hits are 
duplicates. Of course, this has to be explored more fully.

Sorry for any confusion I may have caused.

Take care,
Michael