rtang/ANASEN_analysis
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

modified: Analyzer.C

Browse Source 
modified:   mapping.h
This commit is contained in:
dirac 2024-10-30 09:28:01 -04:00
parent a10081ea81
commit 7a70340b18
2 changed files with 84 additions and 67 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

148
Analyzer.C
View File

@ -33,6 +33,7 @@ TH2F *hAnodeHits;
TH1F *hAnodeHits1d; TH1F *hAnodeHits1d;
TH1F *hPCMultiplicity; TH1F *hPCMultiplicity;
TH1F *hRFtime; TH1F *hRFtime;
TH1F *hSi;
int padID = 0; int padID = 0;
@ -46,9 +47,11 @@ TCutG *PCCoinc;
TCutG *alpha_cut_up; TCutG *alpha_cut_up;
TCutG *alpha_cut_down; TCutG *alpha_cut_down;
TCutG *cutg;
bool inCut; bool inCut;
bool inCutUp; bool inCutUp;
bool inCutDown; bool inCutDown;
bool inCutG;
void Analyzer::Begin(TTree * /*tree*/) { void Analyzer::Begin(TTree * /*tree*/) {
TString option = GetOption(); TString option = GetOption();
@ -87,7 +90,8 @@ void Analyzer::Begin(TTree * /*tree*/) {
hPCMultiplicity = new TH1F("hPCMultiplicity", "Number of PC/Event", 40, 0, 40); hPCMultiplicity = new TH1F("hPCMultiplicity", "Number of PC/Event", 40, 0, 40);
hanVScatsum = new TH2F("hanVScatsum", "Anode vs Cathode Sum; Anode E; Cathode E", 6400, 0, 20000, 6400, 0, 20000); hanVScatsum = new TH2F("hanVScatsum", "Anode vs Cathode Sum; Anode E; Cathode E", 6400, 0, 20000, 6400, 0, 20000);
hICvsSi = new TH2F("hICvsSi", "IC vs Si; Si E; IC E", 800, 0, 20000, 400, 0, 8000); hICvsSi = new TH2F("hICvsSi", "IC vs Si; Si E; IC E", 800, 0, 20000, 400, 0, 8000);
hRFtime = new TH1F("hRFtime", "RF time", 1000, 0, 200000); hSi = new TH1F("hSi", "Si E", 800, 0, 20000);
hRFtime = new TH1F("hRFtime", "RF time (ns)", 500, 0, 3000);
sx3_contr.ConstructGeo(); sx3_contr.ConstructGeo();
pw_contr.ConstructGeo(); pw_contr.ConstructGeo();
@ -97,6 +101,8 @@ void Analyzer::Begin(TTree * /*tree*/) {
alpha_cut_up = (TCutG *)f2->Get("alpha_cut_up"); alpha_cut_up = (TCutG *)f2->Get("alpha_cut_up");
TFile *f3 = new TFile("alpha_cut_down.root"); TFile *f3 = new TFile("alpha_cut_down.root");
alpha_cut_down = (TCutG *)f3->Get("alpha_cut_down"); alpha_cut_down = (TCutG *)f3->Get("alpha_cut_down");
TFile *f4 = new TFile("CUTG.root");
cutg = (TCutG *)f4->Get("CUTG");
// TFile *f1 = new TFile("AnCatSum.root"); // TFile *f1 = new TFile("AnCatSum.root");
// AnCatSum= (TCutG*)f1->Get("AnCatSum"); // AnCatSum= (TCutG*)f1->Get("AnCatSum");
@ -395,22 +401,34 @@ Bool_t Analyzer::Process(Long64_t entry) {
} }
// Miscellaneous channels including the Lollipop IC and Si detectors and hot needle IC // Miscellaneous channels including the Lollipop IC and Si detectors and hot needle IC
bool timing = false;
inCutG = false;
for (int i = 0; i < misc.multi; i++) { for (int i = 0; i < misc.multi; i++) {
if (misc.ch[i] == 1) {
if(misc.e[i] > 7500 && misc.e[i]<15000) hSi->Fill(misc.e[i]);
inCutG = true;
}
for (int j = i + 1; j < misc.multi; j++) { for (int j = i + 1; j < misc.multi; j++) {
// if (misc.ch[i] == 0 || misc.ch[j] == 1) { if (cutg->IsInside(misc.e[i], misc.e[j])) {
hICvsSi->Fill(misc.e[i], misc.e[j]); inCutG = true;
// } }
if (misc.ch[j] == 2 && inCutG ) {
hRFtime->Fill(misc.t[j] + misc.tf[j] * 4 / 1000 - (misc.t[i] + misc.tf[i] * 4 / 1000));
// if (misc.t[j] + misc.tf[j] * 4 / 1000 - (misc.t[i] + misc.tf[i] * 4 / 1000) > 1000 && misc.t[j] + misc.tf[j] * 4 / 1000 - (misc.t[i] + misc.tf[i] * 4 / 1000) < 1100) {
timing = true;
// }
// printf("RF time : %lld %lld %lld %lld %lld\n", misc.t[i], misc.t[j], misc.tf[i], misc.tf[j], (misc.t[j]*1000 + misc.tf[j]*4 - (misc.t[i]*1000 + misc.tf[i]*4)));
}
} }
} }
for (int i = 0; i < misc.multi; i++) { for (int i = 0; i < misc.multi; i++) {
if(misc.ch[i] == 1) hSi->Fill(misc.e[i]);
for (int j = i + 1; j < misc.multi; j++) { for (int j = i + 1; j < misc.multi; j++) {
// if (misc.ch[i] == 2 || misc.ch[j] == 4) { if (timing == true) {
hRFtime->Fill(misc.t[j]*1000 + misc.tf[j]*4 - (misc.t[i]*1000 + misc.tf[i]*4)); hICvsSi->Fill(misc.e[i], misc.e[j]);
// printf("RF time : %lld %lld %lld %lld %lld\n", misc.t[i], misc.t[j], misc.tf[i], misc.tf[j], (misc.t[j]*1000 + misc.tf[j]*4 - (misc.t[i]*1000 + misc.tf[i]*4))); }
// }
} }
} }
@ -431,84 +449,84 @@ Bool_t Analyzer::Process(Long64_t entry) {
void Analyzer::Terminate() { void Analyzer::Terminate() {
gStyle->SetOptStat("neiou"); gStyle->SetOptStat("neiou");
TCanvas *a = new TCanvas("cANASEN", "ANASEN", 800, 600); TCanvas *canvas = new TCanvas("cANASEN", "ANASEN", 2000, 2000);
// canvas->Divide(3, 3); canvas->Divide(3, 3);
hRFtime->Draw(); // hRFtime->Draw();
TCanvas *b = new TCanvas("ANASEN", "ANASEN", 800, 600); // TCanvas *b = new TCanvas("ANASEN", "ANASEN", 800, 600);
hICvsSi->Draw("colz"); // hICvsSi->Draw("colz");
// hsx3VpcIndex->Draw("colz"); // hsx3VpcIndex->Draw("colz");
//=============================================== pad-1 // =============================================== pad-1
// padID++; padID++;
// canvas->cd(padID); canvas->cd(padID);
// canvas->cd(padID)->SetGrid(1); canvas->cd(padID)->SetGrid(1);
// hsx3IndexVE->Draw("colz"); hsx3IndexVE->Draw("colz");
// //=============================================== pad-2 //=============================================== pad-2
// padID++; padID++;
// canvas->cd(padID); canvas->cd(padID);
// canvas->cd(padID)->SetGrid(1); canvas->cd(padID)->SetGrid(1);
// hqqqIndexVE->Draw("colz"); hqqqIndexVE->Draw("colz");
// //=============================================== pad-3 //=============================================== pad-3
// padID++; padID++;
// canvas->cd(padID); canvas->cd(padID);
// canvas->cd(padID)->SetGrid(1); canvas->cd(padID)->SetGrid(1);
// hpcIndexVE->Draw("colz"); hpcIndexVE->Draw("colz");
// //=============================================== pad-4 //=============================================== pad-4
// padID++; padID++;
// canvas->cd(padID); canvas->cd(padID);
// canvas->cd(padID)->SetGrid(1); canvas->cd(padID)->SetGrid(1);
// hsx3Coin->Draw("colz"); hsx3Coin->Draw("colz");
// //=============================================== pad-5 //=============================================== pad-5
// padID++; padID++;
// canvas->cd(padID); canvas->cd(padID);
// canvas->cd(padID)->SetGrid(1); canvas->cd(padID)->SetGrid(1);
// canvas->cd(padID)->SetLogz(true); canvas->cd(padID)->SetLogz(true);
// hqqqCoin->Draw("colz"); hqqqCoin->Draw("colz");
// //=============================================== pad-6 //=============================================== pad-6
// padID++; padID++;
// canvas->cd(padID); canvas->cd(padID);
// canvas->cd(padID)->SetGrid(1); canvas->cd(padID)->SetGrid(1);
// hpcCoin->Draw("colz"); hpcCoin->Draw("colz");
// //=============================================== pad-7 //=============================================== pad-7
// padID++; padID++;
// canvas->cd(padID); canvas->cd(padID);
// canvas->cd(padID)->SetGrid(1); canvas->cd(padID)->SetGrid(1);
// // hsx3VpcIndex ->Draw("colz"); // hsx3VpcIndex ->Draw("colz");
// hsx3VpcE->Draw("colz"); hsx3VpcE->Draw("colz");
// //=============================================== pad-8 //=============================================== pad-8
// padID++; padID++;
// canvas->cd(padID); canvas->cd(padID);
// canvas->cd(padID)->SetGrid(1); canvas->cd(padID)->SetGrid(1);
// // hqqqVpcIndex ->Draw("colz"); // hqqqVpcIndex ->Draw("colz");
// hqqqVpcE->Draw("colz"); hqqqVpcE->Draw("colz");
// //=============================================== pad-9 //=============================================== pad-9
// padID++; padID++;
// // canvas->cd(padID)->DrawFrame(-50, -50, 50, 50); // canvas->cd(padID)->DrawFrame(-50, -50, 50, 50);
// // hqqqPolar->Draw("same colz pol"); // hqqqPolar->Draw("same colz pol");
// canvas->cd(padID); canvas->cd(padID);
// canvas->cd(padID)->SetGrid(1); canvas->cd(padID)->SetGrid(1);
// // hZProj->Draw(); // hZProj->Draw();
// hanVScatsum->Draw("colz"); hanVScatsum->Draw("colz");
// hAnodeHits->Draw("colz"); // hAnodeHits->Draw("colz");
// hAnodeMultiplicity->Draw(); // hAnodeMultiplicity->Draw();
} }

3
mapping.h
View File

@ -85,8 +85,7 @@ const std::vector<int> mapping = {
20100, 20101, 20102, 20103, 20104, 20105, 20106, 20107, 20108, 20109, 20110, 20111, 20112, -1, 20114, 20115, 20100, 20101, 20102, 20103, 20104, 20105, 20106, 20107, 20108, 20109, 20110, 20111, 20112, -1, 20114, 20115,
//================== 405 //================== 405
// 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 // 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
20006, -1, -1, 20009, -1, 20120, 20000, 20019, 20113, 30000, -1, 30001, 30002, 30003, 30004, -1 20006, -1, -1, 20009, -1, 20120, 20000, 20019, 20113, 30000, -1, 30001, 30002, -1, 30003, 30004
}; };
void PrintMapping(){ void PrintMapping(){