EXP1803

#include "TSystem.h"

#include "TFile.h"

#include "TTree.h"

#include "TChain.h"

#include "TCanvas.h"

#include "TBox.h"

#include "TCut.h"

#include "TCutG.h"

using std::cout;

using std::endl;

void showBananas2(const Long64_t drawEntries = 100000) {

//        TFile *fr = new TFile("~/data/exp1804/h5_11_00.root");

        TChain *tr = new TChain("AnalysisxTree");

//        tr->Add("~/data/exp1804/h5_14_0?.root");

        tr->Add("~/data/exp1804/be10_0?_??.root");

        TChain *trCal = new TChain("cal");

//        trCal->Add("~/data/exp1804/h5_11_0?_calib.root");

//        trCal->Add("~/data/exp1804/h5_14_0?_calib.root");

        trCal->Add("~/data/exp1804/h5_14_0?_calib.root");

        tr->AddFriend(trCal);

        TFile *fcal = new TFile("~/data/exp1804/calib/si_20_03_calib.root");

        TTree *tAlphaCal = (TTree*)fcal->Get("cal");

        tAlphaCal->SetMarkerColor(kRed);

        tAlphaCal->SetMarkerStyle(20);

        tAlphaCal->SetMarkerSize(.6);

        TFile *fc = new TFile("cutsIdentification.root", "READ");

        TCutG *cUS = (TCutG*)fc->Get("c7UpperShadow");

        TCutG *cLS = (TCutG*)fc->Get("c7LowerShadow");

        cLS->SetLineColor(kMagenta);

        TCutG *cMA = (TCutG*)fc->Get("c7MainAlpha");

        TCutG *cA = (TCutG*)fc->Get("c7All");

        TCutG *cBL = (TCutG*)fc->Get("cBeamLeft");

        TCutG *cBR = (TCutG*)fc->Get("cBeamRight");

//        const Int_t drawEntries = 6000000;

//        const Int_t drawEntries = tr->GetEntries();

        tr->GetListOfFiles()->Print();

        trCal->GetListOfFiles()->Print();

        cout << tr->GetEntries() << " events." << endl;

        cout << trCal->GetEntries() << " calibrated events." << endl;

        TString varName;

        TString condition;

        cout << drawEntries << " entries in input chain will be processed." << endl;

        const Int_t firstThinStrip = 5;

/////////////////////////////////////////////////////////////////////

//        c1

        /*TCanvas *c1 = new TCanvas("c1", "time issues", 1600, 800);

        c1->Divide(4,2);

        for (Int_t i = 0; i < 8; i++) {

                c1->cd(i+1);

//                varName.Form("SQ20E[%d]:SQLYEsum", i+firstThinStrip);

                varName.Form("SQ20E[%d]:SQ20E[%d]+SQLYEsum", i+firstThinStrip, i+firstThinStrip);

                condition.Form("SQ20E[%d]>0 && SQ20E[%d]<14 && SQLYEsum>0. && SQLYEsum<120. && SQLYE[0]<1. && SQLYE[15]<1.", i+firstThinStrip, i+firstThinStrip);

                varName.Form("SQX_L[%d]:tSQX_L[%d]", i+firstThinStrip, i+firstThinStrip);

                condition.Form("tSQX_L[%d]>0 && SQX_L[%d]>0", i+firstThinStrip, i+firstThinStrip);

                tr->SetMarkerColor(kBlack);

                tr->SetMarkerStyle(20);

                tr->SetMarkerSize(0.1);

                tr->Draw(varName, condition, "", drawEntries);

                c1->Update();

        }

        return;*/

/////////////////////////////////////////////////////////////////////

//        c2

        TCanvas *c2 = new TCanvas("c2", "1 mm Y: 8-15 zoomed", 1600, 800);

        c2->Divide(4,2);

        for (Int_t i = 0; i < 8; i++) {

                c2->cd(i+1);

                varName.Form("SQ20E[%d]:SQLYEsum", i+firstThinStrip);

                varName.Form("SQ20E[%d]:SQ20E[%d]+SQLYEsum", i+firstThinStrip, i+firstThinStrip);

                condition.Form("SQ20E[%d]>1.2 && SQ20E[%d]<5 "

                                "&& SQLYEsum>1.1 && SQLYEsum<30 "

                                "&& SQLYE[0]<1. && SQLYE[15]<1."

                                "&& SQLYmult==1 && SQLXmult==1"

                                "&& trigger==3",

                                i+firstThinStrip, i+firstThinStrip);

                tr->SetMarkerColor(kBlack);

                tr->Draw(varName, condition, "col", drawEntries);

//                tr->Draw(varName, condition, "", drawEntries);

                tAlphaCal->Draw(varName, "", "same");

                c2->Update();

        }//for */

//        return;

/////////////////////////////////////////////////////////////////////

//        c3

        TCanvas *c3 = new TCanvas("c3", "dE-E corrected", 1600, 800);

        c3->Divide(4,2);

        for (Int_t i = 0; i < 8; i++) {

                c3->cd(i+1);

//                varName.Form("SQ20Ecorr[%d]:SQLYEsum", i+firstThinStrip);

                varName.Form("SQ20Ecorr[%d]:SQ20Ecorr[%d]+SQLYEsum", i+firstThinStrip, i+firstThinStrip);

                varName.Form("SQ20Ecorr[%d]:SQ20Ecorr[%d]+SQLXEsum", i+firstThinStrip, i+firstThinStrip);

                condition.Form("SQ20Ecorr[%d]>1.2 && SQ20Ecorr[%d]<5 "

                                "&& SQLYEsum>1.1 && SQLYEsum<30 "

                                "&& SQLXEsum>1.1 && SQLXEsum<30 "

                                "&& SQLYE[0]<1. && SQLYE[15]<1. "

                                "&& SQLYmult==1 && SQLXmult==1"

                                "&& trigger==3",

                                i+firstThinStrip, i+firstThinStrip);

                tr->SetMarkerColor(kBlack);

                tr->Draw(varName, condition, "col", drawEntries);

//                tr->Draw(varName, condition, "", drawEntries);

                tAlphaCal->Draw(varName, "", "same");

                c3->Update();

        }

        return;

/////////////////////////////////////////////////////////////////////

//        c4

        TCanvas *c4 = new TCanvas("c5", "MWPC: different dE-E regions", 1000, 1000);

        c4->Divide(2,2);

        c4->cd(1);

        tr->Draw("y1p:x1p>>hXY1dr(32, -20., 20., 32, -20., 20.)",

                        "TOF<185 && TOF>165 && dEbeam>600 && dEbeam<1150 && trigger==3 && y1p>-50 && x1p>-50",

                        "col", 1000000);

        c4->Update();

        tr->SetMarkerStyle(20);

        tr->SetMarkerSize(0.4);

        condition.Form("TOF<185 && TOF>165 && dEbeam>600 && dEbeam<1150 && trigger==3  && yt>-50 && xt>-50 && c7MainAlpha");

        tr->SetMarkerColor(kGreen);

        tr->SetMarkerSize(1.);

        tr->Draw(varName, condition, "same", drawEntries);

        c4->Update();

        condition.Form("TOF<185 && TOF>165 && dEbeam>600 && dEbeam<1150 && trigger==3  && yt>-50 && xt>-50 && c7UpperShadow");

        tr->SetMarkerColor(kRed);

        tr->SetMarkerStyle(3);

        tr->SetMarkerSize(1.7);

        tr->Draw(varName, condition, "same", drawEntries);

        c4->Update();

        condition.Form("TOF<185 && TOF>165 && dEbeam>600 && dEbeam<1150 && trigger==3  && yt>-50 && xt>-50 && c7LowerShadow");

        tr->SetMarkerColor(kMagenta);

        tr->SetMarkerStyle(20);

        tr->SetMarkerSize(1.);

        tr->Draw(varName, condition, "same", drawEntries);

        c4->Update();

        c4->cd(2);

        tr->Draw("y2p:x2p>>hXY2dr(32, -20., 20., 32, -20., 20.)",

                        "TOF<185 && TOF>165 && dEbeam>600 && dEbeam<1150 && trigger==3 && y2p>-50 && x2p>-50",

                        "col", 1000000);

        c4->Update();

        c4->cd(4);

        condition.Form("TOF<185 && TOF>165 && dEbeam>600 && dEbeam<1150 && trigger==3  && yt>-50 && xt>-50");

        varName.Form("yt:xt>>hXYtDR(32, -40., 40., 32, -40., 40.)");

        tr->Draw(varName, condition, "col", 100000);

        TBox *targetBoarder = new TBox(-25., 20., 25., -20.);

        targetBoarder->SetFillStyle(0);

        targetBoarder->SetLineWidth(3);

        targetBoarder->SetLineColor(kRed);

        targetBoarder->Draw();

        c4->Update();

        tr->SetMarkerStyle(20);

        tr->SetMarkerSize(0.4);

        varName.Form("yt:xt");

        condition.Form("TOF<185 && TOF>165 && dEbeam>600 && dEbeam<1150 && trigger==3  && yt>-50 && xt>-50 && c7MainAlpha");

        tr->SetMarkerColor(kGreen);

        tr->SetMarkerSize(1.);

        tr->Draw(varName, condition, "same", drawEntries);

        c4->Update();

        condition.Form("TOF<185 && TOF>165 && dEbeam>600 && dEbeam<1150 && trigger==3  && yt>-50 && xt>-50 && c7UpperShadow");

        tr->SetMarkerColor(kRed);

        tr->SetMarkerStyle(3);

        tr->SetMarkerSize(1.7);

        tr->Draw(varName, condition, "same", drawEntries);

        c4->Update();

        condition.Form("TOF<185 && TOF>165 && dEbeam>600 && dEbeam<1150 && trigger==3  && yt>-50 && xt>-50 && c7LowerShadow");

        tr->SetMarkerColor(kMagenta);

        tr->SetMarkerStyle(20);

        tr->SetMarkerSize(1.);

        tr->Draw(varName, condition, "same", drawEntries);

        c4->Update();

/////////////////////////////////////////////////////////////////////

//        c5

        TCanvas *c5 = new TCanvas("c4", "MWPC", 1000, 1000);

        c5->Divide(2,2);

        c5->cd(4);

        condition.Form("TOF<185 && TOF>165 && trigger==3  && yt>-50 && xt>-50");

        varName.Form("yt:xt>>hXYt(32, -40., 40., 32, -40., 40.)");

        tr->Draw(varName, condition, "col", drawEntries);

        targetBoarder->Draw();

        c5->Update();

        c5->cd(1);

        tr->Draw("y1p:x1p>>hXY1(32, -20., 20., 32, -20., 20.)",

                        "TOF<185 && TOF>165 && dEbeam>600 && dEbeam<1150 && trigger==3 && y1p>-50 && x1p>-50", "col", 1000000);

        cBL->Draw("same");

        cBR->Draw("same");

        c5->Update();

        c5->cd(2);

        tr->Draw("y2p:x2p>>hXY2(32, -20., 20., 32, -20., 20.)",

                        "TOF<185 && TOF>165 && dEbeam>600 && dEbeam<1150 && trigger==3 && y2p>-50 && x2p>-50", "col", 1000000);

        c5->Update();

        c5->cd(3);

        tr->Draw("dEbeam:TOF>>IDplot(200,,,200,,)", "TOF<200 && TOF>100 && dEbeam < 2500 && trigger==1", "col", 1000000);

        TBox *beamBox = new TBox(165, 1150, 185, 600);

        beamBox->SetFillStyle(0);

        beamBox->SetLineWidth(3);

        beamBox->SetLineColor(kRed);

        beamBox->Draw();

        c5->Update();

/////////////////////////////////////////////////////////////////////

//        c6

        TCanvas *c6 = new TCanvas("c6", "1 mm Y: 8-15 - target locus", 1600, 800);

        c6->Divide(4,2);

        TString cTarget("xt>-10. && xt<10. && yt>-15. && yt<15.");

        c6->cd(1);

        tr->Draw("SQ20E[0]:SQ20E[0]+SQLYEsum>>(200,0,5,200,0,5)", "", "", 100);

        for (Int_t i = 0; i < 8; i++) {

                c6->cd(i+1);

                varName.Form("SQ20E[%d]:SQ20E[%d]+SQLYEsum", i+firstThinStrip, i+firstThinStrip);

                condition.Form("SQ20E[%d]>1.2 && SQ20E[%d]<5.5 && SQLYEsum>1.1 && SQLYEsum<30 && SQLYE[0]<1. && SQLYE[15]<1.",

                                i+firstThinStrip, i+firstThinStrip);

                tr->Draw(varName, condition, "col", drawEntries);

                c6->Update();

                condition.Form("SQ20E[%d]>1.2 && SQ20E[%d]<5.5 && SQLYEsum>1.1 && SQLYEsum<30 && SQLYE[0]<1. && SQLYE[15]<1. && %s",

                                i+firstThinStrip, i+firstThinStrip, cTarget.Data());

                tr->SetMarkerColor(kRed);

                tr->SetMarkerStyle(20);

                tr->SetMarkerSize(.5);

                tr->Draw(varName, condition, "same", drawEntries);

                c6->Update();

        }

/////////////////////////////////////////////////////////////////////

//        c7

        TCanvas *c7 = new TCanvas("c7", "1 mm Y: 8-15", 1600, 800);

        c7->Divide(4,2);

        TString cStripsLow("SQLXE[0]<1. && SQLXE[1]<1. && SQLXE[2]<1. && SQLXE[3]<1. && SQLXE[4]<1. && SQLXE[5]<1. && SQLXE[6]<1. && SQLXE[7]<1. && SQLXE[8]<1. && SQLXE[9]<1. && SQLXE[10]<1.");

        TString cStripsHigh("SQLXE[31]<1. && SQLXE[30]<1. && SQLXE[29]<1. && SQLXE[28]<1. && SQLXE[27]<1. && SQLXE[26]<1. && SQLXE[25]<1. && SQLXE[24]<1. && SQLXE[23]<1. && SQLXE[22]<1. && SQLXE[21]<1. && SQLXE[20]<1.");

        for (Int_t i = 0; i < 8; i++) {

                c7->cd(i+1);

                varName.Form("SQ20E[%d]:SQ20E[%d]+SQLXEsum", i+firstThinStrip, i+firstThinStrip);

                condition.Form("SQ20E[%d]>1.2 && SQ20E[%d]<5.5 "

                                "&& SQLXEsum>1.1 && SQLXEsum<30 "

                                "&& %s && %s"

                                , i+firstThinStrip, i+firstThinStrip, cStripsLow.Data(), cStripsHigh.Data());

                tr->Draw(varName, condition, "col", drawEntries);

                c7->Update();

        }

        return;

}