diff --git a/macros/BeCorrPRC/ES1_6075.cxx b/macros/BeCorrPRC/ES1_6075.cxx index a3aef697a01110ef8605de74ae74510f2760ca62..518e6085d010c524ee1bdd4eccaf2e0612c11370 100644 --- a/macros/BeCorrPRC/ES1_6075.cxx +++ b/macros/BeCorrPRC/ES1_6075.cxx @@ -1,210 +1,124 @@ -//void ES1_6075() +void ES1_6075() { - //open chain with Be events only and draw E_Be - - //order: al0, noal0, al180, noal180 TStopwatch sw; sw.Start(); - const Int_t kMinAngle = 60; - const Int_t kMaxAngle = 75; - - const Int_t lowfile = 0; - const Int_t upfile = 40; - - const Int_t lowSimFile[6] = {0, 0, 0, 0, 0, 0}; - const Int_t upSimFile[6] = {14, 14, 14, 14, 14, 14}; - - - const Bool_t shapeFit = 0; - - TCut cAngles = "fBeThetaCM1>60*TMath::DegToRad() && fBeThetaCM1<75*TMath::DegToRad()"; - TCut crAngles = "f6BeThetaCM1>60*TMath::DegToRad() && f6BeThetaCM1<75*TMath::DegToRad()"; - - //simulations - //aligned, 0 degrees - TString chainAl0name = "../../../be/rootdata/correlations/v5_6075/Sim_mix_br_al_0_50-85_"; - TString inputTreeAl0name = "../../../be/rootdata/correlations/v5_6075/Sim_mix_br_al_0_50-85_000.root"; - //isotropic, 0 degrees - TString chainNoAl0name = "../../../be/rootdata/correlations/v5_6075/Sim_mix_br_noal_0_50-85_"; - TString inputTreeNoAl0name = "../../../be/rootdata/correlations/v5_6075/Sim_mix_br_noal_0_50-85_000.root"; - //aligned, 180 degrees - TString chainAl180name = "../../../be/rootdata/correlations/v5_6075/Sim_mix_br_al_180_50-85_"; - TString inputTreeAl180name = "../../../be/rootdata/correlations/v5_6075/Sim_mix_br_al_180_50-85_000.root"; - //isotropic, 180 degrees - TString chainNoAl180name = "../../../be/rootdata/correlations/v5_6075/Sim_mix_br_noal_180_50-85_"; - TString inputTreeNoAl180name = "../../../be/rootdata/correlations/v5_6075/Sim_mix_br_noal_180_50-85_000.root"; - //isotropic, 90 degrees - TString chainAl90name = "../../../be/rootdata/correlations/v5_6075/Sim_mix_br_al_90_50-85_"; - TString inputTreeAl90name = "../../../be/rootdata/correlations/v5_6075/Sim_mix_br_al_90_50-85_000.root"; - //isotropic, 90 degrees - TString chainNoAl90name = "../../../be/rootdata/correlations/v5_6075/Sim_mix_br_noal_90_50-85_"; - TString inputTreeNoAl90name = "../../../be/rootdata/correlations/v5_6075/Sim_mix_br_noal_90_50-85_000.root"; - -// TChain *che = BeWork::OpenChain("/data2/be/rootdata/correlations/Be.", -// lowfile, upfile, "beonly"); //original file -// Info("statesRatioFitting.cxx", "%d events in chain \"%s\" containing experimental data", -// che->GetEntries(), che->GetName()); - - //open tree with MC simulations - gROOT->ProcessLine(".x ../corr_templates/sfAngInt_opensimtree.cxx"); - - TChain *che = BeWork::OpenChain("../../../be/rootdata/correlations/v5_exp/Be.", - lowfile, upfile, "beonly"); //original file - Info("statesRatioFitting.cxx", "%d events in chain \"%s\" containing experimental data", - che->GetEntries(), che->GetName()); - - gROOT->ProcessLine(".x ../corr_templates/sfAngInt_cuts.cxx"); - - //auxiliary strings - TString drawCommand; - TString hsName; - TString heName; - TString hiName; - - sw.Print(); - sw.Continue(); - - const Bool_t savePictures = 0; - TString ppath = "pictures/correlations/v5_6075/"; - TString configuration = "6075"; -// TString pictFormat = ".ps"; -// gStyle->SetHistLineWidth(1); - TString pictFormat = ".eps"; -// TString pictFormat = ".gif"; + BeAnalysis ana; - const Bool_t kVerbose = 1; - const Bool_t kAutoRange = 1; - // const Float_t kRangeProportion = 1.2; - const Float_t kRangeProportion = 1.1; - // const Float_t kRangeProportionIn = 1.1; - const Float_t kRangeProportionIn = 2.0; + ana.SetCMAngularRange(60, 75); -// const Bool_t kChains[6] = {0,0,0,1,0,1}; - const Bool_t kChains[6] = {1,1,1,1,1,1}; + // ana.SetFigures("figures/", ".eps"); - const Bool_t spectra = 0; + // ana.SetChainsToDraw(0, 0, 0, 0, 1, 1); + ana.SetChainsToDraw(); + Bool_t epsilonT[6] = {0,0,1,0,0}; + Bool_t cosThetaT[6] = {0,0,1,0,0}; - const Int_t nointervals = 5; + Bool_t epsilonY[6] = {0,0,1,0,0}; + Bool_t cosThetaY[6] = {0,0,1,0,0}; - const Bool_t cosThetaT[nointervals] = {0,0,0,0,0}; -// const Bool_t cosThetaT[nointervals] = {1,1,1,1,1}; + Bool_t thetaAT[6] = {0,0,1,0,0}; - const Bool_t epsilonT[nointervals] = {1,1,1,1,1}; -// const Bool_t epsilonT[nointervals] = {1,1,1,1,1}; - const Bool_t epsilonY[nointervals] = {0,0,0,0,0}; -// const Bool_t epsilonY[nointervals] = {1,1,1,1,1}; - const Bool_t cosThetaY[nointervals] = {0,1,0,1,0}; -// const Bool_t cosThetaY[nointervals] = {1,1,1,1,1}; + ana.SetExpChain("../../../be/rootdata/correlations/v5_exp/Be.", 0, 40); - const Bool_t thetaAT[nointervals] = {0,0,0,0,0}; -// const Bool_t thetaAT[nointervals] = {1,1,1,1,1}; - - const Bool_t thetaP1T[nointervals] = {0,0,1,0,0}; -// const Bool_t thetaP1T[nointervals] = {1,1,1,1,1}; - - TCut sRatioAl0 = "sRatio>0.051 && sRatio<0.101"; -// TCut sRatioAl0 = "sRatio>0.060 && sRatio<0.110"; - TCut sRatioNoAl0 = "sRatio>0.050 && sRatio<0.100"; - TCut sRatioAl180 = "sRatio>0.047 && sRatio<0.097"; - TCut sRatioNoAl180 = "sRatio>0.048 && sRatio<0.098"; - TCut sRatioAl90 = "sRatio>0.051 && sRatio<0.101"; - TCut sRatioNoAl90 = "sRatio>0.051 && sRatio<0.101"; - TCut sRatio[6] = {sRatioAl0, sRatioNoAl0, sRatioAl180, sRatioNoAl180, sRatioAl90, sRatioNoAl90}; - - //integral fit const Long64_t eMaxEvents = 5000000; - Long64_t eEventsAl0 = eMaxEvents; - Long64_t eEventsNoAl0 = eMaxEvents; - Long64_t eEventsAl180 = eMaxEvents; - Long64_t eEventsNoAl180 = eMaxEvents; - Long64_t eEventsAl90 = eMaxEvents; - Long64_t eEventsNoAl90 = eMaxEvents; + Long64_t eEvents[6] = {eMaxEvents, eMaxEvents, eMaxEvents, eMaxEvents, eMaxEvents, eMaxEvents}; + ana.SetNoExpEvents(eEvents); + + const Long64_t sMaxEvents = 3000000; + Long64_t sEvents[6] = {sMaxEvents, 2700000, sMaxEvents, 2700000, 2690000, 2710000}; + ana.SetNoSimEvents(sEvents); + + Int_t lowSimFile[6] = {0, 0, 0, 0, 0, 0}; + Int_t upSimFile[6] = {14, 14, 14, 14, 14, 14}; + ana.SetNoSimFiles(lowSimFile, upSimFile); + ana.SetSimChains("../../../be/rootdata/correlations/v5_6075/Sim_mix_br_al_0_50-85_", + "../../../be/rootdata/correlations/v5_6075/Sim_mix_br_noal_0_50-85_", + "../../../be/rootdata/correlations/v5_6075/Sim_mix_br_al_180_50-85_", + "../../../be/rootdata/correlations/v5_6075/Sim_mix_br_noal_180_50-85_", + "../../../be/rootdata/correlations/v5_6075/Sim_mix_br_al_90_50-85_", + "../../../be/rootdata/correlations/v5_6075/Sim_mix_br_noal_90_50-85_"); + + ana.SetInputChains("../../../be/rootdata/correlations/v5_6075/Sim_mix_br_al_0_50-85_000.root", + "../../../be/rootdata/correlations/v5_6075/Sim_mix_br_noal_0_50-85_000.root", + "../../../be/rootdata/correlations/v5_6075/Sim_mix_br_al_180_50-85_000.root", + "../../../be/rootdata/correlations/v5_6075/Sim_mix_br_noal_180_50-85_000.root", + "../../../be/rootdata/correlations/v5_6075/Sim_mix_br_al_90_50-85_000.root", + "../../../be/rootdata/correlations/v5_6075/Sim_mix_br_noal_90_50-85_000.root"); + + + ana.OpenExpChain(); + ana.OpenSimChains(); + + ana.SetRatiosGStoEX("sRatio>0.051 && sRatio<0.101", + "sRatio>0.050 && sRatio<0.100", + "sRatio>0.047 && sRatio<0.097", + "sRatio>0.048 && sRatio<0.098", + "sRatio>0.051 && sRatio<0.101", + "sRatio>0.051 && sRatio<0.101"); + ana.SetCuts(); Long64_t eEventsECuts[5][6] = {3950000, 3950000, 3800000, 3850000, 3950000, 3900000, eMaxEvents, eMaxEvents, 4100000, eMaxEvents, eMaxEvents, eMaxEvents, 3050000, 3150000, 3200000, 3200000, 3250000, 3350000, 3850000, 3950000, 3900000, 3900000, 4000000, 4000000, 4100000, 4100000, 4000000, 4100000, eMaxEvents, eMaxEvents}; - - const Long64_t eEvents[6] = {eEventsAl0, eEventsNoAl0, eEventsAl180, eEventsNoAl180, eEventsAl90, eEventsNoAl90}; - const TChain *chs[6] = {chsAl0, chsNoAl0, chsAl180, chsNoAl180, chsAl90, chsNoAl90}; - const TTree *ti[6] = {tiAl0, tiNoAl0, tiAl180, tiNoAl180, tiAl90, tiNoAl90}; + ana.ExpEventsECuts(eEventsECuts); const Long64_t sMaxEvents = 3000000; - Long64_t sEventsAl0 = sMaxEvents; - Long64_t sEventsNoAl0 = 2700000; - Long64_t sEventsAl180 = sMaxEvents; - Long64_t sEventsNoAl180 = 2700000; - Long64_t sEventsAl90 = 2690000; - Long64_t sEventsNoAl90 = 2710000; - const Long64_t sEvents[6] = {sEventsAl0, sEventsNoAl0, sEventsAl180, sEventsNoAl180, sEventsAl90, sEventsNoAl90}; Long64_t sEventsECuts[5][6] = {sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents, - sMaxEvents, 2630000, sMaxEvents, sMaxEvents, 2660000, 2660000, - sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents, - sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents, - sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents}; - - TString canvasTitle; - TString canvasName; - - const Double_t sRange[6] = {2000, 2000, 2000, 2000, 2000, 2000}; - if (spectra) { - gROOT->ProcessLine(".x ../corr_templates/sfAngInt_spectra_v5.cxx"); - } - - // "T" system - const Double_t eTRange[5][6] = {300, 300, 300, 300, 300, 300, - 380, 380, 380, 380, 380, 380, - 300, 300, 300, 300, 300, 300, - 500, 500, 500, 500, 500, 500, - 500, 500, 500, 500, 500, 500}; - gROOT->ProcessLine(".x ../corr_templates/sfAngInt_epsilonT.cxx"); - -// return; - - const Double_t coskTRange[5][6] = {350, 350, 350, 350, 350, 350, - 400, 400, 400, 400, 400, 400, - 350, 350, 350, 350, 350, 350, - 700, 700, 700, 700, 700, 700, - 900, 900, 900, 900, 900, 900}; - gROOT->ProcessLine(".x ../corr_templates/sfAngInt_cosThetaTk.cxx"); - - // "Y" system - const Double_t eYRange[5][6] = {400, 400, 400, 400, 400, 400, - 450, 450, 450, 450, 450, 450, - 400, 400, 400, 400, 400, 400, - 800, 800, 800, 800, 800, 800, - 1000, 1000, 1000, 1000, 1000, 1000}; - gROOT->ProcessLine(".x ../corr_templates/sfAngInt_epsilonY.cxx"); - - const Double_t coskYRange[5][6] = {250, 250, 250, 250, 250, 250, - 350, 350, 350, 350, 350, 350, - 300, 300, 300, 300, 300, 300, - 450, 450, 450, 450, 450, 450, - 500, 500, 500, 500, 500, 500}; - gROOT->ProcessLine(".x ../corr_templates/sfAngInt_cosThetaYk.cxx"); - - - - const Double_t tP1Range[5][6] = {90, 90, 90, 90, 90, 90, - 140, 140, 140, 140, 140, 140, - 130, 130, 130, 130, 130, 130, - 190, 190, 190, 190, 190, 190, - 200, 200, 200, 200, 200, 200}; - gROOT->ProcessLine(".x ../corr_templates/sfAngInt_thetaP1T.cxx"); - - const Double_t tARange[5][6] = {90, 90, 90, 90, 90, 90, - 140, 140, 140, 140, 140, 140, - 160, 120, 160, 120, 160, 120, - 230, 180, 230, 180, 230, 180, - 210, 210, 210, 210, 210, 210}; - gROOT->ProcessLine(".x ../corr_templates/sfAngInt_thetaAT.cxx"); - - Info("correlations_v5_6075.cxx", "Finished."); + sMaxEvents, 2630000, sMaxEvents, sMaxEvents, 2660000, 2660000, + sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents, + sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents, + sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents}; + ana.SimEventsECuts(sEventsECuts); + + ana.SetRangeProportion(1.1, 1); + + /////////////////////////////////////////////////////////////////////////////////// + + ana.SetSpectra(); + ana.Spectra(); + + ana.SetEpsilonTintervals(epsilonT); + // ana.SetEpsilonTintervals(); + ana.EpsilonT(); + + sw.Print(); + sw.Continue(); + + ana.SetCosThetaTkIntervals(cosThetaT); + // ana.CosThetaTk(); + + sw.Print(); + sw.Continue(); + + ana.SetEpsilonYintervals(epsilonY); + // ana.EpsilonY(); + + sw.Print(); + sw.Continue(); + + ana.SetCosThetaYkIntervals(cosThetaY); + // ana.CosThetaYk(); + + sw.Print(); + sw.Continue(); + + ana.SetThetaATintervals(thetaAT); + //// ana.SetThetaATintervals(); + // ana.ThetaAT(); + + sw.Print(); + sw.Continue(); + + Info("ES1_90120.cxx", "Finished."); + + return; + }