Commit 770d69d6 authored by Valentina Fioretti's avatar Valentina Fioretti
Browse files

spo coating 0 = only Si

parent 14baca54
Loading
Loading
Loading
Loading
+10 −42
Original line number Original line Diff line number Diff line
@@ -295,27 +295,19 @@ G4VPhysicalVolume* EXACRAD::ConstructGeometry(G4VPhysicalVolume* World_phys) {
            
            
            G4double spo_plate_t = 0.140*mm; //
            G4double spo_plate_t = 0.140*mm; //
            G4double spo_plate_t_edge = 0.775*mm; //
            G4double spo_plate_t_edge = 0.775*mm; //
            G4double spo_coating1_t_mean = 1*um; //
            G4double spo_coating1_t_min = 0.6*um; //
            G4double spo_coating1_t_max = 1.2*um; //
            G4double spo_coating2_t = 10*nm; //
            G4double spo_coating3_t = 7*nm; //


            G4double plate_yside = 49.0*mm;
            G4double plate_yside = 49.0*mm;
            G4double plate_xside = 110.0*mm;
            G4double plate_xside = 110.0*mm;
            
            
            G4int spo_mother_copyno = 30;
            G4int spo_mother_copyno = 30;
            G4int spo_plate_copyno = 31;
            G4int spo_plate_copyno = 31;
            G4int spo_coating1_copyno = 32;
            G4int spo_coating2_copyno = 33;
            G4int spo_coating3_copyno = 34;
            
            
            
            
            if (coating1_type == 0){
            if (coating1_type == 0){
                
                
                // --------- mother target
                // --------- mother target
                // mother height
                // mother height
                G4double mother_target_t = spo_plate_t_edge + spo_coating1_t_mean + spo_coating2_t;
                G4double mother_target_t = spo_plate_t_edge;
                G4Box* mother_target = new G4Box("mother_target",
                G4Box* mother_target = new G4Box("mother_target",
                                                 plate_xside/2.,
                                                 plate_xside/2.,
                                                 plate_yside/2.,
                                                 plate_yside/2.,
@@ -555,42 +547,18 @@ G4VPhysicalVolume* EXACRAD::ConstructGeometry(G4VPhysicalVolume* World_phys) {
                log_rib10->SetUserLimits(target_limit);
                log_rib10->SetUserLimits(target_limit);
                region_target->AddRootLogicalVolume(log_rib10);
                region_target->AddRootLogicalVolume(log_rib10);
            
            
                // ------------------------------------------------
                // ----- coating

                // ----- coating 1
                G4Box* coating1 = new G4Box("coating1",
                                             plate_xside/2.,
                                             plate_yside/2.,
                                             spo_coating1_t_mean/2.);
                
                G4LogicalVolume* log_coating1 = new G4LogicalVolume(coating1, spo_coating1_mat, "log_coating1");
                G4VPhysicalVolume* phys_coating1 = new G4PVPlacement(0, G4ThreeVector(0.0,0.0, -(mother_target_t/2.) + spo_plate_t_edge + spo_coating1_t_mean/2.), log_coating1, "phys_coating1", log_mother_target, false, spo_coating1_copyno);
                
                G4VisAttributes* visCoating1 = new G4VisAttributes(G4Colour::Yellow());
                log_coating1->SetVisAttributes(visCoating1);
                gm.AddXYZDetector( log_coating1);
                log_coating1->SetUserLimits(target_limit);
                region_target->AddRootLogicalVolume(log_coating1);

                // ----- coating 2
                G4Box* coating2 = new G4Box("coating2",
                                            plate_xside/2.,
                                            plate_yside/2.,
                                            spo_coating2_t/2.);
                
                G4LogicalVolume* log_coating2 = new G4LogicalVolume(coating2, spo_coating2_mat, "log_coating2");
                G4VPhysicalVolume* phys_coating2 = new G4PVPlacement(0, G4ThreeVector(0.0,0.0, -(mother_target_t/2.) + spo_plate_t_edge + spo_coating1_t_mean + spo_coating2_t/2.), log_coating2, "phys_coating2", log_mother_target, false, spo_coating2_copyno);
                
                G4VisAttributes* visCoating2 = new G4VisAttributes(G4Colour::Brown());
                log_coating2->SetVisAttributes(visCoating2);
                gm.AddXYZDetector( log_coating2);
                log_coating2->SetUserLimits(target_limit);
                region_target->AddRootLogicalVolume(log_coating2);

            }
            }
            if (coating1_type == 1){
            if (coating1_type == 1){
                
                
                G4double spo_coating1_t_mean = 1*um; //
                G4double spo_coating1_t_min = 0.6*um; //
                G4double spo_coating1_t_max = 1.2*um; //
                G4double spo_coating2_t = 10*nm; //
                G4double spo_coating3_t = 7*nm; //
                
                G4int spo_coating1_copyno = 32;
                G4int spo_coating2_copyno = 33;
                G4int spo_coating3_copyno = 34;


                // --------- mother target
                // --------- mother target
                // mother height
                // mother height