コード例 #1
0
ファイル: alexsim.cpp プロジェクト: jkriege2/QuickFit3
void AlexSim::simScenario1()
{
    simInit();
    double muTimestep=2*burstDuration;

    background(muTimestep);
    burst();
    background(muTimestep);
    burst();
    background(muTimestep);
    burstBleachingDonor();
    background(muTimestep);
    burstBlinkingAcceptor();
    background(muTimestep);
    burst();
    background(muTimestep);
    burst();
    background(muTimestep);
    burstCoincidenceDonor();
    background(muTimestep);
    burst();
    background(muTimestep);

    writeHist("AlexSimHist.txt");
    writeArrivalTimes("./AlexSimPhotonTimes.txt");
}
コード例 #2
0
ファイル: alexsim.cpp プロジェクト: jkriege2/QuickFit3
void AlexSim::simScenario2(bool writeToFile)
{
    simInit();
    double muTimestep=5*burstDuration;

    background(muTimestep);    burst();
    background(muTimestep);    burst();
    background(muTimestep);    burstFRET(0.5);
    background(muTimestep);    burstFRET(0.5);
    background(muTimestep);    burst();
    background(muTimestep);    burstAcceptorOnly();
    background(muTimestep);    burstAcceptorOnly();
    background(muTimestep);    burstDonorOnly();
    background(muTimestep);    burstDonorOnly();

    background(muTimestep);    burstDelayedDonor();
    background(muTimestep);    burstDelayedDonor();
    background(muTimestep);    burstDelayedAcceptor();
    background(muTimestep);    burstDelayedAcceptor();

    background(muTimestep);    burstBleachingDonor();
    background(muTimestep);    burstBleachingDonor();
    background(muTimestep);    burstBleachingAcceptor();
    background(muTimestep);    burstBleachingAcceptor();

    background(muTimestep);    burstBlinkingDonor();
    background(muTimestep);    burstBlinkingDonor();
    background(muTimestep);    burstBlinkingAcceptor();
    background(muTimestep);    burstBlinkingAcceptor();

    background(muTimestep);    burstDelayedDonor();     burstAcceptorOnly();
    background(muTimestep);    burstDelayedDonor();     burstAcceptorOnly();
    background(muTimestep);    burstDelayedAcceptor();  burstDonorOnly();
    background(muTimestep);    burstDelayedAcceptor();  burstDonorOnly();

    background(muTimestep);    burstDelayedDonor();     burstDonorOnly();
    background(muTimestep);    burstDelayedDonor();     burstDonorOnly();
    background(muTimestep);    burstDelayedAcceptor();  burstAcceptorOnly();
    background(muTimestep);    burstDelayedAcceptor();  burstAcceptorOnly();

    background(muTimestep);    burstCoincidenceDonor();
    background(muTimestep);    burstCoincidenceDonor();
    background(muTimestep);    burstCoincidenceAcceptor();
    background(muTimestep);    burstCoincidenceAcceptor();

    background(muTimestep);

    writeHist("AlexSimHist.txt");
    if (writeToFile) writeArrivalTimes("./AlexSimPhotonTimes.csv");
}
コード例 #3
0
ファイル: Zprime_Py8.C プロジェクト: kalanand/UserCode
bool Zprime_Py8::write(const std::string& outName)
{
  TFile* file = new TFile(outName.c_str(),"RECREATE");
  
  writeHist(h_nsigparts);
  writeHist(h_etasigparts);
  writeHist(h_ptsigparts);

  writeHist(h_nallparts);
  writeHist(h_etaallparts);
  writeHist(h_ptallparts);

  writeHist(h_npupparts);
  writeHist(h_etapupparts);
  writeHist(h_ptpupparts);

  for(std::map<TString, TH1*>::iterator 
	iter = m_HistNames1D.begin(); 
      iter != m_HistNames1D.end(); ++iter) {
    writeHist(iter->second);
    m_HistNames1D.erase(iter);
  }

  for(std::map<TString, TProfile*>::iterator 
	iter2 = m_HistNames2D.begin(); 
      iter2 != m_HistNames2D.end(); ++iter2) {
    writeHist(iter2->second);
    m_HistNames2D.erase(iter2);
  }


  file->Close();
  return true;

}