ROOT logo
//=======================================================================//
//Macro to draw the main results of the charge fluctuation analysis:
//=======================================================================//

//+++++++++++++++++++++GLOBAL VARIABLES+++++++++++++++++++++//
const Int_t nCentralityBins = 20;

Double_t gNParticipants20[nCentralityBins] = {382.646,329.213,280.911,238.617,201.81,169.339,141.067,116.208,94.2515,75.4558,59.4054,45.7565,34.4839,25.3127,18.0293,12.6535,8.84139,6.16348,4.37454,3.06182};
Double_t gNParticipants20Error[nCentralityBins] = {16.9428,18.0505,17.0971,16.2076,15.5418,14.9458,14.3174,13.9067,13.2661,12.6134,11.8133,11.0495,10.0939,8.99737,7.7884,6.48725,5.21602,3.91988,2.78741,1.75066};
Double_t gNParticipants10[nCentralityBins] = {356.132,260.175,185.76,128.428,84.5666,52.3432,29.7072,15.2207,7.47192,3.71973};//10 bins
Double_t gNParticipants10Error[nCentralityBins] = {31.8228,26.9093,22.3767,18.8802,15.9735,13.3103,10.5806,7.62745,4.78538,2.41942};//10 bins

//================================ALICE================================//
Double_t gNuDynALICEData20[nCentralityBins] = {10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.};
Double_t gNuDynALICEData20Error[nCentralityBins] = {0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.};
Double_t gNetChargeALICEData20[nCentralityBins] = {10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.};
Double_t gNetChargeALICEData20Error[nCentralityBins] = {0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.};
Double_t gNetChargeRmsALICEData20[nCentralityBins] = {10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.};
Double_t gNetChargeRmsALICEData20Error[nCentralityBins] = {0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.};

Double_t gNuDynALICEData10[nCentralityBins] = {10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.};
Double_t gNuDynALICEData10Error[nCentralityBins] = {0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.};
Double_t gNetChargeALICEData10[nCentralityBins] = {10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.};
Double_t gNetChargeALICEData10Error[nCentralityBins] = {0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.};
Double_t gNetChargeRmsALICEData10[nCentralityBins] = {10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.};
Double_t gNetChargeRmsALICEData10Error[nCentralityBins] = {0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.};
//================================ALICE================================//
//+++++++++++++++++++++END OF VARIABLES+++++++++++++++++++++//

//_____________________________________________________//
void compareCentralitySets(const char*resultPath1 = "../LHC10h/Pass1_4Plus/20CentralityBins/", const char*resultPath2 = "../LHC10h/Pass1_4Plus/10CentralityBins/") {
  //Draws the nu_dyn vs centrality percentile
  SetDataPoints(resultPath1,resultPath2);
  //================================================//
  //ALICE PbPb @ 2.76 TeV
  TGraphErrors *grALICEData20Nudyn = new TGraphErrors(nCentralityBins,
						    gNParticipants20,
						    gNuDynALICEData20,
						    gNParticipants20Error,
						    gNuDynALICEData20Error);
  grALICEData20Nudyn->SetMarkerStyle(20);
  grALICEData20Nudyn->SetMarkerColor(2);

  TGraphErrors *grALICEData20NetCharge = new TGraphErrors(nCentralityBins,
							gNParticipants20,
							gNetChargeALICEData20,
							gNParticipants20Error,
							gNetChargeALICEData20Error);
  grALICEData20NetCharge->SetMarkerStyle(20);
  grALICEData20NetCharge->SetMarkerColor(2);

  TGraphErrors *grALICEData20NetChargeRms = new TGraphErrors(nCentralityBins,
							   gNParticipants20,
							   gNetChargeRmsALICEData20,
							   gNParticipants20Error,
							   gNetChargeRmsALICEData20Error);
  grALICEData20NetChargeRms->SetMarkerStyle(20);
  grALICEData20NetChargeRms->SetMarkerColor(2);

  //10 bins
  TGraphErrors *grALICEData10NetCharge = new TGraphErrors(nCentralityBins,
							gNParticipants10,
							gNetChargeALICEData10,
							gNParticipants10Error,
							gNetChargeALICEData10Error);
  grALICEData10NetCharge->SetMarkerStyle(24);
  grALICEData10NetCharge->SetMarkerColor(2);

  TGraphErrors *grALICEData10NetChargeRms = new TGraphErrors(nCentralityBins,
							   gNParticipants10,
							   gNetChargeRmsALICEData10,
							   gNParticipants10Error,
							   gNetChargeRmsALICEData10Error);
  grALICEData10NetChargeRms->SetMarkerStyle(24);
  grALICEData10NetChargeRms->SetMarkerColor(2);

  TGraphErrors *grALICEData10Nudyn = new TGraphErrors(nCentralityBins,
						    gNParticipants10,
						    gNuDynALICEData10,
						    gNParticipants10Error,
						    gNuDynALICEData10Error);
  grALICEData10Nudyn->SetMarkerStyle(24);
  grALICEData10Nudyn->SetMarkerColor(2);
  
  //_____________________________________________________//
  //Draw the results
  //_____________________________________________________//
  TLatex *latex = new TLatex();
  latex->SetTextSize(0.035);

  //====================================//
  //Results vs centrality
  TH2F *gEmpty1 = new TH2F("gEmpty1",
			   ";N_{part.};",
			   100,0,400,10000,-0.1,50.);
  gEmpty1->SetStats(kFALSE);
  gEmpty1->GetYaxis()->SetTitleOffset(1.5);
  gEmpty1->GetXaxis()->SetTitleOffset(1.5);
  gEmpty1->GetYaxis()->SetNdivisions(10);
  gEmpty1->GetXaxis()->SetNdivisions(10);

  TF1 *f1 = new TF1("f1","0",0,1000);
  f1->SetLineColor(1); f1->SetLineStyle(3); f1->SetLineWidth(2);

  //============================================================//
  //nu_{dyn.}
  TCanvas *c1 = new TCanvas("c1","Centrality dependence: nu_dyn",
			    0,0,500,500);
  c1->SetFillColor(10); c1->SetHighLightColor(10);
  c1->SetLeftMargin(0.15); c1->SetBottomMargin(0.15);
  c1->SetGridx(); c1->SetGridy();
  gEmpty1->GetYaxis()->SetRangeUser(-0.08,0.01);
  gEmpty1->GetYaxis()->SetTitle("#nu_{dyn.}");
  gEmpty1->DrawCopy();
  f1->Draw("same");
  grALICEData20Nudyn->Draw("P");
  grALICEData10Nudyn->Draw("P");

  DrawMarker(140., -0.043, 20, 1.6, 2);
  latex->DrawLatex(150.,-0.0445,"20 centrality bins");
  DrawMarker(140., -0.05, 24, 1.6, 2);
  latex->DrawLatex(150.,-0.0515,"10 centrality bins");

  c1->SaveAs("comparissonCentralitySetsNparticipantsDependenceNuDyn.png");

  //============================================================//
  //net charge
  TCanvas *c2 = new TCanvas("c2","Centrality dependence: net charge",
			    100,100,500,500);
  c2->SetFillColor(10); c2->SetHighLightColor(10);
  c2->SetLeftMargin(0.15); c2->SetBottomMargin(0.15);
  c2->SetGridx(); c2->SetGridy();
  gEmpty1->GetYaxis()->SetRangeUser(0.0,10.0);
  gEmpty1->GetYaxis()->SetTitle("#LT #Delta Q = N_{+} - N_{-}#GT");
  gEmpty1->DrawCopy();
  grALICEData20NetCharge->Draw("P");
  grALICEData10NetCharge->Draw("P");

  DrawMarker(40., 8.3, 20, 1.6, 2);
  latex->DrawLatex(50.,8.1,"20 centrality bins");
  DrawMarker(40., 7.3, 24, 1.6, 2);
  latex->DrawLatex(50.,7.1,"10 centrality bins");

  c2->SaveAs("comparissonCentralitySetsNparticipantsDependenceNetCharge.png");

  //============================================================//
  //rms net charge
  TCanvas *c3 = new TCanvas("c3","Centrality dependence: net charge rms",
			    200,200,500,500);
  c3->SetFillColor(10); c3->SetHighLightColor(10);
  c3->SetLeftMargin(0.15); c3->SetBottomMargin(0.15);
  c3->SetGridx(); c3->SetGridy();
  gEmpty1->GetYaxis()->SetRangeUser(0.0,40.0);
  gEmpty1->GetYaxis()->SetTitle("#sigma_{#Delta Q}");
  gEmpty1->DrawCopy();
  grALICEData20NetChargeRms->Draw("P");
  grALICEData10NetChargeRms->Draw("P");

  DrawMarker(40., 38., 20, 1.6, 2);
  latex->DrawLatex(50.,37.5,"20 centrality bins");
  DrawMarker(40., 34., 24, 1.6, 2);
  latex->DrawLatex(50.,33.5,"10 centrality bins");

  c3->SaveAs("comparissonCentralitySetsNparticipantsDependenceNetChargeRms.png");
}

//_______________________________________________________________//
void SetDataPoints(const char* resultPath1,
		   const char* resultPath2) {
  //Calculate the mean and the statistical error of the data points
  const Int_t nRuns = 7;
  Double_t nuDyn[nRuns][nCentralityBins];
  Double_t gNetCharge[nRuns][nCentralityBins];
  Double_t gNetChargeRms[nRuns][nCentralityBins];
  Double_t nuDyn[nRuns][nCentralityBins];
  Int_t nEvents[nRuns][nCentralityBins];
  Int_t nRunNumbers[nRuns] = {137161,137431,137549,
			      137595,137638,137639,137693};
  Int_t gCentrality;
  Double_t netChargeValues[nCentralityBins] = {0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,
					   0.,0.,0.,0.,0.,0.,0.,0.,0.,0.};
  Double_t netChargeValuesError[nCentralityBins] = {0.,0.,0.,0.,0.,0.,0.,0.,
						    0.,0.,0.,0.,0.,0.,0.,0.,
						    0.,0.,0.,0.};
  Double_t netChargeRmsValues[nCentralityBins] = {0.,0.,0.,0.,0.,0.,0.,
						  0.,0.,0.,0.,0.,0.,0.,
						  0.,0.,0.,0.,0.,0.};
  Double_t netChargeRmsValuesError[nCentralityBins] = {0.,0.,0.,0.,0.,0.,0.,0.,
						       0.,0.,0.,0.,0.,0.,0.,0.,
						    0.,0.,0.,0.};
  Double_t nuDynValues[nCentralityBins] = {0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,
					   0.,0.,0.,0.,0.,0.,0.,0.,0.,0.};
  Double_t nuDynValuesError[nCentralityBins] = {0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,
						0.,0.,0.,0.,0.,0.,0.,0.,0.,0.};
  Double_t nEventsValues[nCentralityBins] = {0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,
					     0.,0.,0.,0.,0.,0.,0.,0.,0.,0.};

  //20 centrality bins
  TString inputFileName;
  for(Int_t iRun = 0; iRun < nRuns; iRun++) {
    ifstream inputAscii20;
    Printf("Adding run %d",nRunNumbers[iRun]);
    inputFileName = resultPath1;
    inputFileName += "/output."; inputFileName += nRunNumbers[iRun];
    inputFileName += ".txt";

    Printf("Filename: %s",inputFileName.Data());
    inputAscii20.open(inputFileName.Data());
    for(Int_t iCentrality = 0; iCentrality < nCentralityBins; iCentrality++) {
      inputAscii20>>gCentrality>>nuDyn[iRun][iCentrality]>>gNetCharge[iRun][iCentrality]>>gNetChargeRms[iRun][iCentrality]>>nEvents[iRun][iCentrality];
      cout<<nuDyn[iRun][iCentrality]<<" "<<nEvents[iRun][iCentrality]<<endl;
    }
    inputAscii20.close();
  }

  Int_t nRunCounter = 0;
  for(Int_t iCentrality = 0; iCentrality < nCentralityBins; iCentrality++) {
    for(Int_t iRun = 0; iRun < nRuns; iRun++) {
      nuDynValues[iCentrality] += nuDyn[iRun][iCentrality]*nEvents[iRun][iCentrality];
      netChargeValues[iCentrality] += gNetCharge[iRun][iCentrality]*nEvents[iRun][iCentrality];
      netChargeRmsValues[iCentrality] += gNetChargeRms[iRun][iCentrality]*nEvents[iRun][iCentrality];
      nEventsValues[iCentrality] += nEvents[iRun][iCentrality];
    }

    if(nEventsValues[iCentrality] != 0) {
      nRunCounter += 1;

      netChargeValues[iCentrality] /= nEventsValues[iCentrality];
      netChargeRmsValues[iCentrality] /= nEventsValues[iCentrality];
      nuDynValues[iCentrality] /= nEventsValues[iCentrality];
    }
    else {
      netChargeValues[iCentrality] = 999.;
      netChargeRmsValues[iCentrality] = 999.;
      nuDynValues[iCentrality] = 999.;
    }

    gNetChargeALICEData20[iCentrality] = netChargeValues[iCentrality];
    gNetChargeRmsALICEData20[iCentrality] = netChargeRmsValues[iCentrality];
    gNuDynALICEData20[iCentrality] = nuDynValues[iCentrality];
  }

  for(Int_t iCentrality = 0; iCentrality < nCentralityBins; iCentrality++) {
    for(Int_t iRun = 0; iRun < nRuns; iRun++) {
      netChargeValuesError[iCentrality] += TMath::Power((netChargeValues[iCentrality]-gNetCharge[iRun][iCentrality]),2);
      netChargeRmsValuesError[iCentrality] += TMath::Power((netChargeRmsValues[iCentrality]-gNetChargeRms[iRun][iCentrality]),2);
      nuDynValuesError[iCentrality] += TMath::Power((nuDynValues[iCentrality]-nuDyn[iRun][iCentrality]),2);
    }
    gNetChargeALICEData20Error[iCentrality] = TMath::Sqrt(netChargeValuesError[iCentrality]/(nRunCounter*(nRunCounter-1)));
    gNetChargeRmsALICEData20Error[iCentrality] = TMath::Sqrt(netChargeRmsValuesError[iCentrality]/(nRunCounter*(nRunCounter-1)));
    gNuDynALICEData20Error[iCentrality] = TMath::Sqrt(nuDynValuesError[iCentrality]/(nRunCounter*(nRunCounter-1)));
    Printf("Centrality: %d - nu_dyn: %lf +- %lf - Net charge: %lf +- %lf - RMS: %lf - %lf",iCentrality+1,
	   gNuDynALICEData20[iCentrality],
	   gNuDynALICEData20Error[iCentrality],
	   gNetChargeALICEData20[iCentrality],
	   gNetChargeALICEData20Error[iCentrality],
	   gNetChargeRmsALICEData20[iCentrality],
	   gNetChargeRmsALICEData20Error[iCentrality]);
  }

  //10 centrality bins
  nRunCounter = 0;
  for(Int_t iRun = 0; iRun < nRuns; iRun++) {
    for(Int_t iCentrality = 0; iCentrality < nCentralityBins; iCentrality++) {
      nuDyn[iRun][iCentrality] = 0.;
      gNetCharge[iRun][iCentrality] = 0.;
      gNetChargeRms[iRun][iCentrality] = 0.;
      nEvents[iRun][iCentrality] = 0.;
      nuDynValues[iCentrality] = 0.;
      netChargeValues[iCentrality] = 0.;
      netChargeRmsValues[iCentrality] = 0.;
      nEventsValues[iCentrality] = 0.;
      netChargeValuesError[iCentrality] = 0.;
      netChargeRmsValuesError[iCentrality] = 0.;
      nuDynValuesError[iCentrality] = 0.;
    }
  }

  for(Int_t iRun = 0; iRun < nRuns; iRun++) {
    ifstream inputAscii10;
    Printf("Adding run %d",nRunNumbers[iRun]);
    inputFileName = resultPath2;
    inputFileName += "/output."; inputFileName += nRunNumbers[iRun];
    inputFileName += ".txt";

    Printf("Filename: %s",inputFileName.Data());
    inputAscii10.open(inputFileName.Data());
    for(Int_t iCentrality = 0; iCentrality < nCentralityBins; iCentrality++) {
      inputAscii10>>gCentrality>>nuDyn[iRun][iCentrality]>>gNetCharge[iRun][iCentrality]>>gNetChargeRms[iRun][iCentrality]>>nEvents[iRun][iCentrality];
      //cout<<nuDyn[iRun][iCentrality]<<" "<<nEvents[iRun][iCentrality]<<endl;
    }
    inputAscii10.close();
  }

  Int_t nRunCounter = 0;
  for(Int_t iCentrality = 0; iCentrality < nCentralityBins; iCentrality++) {
    for(Int_t iRun = 0; iRun < nRuns; iRun++) {
      nuDynValues[iCentrality] += nuDyn[iRun][iCentrality]*nEvents[iRun][iCentrality];
      netChargeValues[iCentrality] += gNetCharge[iRun][iCentrality]*nEvents[iRun][iCentrality];
      netChargeRmsValues[iCentrality] += gNetChargeRms[iRun][iCentrality]*nEvents[iRun][iCentrality];
      nEventsValues[iCentrality] += nEvents[iRun][iCentrality];
    }

    if(nEventsValues[iCentrality] != 0) {
      nRunCounter += 1;

      netChargeValues[iCentrality] /= nEventsValues[iCentrality];
      netChargeRmsValues[iCentrality] /= nEventsValues[iCentrality];
      nuDynValues[iCentrality] /= nEventsValues[iCentrality];
    }
    else {
      netChargeValues[iCentrality] = 999.;
      netChargeRmsValues[iCentrality] = 999.;
      nuDynValues[iCentrality] = 999.;
    }

    gNetChargeALICEData10[iCentrality] = netChargeValues[iCentrality];
    gNetChargeRmsALICEData10[iCentrality] = netChargeRmsValues[iCentrality];
    gNuDynALICEData10[iCentrality] = nuDynValues[iCentrality];
  }

  for(Int_t iCentrality = 0; iCentrality < nCentralityBins; iCentrality++) {
    for(Int_t iRun = 0; iRun < nRuns; iRun++) {
      netChargeValuesError[iCentrality] += TMath::Power((netChargeValues[iCentrality]-gNetCharge[iRun][iCentrality]),2);
      netChargeRmsValuesError[iCentrality] += TMath::Power((netChargeRmsValues[iCentrality]-gNetChargeRms[iRun][iCentrality]),2);
      nuDynValuesError[iCentrality] += TMath::Power((nuDynValues[iCentrality]-nuDyn[iRun][iCentrality]),2);
    }
    gNetChargeALICEData10Error[iCentrality] = TMath::Sqrt(netChargeValuesError[iCentrality]/(nRunCounter*(nRunCounter-1)));
    gNetChargeRmsALICEData10Error[iCentrality] = TMath::Sqrt(netChargeRmsValuesError[iCentrality]/(nRunCounter*(nRunCounter-1)));
    gNuDynALICEData10Error[iCentrality] = TMath::Sqrt(nuDynValuesError[iCentrality]/(nRunCounter*(nRunCounter-1)));
    Printf("Centrality: %d - nu_dyn: %lf +- %lf - Net charge: %lf +- %lf - RMS: %lf - %lf",iCentrality+1,
	   gNuDynALICEData10[iCentrality],
	   gNuDynALICEData10Error[iCentrality],
	   gNetChargeALICEData10[iCentrality],
	   gNetChargeALICEData10Error[iCentrality],
	   gNetChargeRmsALICEData10[iCentrality],
	   gNetChargeRmsALICEData10Error[iCentrality]);
  }

}


//_______________________________________________________________//
void DrawMarker(Double_t x, Double_t y, Int_t style, 
		Double_t size, Int_t color) {
  TMarker *m = new TMarker(x,y,style);
  m->SetMarkerSize(size);
  m->SetMarkerColor(color);
  m->Draw();
}
 compareCentralitySets.C:1
 compareCentralitySets.C:2
 compareCentralitySets.C:3
 compareCentralitySets.C:4
 compareCentralitySets.C:5
 compareCentralitySets.C:6
 compareCentralitySets.C:7
 compareCentralitySets.C:8
 compareCentralitySets.C:9
 compareCentralitySets.C:10
 compareCentralitySets.C:11
 compareCentralitySets.C:12
 compareCentralitySets.C:13
 compareCentralitySets.C:14
 compareCentralitySets.C:15
 compareCentralitySets.C:16
 compareCentralitySets.C:17
 compareCentralitySets.C:18
 compareCentralitySets.C:19
 compareCentralitySets.C:20
 compareCentralitySets.C:21
 compareCentralitySets.C:22
 compareCentralitySets.C:23
 compareCentralitySets.C:24
 compareCentralitySets.C:25
 compareCentralitySets.C:26
 compareCentralitySets.C:27
 compareCentralitySets.C:28
 compareCentralitySets.C:29
 compareCentralitySets.C:30
 compareCentralitySets.C:31
 compareCentralitySets.C:32
 compareCentralitySets.C:33
 compareCentralitySets.C:34
 compareCentralitySets.C:35
 compareCentralitySets.C:36
 compareCentralitySets.C:37
 compareCentralitySets.C:38
 compareCentralitySets.C:39
 compareCentralitySets.C:40
 compareCentralitySets.C:41
 compareCentralitySets.C:42
 compareCentralitySets.C:43
 compareCentralitySets.C:44
 compareCentralitySets.C:45
 compareCentralitySets.C:46
 compareCentralitySets.C:47
 compareCentralitySets.C:48
 compareCentralitySets.C:49
 compareCentralitySets.C:50
 compareCentralitySets.C:51
 compareCentralitySets.C:52
 compareCentralitySets.C:53
 compareCentralitySets.C:54
 compareCentralitySets.C:55
 compareCentralitySets.C:56
 compareCentralitySets.C:57
 compareCentralitySets.C:58
 compareCentralitySets.C:59
 compareCentralitySets.C:60
 compareCentralitySets.C:61
 compareCentralitySets.C:62
 compareCentralitySets.C:63
 compareCentralitySets.C:64
 compareCentralitySets.C:65
 compareCentralitySets.C:66
 compareCentralitySets.C:67
 compareCentralitySets.C:68
 compareCentralitySets.C:69
 compareCentralitySets.C:70
 compareCentralitySets.C:71
 compareCentralitySets.C:72
 compareCentralitySets.C:73
 compareCentralitySets.C:74
 compareCentralitySets.C:75
 compareCentralitySets.C:76
 compareCentralitySets.C:77
 compareCentralitySets.C:78
 compareCentralitySets.C:79
 compareCentralitySets.C:80
 compareCentralitySets.C:81
 compareCentralitySets.C:82
 compareCentralitySets.C:83
 compareCentralitySets.C:84
 compareCentralitySets.C:85
 compareCentralitySets.C:86
 compareCentralitySets.C:87
 compareCentralitySets.C:88
 compareCentralitySets.C:89
 compareCentralitySets.C:90
 compareCentralitySets.C:91
 compareCentralitySets.C:92
 compareCentralitySets.C:93
 compareCentralitySets.C:94
 compareCentralitySets.C:95
 compareCentralitySets.C:96
 compareCentralitySets.C:97
 compareCentralitySets.C:98
 compareCentralitySets.C:99
 compareCentralitySets.C:100
 compareCentralitySets.C:101
 compareCentralitySets.C:102
 compareCentralitySets.C:103
 compareCentralitySets.C:104
 compareCentralitySets.C:105
 compareCentralitySets.C:106
 compareCentralitySets.C:107
 compareCentralitySets.C:108
 compareCentralitySets.C:109
 compareCentralitySets.C:110
 compareCentralitySets.C:111
 compareCentralitySets.C:112
 compareCentralitySets.C:113
 compareCentralitySets.C:114
 compareCentralitySets.C:115
 compareCentralitySets.C:116
 compareCentralitySets.C:117
 compareCentralitySets.C:118
 compareCentralitySets.C:119
 compareCentralitySets.C:120
 compareCentralitySets.C:121
 compareCentralitySets.C:122
 compareCentralitySets.C:123
 compareCentralitySets.C:124
 compareCentralitySets.C:125
 compareCentralitySets.C:126
 compareCentralitySets.C:127
 compareCentralitySets.C:128
 compareCentralitySets.C:129
 compareCentralitySets.C:130
 compareCentralitySets.C:131
 compareCentralitySets.C:132
 compareCentralitySets.C:133
 compareCentralitySets.C:134
 compareCentralitySets.C:135
 compareCentralitySets.C:136
 compareCentralitySets.C:137
 compareCentralitySets.C:138
 compareCentralitySets.C:139
 compareCentralitySets.C:140
 compareCentralitySets.C:141
 compareCentralitySets.C:142
 compareCentralitySets.C:143
 compareCentralitySets.C:144
 compareCentralitySets.C:145
 compareCentralitySets.C:146
 compareCentralitySets.C:147
 compareCentralitySets.C:148
 compareCentralitySets.C:149
 compareCentralitySets.C:150
 compareCentralitySets.C:151
 compareCentralitySets.C:152
 compareCentralitySets.C:153
 compareCentralitySets.C:154
 compareCentralitySets.C:155
 compareCentralitySets.C:156
 compareCentralitySets.C:157
 compareCentralitySets.C:158
 compareCentralitySets.C:159
 compareCentralitySets.C:160
 compareCentralitySets.C:161
 compareCentralitySets.C:162
 compareCentralitySets.C:163
 compareCentralitySets.C:164
 compareCentralitySets.C:165
 compareCentralitySets.C:166
 compareCentralitySets.C:167
 compareCentralitySets.C:168
 compareCentralitySets.C:169
 compareCentralitySets.C:170
 compareCentralitySets.C:171
 compareCentralitySets.C:172
 compareCentralitySets.C:173
 compareCentralitySets.C:174
 compareCentralitySets.C:175
 compareCentralitySets.C:176
 compareCentralitySets.C:177
 compareCentralitySets.C:178
 compareCentralitySets.C:179
 compareCentralitySets.C:180
 compareCentralitySets.C:181
 compareCentralitySets.C:182
 compareCentralitySets.C:183
 compareCentralitySets.C:184
 compareCentralitySets.C:185
 compareCentralitySets.C:186
 compareCentralitySets.C:187
 compareCentralitySets.C:188
 compareCentralitySets.C:189
 compareCentralitySets.C:190
 compareCentralitySets.C:191
 compareCentralitySets.C:192
 compareCentralitySets.C:193
 compareCentralitySets.C:194
 compareCentralitySets.C:195
 compareCentralitySets.C:196
 compareCentralitySets.C:197
 compareCentralitySets.C:198
 compareCentralitySets.C:199
 compareCentralitySets.C:200
 compareCentralitySets.C:201
 compareCentralitySets.C:202
 compareCentralitySets.C:203
 compareCentralitySets.C:204
 compareCentralitySets.C:205
 compareCentralitySets.C:206
 compareCentralitySets.C:207
 compareCentralitySets.C:208
 compareCentralitySets.C:209
 compareCentralitySets.C:210
 compareCentralitySets.C:211
 compareCentralitySets.C:212
 compareCentralitySets.C:213
 compareCentralitySets.C:214
 compareCentralitySets.C:215
 compareCentralitySets.C:216
 compareCentralitySets.C:217
 compareCentralitySets.C:218
 compareCentralitySets.C:219
 compareCentralitySets.C:220
 compareCentralitySets.C:221
 compareCentralitySets.C:222
 compareCentralitySets.C:223
 compareCentralitySets.C:224
 compareCentralitySets.C:225
 compareCentralitySets.C:226
 compareCentralitySets.C:227
 compareCentralitySets.C:228
 compareCentralitySets.C:229
 compareCentralitySets.C:230
 compareCentralitySets.C:231
 compareCentralitySets.C:232
 compareCentralitySets.C:233
 compareCentralitySets.C:234
 compareCentralitySets.C:235
 compareCentralitySets.C:236
 compareCentralitySets.C:237
 compareCentralitySets.C:238
 compareCentralitySets.C:239
 compareCentralitySets.C:240
 compareCentralitySets.C:241
 compareCentralitySets.C:242
 compareCentralitySets.C:243
 compareCentralitySets.C:244
 compareCentralitySets.C:245
 compareCentralitySets.C:246
 compareCentralitySets.C:247
 compareCentralitySets.C:248
 compareCentralitySets.C:249
 compareCentralitySets.C:250
 compareCentralitySets.C:251
 compareCentralitySets.C:252
 compareCentralitySets.C:253
 compareCentralitySets.C:254
 compareCentralitySets.C:255
 compareCentralitySets.C:256
 compareCentralitySets.C:257
 compareCentralitySets.C:258
 compareCentralitySets.C:259
 compareCentralitySets.C:260
 compareCentralitySets.C:261
 compareCentralitySets.C:262
 compareCentralitySets.C:263
 compareCentralitySets.C:264
 compareCentralitySets.C:265
 compareCentralitySets.C:266
 compareCentralitySets.C:267
 compareCentralitySets.C:268
 compareCentralitySets.C:269
 compareCentralitySets.C:270
 compareCentralitySets.C:271
 compareCentralitySets.C:272
 compareCentralitySets.C:273
 compareCentralitySets.C:274
 compareCentralitySets.C:275
 compareCentralitySets.C:276
 compareCentralitySets.C:277
 compareCentralitySets.C:278
 compareCentralitySets.C:279
 compareCentralitySets.C:280
 compareCentralitySets.C:281
 compareCentralitySets.C:282
 compareCentralitySets.C:283
 compareCentralitySets.C:284
 compareCentralitySets.C:285
 compareCentralitySets.C:286
 compareCentralitySets.C:287
 compareCentralitySets.C:288
 compareCentralitySets.C:289
 compareCentralitySets.C:290
 compareCentralitySets.C:291
 compareCentralitySets.C:292
 compareCentralitySets.C:293
 compareCentralitySets.C:294
 compareCentralitySets.C:295
 compareCentralitySets.C:296
 compareCentralitySets.C:297
 compareCentralitySets.C:298
 compareCentralitySets.C:299
 compareCentralitySets.C:300
 compareCentralitySets.C:301
 compareCentralitySets.C:302
 compareCentralitySets.C:303
 compareCentralitySets.C:304
 compareCentralitySets.C:305
 compareCentralitySets.C:306
 compareCentralitySets.C:307
 compareCentralitySets.C:308
 compareCentralitySets.C:309
 compareCentralitySets.C:310
 compareCentralitySets.C:311
 compareCentralitySets.C:312
 compareCentralitySets.C:313
 compareCentralitySets.C:314
 compareCentralitySets.C:315
 compareCentralitySets.C:316
 compareCentralitySets.C:317
 compareCentralitySets.C:318
 compareCentralitySets.C:319
 compareCentralitySets.C:320
 compareCentralitySets.C:321
 compareCentralitySets.C:322
 compareCentralitySets.C:323
 compareCentralitySets.C:324
 compareCentralitySets.C:325
 compareCentralitySets.C:326
 compareCentralitySets.C:327
 compareCentralitySets.C:328
 compareCentralitySets.C:329
 compareCentralitySets.C:330
 compareCentralitySets.C:331
 compareCentralitySets.C:332
 compareCentralitySets.C:333
 compareCentralitySets.C:334
 compareCentralitySets.C:335
 compareCentralitySets.C:336
 compareCentralitySets.C:337
 compareCentralitySets.C:338
 compareCentralitySets.C:339
 compareCentralitySets.C:340
 compareCentralitySets.C:341
 compareCentralitySets.C:342
 compareCentralitySets.C:343
 compareCentralitySets.C:344
 compareCentralitySets.C:345
 compareCentralitySets.C:346
 compareCentralitySets.C:347
 compareCentralitySets.C:348
 compareCentralitySets.C:349
 compareCentralitySets.C:350
 compareCentralitySets.C:351