void plotallchchipx(int iEntry = 0, int chip = 0,const char *title = "",int runnb = 0) { char name[100]; Char_t iChan; // Channel number Char_t iChip; // Channel number UShort_t iClLength; //number of samples UInt_t iEventnb; Int_t iSamples[1024]; // Samples Int_t iMax = 0; TH1F* hSignal[16]; int vschip[16]; int vschan[16]; TCanvas *c1 = new TCanvas("c1","Single cluster",600,600); // TChain ch("tree"); sprintf(name, "readout-%i_0.root", runnb); TFile *inFile = new TFile(name); TTree *ch = (TTree*) inFile->Get("tree"); if(!ch) { cerr << "Tree SALTRO events was not found: " << name << "\n"; return; } printf("Run: %s",name); // Assign variables to branches ch->SetBranchAddress("iEventnb",&iEventnb); ch->SetBranchAddress("iSamples",&iSamples); ch->SetBranchAddress("iChan",&iChan); ch->SetBranchAddress("iChip",&iChip); ch->SetBranchAddress("iClLength",&iClLength); // c1->Divide(1,16); Int_t nentries = ch->GetEntries(); printf("ENTRIES %d\n",nentries); for (Int_t j = 0; j < 16; j++) { vschip[j] = -1; vschan[j] = -1; sprintf(name, "hSignal_%d", j); hSignal[j] = ((TH1F *)(gROOT->FindObject(name))); if (hSignal[j]) delete hSignal[j]; hSignal[j] = new TH1F(name, "",1024,-0.5,1023.5); hSignal[j]->SetMinimum(0); // hSignal[j]->SetMaximum(1024); } gStyle->SetOptStat(0); int iRun = 0; for (Int_t j = 0; j < nentries; j++) { ch->GetEntry(j); if (iEventnb == iEntry) { for(Int_t l = 0; l < iClLength; l++) { hSignal[iRun]->SetBinContent(l, iSamples[iClLength-l]); } vschip[iRun] = iChip; vschan[iRun] = iChan; iRun++; } } for (Int_t j = 0; j < 16; j++) { // c1->cd(j+1); // c1->cd(); sprintf(name, "Run %d Evt %d Time mcm %d chip %d chan %d %s",runnb,iEntry,chip,vschip[j],vschan[j],title); hSignal[j]->SetTitle(name); hSignal[j]->GetYaxis()->SetTitle("Time slice"); hSignal[j]->GetXaxis()->SetTitle("ADC value"); hSignal[j]->SetAxisRange(0, 150, "X"); hSignal[j]->SetAxisRange(0, 1030, "Y"); hSignal[j]->Draw(); if (j == 0) sprintf(name,"chipx-all-chip%d-evt%d-%d-%s.pdf(",chip,iEntry,runnb,title); else if ( j == 15) sprintf(name,"chipx-all-chip%d-evt%d-%d-%s.pdf)",chip,iEntry,runnb,title); else sprintf(name,"chipx-all-chip%d-evt%d-%d-%s.pdf",chip,iEntry,runnb,title); c1->SaveAs(name,"pdf"); // printf("%d\t%5.2f\t%5.3f\t%5.1f\n",j,hSignal[j]->GetMean(),hSignal[j]->GetRMS(),hSignal[j]->GetEntries()); } // sprintf(name,"ch-all-chip%d-time-%s.pdf",chip,title); // c1->SaveAs(name); }