コード例 #1
0
ファイル: sc_io.cpp プロジェクト: DowerChest/codeblocks
 bool RenameFile(const wxString& src, const wxString& dst)
 {
     wxFileName fname1(Manager::Get()->GetMacrosManager()->ReplaceMacros(src));
     wxFileName fname2(Manager::Get()->GetMacrosManager()->ReplaceMacros(dst));
     NormalizePath(fname1, wxEmptyString);
     NormalizePath(fname2, wxEmptyString);
     if (!SecurityAllows(_T("RenameFile"), wxString::Format(_T("%s -> %s"),
                                     fname1.GetFullPath().c_str(), fname2.GetFullPath().c_str())))
         return false;
     if (!wxFileExists(fname1.GetFullPath())) return false;
     return wxRenameFile(fname1.GetFullPath(),
                         fname2.GetFullPath());
 }
コード例 #2
0
ファイル: main.cpp プロジェクト: Jussmith01/Cpp2011Testing
/* --------------------------------
               Main
---------------------------------- */
int main() {
    std::string fname1("screenwriter:fontone.dat");
    std::string fname2("screenwriter:fonttwo.dat");
    std::string fname3("imageloader:imageone.dat");
    std::string fname4("imageloader:imagetwo.dat");
    std::string fname5("imageloader:imagethree.dat");

    MenuResourceManager mrm;
    std::shared_ptr<MenuResourceBase> p1 = mrm.requestResource(fname1);
    std::shared_ptr<MenuResourceBase> p2 = mrm.requestResource(fname1);
    std::shared_ptr<MenuResourceBase> p3 = mrm.requestResource(fname2);
    std::shared_ptr<MenuResourceBase> p4 = mrm.requestResource(fname3);
    std::shared_ptr<MenuResourceBase> p5 = mrm.requestResource(fname4);
    std::shared_ptr<MenuResourceBase> p6 = mrm.requestResource(fname5);

    unsigned int n( std::thread::hardware_concurrency() );
    std::cout << n << " concurrent threads are supported.\n";

    unsigned N(0);
    while (true) { // Simulate a graphix engine loop
        mrm.manageResources();

        std::cout << "Hi, I'm the main thread, waiting for my children threads to finish!\n";
        std::cout << "COUNT: " << N << "\n";
        std::this_thread::sleep_for( std::chrono::milliseconds(100) );

        mrm.printLoad();

        if ( p1 && p1.get()->is_loaded() )
            p1.reset();

        if ( p2 && p2.get()->is_loaded() )
            p2.reset();

        ++N;
    }

    std::cout << "Hi, I'm Main after the threads finished loading!\n";
}
コード例 #3
0
ファイル: MakeSkimQAPlots.C プロジェクト: aashaqshah/cmssw-1
void MakeSkimQAPlots(TString dataset = "")
{
  gStyle->SetMarkerStyle(20);
  gStyle->SetLineWidth(3);
  gStyle->SetPadColor(10);
  gStyle->SetCanvasColor(10);
  gStyle->SetTitleFont(42, "XYZ");
  gStyle->SetTitleSize(0.06, "XYZ");
  gStyle->SetTitleXOffset(0.7);
  gStyle->SetTitleYOffset(1.05);


  TString fname1("mumu.recolevel." + dataset + ".root");
  TString fname2("ee.recolevel." + dataset + ".root");

  //  TFile *f1 = new TFile(fname1);
  //  TFile *f2 = new TFile(fname2);

  //  TTree *t1 = f1->Get("ntp1");
  //  TTree *t2 = f2->Get("ntp1");

  TChain *t1 = new TChain("ntp1");
  t1->Add("crab_0_071121_090253/res/mumu.recolevel_*.root");
  TChain *t2 = new TChain("ntp1");
  t2->Add("crab_0_071121_102409/res/ee.recolevel_*.root");

  TH1F *hpsim = new TH1F("hpsim","hpsim",50,2.5,4.5);
  TH1F *hpsie = new TH1F("hpsie","hpsie",50,2.5,4.5);
  TH1F *hupsm = new TH1F("hupsm","hupsm",50,8,12);
  TH1F *hupse = new TH1F("hupse","hupse",50,8,12);
  TH1F *hzm = new TH1F("hzm","hzm",50,60,120);
  TH1F *hze = new TH1F("hze","hze",50,60,120);

  TCanvas *c1 = new TCanvas("c1","#mu^{+}#mu^{-}",1400,800);
  c1->Divide(3,3);
  c1->cd(1);
  hpsim->SetXTitle("#psi mass");
  t1->Draw("MuMu_mass >> hpsim");
  c1->cd(2);
  hupsm->SetXTitle("#Upsilon mass");
  t1->Draw("MuMu_mass >> hupsm");
  c1->cd(3);
  hzm->SetXTitle("Z mass");
  t1->Draw("MuMu_mass >> hzm");
  c1->cd(4);
  t1->Draw("MuonCand_pt");
  htemp->SetXTitle("p_{T} (#mu)");
  htemp->Draw();
  c1->cd(5);
  t1->Draw("MuonCand_pt[0]-MuonCand_pt[1]");
  htemp->SetXTitle("#Delta p_{T}");
  htemp->Draw();
  c1->cd(6);
  t1->Draw("MuMu_dphi");
  htemp->SetXTitle("#Delta #phi");
  htemp->Draw();
  c1->cd(7);
  t1->Draw("CaloTower_eta");
  htemp->SetXTitle("#eta (towers)");
  htemp->Draw();
  c1->cd(8);
  c1->SetLogy();
  t1->Draw("CaloTower_e");
  htemp->SetXTitle("E (towers)");
  htemp->Draw();
  c1->SetLogy(0);
  c1->cd(9);
  t1->Draw("nExtraCaloTowersE5");
  htemp->SetXTitle("Extra towers (E > 5 GeV)");
  htemp->Draw();
  TString muname = "dimuons.skimqa." + dataset + ".png";
  c1->SaveAs(muname);

  TCanvas *c2 = new TCanvas("c2","e^{+}e^{-}",1400,800);
  c2->Divide(3,3);
  c2->cd(1);
  hpsie->SetXTitle("#psi mass");
  t2->Draw("ElEl_mass >> hpsie");
  c2->cd(2);
  hupse->SetXTitle("#Upsilon mass");
  t2->Draw("ElEl_mass >> hupse");
  c2->cd(3);
  hze->SetXTitle("Z mass");
  t2->Draw("ElEl_mass >> hze");
  c2->cd(4);
  t2->Draw("EleCand_et");
  htemp->SetXTitle("E_{T} (e)");
  htemp->Draw();
  c2->cd(5);
  t2->Draw("EleCand_et[0]-EleCand_et[1]");
  htemp->SetXTitle("#Delta E_{T}");
  htemp->Draw();
  c2->cd(6);
  t2->Draw("ElEl_dphi");
  htemp->SetXTitle("#Delta #phi");
  htemp->Draw();
  c2->cd(7);
  t2->Draw("CaloTower_eta");
  htemp->SetXTitle("#eta (towers)");
  htemp->Draw();
  c2->cd(8);
  c2->SetLogy();
  t2->Draw("CaloTower_e");
  htemp->SetXTitle("E (towers)");
  htemp->Draw();
  c2->SetLogy(0);
  c2->cd(9);
  t2->Draw("nExtraCaloTowersE5");
  htemp->SetXTitle("Extra towers (E > 5 GeV)");
  htemp->Draw();
  TString elname = "dielectrons.skimqa." + dataset + ".png";
  c2->SaveAs(elname);


}
コード例 #4
0
ファイル: main.cpp プロジェクト: marras/fluor
/******************************
 * Cross-correlation analysis *
 ******************************/
void CalcCrossCorrelation (char * filename1 = "test.flu",char * filename2 = "test2.flu") {
    FILE * fp_params = fopen ("params.txt","rt");
    if (fp_params == NULL) {LOG("!Could not find parameters file!"); return;}

    std::string fname1 (filename1);
    std::string fname2 (filename2);
    std::string out_name ("Xcorr");
    
    if (fname1 != "test.flu" || fname2 != "test2.flu") //in case of custom filename, change the output name too
	out_name += fname1.substr(8,100) += fname2.substr(8,100);

    out_name += ".dat";
    LOG ("*Output correlation to %s",out_name.c_str());

    int num_steps = -1;
    double dT = -1.0, TAU_MAX = -1.0;
    int STEPS_PER_BIN = -1;

    if (fscanf (fp_params,"NUM_FRAMES %d\ndT %lf\nTAU_MAX %lf\nBIN_SIZE %d\n",&num_steps, &dT, &TAU_MAX, &STEPS_PER_BIN) != 4)
     { LOG ("!Error reading parameters file (params.txt)."); exit (2); }
    
    LOG ("*Calculating correlation for %s & %s.\nNUM_FRAMES %d, dT %f, TAU_MAX %f, BIN_SIZE %d\n",filename1,filename2,num_steps, dT, TAU_MAX, STEPS_PER_BIN);

    FILE * fp;

    fp = fopen (out_name.c_str(),"wt");
    LOG ("Calculating cross correlation (Logarithmic).");
    fprintf (fp, "@title \"G(tau), dT = %f\"\n",dT);

    int counter = 0; //counter=1;

    if (STEPS_PER_BIN <= 0)	{	//STEPS_PER_BIN undefined = varying bin size
      STEPS_PER_BIN = 1;		//start form bin size = frame size, then increase
      double time = STEPS_PER_BIN*dT;
      Correlator * corr1;

      while (time <= TAU_MAX) {
	  corr1 = new CrossCorrelator (filename1, filename2, dT, num_steps, STEPS_PER_BIN);
	  printf ("STEPS_PER_BIN=%d.\n",STEPS_PER_BIN);

	  double result = corr1->Calculate_correlation (time);
	  fprintf (fp,"%5.9le %5.9le\n",time/1e12, result);
	  printf ("g(%5.9le)=%5.9le\n", time/1e12, result);

	  delete corr1;

	  counter++;
	  if (counter >= 4) {counter=0; STEPS_PER_BIN*=2;}	//every 4 iterations, the bin size (and time change) is doubled
	  time += STEPS_PER_BIN*dT;
      }
    }
    else {				//STEPS_PER_BIN constant, as defined in params.txt
      LOG ("Calculating CROSS-correlation for BIN_SIZE = %d.",STEPS_PER_BIN);

      double time = STEPS_PER_BIN*dT;
      double TIME_STEP = STEPS_PER_BIN;
      Correlator * corr1 = new CrossCorrelator (filename1, filename2, dT, num_steps, STEPS_PER_BIN);

      printf ("STEPS_PER_BIN=%d.\n",STEPS_PER_BIN);

      while (time <= TAU_MAX) {
	  double result = corr1->Calculate_correlation (time);
	  fprintf (fp,"%5.9le %5.9le\n",time/1e12, result);
	  printf ("g(%5.9le)=%5.9le\n", time/1e12, result);

	  counter++;
	  if (counter >= 4) {counter=0; TIME_STEP*=2;}	//every 4 iterations, the correlation time increase is doubled
	  time += TIME_STEP*dT;
      }
      delete corr1;
    }

    fclose (fp);
    LOG ("Finished calculating cross-correlation.");
}