#include "AADetectorConstruction.hh"
#include "G4Material.hh"
#include "G4Box.hh"
#include "G4Tubs.hh"
#include "G4LogicalVolume.hh"
#include "G4ThreeVector.hh"
#include "G4PVPlacement.hh"


#include <iostream>
using namespace std;
using namespace CLHEP;


G4VPhysicalVolume* AADetectorConstruction::Construct()
{
  // todo: add salt etc.

  cout << "AADetectorConstruction::Construct()" << endl;

  G4Material* SeaWater = new G4Material("SeaWater",
                                        10.0,          
                                        18*g/mole, 
                                        1.04*g/cm3 );

  G4Tubs* can = new G4Tubs( "Can", 0.0, 1000 *m, 
			    2000 *m ,0.*deg, 360. *deg);
  
  can_log  = new G4LogicalVolume( can, SeaWater, "can_log", 0,0,0 );

  can_phys = new G4PVPlacement( 0, G4ThreeVector(), can_log, "can", 0, false, 0 );

  return can_phys;
}