quando gli gira va (dadada08) · Commits · Alfonso Pisapia / G4Alfo

include/DetectorConstruction.hh

+2 −2

Original line number	Original line	Diff line number	Diff line
	@@ -286,7 +286,7 @@ class DetectorConstruction : public G4VUserDetectorConstruction
	G4Material* polypropMaterial;		G4Material* polypropMaterial;
	G4Material* kapMaterial;		G4Material* kapMaterial;
	G4Material* alumMaterial;		G4Material* alumMaterial;
	G4Material* FR4;		G4Material* fr4Material;


	};		};

include/PCBParameterisation.hh

+2 −0

Original line number	Original line	Diff line number	Diff line
	@@ -53,6 +53,8 @@ class PCBParameterisation : public G4VPVParameterisation
	G4double fHeight;		G4double fHeight;
	G4double fThick;		G4double fThick;
	G4int fSide;		G4int fSide;
			G4Material* fr4Material;

	};		};

src/DetectorConstruction.cc

+18 −6

Original line number	Original line	Diff line number	Diff line
	@@ -348,13 +348,13 @@ void DetectorConstruction::DefineMaterials()
	G4Material* FR4 = new G4Material("FR4", density = 1.8*g/cm3, ncomponents=2);		G4Material* FR4 = new G4Material("FR4", density = 1.8*g/cm3, ncomponents=2);
	FR4 -> AddMaterial(G4_SILICON_DIOXIDE, fractionmass=60*perCent);		FR4 -> AddMaterial(G4_SILICON_DIOXIDE, fractionmass=60*perCent);
	FR4 -> AddMaterial(Epoxy_Resin, fractionmass=40*perCent);		FR4 -> AddMaterial(Epoxy_Resin, fractionmass=40*perCent);
			fr4Material = FR4; // G4Material::GetMaterial("FR4");

	// Default materials of the World: vacuum		// Default materials of the World: vacuum
	defaultMaterial = Vacuum;		defaultMaterial = Vacuum;
	}		}



	// Detector construction		// Detector construction
	G4VPhysicalVolume* DetectorConstruction::Construct()		G4VPhysicalVolume* DetectorConstruction::Construct()
	{		{
	@@ -365,6 +365,7 @@ G4VPhysicalVolume* DetectorConstruction::Construct()
	G4SolidStore::GetInstance()->Clean();		G4SolidStore::GetInstance()->Clean();



	// ! Experimental hall (world volume)		// ! Experimental hall (world volume)

	// Solid		// Solid
	@@ -587,10 +588,10 @@ G4VPhysicalVolume* DetectorConstruction::Construct()
	// ! PCB		// ! PCB

	// Solid		// Solid
	G4double pcb_side_x = coll_side_x *cm;		G4double pcb_side_x = coll_side_x;
	G4double pcb_side_y = coll_side_y *cm;		G4double pcb_side_y = coll_side_y;
	G4double pcb_thick = 1.6 *mm; // ! si potrebbe mettere come parametro		G4double pcb_thick = 1.6 *mm; // ! si potrebbe mettere come parametro
	G4double pcb_height = optCoupler_thick+sdd_thick+pcb_thick/2. *cm;		G4double pcb_height = optCoupler_thick+sdd_thick+pcb_thick;

	G4Box* pcb_container_box = new G4Box("pcb_container_box", pcb_side_x/2., pcb_side_y/2., pcb_thick/2.);		G4Box* pcb_container_box = new G4Box("pcb_container_box", pcb_side_x/2., pcb_side_y/2., pcb_thick/2.);

	@@ -606,20 +607,31 @@ G4VPhysicalVolume* DetectorConstruction::Construct()
	false,		false,
	0);		0);


	// Solid		// Solid
	G4double pcb_elem_side_x = 3.*mm;		G4double pcb_elem_side_x = 3.*mm;
	G4double pcb_elem_side_y = pcb_side_y;		G4double pcb_elem_side_y = pcb_side_y;
			G4cout << "Materiale FR4: " << fr4Material->GetName() << G4endl;

	G4Box* pcb_box = new G4Box("pcb_box", pcb_elem_side_x/2., pcb_elem_side_y/2., pcb_thick/2.);		G4Box* pcb_box = new G4Box("pcb_box", pcb_elem_side_x/2., pcb_elem_side_y/2., pcb_thick/2.);

	// Logical		// Logical
	pcb_log = new G4LogicalVolume(pcb_box, FR4, "pcb_log", 0, 0, 0);		pcb_log = new G4LogicalVolume(pcb_box, fr4Material, "pcb_log", 0, 0, 0);

	// Physical		// Physical
	G4int n_rep = (n_side-1)/2; // repliche di PCB usando il numero di scintillatori per side		G4int n_rep = (n_side-1)/2; // repliche di PCB usando il numero di scintillatori per side

	pcbParam = new PCBParameterisation(pcb_elem_side_x, pcb_elem_side_y, pcb_height, pcb_thick, n_rep);		pcbParam = new PCBParameterisation(pcb_elem_side_x, pcb_elem_side_y, pcb_height, pcb_thick, n_rep);

			G4cout << " pcb_elem_side_x " << pcb_elem_side_x/mm << G4endl;
			G4cout << " pcb_elem_side_y " << pcb_elem_side_y/mm << G4endl;
			G4cout << " pcb_height " << pcb_height << G4endl;
			G4cout << " pcb_thick " << pcb_thick/mm << G4endl;

			G4cout << " pcb_sie_x " << pcb_side_x << G4endl;
			G4cout << " n_rep " << n_rep << G4endl;


	pcb_phys = new G4PVParameterised("pcb_phys", // their name		pcb_phys = new G4PVParameterised("pcb_phys", // their name
	pcb_log, // their logical volume		pcb_log, // their logical volume
	pcb_container_log, // Mother logical volume		pcb_container_log, // Mother logical volume

src/PCBParameterisation.cc

+5 −4

Original line number	Original line	Diff line number	Diff line
	@@ -30,10 +30,11 @@ void PCBParameterisation::ComputeTransformation (const G4int copyNo, G4VPhysical

	xPos = (-fLength)/2. +1.cm + copyNo 1.*cm; //sono al centro della barra capisci se indietreggiare di mezza barra		xPos = (-fLength)/2. +1.cm + copyNo 1.*cm; //sono al centro della barra capisci se indietreggiare di mezza barra

			G4cout << " copyNo " << copyNo << G4endl;
	G4cout << "* DEBUG *: Replica number " << copyNo << G4endl;		G4cout << " fSpacing " << fSpacing << G4endl;
	G4cout << "* DEBUG *: " << xPos << " " << yPos << " " << zPos << G4endl;		G4cout << " zPos " << zPos << G4endl;
			G4cout << " fLength " << fLength << G4endl;
			G4cout << " fSpacing " << fSpacing << G4endl;

	physVol->SetTranslation(G4ThreeVector(xPos, yPos, zPos));		physVol->SetTranslation(G4ThreeVector(xPos, yPos, zPos));

xgis_M7

(471 KiB)

File changed.

No diff preview for this file type.

View original file

View changed file