|
|
|
#include <comparisonPlots.h>
Public Member Functions | |
| comparisonPlots (std::string filename, std::string outputDir, std::string outputFilename="OUTPUT_comparison.root") | |
| void | plot3x5 (TCut Cut, char *dirName, bool savePlot=false, std::string plotName="plot3x5.eps", bool autolimits=false, int ColorCode=1) |
| void | plot3x5Profile (TCut Cut, char *dirName, int nBins, bool savePlot=false, std::string plotName="plot3x5Profile.eps", bool autolimits=false, int ColorCode=1) |
| void | plotTwist (TCut Cut, char *dirName, bool savePlot=false, std::string plotName="plot3x5.eps", bool autolimits=false, int ColorCode=1) |
| void | readTree () |
| void | Write () |
Public Attributes | |
| float | arrowSize |
| TTree * | data |
| TFile * | fin |
| TFile * | output |
Private Member Functions | |
| void | getHistMaxMin (TH1 *hist, double &max, double &min, int flag) |
| void | getMaxMin () |
Private Attributes | |
| std::string | _outputDir |
| float | alpha_ |
| float | beta_ |
| float | dalpha_ |
| float | dbeta_ |
| int | detDim_ |
| float | dgamma_ |
| float | dphi_ |
| float | dr_ |
| float | dx_ |
| float | dy_ |
| float | dz_ |
| float | eta_ |
| float | gamma_ |
| int | id_ |
| int | level_ |
| float | maxDR |
| float | maxDX |
| float | maxDY |
| float | maxDZ |
| float | maxPhi |
| float | maxR |
| float | maxRDPhi |
| float | maxZ |
| int | mid_ |
| float | minDR |
| float | minDX |
| float | minDY |
| float | minDZ |
| float | minPhi |
| float | minR |
| float | minRDPhi |
| float | minZ |
| int | mlevel_ |
| float | phi_ |
| float | r_ |
| int | sublevel_ |
| int | useDetId_ |
| float | x_ |
| float | y_ |
| float | z_ |
Definition at line 24 of file comparisonPlots.h.
| comparisonPlots::comparisonPlots | ( | std::string | filename, |
| std::string | outputDir, | ||
| std::string | outputFilename = "OUTPUT_comparison.root" |
||
| ) |
Definition at line 16 of file comparisonPlots.cc.
References _outputDir, fin, output, indexGen::outputDir, and readTree().
| void comparisonPlots::getHistMaxMin | ( | TH1 * | hist, |
| double & | max, | ||
| double & | min, | ||
| int | flag | ||
| ) | [private] |
Definition at line 1137 of file comparisonPlots.cc.
References i.
Referenced by plot3x5(), plot3x5Profile(), and plotTwist().
{
int nBins = hist->GetNbinsX();
for (int i = 0; i < nBins; ++i){
double binContent = hist->GetBinContent(i);
if (binContent > 0){
//double binWidth = hist->GetBinLowEdge(i) - hist->GetBinLowEdge(i-1);
//std::cout << "bin width1: " << hist->GetBinWidth(i) << ", bin width2: " << binWidth << std::endl;
if (flag == 0) max = hist->GetBinLowEdge(i) + 2.*hist->GetBinWidth(i) + 0.1*hist->GetBinLowEdge(i);
if (flag == 1) max = hist->GetBinLowEdge(i) + hist->GetBinWidth(i) + 0.1*hist->GetBinLowEdge(i);
}
}
for (int i = (nBins-1); i >= 0; i--){
double binContent = hist->GetBinContent(i);
if (binContent > 0) min = hist->GetBinLowEdge(i) - 0.1*hist->GetBinLowEdge(i);
}
//std::cout << "max: " << max << ", min: " << min << std::endl;
}
| void comparisonPlots::getMaxMin | ( | ) | [private] |
Definition at line 1101 of file comparisonPlots.cc.
References data, dphi_, dr_, dx_, dy_, dz_, i, maxDR, maxDX, maxDY, maxDZ, maxPhi, maxR, maxRDPhi, maxZ, minDR, minDX, minDY, minDZ, minPhi, minR, minRDPhi, minZ, phi_, r_, and z_.
{
data->GetEntry(0);
minR = r_; maxR = r_;
minZ = z_; maxZ = z_;
minPhi = phi_; maxPhi = phi_;
minDR = dr_; maxDR = dr_;
minDZ = dz_; maxDZ = dz_;
minRDPhi = r_*dphi_; maxRDPhi = r_*dphi_;
minDX = dx_; maxDX = dx_;
minDY = dy_; maxDY = dy_;
int nEntries = data->GetEntries();
for (int i = 1; i < nEntries; ++i){
data->GetEntry(i);
if (r_ < minR) minR = r_;
if (r_ > maxR) maxR = r_;
if (z_ < minZ) minZ = z_;
if (z_ > maxZ) maxZ = z_;
if (phi_ < minPhi) minPhi = phi_;
if (phi_ > maxPhi) maxPhi = phi_;
if (dr_ < minDR) minDR = dr_;
if (dr_ > maxDR) maxDR = dr_;
if (dz_ < minDZ) minDZ = dz_;
if (dz_ > maxDZ) maxDZ = dz_;
if (r_*dphi_ < minRDPhi) minRDPhi = r_*dphi_;
if (r_*dphi_ > maxRDPhi) maxRDPhi = r_*dphi_;
if (dx_ < minDX) minDX = dx_;
if (dx_ > maxDX) maxDX = dx_;
if (dy_ < minDY) minDY = dy_;
if (dy_ > maxDY) maxDY = dy_;
}
}
| void comparisonPlots::plot3x5 | ( | TCut | Cut, |
| char * | dirName, | ||
| bool | savePlot = false, |
||
| std::string | plotName = "plot3x5.eps", |
||
| bool | autolimits = false, |
||
| int | ColorCode = 1 |
||
| ) |
Definition at line 149 of file comparisonPlots.cc.
References _outputDir, trackerHits::c, data, getHistMaxMin(), i, j, and output.
{
//int ColorCodeHalfBarrels = 1; //color seperation of moodules corresponding to the different subdets
//if you want to seperate corresponding to the halfbarrels (+z or -z) set colorcode !=1
gStyle->SetTitleAlign(22);
gStyle->SetTitleX(0.5);
gStyle->SetTitleY(0.97);
//gStyle->SetShadowColor(0);
//gStyle->SetTitleSize(.5);
gStyle->SetTitleFont(62);
// --------- create directory for histograms ---------
//const char* dirName = Cut;
TDirectory* plotDir = output->mkdir( dirName );
// --------- get right limits for histogram ---------
double minimumR, maximumR;
double minimumZ, maximumZ;
double minimumPhi, maximumPhi;
double minimumDR, maximumDR;
double minimumDZ, maximumDZ;
double minimumRDPhi, maximumRDPhi;
double minimumDX, maximumDX;
double minimumDY, maximumDY;
if (autolimits){
// --------- get right limits for histogram ---------
TH1F* phr = new TH1F("phr", "phr", 200, 0, 200);
TH1F* phz = new TH1F("phz", "phz", 400, -300, 300);
TH1F* phphi = new TH1F("phphi", "phphi", 200, -3.15, 3.15);
TH1F* phdr = new TH1F("phdr", "phdr", 2000, -10, 10);
TH1F* phdz = new TH1F("phdz", "phdz", 2000, -10, 10);
TH1F* phrdphi = new TH1F("phrdphi", "phrdphi", 200, -10, 10);
TH1F* phdx = new TH1F("phdx", "phy", 2000, -10, 10);
TH1F* phdy = new TH1F("phdy", "phy", 2000, -10, 10);
data->Project("phr","r",Cut);
data->Project("phz","z",Cut);
data->Project("phphi","phi",Cut);
data->Project("phdr","dr",Cut);
data->Project("phdz","dz",Cut);
data->Project("phrdphi","r*dphi",Cut);
data->Project("phdx","dx",Cut);
data->Project("phdy","dy",Cut);
getHistMaxMin(phr, maximumR, minimumR, 0);
getHistMaxMin(phz, maximumZ, minimumZ, 0);
getHistMaxMin(phphi, maximumPhi, minimumPhi, 0);
getHistMaxMin(phdr, maximumDR, minimumDR, 1);
getHistMaxMin(phdz, maximumDZ, minimumDZ, 1);
getHistMaxMin(phrdphi, maximumRDPhi, minimumRDPhi, 1);
getHistMaxMin(phdx, maximumDX, minimumDX, 1);
getHistMaxMin(phdy, maximumDY, minimumDY, 1);
}
else{
minimumR = 0., maximumR = 200.;
minimumZ = -300., maximumZ = 300.;
minimumPhi = -3.15, maximumPhi = 3.15;
minimumDR = -1, maximumDR = 1;
minimumDZ = -1, maximumDZ = 1;
minimumRDPhi = -1, maximumRDPhi = 1;
minimumDX = -1, maximumDX = 1;
minimumDY = -1, maximumDY = 1;
}
// --------- declare histograms ---------
TH1F* h_dr = new TH1F("h_dr", "#Delta r", 2000, minimumDR, maximumDR);
TH1F* h_dz = new TH1F("h_dz", "#Delta z", 2000, minimumDZ, maximumDZ);
TH1F* h_rdphi = new TH1F("h_rdphi", "r* #Delta #phi", 2000, minimumRDPhi, maximumRDPhi);
TH2F* h_drVr = new TH2F("h_drVr","#Delta r vs. r",200,minimumR,maximumR,200,minimumDR,maximumDR);
TH2F* h_dzVr = new TH2F("h_dzVr","#Delta z vs. r",200,minimumR,maximumR,200,minimumDZ,maximumDZ);
TH2F* h_rdphiVr = new TH2F("h_rdphiVr","r#Delta #phi vs. r",200,minimumR,maximumR,200,minimumRDPhi,maximumRDPhi);
TH2F* h_dxVr = new TH2F("h_dxVr","#Delta x vs. r", 200,minimumR,maximumR, 200,minimumDX,maximumDX);
TH2F* h_dyVr = new TH2F("h_dyVr","#Delta y vs. r", 200,minimumR,maximumR, 200,minimumDY,maximumDY);
TH2F* h_drVz = new TH2F("h_drVz","#Delta r vs. z", 200,minimumZ,maximumZ, 200,minimumDR,maximumDR);
TH2F* h_dzVz = new TH2F("h_dzVz","#Delta z vs. z", 200,minimumZ,maximumZ, 200,minimumDZ,maximumDZ);
TH2F* h_rdphiVz = new TH2F("h_rdphiVz","r#Delta #phi vs. z", 200,minimumZ,maximumZ, 200,minimumRDPhi,maximumRDPhi);
TH2F* h_dxVz = new TH2F("h_dxVz","#Delta x vs. z", 200,minimumZ,maximumZ, 200,minimumDX,maximumDX);
TH2F* h_dyVz = new TH2F("h_dyVz","#Delta y vs. z", 200,minimumZ,maximumZ, 200,minimumDY,maximumDY);
TH2F* h_drVphi = new TH2F("h_drVphi","#Delta r vs. #phi", 200,minimumPhi,maximumPhi,200,minimumDR,maximumDR);
TH2F* h_dzVphi = new TH2F("h_dzVphi","#Delta z vs. #phi", 200,minimumPhi,maximumPhi, 200,minimumDZ,maximumDZ);
TH2F* h_rdphiVphi = new TH2F("h_rdphiVphi","r#Delta #phi vs. #phi", 200,minimumPhi,maximumPhi,200,minimumRDPhi,maximumRDPhi);
TH2F* h_dxVphi = new TH2F("h_dxVphi","#Delta x vs. #phi", 200,minimumPhi,maximumPhi, 200,minimumDX,maximumDX);
TH2F* h_dyVphi = new TH2F("h_dyVphi","#Delta y vs. #phi", 200,minimumPhi,maximumPhi, 200,minimumDY,maximumDY);
h_drVr->SetMarkerStyle(6);
h_dzVr->SetMarkerStyle(6);
h_rdphiVr->SetMarkerStyle(6);
h_dxVr->SetMarkerStyle(6);
h_dyVr->SetMarkerStyle(6);
h_drVz->SetMarkerStyle(6);
h_dzVz->SetMarkerStyle(6);
h_rdphiVz->SetMarkerStyle(6);
h_dxVz->SetMarkerStyle(6);
h_dyVz->SetMarkerStyle(6);
h_drVphi->SetMarkerStyle(6);
h_dzVphi->SetMarkerStyle(6);
h_rdphiVphi->SetMarkerStyle(6);
h_dxVphi->SetMarkerStyle(6);
h_dyVphi->SetMarkerStyle(6);
int j = 2;
TCut zCut[6];
zCut[0] = "z < 0";
zCut[1] = "z >= 0";
if (ColorCode == 1) {
j = 6;
// zCut[6];
zCut[0] = "sublevel==1";
zCut[1] = "sublevel==2";
zCut[2] = "sublevel==3";
zCut[3] = "sublevel==4";
zCut[4] = "sublevel==5";
zCut[5] = "sublevel==6";
}
// --------- project tree onto histograms ---------
data->Project("h_dr","dr",Cut);
data->Project("h_dz","dz",Cut);
data->Project("h_rdphi","r*dphi",Cut);
data->Project("h_drVr", "dr:r",Cut);
TGraph* gr_drVr_Array[j];
TMultiGraph* mgr_drVr=new TMultiGraph();
for ( int i = 0; i < j; i++) {
data->Draw("dr:r",Cut+zCut[i]);
gr_drVr_Array[i] = new TGraph(data->GetSelectedRows(),data->GetV2(),data->GetV1());
gr_drVr_Array[i]->SetMarkerColor(int(i/4)+i+1);
gr_drVr_Array[i]->SetMarkerStyle(6);
mgr_drVr->Add(gr_drVr_Array[i],"p");
}
data->Project("h_dzVr", "dz:r",Cut);
TGraph* gr_dzVr_Array[j];
TMultiGraph* mgr_dzVr=new TMultiGraph();
for ( int i = 0; i < j; i++) {
data->Draw("dz:r",Cut+zCut[i]);
gr_dzVr_Array[i] = new TGraph(data->GetSelectedRows(),data->GetV2(),data->GetV1());
gr_dzVr_Array[i]->SetMarkerColor(int(i/4)+i+1);
gr_dzVr_Array[i]->SetMarkerStyle(6);
mgr_dzVr->Add(gr_dzVr_Array[i],"p");
}
data->Project("h_rdphiVr", "r*dphi:r",Cut);
TGraph* gr_rdphiVr_Array[j];
TMultiGraph* mgr_rdphiVr=new TMultiGraph();
for ( int i = 0; i < j; i++) {
data->Draw("r*dphi:r",Cut+zCut[i]);
gr_rdphiVr_Array[i] = new TGraph(data->GetSelectedRows(),data->GetV2(),data->GetV1());
gr_rdphiVr_Array[i]->SetMarkerColor(int(i/4)+i+1);
gr_rdphiVr_Array[i]->SetMarkerStyle(6);
mgr_rdphiVr->Add(gr_rdphiVr_Array[i],"p");
}
data->Project("h_dxVr", "dx:r",Cut);
TGraph* gr_dxVr_Array[j];
TMultiGraph* mgr_dxVr=new TMultiGraph();
for ( int i = 0; i < j; i++) {
data->Draw("dx:r",Cut+zCut[i]);
gr_dxVr_Array[i] = new TGraph(data->GetSelectedRows(),data->GetV2(),data->GetV1());
gr_dxVr_Array[i]->SetMarkerColor(int(i/4)+i+1);
gr_dxVr_Array[i]->SetMarkerStyle(6);
mgr_dxVr->Add(gr_dxVr_Array[i],"p");
}
data->Project("h_dyVr", "dy:r",Cut);
TGraph* gr_dyVr_Array[j];
TMultiGraph* mgr_dyVr=new TMultiGraph();
for ( int i = 0; i < j; i++) {
data->Draw("dy:r",Cut+zCut[i]);
gr_dyVr_Array[i] = new TGraph(data->GetSelectedRows(),data->GetV2(),data->GetV1());
gr_dyVr_Array[i]->SetMarkerColor(int(i/4)+i+1);
gr_dyVr_Array[i]->SetMarkerStyle(6);
mgr_dyVr->Add(gr_dyVr_Array[i],"p");
}
data->Project("h_drVz", "dr:z",Cut);
TGraph* gr_drVz_Array[j];
TMultiGraph* mgr_drVz=new TMultiGraph();
for ( int i = 0; i < j; i++) {
data->Draw("dr:z",Cut+zCut[i]);
gr_drVz_Array[i] = new TGraph(data->GetSelectedRows(),data->GetV2(),data->GetV1());
gr_drVz_Array[i]->SetMarkerColor(int(i/4)+i+1);
gr_drVz_Array[i]->SetMarkerStyle(6);
mgr_drVz->Add(gr_drVz_Array[i],"p");
}
data->Project("h_dzVz", "dz:z",Cut);
TGraph* gr_dzVz_Array[j];
TMultiGraph* mgr_dzVz=new TMultiGraph();
for ( int i = 0; i < j; i++) {
data->Draw("dz:z",Cut+zCut[i]);
gr_dzVz_Array[i] = new TGraph(data->GetSelectedRows(),data->GetV2(),data->GetV1());
gr_dzVz_Array[i]->SetMarkerColor(int(i/4)+i+1);
gr_dzVz_Array[i]->SetMarkerStyle(6);
mgr_dzVz->Add(gr_dzVz_Array[i],"p");
}
data->Project("h_rdphiVz", "r*dphi:z",Cut);
TGraph* gr_rdphiVz_Array[j];
TMultiGraph* mgr_rdphiVz=new TMultiGraph();
for ( int i = 0; i < j; i++) {
data->Draw("r*dphi:z",Cut+zCut[i]);
gr_rdphiVz_Array[i] = new TGraph(data->GetSelectedRows(),data->GetV2(),data->GetV1());
gr_rdphiVz_Array[i]->SetMarkerColor(int(i/4)+i+1);
gr_rdphiVz_Array[i]->SetMarkerStyle(6);
mgr_rdphiVz->Add(gr_rdphiVz_Array[i],"p");
}
data->Project("h_dxVz", "dx:z",Cut);
TGraph* gr_dxVz_Array[j];
TMultiGraph* mgr_dxVz=new TMultiGraph();
for ( int i = 0; i < j; i++) {
data->Draw("dx:z",Cut+zCut[i]);
gr_dxVz_Array[i] = new TGraph(data->GetSelectedRows(),data->GetV2(),data->GetV1());
gr_dxVz_Array[i]->SetMarkerColor(int(i/4)+i+1);
gr_dxVz_Array[i]->SetMarkerStyle(6);
mgr_dxVz->Add(gr_dxVz_Array[i],"p");
}
data->Project("h_dyVz", "dy:z",Cut);
TGraph* gr_dyVz_Array[j];
TMultiGraph* mgr_dyVz=new TMultiGraph();
for ( int i = 0; i < j; i++) {
data->Draw("dy:z",Cut+zCut[i]);
gr_dyVz_Array[i] = new TGraph(data->GetSelectedRows(),data->GetV2(),data->GetV1());
gr_dyVz_Array[i]->SetMarkerColor(int(i/4)+i+1);
gr_dyVz_Array[i]->SetMarkerStyle(6);
mgr_dyVz->Add(gr_dyVz_Array[i],"p");
}
data->Project("h_drVphi", "dr:phi",Cut);
TGraph* gr_drVphi_Array[j];
TMultiGraph* mgr_drVphi=new TMultiGraph();
for ( int i = 0; i < j; i++) {
data->Draw("dr:phi",Cut+zCut[i]);
gr_drVphi_Array[i] = new TGraph(data->GetSelectedRows(),data->GetV2(),data->GetV1());
gr_drVphi_Array[i]->SetMarkerColor(int(i/4)+i+1);
gr_drVphi_Array[i]->SetMarkerStyle(6);
mgr_drVphi->Add(gr_drVphi_Array[i],"p");
}
data->Project("h_dzVphi", "dz:phi",Cut);
TGraph* gr_dzVphi_Array[j];
TMultiGraph* mgr_dzVphi=new TMultiGraph();
for ( int i = 0; i < j; i++) {
data->Draw("dz:phi",Cut+zCut[i]);
gr_dzVphi_Array[i] = new TGraph(data->GetSelectedRows(),data->GetV2(),data->GetV1());
gr_dzVphi_Array[i]->SetMarkerColor(int(i/4)+i+1);
gr_dzVphi_Array[i]->SetMarkerStyle(6);
mgr_dzVphi->Add(gr_dzVphi_Array[i],"p");
}
data->Project("h_rdphiVphi", "r*dphi:phi",Cut);
TGraph* gr_rdphiVphi_Array[j];
TMultiGraph* mgr_rdphiVphi=new TMultiGraph();
for ( int i = 0; i < j; i++) {
data->Draw("r*dphi:phi",Cut+zCut[i]);
gr_rdphiVphi_Array[i] = new TGraph(data->GetSelectedRows(),data->GetV2(),data->GetV1());
gr_rdphiVphi_Array[i]->SetMarkerColor(int(i/4)+i+1);
gr_rdphiVphi_Array[i]->SetMarkerStyle(6);
mgr_rdphiVphi->Add(gr_rdphiVphi_Array[i],"p");
}
data->Project("h_dxVphi", "dx:phi",Cut);
TGraph* gr_dxVphi_Array[j];
TMultiGraph* mgr_dxVphi=new TMultiGraph();
for ( int i = 0; i < j; i++) {
data->Draw("dx:phi",Cut+zCut[i]);
gr_dxVphi_Array[i] = new TGraph(data->GetSelectedRows(),data->GetV2(),data->GetV1());
gr_dxVphi_Array[i]->SetMarkerColor(int(i/4)+i+1);
gr_dxVphi_Array[i]->SetMarkerStyle(6);
mgr_dxVphi->Add(gr_dxVphi_Array[i],"p");
}
data->Project("h_dyVphi", "dy:phi",Cut);
TGraph* gr_dyVphi_Array[j];
TMultiGraph* mgr_dyVphi=new TMultiGraph();
for ( int i = 0; i < j; i++) {
data->Draw("dy:phi",Cut+zCut[i]);
gr_dyVphi_Array[i] = new TGraph(data->GetSelectedRows(),data->GetV2(),data->GetV1());
gr_dyVphi_Array[i]->SetMarkerColor(int(i/4)+i+1);
gr_dyVphi_Array[i]->SetMarkerStyle(6);
mgr_dyVphi->Add(gr_dyVphi_Array[i],"p");
}
// --------- draw histograms ---------
TCanvas* c0 = new TCanvas("c0", "c0", 200, 10, 900, 300);
c0->SetFillColor(0);
c0->Divide(3,1);
c0->cd(1);
h_dr->Draw();
c0->cd(2);
h_dz->Draw();
c0->cd(3);
h_rdphi->Draw();
if (savePlot) c0->Print((_outputDir+"plot3x1_"+plotName).c_str());
mgr_drVr->SetTitle("#Delta r vs. r");
//TStyle::SetTitleAlign(23)
mgr_dzVr->SetTitle("#Delta z vs. r");
mgr_rdphiVr->SetTitle("r#Delta #phi vs. r");
mgr_dxVr->SetTitle("#Delta x vs. r");
mgr_dyVr->SetTitle("#Delta y vs. r");
mgr_drVz->SetTitle("#Delta r vs. z");
mgr_dzVz->SetTitle("#Delta z vs. z");
mgr_rdphiVz->SetTitle("r#Delta #phi vs. z");
mgr_dxVz->SetTitle("#Delta x vs. z");
mgr_dyVz->SetTitle("#Delta y vs. z");
mgr_drVphi->SetTitle("#Delta r vs. #phi");
mgr_dzVphi->SetTitle("#Delta z vs. #phi");
mgr_rdphiVphi->SetTitle("r#Delta #phi vs. #phi");
mgr_dxVphi->SetTitle("#Delta x vs. #phi");
mgr_dyVphi->SetTitle("#Delta y vs. #phi");
// --------- draw histograms ---------
TCanvas* c = new TCanvas("c", "c", 200, 10, 1200, 700);
c->SetFillColor(0);
data->SetMarkerSize(0.5);
data->SetMarkerStyle(6);
c->Divide(5,3);
c->cd(1);
if (!autolimits) h_drVr->Draw();
if (autolimits) mgr_drVr->Draw("a");
c->Update();
c->cd(2);
if (!autolimits) h_dzVr->Draw();
if (autolimits) mgr_dzVr->Draw("a");
c->Update();
c->cd(3);
if (!autolimits) h_rdphiVr->Draw();
if (autolimits) mgr_rdphiVr->Draw("a");
c->cd(4);
if (!autolimits) h_dxVr->Draw();
if (autolimits) mgr_dxVr->Draw("a");
c->cd(5);
if (!autolimits) h_dyVr->Draw();
if (autolimits) mgr_dyVr->Draw("a");
c->cd(6);
if (!autolimits) h_drVz->Draw();
if (autolimits) mgr_drVz->Draw("a");
c->cd(7);
if (!autolimits) h_dzVz->Draw();
if (autolimits) mgr_dzVz->Draw("a");
c->cd(8);
if (!autolimits) h_rdphiVz->Draw();
if (autolimits) mgr_rdphiVz->Draw("a");
c->cd(9);
if (!autolimits) h_dxVz->Draw();
if (autolimits) mgr_dxVz->Draw("a");
c->cd(10);
if (!autolimits) h_dyVz->Draw();
if (autolimits) mgr_dyVz->Draw("a");
c->cd(11);
if (!autolimits) h_drVphi->Draw();
if (autolimits) mgr_drVphi->Draw("a");
c->cd(12);
if (!autolimits) h_dzVphi->Draw();
if (autolimits) mgr_dzVphi->Draw("a");
c->cd(13);
if (!autolimits) h_rdphiVphi->Draw();
if (autolimits) mgr_rdphiVphi->Draw("a");
c->cd(14);
if (!autolimits) h_dxVphi->Draw();
if (autolimits) mgr_dxVphi->Draw("a");
c->cd(15);
if (!autolimits) h_dyVphi->Draw();
if (autolimits) mgr_dyVphi->Draw("a");
c->Update();
// --------- set output directory for histograms ---------
plotDir->cd();
h_dr->Write(); h_dz->Write(); h_rdphi->Write();
h_drVr->Write(); h_dzVr->Write(); h_rdphiVr->Write(); h_dxVr->Write(); h_dyVr->Write();
h_drVz->Write(); h_dzVz->Write(); h_rdphiVz->Write(); h_dxVz->Write(); h_dyVz->Write();
h_drVphi->Write(); h_dzVphi->Write(); h_rdphiVphi->Write(); h_dxVphi->Write(); h_dyVphi->Write();
if (savePlot) c->Print((_outputDir+"plot3x5_"+plotName).c_str());
}
| void comparisonPlots::plot3x5Profile | ( | TCut | Cut, |
| char * | dirName, | ||
| int | nBins, | ||
| bool | savePlot = false, |
||
| std::string | plotName = "plot3x5Profile.eps", |
||
| bool | autolimits = false, |
||
| int | ColorCode = 1 |
||
| ) |
Definition at line 535 of file comparisonPlots.cc.
References _outputDir, trackerHits::c, data, TrackerOfflineValidation_Dqm_cff::dirName, getHistMaxMin(), i, j, output, and alignCSCRings::s.
{
//for all histograms set the title centered
gStyle->SetTitleAlign(22);
gStyle->SetTitleX(0.5);
gStyle->SetTitleY(0.97);
//gStyle->SetShadowColor(0);
//gStyle->SetTitleSize(.5);
gStyle->SetTitleFont(62);
//int ColorCodeHalfBarrels = 1; //color seperation of moodules corresponding to the different subdets
//if you want to seperate corresponding to the halfbarrels (+z or -z) set colorcode !=1
// --------- create directory for histograms ---------
//const char* dirName = Cut;
string s;// = "profile";
s = s + dirName;
s.append("_profile");
TDirectory* plotDir = output->mkdir( s.data() );
int j = 2;
TCut zCut[6];
zCut[0] = "z < 0";
zCut[1] = "z >= 0";
if (ColorCode == 1) {
j = 6;
// zCut[6];
zCut[0] = "sublevel==1";
zCut[1] = "sublevel==2";
zCut[2] = "sublevel==3";
zCut[3] = "sublevel==4";
zCut[4] = "sublevel==5";
zCut[5] = "sublevel==6";
}
double minimumR, maximumR;
double minimumZ, maximumZ;
double minimumPhi, maximumPhi;
double minimumDR, maximumDR;
double minimumDZ, maximumDZ;
double minimumRDPhi, maximumRDPhi;
double minimumDX, maximumDX;
double minimumDY, maximumDY;
if (autolimits){
// --------- get right limits for histogram ---------
TH1F* phr = new TH1F("phr", "phr", 200, 0, 200);
TH1F* phz = new TH1F("phz", "phz", 400, -300, 300);
TH1F* phphi = new TH1F("phphi", "phphi", 200, -3.15, 3.15);
TH1F* phdr = new TH1F("phdr", "phdr", 2000, -10, 10);
TH1F* phdz = new TH1F("phdz", "phdz", 2000, -10, 10);
TH1F* phrdphi = new TH1F("phrdphi", "phrdphi", 200, -10, 10);
TH1F* phdx = new TH1F("phdx", "phy", 2000, -10, 10);
TH1F* phdy = new TH1F("phdy", "phy", 2000, -10, 10);
data->Project("phr","r",Cut);
data->Project("phz","z",Cut);
data->Project("phphi","phi",Cut);
data->Project("phdr","dr",Cut);
data->Project("phdz","dz",Cut);
data->Project("phrdphi","r*dphi",Cut );
data->Project("phdx","dx",Cut );
data->Project("phdy","dy",Cut );
getHistMaxMin(phr, maximumR, minimumR, 0);
getHistMaxMin(phz, maximumZ, minimumZ, 0);
getHistMaxMin(phphi, maximumPhi, minimumPhi, 0);
getHistMaxMin(phdr, maximumDR, minimumDR, 1);
getHistMaxMin(phdz, maximumDZ, minimumDZ, 1);
getHistMaxMin(phrdphi, maximumRDPhi, minimumRDPhi, 1);
getHistMaxMin(phdx, maximumDX, minimumDX, 1);
getHistMaxMin(phdy, maximumDY, minimumDY, 1);
}
else{
minimumR = 0., maximumR = 120.;
minimumZ = -300., maximumZ = 300.;
minimumPhi = -3.15, maximumPhi = 3.15;
minimumDR = -1, maximumDR = 1;
minimumDZ = -1, maximumDZ = 1;
minimumRDPhi = -1, maximumRDPhi = 1;
minimumDX = -1, maximumDX = 1;
minimumDY = -1, maximumDY = 1;
}
// --------- declare histograms ---------
TProfile* hprof_drVr = new TProfile("hprof_drVr","#Delta r vs. r",nBins,minimumR,maximumR,minimumDR,maximumDR);
TProfile* hprof_dzVr = new TProfile("hprof_dzVr","#Delta z vs. r",nBins,minimumR,maximumR,minimumDZ,maximumDZ);
TProfile* hprof_rdphiVr = new TProfile("hprof_rdphiVr","r#Delta #phi vs. r",nBins,minimumR,maximumR,minimumRDPhi,maximumRDPhi);
TProfile* hprof_dxVr = new TProfile("hprof_dxVr","#Delta x vs. r", nBins,minimumR,maximumR,minimumDX,maximumDX);
TProfile* hprof_dyVr = new TProfile("hprof_dyVr","#Delta y vs. r", nBins,minimumR,maximumR,minimumDY,maximumDY);
TProfile* hprof_drVz = new TProfile("hprof_drVz","#Delta r vs. z", nBins,minimumZ,maximumZ,minimumDR,maximumDR);
TProfile* hprof_dzVz = new TProfile("hprof_dzVz","#Delta z vs. z", nBins,minimumZ,maximumZ,minimumDZ,maximumDZ);
TProfile* hprof_rdphiVz = new TProfile("hprof_rdphiVz","r#Delta #phi vs. z", nBins,minimumZ,maximumZ,minimumRDPhi,maximumRDPhi);
TProfile* hprof_dxVz = new TProfile("hprof_dxVz","#Delta x vs. z", nBins,minimumZ,maximumZ,minimumDX,maximumDX);
TProfile* hprof_dyVz = new TProfile("hprof_dyVz","#Delta y vs. z", nBins,minimumZ,maximumZ,minimumDY,maximumDY);
TProfile* hprof_drVphi = new TProfile("hprof_drVphi","#Delta r vs. #phi", nBins,minimumPhi,maximumPhi,minimumDR,maximumDR);
TProfile* hprof_dzVphi = new TProfile("hprof_dzVphi","#Delta z vs. #phi", nBins,minimumPhi,maximumPhi,minimumDZ,maximumDZ);
TProfile* hprof_rdphiVphi = new TProfile("hprof_rdphiVphi","r#Delta #phi vs. #phi", nBins,minimumPhi,maximumPhi,minimumRDPhi,maximumRDPhi);
TProfile* hprof_dxVphi = new TProfile("hprof_dxVphi","#Delta x vs. #phi", nBins,minimumPhi,maximumPhi,minimumDX,maximumDX);
TProfile* hprof_dyVphi = new TProfile("hprof_dyVphi","#Delta y vs. #phi", nBins,minimumPhi,maximumPhi,minimumDY,maximumDY);
TProfile* gr_drVr_Array[j];
THStack* mgr_drVr=new THStack("mgr_drVr","");
for ( int i = 0; i < j; i++) {
data->Draw("dr:r>>hprof_drVr",Cut+zCut[i],"prof");
if (hprof_drVr->GetEntries()>0){
gr_drVr_Array[i] =(TProfile*)hprof_drVr->Clone() ;//data->Draw("dr:r",Cut+zCut[i],"prof");
gr_drVr_Array[i]->SetMarkerColor(int(i/4)+i+1);
gr_drVr_Array[i]->SetLineColor(int(i/4)+i+1);
mgr_drVr->Add(gr_drVr_Array[i]);
}
}
//y-axis scaled 2 times the max value for values >5 micron
//for max values<5 micron set fixed range to 10 micron
if ( (mgr_drVr->GetMaximum("nostack")>0.00005)||(mgr_drVr->GetMinimum("nostack")<-0.00005) ){
mgr_drVr->SetMaximum( 2.*mgr_drVr->GetMaximum("nostack") );
mgr_drVr->SetMinimum( 2.*mgr_drVr->GetMinimum("nostack") );
}
else{
mgr_drVr->SetMaximum(100*mgr_drVr->GetMaximum("nostack"));
mgr_drVr->SetMinimum( 100.*mgr_drVr->GetMinimum("nostack"));
}
//data->Project("hprof_dzVr", "dz:r",Cut,"prof");
TProfile* gr_dzVr_Array[j];
THStack* mgr_dzVr=new THStack("mgr_dzVr","");
for ( int i = 0; i < j; i++) {
data->Draw("dz:r>>hprof_dzVr",Cut+zCut[i],"prof");
if (hprof_dzVr->GetEntries()>0){
gr_dzVr_Array[i] =(TProfile*)hprof_dzVr->Clone() ;
gr_dzVr_Array[i]->SetMarkerColor(int(i/4)+i+1);
gr_dzVr_Array[i]->SetLineColor(int(i/4)+i+1);
mgr_dzVr->Add(gr_dzVr_Array[i]);
}
}
if ( (mgr_dzVr->GetMaximum("nostack")>0.00005)||(mgr_dzVr->GetMinimum("nostack")<-0.00005) ){
mgr_dzVr->SetMaximum( 2.*mgr_dzVr->GetMaximum("nostack") );
mgr_dzVr->SetMinimum( 2.*mgr_dzVr->GetMinimum("nostack") );
}
else{
mgr_dzVr->SetMaximum( 100.*mgr_dzVr->GetMaximum("nostack") );
mgr_dzVr->SetMinimum( 100.*mgr_dzVr->GetMinimum("nostack") );
}
//data->Project("hprof_rdphiVr", "r*dphi:r",Cut,"prof");
TProfile* gr_rdphiVr_Array[j];
THStack* mgr_rdphiVr=new THStack();
for ( int i = 0; i < j; i++) {
data->Draw("r*dphi:r>>hprof_rdphiVr",Cut+zCut[i],"prof");
if (hprof_rdphiVr->GetEntries()>0){
gr_rdphiVr_Array[i] = (TProfile*)hprof_rdphiVr->Clone() ;
gr_rdphiVr_Array[i]->SetMarkerColor(int(i/4)+i+1);
gr_rdphiVr_Array[i]->SetLineColor(int(i/4)+i+1);
mgr_rdphiVr->Add(gr_rdphiVr_Array[i]);
}
}
if ( (mgr_rdphiVr->GetMaximum("nostack")>0.00005)||(mgr_rdphiVr->GetMinimum("nostack")<-0.00005) ){
mgr_rdphiVr->SetMaximum( 2.*mgr_rdphiVr->GetMaximum("nostack") );
mgr_rdphiVr->SetMinimum( 2.*mgr_rdphiVr->GetMinimum("nostack") );
}
else{
mgr_rdphiVr->SetMaximum( 100.*mgr_rdphiVr->GetMaximum("nostack") );
mgr_rdphiVr->SetMinimum( 100.*mgr_rdphiVr->GetMinimum("nostack") );
}
//data->Project("hprof_dxVr", "dx:r",Cut,"prof");
TProfile* gr_dxVr_Array[j];
THStack* mgr_dxVr=new THStack();
for ( int i = 0; i < j; i++) {
data->Draw("dx:r>>hprof_dxVr",Cut+zCut[i],"prof");
if (hprof_dxVr->GetEntries()>0){
gr_dxVr_Array[i] = (TProfile*)hprof_dxVr->Clone() ;
gr_dxVr_Array[i]->SetMarkerColor(int(i/4)+i+1);
gr_dxVr_Array[i]->SetLineColor(int(i/4)+i+1);
mgr_dxVr->Add(gr_dxVr_Array[i]);
}
}
if ( (mgr_dxVr->GetMaximum("nostack")>0.00005)||(mgr_dxVr->GetMinimum("nostack")<-0.00005) ){
mgr_dxVr->SetMaximum( 2.*mgr_dxVr->GetMaximum("nostack") );
mgr_dxVr->SetMinimum( 2.*mgr_dxVr->GetMinimum("nostack") );
}
else{
mgr_dxVr->SetMaximum( 100.*mgr_dxVr->GetMaximum("nostack") );
mgr_dxVr->SetMinimum( 100.*mgr_dxVr->GetMinimum("nostack") );
}
//data->Project("hprof_dyVr", "dy:r",Cut,"prof");
TProfile* gr_dyVr_Array[j];
THStack* mgr_dyVr=new THStack();
for ( int i = 0; i < j; i++) {
data->Draw("dy:r>>hprof_dyVr",Cut+zCut[i],"prof");
if (hprof_dyVr->GetEntries()>0){
gr_dyVr_Array[i] = (TProfile*)hprof_dyVr->Clone() ;
gr_dyVr_Array[i]->SetMarkerColor(int(i/4)+i+1);
gr_dyVr_Array[i]->SetLineColor(int(i/4)+i+1);
mgr_dyVr->Add(gr_dyVr_Array[i]);
}
}
if ( (mgr_dyVr->GetMaximum("nostack")>0.00005)||(mgr_dyVr->GetMinimum("nostack")<-0.00005) ){
mgr_dyVr->SetMaximum( 2.*mgr_dyVr->GetMaximum("nostack") );
mgr_dyVr->SetMinimum( 2.*mgr_dyVr->GetMinimum("nostack") );
}
else{
mgr_dyVr->SetMaximum( 100.*mgr_dyVr->GetMaximum("nostack") );
mgr_dyVr->SetMinimum( 100.*mgr_dyVr->GetMinimum("nostack") );
}
//data->Project("hprof_drVz", "dr:z",Cut,"prof");
TProfile* gr_drVz_Array[j];
THStack* mgr_drVz=new THStack();
for ( int i = 0; i < j; i++) {
data->Draw("dr:z>>hprof_drVz",Cut+zCut[i],"prof");
if (hprof_drVz->GetEntries()>0){
gr_drVz_Array[i] = (TProfile*)hprof_drVz->Clone() ;
gr_drVz_Array[i]->SetMarkerColor(int(i/4)+i+1);
gr_drVz_Array[i]->SetLineColor(int(i/4)+i+1);
mgr_drVz->Add(gr_drVz_Array[i]);
}
}
if ( (mgr_drVr->GetMaximum("nostack")>0.00005)||(mgr_drVr->GetMinimum("nostack")<-0.00005) ){
mgr_drVz->SetMaximum( 2.*mgr_drVz->GetMaximum("nostack") );
mgr_drVz->SetMinimum( 2.*mgr_drVz->GetMinimum("nostack") );
}
else{
mgr_drVz->SetMaximum( 100.*mgr_drVz->GetMaximum("nostack") );
mgr_drVz->SetMinimum( 100.*mgr_drVz->GetMinimum("nostack") );
}
//data->Project("hprof_dzVz", "dz:z",Cut,"prof");
TProfile* gr_dzVz_Array[j];
THStack* mgr_dzVz=new THStack();
for ( int i = 0; i < j; i++) {
data->Draw("dz:z>>hprof_dzVz",Cut+zCut[i],"prof");
if (hprof_dzVz->GetEntries()>0){
gr_dzVz_Array[i] = (TProfile*)hprof_dzVz->Clone() ;
gr_dzVz_Array[i]->SetMarkerColor(int(i/4)+i+1);
gr_dzVz_Array[i]->SetLineColor(int(i/4)+i+1);
mgr_dzVz->Add(gr_dzVz_Array[i]);
}
}
if ( (mgr_dzVr->GetMaximum("nostack")>0.00005)||(mgr_dzVr->GetMinimum("nostack")<-0.00005) ){
mgr_dzVz->SetMaximum( 2.*mgr_dzVz->GetMaximum("nostack") );
mgr_dzVz->SetMinimum( 2.*mgr_dzVz->GetMinimum("nostack") );
}
else{
mgr_dzVz->SetMaximum( 100.*mgr_dzVz->GetMaximum("nostack") );
mgr_dzVz->SetMinimum( 100.*mgr_dzVz->GetMinimum("nostack") );
}
//data->Project("hprof_rdphiVz", "r*dphi:z",Cut,"prof");
TProfile* gr_rdphiVz_Array[j];
THStack* mgr_rdphiVz=new THStack();
for ( int i = 0; i < j; i++) {
data->Draw("r*dphi:z>>hprof_rdphiVz",Cut+zCut[i],"prof");
if (hprof_rdphiVz->GetEntries()>0){
gr_rdphiVz_Array[i] = (TProfile*)hprof_rdphiVz->Clone() ;
gr_rdphiVz_Array[i]->SetMarkerColor(int(i/4)+i+1);
gr_rdphiVz_Array[i]->SetLineColor(int(i/4)+i+1);
mgr_rdphiVz->Add(gr_rdphiVz_Array[i]);
}
}
if ( (mgr_rdphiVz->GetMaximum("nostack")>0.00005)||(mgr_rdphiVz->GetMinimum("nostack")<-0.00005) ){
mgr_rdphiVz->SetMaximum( 2.*mgr_rdphiVz->GetMaximum("nostack") );
mgr_rdphiVz->SetMinimum( 2.*mgr_rdphiVz->GetMinimum("nostack") );
}
else{
mgr_rdphiVz->SetMaximum( 100.*mgr_rdphiVz->GetMaximum("nostack") );
mgr_rdphiVz->SetMinimum( 100.*mgr_rdphiVz->GetMinimum("nostack") );
}
//data->Project("hprof_dxVz", "dx:z",Cut,"prof");
TProfile* gr_dxVz_Array[j];
THStack* mgr_dxVz=new THStack();
for ( int i = 0; i < j; i++) {
data->Draw("dx:z>>hprof_dxVz",Cut+zCut[i],"prof");
if (hprof_dxVz->GetEntries()>0){
gr_dxVz_Array[i] = (TProfile*)hprof_dxVz->Clone() ;
gr_dxVz_Array[i]->SetMarkerColor(int(i/4)+i+1);
gr_dxVz_Array[i]->SetLineColor(int(i/4)+i+1);
mgr_dxVz->Add(gr_dxVz_Array[i]);
}
}
if ( (mgr_dxVz->GetMaximum("nostack")>0.00005)||(mgr_dxVz->GetMinimum("nostack")<-0.00005) ){
mgr_dxVz->SetMaximum( 2.*mgr_dxVz->GetMaximum("nostack") );
mgr_dxVz->SetMinimum( 2.*mgr_dxVz->GetMinimum("nostack") );
}
else{
mgr_dxVz->SetMaximum( 100.*mgr_dxVz->GetMaximum("nostack") );
mgr_dxVz->SetMinimum( 100.*mgr_dxVz->GetMinimum("nostack") );
}
//data->Project("hprof_dyVz", "dy:z",Cut,"prof");
TProfile* gr_dyVz_Array[j];
THStack* mgr_dyVz=new THStack();
for ( int i = 0; i < j; i++) {
data->Draw("dy:z>>hprof_dyVz",Cut+zCut[i],"prof");
if (hprof_dyVz->GetEntries()>0){
gr_dyVz_Array[i] = (TProfile*)hprof_dyVz->Clone() ;
gr_dyVz_Array[i]->SetMarkerColor(int(i/4)+i+1);
gr_dyVz_Array[i]->SetLineColor(int(i/4)+i+1);
mgr_dyVz->Add(gr_dyVz_Array[i],"p");
}
}
if ( (mgr_dyVz->GetMaximum("nostack")>0.00005)||(mgr_dyVz->GetMinimum("nostack")<-0.00005) ){
mgr_dyVz->SetMaximum( 2.*mgr_dyVz->GetMaximum("nostack") );
mgr_dyVz->SetMinimum( 2.*mgr_dyVz->GetMinimum("nostack") );
}
else{
mgr_dyVz->SetMaximum( 100.*mgr_dyVz->GetMaximum("nostack") );
mgr_dyVz->SetMinimum( 100.*mgr_dyVz->GetMinimum("nostack") );
}
//data->Project("hprof_drVphi", "dr:phi",Cut,"prof");
TProfile* gr_drVphi_Array[j];
THStack* mgr_drVphi=new THStack();
for ( int i = 0; i < j; i++) {
data->Draw("dr:phi>>hprof_drVphi",Cut+zCut[i],"prof");
if (hprof_drVphi->GetEntries()>0){
gr_drVphi_Array[i] = (TProfile*)hprof_drVphi->Clone() ;
gr_drVphi_Array[i]->SetMarkerColor(int(i/4)+i+1);
gr_drVphi_Array[i]->SetLineColor(int(i/4)+i+1);
mgr_drVphi->Add(gr_drVphi_Array[i],"p");
}
}
if ( (mgr_drVphi->GetMaximum("nostack")>0.00005)||(mgr_drVphi->GetMinimum("nostack")<-0.00005) ){
mgr_drVphi->SetMaximum( 2.*mgr_drVphi->GetMaximum("nostack") );
mgr_drVphi->SetMinimum( 2.*mgr_drVphi->GetMinimum("nostack") );
}
else{
mgr_drVphi->SetMaximum( 100.*mgr_drVphi->GetMaximum("nostack") );
mgr_drVphi->SetMinimum( 100.*mgr_drVphi->GetMinimum("nostack") );
}
//data->Project("hprof_dzVphi", "dz:phi",Cut,"prof");
TProfile* gr_dzVphi_Array[j];
THStack* mgr_dzVphi=new THStack();
for ( int i = 0; i < j; i++) {
data->Draw("dz:phi>>hprof_dzVphi",Cut+zCut[i],"prof");
if (hprof_dzVphi->GetEntries()>0){
gr_dzVphi_Array[i] = (TProfile*)hprof_dzVphi->Clone() ;
gr_dzVphi_Array[i]->SetMarkerColor(int(i/4)+i+1);
gr_dzVphi_Array[i]->SetLineColor(int(i/4)+i+1);
mgr_dzVphi->Add(gr_dzVphi_Array[i],"p");
}
}
if ( (mgr_dzVphi->GetMaximum("nostack")>0.00005)||(mgr_dzVphi->GetMinimum("nostack")<-0.00005) ){
mgr_dzVphi->SetMaximum( 2.*mgr_dzVphi->GetMaximum("nostack") );
mgr_dzVphi->SetMinimum( 2.*mgr_dzVphi->GetMinimum("nostack") );
}
else{
mgr_dzVphi->SetMaximum( 100.*mgr_dzVphi->GetMaximum("nostack") );
mgr_dzVphi->SetMinimum( 100.*mgr_dzVphi->GetMinimum("nostack") );
}
//data->Project("hprof_rdphiVphi", "r*dphi:phi",Cut,"prof");
TProfile* gr_rdphiVphi_Array[j];
THStack* mgr_rdphiVphi=new THStack();
for ( int i = 0; i < j; i++) {
data->Draw("r*dphi:phi>>hprof_rdphiVphi",Cut+zCut[i],"prof");
if (hprof_rdphiVphi->GetEntries()>0){
gr_rdphiVphi_Array[i] = (TProfile*)hprof_rdphiVphi->Clone() ;
gr_rdphiVphi_Array[i]->SetMarkerColor(int(i/4)+i+1);
gr_rdphiVphi_Array[i]->SetLineColor(int(i/4)+i+1);
mgr_rdphiVphi->Add(gr_rdphiVphi_Array[i],"p");
}
}
if ( (mgr_rdphiVphi->GetMaximum("nostack")>0.00005)||(mgr_rdphiVphi->GetMinimum("nostack")<-0.00005) ){
mgr_rdphiVphi->SetMaximum( 2.*mgr_rdphiVphi->GetMaximum("nostack") );
mgr_rdphiVphi->SetMinimum( 2.*mgr_rdphiVphi->GetMinimum("nostack") );
}
else{
mgr_rdphiVphi->SetMaximum( 100.*mgr_rdphiVphi->GetMaximum("nostack") );
mgr_rdphiVphi->SetMinimum( 100.*mgr_rdphiVphi->GetMinimum("nostack") );
}
//data->Project("hprof_dxVphi", "dx:phi",Cut,"prof");
TProfile* gr_dxVphi_Array[j];
THStack* mgr_dxVphi=new THStack();
for ( int i = 0; i < j; i++) {
data->Draw("dx:phi>>hprof_dxVphi",Cut+zCut[i],"prof");
if (hprof_dxVphi->GetEntries()>0){
gr_dxVphi_Array[i] = (TProfile*)hprof_dxVphi->Clone() ;
gr_dxVphi_Array[i]->SetMarkerColor(int(i/4)+i+1);
gr_dxVphi_Array[i]->SetLineColor(int(i/4)+i+1);
mgr_dxVphi->Add(gr_dxVphi_Array[i],"p");
}
}
if ( (mgr_dxVphi->GetMaximum("nostack")>0.00005)||(mgr_dxVphi->GetMinimum("nostack")<-0.00005) ){
mgr_dxVphi->SetMaximum( 2.*mgr_dxVphi->GetMaximum("nostack") );
mgr_dxVphi->SetMinimum( 2.*mgr_dxVphi->GetMinimum("nostack") );
}
else{
mgr_dxVphi->SetMaximum( 100.*mgr_dxVphi->GetMaximum("nostack") );
mgr_dxVphi->SetMinimum( 100.*mgr_dxVphi->GetMinimum("nostack") );
}
//data->Project("hprof_dyVphi", "dy:phi",Cut,"prof");
TProfile* gr_dyVphi_Array[j];
THStack* mgr_dyVphi=new THStack();
for ( int i = 0; i < j; i++) {
data->Draw("dy:phi>>hprof_dyVphi",Cut+zCut[i],"prof");
if (hprof_dyVphi->GetEntries()>0){
gr_dyVphi_Array[i] = (TProfile*)hprof_dyVphi->Clone() ;
gr_dyVphi_Array[i]->SetMarkerColor(int(i/4)+i+1);
gr_dyVphi_Array[i]->SetLineColor(int(i/4)+i+1);
mgr_dyVphi->Add(gr_dyVphi_Array[i],"p");
}
}
if ( (mgr_dyVphi->GetMaximum("nostack")>0.00005)||(mgr_dyVphi->GetMinimum("nostack")<-0.00005) ){
mgr_dyVphi->SetMaximum( 2.*mgr_dyVphi->GetMaximum("nostack") );
mgr_dyVphi->SetMinimum( 2.*mgr_dyVphi->GetMinimum("nostack") );
}
else{
mgr_dyVphi->SetMaximum( 100.*mgr_dyVphi->GetMaximum("nostack") );
mgr_dyVphi->SetMinimum( 100.*mgr_dyVphi->GetMinimum("nostack") );
}
// --------- draw histograms ---------
/*TCanvas* c0 = new TCanvas("c0", "c0", 200, 10, 900, 300);
c0->SetFillColor(0);
c0->Divide(3,1);
c0->cd(1);
hprof_dr->Draw();
c0->cd(2);
hprof_dz->Draw();
c0->cd(3);
hprof_rdphi->Draw();
if (savePlot) c0->Print((_outputDir+"plot3x1_Profile"+plotName).c_str());
*/
mgr_drVr->SetTitle("#Delta r vs. r");
mgr_dzVr->SetTitle("#Delta z vs. r");
mgr_rdphiVr->SetTitle("r#Delta #phi vs. r");
mgr_dxVr->SetTitle("#Delta x vs. r");
mgr_dyVr->SetTitle("#Delta y vs. r");
mgr_drVz->SetTitle("#Delta r vs. z");
mgr_dzVz->SetTitle("#Delta z vs. z");
mgr_rdphiVz->SetTitle("r#Delta #phi vs. z");
mgr_dxVz->SetTitle("#Delta x vs. z");
mgr_dyVz->SetTitle("#Delta y vs. z");
mgr_drVphi->SetTitle("#Delta r vs. #phi");
mgr_dzVphi->SetTitle("#Delta z vs. #phi");
mgr_rdphiVphi->SetTitle("r#Delta #phi vs. #phi");
mgr_dxVphi->SetTitle("#Delta x vs. #phi");
mgr_dyVphi->SetTitle("#Delta y vs. #phi");
// --------- draw histograms ---------
TCanvas* c = new TCanvas("c", "c", 200, 10, 1200, 700);
c->SetFillColor(0);
data->SetMarkerSize(0.5);
data->SetMarkerStyle(6);
c->Divide(5,3);
c->cd(1);
if (!autolimits) mgr_drVr->Draw("nostack");
if (autolimits) mgr_drVr->Draw("nostack");
c->Update();
c->cd(2);
if (!autolimits) mgr_dzVr->Draw("nostack");
if (autolimits) mgr_dzVr->Draw("nostack");
c->Update();
c->cd(3);
if (!autolimits) mgr_rdphiVr->Draw("nostack");
if (autolimits) mgr_rdphiVr->Draw("nostack");
c->cd(4);
if (!autolimits) mgr_dxVr->Draw("nostack");
if (autolimits) mgr_dxVr->Draw("nostack");
c->cd(5);
if (!autolimits) mgr_dyVr->Draw("nostack");
if (autolimits) mgr_dyVr->Draw("nostack");
c->cd(6);
if (!autolimits) mgr_drVz->Draw("nostack");
if (autolimits) mgr_drVz->Draw("nostack");
c->cd(7);
if (!autolimits) mgr_dzVz->Draw("nostack");
if (autolimits) mgr_dzVz->Draw("nostack");
c->cd(8);
if (!autolimits) mgr_rdphiVz->Draw("nostack");
if (autolimits) mgr_rdphiVz->Draw("nostack");
c->cd(9);
if (!autolimits) mgr_dxVz->Draw("nostack");
if (autolimits) mgr_dxVz->Draw("nostack");
c->cd(10);
if (!autolimits) mgr_dyVz->Draw("nostack");
if (autolimits) mgr_dyVz->Draw("nostack");
c->cd(11);
if (!autolimits) mgr_drVphi->Draw("nostack");
if (autolimits) mgr_drVphi->Draw("nostack");
c->cd(12);
if (!autolimits) mgr_dzVphi->Draw("nostack");
if (autolimits) mgr_dzVphi->Draw("nostack");
c->cd(13);
if (!autolimits) mgr_rdphiVphi->Draw("nostack");
if (autolimits) mgr_rdphiVphi->Draw("nostack");
c->cd(14);
if (!autolimits) mgr_dxVphi->Draw("nostack");
if (autolimits) mgr_dxVphi->Draw("nostack");
c->cd(15);
if (!autolimits) mgr_dyVphi->Draw("nostack");
if (autolimits) mgr_dyVphi->Draw("nostack");
c->Update();
// --------- set output directory for histograms ---------
plotDir->cd();
hprof_drVr->Write(); hprof_dzVr->Write(); hprof_rdphiVr->Write(); hprof_dxVr->Write(); hprof_dyVr->Write();
hprof_drVz->Write(); hprof_dzVz->Write(); hprof_rdphiVz->Write(); hprof_dxVz->Write(); hprof_dyVz->Write();
hprof_drVphi->Write(); hprof_dzVphi->Write(); hprof_rdphiVphi->Write(); hprof_dxVphi->Write(); hprof_dyVphi->Write();
/*
// --------- project tree onto histograms ---------
data->Project("hprof_drVr", "dr:r",Cut ,"prof");
data->Project("hprof_dzVr", "dz:r",Cut ,"prof");
data->Project("hprof_rdphiVr", "r*dphi:r",Cut ,"prof");
data->Project("hprof_dxVr", "dx:r",Cut ,"prof");
data->Project("hprof_dyVr", "dy:r",Cut ,"prof");
data->Project("hprof_drVz", "dr:z",Cut ,"prof");
data->Project("hprof_dzVz", "dz:z",Cut ,"prof");
data->Project("hprof_rdphiVz", "r*dphi:z",Cut ,"prof");
data->Project("hprof_dxVz", "dx:z",Cut ,"prof");
data->Project("hprof_dyVz", "dy:z",Cut ,"prof");
data->Project("hprof_drVphi", "dr:phi",Cut ,"prof");
data->Project("hprof_dzVphi", "dz:phi",Cut ,"prof");
data->Project("hprof_rdphiVphi", "r*dphi:phi",Cut ,"prof");
data->Project("hprof_dxVphi", "dx:phi",Cut ,"prof");
data->Project("hprof_dyVphi", "dy:phi",Cut ,"prof");
// --------- draw histograms ---------
TCanvas* cp = new TCanvas("cp", "cp", 200, 10, 1200, 700);
cp->SetFillColor(0);
data->SetMarkerSize(0.5);
data->SetMarkerStyle(6);
cp->Divide(5,3);
cp->cd(1);
hprof_drVr->Draw();
cp->cd(2);
hprof_dzVr->Draw();
cp->cd(3);
hprof_rdphiVr->Draw();
cp->cd(4);
hprof_dxVr->Draw();
cp->cd(5);
hprof_dyVr->Draw();
cp->cd(6);
hprof_drVz->Draw();
cp->cd(7);
hprof_dzVz->Draw();
cp->cd(8);
hprof_rdphiVz->Draw();
cp->cd(9);
hprof_dxVz->Draw();
cp->cd(10);
hprof_dyVz->Draw();
cp->cd(11);
hprof_drVphi->Draw();
cp->cd(12);
hprof_dzVphi->Draw();
cp->cd(13);
hprof_rdphiVphi->Draw();
cp->cd(14);
hprof_dxVphi->Draw();
cp->cd(15);
hprof_dyVphi->Draw();
// --------- set output directory for histograms ---------
plotDir->cd();
hprof_drVr->Write(); hprof_dzVr->Write(); hprof_rdphiVr->Write(); hprof_dxVr->Write(); hprof_dyVr->Write();
hprof_drVz->Write(); hprof_dzVz->Write(); hprof_rdphiVz->Write(); hprof_dxVz->Write(); hprof_dyVz->Write();
hprof_drVphi->Write(); hprof_dzVphi->Write(); hprof_rdphiVphi->Write(); hprof_dxVphi->Write(); hprof_dyVphi->Write();
*/
if (savePlot) c->Print((_outputDir+"plot3x5Prof_"+plotName).c_str());
}
| void comparisonPlots::plotTwist | ( | TCut | Cut, |
| char * | dirName, | ||
| bool | savePlot = false, |
||
| std::string | plotName = "plot3x5.eps", |
||
| bool | autolimits = false, |
||
| int | ColorCode = 1 |
||
| ) |
Definition at line 60 of file comparisonPlots.cc.
References _outputDir, trackerHits::c, data, getHistMaxMin(), i, j, and output.
{
gStyle->SetTitleAlign(22);
gStyle->SetTitleX(0.5);
gStyle->SetTitleY(0.97);
gStyle->SetTitleFont(62);
TDirectory* plotDir = output->mkdir( dirName );
double minimum, maximum;
if (autolimits)
{
TH1F* phdphi = new TH1F("phdphi", "phdphi", 200, -10, 10);
data->Project("phdphi","dphi",Cut);
getHistMaxMin(phdphi, maximum, minimum, 1);
delete phdphi;
}
else
{
minimum = -0.0004, maximum = 0.0004;
}
TH2F* h_dphiVz = new TH2F("h_dphiVz","#Delta #phi vs. z", 200,minimum,maximum, 200,minimum,maximum);
h_dphiVz->SetMarkerStyle(6);
int j = 2;
TCut zCut[6];
zCut[0] = "z < 0";
zCut[1] = "z >= 0";
if (ColorCode == 1) {
j = 6;
// zCut[6];
zCut[0] = "sublevel==1";
zCut[1] = "sublevel==2";
zCut[2] = "sublevel==3";
zCut[3] = "sublevel==4";
zCut[4] = "sublevel==5";
zCut[5] = "sublevel==6";
}
std::stringstream st;
st << "dphi > " << minimum << " && dphi < " << maximum;
TCut minmaxCut(st.str().c_str());
data->Project("h_dphiVz", "dphi:z",Cut+minmaxCut);
TGraph* gr_dphiVz_Array[j];
TMultiGraph* mgr_dphiVz=new TMultiGraph();
for ( int i = 0; i < j; i++) {
data->Draw("dphi:z",Cut+zCut[i]+minmaxCut);
gr_dphiVz_Array[i] = new TGraph(data->GetSelectedRows(),data->GetV2(),data->GetV1());
gr_dphiVz_Array[i]->SetMarkerColor(int(i/4)+i+1);
gr_dphiVz_Array[i]->SetMarkerStyle(6);
gr_dphiVz_Array[i]->SetMaximum(maximum);
gr_dphiVz_Array[i]->SetMinimum(minimum);
mgr_dphiVz->Add(gr_dphiVz_Array[i],"p");
}
// --------- draw histograms ---------
TCanvas* c = new TCanvas("c_twist", "c_twist", 200, 10, 600, 600);
c->SetFillColor(0);
data->SetMarkerSize(0.5);
data->SetMarkerStyle(6);
c->Divide(1,1);
mgr_dphiVz->SetTitle("#Delta #phi vs. z");
// if (!autolimits)
// h_dphiVz->Draw();
// else
mgr_dphiVz->SetMaximum(maximum);
mgr_dphiVz->SetMinimum(minimum);
mgr_dphiVz->Draw("a");
c->Update();
plotDir->cd();
//h_dphiVz->Write();
if (savePlot) c->Print((_outputDir+"plotTwist_"+plotName).c_str());
delete c;
delete h_dphiVz;
}
| void comparisonPlots::readTree | ( | ) |
Definition at line 31 of file comparisonPlots.cc.
References alpha_, beta_, dalpha_, data, dbeta_, detDim_, dgamma_, dphi_, dr_, dx_, dy_, dz_, eta_, fin, gamma_, id_, level_, mid_, mlevel_, phi_, r_, sublevel_, useDetId_, x_, y_, and z_.
Referenced by comparisonPlots().
{
data = (TTree*)fin->Get("alignTree");
data->SetBranchAddress("id",&id_);
data->SetBranchAddress("mid",&mid_);
data->SetBranchAddress("level",&level_);
data->SetBranchAddress("mlevel",&mlevel_);
data->SetBranchAddress("sublevel",&sublevel_);
data->SetBranchAddress("x",&x_);
data->SetBranchAddress("y",&y_);
data->SetBranchAddress("z",&z_);
data->SetBranchAddress("alpha",&alpha_);
data->SetBranchAddress("beta",&beta_);
data->SetBranchAddress("gamma",&gamma_);
data->SetBranchAddress("phi",&phi_);
data->SetBranchAddress("eta",&eta_);
data->SetBranchAddress("r",&r_);
data->SetBranchAddress("dx",&dx_);
data->SetBranchAddress("dy",&dy_);
data->SetBranchAddress("dz",&dz_);
data->SetBranchAddress("dphi",&dphi_);
data->SetBranchAddress("dr",&dr_);
data->SetBranchAddress("dalpha",&dalpha_);
data->SetBranchAddress("dbeta",&dbeta_);
data->SetBranchAddress("dgamma",&dgamma_);
data->SetBranchAddress("useDetId",&useDetId_);
data->SetBranchAddress("detDim",&detDim_);
}
| void comparisonPlots::Write | ( | ) |
std::string comparisonPlots::_outputDir [private] |
Definition at line 49 of file comparisonPlots.h.
Referenced by comparisonPlots(), plot3x5(), plot3x5Profile(), and plotTwist().
float comparisonPlots::alpha_ [private] |
Definition at line 53 of file comparisonPlots.h.
Referenced by readTree().
Definition at line 36 of file comparisonPlots.h.
float comparisonPlots::beta_ [private] |
Definition at line 53 of file comparisonPlots.h.
Referenced by readTree().
float comparisonPlots::dalpha_ [private] |
Definition at line 54 of file comparisonPlots.h.
Referenced by readTree().
| TTree* comparisonPlots::data |
Definition at line 40 of file comparisonPlots.h.
Referenced by getMaxMin(), plot3x5(), plot3x5Profile(), plotTwist(), and readTree().
float comparisonPlots::dbeta_ [private] |
Definition at line 54 of file comparisonPlots.h.
Referenced by readTree().
int comparisonPlots::detDim_ [private] |
Definition at line 52 of file comparisonPlots.h.
Referenced by readTree().
float comparisonPlots::dgamma_ [private] |
Definition at line 54 of file comparisonPlots.h.
Referenced by readTree().
float comparisonPlots::dphi_ [private] |
Definition at line 54 of file comparisonPlots.h.
Referenced by getMaxMin(), and readTree().
float comparisonPlots::dr_ [private] |
Definition at line 54 of file comparisonPlots.h.
Referenced by getMaxMin(), and readTree().
float comparisonPlots::dx_ [private] |
Definition at line 54 of file comparisonPlots.h.
Referenced by getMaxMin(), and readTree().
float comparisonPlots::dy_ [private] |
Definition at line 54 of file comparisonPlots.h.
Referenced by getMaxMin(), and readTree().
float comparisonPlots::dz_ [private] |
Definition at line 54 of file comparisonPlots.h.
Referenced by getMaxMin(), and readTree().
float comparisonPlots::eta_ [private] |
Definition at line 53 of file comparisonPlots.h.
Referenced by readTree().
| TFile* comparisonPlots::fin |
Definition at line 38 of file comparisonPlots.h.
Referenced by comparisonPlots(), and readTree().
float comparisonPlots::gamma_ [private] |
Definition at line 53 of file comparisonPlots.h.
Referenced by readTree().
int comparisonPlots::id_ [private] |
Definition at line 52 of file comparisonPlots.h.
Referenced by readTree().
int comparisonPlots::level_ [private] |
Definition at line 52 of file comparisonPlots.h.
Referenced by readTree().
float comparisonPlots::maxDR [private] |
Definition at line 59 of file comparisonPlots.h.
Referenced by getMaxMin().
float comparisonPlots::maxDX [private] |
Definition at line 62 of file comparisonPlots.h.
Referenced by getMaxMin().
float comparisonPlots::maxDY [private] |
Definition at line 63 of file comparisonPlots.h.
Referenced by getMaxMin().
float comparisonPlots::maxDZ [private] |
Definition at line 60 of file comparisonPlots.h.
Referenced by getMaxMin().
float comparisonPlots::maxPhi [private] |
Definition at line 58 of file comparisonPlots.h.
Referenced by getMaxMin().
float comparisonPlots::maxR [private] |
Definition at line 56 of file comparisonPlots.h.
Referenced by getMaxMin().
float comparisonPlots::maxRDPhi [private] |
Definition at line 61 of file comparisonPlots.h.
Referenced by getMaxMin().
float comparisonPlots::maxZ [private] |
Definition at line 57 of file comparisonPlots.h.
Referenced by getMaxMin().
int comparisonPlots::mid_ [private] |
Definition at line 52 of file comparisonPlots.h.
Referenced by readTree().
float comparisonPlots::minDR [private] |
Definition at line 59 of file comparisonPlots.h.
Referenced by getMaxMin().
float comparisonPlots::minDX [private] |
Definition at line 62 of file comparisonPlots.h.
Referenced by getMaxMin().
float comparisonPlots::minDY [private] |
Definition at line 63 of file comparisonPlots.h.
Referenced by getMaxMin().
float comparisonPlots::minDZ [private] |
Definition at line 60 of file comparisonPlots.h.
Referenced by getMaxMin().
float comparisonPlots::minPhi [private] |
Definition at line 58 of file comparisonPlots.h.
Referenced by getMaxMin().
float comparisonPlots::minR [private] |
Definition at line 56 of file comparisonPlots.h.
Referenced by getMaxMin().
float comparisonPlots::minRDPhi [private] |
Definition at line 61 of file comparisonPlots.h.
Referenced by getMaxMin().
float comparisonPlots::minZ [private] |
Definition at line 57 of file comparisonPlots.h.
Referenced by getMaxMin().
int comparisonPlots::mlevel_ [private] |
Definition at line 52 of file comparisonPlots.h.
Referenced by readTree().
| TFile* comparisonPlots::output |
Definition at line 39 of file comparisonPlots.h.
Referenced by comparisonPlots(), plot3x5(), plot3x5Profile(), plotTwist(), and Write().
float comparisonPlots::phi_ [private] |
Definition at line 53 of file comparisonPlots.h.
Referenced by getMaxMin(), and readTree().
float comparisonPlots::r_ [private] |
Definition at line 53 of file comparisonPlots.h.
Referenced by getMaxMin(), and readTree().
int comparisonPlots::sublevel_ [private] |
Definition at line 52 of file comparisonPlots.h.
Referenced by readTree().
int comparisonPlots::useDetId_ [private] |
Definition at line 52 of file comparisonPlots.h.
Referenced by readTree().
float comparisonPlots::x_ [private] |
Definition at line 53 of file comparisonPlots.h.
Referenced by readTree().
float comparisonPlots::y_ [private] |
Definition at line 53 of file comparisonPlots.h.
Referenced by readTree().
float comparisonPlots::z_ [private] |
Definition at line 53 of file comparisonPlots.h.
Referenced by getMaxMin(), and readTree().
1.7.3