Message: G4CutTubs why some part cannot tracking? Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None G4CutTubs why some part cannot tracking? 

Forum: Geometry
Date: 03 Oct, 2014
From: T.D.Trong <T.D.Trong>

Dear G4 experts,

I am building an ion chamber with circle electrodes title 30deg. I am using G4CutTubs on geant4 4.10.0 p02

   G4double MS_cell_radius = 80.*mm;
    G4double MS_cell_thick =  20.*mm;
    G4double MS_ring_thick =   2.*mm;
    G4double MS_foil_thick =   2.*um;
    G4int nofLayers = 8;

    G4ThreeVector pLowNorm = G4ThreeVector(0.,-std::tan(fMS_tilt_angle),-1.);
    G4ThreeVector pHigNorm = G4ThreeVector(0., std::tan(fMS_tilt_angle), 1.);

        G4VSolid* outter_elipS
        = new G4EllipticalTube("outter_elip",             // its name
                               MS_cell_radius,MS_cell_radius*std::cos(fMS_tilt_angle),MS_cell_thick +1.4 * MS_cell_radius*std::sin(fMS_tilt_angle)); // its size

        G4VSolid* MS_cellrbeS0
        = new G4CutTubs("MS_cell",            // its name
                        0.,MS_cell_radius, MS_cell_thick/2.,0.,360.*deg,pLowNorm,pHigNorm); // its size
        G4VSolid* MS_cellS
        =new G4IntersectionSolid("MS_cell2",MS_cellrbeS0,outter_elipS);

       G4LogicalVolume* fMScellLV
    = new G4LogicalVolume(
                          MS_cellS,        // its solid
                          MScellM, // its material
                          "MS_cell_LV");          // its name

    for (G4int i=0; i<2*nofLayers ; i++) {
                G4double zpos = (i - nofLayers) *(MS_foil_thick)-7.*mm;
                new G4PVPlacement(
                          0,                // no rotation
                          G4ThreeVector(0., 0., zpos), //  its position
                          fMScellLV,       // its logical volume
                          MS_cell,           // its name
                          world_LV,          // its mother  volume
                          false,            // no boolean operation
                          i,                // copy number
                          fCheckOverlaps);  // checking overlaps
    }
world_LV is a  vacuum box 
MScellM is 1amt gas ( air, CO2,...)    

When I looked visualize and also energy loss, the track path has y_axis < 0 there are some event hit to some cells but not hit to neighborhood cell when y move to near the edge in negative side.

 But when I set fMS_tilt_angle = 0, every cells is hits (Please look at the attach file for more detail).

Who can tell me why?

Thank you,

   Attachment:
      http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2014/10/02/20.33-52434-geant.pdf

Inline Depth:
 1 1
 All All
Outline Depth:
 1 1
 2 2
 All All
Add message: (add)

1 None: Re: G4CutTubs why some part cannot tracking?   (Tatiana Nikitina - 05 Oct, 2014)
(_ None: Re: G4CutTubs why some part cannot tracking?   (T.D.Trong - 05 Oct, 2014)
(_ None: Re: G4CutTubs why some part cannot tracking?   (Tatiana Nikitina - 06 Oct, 2014)
(_ None: Re: G4CutTubs why some part cannot tracking?   (T.D.Trong - 06 Oct, 2014)
(_ None: Re: G4CutTubs why some part cannot tracking?   (Tatiana Nikitina - 07 Oct, 2014)
 Add Message Add Message
to: "G4CutTubs why some part cannot tracking?"

 Subscribe Subscribe

This site runs SLAC HyperNews version 1.11-slac-98, derived from the original HyperNews


[ Geant 4 Home | Geant 4 HyperNews | Search | Request New Forum | Feedback ]