| 33 |
using namespace std; |
using namespace std; |
| 34 |
|
|
| 35 |
Q2TH::Q2TH(TString host, TString user, TString psw){ |
Q2TH::Q2TH(TString host, TString user, TString psw){ |
| 36 |
|
this->Open(host,user,psw); |
| 37 |
|
}; |
| 38 |
|
|
| 39 |
|
void Q2TH::Open(TString host, TString user, TString psw){ |
| 40 |
fh = gSystem->ExpandPathName(host.Data()); |
fh = gSystem->ExpandPathName(host.Data()); |
| 41 |
fu = gSystem->ExpandPathName(user.Data()); |
fu = gSystem->ExpandPathName(user.Data()); |
| 42 |
fp = gSystem->ExpandPathName(psw.Data()); |
fp = gSystem->ExpandPathName(psw.Data()); |
| 43 |
|
printf(" Connecting to DB %s \n",fh.Data()); |
| 44 |
dbc = TSQLServer::Connect(fh.Data(),fu.Data(),fp.Data()); |
dbc = TSQLServer::Connect(fh.Data(),fu.Data(),fp.Data()); |
| 45 |
|
if ( dbc && dbc->IsConnected() ){ |
| 46 |
|
printf(" connected! \n"); |
| 47 |
|
} else { |
| 48 |
|
printf(" ERROR! not connected... :( \n"); |
| 49 |
|
}; |
| 50 |
}; |
}; |
| 51 |
|
|
| 52 |
TObject *Q2TH::Draw(TString query, Bool_t verbose, TString hname){ |
TObject *Q2TH::Draw(TString query, Bool_t verbose, TString hname){ |
| 53 |
// |
// |
| 54 |
|
if ( !strcmp(query.Data(),"help") ){ |
| 55 |
|
printf(" USAGE: \n"); |
| 56 |
|
printf(" 1) start root and create Q2TH object with \n"); |
| 57 |
|
printf(" Q2TH *qt = new Q2TH() \n"); |
| 58 |
|
printf(" or \n"); |
| 59 |
|
printf(" Q2TH *qt = new Q2TH(\"mysql://srvg-g2-01.ts.infn.it/pamelaProcessing9_TS\",\"pamelaprod_ro\",\"mypassword\") \n"); |
| 60 |
|
printf(" 2) query the DB with \n"); |
| 61 |
|
printf(" qt->Draw(\"select REAL_TIME_INIT,BAD_PKT_PERCENTAGE from ROOT_TABLE_MERGING;\"); \n"); |
| 62 |
|
printf(" qt->Draw(\"select REAL_TIME_INIT,BAD_PKT_PERCENTAGE from ROOT_TABLE_MERGING;\",true); this will print numbers on screen \n"); |
| 63 |
|
printf(" qt->Draw(\"select REAL_TIME_INIT from ROOT_TABLE_MERGING;\",true,\"myhisto\"); this will print numbers on screen and create histo \"myhisto\"\n"); |
| 64 |
|
printf(" 3) to use your own THxD create it and then query the DB giving as argument the name of histo: \n"); |
| 65 |
|
printf(" TH2D *myhisto=new TH2D(\"myhisto\",\"myhisto\",5000,1140000000.,1240000000.,10000,0.,1.) \n"); |
| 66 |
|
printf(" qt->Draw(\"select REAL_TIME_INIT,BAD_PKT_PERCENTAGE from ROOT_TABLE_MERGING;\",false,\"myhisto\")\n\n\n"); |
| 67 |
|
|
| 68 |
|
return NULL; |
| 69 |
|
}; |
| 70 |
|
// |
| 71 |
pResult = dbc->Query(query.Data()); |
pResult = dbc->Query(query.Data()); |
| 72 |
// |
// |
| 73 |
Row = pResult->Next(); |
Row = pResult->Next(); |
| 101 |
}; |
}; |
| 102 |
pResult->Delete(); |
pResult->Delete(); |
| 103 |
// |
// |
| 104 |
|
|
| 105 |
|
// |
| 106 |
Int_t f1bin = 70; |
Int_t f1bin = 70; |
| 107 |
Int_t f2bin = 70; |
Int_t f2bin = 70; |
| 108 |
if ( dim == 1 ){ |
if ( dim == 1 ){ |
| 109 |
f1bin = int((maxf1-minf1)/1000.); |
f1bin = int((maxf1-minf1)/1000.); |
| 110 |
if ( f1bin < 70 ) f1bin = 70; |
if ( f1bin < 70 ) f1bin = 70; |
| 111 |
if ( f1bin > 1000 ) f1bin = 1000; |
if ( f1bin > 1000 ) f1bin = 1000; |
| 112 |
h1 = new TH1D(Form("%s1",hname.Data()),Form("%s1",hname.Data()),f1bin,minf1*0.98,maxf1*1.02); |
if ( !strcmp(hname.Data(),"q2th") ) hname += "1"; |
| 113 |
|
// h1 = dynamic_cast<TH1D*>(gDirectory->FindObject(hname.Data())); |
| 114 |
|
h1 = (TH1D*)(gDirectory->FindObject(hname.Data())); |
| 115 |
|
if ( !strcmp(hname.Data(),"q2th1") ){ |
| 116 |
|
if ( h1 ) h1->Delete(); |
| 117 |
|
}; |
| 118 |
|
if ( !h1 ) h1 = new TH1D(hname.Data(),hname.Data(),f1bin,minf1*0.98,maxf1*1.02); |
| 119 |
// h1->SetBit(TH1::kCanRebin); |
// h1->SetBit(TH1::kCanRebin); |
| 120 |
if ( verbose ) printf("\n\n Row %s \n",pResult->GetFieldName(0)); |
if ( verbose ) printf("\n\n Row %s \n",pResult->GetFieldName(0)); |
| 121 |
}; |
}; |
| 123 |
f2bin = int((maxf2-minf2)/1000.); |
f2bin = int((maxf2-minf2)/1000.); |
| 124 |
if ( f2bin < 70 ) f2bin = 70; |
if ( f2bin < 70 ) f2bin = 70; |
| 125 |
if ( f2bin > 1000 ) f2bin = 1000; |
if ( f2bin > 1000 ) f2bin = 1000; |
| 126 |
h2 = new TH2D(Form("%s2",hname.Data()),Form("%s2",hname.Data()),f1bin,minf1*0.98,maxf1*1.02,f2bin,minf2*0.98,maxf2*1.02); |
if ( !strcmp(hname.Data(),"q2th") ) hname += "2"; |
| 127 |
|
// h2 = dynamic_cast<TH2D*>(gDirectory->FindObject(hname.Data())); |
| 128 |
|
h2 = (TH2D*)(gDirectory->FindObject(hname.Data())); |
| 129 |
|
if ( !strcmp(hname.Data(),"q2th2") ){ |
| 130 |
|
if ( h2 ) h2->Delete(); |
| 131 |
|
}; |
| 132 |
|
if ( !h2 ) h2 = new TH2D(hname.Data(),hname.Data(),f1bin,minf1*0.98,maxf1*1.02,f2bin,minf2*0.98,maxf2*1.02); |
| 133 |
// h2->SetBit(TH2::kCanRebin); |
// h2->SetBit(TH2::kCanRebin); |
| 134 |
if ( verbose ) printf("\n\n Row %s %s \n",pResult->GetFieldName(0),pResult->GetFieldName(1)); |
if ( verbose ) printf("\n\n Row %s %s \n",pResult->GetFieldName(0),pResult->GetFieldName(1)); |
| 135 |
}; |
}; |
| 143 |
while ( Row ){ |
while ( Row ){ |
| 144 |
f1 = (Double_t)atof(Row->GetField(0)); |
f1 = (Double_t)atof(Row->GetField(0)); |
| 145 |
if ( dim == 1 ){ |
if ( dim == 1 ){ |
| 146 |
if ( verbose ) printf(" %f \n",r,f1); |
if ( verbose ) printf(" %i %f \n",r,f1); |
| 147 |
h1->Fill(f1); |
h1->Fill(f1); |
| 148 |
} else { |
} else { |
| 149 |
f2 = (Double_t)atof(Row->GetField(1)); |
f2 = (Double_t)atof(Row->GetField(1)); |
| 150 |
if ( verbose ) printf(" %f %f \n",r,f1,f2); |
if ( verbose ) printf(" %i %f %f \n",r,f1,f2); |
| 151 |
h2->Fill(f1,f2); |
h2->Fill(f1,f2); |
| 152 |
}; |
}; |
| 153 |
r++; |
r++; |
| 154 |
Row = pResult->Next(); |
Row = pResult->Next(); |
| 155 |
}; |
}; |
| 156 |
// |
// |
| 157 |
TCanvas *c = new TCanvas(Form("%sc",hname.Data())); |
TCanvas *c = NULL; |
| 158 |
|
TString cname = Form("%sc",hname.Data()); |
| 159 |
|
// c = dynamic_cast<TCanvas*>(gDirectory->FindObject(hname.Data())); |
| 160 |
|
c = (TCanvas*)(gDirectory->FindObject(cname.Data())); |
| 161 |
|
if ( !c ) c = new TCanvas(Form("%sc",cname.Data())); |
| 162 |
|
c->Clear(); |
| 163 |
c->cd(); |
c->cd(); |
| 164 |
if ( dim == 1 ) h1->Draw(); |
if ( dim == 1 ) h1->Draw(); |
| 165 |
if ( dim == 2 ) h2->Draw(); |
if ( dim == 2 ) h2->Draw(); |