QweakSimTarget Class Reference

Definition of the Target. More...

#include <QweakSimTarget.hh>

Collaboration diagram for QweakSimTarget:

Public Member Functions
	QweakSimTarget (QweakSimUserInformation *myUI)
	~QweakSimTarget ()
void	ConstructComponent (G4VPhysicalVolume *)
void	DestroyComponent ()
void	SetTargetCenterPositionInZ (G4double)
G4double	GetTargetCenterPositionInZ ()
void	SetTargetMaterial (G4String)
G4String	GetTargetMaterial ()
void	SetTargetCellMaterial (G4String)
G4String	GetTargetCellMaterial ()
void	SetTargetEntranceWindowMaterial (G4String)
G4String	GetTargetEntranceMaterial ()
void	SetTargetExitWindowMaterial (G4String)
G4String	GetTargetExitWindowMaterial ()
void	SetTargetExitWindowNippleMaterial (G4String)
G4String	GetTargetExitWindowNippleMaterial ()
G4LogicalVolume *	getTargetLogicalVolume ()
G4VPhysicalVolume *	getTargetPhysicalVolume ()
void	SetTargetLength (G4double)
G4double	GetTargetLength ()
void	SetTarget (G4String)
void	SetTargetEntranceWindowLength (G4double)
void	SetTargetExitWindowLength (G4double)
void	SetTargetExitWindowNippleLength (G4double)
G4double	CalculateLuminosity (G4double mass, G4double density, G4double length, G4double pweight)

Private Member Functions
void	CalculateTargetPositions ()
void	ConstructTargetContainer ()
void	ConstructScatteringChamberWindow ()
void	ConstructTargetCell ()
void	ConstructTargetMaterial ()
void	ConstructTargetEntranceWindow ()
void	ConstructTargetExitWindow ()
void	ConstructTargetExitWindowNipple ()

Private Attributes
G4VPhysicalVolume *	theMotherPV
QweakSimMaterial *	pMaterial
G4LogicalVolume *	TargetCell_Logical
G4VPhysicalVolume *	TargetCell_Physical
G4Material *	TargetCell_Material
G4LogicalVolume *	TargetEntranceWindow_Logical
G4VPhysicalVolume *	TargetEntranceWindow_Physical
G4Material *	TargetEntranceWindow_Material
G4LogicalVolume *	TargetExitWindow_Logical
G4VPhysicalVolume *	TargetExitWindow_Physical
G4Material *	TargetExitWindow_Material
G4LogicalVolume *	TargetExitWindowNipple_Logical
G4VPhysicalVolume *	TargetExitWindowNipple_Physical
G4Material *	TargetExitWindowNipple_Material
G4LogicalVolume *	ScatteringChamberWindow_Logical
G4VPhysicalVolume *	ScatteringChamberWindow_Physical
G4Material *	ScatteringChamberWindow_Material
G4LogicalVolume *	TargetMaterial_Logical
G4VPhysicalVolume *	TargetMaterial_Physical
G4Material *	TargetMaterial_Material
G4LogicalVolume *	TargetContainer_Logical
G4VPhysicalVolume *	TargetContainer_Physical
G4Material *	TargetContainer_Material
G4double	targetCellEntranceWindowThickness
G4double	targetCellExitWindowThickness
G4double	targetCellExitWindowNippleThickness
G4double	targetCellWallThickness
G4double	targetCellInnerLength
G4double	targetCellOuterLength
G4double	targetCellFrontRadiusMin
G4double	targetCellFrontInnerRadiusMax
G4double	targetCellFrontOuterRadiusMax
G4double	targetCellBackRadiusMin
G4double	targetCellBackInnerRadiusMax
G4double	targetCellBackOuterRadiusMax
G4double	targetCellExitWindowNippleRadius
G4double	ScatteringChamberWindowRadius
G4double	ScatteringChamberWindowThickness
G4double	targetCellStartingPhi
G4double	targetCellDeltaPhi
G4double	targetZPos
QweakSimTargetMessenger *	targetMessenger
QweakSimUserInformation *	myUserInfo
G4ThreeVector	positionTarget
G4ThreeVector	positionTargetEntranceWindow
G4ThreeVector	positionTargetExitWindow
G4ThreeVector	positionScatteringChamberWindow
G4VSensitiveDetector *	TargetSD

Detailed Description

Definition of the Target.

Placeholder for a long explaination

Definition at line 65 of file QweakSimTarget.hh.

Constructor & Destructor Documentation

QweakSimTarget::QweakSimTarget

(

QweakSimUserInformation *

myUI

)

Definition at line 36 of file QweakSimTarget.cc.

References CalculateLuminosity(), QweakSimMaterial::GetInstance(), QweakSimMaterial::GetMaterial(), inch, mil, myUserInfo, pMaterial, ScatteringChamberWindow_Logical, ScatteringChamberWindow_Material, ScatteringChamberWindow_Physical, ScatteringChamberWindowRadius, ScatteringChamberWindowThickness, TargetCell_Logical, TargetCell_Material, TargetCell_Physical, targetCellBackInnerRadiusMax, targetCellBackOuterRadiusMax, targetCellBackRadiusMin, targetCellDeltaPhi, targetCellEntranceWindowThickness, targetCellExitWindowNippleRadius, targetCellExitWindowNippleThickness, targetCellExitWindowThickness, targetCellFrontInnerRadiusMax, targetCellFrontOuterRadiusMax, targetCellFrontRadiusMin, targetCellInnerLength, targetCellOuterLength, targetCellStartingPhi, targetCellWallThickness, QweakSimUserInformation::TargetCenterPositionZ, TargetContainer_Logical, TargetContainer_Material, TargetContainer_Physical, TargetEntranceWindow_Logical, TargetEntranceWindow_Material, TargetEntranceWindow_Physical, QweakSimUserInformation::TargetEntranceWindowThickness, TargetExitWindow_Logical, TargetExitWindow_Material, TargetExitWindow_Physical, TargetExitWindowNipple_Logical, TargetExitWindowNipple_Material, TargetExitWindowNipple_Physical, QweakSimUserInformation::TargetExitWindowNippleThickness, QweakSimUserInformation::TargetExitWindowThickness, QweakSimUserInformation::TargetLength, QweakSimUserInformation::TargetLuminosityDSALDummy2, QweakSimUserInformation::TargetLuminosityDSALDummy4, QweakSimUserInformation::TargetLuminosityDSALDummy4_Al, QweakSimUserInformation::TargetLuminosityDSALDummy4_Cr, QweakSimUserInformation::TargetLuminosityDSALDummy4_Cu, QweakSimUserInformation::TargetLuminosityDSALDummy4_Fe, QweakSimUserInformation::TargetLuminosityDSALDummy4_Mg, QweakSimUserInformation::TargetLuminosityDSALDummy4_Mn, QweakSimUserInformation::TargetLuminosityDSALDummy4_Si, QweakSimUserInformation::TargetLuminosityDSALDummy4_Ti, QweakSimUserInformation::TargetLuminosityDSALDummy4_Zn, QweakSimUserInformation::TargetLuminosityDSALDummy8, QweakSimUserInformation::TargetLuminosityDSALWindow, QweakSimUserInformation::TargetLuminosityDSCDummy, QweakSimUserInformation::TargetLuminosityLH2, QweakSimUserInformation::TargetLuminosityUSALDummy1, QweakSimUserInformation::TargetLuminosityUSALDummy2, QweakSimUserInformation::TargetLuminosityUSALDummy4, QweakSimUserInformation::TargetLuminosityUSALWindow, QweakSimUserInformation::TargetLuminosityUSCDummy, TargetMaterial_Logical, TargetMaterial_Material, TargetMaterial_Physical, targetMessenger, TargetSD, QweakSimUserInformation::TargetThicknessDSALDummy2, QweakSimUserInformation::TargetThicknessDSALDummy4, QweakSimUserInformation::TargetThicknessDSALDummy8, QweakSimUserInformation::TargetThicknessDSCDummy, QweakSimUserInformation::TargetThicknessUSALDummy1, QweakSimUserInformation::TargetThicknessUSALDummy2, QweakSimUserInformation::TargetThicknessUSALDummy4, QweakSimUserInformation::TargetThicknessUSCDummy, targetZPos, and theMotherPV.

{
    G4cout << G4endl << "###### Calling QweakSimTarget::QweakTarget() " << G4endl << G4endl;

    myUserInfo = myUI;

    theMotherPV               = NULL;

    TargetContainer_Logical      = NULL;
    TargetContainer_Physical     = NULL;
    TargetContainer_Material     = NULL;

    TargetCell_Logical      = NULL;
    TargetCell_Physical     = NULL;
    TargetCell_Material     = NULL;

    TargetEntranceWindow_Logical      = NULL;
    TargetEntranceWindow_Physical     = NULL;
    TargetEntranceWindow_Material     = NULL;

    TargetExitWindow_Logical      = NULL;
    TargetExitWindow_Physical     = NULL;
    TargetExitWindow_Material     = NULL;

    TargetExitWindowNipple_Logical  = NULL;
    TargetExitWindowNipple_Physical = NULL;
    TargetExitWindowNipple_Material = NULL;

    ScatteringChamberWindow_Logical  = NULL;
    ScatteringChamberWindow_Physical = NULL;
    ScatteringChamberWindow_Material = NULL;

    TargetMaterial_Logical  = NULL;
    TargetMaterial_Physical = NULL;
    TargetMaterial_Material = NULL;


    pMaterial = QweakSimMaterial::GetInstance();

    TargetContainer_Material      = pMaterial->GetMaterial("Vacuum");
    TargetCell_Material           = pMaterial->GetMaterial("Aluminum");
    TargetEntranceWindow_Material = pMaterial->GetMaterial("Aluminum");
    TargetExitWindow_Material     = pMaterial->GetMaterial("Aluminum");
    TargetExitWindowNipple_Material = pMaterial->GetMaterial("Aluminum");
    ScatteringChamberWindow_Material = pMaterial->GetMaterial("Aluminum");
    TargetMaterial_Material       = pMaterial->GetMaterial("H2Liquid");

    // define target geometry values

    targetCellEntranceWindowThickness = 4.105*mil; //Avg of Run 1 & 2 (Greg's Target Elog 19)
    targetCellExitWindowThickness =  29.06*mil; //Greg's Target Elog 19 
    targetCellExitWindowNippleThickness =  5.165*mil; //Avg of Run 1 & 2 (Greg's Target Elog 19)
    targetCellWallThickness   = 10.0*mil;

    targetCellExitWindowNippleRadius = 0.5*0.622*inch;

    targetCellInnerLength     = myUserInfo->TargetLength;
    targetCellOuterLength     = targetCellInnerLength
                                + targetCellEntranceWindowThickness
                                + targetCellExitWindowThickness;  // Full length of Target

    targetCellFrontRadiusMin        =  0.0*cm;
    targetCellFrontOuterRadiusMax   =  0.5*2.730*inch;
    targetCellFrontInnerRadiusMax   = targetCellFrontOuterRadiusMax - targetCellWallThickness;

    targetCellBackRadiusMin        =  targetCellExitWindowNippleRadius;  //0.0*cm;
    targetCellBackOuterRadiusMax   =  0.5*7.5*inch;
    targetCellBackInnerRadiusMax   = targetCellBackOuterRadiusMax - targetCellWallThickness;

    targetCellStartingPhi =  0.0*deg;
    targetCellDeltaPhi    = 360*deg;

    targetZPos        = myUserInfo->TargetCenterPositionZ; //-650.0*cm;

    ScatteringChamberWindowRadius = 0.5*23.5*inch;
    ScatteringChamberWindowThickness = 35*mil;

    myUserInfo->TargetEntranceWindowThickness = targetCellEntranceWindowThickness;
    myUserInfo->TargetExitWindowThickness = targetCellExitWindowThickness;
    myUserInfo->TargetExitWindowNippleThickness = targetCellExitWindowNippleThickness;

    //Thicknesses taken from 2012 target survey (https://qweak.jlab.org/elog/Target/21)
    myUserInfo->TargetThicknessUSALDummy1    = 0.8812*mm;
    myUserInfo->TargetThicknessUSALDummy2    = 1.7987*mm;
    myUserInfo->TargetThicknessUSALDummy4    = 3.6030*mm;
    myUserInfo->TargetThicknessDSALDummy2    = 1.8637*mm;
    myUserInfo->TargetThicknessDSALDummy4    = 3.6828*mm;
    myUserInfo->TargetThicknessDSALDummy8    = 7.1980*mm;
    myUserInfo->TargetThicknessUSCDummy      = 0.9973*mm;
    myUserInfo->TargetThicknessDSCDummy      = 3.1876*mm;

    G4double densityLH2 = 0.0713 /(cm*cm*cm);  // [g/cm^3] 
    G4double densityAL  = 2.804   /(cm*cm*cm);  // Avg density, Greg's Target Elog 21
    G4double densityUSC = 1.70   /(cm*cm*cm);  //
    G4double densityDSC = 2.205  /(cm*cm*cm); //

    // Molar masses taken from PDG:  
    // J. Beringer et al. (Particle Data Group), Phys. Rev. D86, 010001 (2012).
    G4double massLH2 = 1.00794;    // [g/mol]
    G4double massAL  = 26.9815386; // [g/mol]
    G4double massC   = 12.0107;    // [g/mol]
    //Added mass for the elements in the alloy (PDG - 2015)
    G4double massZn  = 65.38; //[g/mol]
    G4double massMg  = 24.3050; //[g/mol]
    G4double massCu  = 63.546; //[g/mol]
    G4double massCr  = 51.9961; //[g/mol]
    G4double massFe  = 55.845; //[g/mol]
    G4double massSi  = 28.0855; //[g/mol]
    G4double massMn  = 54.938044; //[g/mol]
    G4double massTi  = 47.867; //[g/mol]

    //Weighting factors for calculation of luminosity of elements in the aluminum alloy
    //Only implementing for DS 4% Al dummy target. Weighting factors take from Greg's 
    //Target ELOG 30 - K.Bartlett 2017/7/13
    G4double weight_std = 1.0;//Standard weight for hydrogen and other solid targets.
    G4double weight_Al = 0.8870;//Aluminum's % by weight[unitless]
    G4double weight_Zn = 0.0630;//Zinc's % by weight[unitless]
    G4double weight_Mg = 0.0270;//Magnesium's % by weight
    G4double weight_Cu = 0.0180;//Copper's % by weight
    G4double weight_Cr = 0.0021;//Chromium's % by weight
    G4double weight_Fe = 0.0012;//Iron's % by weight
    G4double weight_Si = 0.0010;//Silicon's % by weight
    G4double weight_Mn = 0.0004;//Manganese's % by weight
    G4double weight_Ti = 0.0003;//Titanium's % by weight


    myUserInfo->TargetLuminosityLH2 = CalculateLuminosity(massLH2, densityLH2, myUserInfo->TargetLength, weight_std);
    myUserInfo->TargetLuminosityUSALWindow = CalculateLuminosity(massAL, densityAL, myUserInfo->TargetEntranceWindowThickness, weight_std);
    myUserInfo->TargetLuminosityDSALWindow = CalculateLuminosity(massAL, densityAL, myUserInfo->TargetExitWindowNippleThickness, weight_std);
    myUserInfo->TargetLuminosityUSALDummy1 = CalculateLuminosity(massAL, densityAL, myUserInfo->TargetThicknessUSALDummy1, weight_std);
    myUserInfo->TargetLuminosityUSALDummy2 = CalculateLuminosity(massAL, densityAL, myUserInfo->TargetThicknessUSALDummy2, weight_std);
    myUserInfo->TargetLuminosityUSALDummy4 = CalculateLuminosity(massAL, densityAL, myUserInfo->TargetThicknessUSALDummy4, weight_std);
    myUserInfo->TargetLuminosityDSALDummy2 = CalculateLuminosity(massAL, densityAL, myUserInfo->TargetThicknessDSALDummy2, weight_std);
    myUserInfo->TargetLuminosityDSALDummy4 = CalculateLuminosity(massAL, densityAL, myUserInfo->TargetThicknessDSALDummy4, weight_std);
    myUserInfo->TargetLuminosityDSALDummy8 = CalculateLuminosity(massAL, densityAL, myUserInfo->TargetThicknessDSALDummy8, weight_std);
    myUserInfo->TargetLuminosityUSCDummy = CalculateLuminosity(massC, densityUSC, myUserInfo->TargetThicknessUSCDummy, weight_std);
    myUserInfo->TargetLuminosityDSCDummy = CalculateLuminosity(massC, densityDSC, myUserInfo->TargetThicknessDSCDummy, weight_std);

    //Calculate luminosities for elements in alloy for DS 4%.
    myUserInfo->TargetLuminosityDSALDummy4_Al = CalculateLuminosity(massAL, densityAL, myUserInfo->TargetThicknessDSALDummy4, weight_Al);
    myUserInfo->TargetLuminosityDSALDummy4_Zn = CalculateLuminosity(massZn, densityAL, myUserInfo->TargetThicknessDSALDummy4, weight_Zn);
    myUserInfo->TargetLuminosityDSALDummy4_Mg = CalculateLuminosity(massMg, densityAL, myUserInfo->TargetThicknessDSALDummy4, weight_Mg);
    myUserInfo->TargetLuminosityDSALDummy4_Cu = CalculateLuminosity(massCu, densityAL, myUserInfo->TargetThicknessDSALDummy4, weight_Cu);
    myUserInfo->TargetLuminosityDSALDummy4_Cr = CalculateLuminosity(massCr, densityAL, myUserInfo->TargetThicknessDSALDummy4, weight_Cr);
    myUserInfo->TargetLuminosityDSALDummy4_Fe = CalculateLuminosity(massFe, densityAL, myUserInfo->TargetThicknessDSALDummy4, weight_Fe);
    myUserInfo->TargetLuminosityDSALDummy4_Si = CalculateLuminosity(massSi, densityAL, myUserInfo->TargetThicknessDSALDummy4, weight_Si);
    myUserInfo->TargetLuminosityDSALDummy4_Mn = CalculateLuminosity(massMn, densityAL, myUserInfo->TargetThicknessDSALDummy4, weight_Mn);
    myUserInfo->TargetLuminosityDSALDummy4_Ti = CalculateLuminosity(massTi, densityAL, myUserInfo->TargetThicknessDSALDummy4, weight_Ti);



    TargetSD = NULL;

    targetMessenger = new QweakSimTargetMessenger(this);
}

Here is the call graph for this function:

QweakSimTarget::~QweakSimTarget

(

)

Definition at line 193 of file QweakSimTarget.cc.

References targetMessenger.

{
    delete targetMessenger;
}

Member Function Documentation

G4double QweakSimTarget::CalculateLuminosity	(	G4double	mass,
		G4double	density,
		G4double	length,
		G4double	pweight
	)

Definition at line 923 of file QweakSimTarget.cc.

Referenced by QweakSimTarget().

{
  G4double Lum = 0.0;            // Luminosity
  G4double N_b = 6.241e12;       // [Hz/uA] # of particles in the beam (from definition of ampere)
  G4double N_A = 6.02214129e23;  // Avagadro's number

  Lum = N_b*length*pweight*density*N_A/mass;  //  units [Hz/(uA*mm^2)]

  Lum *= 1e-31;  //  Conversion from [Hz/(uA*mm^2)] -> [kHz/(uA*ub)]

  return Lum;
}

Here is the caller graph for this function:

void QweakSimTarget::CalculateTargetPositions ( )

private

Definition at line 266 of file QweakSimTarget.cc.

References positionScatteringChamberWindow, positionTarget, positionTargetEntranceWindow, positionTargetExitWindow, targetCellEntranceWindowThickness, targetCellExitWindowThickness, targetCellInnerLength, targetCellOuterLength, and TargetContainer_Physical.

Referenced by ConstructComponent(), SetTargetEntranceWindowLength(), SetTargetExitWindowLength(), SetTargetExitWindowNippleLength(), and SetTargetLength().

                                             {

    // define Target position
    positionTarget = G4ThreeVector(0,0,0);
    positionTargetEntranceWindow = G4ThreeVector(0,0,-0.5*targetCellInnerLength - 0.5*targetCellEntranceWindowThickness);
    positionTargetExitWindow = G4ThreeVector(0,0, 0.5*targetCellInnerLength + 0.5*targetCellExitWindowThickness);
    positionScatteringChamberWindow = G4ThreeVector(0,0,0.5*targetCellInnerLength + 0.5*targetCellExitWindowThickness + 45.0*cm); // Peiqing: should be located at z=-583.41 cm

    targetCellOuterLength     = targetCellInnerLength
                                + targetCellEntranceWindowThickness
                                + targetCellExitWindowThickness;  // Full length of Target

    if (TargetContainer_Physical) {
      G4cout << "TargetContainer_Physical = " << TargetContainer_Physical->GetTranslation()/cm << " cm" << G4endl;
    }

    G4cout << "targetCellEntranceWindowThickness: " << targetCellEntranceWindowThickness/cm << " cm" << G4endl;
    G4cout << "positionTargetEntranceWindow:      " << positionTargetEntranceWindow/cm << " cm" << G4endl;
    G4cout << "positionTarget:                    " << positionTarget/cm << " cm" << G4endl;
    G4cout << "positionTargetExitWindow:          " << positionTargetExitWindow/cm << " cm" << G4endl;
    G4cout << "targetCellExitWindowThickness:     " << targetCellExitWindowThickness/cm << " cm" << G4endl;
    G4cout << "positionScatteringChamberWindow:   " << positionScatteringChamberWindow/cm << " cm" << G4endl;
    G4cout << "targetCellInnerLength = " << targetCellInnerLength/cm << " cm" << G4endl;
    G4cout << "targetCellOuterLength = " << targetCellOuterLength/cm << " cm" << G4endl;
}

Here is the caller graph for this function:

void QweakSimTarget::ConstructComponent

(

G4VPhysicalVolume *

MotherVolume

)

Definition at line 200 of file QweakSimTarget.cc.

References CalculateTargetPositions(), ConstructScatteringChamberWindow(), ConstructTargetCell(), ConstructTargetContainer(), ConstructTargetEntranceWindow(), ConstructTargetExitWindow(), ConstructTargetExitWindowNipple(), ConstructTargetMaterial(), ScatteringChamberWindow_Logical, TargetCell_Logical, TargetContainer_Logical, TargetEntranceWindow_Logical, TargetExitWindow_Logical, TargetMaterial_Logical, TargetSD, and theMotherPV.

{
    G4cout << G4endl << "###### Calling QweakSimTarget::ConstructComponent() " << G4endl << G4endl;

    theMotherPV = MotherVolume;

    CalculateTargetPositions();

    ConstructTargetContainer(); // scattering chamber
    ConstructScatteringChamberWindow();

    ConstructTargetCell();      // Al cell without end caps
    ConstructTargetMaterial();  // LH2 for production target

    ConstructTargetEntranceWindow();
    ConstructTargetExitWindow();
    ConstructTargetExitWindowNipple();

//--------------------------------------

    G4cout << G4endl << "###### QweakSimTarget: Setting Attributes " << G4endl << G4endl;

    G4Colour  blue  (0.,0.,1.);
    G4Colour  red   (1.,0.,0.);

    G4VisAttributes* TargetContainer_VisAtt = new G4VisAttributes(red);
    TargetContainer_VisAtt -> SetVisibility(false);
    //TargetContainer_VisAtt -> SetForceWireframe(true);
    TargetContainer_Logical -> SetVisAttributes(TargetContainer_VisAtt);

    G4VisAttributes* TargetCell_VisAtt = new G4VisAttributes(blue);
    TargetCell_VisAtt -> SetVisibility(true);
    //TargetCell_VisAtt -> SetForceWireframe(true);
    TargetCell_Logical -> SetVisAttributes(TargetCell_VisAtt);

    G4VisAttributes* TargetWindow_VisAtt = new G4VisAttributes(blue);
    TargetWindow_VisAtt -> SetVisibility(true);
    //TargetWindow_VisAtt -> SetForceWireframe(true);
    TargetEntranceWindow_Logical -> SetVisAttributes(TargetWindow_VisAtt);
    TargetExitWindow_Logical -> SetVisAttributes(TargetWindow_VisAtt);

    G4VisAttributes* ScatteringChamberWindow_VisAtt = new G4VisAttributes(red);
    ScatteringChamberWindow_VisAtt -> SetVisibility(true);
    //ScatteringChamberWindow_VisAtt -> SetForceWireframe(true);
    ScatteringChamberWindow_Logical -> SetVisAttributes(ScatteringChamberWindow_VisAtt);

    G4VisAttributes* TargetMaterial_VisAtt = new G4VisAttributes(red);
    TargetMaterial_VisAtt -> SetVisibility(true);
    //TargetVisAtt -> SetForceWireframe(true);
    TargetMaterial_Logical -> SetVisAttributes(TargetMaterial_VisAtt);

    ///////////////////////////////////////////
    // Define Sensitive Detectors to Target  //
    ///////////////////////////////////////////

    G4SDManager* SDman = G4SDManager::GetSDMpointer();

    TargetSD = new QweakSimTarget_DetectorSD("TargetSD");
    SDman->AddNewDetector(TargetSD);
    ScatteringChamberWindow_Logical->SetSensitiveDetector(TargetSD);


    G4cout << G4endl << "###### Leaving QweakSimTarget::ConstructComponent() " << G4endl << G4endl;

} // end of  QweakSimTarget::ConstructComponent()

Here is the call graph for this function:

void QweakSimTarget::ConstructScatteringChamberWindow ( )

private

Definition at line 509 of file QweakSimTarget.cc.

References positionScatteringChamberWindow, pSurfChk, ScatteringChamberWindow_Logical, ScatteringChamberWindow_Material, ScatteringChamberWindow_Physical, ScatteringChamberWindowRadius, ScatteringChamberWindowThickness, targetCellDeltaPhi, targetCellStartingPhi, and TargetContainer_Physical.

Referenced by ConstructComponent().

{
    G4cout << G4endl << "###### Calling QweakSimTarget::ConstructScatteringChamberWindow() " << G4endl << G4endl;

    //*********************** Define Target scattering chamber Vacuum Window ************************
    G4cout << G4endl << "###### QweakSimTarget: Define ScatteringChamberWindow_Solid " << G4endl << G4endl;
    G4Tubs* ScatteringChamberWindow_Solid    = new G4Tubs("ScatteringChamberWindow_Sol",
            0.,
            ScatteringChamberWindowRadius,
            0.5*ScatteringChamberWindowThickness,
            targetCellStartingPhi,
            targetCellDeltaPhi);

    // define ScatteringChamberWindow logical volume
    G4cout << G4endl << "###### QweakSimTarget: Define ScatteringChamberWindow_Logical " << G4endl << G4endl;

    ScatteringChamberWindow_Logical  = new G4LogicalVolume(ScatteringChamberWindow_Solid,
            ScatteringChamberWindow_Material,
            "ScatteringChamberWindow_Log",
            0,0,0);

    // define Target vacuum window physical volume
    G4cout << G4endl << "###### QweakSimTarget: Define ScatteringChamberWindow_Physical " << G4endl << G4endl;
    ScatteringChamberWindow_Physical   = new G4PVPlacement(0,
            positionScatteringChamberWindow,
            "ScatteringChamberWindow",
            ScatteringChamberWindow_Logical,
            TargetContainer_Physical, //MotherVolume,
            false,
            0,
            pSurfChk);
    G4cout << G4endl << "###### Leaving QweakSimTarget::ConstructScatteringChamberWindow() " << G4endl << G4endl;

//*************************End of Target Vacuum Window *****************************
}

Here is the caller graph for this function:

void QweakSimTarget::ConstructTargetCell ( )

private

Definition at line 321 of file QweakSimTarget.cc.

References positionTarget, pSurfChk, TargetCell_Logical, TargetCell_Material, TargetCell_Physical, targetCellBackInnerRadiusMax, targetCellBackOuterRadiusMax, targetCellDeltaPhi, targetCellFrontInnerRadiusMax, targetCellFrontOuterRadiusMax, targetCellInnerLength, targetCellOuterLength, targetCellStartingPhi, TargetContainer_Logical, and TargetContainer_Physical.

Referenced by ConstructComponent(), and SetTargetLength().

{
    G4cout << G4endl << "###### Calling QweakSimTarget::ConstructTargetCell() " << G4endl << G4endl;

    if (TargetContainer_Logical)
          TargetContainer_Logical->RemoveDaughter(TargetCell_Physical);

    delete TargetCell_Physical;
    TargetCell_Physical = NULL;

    // define target solid volume
    G4cout << G4endl << "###### QweakSimTarget: Define TargetCell_Solid " << G4endl << G4endl;
    // front and back inner radius are at targetCellInnerLength, not targetCellOuterLength
    G4double correction
        = 0.5 * (targetCellOuterLength - targetCellInnerLength)
        * (targetCellBackInnerRadiusMax - targetCellFrontInnerRadiusMax)
        / targetCellInnerLength;
    G4Cons* TargetCell_Solid = new G4Cons("QweakTargetCell_Sol",
                                          targetCellFrontInnerRadiusMax+correction,// G4double  pRmin1,
                                          targetCellFrontOuterRadiusMax, // G4double  pRmax1,
                                          targetCellBackInnerRadiusMax+correction, //G4double  pRmin2,
                                          targetCellBackOuterRadiusMax, //G4double  pRmax2,
                                          0.5*targetCellOuterLength,    //G4double  pDz,
                                          targetCellStartingPhi, //G4double  pSPhi,
                                          targetCellDeltaPhi);   //G4double  pDPhi)

    // define Target logical volume
    G4cout << G4endl << "###### QweakSimTarget: Define TargetCell_Logical " << G4endl << G4endl;

    TargetCell_Logical  = new G4LogicalVolume(TargetCell_Solid,
            TargetCell_Material,
            "QweakTargetCell_Log",
            0,0,0);
    // define Target physical volume
    G4cout << G4endl << "###### QweakSimTarget: Define TargetCell_Physical " << G4endl << G4endl;
    TargetCell_Physical   = new G4PVPlacement(0,
            positionTarget,
            "QweakTargetCell",
            TargetCell_Logical,
            TargetContainer_Physical, //MotherVolume,
            false,
            0,
            pSurfChk);
    G4cout << G4endl << "###### Leaving QweakSimTarget::ConstructTargetCell() " << G4endl << G4endl;
    
}

Here is the caller graph for this function:

void QweakSimTarget::ConstructTargetContainer ( )

private

Definition at line 293 of file QweakSimTarget.cc.

References positionTarget, pSurfChk, ScatteringChamberWindowRadius, targetCellDeltaPhi, targetCellOuterLength, targetCellStartingPhi, TargetContainer_Logical, TargetContainer_Material, TargetContainer_Physical, and theMotherPV.

Referenced by ConstructComponent().

{
    G4cout << G4endl << "###### Calling QweakSimTarget::ConstructTargetContainer() " << G4endl << G4endl;

    G4Tubs* TargetContainer_Solid    = new G4Tubs("TargetContainer_Sol",
            0, //targetCellRadiusMin, jpan@nuclear.uwinnipeg.ca
            ScatteringChamberWindowRadius+0.1*cm,
            0.5*targetCellOuterLength+50.0*cm,
            targetCellStartingPhi,
            targetCellDeltaPhi);

    TargetContainer_Logical  = new G4LogicalVolume(TargetContainer_Solid,
            TargetContainer_Material,
            "TargetContainer_Log",
            0,0,0);

    TargetContainer_Physical   = new G4PVPlacement(0,
            positionTarget,
            "TargetContainer",
            TargetContainer_Logical,
            theMotherPV,
            false,
            0,
            pSurfChk);
    G4cout << G4endl << "###### Leaving QweakSimTarget::ConstructTargetContainer() " << G4endl << G4endl;
    
}

Here is the caller graph for this function:

void QweakSimTarget::ConstructTargetEntranceWindow ( )

private

Definition at line 368 of file QweakSimTarget.cc.

References positionTargetEntranceWindow, pSurfChk, targetCellDeltaPhi, targetCellEntranceWindowThickness, targetCellFrontInnerRadiusMax, targetCellFrontRadiusMin, targetCellStartingPhi, TargetContainer_Logical, TargetContainer_Physical, TargetEntranceWindow_Logical, TargetEntranceWindow_Material, and TargetEntranceWindow_Physical.

Referenced by ConstructComponent(), SetTargetEntranceWindowLength(), and SetTargetLength().

{
    G4cout << G4endl << "###### Calling QweakSimTarget::ConstructTargetEntranceWindow() " << G4endl << G4endl;

    if (TargetContainer_Logical)
          TargetContainer_Logical->RemoveDaughter(TargetEntranceWindow_Physical);

    delete TargetEntranceWindow_Physical;
    TargetEntranceWindow_Physical = NULL;

//--------------------------------------
    // define target window solid volume (front, upstream)
    G4cout << G4endl << "###### QweakSimTarget: Define TargetEntranceWindow_Solid " << G4endl << G4endl;

    G4Tubs* TargetEntranceWindow_Solid    = new G4Tubs("TargetEntranceWindow_Sol",
            targetCellFrontRadiusMin,
            targetCellFrontInnerRadiusMax,
            0.5*targetCellEntranceWindowThickness,
            targetCellStartingPhi,
            targetCellDeltaPhi);

    // define Target window logical volume (front, upstream)
    G4cout << G4endl << "###### QweakSimTarget: Define TargetEntranceWindow_Logical " << G4endl << G4endl;

    TargetEntranceWindow_Logical  = new G4LogicalVolume(TargetEntranceWindow_Solid,
            TargetEntranceWindow_Material,
            "QweakTargetEntranceWindow_Log",
            0,0,0);

    //G4double MaxStepInEntranceWindow = 0.1*targetCellEntranceWindowThickness;
    //TargetEntranceWindow_Logical->SetUserLimits(new G4UserLimits(MaxStepInEntranceWindow));

    // define Target window physical volume (front, upstream)
    G4cout << G4endl << "###### QweakSimTarget: Define TargetEntranceWindow_Physical " << G4endl << G4endl;
    TargetEntranceWindow_Physical   = new G4PVPlacement(0,
            positionTargetEntranceWindow,
            "QweakTargetEntranceWindow",
            TargetEntranceWindow_Logical,
            TargetContainer_Physical, //MotherVolume,
            false,
            0,
            pSurfChk);
    G4cout << G4endl << "###### Leaving QweakSimTarget::ConstructTargetEntranceWindow() " << G4endl << G4endl;

}

Here is the caller graph for this function:

void QweakSimTarget::ConstructTargetExitWindow ( )

private

Definition at line 414 of file QweakSimTarget.cc.

References positionTargetExitWindow, pSurfChk, targetCellBackInnerRadiusMax, targetCellBackRadiusMin, targetCellDeltaPhi, targetCellExitWindowThickness, targetCellStartingPhi, TargetContainer_Logical, TargetContainer_Physical, TargetExitWindow_Logical, TargetExitWindow_Material, and TargetExitWindow_Physical.

Referenced by ConstructComponent(), SetTargetExitWindowLength(), and SetTargetLength().

{
    G4cout << G4endl << "###### Calling QweakSimTarget::ConstructTargetExitWindow() " << G4endl << G4endl;

    if (TargetContainer_Logical)
          TargetContainer_Logical->RemoveDaughter(TargetExitWindow_Physical);

    delete TargetExitWindow_Physical;
    TargetExitWindow_Physical = NULL;

    // define target window solid volume (back, downstream)
    G4cout << G4endl << "###### QweakSimTarget: Define TargetExitWindow_Solid " << G4endl << G4endl;

    G4Tubs* TargetExitWindow_Solid    = new G4Tubs("TargetExitWindow_Sol",
            targetCellBackRadiusMin,
            targetCellBackInnerRadiusMax,
            0.5*targetCellExitWindowThickness,
            targetCellStartingPhi,
            targetCellDeltaPhi);

    // define Target window logical volume (back, downstream)
    G4cout << G4endl << "###### QweakSimTarget: Define TargetExitWindow_Logical " << G4endl << G4endl;

    TargetExitWindow_Logical  = new G4LogicalVolume(TargetExitWindow_Solid,
            TargetExitWindow_Material,
            "QweakTargetExitWindow_Log",
            0,0,0);

    // Set max step size for a certain volume, see
//         http://geant4.web.cern.ch/geant4/G4UsersDocuments/UsersGuides/
//         ForApplicationDeveloper/html/TrackingAndPhysics/thresholdVScut.html
    //G4double MaxStepInExitWindow = 0.1*targetCellExitWindowThickness;
    //TargetExitWindow_Logical->SetUserLimits(new G4UserLimits(MaxStepInExitWindow));


    // define Target window physical volume (back, downstream)
    G4cout << G4endl << "###### QweakSimTarget: Define TargetExitWindow_Physical " << G4endl << G4endl;
    TargetExitWindow_Physical   = new G4PVPlacement(0,
            positionTargetExitWindow,
            "QweakTargetExitWindow",
            TargetExitWindow_Logical,
            TargetContainer_Physical, //MotherVolume,
            false,
            0,
            pSurfChk);
    G4cout << G4endl << "###### Leaving QweakSimTarget::ConstructTargetExitWindow() " << G4endl << G4endl;
    
}

Here is the caller graph for this function:

void QweakSimTarget::ConstructTargetExitWindowNipple ( )

private

Definition at line 463 of file QweakSimTarget.cc.

References positionTargetExitWindow, pSurfChk, targetCellDeltaPhi, targetCellExitWindowNippleRadius, targetCellExitWindowNippleThickness, targetCellStartingPhi, TargetContainer_Logical, TargetContainer_Physical, TargetExitWindowNipple_Logical, TargetExitWindowNipple_Material, and TargetExitWindowNipple_Physical.

Referenced by ConstructComponent(), SetTargetExitWindowNippleLength(), and SetTargetLength().

{
    G4cout << G4endl << "###### Calling QweakSimTarget::ConstructTargetExitWindowNipple() " << G4endl << G4endl;

    if (TargetContainer_Logical)
          TargetContainer_Logical->RemoveDaughter(TargetExitWindowNipple_Physical);

    delete TargetExitWindowNipple_Physical;
    TargetExitWindowNipple_Physical = NULL;   

    // define target window Nipple solid volume (back, downstream)
    G4cout << G4endl << "###### QweakSimTarget: Define TargetExitWindowNipple_Solid " << G4endl << G4endl;

    G4Tubs* TargetExitWindowNipple_Solid    = new G4Tubs("TargetExitWindowNipple_Sol",
            0.,
            targetCellExitWindowNippleRadius,
  //            targetCellExitWindowNippleRadius-0.00001*mm,
            0.5*targetCellExitWindowNippleThickness,
            targetCellStartingPhi,
            targetCellDeltaPhi);

    // define Target window logical volume (back, downstream)
    G4cout << G4endl << "###### QweakSimTarget: Define TargetExitWindowNipple_Logical " << G4endl << G4endl;

    TargetExitWindowNipple_Logical  = new G4LogicalVolume(TargetExitWindowNipple_Solid,
            TargetExitWindowNipple_Material,
            "QweakTargetExitWindowNipple_Log",
            0,0,0);

    //G4double MaxStepInExitWindowNipple = 0.1*targetCellExitWindowNippleThickness;
    //TargetExitWindowNipple_Logical->SetUserLimits(new G4UserLimits(MaxStepInExitWindowNipple));

    // define Target window nipple physical volume (back, downstream)
    G4cout << G4endl << "###### QweakSimTarget: Define TargetExitWindowNipple_Physical " << G4endl << G4endl;
    TargetExitWindowNipple_Physical   = new G4PVPlacement(0,
            positionTargetExitWindow,
            "QweakTargetExitWindowNipple_Physical",
            TargetExitWindowNipple_Logical,
            TargetContainer_Physical, //MotherVolume,
            false,
            0,
            pSurfChk);
    G4cout << G4endl << "###### Leaving QweakSimTarget::ConstructTargetExitWindowNipple() " << G4endl << G4endl;

}

Here is the caller graph for this function:

void QweakSimTarget::ConstructTargetMaterial ( )

private

Definition at line 547 of file QweakSimTarget.cc.

References positionTarget, pSurfChk, targetCellBackInnerRadiusMax, targetCellDeltaPhi, targetCellFrontInnerRadiusMax, targetCellInnerLength, targetCellStartingPhi, TargetContainer_Logical, TargetContainer_Physical, TargetMaterial_Logical, TargetMaterial_Material, and TargetMaterial_Physical.

Referenced by ConstructComponent(), and SetTargetLength().

{
      // define target material solid volume
    G4cout << G4endl << "###### Calling QweakSimTarget::ConstructTargetMaterial() " << G4endl << G4endl;

    if (TargetContainer_Logical)
          TargetContainer_Logical->RemoveDaughter(TargetMaterial_Physical);

    delete TargetMaterial_Physical;
    TargetMaterial_Physical = NULL;

    G4Cons* TargetMaterial_Solid = new G4Cons("QweakTargetMaterial_Sol",
            0., //targetCellFrontRadiusMin, // G4double  pRmin1,
            targetCellFrontInnerRadiusMax, // G4double  pRmax1,
            0., //targetCellBackRadiusMin, //G4double  pRmin2,
            targetCellBackInnerRadiusMax, //G4double  pRmax2,
            0.5*targetCellInnerLength,    //G4double  pDz,
            targetCellStartingPhi, //G4double  pSPhi,
            targetCellDeltaPhi);   //G4double  pDPhi)


    // define Target logical volume
    G4cout << G4endl << "###### QweakSimTarget: Define Target_Logical " << G4endl << G4endl;

    TargetMaterial_Logical  = new G4LogicalVolume(TargetMaterial_Solid,
            TargetMaterial_Material,
            "QweakTargetMaterial_Log",
            0,0,0);

    // set max step size for LH2 target
    //G4double MaxStepInTarget = 0.05*targetCellInnerLength; //step size < 20% of target length
    //TargetMaterial_Logical->SetUserLimits(new G4UserLimits(MaxStepInTarget));

    // define Target material physical volume
    G4cout << G4endl << "###### QweakSimTarget: Define TargetMaterial_Physical " << G4endl << G4endl;
    TargetMaterial_Physical   = new G4PVPlacement(0,
            positionTarget,
            "QweakTargetMaterial",
            TargetMaterial_Logical,
            TargetContainer_Physical, //MotherVolume,
            false,
            0,
            pSurfChk);
    G4cout << G4endl << "###### Leaving QweakSimTarget::ConstructTargetMaterial() " << G4endl << G4endl;

}

Here is the caller graph for this function:

void QweakSimTarget::DestroyComponent

(

)

Definition at line 824 of file QweakSimTarget.cc.

{
}

G4String QweakSimTarget::GetTargetCellMaterial

(

)

G4double QweakSimTarget::GetTargetCenterPositionInZ

(

)

Definition at line 840 of file QweakSimTarget.cc.

References targetZPos.

{
    G4cout << G4endl << "###### Calling QweakSimTarget::GetTargetCenterPositionInZ() " << G4endl << G4endl;
    G4cout << "==== Getting Target CenterPositionZ:  Now the Target Center Position in Z is " << targetZPos/cm << " cm" << G4endl;
    return targetZPos;
}

G4String QweakSimTarget::GetTargetEntranceMaterial

(

)

G4String QweakSimTarget::GetTargetExitWindowMaterial

(

)

G4String QweakSimTarget::GetTargetExitWindowNippleMaterial

(

)

G4double QweakSimTarget::GetTargetLength

(

)

Definition at line 870 of file QweakSimTarget.cc.

References targetCellInnerLength.

{
    G4cout << G4endl << "###### Calling QweakSimTarget::GetTargetLength() " << G4endl << G4endl;
    G4cout << "==== Getting Target Length:  Now the Target Length is " << targetCellInnerLength/cm << " cm" << G4endl;
    return targetCellInnerLength;
}

G4LogicalVolume* QweakSimTarget::getTargetLogicalVolume ( )

inline

Definition at line 92 of file QweakSimTarget.hh.

References TargetMaterial_Logical.

{return TargetMaterial_Logical;} 

G4String QweakSimTarget::GetTargetMaterial

(

)

G4VPhysicalVolume* QweakSimTarget::getTargetPhysicalVolume ( )

inline

Definition at line 93 of file QweakSimTarget.hh.

References TargetMaterial_Physical.

Referenced by QweakSimDetectorConstruction::ConstructQweak().

{return TargetMaterial_Physical;} 

Here is the caller graph for this function:

void QweakSimTarget::SetTarget

(

G4String

targName

)

Definition at line 596 of file QweakSimTarget.cc.

References myUserInfo, SetTargetEntranceWindowLength(), SetTargetEntranceWindowMaterial(), SetTargetExitWindowLength(), SetTargetExitWindowMaterial(), SetTargetExitWindowNippleLength(), SetTargetExitWindowNippleMaterial(), SetTargetLength(), SetTargetMaterial(), targetCellInnerLength, QweakSimUserInformation::TargetThicknessDSALDummy2, QweakSimUserInformation::TargetThicknessDSALDummy4, QweakSimUserInformation::TargetThicknessDSALDummy8, QweakSimUserInformation::TargetThicknessDSCDummy, QweakSimUserInformation::TargetThicknessUSALDummy1, QweakSimUserInformation::TargetThicknessUSALDummy2, QweakSimUserInformation::TargetThicknessUSALDummy4, and QweakSimUserInformation::TargetThicknessUSCDummy.

Referenced by QweakSimTargetMessenger::SetNewValue().

{ // NOTE: Can only make one call at a time.  To change targets you must
  // make a new run, otherwise the positions of the dummy target will be wrong.
  
  double zlength = 0.0; // Use for setting dummy position such that
                        // US face's are aligned (see Target ELOG 18).
  // Technically the US targets are on average 1.39 mm upstream of the entrance window
  // and the DS targets are on average 2.97 mm upstream of the exit window.  For simplicity
  // just align the US faces at -(+)17.173 cm for the US(DS) targets.
  // K.E. Mesick 23/September/2014
  
  if (strcmp(targName,"LH2")==0)
    {
      G4cout << "==== Target is LH2 by default!! " << G4endl;
    }
  else if (strcmp(targName,"USAl1p")==0)
    {
      G4cout << "==== Changing Target to " << targName << G4endl;
      zlength = targetCellInnerLength - 2.0*myUserInfo->TargetThicknessUSALDummy1;
      SetTargetEntranceWindowMaterial("Aluminum");
      SetTargetEntranceWindowLength(myUserInfo->TargetThicknessUSALDummy1);
      SetTargetLength(zlength);
      SetTargetMaterial("Vacuum");
      SetTargetExitWindowMaterial("Vacuum");
      SetTargetExitWindowNippleMaterial("Vacuum");
      G4cout << "==== Changing Target :  Now the Target is " << targName << G4endl;
    }
  else if (strcmp(targName,"USAl2p")==0)
    {
      G4cout << "==== Changing Target to " << targName << G4endl;
      zlength = targetCellInnerLength - 2.0*myUserInfo->TargetThicknessUSALDummy2;
      SetTargetEntranceWindowMaterial("Aluminum");
      SetTargetEntranceWindowLength(myUserInfo->TargetThicknessUSALDummy2);
      SetTargetLength(zlength);
      SetTargetMaterial("Vacuum");
      SetTargetExitWindowMaterial("Vacuum");
      SetTargetExitWindowNippleMaterial("Vacuum");
      G4cout << "==== Changing Target :  Now the Target is " << targName << G4endl;
    }
  else if (strcmp(targName,"USAl4p")==0)
    {
      G4cout << "==== Changing Target to " << targName << G4endl;
      zlength = targetCellInnerLength - 2.0*myUserInfo->TargetThicknessUSALDummy4;
      SetTargetEntranceWindowMaterial("Aluminum");
      SetTargetEntranceWindowLength(myUserInfo->TargetThicknessUSALDummy4);
      SetTargetLength(zlength);
      SetTargetMaterial("Vacuum");
      SetTargetExitWindowMaterial("Vacuum");
      SetTargetExitWindowNippleMaterial("Vacuum");
      G4cout << "==== Changing Target :  Now the Target is " << targName << G4endl;
    }  
  else if (strcmp(targName,"USC")==0)
    {
      G4cout << "==== Changing Target to " << targName << G4endl;
      zlength = targetCellInnerLength - 2.0*myUserInfo->TargetThicknessUSCDummy;
      SetTargetEntranceWindowMaterial("USCarbon");
      SetTargetEntranceWindowLength(myUserInfo->TargetThicknessUSCDummy);
      SetTargetLength(zlength);
      SetTargetMaterial("Vacuum");
      SetTargetExitWindowMaterial("Vacuum");
      SetTargetExitWindowNippleMaterial("Vacuum");
      G4cout << "==== Changing Target :  Now the Target is " << targName << G4endl;
    }
  else if (strcmp(targName,"DSAl2p")==0)
    {
      G4cout << "==== Changing Target to " << targName << G4endl;
      zlength = targetCellInnerLength;
      SetTargetEntranceWindowMaterial("Vacuum");
      SetTargetLength(zlength);
      SetTargetMaterial("Vacuum");
      SetTargetExitWindowLength(myUserInfo->TargetThicknessDSALDummy2);
      SetTargetExitWindowNippleLength(myUserInfo->TargetThicknessDSALDummy2);
      SetTargetExitWindowMaterial("Aluminum"); 
      SetTargetExitWindowNippleMaterial("Aluminum");
      G4cout << "==== Changing Target :  Now the Target is " << targName << G4endl;
    }
  else if (strcmp(targName,"DSAl4p")==0)
    {
      G4cout << "==== Changing Target to " << targName << G4endl;
      zlength = targetCellInnerLength;
      SetTargetEntranceWindowMaterial("Vacuum");
      SetTargetLength(zlength);
      SetTargetMaterial("Vacuum");
      SetTargetExitWindowLength(myUserInfo->TargetThicknessDSALDummy4);
      SetTargetExitWindowNippleLength(myUserInfo->TargetThicknessDSALDummy4);
      //Changed to AlAlloy for alloy simulations. K. Bartlett 2016/9/13
      SetTargetExitWindowMaterial("AlAlloy"); 
      SetTargetExitWindowNippleMaterial("AlAlloy");
      G4cout << "==== Changing Target :  Now the Target is " << targName << G4endl;
    }
  else if (strcmp(targName,"DSAl4pVacuum")==0)
    {
      //Special vacuum DS 4% aluminum target for radiative correction studies.
      G4cout << "==== Changing Target to " << targName << G4endl;
      zlength = targetCellInnerLength;
      SetTargetEntranceWindowMaterial("Vacuum");
      SetTargetLength(zlength);
      SetTargetMaterial("Vacuum");
      SetTargetExitWindowLength(myUserInfo->TargetThicknessDSALDummy4);
      SetTargetExitWindowNippleLength(myUserInfo->TargetThicknessDSALDummy4);
      SetTargetExitWindowMaterial("Vacuum"); 
      SetTargetExitWindowNippleMaterial("Vacuum");
      G4cout << "==== Changing Target :  Now the Target is " << targName << G4endl;
    }
  else if (strcmp(targName,"DSAl8p")==0)
    {
      G4cout << "==== Changing Target to " << targName << G4endl;
      zlength = targetCellInnerLength;
      SetTargetEntranceWindowMaterial("Vacuum");
      SetTargetLength(zlength);
      SetTargetMaterial("Vacuum");
      SetTargetExitWindowLength(myUserInfo->TargetThicknessDSALDummy8);
      SetTargetExitWindowNippleLength(myUserInfo->TargetThicknessDSALDummy8);
      SetTargetExitWindowMaterial("Aluminum"); 
      SetTargetExitWindowNippleMaterial("Aluminum");
      G4cout << "==== Changing Target :  Now the Target is " << targName << G4endl;
    }
  else if (strcmp(targName,"DSC")==0)
    {
      G4cout << "==== Changing Target to " << targName << G4endl;
      zlength = targetCellInnerLength;
      SetTargetEntranceWindowMaterial("Vacuum");
      SetTargetLength(zlength);
      SetTargetMaterial("Vacuum");
      SetTargetExitWindowLength(myUserInfo->TargetThicknessDSCDummy);
      SetTargetExitWindowNippleLength(myUserInfo->TargetThicknessDSCDummy);
      SetTargetExitWindowMaterial("DSCarbon"); 
      SetTargetExitWindowNippleMaterial("DSCarbon");
      G4cout << "==== Changing Target :  Now the Target is " << targName << G4endl;
    }
    else {
        G4cerr << "==== ERROR: Changing Target failed" << G4endl;
    }

}

Here is the call graph for this function:

Here is the caller graph for this function:

void QweakSimTarget::SetTargetCellMaterial

(

G4String

materialName

)

Definition at line 752 of file QweakSimTarget.cc.

References TargetCell_Logical.

Referenced by QweakSimTargetMessenger::SetNewValue().

{
    // search the material by its name
    G4Material* pttoMaterial = G4Material::GetMaterial(materialName);
    if (pttoMaterial)
    {
        G4cout << "==== Changing Target Cell Material: Looking up Material  " << G4endl;
        TargetCell_Logical->SetMaterial(pttoMaterial);
        G4cout << "==== Changing Target Cell Material:  Now the Target Cell is made of " << materialName << G4endl;
    }
    else {
        G4cerr << "==== ERROR: Changing Target Cell Material failed" << G4endl;
    }

}

Here is the caller graph for this function:

void QweakSimTarget::SetTargetCenterPositionInZ

(

G4double

zPos

)

Definition at line 830 of file QweakSimTarget.cc.

References myUserInfo, QweakSimUserInformation::TargetCenterPositionZ, TargetContainer_Physical, and targetZPos.

Referenced by QweakSimTargetMessenger::SetNewValue().

{
    G4cout << G4endl << "###### Calling QweakSimTarget::SetTargetCenterPositionInZ() " << G4endl << G4endl;

    targetZPos = zPos;
    myUserInfo->TargetCenterPositionZ = zPos;
    TargetContainer_Physical->SetTranslation(G4ThreeVector(0.,0.,zPos));
    G4cout << "==== Changing Target CenterPositionZ:  Now the Target Center Position in Z is " << zPos/cm << " cm" << G4endl;
}

Here is the caller graph for this function:

void QweakSimTarget::SetTargetEntranceWindowLength

(

G4double

len

)

Definition at line 878 of file QweakSimTarget.cc.

References CalculateTargetPositions(), ConstructTargetEntranceWindow(), myUserInfo, targetCellEntranceWindowThickness, and QweakSimUserInformation::TargetEntranceWindowThickness.

Referenced by SetTarget().

{
  G4cout << G4endl << "###### Calling QweakSimTarget::SetTargetEntranceWindowLength() " << G4endl << G4endl;
  
  myUserInfo->TargetEntranceWindowThickness = len;
  targetCellEntranceWindowThickness = len;

  G4cout << "==== Changing Target Entrance Window Length:  Now the Target Entrance Window Length is " << len/mm << " mm" << G4endl;

  CalculateTargetPositions();
  ConstructTargetEntranceWindow();

  G4cout << G4endl << "###### Leaving QweakSimTarget::SetTargetEntranceWindowLength() " << G4endl << G4endl;
}

Here is the call graph for this function:

Here is the caller graph for this function:

void QweakSimTarget::SetTargetEntranceWindowMaterial

(

G4String

materialName

)

Definition at line 770 of file QweakSimTarget.cc.

References TargetEntranceWindow_Logical.

Referenced by QweakSimTargetMessenger::SetNewValue(), and SetTarget().

{
    // search the material by its name
    G4Material* pttoMaterial = G4Material::GetMaterial(materialName);
    if (pttoMaterial)
    {
        G4cout << "==== Changing Target Entrance Window Material: Looking up Material  " << G4endl;
        TargetEntranceWindow_Logical->SetMaterial(pttoMaterial);
        G4cout << "==== Changing Target Entrance Window Material:  Now the Target Entrance Window is made of " << materialName << G4endl;
    }
    else {
        G4cerr << "==== ERROR: Changing Target Entrance Window Material failed" << G4endl;
    }

}

Here is the caller graph for this function:

void QweakSimTarget::SetTargetExitWindowLength

(

G4double

len

)

Definition at line 893 of file QweakSimTarget.cc.

References CalculateTargetPositions(), ConstructTargetExitWindow(), myUserInfo, targetCellExitWindowThickness, and QweakSimUserInformation::TargetExitWindowThickness.

Referenced by SetTarget().

{
  G4cout << G4endl << "###### Calling QweakSimTarget::SetTargetExitWindowLength() " << G4endl << G4endl;

  myUserInfo->TargetExitWindowThickness = len;
  targetCellExitWindowThickness = len; 
  G4cout << "==== Changing Target Exit Window Length:  Now the Target Exit Window Length is " << len/mm << " mm" << G4endl;

  CalculateTargetPositions();
  ConstructTargetExitWindow();

  G4cout << G4endl << "###### Leaving QweakSimTarget::SetTargetExitWindowLength() " << G4endl << G4endl;
}

Here is the call graph for this function:

Here is the caller graph for this function:

void QweakSimTarget::SetTargetExitWindowMaterial

(

G4String

materialName

)

Definition at line 788 of file QweakSimTarget.cc.

References TargetExitWindow_Logical.

Referenced by QweakSimTargetMessenger::SetNewValue(), and SetTarget().

{
    // search the material by its name
    G4Material* pttoMaterial = G4Material::GetMaterial(materialName);
    if (pttoMaterial)
    {
        G4cout << "==== Changing Target Exit Window Material: Looking up Material  " << G4endl;
        TargetExitWindow_Logical->SetMaterial(pttoMaterial);
        G4cout << "==== Changing Target Exit Window Material:  Now the Target Exit Window is made of " << materialName << G4endl;
    }
    else {
        G4cerr << "==== ERROR: Changing Target Exit Window Material failed" << G4endl;
    }

}

Here is the caller graph for this function:

void QweakSimTarget::SetTargetExitWindowNippleLength

(

G4double

len

)

Definition at line 907 of file QweakSimTarget.cc.

References CalculateTargetPositions(), ConstructTargetExitWindowNipple(), myUserInfo, targetCellExitWindowNippleThickness, and QweakSimUserInformation::TargetExitWindowNippleThickness.

Referenced by SetTarget().

{
  G4cout << G4endl << "###### Calling QweakSimTarget::SetTargetExitWindowNippleLength() " << G4endl << G4endl;

  myUserInfo->TargetExitWindowNippleThickness = len;
  targetCellExitWindowNippleThickness = len; 
  G4cout << "==== Changing Target Exit Window Nipple Length:  Now the Target Exit Window Nipple Length is " << len/mm << " mm" << G4endl;

  CalculateTargetPositions();
  ConstructTargetExitWindowNipple();

  G4cout << G4endl << "###### Leaving QweakSimTarget::SetTargetExitWindowNippleLength() " << G4endl << G4endl;
}

Here is the call graph for this function:

Here is the caller graph for this function:

void QweakSimTarget::SetTargetExitWindowNippleMaterial

(

G4String

materialName

)

Definition at line 806 of file QweakSimTarget.cc.

References TargetExitWindowNipple_Logical.

Referenced by QweakSimTargetMessenger::SetNewValue(), and SetTarget().

{
    // search the material by its name
    G4Material* pttoMaterial = G4Material::GetMaterial(materialName);
    if (pttoMaterial)
    {
        G4cout << "==== Changing Target Exit Window Nipple Material: Looking up Material  " << G4endl;
        TargetExitWindowNipple_Logical->SetMaterial(pttoMaterial);
        G4cout << "==== Changing Target Exit Window Nipple Material:  Now the Target Exit Window Nipple is made of " << materialName << G4endl;
    }
    else {
        G4cerr << "==== ERROR: Changing Target Exit Window Nipple Material failed" << G4endl;
    }

}

Here is the caller graph for this function:

void QweakSimTarget::SetTargetLength

(

G4double

len

)

Definition at line 848 of file QweakSimTarget.cc.

References CalculateTargetPositions(), ConstructTargetCell(), ConstructTargetEntranceWindow(), ConstructTargetExitWindow(), ConstructTargetExitWindowNipple(), ConstructTargetMaterial(), myUserInfo, targetCellInnerLength, and QweakSimUserInformation::TargetLength.

Referenced by QweakSimTargetMessenger::SetNewValue(), and SetTarget().

{
    G4cout << G4endl << "###### Calling QweakSimTarget::SetTargetLength() " << G4endl << G4endl;

    targetCellInnerLength = len;
    myUserInfo->TargetLength = len;
    G4cout << "==== Changing Target Length:  Now the Target Length is " << len/cm << " cm" << G4endl;

    CalculateTargetPositions();

    // construct target cell/material of appropriate length
    ConstructTargetCell();      // Al cell without end caps
    ConstructTargetMaterial();  // LH2 for production target

    // construct the end caps at the correct location
    ConstructTargetEntranceWindow(); 
    ConstructTargetExitWindow();
    ConstructTargetExitWindowNipple();

    G4cout << G4endl << "###### Leaving QweakSimTarget::SetTargetLength() " << G4endl << G4endl;
}

Here is the call graph for this function:

Here is the caller graph for this function:

void QweakSimTarget::SetTargetMaterial

(

G4String

materialName

)

Definition at line 734 of file QweakSimTarget.cc.

References TargetMaterial_Logical.

Referenced by QweakSimTargetMessenger::SetNewValue(), and SetTarget().

{
    // search the material by its name
    G4Material* pttoMaterial = G4Material::GetMaterial(materialName);
    if (pttoMaterial)
    {
        G4cout << "==== Changing Target Material: Looking up Material  " << G4endl;
        TargetMaterial_Logical->SetMaterial(pttoMaterial);
        G4cout << "==== Changing Target Material:  Now the Target is made of " << materialName << G4endl;
    }
    else {
        G4cerr << "==== ERROR: Changing Target Material failed" << G4endl;
    }

}

Here is the caller graph for this function:

Field Documentation

QweakSimUserInformation* QweakSimTarget::myUserInfo

private

Definition at line 165 of file QweakSimTarget.hh.

Referenced by QweakSimTarget(), SetTarget(), SetTargetCenterPositionInZ(), SetTargetEntranceWindowLength(), SetTargetExitWindowLength(), SetTargetExitWindowNippleLength(), and SetTargetLength().

QweakSimMaterial* QweakSimTarget::pMaterial

private

Definition at line 109 of file QweakSimTarget.hh.

Referenced by QweakSimTarget().

G4ThreeVector QweakSimTarget::positionScatteringChamberWindow

private

Definition at line 170 of file QweakSimTarget.hh.

Referenced by CalculateTargetPositions(), and ConstructScatteringChamberWindow().

G4ThreeVector QweakSimTarget::positionTarget

private

Definition at line 167 of file QweakSimTarget.hh.

Referenced by CalculateTargetPositions(), ConstructTargetCell(), ConstructTargetContainer(), and ConstructTargetMaterial().

G4ThreeVector QweakSimTarget::positionTargetEntranceWindow

private

Definition at line 168 of file QweakSimTarget.hh.

Referenced by CalculateTargetPositions(), and ConstructTargetEntranceWindow().

G4ThreeVector QweakSimTarget::positionTargetExitWindow

private

Definition at line 169 of file QweakSimTarget.hh.

Referenced by CalculateTargetPositions(), ConstructTargetExitWindow(), and ConstructTargetExitWindowNipple().

G4LogicalVolume* QweakSimTarget::ScatteringChamberWindow_Logical

private

Definition at line 127 of file QweakSimTarget.hh.

Referenced by ConstructComponent(), ConstructScatteringChamberWindow(), and QweakSimTarget().

G4Material* QweakSimTarget::ScatteringChamberWindow_Material

private

Definition at line 129 of file QweakSimTarget.hh.

Referenced by ConstructScatteringChamberWindow(), and QweakSimTarget().

G4VPhysicalVolume* QweakSimTarget::ScatteringChamberWindow_Physical

private

Definition at line 128 of file QweakSimTarget.hh.

Referenced by ConstructScatteringChamberWindow(), and QweakSimTarget().

G4double QweakSimTarget::ScatteringChamberWindowRadius

private

Definition at line 157 of file QweakSimTarget.hh.

Referenced by ConstructScatteringChamberWindow(), ConstructTargetContainer(), and QweakSimTarget().

G4double QweakSimTarget::ScatteringChamberWindowThickness

private

Definition at line 158 of file QweakSimTarget.hh.

Referenced by ConstructScatteringChamberWindow(), and QweakSimTarget().

G4LogicalVolume* QweakSimTarget::TargetCell_Logical

private

Definition at line 111 of file QweakSimTarget.hh.

Referenced by ConstructComponent(), ConstructTargetCell(), QweakSimTarget(), and SetTargetCellMaterial().

G4Material* QweakSimTarget::TargetCell_Material

private

Definition at line 113 of file QweakSimTarget.hh.

Referenced by ConstructTargetCell(), and QweakSimTarget().

G4VPhysicalVolume* QweakSimTarget::TargetCell_Physical

private

Definition at line 112 of file QweakSimTarget.hh.

Referenced by ConstructTargetCell(), and QweakSimTarget().

G4double QweakSimTarget::targetCellBackInnerRadiusMax

private

Definition at line 152 of file QweakSimTarget.hh.

Referenced by ConstructTargetCell(), ConstructTargetExitWindow(), ConstructTargetMaterial(), and QweakSimTarget().

G4double QweakSimTarget::targetCellBackOuterRadiusMax

private

Definition at line 153 of file QweakSimTarget.hh.

Referenced by ConstructTargetCell(), and QweakSimTarget().

G4double QweakSimTarget::targetCellBackRadiusMin

private

Definition at line 151 of file QweakSimTarget.hh.

Referenced by ConstructTargetExitWindow(), and QweakSimTarget().

G4double QweakSimTarget::targetCellDeltaPhi

private

Definition at line 161 of file QweakSimTarget.hh.

Referenced by ConstructScatteringChamberWindow(), ConstructTargetCell(), ConstructTargetContainer(), ConstructTargetEntranceWindow(), ConstructTargetExitWindow(), ConstructTargetExitWindowNipple(), ConstructTargetMaterial(), and QweakSimTarget().

G4double QweakSimTarget::targetCellEntranceWindowThickness

private

Definition at line 139 of file QweakSimTarget.hh.

Referenced by CalculateTargetPositions(), ConstructTargetEntranceWindow(), QweakSimTarget(), and SetTargetEntranceWindowLength().

G4double QweakSimTarget::targetCellExitWindowNippleRadius

private

Definition at line 155 of file QweakSimTarget.hh.

Referenced by ConstructTargetExitWindowNipple(), and QweakSimTarget().

G4double QweakSimTarget::targetCellExitWindowNippleThickness

private

Definition at line 141 of file QweakSimTarget.hh.

Referenced by ConstructTargetExitWindowNipple(), QweakSimTarget(), and SetTargetExitWindowNippleLength().

G4double QweakSimTarget::targetCellExitWindowThickness

private

Definition at line 140 of file QweakSimTarget.hh.

Referenced by CalculateTargetPositions(), ConstructTargetExitWindow(), QweakSimTarget(), and SetTargetExitWindowLength().

G4double QweakSimTarget::targetCellFrontInnerRadiusMax

private

Definition at line 148 of file QweakSimTarget.hh.

Referenced by ConstructTargetCell(), ConstructTargetEntranceWindow(), ConstructTargetMaterial(), and QweakSimTarget().

G4double QweakSimTarget::targetCellFrontOuterRadiusMax

private

Definition at line 149 of file QweakSimTarget.hh.

Referenced by ConstructTargetCell(), and QweakSimTarget().

G4double QweakSimTarget::targetCellFrontRadiusMin

private

Definition at line 147 of file QweakSimTarget.hh.

Referenced by ConstructTargetEntranceWindow(), and QweakSimTarget().

G4double QweakSimTarget::targetCellInnerLength

private

Definition at line 144 of file QweakSimTarget.hh.

Referenced by CalculateTargetPositions(), ConstructTargetCell(), ConstructTargetMaterial(), GetTargetLength(), QweakSimTarget(), SetTarget(), and SetTargetLength().

G4double QweakSimTarget::targetCellOuterLength

private

Definition at line 145 of file QweakSimTarget.hh.

Referenced by CalculateTargetPositions(), ConstructTargetCell(), ConstructTargetContainer(), and QweakSimTarget().

G4double QweakSimTarget::targetCellStartingPhi

private

Definition at line 160 of file QweakSimTarget.hh.

Referenced by ConstructScatteringChamberWindow(), ConstructTargetCell(), ConstructTargetContainer(), ConstructTargetEntranceWindow(), ConstructTargetExitWindow(), ConstructTargetExitWindowNipple(), ConstructTargetMaterial(), and QweakSimTarget().

G4double QweakSimTarget::targetCellWallThickness

private

Definition at line 142 of file QweakSimTarget.hh.

Referenced by QweakSimTarget().

G4LogicalVolume* QweakSimTarget::TargetContainer_Logical

private

Definition at line 135 of file QweakSimTarget.hh.

Referenced by ConstructComponent(), ConstructTargetCell(), ConstructTargetContainer(), ConstructTargetEntranceWindow(), ConstructTargetExitWindow(), ConstructTargetExitWindowNipple(), ConstructTargetMaterial(), and QweakSimTarget().

G4Material* QweakSimTarget::TargetContainer_Material

private

Definition at line 137 of file QweakSimTarget.hh.

Referenced by ConstructTargetContainer(), and QweakSimTarget().

G4VPhysicalVolume* QweakSimTarget::TargetContainer_Physical

private

Definition at line 136 of file QweakSimTarget.hh.

Referenced by CalculateTargetPositions(), ConstructScatteringChamberWindow(), ConstructTargetCell(), ConstructTargetContainer(), ConstructTargetEntranceWindow(), ConstructTargetExitWindow(), ConstructTargetExitWindowNipple(), ConstructTargetMaterial(), QweakSimTarget(), and SetTargetCenterPositionInZ().

G4LogicalVolume* QweakSimTarget::TargetEntranceWindow_Logical

private

Definition at line 115 of file QweakSimTarget.hh.

Referenced by ConstructComponent(), ConstructTargetEntranceWindow(), QweakSimTarget(), and SetTargetEntranceWindowMaterial().

G4Material* QweakSimTarget::TargetEntranceWindow_Material

private

Definition at line 117 of file QweakSimTarget.hh.

Referenced by ConstructTargetEntranceWindow(), and QweakSimTarget().

G4VPhysicalVolume* QweakSimTarget::TargetEntranceWindow_Physical

private

Definition at line 116 of file QweakSimTarget.hh.

Referenced by ConstructTargetEntranceWindow(), and QweakSimTarget().

G4LogicalVolume* QweakSimTarget::TargetExitWindow_Logical

private

Definition at line 119 of file QweakSimTarget.hh.

Referenced by ConstructComponent(), ConstructTargetExitWindow(), QweakSimTarget(), and SetTargetExitWindowMaterial().

G4Material* QweakSimTarget::TargetExitWindow_Material

private

Definition at line 121 of file QweakSimTarget.hh.

Referenced by ConstructTargetExitWindow(), and QweakSimTarget().

G4VPhysicalVolume* QweakSimTarget::TargetExitWindow_Physical

private

Definition at line 120 of file QweakSimTarget.hh.

Referenced by ConstructTargetExitWindow(), and QweakSimTarget().

G4LogicalVolume* QweakSimTarget::TargetExitWindowNipple_Logical

private

Definition at line 123 of file QweakSimTarget.hh.

Referenced by ConstructTargetExitWindowNipple(), QweakSimTarget(), and SetTargetExitWindowNippleMaterial().

G4Material* QweakSimTarget::TargetExitWindowNipple_Material

private

Definition at line 125 of file QweakSimTarget.hh.

Referenced by ConstructTargetExitWindowNipple(), and QweakSimTarget().

G4VPhysicalVolume* QweakSimTarget::TargetExitWindowNipple_Physical

private

Definition at line 124 of file QweakSimTarget.hh.

Referenced by ConstructTargetExitWindowNipple(), and QweakSimTarget().

G4LogicalVolume* QweakSimTarget::TargetMaterial_Logical

private

Definition at line 131 of file QweakSimTarget.hh.

Referenced by ConstructComponent(), ConstructTargetMaterial(), getTargetLogicalVolume(), QweakSimTarget(), and SetTargetMaterial().

G4Material* QweakSimTarget::TargetMaterial_Material

private

Definition at line 133 of file QweakSimTarget.hh.

Referenced by ConstructTargetMaterial(), and QweakSimTarget().

G4VPhysicalVolume* QweakSimTarget::TargetMaterial_Physical

private

Definition at line 132 of file QweakSimTarget.hh.

Referenced by ConstructTargetMaterial(), getTargetPhysicalVolume(), and QweakSimTarget().

QweakSimTargetMessenger* QweakSimTarget::targetMessenger

private

Definition at line 164 of file QweakSimTarget.hh.

Referenced by QweakSimTarget(), and ~QweakSimTarget().

G4VSensitiveDetector* QweakSimTarget::TargetSD

private

Definition at line 172 of file QweakSimTarget.hh.

Referenced by ConstructComponent(), and QweakSimTarget().

G4double QweakSimTarget::targetZPos

private

Definition at line 162 of file QweakSimTarget.hh.

Referenced by GetTargetCenterPositionInZ(), QweakSimTarget(), and SetTargetCenterPositionInZ().

G4VPhysicalVolume* QweakSimTarget::theMotherPV

private

Definition at line 107 of file QweakSimTarget.hh.

Referenced by ConstructComponent(), ConstructTargetContainer(), and QweakSimTarget().

---

The documentation for this class was generated from the following files:

• QweakSimTarget.hh
• QweakSimTarget.cc