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

modified: PCGainMatch.C

Browse Source
This commit is contained in:
vs19g 2024-09-06 14:32:22 -04:00
parent 68fc36a8f6
commit 43233ceb02
2 changed files with 101 additions and 54 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
Analyzer.C
View File

@ -331,10 +331,10 @@ Bool_t Analyzer::Process(Long64_t entry){
if( AnCatSum->IsInside(aE, cESum)){ if( AnCatSum->IsInside(aE, cESum)){
inCut=true; inCut=true;
} }
if(inCut){ // if(inCut){
hanVScatsum->Fill(aE,cESum); hanVScatsum->Fill(aE,cESum);
hAnodeHits->Fill(aID, aE); hAnodeHits->Fill(aID, aE);
} // }
// } // }
// } // }
} }

151
PCGainMatch.C
View File

@ -8,7 +8,6 @@
#include <TCutG.h> #include <TCutG.h>
#include <utility> #include <utility>
#include <algorithm> #include <algorithm>
#include <iostream>
#include "Armory/ClassSX3.h" #include "Armory/ClassSX3.h"
#include "Armory/ClassPW.h" #include "Armory/ClassPW.h"
@ -28,6 +27,7 @@ TH2F * hqqqVpcIndex;
TH2F * hqqqVpcE; TH2F * hqqqVpcE;
TH2F * hsx3VpcE; TH2F * hsx3VpcE;
TH2F * hanVScatsum; TH2F * hanVScatsum;
TH2F * hanVScatsum_a[24];
TH2F * hanVScatsum_hcut; TH2F * hanVScatsum_hcut;
TH2F * hanVScatsum_lcut; TH2F * hanVScatsum_lcut;
TH2F * hAnodeHits; TH2F * hAnodeHits;
@ -44,8 +44,11 @@ bool HitNonZero;
TH1F * hZProj; TH1F * hZProj;
TCutG *AnCatSum_high; TCutG *AnCatSum_high;
TCutG *AnCatSum_low; TCutG *AnCatSum_low;
TCutG *PCCoinc_cut1;
TCutG *PCCoinc_cut2;
bool inCuth; bool inCuth;
bool inCut; bool inCutl;
bool inPCCut;
void PCGainMatch::Begin(TTree * /*tree*/){ void PCGainMatch::Begin(TTree * /*tree*/){
TString option = GetOption(); TString option = GetOption();
@ -79,16 +82,25 @@ void PCGainMatch::Begin(TTree * /*tree*/){
hZProj = new TH1F("hZProj", "Nos of anodes", 20, 0, 19); 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 , 20000); 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); 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, 800, 0 , 16000);
hanVScatsum_lcut = new TH2F("hanVScatsum_LCUT", "Anode vs Cathode Sum; Anode E; Cathode E", 400,0 , 10000, 400, 0 , 16000); for (int i = 0; i < 24; i++) {
hanVScatsum_hcut = new TH2F("hanVScatsum_HCUT", "Anode vs Cathode Sum; Anode E; Cathode E", 400,0 , 10000, 400, 0 , 16000); TString histName = Form("hAnodeVsCathode_%d", i);
TString histTitle = Form("Anode %d vs Cathode Sum; Anode E; Cathode Sum E", i);
hanVScatsum_a[i] = new TH2F(histName, histTitle, 400, 0, 10000, 400, 0, 16000);
}
hanVScatsum_lcut = new TH2F("hanVScatsum_LCUT", "Anode vs Cathode Sum; Anode E; Cathode E", 400,0 , 16000, 400, 0 , 16000);
hanVScatsum_hcut = new TH2F("hanVScatsum_HCUT", "Anode vs Cathode Sum; Anode E; Cathode E", 400,0 , 16000, 400, 0 , 16000);
sx3_contr.ConstructGeo(); sx3_contr.ConstructGeo();
pw_contr.ConstructGeo(); pw_contr.ConstructGeo();
TFile *f1 = new TFile("AnCatSum_high.root"); TFile *f1 = new TFile("AnCatSum_high.root");
TFile *f2 = new TFile("AnCatSum_low.root"); TFile *f2 = new TFile("AnCatSum_low.root");
TFile *f3 = new TFile("PCCoinc_cut1.root");
TFile *f4 = new TFile("PCCoinc_cut2.root");
AnCatSum_high= (TCutG*)f1->Get("AnCatSum_high"); AnCatSum_high= (TCutG*)f1->Get("AnCatSum_high");
AnCatSum_low= (TCutG*)f1->Get("AnCatSum_low"); AnCatSum_low= (TCutG*)f2->Get("AnCatSum_low");
PCCoinc_cut1= (TCutG*)f3->Get("PCCoinc_cut1");
PCCoinc_cut2= (TCutG*)f4->Get("PCCoinc_cut2");
} }
@ -132,12 +144,12 @@ Bool_t PCGainMatch::Process(Long64_t entry){
// //======================= PC // //======================= PC
std::vector<std::pair<int, double>> anodeHits; std::vector<std::pair<int, double>> anodeHits={};
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; int anodeCount = 0;
int cathodeCount = 0;
float aE = 0; float aE = 0;
float cE = 0; float cE = 0;
@ -145,77 +157,91 @@ Bool_t PCGainMatch::Process(Long64_t entry){
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};
inCuth=false; inCuth=false;
inCut=false; inCutl=false;
for( int i = 0; i < pc.multi; i ++){ for( int i = 0; i < pc.multi; i ++){
for(int j=0; j<pc.multi;j++){
if(pc.id[j]==0){
anodeCount++;
}
}
if(pc.e[i]>50){ if(pc.e[i]>50 && pc.multi<10){
// hpcIndexVE->Fill( pc.index[i], pc.e[i] );
// for( int j = i+1; j < pc.multi; j++){
// hpcCoin->Fill( pc.index[i], pc.index[j]);
// }
// for (int j=0;j<sx3.multi;j++){ float aESum = 0;
// if(excludeSX3.find(sx3.index[j]) == excludeSX3.end()){
// if(pc.e[i]>100){
float aESum = 0;
float cESum = 0; float cESum = 0;
if (pc.index[i] < 24 ) {
anodeHits.push_back(std::pair<int, double>(pc.index[i], pc.e[i]));
anodeCount++;
} else if (pc.index[i] >= 24) {
cathodeHits.push_back(std::pair<int, double>(pc.index[i], pc.e[i]));
cathodeCount++;
}
if (anodeHits.size() == 1 && cathodeHits.size() >= 1) { for(int j=i+1;j<pc.multi;j++){
if(PCCoinc_cut1->IsInside(pc.index[i], pc.index[j]) || PCCoinc_cut2->IsInside(pc.index[i], pc.index[j])){
// hpcCoin->Fill(pc.index[i], pc.index[j]);
inPCCut = true;
}
hpcCoin->Fill(pc.index[i], pc.index[j]);
}
if (anodeHits.size()>=1 && cathodeHits.size() >= 1) {
// Accumulate total energy from anode hits // Accumulate total energy from anode hits
for (const auto& anode : anodeHits) { for (const auto& anode : anodeHits) {
if(inPCCut){
int aID = anode.first; int aID = anode.first;
float aE = anode.second; float aE = anode.second;
aESum += aE; // Sum the anode energy aESum += aE; // Sum the anode energy
hAnodeHits->Fill(aID, aE);
}
} }
// Accumulate total energy from cathode hits // Accumulate total energy from cathode hits
for (const auto& cathode : cathodeHits) { for (const auto& cathode : cathodeHits) {
if(inPCCut){
int cID = cathode.first; int cID = cathode.first;
float cE = cathode.second; float cE = cathode.second;
cESum += cE; // Sum the cathode energy cESum += cE; // Sum the cathode energy
}
} }
inCuth = false; inCuth = false;
inCut = false; inCutl = false;
inPCCut = false;
// for(int j=i+1;j<pc.multi;j++){
// if(PCCoinc_cut1->IsInside(pc.index[i], pc.index[j]) || PCCoinc_cut2->IsInside(pc.index[i], pc.index[j])){
// // hpcCoin->Fill(pc.index[i], pc.index[j]);
// inPCCut = true;
// }
// hpcCoin->Fill(pc.index[i], pc.index[j]);
// }
// Check if the accumulated energies are within the defined ranges // Check if the accumulated energies are within the defined ranges
if (AnCatSum_high && AnCatSum_high->IsInside(aESum, cESum)) { if (AnCatSum_high && AnCatSum_high->IsInside(aESum, cESum)) {
inCuth = true; inCuth = true;
} }
if (AnCatSum_low && AnCatSum_low->IsInside(aESum, cESum)) { if (AnCatSum_low && AnCatSum_low->IsInside(aESum, cESum)) {
inCut = true; inCutl = true;
} }
// Fill histograms based on the cut conditions // Fill histograms based on the cut conditions
if (inCuth) { if (inCuth && inPCCut) {
hanVScatsum_hcut->Fill(aESum, cESum); hanVScatsum_hcut->Fill(aESum, cESum);
} }
if (inCut) { if (inCutl && inPCCut) {
hanVScatsum_lcut->Fill(aESum, cESum); hanVScatsum_lcut->Fill(aESum, cESum);
} }
// for(auto anode : anodeHits){
// Fill additional histograms
hAnodeMultiplicity->Fill(anodeHits.size()); // float aE = anode.second;
hanVScatsum->Fill(aESum, cESum); // aESum += aE;
hanVScatsum->Fill(aESum, cESum);
if (pc.index[i] < 24 && pc.e[i] > 50) {
hanVScatsum_a[pc.index[i]]->Fill(pc.e[i], cESum);
}
for (const auto& anode : anodeHits) {
int aID = anode.first;
float aE = anode.second;
hAnodeHits->Fill(aID, aE);
}
// }
// Fill histograms for the `pc` data // Fill histograms for the `pc` data
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++) { if(inPCCut){
hpcCoin->Fill(pc.index[i], pc.index[j]); hAnodeMultiplicity->Fill(anodeHits.size());
} }
} }
@ -229,6 +255,18 @@ Bool_t PCGainMatch::Process(Long64_t entry){
std::vector<std::pair<int, int>> ID; // first = id, 2nd = index std::vector<std::pair<int, int>> ID; // first = id, 2nd = index
for( int i = 0; i < sx3.multi; i ++){ for( int i = 0; i < sx3.multi; i ++){
ID.push_back(std::pair<int, int>(sx3.id[i], i)); ID.push_back(std::pair<int, int>(sx3.id[i], i));
hsx3IndexVE->Fill( sx3.index[i], sx3.e[i] );
for( int j = i+1; j < sx3.multi; j++){
hsx3Coin->Fill( sx3.index[i], sx3.index[j]);
}
for( int j = 0; j < pc.multi; j++){
hsx3VpcIndex->Fill( sx3.index[i], pc.index[j] );
// if( sx3.ch[index] > 8 ){
// hsx3VpcE->Fill( sx3.e[i], pc.e[j] );
// }
}
} }
@ -281,7 +319,7 @@ Bool_t PCGainMatch::Process(Long64_t entry){
for( int j = 0; j < pc.multi; j++){ for( int j = 0; j < pc.multi; j++){
// hsx3VpcIndex->Fill( sx3.index[i], pc.index[j] ); // hsx3VpcIndex->Fill( sx3.index[i], pc.index[j] );
if( sx3.ch[index] > 8 && pc.index[j]<24 && pc.e[j]>50 && inCuth){ if( sx3.ch[index] > 8 && pc.index[j]<24 && pc.e[j]>50 && inCuth ){
hsx3VpcE->Fill( sx3.e[i], pc.e[j] ); hsx3VpcE->Fill( sx3.e[i], pc.e[j] );
// hpcIndexVE->Fill( pc.index[i], pc.e[i] ); // hpcIndexVE->Fill( pc.index[i], pc.e[i] );
} }
@ -302,15 +340,24 @@ Bool_t PCGainMatch::Process(Long64_t entry){
for( int j = 0; j < pc.multi; j++){ for( int j = 0; j < pc.multi; j++){
if(pc.index[j]<24 && pc.e[j]>50 ){ if(pc.index[j]<24 && pc.e[j]>50 ){
hqqqVpcE->Fill( qqq.e[i], pc.e[j] ); hqqqIndexVE->Fill( qqq.index[i], qqq.e[i] );
// hpcIndexVE->Fill( pc.index[i], pc.e[i] ); for( int j = 0; j < qqq.multi; j++){
if ( j == i ) continue;
hqqqCoin->Fill( qqq.index[i], qqq.index[j]);
}
for( int j = i + 1; j < qqq.multi; j++){
for( int k = 0; k < pc.multi; k++){
if(pc.index[k]<24 && pc.e[k]>50 && inCuth){
hqqqVpcE->Fill( qqq.e[i], pc.e[k] );
hqqqVpcIndex->Fill( qqq.index[i], pc.index[j] ); hqqqVpcIndex->Fill( qqq.index[i], pc.index[j] );
} }
// } }
}
}
} }
}
}
// hanVScatsum->Fill(aE,cE); // hanVScatsum->Fill(aE,cE);