bcal_occupancy.C

Go to the documentation of this file.

// hnamepath: /bcal/bcal_fadc_digi_occ
// hnamepath: /bcal/bcal_fadc_digi_occ_layer1
// hnamepath: /bcal/bcal_fadc_digi_occ_layer2
// hnamepath: /bcal/bcal_fadc_digi_occ_layer3
// hnamepath: /bcal/bcal_fadc_digi_occ_layer4
// hnamepath: /bcal/bcal_tdc_digi_occ
// hnamepath: /bcal/bcal_tdc_digi_occ_layer1
// hnamepath: /bcal/bcal_tdc_digi_occ_layer2
// hnamepath: /bcal/bcal_tdc_digi_occ_layer3



{
   gStyle->SetOptStat(111110);   
   gStyle->SetStatX(0.99);
   gStyle->SetStatY(0.99);
   gStyle->SetStatW(0.25);

   gStyle->SetTitleOffset(1, "Y");
   gStyle->SetTitleSize(0.05,"xyz");
   gStyle->SetTitleSize(0.08,"h");
   gStyle->SetLabelSize(0.05,"xyz");
   gStyle->SetTitleX(0);
   gStyle->SetTitleAlign(13);
   gStyle->SetNdivisions(505,"xy");

   int col_layer1 = 1;
   int col_layer2 = 2;
   int col_layer3 = 3;
   int col_layer4 = 4;

   TDirectory *dir = (TDirectory*)gDirectory->FindObjectAny("bcal");
   if(dir) dir->cd();

   TH1I *bcal_fadc_digi_occ = (TH1I*)gDirectory->FindObjectAny("bcal_fadc_digi_occ");
   TH1I *bcal_tdc_digi_occ = (TH1I*)gDirectory->FindObjectAny("bcal_tdc_digi_occ");
   TH1I *bcal_fadc_digi_occ_layer1 = (TH1I*)gDirectory->FindObjectAny("bcal_fadc_digi_occ_layer1");
   TH1I *bcal_fadc_digi_occ_layer2 = (TH1I*)gDirectory->FindObjectAny("bcal_fadc_digi_occ_layer2");
   TH1I *bcal_fadc_digi_occ_layer3 = (TH1I*)gDirectory->FindObjectAny("bcal_fadc_digi_occ_layer3");
   TH1I *bcal_fadc_digi_occ_layer4 = (TH1I*)gDirectory->FindObjectAny("bcal_fadc_digi_occ_layer4");
   
   // Just for testing
   if(gPad == NULL){
      TCanvas *c1 = new TCanvas("c1");
      c1->cd(0);
      c1->Draw();
      c1->Update();
   }
   
   if(!gPad) return;
   TCanvas *c1 = gPad->GetCanvas();
   c1->Divide(2, 2);
   
   c1->cd(1);
   gPad->SetTicks();
   gPad->SetGrid();
   if(bcal_fadc_digi_occ) bcal_fadc_digi_occ->Draw("colz");

   c1->cd(2);
   gPad->SetTicks();
   gPad->SetGrid();
   if(bcal_tdc_digi_occ) bcal_tdc_digi_occ->Draw("colz");

   c1->cd(3);
   gPad->SetTicks();
   gPad->SetGrid();
   THStack *occ_stack = new THStack("occ_stack","Occupancy by layer (DBCALDigiHit);global sector  (4 x module + sector);hits");
   TLegend *occ_stack_legend = new TLegend(0.85,0.75,0.99,0.99);
   if (bcal_fadc_digi_occ_layer1 != NULL) {
      bcal_fadc_digi_occ_layer1->SetLineColor(col_layer1);
      bcal_fadc_digi_occ_layer1->SetMarkerColor(col_layer1);
      bcal_fadc_digi_occ_layer1->SetMarkerStyle(21);
      occ_stack->Add(bcal_fadc_digi_occ_layer1);
      occ_stack_legend->AddEntry(bcal_fadc_digi_occ_layer1,"Layer 1","p");  
   }
   if (bcal_fadc_digi_occ_layer2 != NULL) {
      bcal_fadc_digi_occ_layer2->SetLineColor(col_layer2);
      bcal_fadc_digi_occ_layer2->SetMarkerColor(col_layer2);
      bcal_fadc_digi_occ_layer2->SetMarkerStyle(21);
      occ_stack->Add(bcal_fadc_digi_occ_layer2);
      occ_stack_legend->AddEntry(bcal_fadc_digi_occ_layer2,"Layer 2","p");  
   }
   if (bcal_fadc_digi_occ_layer3 != NULL) {
      bcal_fadc_digi_occ_layer3->SetLineColor(col_layer3);
      bcal_fadc_digi_occ_layer3->SetMarkerColor(col_layer3);
      bcal_fadc_digi_occ_layer3->SetMarkerStyle(21);
      occ_stack->Add(bcal_fadc_digi_occ_layer3);
      occ_stack_legend->AddEntry(bcal_fadc_digi_occ_layer3,"Layer 3","p");  
   }
   if (bcal_fadc_digi_occ_layer4 != NULL) {
      bcal_fadc_digi_occ_layer4->SetLineColor(col_layer4);
      bcal_fadc_digi_occ_layer4->SetMarkerColor(col_layer4);
      bcal_fadc_digi_occ_layer4->SetMarkerStyle(21);
      occ_stack->Add(bcal_fadc_digi_occ_layer4);
      occ_stack_legend->AddEntry(bcal_fadc_digi_occ_layer4,"Layer 4","p");  
   }
   occ_stack->Draw("nostack");
   occ_stack_legend->Draw();




   c1->cd(4);
   gPad->SetTicks();
   gPad->SetGrid();
   TH1I *bcal_tdc_digi_occ_layer1 = (TH1I*)gROOT->FindObject("bcal_tdc_digi_occ_layer1");
   TH1I *bcal_tdc_digi_occ_layer2 = (TH1I*)gROOT->FindObject("bcal_tdc_digi_occ_layer2");
   TH1I *bcal_tdc_digi_occ_layer3 = (TH1I*)gROOT->FindObject("bcal_tdc_digi_occ_layer3");
   THStack *tdc_occ_stack = new THStack("tdc_occ_stack","TDC Occupancy by layer (DBCALDigiHit);global sector  (4 x module + sector);hits");
   TLegend *tdc_occ_stack_legend = new TLegend(0.85,0.75,0.99,0.99);
   if (bcal_tdc_digi_occ_layer1 != NULL) {
      bcal_tdc_digi_occ_layer1->SetLineColor(col_layer1);
      bcal_tdc_digi_occ_layer1->SetMarkerColor(col_layer1);
      bcal_tdc_digi_occ_layer1->SetMarkerStyle(21);
      tdc_occ_stack->Add(bcal_tdc_digi_occ_layer1);    
      tdc_occ_stack_legend->AddEntry(bcal_tdc_digi_occ_layer1,"Layer 1","p");  
   }
   if (bcal_tdc_digi_occ_layer2 != NULL) {
      bcal_tdc_digi_occ_layer2->SetLineColor(col_layer2);
      bcal_tdc_digi_occ_layer2->SetMarkerColor(col_layer2);
      bcal_tdc_digi_occ_layer2->SetMarkerStyle(21);
      tdc_occ_stack->Add(bcal_tdc_digi_occ_layer2);
      tdc_occ_stack_legend->AddEntry(bcal_tdc_digi_occ_layer2,"Layer 2","p");  
   }
   if (bcal_tdc_digi_occ_layer3 != NULL) {
      bcal_tdc_digi_occ_layer3->SetLineColor(col_layer3);
      bcal_tdc_digi_occ_layer3->SetMarkerColor(col_layer3);
      bcal_tdc_digi_occ_layer3->SetMarkerStyle(21);
      tdc_occ_stack->Add(bcal_tdc_digi_occ_layer3);
      tdc_occ_stack_legend->AddEntry(bcal_tdc_digi_occ_layer3,"Layer 3","p");  
   }
   tdc_occ_stack->Draw("nostack");
   tdc_occ_stack_legend->Draw();

}