CellKG.cxx

Go to the documentation of this file.

// Cluster3D package
//
// CellKG.cxx
//
// CellKG - single cell 
//
// Author:  Katarzyna Grzelak K.Grzelak1@physics.ox.ac.uk
//

#include "Cluster3D/CellKG.h"
#include "MessageService/MsgService.h" // needed for CVS
#include "Cluster3D/StripKG.h"
#include "Cluster3D/StripListKG.h"
#include <math.h>

ClassImp(CellKG)

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

CellKG::CellKG():index(0),u(0.),v(0.),z(0.),plane(0),planeview(0){
//
//  Purpose:  Default constructor
//

  MSG("Cluster3D", Msg::kVerbose)
    << "Default constructor " << endl;
}

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

Bool_t CellKG::Create(UInt_t index1, UInt_t index2, StripListKG &all1, Int_t whichOne, vector<CellKG> &allcells, Int_t &cellnb){
//
//  Purpose:  Create a cell
//

  Bool_t unique = true;
  
  if(all1.GetStripView(index1)==2){
    u=all1.GetStripTPos(index1);  
    v=all1.GetStripTPos(index2);  
  }
  else if(all1.GetStripView(index1)==3){
    u=all1.GetStripTPos(index2);    
    v=all1.GetStripTPos(index1);  
  }

  elem1=index1;
  elem2=index2;
  
  if(whichOne==1){
    z=all1.GetStripZPos(index1);  
    plane=all1.GetStripPlane(index1);  
    if(all1.GetStripView(index1)==2){
      planeview=2;
    }
    else if(all1.GetStripView(index1)==3){
      planeview=3;
    }
  }
  else if(whichOne==2){
    z=all1.GetStripZPos(index2);  
    plane=all1.GetStripPlane(index2);  

    if(all1.GetStripView(index2)==2){
      planeview=2;
    }
    else if(all1.GetStripView(index2)==3){
      planeview=3;
    }
  }

   for(UInt_t ii=0; ii<allcells.size(); ++ii){

     if( (( fabs(u - allcells[ii].u)<0.000001)
          &&  (fabs(v - allcells[ii].v)<0.000001)
         &&  (fabs(z - allcells[ii].z)<0.000001)
        ) ){
       unique = false;
       break;
     } 
     else{
       unique = true;
     }
   }

   if(unique) {
      index=cellnb;
      ++cellnb;
   }

   return unique;

}

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

Bool_t CellKG::Create(UInt_t index, StripListKG &all1, vector<CellKG> &allcells, Int_t &cellnb){
//
//  Purpose:  Create a long cell (from one strip)
//
   Bool_t unique = true;

  if(all1.GetStripView(index)==2){
    u=all1.GetStripTPos(index);    
    v=9999.;
    planeview=2;
  }
  else if(all1.GetStripView(index)==3){
    u=9999.;  
    v=all1.GetStripTPos(index);  
    planeview=3;
  }

  elem1=index;
  elem2=index;
  
  z=all1.GetStripZPos(index);  
  plane=all1.GetStripPlane(index);  

   for(UInt_t ii=0; ii<allcells.size(); ++ii){

     if( ( planeview==2 && fabs(u - allcells[ii].u)<0.000001)
         || (planeview==3 && fabs(v - allcells[ii].v)<0.000001 )
         &&  (fabs(z - allcells[ii].z)<0.000001) ){
       unique = false;
       break;
     } 
     else{
       unique = true;
     }
   }

   if(unique) {
      index=cellnb;
      ++cellnb;
   }

   return unique;

}

---