#include "SolidMaterials.hh"
#include "NemenixMuVeto.hh"
#include "G4Material.hh"
#include "G4MaterialTable.hh"
#include "G4LogicalBorderSurface.hh"
#include "G4OpBoundaryProcess.hh"
#include "G4ThreeVector.hh"
#include "G4VisAttributes.hh"
#include "G4Transform3D.hh"
#include "G4UnionSolid.hh"
#include <fstream>
#include <iostream>
#include "G4Box.hh"
#include "G4SubtractionSolid.hh"
#include "G4Tubs.hh"
#include "G4PVPlacement.hh"
#include "G4UnitsTable.hh"
#include <math.h>
#include "G4LogicalBorderSurface.hh"

#include "G4SystemOfUnits.hh"
#include "G4PhysicalConstants.hh"

//#include "SolidMuVetoSD.hh"
#include "G4SDManager.hh"

NemenixMuVeto::NemenixMuVeto(){
    
}

NemenixMuVeto::~NemenixMuVeto() {
}

G4LogicalVolume *NemenixMuVeto::MuVetoConstruction(G4double xx, G4double yy, G4double zz)
{
    
    G4Material * mate;
    
    // make colours *********************************************************
    G4Colour  white   (1.0, 1.0, 1.0) ;
    G4Colour  white_t (1.0, 1.0, 1.0, .95) ;
    G4Colour  grey    (0.5, 0.5, 0.5) ;
    G4Colour  lgrey   (.75, .75, .75) ;
    G4Colour  red     (1.0, 0.0, 0.0, 0.5) ;
    G4Colour  blue    (0.0, 0.0, 1.0) ;
    G4Colour  blue_t  (0.0, 0.0, 1.0, .75) ;
    G4Colour  cyan    (0.0, 1.0, 1.0, 0.5) ;
    G4Colour  magenta (1.0, 0.0, 1.0) ;
    G4Colour  yellow  (1.0, 1.0, 0.0, 0.5) ;
    G4Colour  lblue   (0.0, 0.0, .75) ;
    G4Colour  black   (0.0, 0.0, 0.0) ;
    G4Colour  green   (0.0, 1.0, 0.0) ;
    G4Colour  lgreen  (0.0, .75, 0.0) ;
    
    G4Box* MuVetoBox = new G4Box("MuVetoBox",0.5*xx,0.5*yy,0.5*zz);
    
    mate = G4Material::GetMaterial("Air");
    logMuVetoMod = new G4LogicalVolume(MuVetoBox, mate,"logDetMod",0,0,0);
    
    return logMuVetoMod;
}