void TPTiming ()
{
TAxis * ax = TPMatchEmul2D->GetZaxis() ;
ax->SetRangeUser(-1,6) ;
TCanvas* canv = new TCanvas("canv", "canv") ;
canv->SetLogz(0) ;
gStyle->SetOptStat(10) ;
int color[10]= {1,10,3,4,5,6,7,8,9,2} ;
gStyle->SetPalette(7, color) ;
TPMatchEmul2D->GetXaxis()->SetTitle("Phi index");
TPMatchEmul2D->GetYaxis()->SetTitle("Eta index");
TPMatchEmul2D->Draw("colz") ;
TH2I * label = new TH2I("label", "",72, 1, 73, 38, -19, 19) ;
label->SetMarkerSize(0.6);
label->SetBit(kCanDelete);
for (int x=3 ; x<73 ; x+=4) {
for (int y=21; y<=37; y++) {
int towernb = 4*(y-21)+1 ;
label->SetBinContent(x-1, y, towernb) ; //EB+
label->SetBinContent(x, 40-y, towernb) ; //EB-
}
}
label->Draw("same text") ;
TLatex txt;
txt.SetTextSize(0.02);
TLine line;
line.SetLineColor(1) ;
line.SetLineStyle(1) ;
line.SetLineWidth(1) ;
TAxis* xAxis = TPMatchEmul2D->GetXaxis();
TAxis* yAxis = TPMatchEmul2D->GetYaxis();
// draw SM borders and numbers
float sm ;
for (int i=0; i<36 ; i++ ) {
if (i<18) {
sm = 4*i+3 ;
line.DrawLine(sm, 1, sm, 18) ;
txt.SetTextAlign(32);
txt.DrawText(sm-1+0.3, -17.7, Form("-%d",i+1));
}
else {
sm = 4*(i-18)+3 ;
line.DrawLine(sm, 0, sm, -17) ;
txt.SetTextAlign(12);
txt.DrawText(sm-2+0.3, 18.5, Form("+%d",i-17));
}
}
line.DrawLine(1, 0, 73, 0) ;
line.DrawLine(1, -17, 73, -17) ;
line.DrawLine(1, 1, 73, 1) ;
line.DrawLine(1, 18, 73, 18) ;
}
void CCDimage14(int firstrun, int lastrun)
{
gStyle->SetPalette(1,0);
gSystem->Load("libMaxCam");
gSystem->Load("libWaveformTools.so");
std::stringstream sstm;
TH2F * sum = 0;
TH2I * count = 0;
TH2F * p=0;
sum = new TH2F("alpha_sum","alpha_sum",256,0,1024,256,0,1024);
count = new TH2I("alpha_count","alpha_count",256,0,1024,256,0,1024);
const int c = 0;
gROOT->cd();
int runnum;
int setnum;
if (firstrun>=889 && lastrun<=954) {setnum=0;}
else if (firstrun>=955 && lastrun<=3868) {setnum=1;}
else if (firstrun>=3876 && lastrun<=5789) {setnum=2;}
else if (firstrun>=5792) {setnum=3;}
else
{
gApplication->Terminate();
exit;}
for (int x = firstrun; x <= lastrun; x++){
runnum=x;
string origfile = "/net/nudsk0001/d00/scratch/dctpc_tmp/bigdctpc_data/BigDCTPC_run_";
string skimfile = "/net/nudsk0001/d00/scratch/dctpc_tmp/bigdctpc_skim/BigDCTPC_run_";
// string origfile = "/net/hisrv0001/home/spitzj/myOutput_";
// string skimfile = "/net/hisrv0001/home/spitzj//myOutput_";
string skimend = "skim.root";
string origend = ".root";
string origfilename;
string skimfilename;
sstm.str("");
if (x<10000){ origfile+="0"; skimfile+="0"; }
if (x<1000){ origfile+="0"; skimfile+="0"; }
if (x<100){ origfile+="0"; skimfile+="0"; }
if (x<10){ origfile+="0"; skimfile+="0"; }
sstm << origfile << x << origend;
origfilename = sstm.str();
sstm.str("");
sstm << skimfile << x << skimend;
skimfilename = sstm.str();
cout << origfilename << endl;
ifstream ifile(origfilename.c_str());
ifstream ifile2(skimfilename.c_str());
if(!ifile)
continue;
if(!ifile2)
continue;
DmtpcSkimDataset d;
d.openRootFile(skimfilename.c_str());
d.loadDmtpcEvent(true,origfilename.c_str());
for (int i = 0; i < d.nevents(); i ++){
if(i%100==0){
cout<<"Event: "<<i<<endl;
}
//cout<<"here"<<endl
d.getEvent(i);
//cout<<"here2"<<endl;
for (int t =0 ; t < d.event()->ntracks(c); t++)
{
if(d.event()->maxpixel(0,t)>150){continue;}
if(d.event()->spark(0)){continue;}
p = (TH2F*)d.event()->cluster(c)->getImage();
vector<int> clust = d.event()->cluster(c)->getCluster(t);
for (vector<int>::iterator it = clust.begin(); it!=clust.end(); it++)
{
sum->SetBinContent(*it, sum->GetArray()[*it] + p->GetArray()[*it]);
count->SetBinContent(*it, count->GetArray()[*it] + 1);
}
}
}
}
TH2F * normalized = (TH2F*) sum->Clone("normalized");
normalized->SetName("normalized");
normalized->Divide(count);
normalized->SetStats(0);
normalized->Draw("COLZ");
c1->Print(Form("CCDimage14_set_%d_runs_%d_%d.pdf",setnum,firstrun,lastrun));
TFile *f=new TFile(Form("CCDimage14_set_%d_runs_%d_%d.root",setnum,firstrun,lastrun),"RECREATE");
normalized->Write();
f->Close();
gApplication->Terminate();
}
void* Monitoring::MonitorThread(void* arg){
//std::cout<<"d1"<<std::endl;
monitor_thread_args* args= static_cast<monitor_thread_args*>(arg);
std::string outpath=args->outputpath;
zmq::socket_t Ireceive (*(args->context), ZMQ_PAIR);
Ireceive.connect("inproc://MonitorThread");
// std::vector<CardData*> carddata;
std::map<int,std::vector<TH1F> > PedTime;
std::map<int,std::vector<TH1F> > PedRMSTime;
std::vector<TH1F> rates;
std::vector<TH1F> averagesize;
std::vector<TH1I> tfreqplots;
std::map<int,std::vector<std::vector<float > > > pedpars;
TCanvas c1("c1","c1",600,400);
bool running=true;
bool init=true;
std::vector<PMT> PMTInfo;
/////////////////// Connect to sql ///////////////////////
//std::cout<<"d2"<<std::endl;
pqxx::connection *C;
std::stringstream tmp;
tmp<<"dbname=annie"<<" hostaddr=127.0.0.1"<<" port=5432" ;
C=new pqxx::connection(tmp.str().c_str());
if (C->is_open()) {
// std::cout << "Opened database successfully: " << C->dbname() << std::endl;
}
else {
std::cout << "Can't open database" << std::endl;
return false;
}
tmp.str("");
pqxx::nontransaction N(*C);
tmp<<"select gx,gy,gz,vmecard,vmechannel from pmtconnections order by channel; ";
/* Execute SQL query */
pqxx::result R( N.exec( tmp.str().c_str() ));
//pqxx::result::const_iterator c = R.begin();
///////// Fill PMT Info////////////////
for ( pqxx::result::const_iterator c = R.begin(); c != R.end(); ++c) {
PMT tmp;
tmp.gx= c[0].as<int>();
tmp.gy= c[1].as<int>();
tmp.gz= c[2].as<int>();
tmp.card= c[3].as<int>();
tmp.channel= c[4].as<int>()-1;
PMTInfo.push_back(tmp);
}
//std::cout<<"d3"<<std::endl;
while (running){
//std::cout<<"d4"<<std::endl;
zmq::message_t comm;
Ireceive.recv(&comm);
std::istringstream iss(static_cast<char*>(comm.data()));
std::string arg1="";
iss>>arg1;
//std::cout<<"d5"<<std::endl;
if(arg1=="Data"){
////////// Setting up plots/////////
std::vector<TGraph2D*> mg;
TH2I EventDisplay ("Event Display", "Event Display", 10, -1, 8, 10, -1, 8);
TH2I RMSDisplay ("RMS Display", "RMS Display", 10, -1, 8, 10, -1, 8);
std::vector<TH1F> temporalplots;
std::vector<TH1I> freqplots;
CardData* carddata;
int size=0;
iss>>size;
//freqplots.clear();
//std::cout<<"d6"<<std::endl;
for(int i=0;i<size;i++){
//std::cout<<"d7"<<std::endl;
long long unsigned int pointer;
iss>>std::hex>>pointer;
carddata=(reinterpret_cast<CardData *>(pointer));
if(init){ ////make initial freq plot and ped vector ped time and ped rms////
for(int j=0;j<carddata->channels;j++){
std::stringstream tmp;
tmp<<"Channel "<<(i*4)+j<<" frequency";
TH1I tmpfreq(tmp.str().c_str(),tmp.str().c_str(),10,0,9);
tfreqplots.push_back(tmpfreq);
tmp.str("");
tmp<<"Channel "<<(i*4)+j<<" Pedistal";
TH1F tmppedtime(tmp.str().c_str(),tmp.str().c_str(),100,0,99);
PedTime[carddata->CardID].push_back(tmppedtime);
tmp.str("");
tmp<<"Channel "<<(i*4)+j<<" Pedistal RMS";
TH1F tmppedrmstime(tmp.str().c_str(),tmp.str().c_str(),100,0,99);
PedRMSTime[carddata->CardID].push_back(tmppedrmstime);
std::vector<float> tmppedpars;
tmppedpars.push_back(0);
tmppedpars.push_back(0);
pedpars[carddata->CardID].push_back(tmppedpars);
}
if(i==size-1)init=false;
}
//std::cout<<"d8"<<std::endl;
// std::cout<<"d1"<<std::endl;
///////Make temporal plot //////////
for(int j=0;j<carddata->channels;j++){
std::stringstream tmp;
tmp<<"Channel "<<(i*4)+j<<" temporal";
TH1F temporal(tmp.str().c_str(),tmp.str().c_str(),carddata->buffersize,0,carddata->buffersize-1);
long sum=0;
//////////Make freq plot///////////////
tmp.str("");
tmp<<"Channel "<<(i*4)+j<<" frequency";
TH1I freq(tmp.str().c_str(),tmp.str().c_str(),200,200,399);
// std::cout<<"d2"<<std::endl;
//std::cout<<"d9"<<std::endl;
///// Calculate sum for event dispkay and fill freq plots /////////
for(int k=0;k<carddata->buffersize;k++){
//std::cout<<"d10"<<std::endl;
// std::cout<<"i="<<i<<" j="<<j<<std::endl;
//std::cout<<"d2.5 "<<(i*4)+j<<" feqplot.size = "<<freqplots.size()<<std::endl;
if(carddata->Data[(j*carddata->buffersize)+k]>pedpars[carddata->CardID].at(j).at(0)+(pedpars[carddata->CardID].at(j).at(1)*5))sum+=carddata->Data[(j*carddata->buffersize)+k];
freq.Fill(carddata->Data[(j*carddata->buffersize)+k]);
//temporal.SetBinContent(k,carddata->Data[(j*carddata->buffersize)+k]);
}
freqplots.push_back(freq);
//////// find pedistall fill ped temporals//////////
freq.Fit("gaus");
TF1 *gaus = freq.GetFunction("gaus");
pedpars[carddata->CardID].at(j).at(0)=(gaus->GetParameter(1));
pedpars[carddata->CardID].at(j).at(1)=(gaus->GetParameter(2));
gaus->SetLineColor(j+1);
//std::cout<<"d11"<<std::endl;
for(int bin=99;bin>0;bin--){
PedTime[carddata->CardID].at(j).SetBinContent(bin,PedTime[carddata->CardID].at(j).GetBinContent(bin-1));
PedRMSTime[carddata->CardID].at(j).SetBinContent(bin,PedRMSTime[carddata->CardID].at(j).GetBinContent(bin-1));
}
PedTime[carddata->CardID].at(j).SetBinContent(0, pedpars[carddata->CardID].at(j).at(0));
PedRMSTime[carddata->CardID].at(j).SetBinContent(0, pedpars[carddata->CardID].at(j).at(1));
//////// fill temporal plot/////////
for(int k=0;k<carddata->buffersize/4;k++){
//std::cout<<"d12"<<std::endl;
//std::cout<<"j*4 = "<<j*4<<std::endl;
//std::cout<<"(i*BufferSize)+(j*4) = "<<(i*BufferSize)+(j*4)<<std::endl;
//std::cout<<"i*BufferSize)+(j*4)+(BufferSize/2) = "<<(i*BufferSize)+(j*4)+(BufferSize/2)<<std::endl;
//std::cout<<"(i*BufferSize)+(j*4)+(BufferSize/2)+1 = "<<(i*BufferSize)+(j*4)+(BufferSize/2)+1<<std::endl;
int offset=pedpars[carddata->CardID].at(j).at(0);
double conversion=2.415/pow(2.0, 12.0);
temporal.SetBinContent(k*4,(carddata->Data[(j*carddata->buffersize)+(k*2)]-offset)*conversion);
temporal.SetBinContent((k*4)+1,(carddata->Data[(j*carddata->buffersize)+(k*2)+1]-offset)*conversion);
temporal.SetBinContent((k*4)+2,(carddata->Data[(j*carddata->buffersize)+(k*2)+(carddata->buffersize/2)]-offset)*conversion);
temporal.SetBinContent((k*4)+3,(carddata->Data[(j*carddata->buffersize)+(k*2)+(carddata->buffersize/2)+1]-offset)*conversion);
}
//std::cout<<"d13"<<std::endl;
//std::cout<<"d3"<<std::endl;
temporalplots.push_back(temporal);
////// find x,y,z fill event display /////////
int x=-10;
int z=-10;
int y=-10;
for(int pmt=0;pmt<PMTInfo.size();pmt++){
//std::cout<<"d4"<<std::endl;
if(PMTInfo.at(pmt).card==carddata->CardID && PMTInfo.at(pmt).channel==j){
x=PMTInfo.at(pmt).gx;
z=PMTInfo.at(pmt).gz;
y=PMTInfo.at(pmt).gy;
//std::cout<<"d15"<<std::endl;
}
}
/*
int x=(((i*4)+j)%8);
int y=(floor(((i*4)+j)/8.0));
if(x==0 && y==0){x=-10;y=-10;}
if(x==7 && y==0){x=-10;y=-10;}
if(x==0 && y==7){x=-10;y=-10;}
if(x==7 && y==7){x=-10;y=-10;}
if (y>7){x=-10;y=-10;}
std::cout<<"i="<<i<<" j="<<j<<" (i*4)+j)="<<((i*4)+j)<<" x="<<x<<" y="<<y<<" sum="<<sum<<std::endl;
EventDisplay.SetBinContent(x+1,y+1,sum);
//EventDisplay.SetBinContent(((i*4)+j),sum);
*/
//std::cout<<"d16"<<std::endl;
if(x!=-10 && z!=-10){
//std::cout<<"d17"<<std::endl;
//std::cout<<"gx = "<<x<<" , gz="<<z<<std::endl;
EventDisplay.SetBinContent(x+2,z+2,sum);
RMSDisplay.SetBinContent(x+2,z+2,gaus->GetParameter(2)*100);
//// Attempted 2ne event display ///
TGraph2D *dt=new TGraph2D(1);
dt->SetPoint(0,x,z,z);
dt->SetMarkerStyle(20);
//dt->GetXaxis()->SetRangeUser(-1,8);
// dt->GetYaxis()->SetRangeUser(-1,8);
//dt->GetZaxis()->SetRangeUser(-1,8);
mg.push_back(dt);
}
}
//std::cout<<"d18"<<std::endl;
///////Find max freq for scaling//////////
int maxplot=0;
long maxvalue=0;
//std::cout<<"i="<<i<<" (i*4)="<<(i*4)<<" (i*4)+4="<<(i*4)+4<<" size="<<freqplots.size()<<std::endl;
for(int j=(i*4);j<(i*4)+4;j++){
if (freqplots.at(j).GetMaximum()>maxvalue){
//std::cout<<"d19"<<std::endl;
maxvalue=freqplots.at(j).GetMaximum();
maxplot=j;
}
}
//std::cout<<"d20"<<std::endl;
////////Find current time and plot frewuency plot
time_t t = time(0); // get time now
struct tm * now = localtime( & t );
std::stringstream title;
title<<"Card "<<carddata->CardID<<" frequency: "<<(now->tm_year + 1900) << '-' << (now->tm_mon + 1) << '-' << now->tm_mday<<','<<now->tm_hour<<':'<<now->tm_min<<':'<<now->tm_sec;
//std::cout<<"d21"<<std::endl;
freqplots.at(maxplot).SetTitle(title.str().c_str());
freqplots.at(maxplot).GetXaxis()->SetTitle("ADC Value");
freqplots.at(maxplot).GetYaxis()->SetTitle("Frequency");
freqplots.at(maxplot).SetLineColor((maxplot%4)+1);
freqplots.at(maxplot).Draw();
TLegend leg(0.8,0.4,1.0,0.7);
//leg.SetHeader("The Legend Title");
//std::cout<<"d22"<<std::endl;
for(int j=(i*4);j<(i*4)+4;j++){
//std::cout<<"d23"<<std::endl;
std::stringstream legend;
legend<<"Channel "<<j-(i*4);
leg.AddEntry(&freqplots.at(j),legend.str().c_str(),"l");
freqplots.at(j).SetLineColor((j%4)+1);
if(j==maxplot){;}//freqplots.at(i).Draw();
else freqplots.at(j).Draw("same");
}
leg.Draw();
//std::cout<<"d24"<<std::endl;
std::stringstream tmp;
tmp<<outpath<<carddata->CardID<<"freq.jpg";
c1.SaveAs(tmp.str().c_str());
//std::cout<<"d25"<<std::endl;
///////find max tmporal plot for scaling
///temporal
maxplot=0;
maxvalue=0;
//std::cout<<"i="<<i<<" (i*4)="<<(i*4)<<" (i*4)+4="<<(i*4)+4<<" size="<<freqplots.size()<<std::endl;
for(int j=(i*4);j<(i*4)+4;j++){
//std::cout<<"d26"<<std::endl;
if (temporalplots.at(j).GetMaximum()>maxvalue){
maxvalue=temporalplots.at(j).GetMaximum();
maxplot=j;
}
}
//std::cout<<"d27"<<std::endl;
//////// Find time and plot temporal plot
t = time(0); // get time now
now = localtime( & t );
std::stringstream title2;
title2<<"Card "<<carddata->CardID<<" Temporal: "<<(now->tm_year + 1900) << '-' << (now->tm_mon + 1) << '-' << now->tm_mday<<','<<now->tm_hour<<':'<<now->tm_min<<':'<<now->tm_sec;
//std::cout<<"d28"<<std::endl;
temporalplots.at(maxplot).SetTitle(title2.str().c_str());
temporalplots.at(maxplot).GetXaxis()->SetTitle("Samples");
temporalplots.at(maxplot).GetYaxis()->SetTitle("Volate (V)");
temporalplots.at(maxplot).SetLineColor((maxplot%4)+1);
temporalplots.at(maxplot).Draw();
TLegend leg2(0.8,0.4,1.0,0.7);
//std::cout<<"d29"<<std::endl;
for(int j=(i*4);j<(i*4)+4;j++){
//std::cout<<"d30"<<std::endl;
std::stringstream legend;
legend<<"Channel "<<j-(i*4);
leg2.AddEntry(&temporalplots.at(j),legend.str().c_str(),"l");
temporalplots.at(j).SetLineColor((j%4)+1);
if(j==0){;}
else temporalplots.at(j).Draw("same");
}
//std::cout<<"d31"<<std::endl;
leg2.Draw();
//std::cout<<"d6"<<std::endl;
std::stringstream tmp2;
tmp2<<outpath<<carddata->CardID<<"temporal.jpg";
c1.SaveAs(tmp2.str().c_str());
//std::cout<<"d32"<<std::endl;
///////plotting PED time and ped rms time //////
maxplot=0;
maxvalue=0;
for(int j=0;j<4;j++){
//std::cout<<"d26"<<std::endl;
if ( PedTime[carddata->CardID].at(j).GetMaximum()>maxvalue){
maxvalue=PedTime[carddata->CardID].at(j).GetMaximum();
maxplot=j;
}
}
t = time(0); // get time now
now = localtime( & t );
title2.str("");
title2<<"Card "<<carddata->CardID<<" Pedistal Variation: "<<(now->tm_year + 1900) << '-' << (now->tm_mon + 1) << '-' << now->tm_mday<<','<<now->tm_hour<<':'<<now->tm_min<<':'<<now->tm_sec;
//std::cout<<"d28"<<std::endl;
PedTime[carddata->CardID].at(maxplot).SetTitle(title2.str().c_str());
PedTime[carddata->CardID].at(maxplot).GetXaxis()->SetTitle("Samples");
PedTime[carddata->CardID].at(maxplot).GetYaxis()->SetTitle("ADC Value");
PedTime[carddata->CardID].at(maxplot).SetLineColor((maxplot%4)+1);
PedTime[carddata->CardID].at(maxplot).Draw();
TLegend leg3(0.8,0.4,1.0,0.7);
//std::cout<<"d29"<<std::endl;
for(int j=0;j<4;j++){
//std::cout<<"d30"<<std::endl;
std::stringstream legend;
legend<<"Channel "<<j;
leg3.AddEntry(&PedTime[carddata->CardID].at(j),legend.str().c_str(),"l");
PedTime[carddata->CardID].at(j).SetLineColor((j%4)+1);
PedTime[carddata->CardID].at(j).Draw("same");
}
//std::cout<<"d31"<<std::endl;
leg3.Draw();
//std::cout<<"d6"<<std::endl;
tmp2.str("");
tmp2<<outpath<<"plots2/"<<carddata->CardID<<"PedTime.jpg";
c1.SaveAs(tmp2.str().c_str());
PedRMSTime[carddata->CardID].at(0).Draw();
for (int channel=1;channel<4;channel++){
PedRMSTime[carddata->CardID].at(channel).Draw("same");
}
tmp2.str("");
tmp2<<outpath<<"plots2/"<<carddata->CardID<<"PedRMSTime.jpg";
c1.SaveAs(tmp2.str().c_str());
delete carddata;
} /// size i
/*
//std::cout<<"d4"<<std::endl;
int maxplot=0;
long maxvalue=0;
for(int i=0;i<freqplots.size();i++){
if (freqplots.at(i).GetMaximum()>maxvalue){
maxvalue=freqplots.at(i).GetMaximum();
maxplot=i;
}
}
freqplots.at(maxplot).SetLineColor(maxplot+1);
freqplots.at(maxplot).Draw();
for(int i=0;i<freqplots.size();i++){
// Double_t scale = 1/freqplots.at(i).GetMaximum();
// freqplots.at(i).Scale(scale);
freqplots.at(i).SetLineColor(i+1);
if(i==maxplot);//freqplots.at(i).Draw();
else freqplots.at(i).Draw("same");
//freqplots.at(i).Scale((1.0/scale));
}
std::cout<<"d5"<<std::endl;
std::stringstream tmp;
tmp<<outpath<<"freq.jpg";
c1.SaveAs(tmp.str().c_str());
for(int i=0;i<temporalplots.size();i++){
temporalplots.at(i).SetLineColor(i+1);
if(i==0)temporalplots.at(i).Draw();
else temporalplots.at(i).Draw("same");
}
std::cout<<"d6"<<std::endl;
std::stringstream tmp2;
tmp2<<outpath<<"temporal.jpg";
c1.SaveAs(tmp2.str().c_str());
*/
//std::cout<<"d33"<<std::endl;
/////////plot event display /////////
EventDisplay.Draw("COLZ");
std::stringstream tmp3;
tmp3<<outpath<<"0EventDisplay.jpg";
c1.SaveAs(tmp3.str().c_str());
//std::cout<<"d34 ="<<mg.size()<<std::endl;
/////////plot RMS display /////////
RMSDisplay.Draw("COLZ");
tmp3.str("");
tmp3<<outpath<<"0RMSDisplay.jpg";
c1.SaveAs(tmp3.str().c_str());
///plot atempted 3d event display///////
tmp3.str("");
if(mg.size()>0) mg.at(0)->Draw();
for(int plots=1;plots<mg.size();plots++){
mg.at(plots)->Draw("same");
//std::cout<<"d35"<<std::endl;
}
tmp3<<outpath<<"0EventDisplay3D.jpg";
//c1.SaveAs(tmp3.str().c_str());
}
else if(arg1=="Quit"){