Fiducial variation - Githubissues

jtenavidal commented 4 months ago

Bool_t GetEPhiLimits(double beam_en, Float_t momentum, Float_t theta, Int_t sector,Float_t *EPhiMin, Float_t *EPhiMax);
    Bool_t EFiducialCut(double beam_en, TVector3 momentum);
    Bool_t PFiducialCut(double beam_en, TVector3 momentum);
    Bool_t PiplFiducialCut(double beam_en, TVector3 momentum,Float_t *philow,Float_t *phiup);
    Bool_t PimiFiducialCut(double beam_en, TVector3 momentum, Float_t *pimi_philow, Float_t *pimi_phiup);

PiplFiducialCut returns philow, and phiup. We can apply a variation there and see if it is still in fiducial

For electron we use GetEPhiLimits

We can add this as a new function with default 0 variation

jtenavidal commented 4 months ago

[x] Add on Configuration, Phi/Theta Shift for fiducial. It has to get a number from it ![Uploading image.png…]()
[x] From the shift, need to add a function in Fiducial to check if it is in the fiducial after the shift
[x] Add shift info in output file name if it is not 0

jtenavidal commented 4 months ago

Angle shift was already included:

    } else if ( param[i] == "SystFidAngleShift" ) {
      fFidAngleShift = stod(value[i]) ;
    }

jtenavidal commented 4 months ago

bool GetFidAngleShift(void) { return fFidAngleShift; }

jtenavidal commented 4 months ago

EDIT THIS AND ADD IN CODE: `Bool_t Fiducial::GetFiducialCut(unsigned int pdg, double beam_en, TVector3 momentum, double shift = 0){

bool in_fid = false ; if( pdg == conf::kPdgElectron ) in_fid = Fiducial::EFiducialCut( beam_en, momentum); if( pdg == conf::kPdgProton ) in_fid = PFiducialCut( beam_en, momentum ) ; if( pdg == conf::kPdgPim ) in_fid = PimiFiducialCut( beam_en, momentum ) ; if( pdg == conf::kPdgPip ) in_fid return in_fid ; / Bool_t PFiducialCut(double beam_en, TVector3 momentum);
Bool_t PiplFiducialCut(double beam_en, TVector3 momentum,Float_t philow,Float_t phiup);
Bool_t PimiFiducialCut(double beam_en, TVector3 momentum, Float_t pimi_philow, Float_t pimi_phiup);
bool Phot_fid(TVector3 V3_phot);
bool Pi_phot_fid_united(double beam_en, TVector3 V3_pi_phot, int q_pi_phot);/ } `

jtenavidal commented 4 months ago

Configured correctly: