pileup

• DC1ではatlsim(Geant3ベース)の環境で行った。(Athena 4.0.1)
• Low = 2.0 x 10^33、High = 10^34
• SignalイベントにPileupイベントを埋め込むことで実現している。
• ファイルの混ぜ方 (ATLAS Japan)
  • /home/atljphys/athena/dc1/simul/dset/002099/data/以下にPileupイベントが500ファイルある。
  • Lowでは25ファイル/400イベント、Highでは60ファイル/200イベント必要。
  • 10回ごとに混ぜるファイルを変更していく。
    500ファイルを如何に有効に使うかということがポイント。
    以下がそのときのサンプルコード。
    ex1) calPartition4minibias 1 l 500 ... "1イベントから400イベントまで"用
    ex2) calPartition4minibias 2 l 500 ... "401イベントから800イベントまで"用
    ただし、例えば200イベントしかなくても同じ計算方法を利用することもある。

    #include <iostream>
    #include <string>
    
    using namespace std;
    
    int
    main(int argc, char * argv[])
    {
      int nPartition = atoi(argv[1]);
      int nSet  = 10;
      int nFile = 60; // High luminosity
      if (string(argv[2]) == string("l") ||
          string(argv[2]) == string("L")) {
        nFile = 25; // Low luminosity
      }
      int nTotal = atoi(argv[3]);
    
      int First = ((nPartition-1)/nSet*nFile)%nTotal+1;
      int Last  = First+nFile-1;
    
      cout << First << " " << Last << endl;
    
      return 0;
    }