changed SM range, splitted plotting section
authorfronchet <fronchet@f7af4fe6-9843-0410-8265-dc069ae4e863>
Sat, 28 May 2011 17:48:17 +0000 (17:48 +0000)
committerfronchet <fronchet@f7af4fe6-9843-0410-8265-dc069ae4e863>
Sat, 28 May 2011 17:48:17 +0000 (17:48 +0000)
HLT/EMCAL/macros/hlt_plot.C [new file with mode: 0644]
HLT/EMCAL/macros/hlt_rec.C

diff --git a/HLT/EMCAL/macros/hlt_plot.C b/HLT/EMCAL/macros/hlt_plot.C
new file mode 100644 (file)
index 0000000..7f5d51c
--- /dev/null
@@ -0,0 +1,79 @@
+void hlt_plot() {
+  
+  gStyle->SetPalette(1);
+  gStyle->SetOptStat("nemruo");
+
+  TFile *f = new TFile("histofile_local.root");  
+  TProfile2D *h[4];
+  TH2F *g[4];
+
+  TCanvas *c[3];
+
+  Char_t outname[256];
+
+  for (int i=0; i<4; i++) {
+    sprintf(outname, "can%d",i);
+    c[i]= new TCanvas(outname,"",0,0,600,400);
+    c[i]->Divide(2,2);    
+ }
+  // TCanvas *c1 = new TCanvas("c1","",0,0,600,400);
+  // c1->Divide(2,2);
+
+  // TCanvas *c2 = new TCanvas("c2","",0,0,600,400);
+  // c2->Divide(2,2);
+
+  // TCanvas *c3 = new TCanvas("c3","",0,0,600,400);
+  // c3->Divide(2,2);
+  
+  
+  for (int i=0; i<4; i++) {
+
+    c[0]->cd(i+1);
+    h[i] = (TProfile2D*)f->Get(Form("fAmp%d", i));
+    h[i]->SetMaximum(1023);
+    cout << "fAmp" << i << "entries: " << h[i]->GetEntries()<<endl; 
+    h[i]->Draw("colz");
+
+    c[1]->cd(i+1);
+    h[i] = (TProfile2D*)f->Get(Form("fTime%d", i));
+    h[i]->SetMaximum(15);
+    cout << "fTime" << i << "entries: " << h[i]->GetEntries()<<endl;
+    h[i]->Draw("colz");
+    
+    c[2]->cd(i+1);
+    g[i] = (TH2F*)f->Get(Form("fAT%d", i));
+    cout << "fAT" << i << "entries: " << g[i]->GetEntries()<<endl;
+    g[i]->Draw("colz");
+
+
+    c[3]->cd(i+1);
+    g[i] = (TH2F*)f->Get(Form("fCellVsEne%d", i));
+    cout << "fCellVsEne" << i << "entries: " << g[i]->GetEntries()<<endl;
+    g[i]->Draw("colz");
+
+
+  }
+  
+  sprintf(outname, "fAmp.gif");
+  c[0]->SaveAs(outname);
+  
+  sprintf(outname, "fTime.gif");
+  c[1]->SaveAs(outname);
+
+  sprintf(outname, "fAT.gif");
+  c[2]->SaveAs(outname);
+
+  sprintf(outname, "fCellVsEne.gif");
+  c[3]->SaveAs(outname);
+
+
+
+ TCanvas *c4 = new TCanvas("c4","",0,0,600,400);
+ fCellEne->Draw();
+ fClusCellEne->SetLineColor(kRed);
+ fClusCellEne->Draw("same");
+ sprintf(outname, "ClusCellEne.gif");
+ c4->SaveAs(outname);
+
+}
index 221936e..ec35194 100644 (file)
@@ -36,7 +36,7 @@ void hlt_rec(const char* input="./") {
 
   // Define which modules and RCUs to use
   int moduleStart = 0;
-  int moduleEnd = 3;
+  int moduleEnd = 9;
   int rcuStart = 0;
   int rcuEnd = 1;
   
@@ -77,7 +77,7 @@ void hlt_rec(const char* input="./") {
           dm.Form("EMC-DM_%02d_%d", module, rcu);
           arg="";
          // HI-GAIN ~ 15/256 ------ CHECK LOW GAIN FACTOR
-          arg.Form("-highgainfactor 0.057 -lowgainfactor 0.08");
+          arg.Form("-highgainfactor 0.0167 -lowgainfactor 0.08");
           AliHLTConfiguration dmConf(dm.Data(), "EmcalDigitMaker", rawanalyzer.Data(), arg.Data());
           if(clInput.Length() > 0) clInput += " ";
           clInput+=dm;
@@ -90,7 +90,7 @@ void hlt_rec(const char* input="./") {
        cl.Form("EMC-CL_%02d", module);
        arg = "";
        // digitthreshold was 0.005
-       arg.Form("-digitthreshold 0.057 -recpointthreshold 0.1");
+       arg.Form("-digitthreshold 0.0167 -recpointthreshold 0.1");
        AliHLTConfiguration clConf(cl.Data(), "EmcalClusterizer", clInput.Data(), arg.Data());
        if(ecInput.Length() > 0) ecInput += " ";
        ecInput += cl;
@@ -139,7 +139,7 @@ void hlt_rec(const char* input="./") {
   //option+="rootFileHisto loglevel=0x5f";
 
   AliReconstruction rec;
-  Char_t outname[256];
+  
 
   // uncomment for simulation
   rec.SetDefaultStorage("local://$ALICE_ROOT/OCDB");
@@ -161,57 +161,5 @@ void hlt_rec(const char* input="./") {
   
   rec.Run();
 
-  TCanvas *c1 = new TCanvas("c1","",0,0,600,400);
-  c1->Divide(2,2);
-
- // for (int i=0; i<3; i++) {
-//  c1->cd(i+1);
-//  fChannelETMap[i]->Draw("colz");
-
-//  }
-
- c1->cd(1);
- fChannelETMap0->Draw("colz");
- c1->cd(2);
- fChannelETMap1->Draw("colz");
- c1->cd(3);
- fChannelETMap2->Draw("colz");
- c1->cd(4);
- fChannelETMap3->Draw("colz");
-
- sprintf(outname, "ChannleETMap.gif");
- c1->SaveAs(outname);
-
-
- TCanvas *c2 = new TCanvas("c2","",0,0,600,400);
- c2->Divide(2,2);
-
-
- c2->cd(1);
- fChannelEMap0->Draw("colz");
- c2->cd(2);
- fChannelEMap1->Draw("colz");
- c2->cd(3);
- fChannelEMap2->Draw("colz");
- c2->cd(4);
- fChannelEMap3->Draw("colz");
-
- sprintf(outname, "ChannleEMap.gif");
- c2->SaveAs(outname);
-
- TCanvas *c3 = new TCanvas("c3","",0,0,600,400);
- c3->Divide(2,2);
-
- c3->cd(1);
- fChannelTMap0->Draw("colz");
- c3->cd(2);
- fChannelTMap1->Draw("colz");
- c3->cd(3);
- fChannelTMap2->Draw("colz");
- c3->cd(4);
- fChannelTMap3->Draw("colz");
-
- sprintf(outname, "ChannleTMap.gif");
- c3->SaveAs(outname);
 
 }