MollerTarget.cc

Go to the documentation of this file.

/**

   \file MollerTarget.cc
   \author Peiqing Wang
   \date Tue Mar 13 01:52:47 EDT 2012
*/

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

#include "MollerTarget.hh"

// user classes
class MollerMaterial;
class MollerTargetMessenger;

// geant4 includes
#include "G4UserLimits.hh"

// user includes
#include "MollerTargetMessenger.hh"

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

MollerTarget::MollerTarget()
{
    G4cout << G4endl << "###### Calling MollerTarget::MollerTarget() " << G4endl << G4endl;

    //myUserInfo = myUI;

    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;

    TargetMaterial_Logical  = NULL;
    TargetMaterial_Physical = NULL;
    Target_Material         = NULL;


    // definition of a mil = inch/1000
    static const G4double mil = 0.001*2.54*cm;
    static const G4double inch = 2.54*cm;

     pMaterial = new MollerMaterial();
     pMaterial->DefineMaterials();

    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");
    Target_Material                 = pMaterial->GetMaterial("H2Liquid");

    // define target geometry values

    targetCellEntranceWindowThickness =  5*mil;;
    targetCellExitWindowThickness =  20*mil;
    targetCellExitWindowNippleThickness =  5*mil;
    targetCellWallThickness   = 10.0*mil;

    targetCellExitWindowNippleRadius = 0.5*0.622*inch;

    targetCellInnerLength     = 34.346*cm;
    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        = 0.0*cm;
    //targetZPos        = myUserInfo->TargetCenterPositionZ; //-650.0*cm;

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

    //myUserInfo->TargetLength = targetCellInnerLength;
    //myUserInfo->TargetEntranceWindowThickness = targetCellEntranceWindowThickness;
    //myUserInfo->TargetExitWindowThickness = targetCellExitWindowThickness;
    targetMessenger = new MollerTargetMessenger(this);

    fTargetLength  = 150.0*cm;  

}

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
MollerTarget::~MollerTarget()
{
    delete pMaterial;
    delete targetMessenger;
}

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

void MollerTarget::ConstructComponent(G4VPhysicalVolume* MotherVolume)
{
    G4cout << G4endl << "###### Calling MollerTarget::ConstructComponent() " << G4endl << G4endl;
   
      //------Target-----//
  G4double inner_rad = 0.0*cm;
  G4double start_angle = 0.0*deg;
  G4double end_angle = 360.0*deg;

  G4double outer_rad_targ = 4.0*cm;
  G4double half_len_targ = fTargetLength/2;

  G4cout<< "target half-length = "<<half_len_targ/cm<<G4endl;

  G4Tubs* target_solid = new G4Tubs("tubeTarget", inner_rad, outer_rad_targ,
           half_len_targ, start_angle, end_angle);

  G4LogicalVolume* target_logic = new G4LogicalVolume(target_solid, Target_Material, "logicTarget");

      G4VPhysicalVolume* target_phys   = new G4PVPlacement(0,
            G4ThreeVector(0, 0, targetZPos),
            "Target",
            target_logic,
            MotherVolume,
            false,
            0);

  G4VisAttributes* targetVisAtt= new G4VisAttributes(G4Colour(1.0,1.0,0.0));
  targetVisAtt->SetVisibility(true);
  target_logic->SetVisAttributes(targetVisAtt);
  
    // define target entrance window (front, upstream)
    G4Tubs* TargetEntranceWindow_Solid    = new G4Tubs("TargetEntranceWindow_Sol",
             inner_rad, outer_rad_targ,
            0.5*targetCellEntranceWindowThickness,
            start_angle, end_angle);

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

   TargetEntranceWindow_Physical   = new G4PVPlacement(0,
            G4ThreeVector(0, 0, targetZPos-(half_len_targ+0.5*targetCellEntranceWindowThickness)),// positionTargetEntranceWindow,
            "MollerTargetEntranceWindow",
            TargetEntranceWindow_Logical,
            MotherVolume,
            false,
            0);

    // define target exit window (back, downstream)
    G4Tubs* TargetExitWindow_Solid    = new G4Tubs("TargetExitWindow_Sol",
             inner_rad, outer_rad_targ,
            0.5*targetCellExitWindowThickness,
            start_angle, end_angle);

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

   TargetExitWindow_Physical   = new G4PVPlacement(0,
            G4ThreeVector(0, 0, targetZPos+(half_len_targ+0.5*targetCellExitWindowThickness)),// positionTargetEntranceWindow,
            "MollerTargetExiteWindow",
            TargetExitWindow_Logical,
            MotherVolume,
            false,
            0);
   
  G4Colour  blue      (   0/255.,   0/255., 255/255.);
  G4VisAttributes* targetWindowVisAtt= new G4VisAttributes(blue);
  targetWindowVisAtt->SetVisibility(true);
  TargetEntranceWindow_Logical->SetVisAttributes(targetWindowVisAtt);
  TargetExitWindow_Logical->SetVisAttributes(targetWindowVisAtt);


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

} // end of  MollerTarget::ConstructComponent()

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

void MollerTarget::SetTargetMaterial(G4String materialName)
{
    // 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;
    }

}

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

void MollerTarget::SetTargetCellMaterial(G4String materialName)
{
    // 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;
    }

}

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

void MollerTarget::SetTargetEntranceWindowMaterial(G4String materialName)
{
    // 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;
    }

}

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

void MollerTarget::SetTargetExitWindowMaterial(G4String materialName)
{
    // 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;
    }

}

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

void MollerTarget::SetTargetExitWindowNippleMaterial(G4String materialName)
{
    // 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;
    }

}

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

void MollerTarget::DestroyComponent()
{
}

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

void MollerTarget::SetTargetCenterPositionInZ(G4double zPos)
{
    G4cout << G4endl << "###### Calling MollerTarget::SetTargetCenterPositionInZ() " << G4endl << G4endl;

    targetZPos = zPos;
    //myUserInfo->TargetCenterPositionZ = zPos;
    //TargetCell_Physical->SetTranslation(G4ThreeVector(0.,0., zPos));
    TargetContainer_Physical->SetTranslation(G4ThreeVector(0.,0., zPos));
}

G4double MollerTarget::GetTargetCenterPositionInZ()
{
    G4cout << G4endl << "###### Calling MollerTarget::GetTargetCenterPositionInZ() " << G4endl << G4endl;
    return targetZPos;
}

void MollerTarget::SetTargetLength(G4double length)
{
    G4cout << G4endl << "###### Calling MollerTarget::SetTargetLength() " << G4endl << G4endl;
    fTargetLength = length;
}
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......