CDCJobCntlParModifierModule.cc

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2012 - Belle II Collaboration                             *
 *                                                                        *
 * Author: The Belle II Collaboration                                     *
 * Contributors: CDC group                                    *
 *                                                                        *
 * This software is provided "as is" without any warranty.                *
 **************************************************************************/
 
#include <cdc/modules/cdcJobCntlParModifier/CDCJobCntlParModifierModule.h>
//#include <framework/logging/Logger.h>
 
using namespace std;
using namespace Belle2;
using namespace CDC;
 
// register module
REG_MODULE(CDCJobCntlParModifier)
CDCJobCntlParModifierModule::CDCJobCntlParModifierModule() : Module(), m_scp(CDCSimControlPar::getInstance()),
  m_gcp(CDCGeoControlPar::getInstance()), m_timeWalk(), m_wireSag(), m_modLeftRightFlag(), m_debug4Sim(), m_thresholdEnergyDeposit(),
  m_minTrackLength(), m_maxSpaceResol(), m_addFudgeFactorForSigmaForData(), m_addFudgeFactorForSigmaForMC(),
  m_mapperGeometry(), m_mapperPhiAngle(), m_debug4Geo(), m_printMaterialTable(),
  m_materialDefinitionMode(), m_senseWireZposMode(),
  m_displacement(),
  m_alignment(),
  m_misalignment(),
  m_displacementInputType(), m_alignmentInputType(), m_misalignmentInputType(), m_xtInputType(), m_sigmaInputType(),
  m_propSpeedInputType(), m_t0InputType(), m_twInputType(), m_bwInputType(), m_chMapInputType(), m_eDepToADCInputType(),
  m_displacementFile(),
  m_alignmentFile(), m_misalignmentFile(), m_xtFile(), m_sigmaFile(), m_propSpeedFile(), m_t0File(), m_twFile(), m_bwFile(),
  m_chMapFile(), m_eDepToADCFile()
 
{
  //  B2INFO("CDCJobCntlParModifierModule::constructor called.");
  // Set description
  setDescription("Change job contorol parameters. Please put this module in the path (before Geometry module) with specified input parameters when you want to change them.");
  setPropertyFlags(c_ParallelProcessingCertified);
 
  //N.B. The following default values must be identical to the ones in xxControlPar objects.
  //For Simulation
  //Switch for debug
  addParam("Debug4Sim", m_debug4Sim, "Switch on/off debug in FullSim.", false);
  //Switch for time walk in translator
  addParam("TimeWalkInTranslator", m_timeWalk,
           "Switch on/off time-walk correction in the TDC-count translator for MC events. Use the CDCDigitizer's flag if you want to add/not add time walk in digitization.",
           true);
  //Switch for wire sag
  addParam("WireSag", m_wireSag,
           "Switch on/off sense wire (gravitational) sag in FullSim. Here, sag means the main part which corresponds to design+displacement in case of wire position. You can control the perturbative part (corresponting to (mis)alignment in case of wire-position) of sag in Digitizer.",
           true);
  //Switch for modified left/right flag
  addParam("ModLeftRightFlag", m_modLeftRightFlag, "Switch on/off calculation of modified left/right flag in FullSim.", false);
  //energy thresh
  addParam("ThresholdEnergyDeposit",  m_thresholdEnergyDeposit,
           "Energy deposite (edep) thresh. for G4 step (GeV). All hits with smaller edep will be dropped at FullSim level. Set this to a negative value if you want to keep simhits with edep=0.",
           0.0);
  //min. track length
  addParam("MinTrackLength",  m_minTrackLength,
           "Minimum track length (cm) required for saving in MCParticle.; to be applied to all particles other than primaries and decay-daughters.",
           15.0);
 
  //For Geometry
  //Switch for debug
  addParam("Debug4Geo", m_debug4Geo, "Switch on/off debug in Geo.", false);
  //Switch for printing material table
  addParam("PrintMaterialTable", m_printMaterialTable,
           "Switch on/off printing the G4 material table at the stage of CDC geometry creation.", false);
  //material definition mode
  addParam("MaterialDefinitionMode",  m_materialDefinitionMode,
           "Material definition mode: =0: define a mixture of gases and wires in the entire tracking volume; =1: dummy; =2: define all sense and field wires explicitly in the volume.",
           0);
  //sense wire z pos mode
  addParam("SenseWireZposMode",  m_senseWireZposMode,
           "Sense wire z position mode: =1: correct for feedthrough effect; =0: not correct.",  1);
  //displacement switch
  addParam("Displacement", m_displacement, "Switch for wire displacement: on/off.",  true);
  //alignmentt switch
  addParam("Alignment", m_alignment, "Switch for wire alignment: on/off.",  true);
  //misalignment switch
  //  addParam("Misalignment", m_misalignment, "Switch for wire misalignment: on/off.",  true);
  m_misalignment = false;
 
 
  //input type for displacement
  addParam("DisplacementInputType", m_displacementInputType, "Input type for wire displacement; db-object (true); text-file (false).",
           true);
  //input type for alignment
  addParam("AlignmentInputType", m_alignmentInputType, "Input type for wire alignment; db-object (true); text-file (false).", true);
  //input type for misalignment
  //  addParam("MisalignmentInputType", m_misalignmentInputType, "Input type for wire misalignment; db-object (true); text-file (false).", true);
  m_misalignmentInputType = true;
  //input type for xt-relation
  addParam("XtInputType", m_xtInputType, "Input type for xt-relations; db-object (true); text-file (false).", true);
  //input type for sigma
  addParam("SigmaInputType", m_sigmaInputType, "Input type for sigmas; db-object (true); text-file (false).", true);
  //input type for prop-speed
  addParam("PropSpeedInputType", m_propSpeedInputType, "Input type for prop-speeds; db-object (true); text-file (false).", true);
  //input type for t0
  addParam("T0InputType", m_t0InputType, "Input type for t0s; db-object (true); text-file (false).", true);
  //input type for time walk
  addParam("TimeWalkInputType", m_twInputType, "Input type for time walks; db-object (true); text-file (false).", true);
  //input type for bad wire
  addParam("BadWireInputType", m_bwInputType, "Input type for bad wires; db-object (true); text-file (false).", true);
  //input type for channel map
  addParam("ChannelMapInputType", m_chMapInputType, "Input type for channel map; db-object (true); text-file (false).", true);
  //input type for edep-to-adc
  addParam("EDepToADCInputType", m_eDepToADCInputType, "Input type for edep-to-adc; db-object (true); text-file (false).", true);
 
 
  //displacement file
  addParam("DisplacementFile", m_displacementFile, "Input file name (on cdc/data) for wire displacement.",
           string("displacement_v2.2.1.dat"));
  //alignment file
  addParam("AlignmentFile", m_alignmentFile, "Input file name (on cdc/data) for wire alignment.",  string("alignment_v2.dat"));
  //misalignment file
  //  addParam("MisalignmentFile", m_misalignmentFile, "Input file name (on cdc/data) for wire misalignment.", string("misalignment_v2.dat"));
  m_misalignmentFile = string("misalignment_v2.dat");
  //xt-relation
  addParam("XtFile", m_xtFile, "Input file name (on cdc/data) for xt-relations. You can specify either an uncompressed or gzip file.",
           string("xt_v3_chebyshev.dat.gz"));
  //sigma
  addParam("SigmaFile", m_sigmaFile, "Input file name (on cdc/data) for sigmas.",  string("sigma_v2.dat"));
  //prop-speed
  addParam("PropSpeedFile", m_propSpeedFile, "Input file name (on cdc/data) for prop-speeds.",  string("propspeed_v0.dat"));
  //t0
  addParam("T0File", m_t0File, "Input file name (on cdc/data) for t0s.",  string("t0_v1.dat"));
  //time walk
  addParam("TimeWalkFile", m_twFile, "Input file name (on cdc/data) for time walks.",  string("tw_off.dat"));
  //bad wire
  addParam("BadWireFile", m_bwFile, "Input file name (on cdc/data) for bad wires.",  string("badwire_v1.dat"));
  //channel map
  addParam("ChannelMapFile", m_chMapFile, "Input file name (on cdc/data) for channel map.",  string("ch_map.dat"));
  //edep-to-adc
  addParam("EDepToADCFile", m_eDepToADCFile, "Input file name (on cdc/data) for edep-to-adc.",  string("edeptoadc.dat"));
 
  //max. space resolution
  addParam("MaxSpaceResol", m_maxSpaceResol,
           "Maximum space resolution (cm) in CDCGeometryPar::getSigma() to avoid a too large value; from 2011 beam test; a bit larger value may be better...",
           double(2.5 * 0.0130));
 
  //additional fudge factor for data
  addParam("AddFudgeFactorForSigmaForData", m_addFudgeFactorForSigmaForData,
           "Additional fudge factor for space resol. used in real data event reconstruction (common to all cells).",
           double(1.));
 
  //additional fudge factor for MC
  addParam("AddFudgeFactorForSigmaForMC", m_addFudgeFactorForSigmaForMC,
           "Additional fudge factor for space resol. used in MC event reconstruction (common to all cells). N.B. This factor is not applied in digitization; only applied in reconstruction. Pplease set the fudge parameter of digitizer instead if you want to change the digitization sigma indep. of sigma used in reconstruction.",
           double(1.));
 
  //mapper geometry flag
  addParam("MapperGeometry", m_mapperGeometry, "Define B-field mapper geometry used in GCR in 2017 summer. Tentative option.",
           bool(false));
 
  //mapper phi-angle
  addParam("MapperPhiAngle", m_mapperPhiAngle, "Phi-angle (deg.) of B-field mapper used in GCR in 2017 summer. Tentative option.",
           double(16.7));
 
}
 
void CDCJobCntlParModifierModule::initialize()
{
  //  B2INFO("CDCJobCntlParModifierModule::initialize() called.");
  //For Simulation
  if (m_scp.getTimeWalk() != m_timeWalk) {
    B2INFO("CDCJobCntlParModifier: timeWalk modified: " << m_scp.getTimeWalk() << " to " << m_timeWalk);
    m_scp.setTimeWalk(m_timeWalk);
  }
 
  if (m_scp.getWireSag() != m_wireSag) {
    B2INFO("CDCJobCntlParModifier: wireSag modified: " << m_scp.getWireSag() << " to " << m_wireSag);
    m_scp.setWireSag(m_wireSag);
  }
 
  if (m_scp.getModLeftRightFlag()  !=  m_modLeftRightFlag) {
    B2INFO("CDCJobCntlParModifier: modLeftRightFlag modified: " << m_scp.getModLeftRightFlag()  << " to " << m_modLeftRightFlag);
    m_scp.setModLeftRightFlag(m_modLeftRightFlag);
  }
 
  if (m_scp.getDebug()  !=  m_debug4Sim) {
    B2INFO("CDCJobCntlParModifier: debug4Sim modified: " << m_scp.getDebug()  << " to " << m_debug4Sim);
    m_scp.setDebug(m_debug4Sim);
  }
 
  if (m_scp.getThresholdEnergyDeposit() != m_thresholdEnergyDeposit) {
    B2INFO("CDCJobCntlParModifier: thresholdEnergyDeposit modified: " << m_scp.getThresholdEnergyDeposit() << " to " <<
           m_thresholdEnergyDeposit);
    m_scp.setThresholdEnergyDeposit(m_thresholdEnergyDeposit);
  }
 
  if (m_scp.getMinTrackLength() != m_minTrackLength) {
    B2INFO("CDCJobCntlParModifier: minTrackLength modified: "  << m_scp.getMinTrackLength() << " to " << m_minTrackLength);
    m_scp.setMinTrackLength(m_minTrackLength);
  }
 
  //For Geometry
  if (m_gcp.getMaterialDefinitionMode() != m_materialDefinitionMode) {
    B2INFO("CDCJobCntlParModifier: materialDefinitionMode modified: " << m_gcp.getMaterialDefinitionMode() << " to " <<
           m_materialDefinitionMode);
    m_gcp.setMaterialDefinitionMode(m_materialDefinitionMode);
  }
 
  if (m_gcp.getSenseWireZposMode() != m_senseWireZposMode) {
    B2INFO("CDCJobCntlParModifier: senseWireZposMode modified: " << m_gcp.getSenseWireZposMode() << " " << m_senseWireZposMode);
    m_gcp.setSenseWireZposMode(m_senseWireZposMode);
  }
 
  if (m_gcp.getPrintMaterialTable() != m_printMaterialTable) {
    B2INFO("CDCJobCntlParModifier: printMaterialTable modified: " << m_gcp.getPrintMaterialTable() << " to " << m_printMaterialTable);
    m_gcp.setPrintMaterialTable(m_printMaterialTable);
  }
 
  if (m_gcp.getDebug() != m_debug4Geo) {
    B2INFO("CDCJobCntlParModifier: debug4Geo modified: " << m_gcp.getDebug() << " to " << m_debug4Geo);
    m_gcp.setDebug(m_debug4Geo);
  }
 
  if (m_gcp.getDisplacement() != m_displacement) {
    B2INFO("CDCJobCntlParModifier: displacement switch modified: " << m_gcp.getDisplacement() << " to " << m_displacement);
    m_gcp.setDisplacement(m_displacement);
  }
 
  if (m_gcp.getAlignment() != m_alignment) {
    B2INFO("CDCJobCntlParModifier: alignment switch modified: " << m_gcp.getAlignment() << " to " << m_alignment);
    m_gcp.setAlignment(m_alignment);
  }
 
  if (m_gcp.getMisalignment() != m_misalignment) {
    B2INFO("CDCJobCntlParModifier: misalignment switch modified: " << m_gcp.getMisalignment() << " to " << m_misalignment);
    m_gcp.setMisalignment(m_misalignment);
  }
 
  if (m_gcp.getDisplacementInputType() != m_displacementInputType) {
    B2INFO("CDCJobCntlParModifier: displacementInputType modified: " << m_gcp.getDisplacementInputType() << " to " <<
           m_displacementInputType);
    m_gcp.setDisplacementInputType(m_displacementInputType);
  }
 
  if (m_gcp.getAlignmentInputType() != m_alignmentInputType) {
    B2INFO("CDCJobCntlParModifier: alignmentInputType modified: " << m_gcp.getAlignmentInputType() << " to " << m_alignmentInputType);
    m_gcp.setAlignmentInputType(m_alignmentInputType);
  }
 
  if (m_gcp.getMisalignmentInputType() != m_misalignmentInputType) {
    B2INFO("CDCJobCntlParModifier: misalignmentInputType modified: " << m_gcp.getMisalignmentInputType() << " to " <<
           m_misalignmentInputType);
    m_gcp.setMisalignmentInputType(m_misalignmentInputType);
  }
 
  if (m_gcp.getXtInputType() != m_xtInputType) {
    B2INFO("CDCJobCntlParModifier: xtInputType modified: " << m_gcp.getXtInputType() << " to " << m_xtInputType);
    m_gcp.setXtInputType(m_xtInputType);
  }
 
  if (m_gcp.getSigmaInputType() != m_sigmaInputType) {
    B2INFO("CDCJobCntlParModifier: sigmaInputType modified: " << m_gcp.getSigmaInputType() << " to " << m_sigmaInputType);
    m_gcp.setSigmaInputType(m_sigmaInputType);
  }
 
  if (m_gcp.getPropSpeedInputType() != m_propSpeedInputType) {
    B2INFO("CDCJobCntlParModifier: propSpeedInputType modified: " << m_gcp.getPropSpeedInputType() << " to " << m_propSpeedInputType);
    m_gcp.setPropSpeedInputType(m_propSpeedInputType);
  }
 
  if (m_gcp.getT0InputType() != m_t0InputType) {
    B2INFO("CDCJobCntlParModifier: t0InputType modified: " << m_gcp.getT0InputType() << " to " << m_t0InputType);
    m_gcp.setT0InputType(m_t0InputType);
  }
 
  if (m_gcp.getTwInputType() != m_twInputType) {
    B2INFO("CDCJobCntlParModifier: twInputType modified: " << m_gcp.getTwInputType() << " to " << m_twInputType);
    m_gcp.setTwInputType(m_twInputType);
  }
 
  if (m_gcp.getBwInputType() != m_bwInputType) {
    B2INFO("CDCJobCntlParModifier: bwInputType modified: " << m_gcp.getBwInputType() << " to " << m_bwInputType);
    m_gcp.setBwInputType(m_bwInputType);
  }
 
  if (m_gcp.getChMapInputType() != m_chMapInputType) {
    B2INFO("CDCJobCntlParModifier: chMapInputType modified: " << m_gcp.getChMapInputType() << " to " << m_chMapInputType);
    m_gcp.setChMapInputType(m_chMapInputType);
  }
 
  if (m_gcp.getEDepToADCInputType() != m_eDepToADCInputType) {
    B2INFO("CDCJobCntlParModifier: edep-to-adc modified: " << m_gcp.getEDepToADCInputType() << " to " << m_eDepToADCInputType);
    m_gcp.setEDepToADCInputType(m_eDepToADCInputType);
  }
 
  if (m_gcp.getDisplacementFile() != m_displacementFile) {
    B2INFO("CDCJobCntlParModifier: displacementFile modified: " << m_gcp.getDisplacementFile() << " to " << m_displacementFile);
    m_gcp.setDisplacementFile(m_displacementFile);
  }
 
  if (m_gcp.getAlignmentFile() != m_alignmentFile) {
    B2INFO("CDCJobCntlParModifier: alignmentFile modified: " << m_gcp.getAlignmentFile() << " to " << m_alignmentFile);
    m_gcp.setAlignmentFile(m_alignmentFile);
  }
 
  if (m_gcp.getMisalignmentFile() != m_misalignmentFile) {
    B2INFO("CDCJobCntlParModifier: misalignmentFile modified: " << m_gcp.getMisalignmentFile() << " to " << m_misalignmentFile);
    m_gcp.setMisalignmentFile(m_misalignmentFile);
  }
 
  if (m_gcp.getXtFile() != m_xtFile) {
    B2INFO("CDCJobCntlParModifier: xtFile modified: " << m_gcp.getXtFile() << " to " << m_xtFile);
    m_gcp.setXtFile(m_xtFile);
  }
 
  if (m_gcp.getSigmaFile() != m_sigmaFile) {
    B2INFO("CDCJobCntlParModifier: sigmaFile modified: " << m_gcp.getSigmaFile() << " to " << m_sigmaFile);
    m_gcp.setSigmaFile(m_sigmaFile);
  }
 
  if (m_gcp.getPropSpeedFile() != m_propSpeedFile) {
    B2INFO("CDCJobCntlParModifier: propSpeedFile modified: " << m_gcp.getPropSpeedFile() << " to " << m_propSpeedFile);
    m_gcp.setPropSpeedFile(m_propSpeedFile);
  }
 
  if (m_gcp.getT0File() != m_t0File) {
    B2INFO("CDCJobCntlParModifier: t0File modified: " << m_gcp.getT0File() << " to " << m_t0File);
    m_gcp.setT0File(m_t0File);
  }
 
  if (m_gcp.getTwFile() != m_twFile) {
    B2INFO("CDCJobCntlParModifier: timeWalkFile modified: " << m_gcp.getTwFile() << " to " << m_twFile);
    m_gcp.setTwFile(m_twFile);
  }
 
  if (m_gcp.getBwFile() != m_bwFile) {
    B2INFO("CDCJobCntlParModifier: badWireFile modified: " << m_gcp.getBwFile() << " to " << m_bwFile);
    m_gcp.setBwFile(m_bwFile);
  }
 
  if (m_gcp.getChMapFile() != m_chMapFile) {
    B2INFO("CDCJobCntlParModifier: chMapFile modified: " << m_gcp.getChMapFile() << " to " << m_chMapFile);
    m_gcp.setChMapFile(m_chMapFile);
  }
 
  if (m_gcp.getEDepToADCFile() != m_eDepToADCFile) {
    B2INFO("CDCJobCntlParModifier: edep-to-adcFile modified: " << m_gcp.getEDepToADCFile() << " to " << m_eDepToADCFile);
    m_gcp.setEDepToADCFile(m_eDepToADCFile);
  }
 
  if (m_gcp.getMaxSpaceResolution() != m_maxSpaceResol) {
    B2INFO("CDCJobCntlParModifier: maxSpaceResol modified: " << m_gcp.getMaxSpaceResolution() << " to " << m_maxSpaceResol);
    m_gcp.setMaxSpaceResolution(m_maxSpaceResol);
  }
 
  if (m_gcp.getAddFudgeFactorForSigmaForData() != m_addFudgeFactorForSigmaForData) {
    B2INFO("CDCJobCntlParModifier: addFudgeFactorForSigmaForData modified: " << m_gcp.getAddFudgeFactorForSigmaForData() << " to "
           << m_addFudgeFactorForSigmaForData);
    m_gcp.setAddFudgeFactorForSigmaForData(m_addFudgeFactorForSigmaForData);
  }
 
  if (m_gcp.getAddFudgeFactorForSigmaForMC() != m_addFudgeFactorForSigmaForMC) {
    B2INFO("CDCJobCntlParModifier: addFudgeFactorForSigmaForMC modified: " << m_gcp.getAddFudgeFactorForSigmaForMC() << " to " <<
           m_addFudgeFactorForSigmaForMC);
    m_gcp.setAddFudgeFactorForSigmaForMC(m_addFudgeFactorForSigmaForMC);
  }
 
  if (m_gcp.getMapperGeometry() != m_mapperGeometry) {
    B2INFO("CDCJobCntlParModifier: mapper geometry flag modified: " << m_gcp.getMapperGeometry() << " to " << m_mapperGeometry);
    m_gcp.setMapperGeometry(m_mapperGeometry);
  }
 
  if (m_gcp.getMapperPhiAngle() != m_mapperPhiAngle) {
    B2INFO("CDCJobCntlParModifier: mapper phi-angle modified: " << m_gcp.getMapperPhiAngle() << " to " << m_mapperPhiAngle);
    m_gcp.setMapperPhiAngle(m_mapperPhiAngle);
  }
}
 
void CDCJobCntlParModifierModule::event()
{
}