test_MinimalLArTPC.fcl

Go to the documentation of this file.

#no experiment specific configurations because SingleGen is detector agnostic

#include "simulationservices.fcl"
#include "spacecharge.fcl"
#include "larproperties.fcl"
######include "detectorproperties.fcl"
#include "detectorclocks.fcl"
#include "larfft.fcl"

standard_singlep:
{
 module_type:           "SingleGen" 
 ParticleSelectionMode: "all"       # 0 = use full list, 1 =  randomly select a single listed particle
 PadOutVectors:        false        # false: require all vectors to be same length
                                    # true:  pad out if a vector is size one    
 PDG:                   [ 13 ]      # list of pdg codes for particles to make
 P0:                    [ 6. ]      # central value of momentum for each particle
 SigmaP:                [ 0. ]      # variation about the central value 
 PDist:                 "Gaussian"  # 0 - uniform, 1 - gaussian distribution                               
 X0:                    [ 0. ]     # in cm in world coordinates, ie x = 0 is at the wire plane          
                                    # and increases away from the wire plane                             
 Y0:                    [ 0. ]      # in cm in world coordinates, ie y = 0 is at the center of the TPC   
 Z0:                    [ -130. ]   # in cm in world coordinates, ie z = 0 is at the upstream edge of    
                                    # the TPC and increases with the beam direction                      
 T0:                    [ 0. ]      # starting time
 SigmaX:                [ 0. ]      # variation in the starting x position
 SigmaY:                [ 0. ]      # variation in the starting y position                            
 SigmaZ:                [ 0.0 ]     # variation in the starting z position                            
 SigmaT:                [ 0.0 ]     # variation in the starting time
 PosDist:               "uniform"   # 0 - uniform, 1 - gaussian
 TDist:                 "uniform"   # 0 - uniform, 1 - gaussian         
 Theta0XZ:              [ 0. ]      #angle in XZ plane (degrees)                                         
 Theta0YZ:              [ 0. ]      #angle in YZ plane (degrees)                                         
 SigmaThetaXZ:          [ 0. ]      #in degrees                                                          
 SigmaThetaYZ:          [ 0. ]      #in degrees                                                          
 AngleDist:             "Gaussian"  # 0 - uniform, 1 - gaussian
}

process_name:processA

services:
{
  message : {
     debugModules : ["*"]
     suppressInfo : []    
     destinations : {
       LogToConsole : {
         type : "cout" 
         threshold : "DEBUG" 
         categories : {
           default : { limit : 50 }
         }
       }
     }
  }

  TFileService : { fileName: "testlarg4.root" }
  RandomNumberGenerator: {}
  NuRandomService:{
    service_type: "NuRandomService" 
    endOfJobSummary: true
    policy: "random" 
  }

  DetectorHolder: {}
  ActionHolder: {}
  PhysicsListHolder: {} 
  PhysicsList: { 
    PhysicsListName: "FTFP_BERT" 
    DumpList: false    
    enableCerenkov: false
    enableScintillation: false
    ScintillationByParticleType: false
    enableAbsorption: false   
    enableRayleigh: false     
    enableMieHG: false       
    enableBoundary: false    
    enableWLS: false
  }
  LArG4Detector : 
  {
    category: "world" 
    gdmlFileName_ : "MinimalLArTPCdetector_nowires.gdml" 
  }   

  MCTruthEventAction: {service_type: "MCTruthEventActionService"}
  ParticleListAction: {service_type: "ParticleListActionService"}

  Geometry: {
     Name: "MinimalLArTPCdetector_g4"
     GDML: "MinimalLArTPCdetector_g4.gdml" 
     ROOT: "MinimalLArTPCdetector_g4.gdml"
     SurfaceY: 0
     DisableWiresInG4: true
     SortingParameters: {}
     service_type: "Geometry"
  }

  ExptGeoHelperInterface: {
     service_provider : StandardGeometryHelper
     service_type: "ExptGeoHelperInterface"
  }
  
  GeometryConfigurationWriter: {}

  # -- in larsim
  LArG4Parameters:              @local::standard_largeantparameters
  # -- in larevt
  SpaceCharge:                  @local::standard_spacecharge
  # -- following 3 are in lardataalg
  LArPropertiesService:            @local::standard_properties
  DetectorPropertiesService:       {
  service_provider: "DetectorPropertiesServiceStandard"

  # Drift properties 

     SternheimerA:     0.1956  # Ar Sternheimer parameter a.
     SternheimerK:     3.0000  # Ar Sternheimer parameter k.
     SternheimerX0:    0.2000  # Ar Sternheimer parameter x0.
     SternheimerX1:    3.0000  # Ar Sternheimer parameter x0.
     SternheimerCbar:  5.2146  # Ar Sternheimer parameter Cbar.
    
     DriftVelFudgeFactor: 1.
     
     ElectronsToADC:    1.208041e-3  # in ADC/e; 6241.5 electrons = 1fC = 7.54 ADC counts for ArgoNeuT
     NumberTimeSamples: 2048         # ticks of the clock per readout frame
     ReadOutWindowSize: 2048
   
     SimpleBoundaryProcess: true  #enable opticalBoundaryProcessSimple instead of G4 default

Temperature:       87            # in K
Electronlifetime:  3.0e3         # in us
Efield:            [0.5,4.0,0.0] # in kV/cm
ElectronsToADC:    1.02534e-03   # 
NumberTimeSamples: 1667
ReadOutWindowSize: 1667
#TimeOffsetX:       0.
TimeOffsetY:       0.
TimeOffsetZ:       0.

  }
  DetectorClocksService:        @local::standard_detectorclocks

  LArFFT:                  @local::standard_larfft
}

source: {
  module_type: EmptyEvent
  maxEvents:  10
}

physics:
{

  producers:
  {
    generator: @local::standard_singlep
    larg4Main: {
      module_type: larg4Main
      enableVisualization: false
      macroPath: ".:./macros" 
      visMacro: "vis.mac" 
    }
    elecDrift: {
      module_type:     "SimDriftElectrons"
      SimulationLabel: "larg4Main:LArG4DetectorServicevolTPCActive"
      StoreDriftedElectronClusters: true
    }
    simWire: {
      module_type:   "SimWire"
      DriftEModuleLabel:  "elecDrift"
      #NoiseFact:    0.0132 # Noise Scale
      NoiseFact:     132. # Noise Scale
      #NoiseWidth:   62.4   # Exponential Noise width (kHz)
      NoiseWidth:    62.4   # Exponential Noise width (kHz)
      LowCutoff:     7.5    # Low frequency filter cutoff (kHz)
      FieldBins:     75
      Col3DCorrection:    2.5
      Ind3DCorrection:    1.5
      ColFieldRespAmp:    0.0354
      IndFieldRespAmp:    0.018
      ShapeTimeConst:        [ 3000., 900. ]
      CompressionType:    "none"
    }
  }
  analyzers:
  {
    CheckSimEnergyDeposit: {
      module_type: CheckSimEnergyDeposit
      hist_dir: "HistoDir" 
    }
    SimDriftedElectronClusterAna: {
      module_type: SimDriftedElectronClusterAna
      EDepModuleLabel: "elecDrift"
    }               
  }

  path1: [ generator,larg4Main,elecDrift,simWire ]
  stream1:  [ out1,CheckSimEnergyDeposit,SimDriftedElectronClusterAna ]
  trigger_paths: [ path1 ]
  end_paths: [ stream1 ]
}

outputs:
{
  out1:
  {
    module_type: RootOutput   
    fileName: "Testingout.root" 
  }
}