Ejemplo n.º 1
0
int
main ()
{
  //findgraph ( 2, 3, .1, .0001, 5,    &fsqrt, "Find square Root of 5");
  //findgraph (10,500, 5, .0000001,    5487, &fsqrt, "Find square Root of 5487");

  findgraph ( 0, 1,   .1,  .0001,  0, &hw2p1, "hw2p1");
  findgraph ( 1, 3.2, .1,  .0001,  0, &hw2p2, "hw2p2");
  findgraph ( 1, 2,   .1,  .0001,  0, &hw2p3, "hw2p3");

  findgraph (-1, 0,   .01, .001,  0, &hw2p4, "hw2p4");
  secant (-1, 0, .001,  0, &hw2p4, "hw2p4");
  falseposition (-1, 0, .001,  0, &hw2p4, "hw2p4");

  findgraph (-1, 0,   .01, .00001, 0, &hw2p5, "hw2p5 from -1 to 0");
  findgraph ( 0, 1,   .01, .00001, 0, &hw2p5, "hw2p5 from 0 to 1");
  findgraph (-1 , 1,  .1,  .00001, 0, &hw2p5, "hw2p5 from 0 to 1");

  findpoint (0, 2, 0, 1, .001, 8, &hw2p6, "hw2p6");

  // testing the point function
  //findpoint (0, 2, 1, 1, .001, 8, &hw2p6, "hw2p6");
  //findpoint (1, 3, 2, 4, .001, 8, &hw2p6, "hw2p6");
  //findpoint (1, 4, 0, 2, .001, 8, &hw2p6, "hw2p6");
  //findpoint (1, 8, 18, 0, .001, 8, &hw2p6, "hw2p6");

  findgraph (0, 1, .1,  .01,  12.4,  &hw2p7a, "hw2p7 trough with 12.4 volume");

///*
  printf ("finding a range of values for testing.\n");
  long double i;
  for (i=0; i< 1; i=i+.1){
    char name[40];
    sprintf(name, "hw2p7 %Lf", i);
    //findgraph (0, 1, .01, .00001, i,    &hw2p7, name);
    hw2p7(i, 0);
    hw2p7a(i, 0);
  }
//*/

  //findgraph (-1, 3, .1, .001, 38, &hw2ptri, "hw2ptri full");
  //findgraph (-1, 3, .1, .001, 10, &hw2ptri, "hw2ptri half full");
  //findgraph (-1, 3, .1, .001, .1, &hw2ptri, "hw2ptri empty");

  //findgraph   ( 0, 2, .1, .00000001, 0, &fex, "fex");
}
Ejemplo n.º 2
0
Archivo: findcent.C Proyecto: XuQiao/HI
void findcent()
{
	const int color[10] = {2,3,4,5,6,7,8,4,3,2};
	const double centbin[] = {0,0.005,0.01,0.05,0.1,0.2,0.3,0.4,0.6,1.0};
	//const double centbin[] = {0,0.001,0.005,0.01,0.05,0.1,0.2,0.4,0.6,1.0};
        //double centbin[101];
        //for(int i=0;i<101;i++){
        //    centbin[i]=i*0.01;
        //}
	const double erf = 1e-5;
	int ncent = (int)sizeof(centbin)/sizeof(double);
	//TFile *f = TFile::Open("merged_AnapAlfvtxor.root");
	//TFile *f = TFile::Open("merged_AnapAlfvtxand.root");
	//TFile *f = TFile::Open("merged_AnapAlfvtxsouth.root");
	TFile *f = TFile::Open("merged_AnaMWGppmb.root");
//	TH2F* pc1hitsbbcdis = (TH2F*)f->Get("pc1hitsbbcdis"); 
//	TH1F* pc1hitsdis = (TH1F*)f->Get("pc1hitsdis"); 
//	TH2F* hbbcntrk = (TH2F*)f->Get("hbbcnbbc");
	TH2F* hbbcntrk = (TH2F*)f->Get("hnfvtxtrkbbc");
//	hbbcntrk->RebinX(3);
//	TH1F* bbcsouthdis = (TH1F*)pc1hitsbbcdis->ProjectionX(0,-1);
	TH1F* bbcsouthdis = (TH1F*)hbbcntrk->ProjectionX(0,-1);
//	TH2F* bbcsouthdis = (TH2F*)hbbcntrk->Clone("bbcsouthdis");
	bbcsouthdis->SetTitle("BBC charge distribution");	
	bbcsouthdis->GetXaxis()->SetTitle("BBC charge");	
	bbcsouthdis->GetYaxis()->SetTitle("# of Events");	
//	pc1hitsdis->SetTitle("PC1 hits distribution");	
//	pc1hitsdis->GetXaxis()->SetTitle("PC1 hits");	
//	pc1hitsdis->GetYaxis()->SetTitle("# of Events");	


	ofstream fstr("centbin_bbcmb.txt");
//	fstr<<"centbin"<<"\t"<<"bbcsouth"<<"\t"<<"pc1hits"<<"\t"<<"bbc ratio"<<"\t"<<"pc1hits ratio"<<endl;
	fstr<<"centbin"<<"\t"<<"bbc"<<"\t"<<"bbc ratio"<<endl;
	
	vector<double> bbcs;
//	vector<double> pc1hs;

	for(int icent = 0; icent < ncent; icent++){
		double bbcskp = findpoint(bbcsouthdis,centbin[icent]);
//		double pc1hskp = findpoint(pc1hitsdis,centbin[icent]);
		if(bbcskp<0) bbcskp = 0;
//		if(pc1hskp<0) pc1hskp = 0;
		bbcs.push_back(bbcskp);
cout<<bbcskp<<",";
//		pc1hs.push_back(pc1hskp);
	}
	
	for(int icent = 0; icent < ncent; icent++){
	//	fstr<<centbin[icent]<<"\t"<<bbcs[icent]<<"\t"<<pc1hs[icent]<<"\t"<<bbcsouthdis->Integral(bbcsouthdis->FindBin(bbcs[icent]+erf),bbcsouthdis->FindBin(bbcs[0]-erf))/bbcsouthdis->Integral(bbcsouthdis->FindBin(bbcs[ncent-1]-erf),bbcsouthdis->FindBin(bbcs[0]+erf))<<"\t"<<pc1hitsdis->Integral(pc1hitsdis->FindBin(pc1hs[icent]+erf),pc1hitsdis->FindBin(pc1hs[0]-erf))/pc1hitsdis->Integral(pc1hitsdis->FindBin(pc1hs[ncent-1]-erf),pc1hitsdis->FindBin(pc1hs[0]+erf))<<endl;
		fstr<<centbin[icent]<<"\t"<<bbcs[icent]<<"\t"<<bbcsouthdis->Integral(bbcsouthdis->FindBin(bbcs[icent]+erf),bbcsouthdis->FindBin(bbcs[0]-erf))/bbcsouthdis->Integral(bbcsouthdis->FindBin(bbcs[ncent-1]-erf),bbcsouthdis->FindBin(bbcs[0]+erf))<<endl;
		//fstr<<centbin[icent]<<"\t"<<bbcs[icent]<<"\t"<<bbcsouthdis->Integral(bbcsouthdis->GetXaxis()->FindBin(bbcs[icent]+erf),bbcsouthdis->GetXaxis()->FindBin(bbcs[0]-erf),bbcsouthdis->GetYaxis()->FindBin(bbcs[icent]+erf),bbcsouthdis->GetYaxis()->FindBin(bbcs[0]-erf))/bbcsouthdis->Integral(bbcsouthdis->FindBin(bbcs[ncent-1]-erf,bbcs[ncent-1]-erf),bbcsouthdis->FindBin(bbcs[0]+erf,bbcs[0]+erf))<<endl;
	}

	TCanvas *c1 = new TCanvas();
	//c1->SetLogx();
	c1->SetLogy();
//	c1->SetLogz();
	bbcsouthdis->GetXaxis()->SetRangeUser(0,500);
//	bbcsouthdis->GetYaxis()->SetRangeUser(0,100);
	bbcsouthdis->Draw();
	for(int icent = 0; icent < ncent; icent++){
		TH1F* bbcs_t = (TH1F*)bbcsouthdis->Clone(Form("bbcs_%d",icent));
		bbcs_t->GetXaxis()->SetRangeUser(bbcs[icent+1]-erf,bbcs[icent]+erf);
		//bbcs_t->GetYaxis()->SetRangeUser(bbcs[icent+1]-erf,bbcs[icent]+erf);
		bbcs_t->SetFillColor(color[icent]);
		bbcs_t->Draw("HIST same");
	}
	c1->Print("bbcaddcent.png");
	/*
 	TCanvas *c2 = new TCanvas();
	c2->SetLogy();
	pc1hitsdis->Draw();
	for(int icent = 0; icent < ncent; icent++){
		TH1F* pc1hs_t = (TH1F*)pc1hitsdis->Clone(Form("pc1hs_%d",icent));
		pc1hs_t->GetXaxis()->SetRangeUser(pc1hs[icent+1]-erf,pc1hs[icent]+erf);
		pc1hs_t->SetFillColor(color[icent]);
		pc1hs_t->Draw("HIST same");
	}
	c2->Print("pc1hits.png");
*/
}