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

modified: Analyzer.C implemented cuts and histograms for the 2D histograms

Browse Source
This commit is contained in:
vs19g 2024-08-26 11:41:06 -04:00
parent f0a393abe2
commit 48ede97992
1 changed files with 29 additions and 14 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

43
Analyzer.C
View File

@ -5,7 +5,7 @@
#include <TStyle.h> #include <TStyle.h>
#include <TCanvas.h> #include <TCanvas.h>
#include <TMath.h> #include <TMath.h>
#include <TCutG.h>
#include <utility> #include <utility>
#include <algorithm> #include <algorithm>
@ -29,6 +29,8 @@ TH2F * hqqqVpcE;
TH2F * hsx3VpcE; TH2F * hsx3VpcE;
TH2F * hanVScatsum; TH2F * hanVScatsum;
TH2F * hAnodeHits; TH2F * hAnodeHits;
TH1F * hAnodeMultiplicity;
int padID = 0; int padID = 0;
@ -38,6 +40,8 @@ TVector3 hitPos;
bool HitNonZero; bool HitNonZero;
TH1F * hZProj; TH1F * hZProj;
TCutG *AnCatSum;
bool inCut;
void Analyzer::Begin(TTree * /*tree*/){ void Analyzer::Begin(TTree * /*tree*/){
TString option = GetOption(); TString option = GetOption();
@ -68,14 +72,15 @@ void Analyzer::Begin(TTree * /*tree*/){
hsx3VpcE->SetNdivisions( -612, "x"); hsx3VpcE->SetNdivisions( -612, "x");
hsx3VpcE->SetNdivisions( -12, "y"); hsx3VpcE->SetNdivisions( -12, "y");
hZProj = new TH1F("hZProj", "Z Projection", 200, -600, 600); hZProj = new TH1F("hZProj", "Nos of anodes", 20, 0, 19);
hAnodeHits = new TH2F("hAnodeHits", "Anode vs Anode Energy, Anode ID; Anode E", 24,0 , 23, 400, 0 , 16000); hAnodeHits = new TH2F("hAnodeHits", "Anode vs Anode Energy, Anode ID; Anode E", 24,0 , 23, 400, 0 , 20000);
hAnodeMultiplicity = new TH1F("hAnodeMultiplicity", "Number of Anodes/Event", 40, 0, 40);
hanVScatsum = new TH2F("hanVScatsum", "Anode vs Cathode Sum; Anode E; Cathode E", 400,0 , 10000, 400, 0 , 16000); hanVScatsum = new TH2F("hanVScatsum", "Anode vs Cathode Sum; Anode E; Cathode E", 400,0 , 10000, 400, 0 , 16000);
sx3_contr.ConstructGeo(); sx3_contr.ConstructGeo();
pw_contr.ConstructGeo(); pw_contr.ConstructGeo();
TFile *f1 = new TFile("AnCatSum.root");
AnCatSum= (TCutG*)f1->Get("AnCatSum");
} }
@ -264,6 +269,7 @@ Bool_t Analyzer::Process(Long64_t entry){
std::vector<std::pair<int, double>> cathodeHits; std::vector<std::pair<int, double>> cathodeHits;
int aID = 0; int aID = 0;
int cID = 0; int cID = 0;
int anodeCount = 0;
float aE = 0; float aE = 0;
float cE = 0; float cE = 0;
@ -271,9 +277,9 @@ Bool_t Analyzer::Process(Long64_t entry){
// Define the excluded SX3 and QQQ channels // Define the excluded SX3 and QQQ channels
std::unordered_set<int> excludeSX3 = {34, 35, 36, 37, 61, 62, 67, 73, 74, 75, 76, 77, 78, 79, 80, 93, 97, 100, 103, 108, 109, 110, 111, 112}; std::unordered_set<int> excludeSX3 = {34, 35, 36, 37, 61, 62, 67, 73, 74, 75, 76, 77, 78, 79, 80, 93, 97, 100, 103, 108, 109, 110, 111, 112};
std::unordered_set<int> excludeQQQ = {0, 17, 109, 110, 111, 112, 113, 119, 127, 128}; std::unordered_set<int> excludeQQQ = {0, 17, 109, 110, 111, 112, 113, 119, 127, 128};
inCut=false;
for( int i = 0; i < pc.multi; i ++){ for( int i = 0; i < pc.multi; i ++){
if(pc.e[i]>50){ if(pc.e[i]>50 && pc.multi<10){
// hpcIndexVE->Fill( pc.index[i], pc.e[i] ); // hpcIndexVE->Fill( pc.index[i], pc.e[i] );
// for( int j = i+1; j < pc.multi; j++){ // for( int j = i+1; j < pc.multi; j++){
// hpcCoin->Fill( pc.index[i], pc.index[j]); // hpcCoin->Fill( pc.index[i], pc.index[j]);
@ -286,23 +292,27 @@ Bool_t Analyzer::Process(Long64_t entry){
for( int j = i+1; j < pc.multi; j++){ for( int j = i+1; j < pc.multi; j++){
hpcCoin->Fill( pc.index[i], pc.index[j]); hpcCoin->Fill( pc.index[i], pc.index[j]);
} }
if(pc.e[i]>100){ // if(pc.e[i]>100){
if (pc.index[i] < 24 ){ if (pc.index[i] < 24 ){
anodeHits.push_back(std::pair<int, double>(pc.index[i], pc.e[i])); anodeHits.push_back(std::pair<int, double>(pc.index[i], pc.e[i]));
// anodeCount++;
} else if (pc.index[i] >= 24){ } else if (pc.index[i] >= 24){
cathodeHits.push_back(std::pair<int, double>(pc.index[i], pc.e[i])); cathodeHits.push_back(std::pair<int, double>(pc.index[i], pc.e[i]));
} }
} // }
} }
} }
// hpcIndexVE->Fill( pc.index[i], pc.e[i] ); // hpcIndexVE->Fill( pc.index[i], pc.e[i] );
anodeCount=pc.multi;
hAnodeMultiplicity->Fill(anodeCount);
float aESum = 0; float aESum = 0;
if (anodeHits.size()>=1 && cathodeHits.size() >= 1 && pc.e[i]>100){ if (anodeHits.size()>=1 && cathodeHits.size() >= 1){
for (const auto& anode : anodeHits) { for (const auto& anode : anodeHits) {
float cESum = 0; float cESum = 0;
// for(int l=0; l<sx3.multi; l++){ // for(int l=0; l<sx3.multi; l++){
// if (sx3.index[l]==80){ // if (sx3.index[l]==80){
int aID = anode.first; int aID = anode.first;
float aE = anode.second; float aE = anode.second;
aESum += aE; aESum += aE;
@ -313,8 +323,13 @@ Bool_t Analyzer::Process(Long64_t entry){
cESum += cE; cESum += cE;
} }
hanVScatsum->Fill(aESum,cESum); if( AnCatSum->IsInside(aE, cESum)){
inCut=true;
}
if(inCut){
hanVScatsum->Fill(aE,cESum);
hAnodeHits->Fill(aID, aE); hAnodeHits->Fill(aID, aE);
}
// } // }
// } // }
} }
@ -403,7 +418,7 @@ void Analyzer::Terminate(){
canvas->cd(padID); canvas->cd(padID)->SetGrid(1); canvas->cd(padID); canvas->cd(padID)->SetGrid(1);
// hZProj->Draw(); // hZProj->Draw();
// hanVScatsum->Draw("colz"); hanVScatsum->Draw("colz");
hAnodeHits->Draw("colz"); // hAnodeHits->Draw("colz");
// hAnodeMultiplicity->Draw();?
} }