TpcSimGeometry.cxx

Go to the documentation of this file.

#include "TCanvas.h"
#include "TGraph2D.h"
#include "TH3F.h"
#include "TRandom.h"

#include "TpcSimGeometry.h"

ClassImp(TpcSimGeometry)

TpcSimGeometry::TpcSimGeometry()
{
  fX = 0;
  fY = 0;
  fZ = 0;

  fGasSizeX = 300;
  fGasSizeY = 300;
  fGasSizeZ = 300;

  fBoxThickness = 10;


  for(Int_t i = 0; i<fkNScinti; i++){
    fScintiX[i] = 0;
    fScintiY[i] = 0;
    fScintiZ[i] = 0;
    fScintiSizeX[i] = 0;
    fScintiSizeY[i] = 0;
    fScintiSizeZ[i] = 0;
  }

}


Int_t TpcSimGeometry::GetVolume(Double_t x, Double_t y, Double_t z)
{

  if(TMath::Abs(x-fX)<fGasSizeX/2 &&
     TMath::Abs(y-fY)<fGasSizeY/2 && 
     TMath::Abs(z-fZ-fGasSizeZ/2)<fGasSizeZ/2)
    return 1; // Gas Volume

  if(TMath::Abs(x-fX)<fGasSizeX/2 + fBoxThickness &&
     TMath::Abs(y-fY)<fGasSizeY/2 + fBoxThickness && 
     TMath::Abs(z-fZ-fGasSizeZ/2)<fGasSizeZ/2 + fBoxThickness)
    return 2; // Box/Field cage


  return 0;

}


TCanvas* TpcSimGeometry::DrawDetector()
{

  TGraph2D* g[fkNScinti+3];
  for(Int_t i = 0; i<fkNScinti+3; i++) 
    g[i] = new TGraph2D();

  for(Int_t i = 0; i<100000; i++){
    Double_t x = 400*gRandom->Rndm() - 200;
    Double_t y = 400*gRandom->Rndm() - 200;
    Double_t z = 400*gRandom->Rndm() - 50;

    Int_t vol = GetVolume(x,y,z);
    g[vol]->SetPoint(g[vol]->GetN(),x,y,z);
  }

  g[0]->SetMarkerColor(kCyan);
  g[1]->SetMarkerColor(kBlue);
  g[2]->SetMarkerColor(kRed);

  TH3F* h = new TH3F("h",";x;y;z;",100,0,200,100,0,200,100,-50,350);
  TCanvas* c = new TCanvas();
  h->Draw();
  for(Int_t i = 0; i<fkNScinti+3; i++) 
    g[i]->Draw("Psame");


  return c;
}