Back to home page

sPhenix code displayed by LXR

 
 

    


File indexing completed on 2025-12-16 09:18:01

0001 void draw_sicalo_corr(){
0002   //TFile *f = TFile::Open("dphi_distribution_em.root");
0003   //TFile *f = TFile::Open("dphi_distribution_jpsi.root");
0004   //TFile *f = TFile::Open("dphi_distribution_pythia.root");
0005   //TFile *f = TFile::Open("dphi_distribution_pythia.root");
0006   TFile *f = TFile::Open("dphi_distribution.root");
0007 
0008 //  TNtuple *ntp_sicalo = (TNtuple*)f->Get("ntp_sicalo");
0009   TNtuple *ntp_sicalo = (TNtuple*)f->Get("ntp_sicalotrk");
0010 
0011   TCanvas *c1 = new TCanvas("c1", "c1", 1000, 750);
0012   c1->Divide(3,2);
0013   c1->cd(1);
0014   ntp_sicalo->Draw("z_emc:z_pemc>>h(100,-150,150,100,-150,150)","e_emc>0.2&&pt>0.5","colz");
0015   c1->cd(2);
0016   ntp_sicalo->Draw("z_emc-z_pemc:z_pemc>>h2(100,-150,150, 100,-50,50)","e_emc>0.2&&pt>0.5","colz");
0017   c1->cd(4);
0018   ntp_sicalo->Draw("phi_emc:phi_pemc>>h3(100,-3.2,3.2,100,-3.2,3.2)","e_emc>0.2&&pt>0.5","colz");
0019   c1->cd(5);
0020   ntp_sicalo->Draw("phi_emc:phi_pemc>>h4(100,-3.2,3.2,100,-3.2,3.2)","e_emc>0.2&&pt>0.5&&abs(z_emc-z_pemc)<6","colz");
0021 
0022   c1->cd(3);
0023   ntp_sicalo->Draw("e_emc/sqrt(pt*pt+pz*pz):sqrt(pt*pt+pz*pz)>>h5(50,0,5,100,0,3)","abs(z_emc-z_pemc)<5","colz");
0024   c1->cd(6);
0025   ntp_sicalo->Draw("phi_emc-phi_pemc:pt>>h6(100,0,5,100,-2,2)","e_emc>0.2&&pt>0.&&abs(z_emc-z_pemc)<5","colz");
0026 
0027 
0028   TCanvas *c2 = new TCanvas("c2", "c2", 1000, 750);
0029   c2->Divide(2,2);
0030   c2->cd(1);
0031   ntp_sicalo->Draw("phi_emc-phi_pemc:pt>>h7(100,0,5,100,-2,2)","e_emc>0.2&&pt>0.&&abs(z_emc-z_pemc)<5&&abs(phi_pemc)<1.6&&c>0","colz");
0032   c2->cd(2);
0033   ntp_sicalo->Draw("phi_emc-phi_pemc:pt>>h8(100,0,5,100,-2,2)","e_emc>0.2&&pt>0.&&abs(z_emc-z_pemc)<5&&abs(phi_pemc)<1.6&&c<0","colz");
0034   c2->cd(3);
0035   ntp_sicalo->Draw("phi_emc-phi_pemc:pt>>h9(100,0,5,100,-2,2)","e_emc>0.2&&pt>0.&&abs(z_emc-z_pemc)<5&&abs(phi_pemc)>1.6&&c>0","colz");
0036   c2->cd(4);
0037   ntp_sicalo->Draw("phi_emc-phi_pemc:pt>>h10(100,0,5,100,-2,2)","e_emc>0.2&&pt>0.&&abs(z_emc-z_pemc)<5&&abs(phi_pemc)>1.6&&c<0","colz");
0038 
0039   TCanvas *c3 = new TCanvas("c3", "c3", 1000, 750);
0040   c3->Divide(2,2);
0041   c3->cd(1);
0042   ntp_sicalo->Draw("c*(phi_emc-phi_pemc):pt>>h11(100,0,5,100,-2,2)","e_emc>0.2&&pt>0.&&abs(z_emc-z_pemc)<5","colz");
0043   c3->cd(2);
0044   gPad->SetLogx();
0045   gPad->SetLogy();
0046   TH2* hh = (TH2*)gROOT->FindObject("h11");
0047   hh->Draw("colz");
0048 
0049   TF1* ff = new TF1("ff", "[0]*pow(x,[1])");
0050   //ff->SetParameters(0.18, -0.986);
0051   ff->SetParameters(0.18, -0.988);
0052   ff->Draw("same");
0053 
0054   c3->cd(3);
0055   ntp_sicalo->Draw("c*(phi_emc-phi_pemc)-(0.18*pow(pt,-0.986)):pt>>h12(100,0,5,100,-2,2)","e_emc>0.2&&pt>0.&&abs(z_emc-z_pemc)<5","colz");
0056 
0057   c3->cd(4);
0058 ntp_sicalo->Draw("c*(phi_emc-phi_pemc):pt>>h13(100,0,5,100,-2,2)","e_emc>0.2&&pt>0.&&abs(z_emc-z_pemc)<5&&abs(c*(phi_emc-phi_pemc)-(0.21*pow(pt,-0.986)))<0.2","colz");
0059 
0060 //
0061 //--  ntp_sicalo->Draw("phi_emc-phi_pemc:phi_pemc>>h(100,-3.2,3.2,100,-1,1)","e_emc>0.2&&pt>0.5&&abs(z_emc-z_pemc)<5","colz")
0062 //--  ntp_sicalo->Draw("phi_emc-phi_pemc:phi_pemc>>h(100,-3.2,3.2,100,-1,1)","e_emc>0.2&&pt>0.5&&pt<2&&abs(z_emc-z_pemc)<5","colz")
0063 //--  ntp_sicalo->Draw("phi_emc-phi_pemc>>h(100,-2,2)","e_emc>0.2&&pt>0.5&&pt<2&&abs(z_emc-z_pemc)<5","colz")
0064 //--  ntp_sicalo->Draw("phi_emc-phi_pemc+(-1.6<phi_pemc&&phi_pemc<1.6?-0.2:0.2):phi_pemc>>h(100,-3.2,3.2,100,-2,2)","e_emc>0.2&&pt>0.5&&pt<2&&abs(z_emc-z_pemc)<5","colz")
0065 //--  ntp_sicalo->Draw("phi_emc-phi_pemc+(-1.6<phi_pemc&&phi_pemc<1.6?-0.2:0.2)>>h(100,-2,2)","e_emc>0.2&&pt>0.5&&pt<2&&abs(z_emc-z_pemc)<5","colz")
0066 //--  ntp_sicalo->Draw("phi_emc-phi_pemc+(-1.6<phi_pemc&&phi_pemc<1.6?-0.2:0.2):pt>>h(100,0,5,100,-2,2)","e_emc>0.2&&pt>0.&&abs(z_emc-z_pemc)<5","colz")
0067 //--  ntp_sicalo->Draw("phi_emc-phi_pemc>>h(100,-2,2)","e_emc>0.2&&pt>0.5&&pt<2&&abs(z_emc-z_pemc)<5","colz")
0068 //--  ntp_sicalo->Draw("phi_emc-phi_pemc:phi_pemc>>h(100,-3.2,3.2,100,-2,2)","e_emc>0.2&&pt>0.5&&pt<2&&abs(z_emc-z_pemc)<5","colz")
0069 //--  ntp_sicalo->Draw("phi_emc-phi_pemc>>h(100,-2,2)","e_emc>0.2&&pt>0.5&&pt<2&&abs(z_emc-z_pemc)<5","colz")
0070 //--  ntp_sicalo->Draw("phi_emc-phi_pemc+(-1.6<phi_pemc&&phi_pemc<1.6?-0.15:0.25):pt>>h(100,0,5,100,-2,2)","e_emc>0.2&&pt>0.&&abs(z_emc-z_pemc)<5","colz")
0071 //--  ntp_sicalo->Draw("phi_emc-phi_pemc+(-1.6<phi_pemc&&phi_pemc<1.6?-0.15:0.25):phi_pemc>>h(100,-3.2,3.2,100,-2,2)","e_emc>0.2&&pt>0.&&abs(z_emc-z_pemc)<5","colz")
0072 //--  ntp_sicalo->Draw("phi_emc-phi_pemc+(-1.6<phi_pemc&&phi_pemc<1.6?-0.:0.):phi_pemc>>h(100,-3.2,3.2,100,-2,2)","e_emc>0.2&&pt>0.&&abs(z_emc-z_pemc)<5","colz")
0073 //--  ntp_sicalo->Draw("phi_emc-phi_pemc+(-1.6<phi_pemc&&phi_pemc<1.6?-0.15:0.25):phi_pemc>>h(100,-3.2,3.2,100,-2,2)","e_emc>0.2&&pt>0.&&abs(z_emc-z_pemc)<5","colz")
0074 //--  ntp_sicalo->Draw("phi_emc-phi_pemc+(-1.6<phi_pemc&&phi_pemc<1.6?-0.15:0.25)>>h(100,-2,2)","e_emc>0.2&&pt>0.&&abs(z_emc-z_pemc)<5","colz")
0075 //--  ntp_sicalo->Draw("phi_emc-phi_pemc+(-1.6<phi_pemc&&phi_pemc<1.6?-0.15:0.25):phi_pemc>h(100,-3.2,3.2,100,-2,2)","e_emc>0.2&&pt>0.&&abs(z_emc-z_pemc)<5","colz")
0076 //--  ntp_sicalo->Draw("phi_emc-phi_pemc+(-1.6<phi_pemc&&phi_pemc<1.6?-0.15:0.25):phi_pemc>h(100,-3.2,3.2,100,-2,2)","e_emc>0.2&&pt>0.&&abs(z_emc-z_pemc)<5","colz")
0077 //--  ntp_sicalo->Draw("phi_emc-phi_pemc+(-1.6<phi_pemc&&phi_pemc<1.6?-0.15:0.25):phi_pemc>>h(100,-3.2,3.2,100,-2,2)","e_emc>0.2&&pt>0.&&abs(z_emc-z_pemc)<5","colz")
0078 //--  ntp_sicalo->Draw("phi_emc-phi_pemc+(-1.6<phi_pemc&&phi_pemc<1.6?-0.15:0.25):pt>>h(100,0,5,100,-2,2)","e_emc>0.2&&pt>0.&&abs(z_emc-z_pemc)<5","colz")
0079 //--  ntp_sicalo->Draw("c*(phi_emc-phi_pemc+(-1.6<phi_pemc&&phi_pemc<1.6?-0.15:0.25)):pt>>h(100,0,5,100,-2,2)","e_emc>0.2&&pt>0.&&abs(z_emc-z_pemc)<5","colz")
0080 //--  ntp_sicalo->Draw("(phi_emc-phi_pemc+(-1.6<phi_pemc&&phi_pemc<1.6?-0.15:0.25)):pt>>h(100,0,5,100,-2,2)","e_emc>0.2&&pt>0.&&abs(z_emc-z_pemc)<5&&c>0","colz")
0081 //--  ntp_sicalo->Draw("(phi_emc-phi_pemc+(-1.6<phi_pemc&&phi_pemc<1.6?-0.15:0.25)):pt>>h(100,0,5,100,-2,2)","e_emc>0.2&&pt>0.&&abs(z_emc-z_pemc)<5&&c<0","colz")
0082 //--  ntp_sicalo->Draw("(phi_emc-phi_pemc+(-1.6<phi_pemc&&phi_pemc<1.6?-0.15:0.25)):pt>>h(100,0,5,100,-2,2)","e_emc>0.2&&pt>0.&&abs(z_emc-z_pemc)<5&&c>0","colz")
0083 //--  ntp_sicalo->Draw("(phi_emc-phi_pemc+(-1.6<phi_pemc&&phi_pemc<1.6?-0.15:0.25)):pt>>h(100,0,5,100,-2,2)","e_emc>0.2&&pt>0.&&abs(z_emc-z_pemc)<5&&c<0","colz")
0084 //--  ntp_sicalo->Draw("(phi_emc-phi_pemc+(-1.6<phi_pemc&&phi_pemc<1.6?-0.15:0.25)):pt>>h(100,0,5,100,-2,2)","e_emc>0.2&&pt>0.&&abs(z_emc-z_pemc)<5&&c>0","colz")
0085 //--  ntp_sicalo->Draw("(phi_emc-phi_pemc+(-1.6<phi_pemc&&phi_pemc<1.6?-0.15:0.25)):pt>>h(100,0,5,100,-2,2)","e_emc>0.2&&pt>0.&&abs(z_emc-z_pemc)<5&&c<0","colz")
0086 //--  ntp_sicalo->Draw("(phi_emc-phi_pemc+(-1.6<phi_pemc&&phi_pemc<1.6?-0.15:0.25)):pt>>h(100,0,5,100,-2,2)","e_emc>0.2&&pt>0.&&abs(z_emc-z_pemc)<5&&c>0","colz")
0087 //--  ntp_sicalo->Draw("(phi_emc-phi_pemc+(-1.6<phi_pemc&&phi_pemc<1.6?-0.15:0.25)):pt>>h(100,0,5,100,-2,2)","e_emc>0.2&&pt>0.&&abs(z_emc-z_pemc)<5&&c<0","colz")
0088 //--  ntp_sicalo->Draw("(phi_emc-phi_pemc+(-1.6<phi_pemc&&phi_pemc<1.6?-0.15:0.25)):pt>>h(100,0,5,100,-2,2)","e_emc>0.2&&pt>0.&&abs(z_emc-z_pemc)<5&&c<0&&abs(phi_pemc)<1.6","colz")
0089 //--  ntp_sicalo->Draw("(phi_emc-phi_pemc+(-1.6<phi_pemc&&phi_pemc<1.6?-0.15:0.25)):pt>>h(100,0,5,100,-2,2)","e_emc>0.2&&pt>0.&&abs(z_emc-z_pemc)<5&&c>0&&abs(phi_pemc)<1.6","colz")
0090 //--  ntp_sicalo->Draw("(phi_emc-phi_pemc+(-1.6<phi_pemc&&phi_pemc<1.6?-0.15:0.25)):pt>>h(100,0,5,100,-2,2)","e_emc>0.2&&pt>0.&&abs(z_emc-z_pemc)<5&&c>0&&abs(phi_pemc)>1.6","colz")
0091 //--  ntp_sicalo->Draw("(phi_emc-phi_pemc+(-1.6<phi_pemc&&phi_pemc<1.6?-0.15:0.25)):pt>>h(100,0,5,100,-2,2)","e_emc>0.2&&pt>0.&&abs(z_emc-z_pemc)<5&&c<0&&abs(phi_pemc)>1.6","colz")
0092 //--  ntp_sicalo->Draw("c(phi_emc-phi_pemc+(-1.6<phi_pemc&&phi_pemc<1.6?-0.15:0.25)):pt>>h(100,0,5,100,-2,2)","e_emc>0.2&&pt>0.&&abs(z_emc-z_pemc)<5&&c<0&&abs(phi_pemc)>1.6","colz")
0093 //--  ntp_sicalo->Draw("c*(phi_emc-phi_pemc+(-1.6<phi_pemc&&phi_pemc<1.6?-0.15:0.25)):pt>>h(100,0,5,100,-2,2)","e_emc>0.2&&pt>0.&&abs(z_emc-z_pemc)<5&&c<0&&abs(phi_pemc)>1.6","colz")
0094 //--  ntp_sicalo->Draw("c*(phi_emc-phi_pemc+(-1.6<phi_pemc&&phi_pemc<1.6?-0.15:0.25)):pt>>h(100,0,5,100,-2,2)","e_emc>0.2&&pt>0.&&abs(z_emc-z_pemc)<5&&abs(phi_pemc)>1.6","colz")
0095 //--  ntp_sicalo->Draw("c*(phi_emc-phi_pemc+(-1.6<phi_pemc&&phi_pemc<1.6?-0.15:0.25)):pt>>h(100,0,5,100,-2,2)","e_emc>0.2&&pt>0.&&abs(z_emc-z_pemc)<5&&c<0&&abs(phi_pemc)<1.6","colz")
0096 //--  ntp_sicalo->Draw("c*(phi_emc-phi_pemc+(-1.6<phi_pemc&&phi_pemc<1.6?-0.15:0.25)):pt>>h(100,0,5,100,-2,2)","e_emc>0.2&&pt>0.&&abs(z_emc-z_pemc)<5&&c<0&&abs(phi_pemc)>1.6","colz")
0097 //--  ntp_sicalo->Draw("(phi_emc-phi_pemc+(-1.6<phi_pemc&&phi_pemc<1.6?-0.15:0.25)):pt>>h(100,0,5,100,-2,2)","e_emc>0.2&&pt>0.&&abs(z_emc-z_pemc)<5&&c<0&&abs(phi_pemc)>1.6&&c>0","colz")
0098 //--  ntp_sicalo->Draw("(phi_emc-phi_pemc+(-1.6<phi_pemc&&phi_pemc<1.6?-0.15:0.25)):pt>>h(100,0,5,100,-2,2)","e_emc>0.2&&pt>0.&&abs(z_emc-z_pemc)<5&&abs(phi_pemc)>1.6&&c>0","colz")
0099 //--  ntp_sicalo->Draw("(phi_emc-phi_pemc+(-1.6<phi_pemc&&phi_pemc<1.6?-0.15:0.25)):pt>>h(100,0,5,100,-2,2)","e_emc>0.2&&pt>0.&&abs(z_emc-z_pemc)<5&&abs(phi_pemc)>1.6&&c<0","colz")
0100 //--  ntp_sicalo->Draw("(phi_emc-phi_pemc+(-1.6<phi_pemc&&phi_pemc<1.6?-0.15:0.25)):pt>>h(100,0,5,100,-2,2)","e_emc>0.2&&pt>0.&&abs(z_emc-z_pemc)<5&&abs(phi_pemc)<1.6&&c<0","colz")
0101 //--  ntp_sicalo->Draw("(phi_emc-phi_pemc+(-1.6<phi_pemc&&phi_pemc<1.6?-0.15:0.25)):pt>>h(100,0,5,100,-2,2)","e_emc>0.2&&pt>0.&&abs(z_emc-z_pemc)<5&&abs(phi_pemc)<3.2&&c<0","colz")
0102 //--  ntp_sicalo->Draw("(phi_emc-phi_pemc+(-1.6<phi_pemc&&phi_pemc<1.6?-0.15:0.25)):pt>>h(100,0,5,100,-2,2)","e_emc>0.2&&pt>0.&&abs(z_emc-z_pemc)<5&&abs(phi_pemc)<1.6&&c<0","colz")
0103 //--  ntp_sicalo->Draw("(phi_emc-phi_pemc+(-1.6<phi_pemc&&phi_pemc<1.6?-0.15:0.25)):pt>>h(100,0,5,100,-2,2)","e_emc>0.2&&pt>0.&&abs(z_emc-z_pemc)<5&&abs(phi_pemc)>1.6&&c<0","colz")
0104 //--  ntp_sicalo->Draw("(phi_emc-phi_pemc+(-1.6<phi_pemc&&phi_pemc<1.6?-0.15:0.25)):pt>>h(100,0,5,100,-2,2)","e_emc>0.2&&pt>0.&&abs(z_emc-z_pemc)<5&&abs(phi_pemc)<1.6&&c<0","colz")
0105 //--  ntp_sicalo->Draw("(phi_emc-phi_pemc+(-1.6<phi_pemc&&phi_pemc<1.6?-0.15:0.25)):pt>>h(100,0,5,100,-2,2)","e_emc>0.2&&pt>0.&&abs(z_emc-z_pemc)<5&&abs(phi_pemc)>1.6&&c<0","colz")
0106 }
0107 
0108 //--.x Fun4All_DataDST_SiliconSeedAna.C
0109 //--TFile *_file0 = TFile::Open("test_ana.root")
0110 //--evtTree->Draw("nmaps:nintt50","", "colz")
0111 //--.q
0112 //--.x Fun4All_DataDST_SiliconSeedAna.C
0113 //--TFile *_file0 = TFile::Open("test_ana.root")
0114 //--evtTree->Draw("nmaps:nintt50","", "colz")
0115 //--.ls
0116 //--evtTree->Print()
0117 //--evtTree->Draw("xvtx")
0118 //--evtTree->Draw("yvtx")
0119 //--evtTree->Draw("zvtx")
0120 //--evtTree->Draw("yvtx:xvtx>>h(100,-1,1,100,-1,1)","","colz")
0121 //--evtTree->Draw("yvtx:xvtx>>h(200,-1,1,200,-1,1)","","colz")
0122 //--evtTree->Draw("yvtx:xvtx>>h(200,-0.5,0.5,200,-0.5,0.5)","","colz")
0123 //--evtTree->Draw("xvtx>>h(200,-0.5,0.5)","","colz")
0124 //--evtTree->Draw("xvtx>>h(400,-0.5,0.5)","","colz")
0125 //--evtTree->Draw("yvtx>>h(400,-0.5,0.5)","","colz")
0126 //--evtTree->Draw("zvtx>>h(400,-20,20)","","colz")
0127 //--evtTree->Draw("zvtx>>h(100,-20,20)","","colz")
0128 //--evtTree->Draw("nmaps:nintt50","abs(zvtx)<10", "colz")
0129 //--evtTree->Draw("nmaps:nintt50>>h(100,0,500,100,0,500)","abs(zvtx)<10", "colz")
0130 //--evtTree->Draw("nmaps:nintt50>>h(100,0,1000,100,0,1000)","abs(zvtx)<10", "colz")
0131 //--evtTree->Draw("nmaps:nintt50>>h(100,0,1000,100,0,1000)","abs(zvtx)<5", "colz")
0132 //--evtTree->Draw("nmaps:nintt50>>h(100,0,1000,100,0,1000)","abs(zvtx)<10", "colz")
0133 //--evtTree->Draw("nmaps:nintt50>>h(100,0,1000,100,0,1000)","abs(zvtx)<30", "colz")
0134 //--.q
0135 //--.x Fun4All_DataDST_SiliconSeedAna.C
0136 //--.x analyze_SiSeedPair.C
0137 //--.x analyze_SiSeedPair.C
0138 //--TFile *_file0 = TFile::Open("dphi_distribution.root")
0139 //--.ls
0140 //--ntp_sicalo->Print()
0141 //--ntp_sicalo->Draw("phi_emc:phi_pemc")
0142 //--ntp_sicalo->Draw("phi_emc:phi_pemc>>h(100,-3.2,3.2,100,-3.2,3.2)","","colz")
0143 //--ntp_sicalo->Draw("phi","","colz")
0144 //--ntp_sicalo->Draw("phi_emc:phi_pemc>>h(100,-3.2,3.2,100,-3.2,3.2)","","colz")
0145 //--ntp_sicalo->Draw("phi_emc:phi_pemc>>h(100,-3.2,3.2,100,-3.2,3.2)","e_emc>0.2","colz")
0146 //--ntp_sicalo->Draw("e_emc","","colz")
0147 //--ntp_sicalo->Draw("phi_emc:phi_pemc>>h(100,-3.2,3.2,100,-3.2,3.2)","e_emc>0.2","colz")
0148 //--ntp_sicalo->Draw("phi_emc:phi_pemc>>h(320,-3.2,3.2,320,-3.2,3.2)","e_emc>0.2","colz")
0149 //--ntp_sicalo->Draw("phi_emc:phi_pemc>>h(160,-3.2,3.2,160,-3.2,3.2)","e_emc>0.2","colz")
0150 //--ntp_sicalo->Draw("phi_emc>>h(160,-3.2,3.2,160,-3.2,3.2)","e_emc>0.2","colz")
0151 //--.ls
0152 //--.q
0153 //--TFile *_file0 = TFile::Open("../macro/data//test_ana2k.root")
0154 //--.ls
0155 //--caloTree->Print()
0156 //--caloTree->Draw("phi")
0157 //--caloTree->Draw("phi>>h(160,-3,2,3,2)")
0158 //--caloTree->Draw("phi>>h(175,-3,5,3,5)")
0159 //--caloTree->Draw("phi>>h(175,-3.5,3.5)")
0160 //--caloTree->Draw("phi:z>>h(100,-200,200,175,-3.5,3.5)","","colz")
0161 //--caloTree->Draw("z>>h(100,-200,200)","","colz")
0162 //--.q
0163 //--TFile *_file0 = TFile::Open("dphi_distribution.root")
0164 //--ntp_sicalo->Print()
0165 //--ntp_sicalo->Draw("phi_emc:phi_pemc>>h(320,-3.2,3.2,320,-3.2,3.2)","e_emc>0.2&&pt>0.5","colz")
0166 //--ntp_sicalo->Draw("phi_emc:phi_pemc>>h(160,-3.2,3.2,160,-3.2,3.2)","e_emc>0.2&&pt>0.5","colz")
0167 //--ntp_sicalo->Draw("phi_emc-phi_pemc","e_emc>0.2&&pt>0.5","colz")
0168 //--ntp_sicalo->Draw("z_emc-z_pemc","e_emc>0.2&&pt>0.5","colz")
0169 //--ntp_sicalo->Draw("z_emc:z_pemc","e_emc>0.2&&pt>0.5","colz")
0170 //--ntp_sicalo->Draw("z_emc:z_pemc>>h(100,-150,150,100,-150,150)","e_emc>0.2&&pt>0.5","colz")
0171 //--ntp_sicalo->Draw("z_emc:z_pemc>>h(150,-150,150,150,-150,150)","e_emc>0.2&&pt>0.5","colz")
0172 //--ntp_sicalo->Draw("z_emc:z_pemc>>h(100,-150,150,100,-150,150)","e_emc>0.2&&pt>0.5","colz")
0173 //--ntp_sicalo->Draw("phi_emc:phi_pemc>>h(160,-3.2,3.2,160,-3.2,3.2)","e_emc>0.2&&pt>0.5","colz")
0174 //--ntp_sicalo->Draw("phi_emc:phi_pemc>>h(80,-3.2,3.2,80,-3.2,3.2)","e_emc>0.2&&pt>0.5","colz")
0175 //--ntp_sicalo->Draw("phi_emc:phi_pemc>>h(160,-3.2,3.2,160,-3.2,3.2)","e_emc>0.2&&pt>0.5","colz")
0176 //--ntp_sicalo->Draw("z_emc-z_pemc:z_pemc>>h(100,-150,150,100,-150,150)","e_emc>0.2&&pt>0.5","colz")
0177 //--ntp_sicalo->Draw("z_emc-z_pemc:z_pemc>>h(100,-150,150,100,-20,20)","e_emc>0.2&&pt>0.5","colz")
0178 //--ntp_sicalo->Draw("z_emc-z_pemc:z_pemc>>h(100,-150,150,100,-50,50)","e_emc>0.2&&pt>0.5","colz")
0179 //--ntp_sicalo->Draw("z_emc-z_pemc>>h(100,-50,50)","e_emc>0.2&&pt>0.5","colz")
0180 //--ntp_sicalo->Draw("phi_emc:phi_pemc>>h(160,-320,320,160,-320,320)","e_emc>0.2&&pt>0.5","colz")
0181 //--ntp_sicalo->Draw("phi_emc:phi_pemc>>h(160,-3.20,3.20,160,-3.20,3.20)","e_emc>0.2&&pt>0.5","colz")
0182 //--ntp_sicalo->Draw("phi_emc:phi_pemc>>h(160,-3.20,3.20,160,-3.20,3.20)","e_emc>0.2&&pt>0.5&&abs(z_emc-z_pemc)<5","colz")
0183 //--ntp_sicalo->Draw("phi_emc-phi_pemc>>h(100,-1,1)","e_emc>0.2&&pt>0.5&&abs(z_emc-z_pemc)<5","colz")
0184 //--ntp_sicalo->Draw("phi_emc-phi_pemc>>h(100,-2,2)","e_emc>0.2&&pt>0.5&&abs(z_emc-z_pemc)<5","colz")
0185 //--ntp_sicalo->Draw("phi_emc-phi_pemc:phi_pemc>>h(100,-3.2,3.2,100,-2,2)","e_emc>0.2&&pt>0.5&&abs(z_emc-z_pemc)<5","colz")
0186 //--ntp_sicalo->Print()
0187 //--ntp_sicalo->Draw("e_emc:sqrt(pt*pt+pz*pz)","","colz"
0188 //--.@
0189 //--ntp_sicalo->Draw("e_emc:sqrt(pt*pt+pz*pz)","","colz")
0190 //--ntp_sicalo->Draw("e_emc:sqrt(pt*pt+pz*pz)","pt<10","colz")
0191 //--ntp_sicalo->Draw("e_emc:sqrt(pt*pt+pz*pz)","pt<10&&pz<10","colz")
0192 //--ntp_sicalo->Draw("sqrt(pt*pt+pz*pz)","pt<10&&pz<10","colz")
0193 //--ntp_sicalo->Draw("e_emc:sqrt(pt*pt+pz*pz)","pt<10&&pz<10","colz")
0194 //--ntp_sicalo->Draw("e_emc/sqrt(pt*pt+pz*pz)","pt<10&&pz<10","colz")
0195 //--ntp_sicalo->Draw("e_emc/sqrt(pt*pt+pz*pz)>>h(100,0,3)","pt<10&&pz<10","colz")
0196 //--ntp_sicalo->Draw("e_emc/sqrt(pt*pt+pz*pz):pt>>h(100,0,10,100,0,3)","pt<10&&pz<10","colz")
0197 //--ntp_sicalo->Draw("e_emc/sqrt(pt*pt+pz*pz):sqrt(pt*pt+pz*pz)>>h(100,0,10,100,0,3)","pt<10&&pz<10","colz")
0198 //--ntp_sicalo->Draw("e_emc/sqrt(pt*pt+pz*pz):sqrt(pt*pt+pz*pz)>>h(500,0,50,100,0,3)","pt<10&&pz<10","colz")
0199 //--ntp_sicalo->Draw("e_emc/sqrt(pt*pt+pz*pz):sqrt(pt*pt+pz*pz)>>h(50,0,50,100,0,3)","pt<10&&pz<10","colz")
0200 //--ntp_sicalo->Draw("e_emc/sqrt(pt*pt+pz*pz):sqrt(pt*pt+pz*pz)>>h(50,0,5,100,0,3)","pt<10&&pz<10","colz")
0201 //--ntp_sicalo->Draw("e_emc/sqrt(pt*pt+pz*pz):sqrt(pt*pt+pz*pz)>>h(50,0,50,100,0,3)","pt<10&&pz<10&&abs(z_emc-z_pemc)<5","colz")
0202 //--ntp_sicalo->Draw("e_emc/sqrt(pt*pt+pz*pz):sqrt(pt*pt+pz*pz)>>h(50,0,5,100,0,3)","pt<10&&pz<10&&abs(z_emc-z_pemc)<5","colz")
0203 //--.q
0204 //--TFile *_file0 = TFile::Open("../macro/data/test_ana.root")
0205 //--.ls
0206 //--SiClusTree->Print()
0207 //--SiClusTree->Draw("Siclus_z:layer","","colz")
0208 //--SiClusTree->Draw("Siclus_z:Siclus_layer","","colz")
0209 //--SiClusTree->Draw("Siclus_z:Siclus_layer>>h(8,0,8,100,-30,30)","","colz")
0210 //--SiClusTree->Draw("Siclus_z:Siclus_layer>>h(8,0,8,300,-30,30)","","colz")
0211 //--SiClusTree->Draw("Siclus_y:Siclus_x>>h(300,-30,30,300,-30,30)","","colz")
0212 //--SiClusTree->Draw("Siclus_y:Siclus_x>>h(300,-15,15,300,-15,15)","","colz")
0213 //--SiClusTree->Draw("Siclus_z:Siclus_layer>>h(8,0,8,300,-30,30)","","colz")
0214 //--.q
0215 //--TFile *_file0 = TFile::Open("/sphenix/lustre01/sphnxpro/production2/run2pp/physics/new_nocdbtag_v007/run_00052000_00052100/DST_TRIGGERED_EVENT_seb00_run2pp_new_nocdbtag_v007-00052077-00000.root")
0216 //--.q
0217 //--TFile *_file0 = TFile::Open("test_ana.root")
0218 //--ls
0219 //--.ls
0220 //--evtTree->Print()
0221 //--.q
0222 //--.x Fun4All_DataDST_SiliconSeedAna.C
0223 //--.x analyze_SiSeedPair.C
0224 //--TFile *_file0 = TFile::Open("dphi_distribution.root")
0225 //--SiClusTree->Draw("Siclus_y:Siclus_x>>h(300,-15,15,300,-15,15)","","colz")
0226 //--.q