void checkOpenHLT::fillAllHistograms(const TString & name, const int arraySize, const bool * selectionArray)
{
fillHistograms("pt"+name, ohMuL2NoVtxPt, arraySize, selectionArray);
fillHistograms("eta"+name, ohMuL2NoVtxEta, arraySize, selectionArray);
fillHistograms("phi"+name, ohMuL2NoVtxPhi, arraySize, selectionArray);
// fillHistograms("chg"+name, ohMuL2NoVtxChg, arraySize, selectionArray);
fillHistograms("nhits"+name, ohMuL2NoVtxNhits, arraySize, selectionArray);
fillHistograms("nchambers"+name, ohMuL2NoVtxNchambers, arraySize, selectionArray);
if( arraySize > 1 && selectionArray[0] && selectionArray[1] ) histoMap_["parallelism"+name]->Fill(parallelDiff_);
}
void Fit3D::ntupleLoopCore(const int& entry_id)
{
x_value_ = l0_pcm + l1_pcm;
y_value_ = cos_thta;
if (evt_sign != 0) {
z_value_ = abs(l0_ip_dz - l1_ip_dz);
} else {
z_value_ = l0_ip_dz - l1_ip_dz;
}
AnalysisSelectors cuts(*this);
fillDataSet(cuts.signal_type());
fillHistograms(cuts.signal_type());
}
void CombineSpectra_minbias()
{
gROOT->Reset();
gROOT->ProcessLine(".x rootlogonChristof.C");
gROOT->ForceStyle();
gStyle->SetPalette(1);
gStyle->SetTitleYOffset(1.30);
gStyle->SetPadRightMargin(0.16);
bool doSave = true;
TFile * f_output = new TFile("CombineSpectra_minbias.root","recreate");
// Define the input file and HiForest
char *infName_1 = "root://eoscms//eos/cms/store/group/phys_heavyions/yjlee/pPb2013/promptReco/PA2013_HiForest_PromptReco_HLT_SingleTrack_JSonPPb_forestv84.root";
// char *infName_2 = "root://eoscms//eos/cms/store/group/phys_heavyions/kjung/MinBiasUPCForest_v71/MergedForest_withCones_MinBiasUPC_v71_2.root";
HiForest *c_minbias_1 = new HiForest(infName_1,"",cPPb);
// HiForest *c_minbias_2 = new HiForest(infName_2,"",cPPb);
// c_minbias_1->doTrackCorrections=1;
// c_minbias_1->InitTree();
std::vector<Double_t> ptBins;
const Double_t small = 1e-3;
Double_t pt;
for(pt = 0. ; pt < 48.-small; pt += 2./2.) ptBins.push_back(pt);
for(pt = 48.; pt < 60.-small; pt += 4./2. ) ptBins.push_back(pt);
for(pt = 60.; pt < 76.-small; pt += 8./2. ) ptBins.push_back(pt);
for(pt = 76.; pt < 140.-small; pt += 16./2. ) ptBins.push_back(pt);
for(pt = 140.; pt < 364.-small; pt += 32./2. ) ptBins.push_back(pt);
ptBins.push_back(364.0);
std::vector<Double_t> ptBins_smaller;
Double_t pt_smaller;
for(pt_smaller = 0.2; pt_smaller < 20.0-small; pt_smaller += 0.1 ) ptBins_smaller.push_back(pt_smaller);
for(pt_smaller = 20.0; pt_smaller < 40.0-small; pt_smaller += 0.2 ) ptBins_smaller.push_back(pt_smaller);
for(pt_smaller = 40.0; pt_smaller < 103.6-small; pt_smaller += 0.4 ) ptBins_smaller.push_back(pt_smaller);
ptBins_smaller.push_back(103.6);
/*
std::vector<Double_t> ptBins_part;
Double_t pt_part;
for(pt_part = 0.2 ; pt_part < 1.2-small; pt_part += 0.05) ptBins_part.push_back(pt_part);
for(pt_part = 1.2; pt_part < 2.4-small; pt_part += 0.1 ) ptBins_part.push_back(pt_part); //12 bins
for(pt_part = 2.4; pt_part < 7.2-small; pt_part += 0.4 ) ptBins_part.push_back(pt_part); //12 bins
for(pt_part = 7.2; pt_part < 16.8-small; pt_part += 1.6 ) ptBins_part.push_back(pt_part);//it was 3.2
for(pt_part = 16.8; pt_part < 30.0-small; pt_part += 6.6 ) ptBins_part.push_back(pt_part);
for(pt_part = 30.0; pt_part < 100.0-small; pt_part += 17.5 ) ptBins_part.push_back(pt_part);
for(pt_part = 100.0; pt_part < 300.0-small; pt_part += 25. ) ptBins_part.push_back(pt_part);
ptBins_part.push_back(300.0);
*/
//orig
// Int_t NumOfPtBins_part = 34;
// double ptBins_part[] = {0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0, 1.1, 1.2, 1.4, 1.6, 1.8, 2.0, 2.2, 2.4, 3.2, 4.0, 4.8, 5.6, 6.4, 7.2, 9.6, 12.0, 14.4, 19.2, 24.0, 28.8, 35.2, 41.6, 48.0, 60.8, 73.6, 86.4, 103.6};//, 120.8, 138., 155.2, 172.4};
Int_t NumOfPtBins_part = 40;
double ptBins_part[] = {0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0, 1.1, 1.2, 1.4, 1.6, 1.8, 2.0, 2.2, 2.4, 3.2, 4.0, 4.8, 5.6, 6.4, 7.2, 9.6, 12.0, 14.4, 19.2, 24.0, 28.8, 35.2, 41.6, 48.0, 60.8, 73.6, 86.4, 103.6, 120.8, 138., 155.2, 172.4, 189.6, 206.8};
TTree * t = new TTree("t","whatever");
JetData data(t,1);
TH1D * hPartPt_minbias = new TH1D("hPartPt_minbias","hPartPt_minbias",NumOfPtBins_part-1,&ptBins_part[0]);
TH1D * hPartPt_minbias_trkCorr = new TH1D("hPartPt_minbias_trkCorr","hPartPt_minbias_trkCorr",NumOfPtBins_part-1,&ptBins_part[0]);
TH1D * hPartPt_minbias_trkCorr_trigCorr = new TH1D("hPartPt_minbias_trkCorr_trigCorr","hPartPt_minbias_trkCorr_trigCorr",NumOfPtBins_part-1,&ptBins_part[0]);
TH1D * hPartPt_minbias_trkCorr_trigCorr_smallerBins = new TH1D("hPartPt_minbias_trkCorr_trigCorr_smallerBins","hPartPt_minbias_trkCorr_trigCorr_smallerBins",ptBins_smaller.size()-1,&ptBins_smaller[0]);
TH1D * hPartPt_0_14_minbias = new TH1D("hPartPt_0_14_minbias","hPartPt_0_14_minbias",NumOfPtBins_part-1,&ptBins_part[0]);
TH1D * hPartPt_0_14_minbias_trkCorr = new TH1D("hPartPt_0_14_minbias_trkCorr","hPartPt_0_14_minbias_trkCorr",NumOfPtBins_part-1,&ptBins_part[0]);
TH1D * hPartPt_0_14_minbias_trkCorr_trigCorr = new TH1D("hPartPt_0_14_minbias_trkCorr_trigCorr","hPartPt_0_14_minbias_trkCorr_trigCorr",NumOfPtBins_part-1,&ptBins_part[0]);
TH1D * hNumEv_minbias = new TH1D("hNumEv_minbias","hNumEv_minbias",1,0.,1.);
TH1D * hNumEv_minbias_trigCorr = new TH1D("hNumEv_minbias_trigCorr","hNumEv_minbias_trigCorr",1,0.,1.);
TH1D * hNumEv_0_14_minbias = new TH1D("hNumEv_0_14_minbias","hNumEv_0_14_minbias",1,0.,1.);
TH1D * hNumEv_0_14_minbias_trigCorr = new TH1D("hNumEv_0_14_minbias_trigCorr","hNumEv_0_14_minbias_trigCorr",1,0.,1.);
TH1D * hNumEv_14_22_minbias = new TH1D("hNumEv_14_22_minbias","hNumEv_14_22_minbias",1,0.,1.);
TH1D * hNumEv_14_22_minbias_trigCorr = new TH1D("hNumEv_14_22_minbias_trigCorr","hNumEv_14_22_minbias_trigCorr",1,0.,1.);
TH1D * hNumEv_22_32_minbias = new TH1D("hNumEv_22_32_minbias","hNumEv_22_32_minbias",1,0.,1.);
TH1D * hNumEv_22_32_minbias_trigCorr = new TH1D("hNumEv_22_32_minbias_trigCorr","hNumEv_22_32_minbias_trigCorr",1,0.,1.);
TH1D * hNumEv_32_X_minbias = new TH1D("hNumEv_32_X_minbias","hNumEv_32_X_minbias",1,0.,1.);
TH1D * hRecMultiplicity = new TH1D("hRecMultiplicity","hRecMultiplicity",300.,0.,300.);
hPartPt_minbias->Sumw2();
hPartPt_minbias_trkCorr->Sumw2();
hPartPt_minbias_trkCorr_trigCorr->Sumw2();
hPartPt_minbias_trkCorr_trigCorr_smallerBins->Sumw2();
hPartPt_0_14_minbias->Sumw2();
hPartPt_0_14_minbias_trkCorr->Sumw2();
hPartPt_0_14_minbias_trkCorr_trigCorr->Sumw2();
hRecMultiplicity->Sumw2();
int numev_minbias = 0;
float numev_minbias_trigCorr = 0;
int numev_0_14_minbias = 0;
float numev_0_14_minbias_trigCorr = 0;
int numev_14_22_minbias = 0;
float numev_14_22_minbias_trigCorr = 0;
int numev_22_32_minbias = 0;
float numev_22_32_minbias_trigCorr = 0;
int numev_32_X_minbias = 0;
fillHistograms(c_minbias_1, data, numev_minbias, numev_minbias_trigCorr, numev_0_14_minbias, numev_0_14_minbias_trigCorr, numev_14_22_minbias, numev_14_22_minbias_trigCorr, numev_22_32_minbias, numev_22_32_minbias_trigCorr, numev_32_X_minbias, hPartPt_minbias, hPartPt_minbias_trkCorr, hPartPt_minbias_trkCorr_trigCorr, hPartPt_minbias_trkCorr_trigCorr_smallerBins, hPartPt_0_14_minbias, hPartPt_0_14_minbias_trkCorr, hPartPt_0_14_minbias_trkCorr_trigCorr,hRecMultiplicity);
// fillHistograms(c_minbias_2, data, numev_minbias, numev_minbias_trigCorr, numev_0_14_minbias, numev_0_14_minbias_trigCorr, numev_14_22_minbias, numev_14_22_minbias_trigCorr, numev_22_32_minbias, numev_22_32_minbias_trigCorr, numev_32_X_minbias, hPartPt_minbias, hPartPt_minbias_trkCorr, hPartPt_minbias_trkCorr_trigCorr, hPartPt_minbias_trkCorr_trigCorr_smallerBins, hPartPt_0_14_minbias, hPartPt_0_14_minbias_trkCorr, hPartPt_0_14_minbias_trkCorr_trigCorr,hRecMultiplicity);
hNumEv_minbias->SetBinContent(1,numev_minbias);
hNumEv_minbias_trigCorr->SetBinContent(1,numev_minbias_trigCorr);
hNumEv_0_14_minbias->SetBinContent(1,numev_0_14_minbias);
hNumEv_0_14_minbias_trigCorr->SetBinContent(1,numev_0_14_minbias_trigCorr);
hNumEv_14_22_minbias->SetBinContent(1,numev_14_22_minbias);
hNumEv_14_22_minbias_trigCorr->SetBinContent(1,numev_14_22_minbias_trigCorr);
hNumEv_22_32_minbias->SetBinContent(1,numev_22_32_minbias);
hNumEv_22_32_minbias_trigCorr->SetBinContent(1,numev_22_32_minbias_trigCorr);
hNumEv_32_X_minbias->SetBinContent(1,numev_32_X_minbias);
normalizeByBinWidth(hPartPt_minbias);
normalizeByBinWidth(hPartPt_minbias_trkCorr);
normalizeByBinWidth(hPartPt_minbias_trkCorr_trigCorr);
normalizeByBinWidth(hPartPt_minbias_trkCorr_trigCorr_smallerBins);
normalizeByBinWidth(hPartPt_0_14_minbias);
normalizeByBinWidth(hPartPt_0_14_minbias_trkCorr);
normalizeByBinWidth(hPartPt_0_14_minbias_trkCorr_trigCorr);
f_output->cd();
hPartPt_minbias->Write();
hPartPt_minbias_trkCorr->Write();
hPartPt_minbias_trkCorr_trigCorr->Write();
hPartPt_minbias_trkCorr_trigCorr_smallerBins->Write();
hPartPt_0_14_minbias->Write();
hPartPt_0_14_minbias_trkCorr->Write();
hPartPt_0_14_minbias_trkCorr_trigCorr->Write();
hNumEv_minbias->Write();
hNumEv_minbias_trigCorr->Write();
hNumEv_0_14_minbias->Write();
hNumEv_0_14_minbias_trigCorr->Write();
hNumEv_14_22_minbias->Write();
hNumEv_14_22_minbias_trigCorr->Write();
hNumEv_22_32_minbias->Write();
hNumEv_22_32_minbias_trigCorr->Write();
hNumEv_32_X_minbias->Write();
hRecMultiplicity->Write();
f_output->Close();
}
