-
hallc-online authored
1) Add tree variables for the three BPMs and projections to target in the EPICs coordinate system (+X beam right, +Y up) rb.raster.fr_xbpm_tar rb.raster.fr_ybpm_tar rb.raster.fr_xbpmA rb.raster.fr_ybpmA rb.raster.fr_xbpmB rb.raster.fr_ybpmB rb.raster.fr_xbpmC rb.raster.fr_ybpmC 2) Modified code so that it only modifies the BPM variables if they are all !=0 ( beam on) . Otherwise it sets them to the previous EPICS read. This gets rid of all zero EPICS reads,except if the first EPICs read the beam is off. 3) Removed the reading of gpbeam from the global parameter list in the Process method and just have it in ReadDatabase. 4) Change the calculation of beam at target to use the projection of BPM_A and BPM_C instead of BPM_A and BPM_B. 5) Remove setting the raster position to the hcana beam coordinate system with (+X beam left). The raster is in the EPICS coordinate system. This makes it easier to interpret the carbon hole runs.hallc-online authored1) Add tree variables for the three BPMs and projections to target in the EPICs coordinate system (+X beam right, +Y up) rb.raster.fr_xbpm_tar rb.raster.fr_ybpm_tar rb.raster.fr_xbpmA rb.raster.fr_ybpmA rb.raster.fr_xbpmB rb.raster.fr_ybpmB rb.raster.fr_xbpmC rb.raster.fr_ybpmC 2) Modified code so that it only modifies the BPM variables if they are all !=0 ( beam on) . Otherwise it sets them to the previous EPICS read. This gets rid of all zero EPICS reads,except if the first EPICs read the beam is off. 3) Removed the reading of gpbeam from the global parameter list in the Process method and just have it in ReadDatabase. 4) Change the calculation of beam at target to use the projection of BPM_A and BPM_C instead of BPM_A and BPM_B. 5) Remove setting the raster position to the hcana beam coordinate system with (+X beam left). The raster is in the EPICS coordinate system. This makes it easier to interpret the carbon hole runs.
THcRaster.h 5.33 KiB
//author Burcu Duran - Melanie Rehfuss (2017)
#ifndef ROOT_THcRaster
#define ROOT_THcRaster
///////////////////////////////////////////////////////////////////////////////
// //
// THcRaster //
// //
///////////////////////////////////////////////////////////////////////////////
#include "THaBeamDet.h"
#include "TVector.h"
#include "TClonesArray.h"
#include "THcHitList.h"
#include "THcDetectorMap.h"
#include "THcRasterRawHit.h"
#include "THaCutList.h"
#include "THaOutput.h"
#include "THaEpicsEvtHandler.h"
class THcRaster : public THaBeamDet, public THcHitList {
public:
THcRaster(const char* name, const char* description ="",THaApparatus* a = NULL );
~THcRaster();
void Clear(Option_t* opt="");
void AccumulatePedestals(TClonesArray* rawhits);
void CalculatePedestals();
Int_t Decode( const THaEvData& );
Int_t ReadDatabase( const TDatime& date );
Int_t DefineVariables( EMode mode );
EStatus Init( const TDatime& run_time );
void InitArrays() { /* do nothing */;}
void DeleteArrays() { /* do nothing */;}
Int_t Process();
TVector3 GetPosition() const { return fPosition[1]; }
TVector3 GetDirection() const { return fDirection; }
Double_t GetCurrentX() { return fRawPos[0]; }
Double_t GetCurrentY() { return fRawPos[1]; }
THcRaster();
protected:
//event information
Int_t fNhits;
/* void InitializeReconstruction(); */
Double_t fgpbeam; //beam momentum
Double_t fgfrx_dist; //Distance of raster to target
Double_t fgfry_dist;
Double_t fgbeam_xoff; // Beam offsets
Double_t fgbeam_xpoff; // Beam offsets
Double_t fgbeam_yoff; // Beam offsets
Double_t fgbeam_ypoff; // Beam offsets
Double_t fgbpmxa_slope; // BPM slope/offsets
Double_t fgbpmxa_off; // BPM slope/offsets Double_t fgbpmxb_slope; // BPM slope/offsets
Double_t fgbpmxb_off; // BPM slope/offsets
Double_t fgbpmxc_slope; // BPM slope/offsets
Double_t fgbpmxc_off; // BPM slope/offsets
Double_t fgbpmya_slope; // BPM slope/offsets
Double_t fgbpmya_off; // BPM slope/offsets
Double_t fgbpmyb_slope; // BPM slope/offsets
Double_t fgbpmyb_off; // BPM slope/offsets
Double_t fgbpmyc_slope; // BPM slope/offsets
Double_t fgbpmyc_off; // BPM slope/offsets
Double_t fgbpma_zpos; // BPM A z position
Double_t fgbpmb_zpos; // BPM B z position
Double_t fgbpmc_zpos; // BPM C x position
Int_t fgusefr; /* Use Raster for beam position */
Double_t FRXA_rawadc; // XA raw ADC
Double_t FRYA_rawadc; // YA raw ADC
Double_t FRXB_rawadc; // XB raw ADC
Double_t FRYB_rawadc; // YB raw ADC
Double_t BPMXA_raw; // BPM A Raw X Position
Double_t BPMYA_raw; // BPM A Raw Y Position
Double_t BPMXB_raw; // BPM B Raw X Position
Double_t BPMYB_raw; // BPM B Raw Y Position
Double_t BPMXC_raw; // BPM C Raw X Position
Double_t BPMYC_raw; // BPM C Raw Y Position
Double_t BPMXA_pos; // BPM A Corrected X Position
Double_t BPMYA_pos; // BPM A Corrected Y Position
Double_t BPMXB_pos; // BPM B Corrected X Position
Double_t BPMYB_pos; // BPM B Corrected Y Position
Double_t BPMXC_pos; // BPM C Corrected X Position
Double_t BPMYC_pos; // BPM C Corrected Y Position
Double_t fXA_ADC; // XA ADC
Double_t fYA_ADC; // YA ADC
Double_t fXB_ADC; // XB ADC
Double_t fYB_ADC; // YB ADC
Double_t fXA_pos; // XA position
Double_t fYA_pos; // YA position
Double_t fXB_pos; // XB position
Double_t fYB_pos; // YB position
Double_t fXbpm_tar; // X BPM at target (+X is beam right)
Double_t fYbpm_tar; // Y BPM at target (+Y is up)
Double_t fXbpm_A; // X BPM at BPMA (+X is beam right)
Double_t fYbpm_A; // Y BPM at BPMA (+Y is up)
Double_t fXbpm_B; // X BPM at BPMB (+X is beam right)
Double_t fYbpm_B; // Y BPM at BPMB (+Y is up)
Double_t fXbpm_C; // X BPM at BPMC (+X is beam right)
Double_t fYbpm_C; // Y BPM at BPMC (+Y is up)
Double_t fXbeam_prev[4]; //
Double_t fYbeam_prev[4]; //
Double_t fFrXA_ADC_zero_offset;
Double_t fFrYA_ADC_zero_offset;
Double_t fFrXB_ADC_zero_offset;
Double_t fFrYB_ADC_zero_offset;
Double_t fPedADC[4]; // ADC poedestals
//Double_t fAvgPedADC[4]; // Avergage ADC poedestals
Double_t fRawPos[2]; // current in Raster ADCs for position
TVector3 fPosition[3]; // Beam position at 1st, 2nd BPM or at the target (meters)
TVector3 fDirection;
TClonesArray* frPosAdcPulseIntRaw;
THaEpicsEvtHandler* fEpicsHandler;
private:
Bool_t fAnalyzePedestals;
Int_t fNPedestalEvents;
Double_t fFrCalMom; Double_t fFrXA_ADCperCM;
Double_t fFrYA_ADCperCM;
Double_t fFrXB_ADCperCM;
Double_t fFrYB_ADCperCM;
ClassDef(THcRaster, 0); // add THcRaster to ROOT library
};
////////////////////////////////////////////////////////////////////////////////
#endif