Commit b6d0b970 authored by Vratislav Chudoba's avatar Vratislav Chudoba

Interval 60-75 done.

parent 3a99821c
//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};
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;
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.");
}
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment