merge_configure.h

Go to the documentation of this file.

void merge_configure(JobC& jc, JobCPath& path, int seed, int runn)
{
  // expect all file to be processed by this
  ReleaseType::Release_t official_rel = ReleaseType::kDogwood |
                                        ReleaseType::kGarlic |
                                        ReleaseType::k00;

  TString alt_rel = gSystem->Getenv("ENV_GEVPERMIP_SET");
  if(alt_rel != ""){
    // user externally chose an alternative, so be it..
    choose_gevpermip_release(alt_rel.Data());
  }else{
    // go with the official choice
    choose_gevpermip_release(official_rel);
  }


  path.SetAllFilters(0);

  // explicitly set the value of GeVPerMip
  JobCModule& photon = path.Mod("PhotonTransport");

  // If we ever want to use this for FD then we'll need to check everywhere
  // that assumes we're using ND.
  double gpm = get_gevpermip(Detector::kNear);
  photon.Set(TString::Format("GeVPerMip=%lf", gpm));


  // Defaults in photon transport as of 24/1/08
  // Total noise rate per second in the whole detector
  photon.Set("FDNoiseRate=3762138");
  // Set non-zero to use the old cedar FD noise model
  photon.Set("UseSimpleNoiseModel=0");

  /*
  if(det == Detector::kFar){
    JobCModule& detsim = path.Mod("DetSim");
    photon.Set("FibreVelocityFudge=1.02");
    detsim.Set("vaTimingWidth=0.12");
  }
  */


  JobCModule& mcmerge = path.Mod("MCMerge");
  // Set the RNG seed (0 selects current time in sec)
  mcmerge.Set(TString::Format("RandomSeed=%d", seed));
  mcmerge.Set(TString::Format("Runn=%d", runn));


  JobCEnv& env=JobCEnv::Instance();


  // Configure detector path
  const TString det_file = env.GetFileName(0);
  cout << "Adding file " << det_file << " to \"detector\" path" << endl;
  jc.Input.AddFile(det_file,"det");
  jc.Input.DefineStream("det","SimSnarl");

  jc.Input.SetSequenceMode("det",Per::kSequential);
  jc.Input.SetPushRandom("det",true);
  jc.Input.SetMeanMom("det",18.64); // value used in carrot processing
  jc.Input.SetMaxFileRepeat("det", 0);


  // Configure rock path
  const TString rock_file = env.GetFileName(1);
  cout << "Adding file " << rock_file << " to \"rock\" path" << endl;
  jc.Input.AddFile(rock_file,"rock");
  jc.Input.DefineStream("rock","SimSnarl");

  jc.Input.SetSequenceMode("rock",Per::kRandom);
  jc.Input.SetPushRandom("rock",true);
  jc.Input.SetMeanMom("rock",21.15); // value used in carrot processing
  jc.Input.SetMaxFileRepeat("rock", 10); // reuse each input rock file 10x


  // If a third file is specified on the command line, it's a special
  // signal file: put it in the signal path
  if(env.GetNfile() > 2){
    const TString sig_file = env.GetFileName(2);
    cout << "Adding file " << sig_file << " to special \"signal\" path" << endl;
    jc.Input.AddFile(sig_file,"signal");
    jc.Input.DefineStream("signal","SimSnarl");

    // TODO: What are these supposed to actually be? I just copied them
    // from the "det" stream...
    jc.Input.SetSequenceMode("signal",Per::kSequential);
    jc.Input.SetPushRandom("signal",true);
    jc.Input.SetMeanMom("signal",18.64); // value used in carrot processing
    jc.Input.SetMaxFileRepeat("signal", 0);
  }
}

---