ClusterKG.cxx

Go to the documentation of this file.

// Cluster3D package
//
// ClusterKG.cxx
//
// ClusterKG - a single cluster
//
// Author:  Katarzyna Grzelak K.Grzelak1@physics.ox.ac.uk
//

#include "Cluster3D/ClusterKG.h"
#include "Algorithm/AlgConfig.h"
#include "MessageService/MsgService.h" // needed for CVS
#include "Cluster3D/StripKG.h"
#include "Cluster3D/CellKG.h"
#include "Cluster3D/StripListKG.h"
#include "Cluster3D/CellListKG.h"

ClassImp(ClusterKG)

CVSID("$Id: ClusterKG.cxx,v 1.3 2007/03/01 17:44:17 rhatcher Exp $");

ClusterKG::ClusterKG():energyRaw(0.),energyPe(0.),energyGeV(0.),nbOfAmbStrips(0){
//
//  Purpose:  Default constructor
//
  MSG("Cluster3D", Msg::kVerbose)
    << "Default constructor " << endl;
}

//............................................................................

Bool_t ClusterKG::Create(AlgConfig &ac, StripListKG &all1,  CellListKG &all2, const std::vector<Int_t> &row){
//
//  Purpose: Create a single cluster 
//
//  Arguments:
//
//  Return: 
//
//  Contact:   K.Grzelak1@physics.ox.ac.uk
//
// //  Program Notes:-
//  =============
// Finally cluster should consist of strips, not cells  

  Bool_t bigcluster = true;
  stpindexes.clear();

  //
  // Make sure that cluster consists of unique strips only
  //
  for(UInt_t jj=0; jj<row.size(); ++jj){
      Int_t imult1=0;
      Int_t imult2=0;

      Int_t istp1 = all1.GetStripIndex(all2.GetCellElem1(row[jj])); 
      Int_t istp2 = all1.GetStripIndex(all2.GetCellElem2(row[jj]));          

      // "Long" cells have elem1 = elem2

      for(UInt_t ll=0; ll<stpindexes.size(); ++ll){
        Int_t istp = all1.GetStripIndex(stpindexes[ll]); 
        if(istp1 == istp){
          ++imult1;
        }
        if(istp2 == istp){
          ++imult2;
        }
      }
      if(imult1==0){
        Int_t ind = all2.GetCellElem1(row[jj]) ;
        stpindexes.push_back(ind);
        all1.SetStripAmbiguity(ind); 
      }
      // Check if it was not a "long" cell
      if(imult2==0 && istp1 != istp2){
        Int_t ind = all2.GetCellElem2(row[jj]) ;
        stpindexes.push_back(ind);
        all1.SetStripAmbiguity(ind); 
      }          
  } 

  // Set cluster energy
  SetEnergy(all1);
  // Check how many strips belong to other clusters as well
  FindNbOfAmbStrips(all1);

  MSG("Cluster3D", Msg::kDebug)
      << " Cluster3D create "
      << "energyRaw= " << energyRaw 
      << "energyPe= " << energyPe << endl;

  const Int_t minStp = ac.GetInt("minClustStp");
  const Int_t minAmbStp = ac.GetInt("minClustAmbStp");
  const Double_t minClustPe = ac.GetDouble("minClustPe");

  MSG("Cluster3D", Msg::kDebug)
      << " Cluster3D create "
      << " nbOfAmbStrips= " << nbOfAmbStrips
      << "stpindexes.size()= " << stpindexes.size() << endl;

  //
  // How many strips belongs ONLY to this cluster ?
  // If less than minAmbStp, than it is not an interesting  cluster 
  //

  if(((Int_t)stpindexes.size() - nbOfAmbStrips) < minAmbStp){
    bigcluster=false;
  }

  //
  // Is it simply a big or a small cluster ?
  // If total number of strips is less than  minSt
  // or total pulse height is less than minClustPe
  // than again -> it is not an interesting cluster  
  //

  if((Int_t)stpindexes.size() < minStp || energyPe < minClustPe){
    bigcluster=false; 
  } 

  for(UInt_t ll=0; ll<stpindexes.size(); ++ll){
      MSG("Cluster3D", Msg::kDebug)
        << " Cluster3D create "
        << "all1.GetStripXTalk(stpindexes[ll])= " 
        << all1.GetStripXTalk(stpindexes[ll])
        << "all1.GetStripAmbiguity(stpindexes[ll])= " 
        << all1.GetStripAmbiguity(stpindexes[ll])
        << "all1.GetStripStrip(stpindexes[ll])= " 
        << all1.GetStripStrip(stpindexes[ll])
        << "all1.GetStripPlane(stpindexes[ll])= " 
        << all1.GetStripPlane(stpindexes[ll]) << endl;
  }
  return bigcluster;

}

//............................................................................

void ClusterKG::SetEnergy(StripListKG &all1){
//
//  Purpose: Calculate energy of a cluster  
//
//  Arguments:
//
//  Return: 
//
//  Contact:   K.Grzelak1@physics.ox.ac.uk
//
// //  Program Notes:-
//  =============

  energyRaw=0.0;
  energyPe=0.0;
  
  for(UInt_t i=0; i<stpindexes.size(); ++i){
    energyRaw=energyRaw+all1.GetStripRawE(stpindexes[i])
      +all1.GetStripRawW(stpindexes[i]);
    energyPe=energyPe+all1.GetStripPeE(stpindexes[i])
      +all1.GetStripPeW(stpindexes[i]);
  }
}

//............................................................................

Double_t ClusterKG::GetEnergy(Int_t whatKind) const{
//
//  Purpose: Returns energy of a cluster 
//
//  Arguments:
//
//  Return: 
//
//  Contact:   K.Grzelak1@physics.ox.ac.uk
//
// //  Program Notes:-
//  =============
// 

  if(whatKind==0){
    return energyRaw;
  }
  else if(whatKind==1){
    return energyPe;
  }
  else if(whatKind==2){
    return energyGeV;
  }
  else{
    return energyPe;
  }

  return energyPe;
}

//............................................................................

void ClusterKG::FindNbOfAmbStrips(StripListKG &all1){
//
//  Purpose: Find number of strips which belongs not only to this cluster  
//
//  Arguments: 
//
//  Return: 
//
//  Contact:   K.Grzelak1@physics.ox.ac.uk
//
//  Program Notes:-
//  =============

  for(UInt_t i=0; i<stpindexes.size(); ++i){
    if(all1.GetStripAmbiguity(stpindexes[i])>1){
      nbOfAmbStrips=nbOfAmbStrips+1;
    } 
  }
}

//............................................................................

Int_t ClusterKG::GetNbOfAmbStrips() const{
//
//  Purpose: Return number of strips which belongs not only to this cluster 
//
//  Arguments:
//
//  Return: 
//
//  Contact:   K.Grzelak1@physics.ox.ac.uk
//
//  Program Notes:-
//  =============
// 

  return nbOfAmbStrips;
}

---