diff --git a/CODE/construction.cc b/CODE/construction.cc
deleted file mode 100644
index 14881326335b14176ac50e59256aa3e134dce8d3..0000000000000000000000000000000000000000
--- a/CODE/construction.cc
+++ /dev/null
@@ -1,144 +0,0 @@
-#include "construction.hh"
-
-MyDetectorConstruction::MyDetectorConstruction()
-{
-//************************//
-// Number of layer //
- nCols = 10;
- nRows = 10;
-
- fMessenger = new G4GenericMessenger(this, "/detector/", "Detector Construction");
-
- fMessenger->DeclareProperty("nCols", nCols, "Number of cols");
- fMessenger->DeclareProperty("nRows", nRows, "Number of rows");
- fMessenger->DeclareProperty("atmosphere", isAtmosphere, "Construct Atmosphere");
-
- DefineMaterials();
-
- isAtmosphere = true;
-//************************//
-//Dimensions of Atmosphere//
-
- xWorld = 1000.*m;
- yWorld = 1000.*m;
- zWorld = 1000.*m;
-}
-
-MyDetectorConstruction::~MyDetectorConstruction()
-{}
-
-
-//***********************************//
-// Materials //
-
-
-void MyDetectorConstruction::DefineMaterials()
-{
- G4NistManager* nist = G4NistManager::Instance();
-
- worldMat = nist->FindOrBuildMaterial("G4_AIR");
-
- //**********************//
-
- G4double density0 = 1.29*kg/m3;
-
- //Set it molar mass of elemets in the Atmosphere//
-
- G4double aN = 14.01*g/mole;
- G4double aO = 16.00*g/mole;
- G4double aAr = 39.95*g/mole;
-
- //************************//
- //G4Element("Name of element", "Abreviation", atomic number, atomic masa) //
-
- G4Element *elN = new G4Element("Nitrogen", "N", 7, aN);
- G4Element *elO = new G4Element("Oxygen", "O", 8, aO);
- G4Element* elAr = new G4Element("Argon", "Ar", 18, aAr);
-
- //**********************//
- //Set it Constants //
-
- G4double f = 3; //degrees of freedom
- G4double R = 8.3144626181532; // Universal constant of gas
- G4double g0 = 9.81; // Gravity
- G4double kappa = (f + 2) / f; // parametrization
- G4double T = 293.15; // Standard temperatur at sea level
- G4double M = (0.23 * 14.01 + 0.76 * 16.0 + 0.01 * 39.95 ) / 1000.; // Average molar mass of atmosphere
-
-
- //*****************************//
-
- // Creating the air
-
- for (G4int i = 0; i < 10; i++)
- {
- std::stringstream stri;
- stri << i;
- G4double h = 2.e3 / 10. * i;
- //h is the heught of each layer
- //Divide by 5 because is the number of layers
-
- G4double density = density0 * pow((1 - (kappa - 1) / kappa * M * g0 * h / (R * T)), (1 / (kappa - 1)));
-
- Air[i] = new G4Material("G4_AIR_" + stri.str(), density, 3);
- //("G4_AIR_" + stri.str() is the name of the material. It uses a string that combines "G4_AIR_" with the value of stri.str() to create unique names for each material.)
-
-
- // Add the elements and their amount
-
- Air[i]->AddElement(elN, 76 * perCent);
- Air[i]->AddElement(elO, 23 * perCent);
- Air[i]->AddElement(elAr, 1 * perCent);
- }
-}
- //MATERIAL SESSION END
-
-//********************************//
-
-//ATMOSPHERE SET UP
-
-void MyDetectorConstruction::ConstructAtmosphere()
-{
- solidAir = new G4Box("solidAir", xWorld, yWorld, zWorld/10.);
-
-//*********************************//
-// IT create layer : Air[0],Air[1]..., Air[4]
-
- for(G4int i = 0; i < 10; i++)
- {
- logicAir[i] = new G4LogicalVolume(solidAir, Air[i], "logicAir");
-
- physAir[i] = new G4PVPlacement(0, G4ThreeVector(0, 0, zWorld/10.*2.*i - zWorld + zWorld/10.), logicAir[i], "physAir", logicWorld, false, i, true);
- }
-}
-//********************************//
-//END ATMOSPHERE
-
-//
-
-//************************************//
-
-
-G4VPhysicalVolume *MyDetectorConstruction::Construct()
-{
-
-// CREATE THE WORLD
-
- solidWorld = new G4Box("solidWorld", xWorld, yWorld, zWorld);
-
- logicWorld = new G4LogicalVolume(solidWorld, worldMat, "logicWorld");
-
- physWorld = new G4PVPlacement(0, G4ThreeVector(0., 0., 0.), logicWorld, "physWorld", 0, false, 0, true);
-
-// CREATE ATMOSPHERE
-
- if(isAtmosphere)
- ConstructAtmosphere();
-
- return physWorld;
-
-}
-
-void MyDetectorConstruction::ConstructSDandField()
-{}
-