Skip to content
Snippets Groups Projects
Commit 5dda4eef authored by Theodore Efremov's avatar Theodore Efremov :hibiscus:
Browse files

Modified TTime and filled it in the convert to add the section of the MWPC and did some analysis

parent 816060ba
No related branches found
No related tags found
1 merge request!27Draft: [Epic] Preparation of the environement for the new GaseousDetectorScorers...
Hide whitespace changes
Inline Side-by-side
NPLib/Detectors/PISTA/TTimeData.cxx
+ 4
− 0
View file @ 5dda4eef
......@@ -53,6 +53,10 @@ void TTimeData::Clear() {
fTime_MWPC14.clear();
fTime_MWPC23.clear();
fTime_MWPC24.clear();
fSection_MWPC3.clear();
fSection_MWPC4.clear();
}
......
NPLib/Detectors/PISTA/TTimeData.h
+ 13
− 2
View file @ 5dda4eef
......@@ -44,6 +44,8 @@ class TTimeData : public TObject {
vector<float> fTime_MWPC23;
vector<float> fTime_MWPC24;
vector<int> fSection_MWPC3;
vector<int> fSection_MWPC4;
......@@ -80,6 +82,10 @@ class TTimeData : public TObject {
inline void SetTime_MWPC23(float time ){fTime_MWPC23.push_back(time);};//!
inline void SetTime_MWPC24(float time ){fTime_MWPC24.push_back(time);};//!
inline void SetSection_MWPC3(int section ){fSection_MWPC3.push_back(section);};//!
inline void SetSection_MWPC4(int section ){fSection_MWPC4.push_back(section);};//!
////////////////////// GETTERS ////////////////////////
......@@ -100,7 +106,12 @@ class TTimeData : public TObject {
{return fTime_MWPC23.at(i) ;}//!
inline float GetTime_MWPC24(const unsigned int &i) const
{return fTime_MWPC24.at(i) ;}//!
inline int GetSection_MWPC3(const unsigned int &i) const
{return fSection_MWPC3.at(i) ;}//!
inline int GetSection_MWPC4(const unsigned int &i) const
{return fSection_MWPC4.at(i) ;}//!
inline Int_t GetMWPC13Mult() const
{return static_cast<int>(fTime_MWPC13.size());}
inline Int_t GetMWPC14Mult() const
......@@ -112,7 +123,7 @@ class TTimeData : public TObject {
//////////////////////////////////////////////////////////////
// Required for ROOT dictionnary
ClassDef(TTimeData,1) // TimeData structure
ClassDef(TTimeData,2) // TimeData structure
};
#endif
Projects/AlPhaPha/2024/Analysis.cxx
+ 3
− 2
View file @ 5dda4eef
......@@ -505,10 +505,10 @@ void Analysis::VamosAnalysis(){
// To know time of electron drift in a section of the chio
for (int seg = 0 ; seg < IC->fIC_TS.size() ; seg++){
if (FPMW_Section != 0){
FF_DriftTime.push_back(10* (IC->fIC_TS.at(seg) - Time->GetTS_MWPC13(0) - ((Time->GetTime_MWPC14(0)-Toff14.at(FPMW_Section-1)) - (Time->GetTime_MWPC13(0)- Toff13.at(FPMW_Section-1)))));
FF_DriftTime.push_back(10* (IC->fIC_TS.at(seg) - Time->GetTS_MWPC13(0)) - ((Time->GetTime_MWPC13(0)+Toff13.at(FPMW_Section-1))));
}
else {
FF_DriftTime.push_back(10* (IC->fIC_TS.at(seg) - Time->GetTS_MWPC13(0) - ((Time->GetTime_MWPC14(0)) - (Time->GetTime_MWPC13(0)))));
FF_DriftTime.push_back(10* (IC->fIC_TS.at(seg) - Time->GetTS_MWPC13(0)) - ((Time->GetTime_MWPC13(0))));
}
}
......@@ -516,6 +516,7 @@ void Analysis::VamosAnalysis(){
FF_T13 = - MTOF_FP0_T0V[0] + Toff[FPMWPat];
//FF_T13 = -Time->GetTime_MWPC13(0) - Toff13[FPMW_Section-1] + 589;
FF_D13 = FF_Path - path1 + path2;
FF_V13 = FF_D13/FF_T13;
FF_Beta13 = FF_V13/29.9792458;
......
Projects/AlPhaPha/2024/ReadMe.md
+ 1
− 1
View file @ 5dda4eef
......@@ -39,7 +39,7 @@ Don't forget to modify the input file .
Your project.config must point to the root dir of your computer
# Compatibility AlmaLinux
Never import a detector class in a macro, if you do be sure to put a
Never import a detector class in a macro, if you do, be sure to put a
```
#ifdef 0
#include DETECTOR.h
......
Projects/AlPhaPha/2024/macro/mwpc/Toff/ToffGenerator.cxx
+ 49
− 7
View file @ 5dda4eef
#include "TICPhysics.h"
#include "TTimeData.h"
#include <TCanvas.h>
#include <TChain.h>
......@@ -27,7 +28,7 @@ void ToffGenerator(bool Create = true ){
UShort_t S13[20] , S14[20] , S23[20] , S24[20];
TTimeData *Time = new TTimeData();
TICPhysics *IC = new TICPhysics();
chain->SetBranchStatus("MTOF_FP0_T0VN", true);
chain->SetBranchAddress("MTOF_FP0_T0VN", &S13);
......@@ -44,6 +45,8 @@ void ToffGenerator(bool Create = true ){
chain->SetBranchStatus("Time", true);
chain->SetBranchAddress("Time", &Time);
chain->SetBranchStatus("IC", true);
chain->SetBranchAddress("IC", &IC);
//===========================================================================================================
// Histograms
//===========================================================================================================
......@@ -53,6 +56,12 @@ void ToffGenerator(bool Create = true ){
TH2F *hToff23;
TH2F *hToff24;
TH2F *hToffOnline13;
TH2F *hToffOnline14;
TH2F *hToffOnline23;
TH2F *hToffOnline24;
double xmin[4], xmax[4];
double ymin[4], ymax[4];
int binx = 20 , biny = 400;
......@@ -62,6 +71,9 @@ void ToffGenerator(bool Create = true ){
xmin[2] = 0 ; xmax[2] = 20 ; ymin[2] = -10 ; ymax[2] = 10 ;
xmin[3] = 0 ; xmax[3] = 20 ; ymin[3] = -10 ; ymax[3] = 10 ;
// Toff from aligning AoQ
double Toff[20] = {0, 588.0, 588.5, 587.95, 588.04, 587.72, 587.92, 587.9, 587.9, 588.66, 588.80, 588.67, 588.64, 588.75, 588.47, 588.65, 588.65, 588.67, 589.05, 590.3};
if (Create == true){
hToff13 = new TH2F("Toff13","Toff13",binx ,xmin[0] , xmax[0], biny, ymin[0], ymax[0]);
......@@ -69,6 +81,12 @@ void ToffGenerator(bool Create = true ){
hToff23 = new TH2F("Toff23","Toff23",binx ,xmin[2] , xmax[2], biny, ymin[2], ymax[2]);
hToff24 = new TH2F("Toff24","Toff24",binx ,xmin[3] , xmax[3], biny, ymin[3], ymax[3]);
hToffOnline13 = new TH2F("ToffOnline13","ToffOnline13",binx ,xmin[0] , xmax[0], biny, ymin[0], ymax[0]);
hToffOnline14 = new TH2F("ToffOnline14","ToffOnline14",binx ,xmin[1] , xmax[1], biny, ymin[1], ymax[1]);
hToffOnline23 = new TH2F("ToffOnline23","ToffOnline23",binx ,xmin[2] , xmax[2], biny, ymin[2], ymax[2]);
hToffOnline24 = new TH2F("ToffOnline24","ToffOnline24",binx ,xmin[3] , xmax[3], biny, ymin[3], ymax[3]);
//===========================================================================================================
// Loop on entries
//===========================================================================================================
......@@ -94,26 +112,29 @@ void ToffGenerator(bool Create = true ){
M14 = Time->GetMWPC14Mult();
M23 = Time->GetMWPC23Mult();
M24 = Time->GetMWPC24Mult();
if (M13 == 2 && (abs(S13[0]-S13[1]) == 1)){
vector<int> Pos = SortSection(S13);
hToff13->Fill(S13[0],Time->GetTime_MWPC13(Pos[0]) - Time->GetTime_MWPC13(Pos[1]));
hToffOnline13->Fill(S13[0],(Time->GetTime_MWPC13(Pos[0]) - Toff[S13[0]]) -( Time->GetTime_MWPC13(Pos[1]) -Toff[S13[1]] ) );
}
if (M23 == 2 && (abs(S23[0]-S23[1]) == 1)){
vector<int> Pos = SortSection(S23);
hToff23->Fill(S23[0],Time->GetTime_MWPC23(Pos[0]) - Time->GetTime_MWPC23(Pos[1]));
hToffOnline23->Fill(S23[0],(Time->GetTime_MWPC23(Pos[0]) - Toff[S23[0]]) -( Time->GetTime_MWPC23(Pos[1]) -Toff[S23[1]] ) );
}
if (M14 == 2 && (abs(S14[0]-S14[1]) == 1)){
vector<int> Pos = SortSection(S14);
hToff14->Fill(S14[0],Time->GetTime_MWPC14(Pos[0]) - Time->GetTime_MWPC14(Pos[1]));
hToffOnline14->Fill(S14[0],(Time->GetTime_MWPC14(Pos[0]) - Toff[S14[0]]) -( Time->GetTime_MWPC14(Pos[1]) -Toff[S14[1]] ) );
}
if (M24 == 2 && (abs(S24[0]-S24[1]) == 1)){
vector<int> Pos = SortSection(S24);
hToff24->Fill(S24[0],Time->GetTime_MWPC24(Pos[0]) - Time->GetTime_MWPC24(Pos[1]));
hToffOnline24->Fill(S24[0],(Time->GetTime_MWPC24(Pos[0]) - Toff[S24[0]]) -( Time->GetTime_MWPC24(Pos[1]) -Toff[S24[1]] ) );
}
} // End loop event
......@@ -122,6 +143,11 @@ void ToffGenerator(bool Create = true ){
hToff14->Write();
hToff23->Write();
hToff24->Write();
hToffOnline13->Write();
hToffOnline14->Write();
hToffOnline23->Write();
hToffOnline24->Write();
}// end if create
else {
......@@ -131,6 +157,12 @@ void ToffGenerator(bool Create = true ){
hToff14 = (TH2F*)in->Get("Toff14");
hToff23 = (TH2F*)in->Get("Toff23");
hToff24 = (TH2F*)in->Get("Toff24");
hToffOnline13 = (TH2F*)in->Get("ToffOnline13");
hToffOnline14 = (TH2F*)in->Get("ToffOnline14");
hToffOnline23 = (TH2F*)in->Get("ToffOnline23");
hToffOnline24 = (TH2F*)in->Get("ToffOnline24");
}
int MinX = 2;
......@@ -419,39 +451,49 @@ void ToffGenerator(bool Create = true ){
//===========================================================================================================
TCanvas *c1 = new TCanvas("c1","c1");
c1->Divide(2);
c1->Divide(3);
c1->cd(1);
hToff13->Draw("colz");
c1->cd(2);
hCorrToff13->Draw("colz");
c1->cd(3);
hToffOnline13->Draw("colz");
TCanvas *c2 = new TCanvas("c2","c2");
c2->Divide(2);
c2->Divide(3);
c2->cd(1);
hToff23->Draw("colz");
c2->cd(2);
hCorrToff23->Draw("colz");
c2->cd(3);
hToffOnline23->Draw("colz");
TCanvas *c3 = new TCanvas("c3","c3");
c3->Divide(2);
c3->Divide(3);
c3->cd(1);
hToff14->Draw("colz");
c3->cd(2);
hCorrToff14->Draw("colz");
c3->cd(3);
hToffOnline14->Draw("colz");
TCanvas *c4 = new TCanvas("c4","c4");
c4->Divide(2);
c4->Divide(3);
c4->cd(1);
hToff24->Draw("colz");
c4->cd(2);
hCorrToff24->Draw("colz");
c4->cd(3);
hToffOnline24->Draw("colz");
}
vector<int> SortSection(UShort_t STargFocal[2]){
......
Projects/AlPhaPha/conversion/2024/convert.C
+ 2
− 0
View file @ 5dda4eef
......@@ -86,6 +86,7 @@ void convert(int run=204){
for (int i=0; i<MTOF_FP0_T0VM;i++){
m_time->SetTS_MWPC13(MTOF_FP0_T0VTS[i]);
m_time->SetTime_MWPC13(MTOF_FP0_T0V[i]);
m_time->SetSection_MWPC3(MTOF_FP0_T0VN[i]);
}
for (int i=0; i<MTOF_FP0_T1VM;i++){
m_time->SetTS_MWPC14(MTOF_FP0_T1VTS[i]);
......@@ -95,6 +96,7 @@ void convert(int run=204){
for (int i=0; i<MTOF_FP1_T0VM;i++){
m_time->SetTS_MWPC23(MTOF_FP1_T0VTS[i]);
m_time->SetTime_MWPC23(MTOF_FP1_T0V[i]);
m_time->SetSection_MWPC3(MTOF_FP1_T0VN[i]);
}
for (int i=0; i<MTOF_FP1_T1VM;i++){
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment