TMBTRefs.cpp

Go to the documentation of this file.

//                                                                        //
// File: TMBTRefs.cpp                                                     //
//                                                                        //
// Purpose:  Container of TRefArray's for related TMB physics objects     //
//                                                                        //
// Created: 22-MAY-2002  Serban Protopopescu                              //
//                                                                        //

#include "tmb_tree/TMBTRefs.hpp"

ClassImp(TMBTRefs)

TMBTRefs::TMBTRefs(TRefArray* ems, TRefArray* muons, TRefArray* taus, TRefArray* jets,
                   TRefArray* mcparts)
{
  _ems.Clear();
  if (ems) {
    for (Int_t i=0; i<ems->GetLast()+1; i++) {
      _ems.Add( ems->At(i) );
    }
  }
 
  _muons.Clear();
  if (muons) {
    for (Int_t i=0; i<muons->GetLast()+1; i++) {
      _muons.Add( muons->At(i) );
    }
  }

  _taus.Clear();
  if (taus) {
    for (Int_t i=0; i<taus->GetLast()+1; i++) {
      _taus.Add( taus->At(i) );
    }
  }

  _jets.Clear();
  if (jets) {
    for (Int_t i=0; i<jets->GetLast()+1; i++) {
      _jets.Add( jets->At(i) );
    }
  }

  _mcparts.Clear();
  if (mcparts) {
    for (Int_t i=0; i<mcparts->GetLast()+1; i++) {
      _mcparts.Add( mcparts->At(i) );
    }
  }
}

const TMBMCpart* TMBTRefs::mcmatch(const TMBEMCluster* emcl) const{
  int nparts=NMCparts();
  const TMBMCpart* mcpart=0;
  if(nparts==0) return mcpart;
  int nems=NEmcls();
  if(nems==0) return mcpart;

  for(int ip=0; ip<nparts; ++ip){
    mcpart=getTMBMCpart(ip);
    for(int iem=0; iem<nems && abs(mcpart->pdgid())==11 ; ++iem){
      if(emcl==getEmcl(iem)) return mcpart;
    }
  }
  mcpart=0;
  return mcpart;
}

const TMBMCpart* TMBTRefs::mcmatch(const TMBMuon* muon) const{

  int nparts=NMCparts();
  const TMBMCpart* mcpart=0;
  if(nparts==0) return mcpart;
  int nmuons=NMuons();
  if(nmuons==0) return mcpart;

  for(int ip=0; ip<nparts; ++ip){
    mcpart=getTMBMCpart(ip);
    for(int imu=0; imu<nmuons && abs(mcpart->pdgid())==13 ; ++imu){
      if(muon==getMuon(imu)) return mcpart;
    }
  }
  mcpart=0;
  return mcpart;
}

const TMBMCpart* TMBTRefs::mcmatch(const TMBTau* tau) const{
  int nparts=NMCparts();
  const TMBMCpart* mcpart=0;
  if(nparts==0) return mcpart;
  int ntaus=NTaus();
  if(ntaus==0) return mcpart;

  for(int ip=0; ip<nparts; ++ip){
    mcpart=getTMBMCpart(ip);
    for(int itau=0; itau<ntaus && abs(mcpart->pdgid())==15 ; ++itau){
      if(tau==getTau(itau)) return mcpart;
    }
  }
  mcpart=0;
  return mcpart;
}

bool TMBTRefs::isInThis(TObject* objptr) const{

  int n=NEmcls();
  for(int i=0; i<n; ++i){
    if(objptr == _ems.At(i)) return true;
  }

  n=NMuons();
  for(int i=0; i<n; ++i){
    if(objptr == _muons.At(i)) return true;
  }

  n=NTaus();
  for(int i=0; i<n; ++i){
    if(objptr == _taus.At(i)) return true;
  }

  n=NJets();
  for(int i=0; i<n; ++i){
    if(objptr == _jets.At(i)) return true;
  }

  n=NMuons();
  for(int i=0; i<n; ++i){
    if(objptr == _muons.At(i)) return true;
  }

  n=NMCparts();
  for(int i=0; i<n; ++i){
    if(objptr == _mcparts.At(i)) return true;
  }
  return false;
}
        

---