-
Notifications
You must be signed in to change notification settings - Fork 1
Commit
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
Implemented Veto Wall position and energy calibration
- Loading branch information
Daniele DellAquila
committed
Jan 25, 2019
1 parent
8db8fe5
commit d8e18dc
Showing
166 changed files
with
18,325 additions
and
18,132 deletions.
There are no files selected for viewing
Large diffs are not rendered by default.
Oops, something went wrong.
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -1,181 +1,181 @@ | ||
| #include "include/E15190Reader.h" | ||
|
|
||
| //____________________________________________________ | ||
| void E15190Reader::CreateUsefulVetoWallHistograms(const char * file_name, Long64_t evt_amount) | ||
| { | ||
| // this method constructs some useful VetoWall histograms | ||
|
|
||
| if(!fIsNWB || !fIsFA || !fIsVW) return; | ||
|
|
||
| TFile * FileOut = new TFile(file_name, "RECREATE"); | ||
| if(FileOut->IsZombie()) { | ||
| printf("Error: Cannot create file %s\n", file_name); | ||
| return; | ||
| } | ||
|
|
||
| // Initializing output histograms | ||
| TH2D * VWDETOF[NUM_BARS_VW]; | ||
| TH2D * VWDETOFMatched[NUM_BARS_VW]; | ||
| for(int i=0; i<NUM_BARS_VW; i++) { | ||
| VWDETOF[i] = new TH2D (Form("VWDETOF%02d",i), Form("VWDETOF%02d",i), 600,0,300,1024,0,4096); | ||
| VWDETOFMatched[i] = new TH2D (Form("VWDETOF%02dMatched",i), Form("VWDETOF%02dMatched",i), 600,0,300,1024,0,4096); | ||
| } | ||
|
|
||
| Long64_t nentries=fChain->GetEntries(); | ||
| if(evt_amount!=0) { | ||
| nentries=evt_amount; | ||
| } | ||
| Long64_t jentry=0; | ||
| std::cout << "found " << nentries << " entries\n"; | ||
| fStartTime = clock(); // initializing start time | ||
| for(;fE15190Reader->Next() && jentry<nentries; jentry++) | ||
| { | ||
| if(jentry%100000==0) { | ||
| PrintPercentage(jentry,nentries); | ||
| } | ||
|
|
||
| NeutronWallCalibratedData * NWB = fNWBCal->Get(); | ||
| ForwardArrayCalibratedData * ForwardArray = fForwardArrayCal->Get(); | ||
| VetoWallCalibratedData * VetoWall = fVetoWallCal->Get(); | ||
|
|
||
| for(int i=0; i<NWB->fmulti; i++) { | ||
| double TOF = NWB->fTimeMean[i]-ForwardArray->fTimeMin; | ||
| for(int j=0; j<VetoWall->fmulti; j++) { | ||
| VWDETOF[VetoWall->fnumbar[j]]->Fill(TOF,VetoWall->fGeoMean[j]); | ||
| VWDETOFMatched[VetoWall->fnumbar[j]]->Fill(TOF,VetoWall->fGeoMeanMatched[j]); | ||
| } | ||
| } | ||
| } | ||
|
|
||
| //Writing histograms to file | ||
| for(int i=0; i<NUM_BARS_VW; i++) { | ||
| FileOut->WriteTObject(VWDETOF[i], VWDETOF[i]->GetName()); | ||
| FileOut->WriteTObject(VWDETOFMatched[i], VWDETOFMatched[i]->GetName()); | ||
| } | ||
|
|
||
| // closing output file | ||
| FileOut->Close(); | ||
| } | ||
|
|
||
| //____________________________________________________ | ||
| void E15190Reader::CreateUsefulForwardArrayHistograms(const char * file_name, Long64_t evt_amount) | ||
| { | ||
| // this method constructs some useful VetoWall histograms | ||
|
|
||
| if(!fIsNWB || !fIsFA) return; | ||
|
|
||
| TFile * FileOut = new TFile(file_name, "RECREATE"); | ||
| if(FileOut->IsZombie()) { | ||
| printf("Error: Cannot create file %s\n", file_name); | ||
| return; | ||
| } | ||
|
|
||
| // Initializing output histograms | ||
| TH1D * FANWBToF[NUM_DETECTORS_FA]; | ||
| for(int i=0; i<NUM_DETECTORS_FA; i++) { | ||
| FANWBToF[i] = new TH1D(Form("FA%02dNWBToF", i+1),Form("FA%02dNWBToF", i+1),600,-150,150); | ||
| } | ||
|
|
||
| Long64_t nentries=fChain->GetEntries(); | ||
| if(evt_amount!=0) { | ||
| nentries=evt_amount; | ||
| } | ||
| Long64_t jentry=0; | ||
| std::cout << "found " << nentries << " entries\n"; | ||
| fStartTime = clock(); // initializing start time | ||
| for(;fE15190Reader->Next() && jentry<nentries; jentry++) | ||
| { | ||
| if(jentry%100000==0) { | ||
| PrintPercentage(jentry,nentries); | ||
| } | ||
|
|
||
| NeutronWallCalibratedData * NWB = fNWBCal->Get(); | ||
| ForwardArrayCalibratedData * ForwardArray = fForwardArrayCal->Get(); | ||
|
|
||
| for(int FAhit=0; FAhit<ForwardArray->fmulti; FAhit++) | ||
| { | ||
| for(int NWhit=0; NWhit<NWB->fmulti; NWhit++) { | ||
| // I record the ToF spectrum only if the FA element has a significant energy (more than 300) | ||
| if(ForwardArray->fE[FAhit]>300) FANWBToF[ForwardArray->fnumdet[FAhit]-1]->Fill(NWB->fTimeMean[NWhit]-ForwardArray->fTime[FAhit]); | ||
| } | ||
| } | ||
|
|
||
| } | ||
|
|
||
| //Writing histograms to file | ||
| for(int i=0; i<NUM_DETECTORS_FA; i++) { | ||
| FileOut->WriteTObject(FANWBToF[i], FANWBToF[i]->GetName()); | ||
| } | ||
|
|
||
| // closing output file | ||
| FileOut->Close(); | ||
| } | ||
|
|
||
|
|
||
| //____________________________________________________ | ||
| void E15190Reader::CreateUsefulHiRAHistograms(const char * file_name, Long64_t evt_amount) | ||
| { | ||
| // this method constructs some useful HiRA histograms | ||
|
|
||
| if(!fIsHiRA) return; | ||
|
|
||
| TFile * FileOut = new TFile(file_name, "RECREATE"); | ||
| if(FileOut->IsZombie()) { | ||
| printf("Error: Cannot create file %s\n", file_name); | ||
| return; | ||
| } | ||
|
|
||
| // Initializing output histograms | ||
| //TH2D * HiRADERAWCSIRAW_SINGLESTRIPCSI=new TH2D("DE_Raw_ERaw_00_05_00", "DE_Raw_ERaw_00_05_00", 3000,0,3000,4000,0,8000); | ||
| TH2D * HiRADECALCSIRAW[NUM_CSI_TEL*NUM_TEL]; | ||
| TH2D * HiRADECALCSICALPROTONS[NUM_TEL]; | ||
| TH2D * HiRADECALCSICALALLPROTONS = new TH2D(Form("DECal_ECal_Protons"),Form("DECal_ECal_Protons"),2000,0,900,1000,0,200); | ||
| TH2D * HiRADECALCSICALALL = new TH2D(Form("DECal_ECal"),Form("DECal_ECal"),2000,0,500,900,0,200); | ||
| for(int i=0; i<NUM_TEL; i++) { | ||
| HiRADECALCSICALPROTONS[i] = new TH2D(Form("DECal_ECal_Protons_%02d", i),Form("DECal_ECal_Protons_%02d", i),2000,0,500,1000,0,200); | ||
| for(int j=0; j<NUM_CSI_TEL; j++) { | ||
| HiRADECALCSIRAW[NUM_CSI_TEL*i+j] = new TH2D(Form("DECal_ERaw_%02d_%02d", i, j ),Form("DECal_ERaw_%02d_%02d", i, j ),3000,0,3000,3000,0,200); | ||
| } | ||
| } | ||
|
|
||
| Long64_t nentries=fChain->GetEntries(); | ||
| if(evt_amount!=0) { | ||
| nentries=evt_amount; | ||
| } | ||
| Long64_t jentry=0; | ||
| std::cout << "found " << nentries << " entries\n"; | ||
| fStartTime = clock(); // initializing start time | ||
| for(;fE15190Reader->Next() && jentry<nentries; jentry++) | ||
| { | ||
| if(jentry%100000==0) { | ||
| PrintPercentage(jentry,nentries); | ||
| } | ||
|
|
||
| HiRACalibratedData * HiRA = fHiRACal->Get(); | ||
|
|
||
| for(int i=0; i< HiRA->fmulti; i++) { | ||
|
|
||
| //if(HiRA->fnumtel[i]==0 && HiRA->fnumcsi[i]==0 && HiRA->fnumstripf[i]==5) HiRADERAWCSIRAW_SINGLESTRIPCSI->Fill(HiRA->fEnergycsi[i], HiRA->fEnergySifHi[i]); | ||
|
|
||
| HiRADECALCSICALALLPROTONS->Fill(HiRA->fEnergycsiCalProtons[i], HiRA->fEnergySifCal[i]); | ||
| HiRADECALCSICALPROTONS[HiRA->fnumtel[i]]->Fill(HiRA->fEnergycsiCalProtons[i], HiRA->fEnergySifCal[i]); | ||
| HiRADECALCSICALALL->Fill(HiRA->fEnergycsiCal[i], HiRA->fEnergySifCal[i]); | ||
| HiRADECALCSIRAW[NUM_CSI_TEL*HiRA->fnumtel[i]+HiRA->fnumcsi[i]]->Fill(HiRA->fEnergycsi[i], HiRA->fEnergySifCal[i]); | ||
| } | ||
|
|
||
| } | ||
|
|
||
| //Writing histograms to file | ||
| //FileOut->WriteTObject(HiRADERAWCSIRAW_SINGLESTRIPCSI, HiRADERAWCSIRAW_SINGLESTRIPCSI->GetName()); | ||
| FileOut->WriteTObject(HiRADECALCSICALALLPROTONS, HiRADECALCSICALALLPROTONS->GetName()); | ||
| //FileOut->WriteTObject(HiRADECALCSICALALL, HiRADECALCSICALALL->GetName()); | ||
| for(int i=0; i<NUM_TEL; i++) { | ||
| //FileOut->WriteTObject(HiRADECALCSICALPROTONS[i], HiRADECALCSICALPROTONS[i]->GetName()); | ||
| for(int j=0; j<NUM_CSI_TEL; j++) { | ||
| FileOut->WriteTObject(HiRADECALCSIRAW[NUM_CSI_TEL*i+j], HiRADECALCSIRAW[NUM_CSI_TEL*i+j]->GetName()); | ||
| } | ||
| } | ||
|
|
||
| // closing output file | ||
| FileOut->Close(); | ||
| } | ||
| #include "include/E15190Reader.h" | ||
|
|
||
| //____________________________________________________ | ||
| void E15190Reader::CreateUsefulVetoWallHistograms(const char * file_name, Long64_t evt_amount) | ||
| { | ||
| // this method constructs some useful VetoWall histograms | ||
|
|
||
| if(!fIsNWB || !fIsFA || !fIsVW) return; | ||
|
|
||
| TFile * FileOut = new TFile(file_name, "RECREATE"); | ||
| if(FileOut->IsZombie()) { | ||
| printf("Error: Cannot create file %s\n", file_name); | ||
| return; | ||
| } | ||
|
|
||
| // Initializing output histograms | ||
| TH2D * VWDETOF[NUM_BARS_VW]; | ||
| TH2D * VWDETOFMatched[NUM_BARS_VW]; | ||
| for(int i=0; i<NUM_BARS_VW; i++) { | ||
| VWDETOF[i] = new TH2D (Form("VWDETOF%02d",i), Form("VWDETOF%02d",i), 600,0,300,1024,0,4096); | ||
| VWDETOFMatched[i] = new TH2D (Form("VWDETOF%02dMatched",i), Form("VWDETOF%02dMatched",i), 600,0,300,1024,0,4096); | ||
| } | ||
|
|
||
| Long64_t nentries=fChain->GetEntries(); | ||
| if(evt_amount!=0) { | ||
| nentries=evt_amount; | ||
| } | ||
| Long64_t jentry=0; | ||
| std::cout << "found " << nentries << " entries\n"; | ||
| fStartTime = clock(); // initializing start time | ||
| for(;fE15190Reader->Next() && jentry<nentries; jentry++) | ||
| { | ||
| if(jentry%100000==0) { | ||
| PrintPercentage(jentry,nentries); | ||
| } | ||
|
|
||
| NeutronWallCalibratedData * NWB = fNWBCal->Get(); | ||
| ForwardArrayCalibratedData * ForwardArray = fForwardArrayCal->Get(); | ||
| VetoWallCalibratedData * VetoWall = fVetoWallCal->Get(); | ||
|
|
||
| for(int i=0; i<NWB->fmulti; i++) { | ||
| double TOF = NWB->fTimeMean[i]-ForwardArray->fTimeMin; | ||
| for(int j=0; j<VetoWall->fmulti; j++) { | ||
| VWDETOF[VetoWall->fnumbar[j]]->Fill(TOF,VetoWall->fGeoMean[j]); | ||
| VWDETOFMatched[VetoWall->fnumbar[j]]->Fill(TOF,VetoWall->fGeoMeanMatched[j]); | ||
| } | ||
| } | ||
| } | ||
|
|
||
| //Writing histograms to file | ||
| for(int i=0; i<NUM_BARS_VW; i++) { | ||
| FileOut->WriteTObject(VWDETOF[i], VWDETOF[i]->GetName()); | ||
| FileOut->WriteTObject(VWDETOFMatched[i], VWDETOFMatched[i]->GetName()); | ||
| } | ||
|
|
||
| // closing output file | ||
| FileOut->Close(); | ||
| } | ||
|
|
||
| //____________________________________________________ | ||
| void E15190Reader::CreateUsefulForwardArrayHistograms(const char * file_name, Long64_t evt_amount) | ||
| { | ||
| // this method constructs some useful VetoWall histograms | ||
|
|
||
| if(!fIsNWB || !fIsFA) return; | ||
|
|
||
| TFile * FileOut = new TFile(file_name, "RECREATE"); | ||
| if(FileOut->IsZombie()) { | ||
| printf("Error: Cannot create file %s\n", file_name); | ||
| return; | ||
| } | ||
|
|
||
| // Initializing output histograms | ||
| TH1D * FANWBToF[NUM_DETECTORS_FA]; | ||
| for(int i=0; i<NUM_DETECTORS_FA; i++) { | ||
| FANWBToF[i] = new TH1D(Form("FA%02dNWBToF", i+1),Form("FA%02dNWBToF", i+1),600,-150,150); | ||
| } | ||
|
|
||
| Long64_t nentries=fChain->GetEntries(); | ||
| if(evt_amount!=0) { | ||
| nentries=evt_amount; | ||
| } | ||
| Long64_t jentry=0; | ||
| std::cout << "found " << nentries << " entries\n"; | ||
| fStartTime = clock(); // initializing start time | ||
| for(;fE15190Reader->Next() && jentry<nentries; jentry++) | ||
| { | ||
| if(jentry%100000==0) { | ||
| PrintPercentage(jentry,nentries); | ||
| } | ||
|
|
||
| NeutronWallCalibratedData * NWB = fNWBCal->Get(); | ||
| ForwardArrayCalibratedData * ForwardArray = fForwardArrayCal->Get(); | ||
|
|
||
| for(int FAhit=0; FAhit<ForwardArray->fmulti; FAhit++) | ||
| { | ||
| for(int NWhit=0; NWhit<NWB->fmulti; NWhit++) { | ||
| // I record the ToF spectrum only if the FA element has a significant energy (more than 300) | ||
| if(ForwardArray->fE[FAhit]>300) FANWBToF[ForwardArray->fnumdet[FAhit]-1]->Fill(NWB->fTimeMean[NWhit]-ForwardArray->fTime[FAhit]); | ||
| } | ||
| } | ||
|
|
||
| } | ||
|
|
||
| //Writing histograms to file | ||
| for(int i=0; i<NUM_DETECTORS_FA; i++) { | ||
| FileOut->WriteTObject(FANWBToF[i], FANWBToF[i]->GetName()); | ||
| } | ||
|
|
||
| // closing output file | ||
| FileOut->Close(); | ||
| } | ||
|
|
||
|
|
||
| //____________________________________________________ | ||
| void E15190Reader::CreateUsefulHiRAHistograms(const char * file_name, Long64_t evt_amount) | ||
| { | ||
| // this method constructs some useful HiRA histograms | ||
|
|
||
| if(!fIsHiRA) return; | ||
|
|
||
| TFile * FileOut = new TFile(file_name, "RECREATE"); | ||
| if(FileOut->IsZombie()) { | ||
| printf("Error: Cannot create file %s\n", file_name); | ||
| return; | ||
| } | ||
|
|
||
| // Initializing output histograms | ||
| //TH2D * HiRADERAWCSIRAW_SINGLESTRIPCSI=new TH2D("DE_Raw_ERaw_00_05_00", "DE_Raw_ERaw_00_05_00", 3000,0,3000,4000,0,8000); | ||
| TH2D * HiRADECALCSIRAW[NUM_CSI_TEL*NUM_TEL]; | ||
| TH2D * HiRADECALCSICALPROTONS[NUM_TEL]; | ||
| TH2D * HiRADECALCSICALALLPROTONS = new TH2D(Form("DECal_ECal_Protons"),Form("DECal_ECal_Protons"),2000,0,900,1000,0,200); | ||
| TH2D * HiRADECALCSICALALL = new TH2D(Form("DECal_ECal"),Form("DECal_ECal"),2000,0,500,900,0,200); | ||
| for(int i=0; i<NUM_TEL; i++) { | ||
| HiRADECALCSICALPROTONS[i] = new TH2D(Form("DECal_ECal_Protons_%02d", i),Form("DECal_ECal_Protons_%02d", i),2000,0,500,1000,0,200); | ||
| for(int j=0; j<NUM_CSI_TEL; j++) { | ||
| HiRADECALCSIRAW[NUM_CSI_TEL*i+j] = new TH2D(Form("DECal_ERaw_%02d_%02d", i, j ),Form("DECal_ERaw_%02d_%02d", i, j ),3000,0,3000,3000,0,200); | ||
| } | ||
| } | ||
|
|
||
| Long64_t nentries=fChain->GetEntries(); | ||
| if(evt_amount!=0) { | ||
| nentries=evt_amount; | ||
| } | ||
| Long64_t jentry=0; | ||
| std::cout << "found " << nentries << " entries\n"; | ||
| fStartTime = clock(); // initializing start time | ||
| for(;fE15190Reader->Next() && jentry<nentries; jentry++) | ||
| { | ||
| if(jentry%100000==0) { | ||
| PrintPercentage(jentry,nentries); | ||
| } | ||
|
|
||
| HiRACalibratedData * HiRA = fHiRACal->Get(); | ||
|
|
||
| for(int i=0; i< HiRA->fmulti; i++) { | ||
|
|
||
| //if(HiRA->fnumtel[i]==0 && HiRA->fnumcsi[i]==0 && HiRA->fnumstripf[i]==5) HiRADERAWCSIRAW_SINGLESTRIPCSI->Fill(HiRA->fEnergycsi[i], HiRA->fEnergySifHi[i]); | ||
|
|
||
| HiRADECALCSICALALLPROTONS->Fill(HiRA->fEnergycsiCalProtons[i], HiRA->fEnergySifCal[i]); | ||
| HiRADECALCSICALPROTONS[HiRA->fnumtel[i]]->Fill(HiRA->fEnergycsiCalProtons[i], HiRA->fEnergySifCal[i]); | ||
| HiRADECALCSICALALL->Fill(HiRA->fEnergycsiCal[i], HiRA->fEnergySifCal[i]); | ||
| HiRADECALCSIRAW[NUM_CSI_TEL*HiRA->fnumtel[i]+HiRA->fnumcsi[i]]->Fill(HiRA->fEnergycsi[i], HiRA->fEnergySifCal[i]); | ||
| } | ||
|
|
||
| } | ||
|
|
||
| //Writing histograms to file | ||
| //FileOut->WriteTObject(HiRADERAWCSIRAW_SINGLESTRIPCSI, HiRADERAWCSIRAW_SINGLESTRIPCSI->GetName()); | ||
| FileOut->WriteTObject(HiRADECALCSICALALLPROTONS, HiRADECALCSICALALLPROTONS->GetName()); | ||
| //FileOut->WriteTObject(HiRADECALCSICALALL, HiRADECALCSICALALL->GetName()); | ||
| for(int i=0; i<NUM_TEL; i++) { | ||
| //FileOut->WriteTObject(HiRADECALCSICALPROTONS[i], HiRADECALCSICALPROTONS[i]->GetName()); | ||
| for(int j=0; j<NUM_CSI_TEL; j++) { | ||
| FileOut->WriteTObject(HiRADECALCSIRAW[NUM_CSI_TEL*i+j], HiRADECALCSIRAW[NUM_CSI_TEL*i+j]->GetName()); | ||
| } | ||
| } | ||
|
|
||
| // closing output file | ||
| FileOut->Close(); | ||
| } |
Oops, something went wrong.