WCSimDetectorMessenger.cc

Go to the documentation of this file.

#include "WCSimDetectorMessenger.hh"

#include "WCSimDetectorConstruction.hh"
#include "G4UIdirectory.hh"
#include "G4UIcommand.hh"
#include "G4UIparameter.hh"
#include "G4UIcmdWithAString.hh"
#include "G4UIcmdWithADoubleAndUnit.hh"

WCSimDetectorMessenger::WCSimDetectorMessenger(WCSimDetectorConstruction* WCSimDet)
:WCSimDetector(WCSimDet)
{ 
  WCSimDir = new G4UIdirectory("/WCSim/");
  WCSimDir->SetGuidance("Commands to change the geometry of the simulation");


  PMTConfig = new G4UIcmdWithAString("/WCSim/WCgeom",this);
  PMTConfig->SetGuidance("Set the geometry configuration for the WC.");
  PMTConfig->SetGuidance("Available options are:\n"
                          "SuperK\n"
                          "SuperK_20inchPMT_20perCent\n"
                          "SuperK_20inchBandL_20perCent\n"
                          "SuperK_12inchBandL_15perCent\n"
                          "SuperK_20inchBandL_14perCent\n"
                          "Cylinder_60x74_20inchBandL_14perCent\n"
                          "Cylinder_60x74_20inchBandL_40perCent\n"
                          "Cylinder_12inchHPD_15perCent\n"
                          "HyperK\n"
                          "HyperK_20perCent\n"
                          "HyperKWithOD\n"
                          "EggShapedHyperK\n"
                          "EggShapedHyperK_withHPD\n"
                          );
  PMTConfig->SetParameterName("PMTConfig", false);
  PMTConfig->SetCandidates("SuperK "
                          "SuperK_20inchPMT_20perCent "
                          "SuperK_20inchBandL_20perCent "
                          "SuperK_12inchBandL_15perCent "
                          "SuperK_20inchBandL_14perCent "
                          "Cylinder_60x74_20inchBandL_14perCent "
                          "Cylinder_60x74_20inchBandL_40perCent "
                          "Cylinder_12inchHPD_15perCent "
                          "HyperK "
                          "HyperK_20perCent "
                          "HyperKWithOD "
                          "EggShapedHyperK "
                          "EggShapedHyperK_withHPD "
                          );
  PMTConfig->AvailableForStates(G4State_PreInit, G4State_Idle);

  WCVisChoice = new G4UIcmdWithAString("/WCSim/WCVisChoice",this);
  WCVisChoice->SetGuidance("Set the visualization style for the WC.");
  WCVisChoice->SetGuidance("Available options are:\n"
                          "OGLSX\n"
                          "RayTracer\n"
                          );
  WCVisChoice->SetParameterName("WCVisChoice", false);
  WCVisChoice->SetCandidates("OGLSX "
                          "RayTracer "
                          );
  WCVisChoice->AvailableForStates(G4State_PreInit, G4State_Idle);


  PMTSize = new G4UIcmdWithAString("/WCSim/WCPMTsize",this);
  PMTSize->SetGuidance("Set alternate PMT size for the WC (Must be entered after geometry details are set).");
  PMTSize->SetGuidance("Available options 20inch 10inch");
  PMTSize->SetParameterName("PMTSize", false);
  PMTSize->SetCandidates("20inch 10inch");
  PMTSize->AvailableForStates(G4State_PreInit, G4State_Idle);

  SavePi0 = new G4UIcmdWithAString("/WCSim/SavePi0", this);
  SavePi0->SetGuidance("true or false");
  SavePi0->SetParameterName("SavePi0",false);
  SavePi0->SetCandidates("true false");
  SavePi0->AvailableForStates(G4State_PreInit, G4State_Idle);
  
  
  PMTQEMethod = new G4UIcmdWithAString("/WCSim/PMTQEMethod", this);
  PMTQEMethod->SetGuidance("Set the PMT configuration.");
  PMTQEMethod->SetGuidance("Available options are:\n"
                             "Stacking_Only\n"
                             "Stacking_And_SensitiveDetector\n"
                             "SensitiveDetector_Only\n"
                           "DoNotApplyQE");
  PMTQEMethod->SetParameterName("PMTQEMethod", false);
  PMTQEMethod->SetCandidates("Stacking_Only "
                             "Stacking_And_SensitiveDetector "
                             "SensitiveDetector_Only "
                             "DoNotApplyQE");
  PMTQEMethod->AvailableForStates(G4State_PreInit, G4State_Idle);

  PMTCollEff = new G4UIcmdWithAString("/WCSim/PMTCollEff", this);
  PMTCollEff->SetGuidance("Set the PMT configuration.");
  PMTCollEff->SetGuidance("Available options are:\n"
                          "on\n"
                          "off\n");
  PMTCollEff->SetParameterName("PMTCollEff", false);
  PMTCollEff->SetCandidates("on "
                            "off ");
  PMTCollEff->AvailableForStates(G4State_PreInit, G4State_Idle);

  waterTank_Length = new G4UIcmdWithADoubleAndUnit("/WCSim/EggShapedHyperK/waterTank_Length", this);
  waterTank_Length->SetGuidance("Set the Length of the egg-shaped HyperK detector (unit: mm cm m).");
  waterTank_Length->SetParameterName("waterTank_length", true);
  waterTank_Length->SetDefaultValue(49500.);
  waterTank_Length->SetUnitCategory("Length");
  waterTank_Length->SetDefaultUnit("mm");
  waterTank_Length->SetUnitCandidates("mm cm m");

  LCConfig = new G4UIcmdWithAString("/WCSim/LCConfig",this);
  LCConfig->SetGuidance("Set the geometry configuration for the light collecting mirror.");
  LCConfig->SetGuidance("**For 20 inch PMT Only**");
  LCConfig->SetGuidance("Available options are:\n"
                          "No_Mirror\n"
                          "Mirror_OldLC\n"
                          "Mirror_2018Oct\n"
                          );
  LCConfig->SetParameterName("LCConfig", true);
  LCConfig->SetCandidates("No_Mirror "
                          "Mirror_OldLC "
                          "Mirror_2018Oct "
                          );
  LCConfig->SetDefaultValue("No_Mirror");
  LCConfig->AvailableForStates(G4State_PreInit, G4State_Idle);

  WCConstruct = new G4UIcmdWithoutParameter("/WCSim/Construct", this);
  WCConstruct->SetGuidance("Update detector construction with new settings.");
}

WCSimDetectorMessenger::~WCSimDetectorMessenger()
{
  delete PMTConfig;
  delete LCConfig;
  delete SavePi0;
  delete PMTQEMethod;
  delete PMTCollEff;
  delete waterTank_Length;
  delete WCVisChoice;

  delete tubeCmd;
  delete distortionCmd;
  delete WCSimDir;
}

void WCSimDetectorMessenger::SetNewValue(G4UIcommand* command,G4String newValue)
{    
        if( command == PMTConfig ) { 
                WCSimDetector->SetIsUpright(false);
                WCSimDetector->SetIsEggShapedHyperK(false);
                if ( newValue == "SuperK") {
                  WCSimDetector->SetSuperKGeometry();
                } else if ( newValue == "SuperK_20inchPMT_20perCent" ){
                  WCSimDetector->SuperK_20inchPMT_20perCent();
                } else if ( newValue == "SuperK_20inchBandL_20perCent" ){
                  WCSimDetector->SuperK_20inchBandL_20perCent();
                } else if ( newValue == "SuperK_12inchBandL_15perCent" ) {
                  WCSimDetector->SuperK_12inchBandL_15perCent();
                } else if ( newValue == "SuperK_20inchBandL_14perCent" ) {
                  WCSimDetector->SuperK_20inchBandL_14perCent();
                } else if ( newValue == "Cylinder_60x74_20inchBandL_14perCent" ) {
                  WCSimDetector->Cylinder_60x74_20inchBandL_14perCent();
                } else if ( newValue == "Cylinder_60x74_20inchBandL_40perCent" ) {
                  WCSimDetector->Cylinder_60x74_20inchBandL_40perCent();
                } else if ( newValue == "Cylinder_12inchHPD_15perCent" ){
                  WCSimDetector->Cylinder_12inchHPD_15perCent();
                } else if ( newValue == "HyperK" ){
                  WCSimDetector->SetHyperKGeometry();
                } else if ( newValue == "HyperK_20perCent" ){
                  WCSimDetector->SetHyperKGeometry_20perCent();
                } else if ( newValue == "HyperKWithOD" ){
                  WCSimDetector->SetHyperKWithODGeometry();
                } else if ( newValue == "EggShapedHyperK") {
                  WCSimDetector->SetIsEggShapedHyperK(true);
                  WCSimDetector->SetEggShapedHyperKGeometry();
                } else if ( newValue == "EggShapedHyperK_withHPD") {
                  WCSimDetector->SetIsEggShapedHyperK(true);
                  WCSimDetector->SetEggShapedHyperKGeometry_withHPD();
                } else
                  G4cout << "That geometry choice is not defined!" << G4endl;
        }
  
        if (command == SavePi0){
          G4cout << "Set the flag for saving pi0 info " << newValue << G4endl;
          if (newValue=="true"){
            WCSimDetector->SavePi0Info(true);
          }else if (newValue == "false"){
            WCSimDetector->SavePi0Info(false);
          }else{
            
          }
        }

        if (command == PMTQEMethod){
          G4cout << "Set PMT QE Method " << newValue << " ";
          if (newValue == "Stacking_Only"){
            WCSimDetector->SetPMT_QE_Method(1);
            G4cout << "1";
          }else if (newValue == "Stacking_And_SensitiveDetector"){
            WCSimDetector->SetPMT_QE_Method(2);
            G4cout << "2";
          }else if (newValue == "SensitiveDetector_Only"){
            WCSimDetector->SetPMT_QE_Method(3);
            G4cout << "3";
          }else if (newValue == "DoNotApplyQE"){
            WCSimDetector->SetPMT_QE_Method(4);
            G4cout << "4";
          }else{
            
          }
          G4cout << G4endl;
        }

        if (command == PMTCollEff){
          G4cout << "Set PMT Collection Efficiency " << newValue << " ";
          if (newValue == "on"){
            WCSimDetector->SetPMT_Coll_Eff(1);
            G4cout << "1";
          }else if (newValue == "off"){
            WCSimDetector->SetPMT_Coll_Eff(0);
            G4cout << "0";
          }
          G4cout << G4endl;
        }
        
        if (command == WCVisChoice){
          G4cout << "Set Vis Choice " << newValue << " ";
          if (newValue == "OGLSX"){
            WCSimDetector->SetVis_Choice("OGLSX");
          }else if (newValue == "RayTracer"){
            WCSimDetector->SetVis_Choice("RayTracer");
          }
          G4cout << G4endl;
        }

        if (command == waterTank_Length){
          bool isEggShapedHyperK = WCSimDetector->GetIsEggShapedHyperK();
          if(isEggShapedHyperK == true){
            G4cout << "Set Partition Length in a cylinder " << newValue << " " << G4endl;
            WCSimDetector->SetwaterTank_Length(waterTank_Length->GetNewDoubleValue(newValue));
            WCSimDetector->SetWaterTubeLength(waterTank_Length->GetNewDoubleValue(newValue));
            // WCSimDetector->SetIsEggShapedHyperK(true);
          }
          else {
            G4cout << "Not egg-shaped HyperK Geometry. Detector length unchanged." << G4endl;
          }
        }


        if(command == PMTSize) {
                G4cout << "SET PMT SIZE" << G4endl;
                if ( newValue == "20inch"){
//                      WCSimDetector->Set20inchPMTs();
                }else if (newValue == "10inch"){
//                      WCSimDetector->Set10inchPMTs();
                }else
                        G4cout << "That PMT size is not defined!" << G4endl;    
        }

        if( command == LCConfig ) { 
                // LC Type is defined in WCSimDetectorConstruction.hh
                if ( newValue == "No_Mirror") {
                  WCSimDetector->SetLCType(0);
                } else if ( newValue == "Mirror_OldLC" ){
                  WCSimDetector->SetLCType(1);
                } else if ( newValue == "Mirror_2018Oct" ){
                  WCSimDetector->SetLCType(2);
                }
        }

        if(command == WCConstruct) {
                WCSimDetector->UpdateGeometry();
        }

}