+ if(fStudyExotic)
+ {
+ fhExoL0ECross = new TH2F("hExoL0_ECross",
+ "#lambda^{2}_{0} vs 1-E_{+}/E_{max} for E > 5 GeV",
+ 400,0,1,ssbins,ssmin,ssmax);
+ fhExoL0ECross ->SetXTitle("1-E_{+}/E_{cell max}");
+ fhExoL0ECross ->SetYTitle("#lambda^{2}_{0}");
+ outputContainer->Add(fhExoL0ECross) ;
+
+ fhExoL1ECross = new TH2F("hExoL1_ECross",
+ "#lambda^{2}_{1} vs 1-E_{+}/E_{max} for E > 5 GeV",
+ 400,0,1,ssbins,ssmin,ssmax);
+ fhExoL1ECross ->SetXTitle("1-E_{+}/E_{cell max}");
+ fhExoL1ECross ->SetYTitle("#lambda^{2}_{1}");
+ outputContainer->Add(fhExoL1ECross) ;
+
+ for(Int_t ie = 0; ie <fExoNECrossCuts; ie++)
+ {
+
+ fhExoDTime[ie] = new TH2F(Form("hExoDTime_ECross%d",ie),
+ Form("#Delta time = t_{max}-t_{cells} vs E_{cluster} for exotic, 1-E_{+}/E_{max} < %2.2f",fExoECrossCuts[ie]),
+ nptbins,ptmin,ptmax,tdbins,tdmin,tdmax);
+ fhExoDTime[ie] ->SetYTitle("#Delta t (ns)");
+ fhExoDTime[ie] ->SetXTitle("E (GeV)");
+ outputContainer->Add(fhExoDTime[ie]) ;
+
+ for(Int_t idt = 0; idt < fExoNDTimeCuts; idt++)
+ {
+ fhExoNCell[ie][idt] = new TH2F(Form("hExoNCell_ECross%d_DT%d",ie,idt),
+ Form("N cells per cluster vs E cluster, 1-E_{+}/E_{max} < %2.2f, #Delta t < %2.0f",fExoECrossCuts[ie],fExoDTimeCuts[idt]),
+ nptbins,ptmin,ptmax,nceclbins,nceclmin,nceclmax);
+ fhExoNCell[ie][idt] ->SetYTitle("N cells");
+ fhExoNCell[ie][idt] ->SetXTitle("E (GeV)");
+ outputContainer->Add(fhExoNCell[ie][idt]) ;
+
+ fhExoL0 [ie][idt] = new TH2F(Form("hExoL0_ECross%d_DT%d",ie,idt),
+ Form("#lambda^{2}_{0} vs E cluster for exotic, 1-E_{+}/E_{max} < %2.2f, #Delta t = %2.0f",fExoECrossCuts[ie],fExoDTimeCuts[idt]),
+ nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
+ fhExoL0 [ie][idt] ->SetYTitle("#lambda^{2}_{0}");
+ fhExoL0 [ie][idt] ->SetXTitle("E (GeV)");
+ outputContainer->Add(fhExoL0[ie][idt]) ;
+
+ fhExoL1 [ie][idt] = new TH2F(Form("hExoL1_ECross%d_DT%d",ie,idt),
+ Form("#lambda^{2}_{1} vs E cluster for exotic, 1-E_{+}/E_{max} < %2.2f, #Delta t = %2.0f",fExoECrossCuts[ie],fExoDTimeCuts[idt]),
+ nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
+ fhExoL1 [ie][idt] ->SetYTitle("#lambda^{2}_{1}");
+ fhExoL1 [ie][idt] ->SetXTitle("E (GeV)");
+ outputContainer->Add(fhExoL1[ie][idt]) ;
+
+ fhExoECross[ie][idt] = new TH2F(Form("hExoECross_ECross%d_DT%d",ie,idt),
+ Form("E cross for cells vs E cell, 1-E_{+}/E_{max} < %2.2f, #Delta t < %2.0f",fExoECrossCuts[ie],fExoDTimeCuts[idt]),
+ nptbins,ptmin,ptmax,400,0,1);
+ fhExoECross[ie][idt] ->SetYTitle("1-E_{+}/E_{cell max}");
+ fhExoECross[ie][idt] ->SetXTitle("E_{cell} (GeV)");
+ outputContainer->Add(fhExoECross[ie][idt]) ;
+
+ fhExoTime [ie][idt] = new TH2F(Form("hExoTime_ECross%d_DT%d",ie,idt),
+ Form("Time of cluster (max cell) vs E cluster for exotic, 1-E_{+}/E_{max} < %2.2f, #Delta t = %2.0f",fExoECrossCuts[ie],fExoDTimeCuts[idt]),
+ nptbins,ptmin,ptmax,ntimebins,timemin,timemax);
+ fhExoTime [ie][idt] ->SetYTitle("time_{max} (ns)");
+ fhExoTime [ie][idt] ->SetXTitle("E (GeV)");
+ outputContainer->Add(fhExoTime[ie][idt]) ;
+
+ fhExoL0NCell[ie][idt] = new TH2F(Form("hExoL0_NCell%d_DT%d",ie,idt),
+ Form("#lambda^{2}_{0} vs N cells per clusters for E > 5 GeV, 1-E_{+}/E_{max} < %2.2f, #Delta t = %2.0f",fExoECrossCuts[ie],fExoDTimeCuts[idt]),
+ nptbins,ptmin,ptmax,ntimebins,timemin,timemax);
+ fhExoL0NCell[ie][idt] ->SetYTitle("N cells");
+ fhExoL0NCell[ie][idt] ->SetXTitle("#lambda^{2}_{0}");
+ outputContainer->Add(fhExoL0NCell[ie][idt]) ;
+
+ fhExoL1NCell[ie][idt] = new TH2F(Form("hExoL1_NCell%d_DT%d",ie,idt),
+ Form("#lambda^{2}_{1} vs N cells per clusters for E > 5 GeV, 1-E_{+}/E_{max} < %2.2f, #Delta t = %2.0f",fExoECrossCuts[ie],fExoDTimeCuts[idt]),
+ nptbins,ptmin,ptmax,ntimebins,timemin,timemax);
+ fhExoL1NCell[ie][idt] ->SetYTitle("N cells");
+ fhExoL1NCell[ie][idt] ->SetXTitle("#lambda^{2}_{1}");
+ outputContainer->Add(fhExoL1NCell[ie][idt]) ;
+
+ }
+ }
+ }
+