[u/mrichter/AliRoot.git] / HLT / EVE / AliHLTEveMultCorr.cxx

#include "AliHLTEveMultCorr.h"
#include "TH1F.h"
#include "TH2F.h"
#include <iostream>
#include <TCanvas.h>
#include "AliHLTHOMERBlockDesc.h"
#include "AliHLTEveHistoMerger.h"
#include "TList.h"
#include "AliLog.h"

AliHLTEveMultCorr::AliHLTEveMultCorr(const char* name): AliHLTEveBase(name)
        ,fVzeroMultCanvas(0)
        ,fZdcMultCanvas(0)
        ,fTrackMultCanvas(0)
        ,fCorrCanvas(0)
        ,fEtCorrCanvas(0)
        ,fZdcVzeroSpdCorrCanvas(0)
        ,fMerger(0)
        ,fMyList(0)
{
  fMerger = new AliHLTEveHistoMerger();
}

AliHLTEveMultCorr::~AliHLTEveMultCorr()
{
  if(fMerger)
  {
    delete fMerger;
  }
  if(fMyList)
  {
    delete fMyList;
  }
}

void AliHLTEveMultCorr::ResetElements()
{

}

void AliHLTEveMultCorr::UpdateElements()
{
    fVzeroMultCanvas->Update();
    fZdcMultCanvas->Update();
    fTrackMultCanvas->Update();
    fCorrCanvas->Update();
    fEtCorrCanvas->Update();
    fZdcVzeroSpdCorrCanvas->Update();

}

void AliHLTEveMultCorr::ProcessBlock(AliHLTHOMERBlockDesc* block)
{
    TString msg;
    
    TList *hlist = dynamic_cast<TList*>(block->GetTObject());
    if(kFALSE)
    {
      hlist = dynamic_cast<TList*>(fMerger->Process(hlist, block->GetSpecification()));
    }
    else
    {
      if(fMyList)
      {
	delete fMyList;
	fMyList = 0;
      }
      fMyList = dynamic_cast<TList*>(hlist->Clone());
      hlist = fMyList;
    }
    if (hlist)
    {

        if (hlist->Contains("fVzeroMult")) // These are the correlation histograms
        {
            TH1F *oneDf = 0;
            TH2F *twoDf = 0;
	    int cd = 1;

            //VZERO multiplicity hists
            if (!fVzeroMultCanvas) 
	    {
	      fVzeroMultCanvas = CreateCanvas("V0 M", "V0 Multiplicities");
	      fVzeroMultCanvas->Divide(2, 2);
	    }
            oneDf = dynamic_cast<TH1F*>(FindHistogram(hlist, "fVzeroMult"));
            AddHistogramToCanvas(oneDf, fVzeroMultCanvas, cd);
	    cd++;
            twoDf = dynamic_cast<TH2F*>(FindHistogram(hlist, "fVzeroMultAC"));
            AddHistogramToCanvas(twoDf, fVzeroMultCanvas, cd);
	    cd++;
            oneDf = dynamic_cast<TH1F*>(FindHistogram(hlist, "fVzeroFlaggedMult"));
            AddHistogramToCanvas(oneDf, fVzeroMultCanvas, cd);
	    cd++;
            twoDf = dynamic_cast<TH2F*>(FindHistogram(hlist, "fVzeroFlaggedMultAC"));
            AddHistogramToCanvas(twoDf, fVzeroMultCanvas, cd);

	    cd = 1;
            //ZDC multiplicity hists
            if (!fZdcMultCanvas) 
	    {
	      fZdcMultCanvas = CreateCanvas("ZDC M", "ZDC Multiplicities");
	      fZdcMultCanvas->Divide(3, 3);
	    }
            oneDf = dynamic_cast<TH1F*>(FindHistogram(hlist, "fZdcEzdc"));
            AddHistogramToCanvas(oneDf, fZdcMultCanvas, cd);
	    cd++;
            oneDf = dynamic_cast<TH1F*>(FindHistogram(hlist, "fZdcEzn"));
            AddHistogramToCanvas(oneDf, fZdcMultCanvas, cd);
	    cd++;
            oneDf = dynamic_cast<TH1F*>(FindHistogram(hlist, "fZdcEzp"));
            AddHistogramToCanvas(oneDf, fZdcMultCanvas, cd);
	    cd++;
            oneDf = dynamic_cast<TH1F*>(FindHistogram(hlist, "fZdcEzem"));
            AddHistogramToCanvas(oneDf, fZdcMultCanvas, cd);
	    cd++;
            oneDf = dynamic_cast<TH1F*>(FindHistogram(hlist, "fZdcNpart"));
            AddHistogramToCanvas(oneDf, fZdcMultCanvas, cd);
	    cd++;
            oneDf = dynamic_cast<TH1F*>(FindHistogram(hlist, "fZdcB"));
            AddHistogramToCanvas(oneDf, fZdcMultCanvas, cd);
	    cd++;
            twoDf = dynamic_cast<TH2F*>(FindHistogram(hlist, "fZdcEzemEzdc"));
            AddHistogramToCanvas(twoDf, fZdcMultCanvas, cd);

	    cd = 1;
            // TPC multiplicity
            if (!fTrackMultCanvas) 
	    {
	      fTrackMultCanvas = CreateCanvas("Track M", "Track Multiplicities");
	      fTrackMultCanvas->Divide(2, 2);
	    }
            oneDf = dynamic_cast<TH1F*>(FindHistogram(hlist, "fTpcNch2"));
            AddHistogramToCanvas(oneDf, fTrackMultCanvas, cd);
	    cd++;
            oneDf = dynamic_cast<TH1F*>(FindHistogram(hlist, "fTpcNch3"));
            AddHistogramToCanvas(oneDf, fTrackMultCanvas, cd);
	    cd++;
	    oneDf = dynamic_cast<TH1F*>(FindHistogram(hlist, "fSpdNClusters"));
	    
            AddHistogramToCanvas(oneDf, fTrackMultCanvas, cd);

	    cd = 1;
            // Correlations
            if (!fCorrCanvas) 
	    {
	      fCorrCanvas = CreateCanvas("Corr", "Multiplicity Correlations");
	      fCorrCanvas->Divide(3, 2);
	    }
            twoDf = dynamic_cast<TH2F*>(FindHistogram(hlist, "fCorrEzdcNch"));
            AddHistogramToCanvas(twoDf, fCorrCanvas, cd);
	    cd++;
            twoDf = dynamic_cast<TH2F*>(FindHistogram(hlist, "fCorrVzeroNch"));
            AddHistogramToCanvas(twoDf, fCorrCanvas, cd);
	    cd++;
            twoDf = dynamic_cast<TH2F*>(FindHistogram(hlist, "fCorrSpdTpcNch"));
            AddHistogramToCanvas(twoDf, fCorrCanvas, cd);
	    cd++;
            twoDf = dynamic_cast<TH2F*>(FindHistogram(hlist, "fCorrEzemNch"));
            AddHistogramToCanvas(twoDf, fCorrCanvas, cd);
	    cd++;
            twoDf = dynamic_cast<TH2F*>(FindHistogram(hlist, "fCorrEzdcVzero"));
            AddHistogramToCanvas(twoDf, fCorrCanvas, cd);
	    cd++;
            twoDf = dynamic_cast<TH2F*>(FindHistogram(hlist, "fCorrEzemVzero"));
            AddHistogramToCanvas(twoDf, fCorrCanvas, cd);
	    
	    cd = 1;
            // ET Correlations
            if (!fEtCorrCanvas) 
	    {
	      fEtCorrCanvas = CreateCanvas("ET", "E_{T} Correlations");
	      fEtCorrCanvas->Divide(3, 2);
	    }
            twoDf = dynamic_cast<TH2F*>(FindHistogram(hlist, "fCorrZdcTotEvsPhosTotEt"));
            AddHistogramToCanvas(twoDf, fEtCorrCanvas, cd);
	    cd++;
            twoDf = dynamic_cast<TH2F*>(FindHistogram(hlist, "fCorrZdcTotEvsEmcalTotEt"));
            AddHistogramToCanvas(twoDf, fEtCorrCanvas, cd);
	    cd++;
            twoDf = dynamic_cast<TH2F*>(FindHistogram(hlist, "fCorrZdcTotEvsTotEt"));
            AddHistogramToCanvas(twoDf, fEtCorrCanvas, cd);
	    cd++;
            twoDf = dynamic_cast<TH2F*>(FindHistogram(hlist, "fCorrVzerovsPhosTotEt"));
            AddHistogramToCanvas(twoDf, fEtCorrCanvas, cd);
	    cd++;
            twoDf = dynamic_cast<TH2F*>(FindHistogram(hlist, "fCorrVzerovsEmcalTotEt"));
            AddHistogramToCanvas(twoDf, fEtCorrCanvas, cd);
	    cd++;
            twoDf = dynamic_cast<TH2F*>(FindHistogram(hlist, "fCorrVzerovsTotEt"));
            AddHistogramToCanvas(twoDf, fEtCorrCanvas, cd);

	    cd = 1;
            // ZDC,V0 vs SPD
            if (!fZdcVzeroSpdCorrCanvas) 
	    {
	      fZdcVzeroSpdCorrCanvas = CreateCanvas("ZDC/V0 vs SPD", "ZDC/V0 vs. SPD");
	      fZdcVzeroSpdCorrCanvas->Divide(3, 2);
	      
	    }
            twoDf = dynamic_cast<TH2F*>(FindHistogram(hlist, "fCorrVzeroSpd"));
            AddHistogramToCanvas(twoDf, fZdcVzeroSpdCorrCanvas, cd);
	    cd++;
            twoDf = dynamic_cast<TH2F*>(FindHistogram(hlist, "fCorrVzeroASpd"));
            AddHistogramToCanvas(twoDf, fZdcVzeroSpdCorrCanvas, cd);
	    cd++;
            twoDf = dynamic_cast<TH2F*>(FindHistogram(hlist, "fCorrVzeroCSpd"));
            AddHistogramToCanvas(twoDf, fZdcVzeroSpdCorrCanvas, cd);
	    cd++;
            twoDf = dynamic_cast<TH2F*>(FindHistogram(hlist, "fCorrEzdcSpd"));
            AddHistogramToCanvas(twoDf, fZdcVzeroSpdCorrCanvas, cd);
	    cd++;
            twoDf = dynamic_cast<TH2F*>(FindHistogram(hlist, "fCorrEzdcASpd"));
            AddHistogramToCanvas(twoDf, fZdcVzeroSpdCorrCanvas, cd);
	    cd++;
            twoDf = dynamic_cast<TH2F*>(FindHistogram(hlist, "fCorrEzdcCSpd"));
            AddHistogramToCanvas(twoDf, fZdcVzeroSpdCorrCanvas, cd);
        }
        else
        {
            AliWarning("This block does not contain what you think it contains!");
        }
    }
    else
    {
      AliWarning("There is no TList object in the block");
    }


}

void AliHLTEveMultCorr::AddHistogramsToCanvas(AliHLTHOMERBlockDesc* /*block*/, TCanvas* /*canvas*/, Int_t& /*cdCount*/)
{
}

void AliHLTEveMultCorr::AddHistogramToCanvas(TH1* hist, TCanvas* canvas, Int_t& cdCount, Bool_t zoom)
{
  
  TString msg;
    if (hist)
    {
        if (!strcmp(hist->ClassName(), "TH1F"))
        {
	    canvas->cd(cdCount);
	    TPad * pad = dynamic_cast<TPad*>(canvas->cd(cdCount));
	    pad->SetLogy();
	    if(zoom) 
	    {
	      TH1F *h = dynamic_cast<TH1F*>(hist);
	      h->GetXaxis()->SetRange(0, h->GetMaximumBin() + h->GetMaximumBin()*0.2);
	    }
            dynamic_cast<TH1F*>(hist)->Draw();
        }
        else if (!strcmp(hist->ClassName(), "TH2F"))
        {
	    canvas->cd(cdCount);
	    if(zoom)
	    {
	      TH2F *h = dynamic_cast<TH2F*>(hist);
	      h->GetXaxis()->SetRange(0, h->GetMaximumBin() + h->GetMaximumBin()*0.2);
	      h->GetYaxis()->SetRange(0, h->GetMaximumBin() + h->GetMaximumBin()*0.2);
	    }
            dynamic_cast<TH2F*>(hist)->Draw("COLZ");
        }
        else
        {	
	  msg.Form("I don't want histograms of type %s", hist->ClassName());
	  AliWarning(msg.Data());
        }
    }
    else
    {
      AliWarning("You gave me a null pointer!");
    }
}

TH1* AliHLTEveMultCorr::FindHistogram(TCollection* coll, const char* name)
{
  TString msg;
  TH1 *hist = dynamic_cast<TH1*>(coll->FindObject(name));
  if(!hist)
  {
    msg.Form("Could not find object %s", name);
    AliWarning(msg);
  }
  return hist;
}
Commit	Line	Data
2127802b	1	#include "AliHLTEveMultCorr.h"
	2	#include "TH1F.h"
	3	#include "TH2F.h"
	4	#include <iostream>
	5	#include <TCanvas.h>
	6	#include "AliHLTHOMERBlockDesc.h"
0ff4c5eb	7	#include "AliHLTEveHistoMerger.h"
2127802b	8	#include "TList.h"
0ff4c5eb	9	#include "AliLog.h"
2127802b	10
	11	AliHLTEveMultCorr::AliHLTEveMultCorr(const char* name): AliHLTEveBase(name)
	12	,fVzeroMultCanvas(0)
	13	,fZdcMultCanvas(0)
0ff4c5eb	14	,fTrackMultCanvas(0)
2127802b	15	,fCorrCanvas(0)
	16	,fEtCorrCanvas(0)
	17	,fZdcVzeroSpdCorrCanvas(0)
0ff4c5eb	18	,fMerger(0)
0ff4c5eb	19	,fMyList(0)
2127802b	20	{
0ff4c5eb	21	fMerger = new AliHLTEveHistoMerger();
2127802b	22	}
	23
	24	AliHLTEveMultCorr::~AliHLTEveMultCorr()
	25	{
0ff4c5eb	26	if(fMerger)
	27	{
	28	delete fMerger;
	29	}
	30	if(fMyList)
	31	{
	32	delete fMyList;
	33	}
2127802b	34	}
	35
	36	void AliHLTEveMultCorr::ResetElements()
	37	{
	38
	39	}
	40
	41	void AliHLTEveMultCorr::UpdateElements()
	42	{
	43	fVzeroMultCanvas->Update();
	44	fZdcMultCanvas->Update();
0ff4c5eb	45	fTrackMultCanvas->Update();
2127802b	46	fCorrCanvas->Update();
	47	fEtCorrCanvas->Update();
	48	fZdcVzeroSpdCorrCanvas->Update();
	49
	50	}
	51
	52	void AliHLTEveMultCorr::ProcessBlock(AliHLTHOMERBlockDesc* block)
	53	{
0ff4c5eb	54	TString msg;
0ff4c5eb	55
2127802b	56	TList hlist = dynamic_cast<TList>(block->GetTObject());
0ff4c5eb	57	if(kFALSE)
	58	{
	59	hlist = dynamic_cast<TList*>(fMerger->Process(hlist, block->GetSpecification()));
	60	}
	61	else
	62	{
	63	if(fMyList)
	64	{
	65	delete fMyList;
	66	fMyList = 0;
	67	}
	68	fMyList = dynamic_cast<TList*>(hlist->Clone());
	69	hlist = fMyList;
	70	}
2127802b	71	if (hlist)
	72	{
	73
	74	if (hlist->Contains("fVzeroMult")) // These are the correlation histograms
	75	{
	76	TH1F *oneDf = 0;
	77	TH2F *twoDf = 0;
	78	int cd = 1;
	79
	80	//VZERO multiplicity hists
	81	if (!fVzeroMultCanvas)
	82	{
0ff4c5eb	83	fVzeroMultCanvas = CreateCanvas("V0 M", "V0 Multiplicities");
2127802b	84	fVzeroMultCanvas->Divide(2, 2);
2127802b	85	}
0ff4c5eb	86	oneDf = dynamic_cast<TH1F*>(FindHistogram(hlist, "fVzeroMult"));
2127802b	87	AddHistogramToCanvas(oneDf, fVzeroMultCanvas, cd);
2127802b	88	cd++;
0ff4c5eb	89	twoDf = dynamic_cast<TH2F*>(FindHistogram(hlist, "fVzeroMultAC"));
2127802b	90	AddHistogramToCanvas(twoDf, fVzeroMultCanvas, cd);
2127802b	91	cd++;
0ff4c5eb	92	oneDf = dynamic_cast<TH1F*>(FindHistogram(hlist, "fVzeroFlaggedMult"));
2127802b	93	AddHistogramToCanvas(oneDf, fVzeroMultCanvas, cd);
2127802b	94	cd++;
0ff4c5eb	95	twoDf = dynamic_cast<TH2F*>(FindHistogram(hlist, "fVzeroFlaggedMultAC"));
2127802b	96	AddHistogramToCanvas(twoDf, fVzeroMultCanvas, cd);
	97
	98	cd = 1;
	99	//ZDC multiplicity hists
	100	if (!fZdcMultCanvas)
	101	{
0ff4c5eb	102	fZdcMultCanvas = CreateCanvas("ZDC M", "ZDC Multiplicities");
0ff4c5eb	103	fZdcMultCanvas->Divide(3, 3);
2127802b	104	}
0ff4c5eb	105	oneDf = dynamic_cast<TH1F*>(FindHistogram(hlist, "fZdcEzdc"));
2127802b	106	AddHistogramToCanvas(oneDf, fZdcMultCanvas, cd);
2127802b	107	cd++;
0ff4c5eb	108	oneDf = dynamic_cast<TH1F*>(FindHistogram(hlist, "fZdcEzn"));
2127802b	109	AddHistogramToCanvas(oneDf, fZdcMultCanvas, cd);
2127802b	110	cd++;
0ff4c5eb	111	oneDf = dynamic_cast<TH1F*>(FindHistogram(hlist, "fZdcEzp"));
2127802b	112	AddHistogramToCanvas(oneDf, fZdcMultCanvas, cd);
2127802b	113	cd++;
0ff4c5eb	114	oneDf = dynamic_cast<TH1F*>(FindHistogram(hlist, "fZdcEzem"));
2127802b	115	AddHistogramToCanvas(oneDf, fZdcMultCanvas, cd);
2127802b	116	cd++;
0ff4c5eb	117	oneDf = dynamic_cast<TH1F*>(FindHistogram(hlist, "fZdcNpart"));
2127802b	118	AddHistogramToCanvas(oneDf, fZdcMultCanvas, cd);
2127802b	119	cd++;
0ff4c5eb	120	oneDf = dynamic_cast<TH1F*>(FindHistogram(hlist, "fZdcB"));
2127802b	121	AddHistogramToCanvas(oneDf, fZdcMultCanvas, cd);
2127802b	122	cd++;
0ff4c5eb	123	twoDf = dynamic_cast<TH2F*>(FindHistogram(hlist, "fZdcEzemEzdc"));
2127802b	124	AddHistogramToCanvas(twoDf, fZdcMultCanvas, cd);
	125
	126	cd = 1;
	127	// TPC multiplicity
0ff4c5eb	128	if (!fTrackMultCanvas)
2127802b	129	{
0ff4c5eb	130	fTrackMultCanvas = CreateCanvas("Track M", "Track Multiplicities");
0ff4c5eb	131	fTrackMultCanvas->Divide(2, 2);
2127802b	132	}
0ff4c5eb	133	oneDf = dynamic_cast<TH1F*>(FindHistogram(hlist, "fTpcNch2"));
0ff4c5eb	134	AddHistogramToCanvas(oneDf, fTrackMultCanvas, cd);
2127802b	135	cd++;
0ff4c5eb	136	oneDf = dynamic_cast<TH1F*>(FindHistogram(hlist, "fTpcNch3"));
0ff4c5eb	137	AddHistogramToCanvas(oneDf, fTrackMultCanvas, cd);
9e34efc3	138	cd++;
0ff4c5eb	139	oneDf = dynamic_cast<TH1F*>(FindHistogram(hlist, "fSpdNClusters"));
9e34efc3	140
0ff4c5eb	141	AddHistogramToCanvas(oneDf, fTrackMultCanvas, cd);
2127802b	142
	143	cd = 1;
	144	// Correlations
	145	if (!fCorrCanvas)
	146	{
0ff4c5eb	147	fCorrCanvas = CreateCanvas("Corr", "Multiplicity Correlations");
0ff4c5eb	148	fCorrCanvas->Divide(3, 2);
2127802b	149	}
0ff4c5eb	150	twoDf = dynamic_cast<TH2F*>(FindHistogram(hlist, "fCorrEzdcNch"));
2127802b	151	AddHistogramToCanvas(twoDf, fCorrCanvas, cd);
2127802b	152	cd++;
0ff4c5eb	153	twoDf = dynamic_cast<TH2F*>(FindHistogram(hlist, "fCorrVzeroNch"));
2127802b	154	AddHistogramToCanvas(twoDf, fCorrCanvas, cd);
2127802b	155	cd++;
0ff4c5eb	156	twoDf = dynamic_cast<TH2F*>(FindHistogram(hlist, "fCorrSpdTpcNch"));
2127802b	157	AddHistogramToCanvas(twoDf, fCorrCanvas, cd);
2127802b	158	cd++;
0ff4c5eb	159	twoDf = dynamic_cast<TH2F*>(FindHistogram(hlist, "fCorrEzemNch"));
2127802b	160	AddHistogramToCanvas(twoDf, fCorrCanvas, cd);
2127802b	161	cd++;
0ff4c5eb	162	twoDf = dynamic_cast<TH2F*>(FindHistogram(hlist, "fCorrEzdcVzero"));
2127802b	163	AddHistogramToCanvas(twoDf, fCorrCanvas, cd);
2127802b	164	cd++;
0ff4c5eb	165	twoDf = dynamic_cast<TH2F*>(FindHistogram(hlist, "fCorrEzemVzero"));
2127802b	166	AddHistogramToCanvas(twoDf, fCorrCanvas, cd);
	167
	168	cd = 1;
	169	// ET Correlations
	170	if (!fEtCorrCanvas)
	171	{
0ff4c5eb	172	fEtCorrCanvas = CreateCanvas("ET", "E_{T} Correlations");
0ff4c5eb	173	fEtCorrCanvas->Divide(3, 2);
2127802b	174	}
0ff4c5eb	175	twoDf = dynamic_cast<TH2F*>(FindHistogram(hlist, "fCorrZdcTotEvsPhosTotEt"));
2127802b	176	AddHistogramToCanvas(twoDf, fEtCorrCanvas, cd);
2127802b	177	cd++;
0ff4c5eb	178	twoDf = dynamic_cast<TH2F*>(FindHistogram(hlist, "fCorrZdcTotEvsEmcalTotEt"));
2127802b	179	AddHistogramToCanvas(twoDf, fEtCorrCanvas, cd);
2127802b	180	cd++;
0ff4c5eb	181	twoDf = dynamic_cast<TH2F*>(FindHistogram(hlist, "fCorrZdcTotEvsTotEt"));
2127802b	182	AddHistogramToCanvas(twoDf, fEtCorrCanvas, cd);
2127802b	183	cd++;
0ff4c5eb	184	twoDf = dynamic_cast<TH2F*>(FindHistogram(hlist, "fCorrVzerovsPhosTotEt"));
2127802b	185	AddHistogramToCanvas(twoDf, fEtCorrCanvas, cd);
2127802b	186	cd++;
0ff4c5eb	187	twoDf = dynamic_cast<TH2F*>(FindHistogram(hlist, "fCorrVzerovsEmcalTotEt"));
2127802b	188	AddHistogramToCanvas(twoDf, fEtCorrCanvas, cd);
2127802b	189	cd++;
0ff4c5eb	190	twoDf = dynamic_cast<TH2F*>(FindHistogram(hlist, "fCorrVzerovsTotEt"));
2127802b	191	AddHistogramToCanvas(twoDf, fEtCorrCanvas, cd);
	192
	193	cd = 1;
	194	// ZDC,V0 vs SPD
	195	if (!fZdcVzeroSpdCorrCanvas)
	196	{
	197	fZdcVzeroSpdCorrCanvas = CreateCanvas("ZDC/V0 vs SPD", "ZDC/V0 vs. SPD");
0ff4c5eb	198	fZdcVzeroSpdCorrCanvas->Divide(3, 2);
2127802b	199
2127802b	200	}
0ff4c5eb	201	twoDf = dynamic_cast<TH2F*>(FindHistogram(hlist, "fCorrVzeroSpd"));
2127802b	202	AddHistogramToCanvas(twoDf, fZdcVzeroSpdCorrCanvas, cd);
2127802b	203	cd++;
0ff4c5eb	204	twoDf = dynamic_cast<TH2F*>(FindHistogram(hlist, "fCorrVzeroASpd"));
2127802b	205	AddHistogramToCanvas(twoDf, fZdcVzeroSpdCorrCanvas, cd);
2127802b	206	cd++;
0ff4c5eb	207	twoDf = dynamic_cast<TH2F*>(FindHistogram(hlist, "fCorrVzeroCSpd"));
2127802b	208	AddHistogramToCanvas(twoDf, fZdcVzeroSpdCorrCanvas, cd);
2127802b	209	cd++;
0ff4c5eb	210	twoDf = dynamic_cast<TH2F*>(FindHistogram(hlist, "fCorrEzdcSpd"));
2127802b	211	AddHistogramToCanvas(twoDf, fZdcVzeroSpdCorrCanvas, cd);
2127802b	212	cd++;
0ff4c5eb	213	twoDf = dynamic_cast<TH2F*>(FindHistogram(hlist, "fCorrEzdcASpd"));
2127802b	214	AddHistogramToCanvas(twoDf, fZdcVzeroSpdCorrCanvas, cd);
2127802b	215	cd++;
0ff4c5eb	216	twoDf = dynamic_cast<TH2F*>(FindHistogram(hlist, "fCorrEzdcCSpd"));
2127802b	217	AddHistogramToCanvas(twoDf, fZdcVzeroSpdCorrCanvas, cd);
	218	}
	219	else
	220	{
0ff4c5eb	221	AliWarning("This block does not contain what you think it contains!");
2127802b	222	}
	223	}
	224	else
	225	{
0ff4c5eb	226	AliWarning("There is no TList object in the block");
2127802b	227	}
	228
	229
	230	}
	231
	232	void AliHLTEveMultCorr::AddHistogramsToCanvas(AliHLTHOMERBlockDesc* /block/, TCanvas* /canvas/, Int_t& /cdCount/)
	233	{
	234	}
	235
0ff4c5eb	236	void AliHLTEveMultCorr::AddHistogramToCanvas(TH1* hist, TCanvas* canvas, Int_t& cdCount, Bool_t zoom)
2127802b	237	{
0ff4c5eb	238
0ff4c5eb	239	TString msg;
2127802b	240	if (hist)
	241	{
	242	if (!strcmp(hist->ClassName(), "TH1F"))
	243	{
	244	canvas->cd(cdCount);
9e34efc3	245	TPad * pad = dynamic_cast<TPad*>(canvas->cd(cdCount));
9e34efc3	246	pad->SetLogy();
0ff4c5eb	247	if(zoom)
	248	{
	249	TH1F h = dynamic_cast<TH1F>(hist);
	250	h->GetXaxis()->SetRange(0, h->GetMaximumBin() + h->GetMaximumBin()*0.2);
	251	}
2127802b	252	dynamic_cast<TH1F*>(hist)->Draw();
	253	}
	254	else if (!strcmp(hist->ClassName(), "TH2F"))
	255	{
	256	canvas->cd(cdCount);
0ff4c5eb	257	if(zoom)
	258	{
	259	TH2F h = dynamic_cast<TH2F>(hist);
	260	h->GetXaxis()->SetRange(0, h->GetMaximumBin() + h->GetMaximumBin()*0.2);
	261	h->GetYaxis()->SetRange(0, h->GetMaximumBin() + h->GetMaximumBin()*0.2);
	262	}
2127802b	263	dynamic_cast<TH2F*>(hist)->Draw("COLZ");
	264	}
	265	else
0ff4c5eb	266	{
	267	msg.Form("I don't want histograms of type %s", hist->ClassName());
	268	AliWarning(msg.Data());
2127802b	269	}
	270	}
	271	else
	272	{
0ff4c5eb	273	AliWarning("You gave me a null pointer!");
2127802b	274	}
	275	}
	276
0ff4c5eb	277	TH1* AliHLTEveMultCorr::FindHistogram(TCollection* coll, const char* name)
	278	{
	279	TString msg;
	280	TH1 hist = dynamic_cast<TH1>(coll->FindObject(name));
	281	if(!hist)
	282	{
	283	msg.Form("Could not find object %s", name);
	284	AliWarning(msg);
	285	}
	286	return hist;
	287	}