Closed jtenavidal closed 11 months ago
TF1 *myElectronFit = new TF1("myElectronFit","[0]+[1]/x",0.,5.);
if (en_beam[fbeam_en] == 1.161) { myElectronFit->SetParameters(17,7); }
if (en_beam[fbeam_en] == 2.261) { myElectronFit->SetParameters(16,10.5); }
if (en_beam[fbeam_en] == 4.461) { myElectronFit->SetParameters(13.5,15); }
...
double theta_min = myElectronFit->Eval(el_momentum); // in deg
if (el_theta*180./TMath::Pi() < theta_min) { continue; }
// Explicit cuts on electron momentum
if (fbeam_en=="1161" && el_momentum < 0.4) { continue; }
if (fbeam_en=="2261" && el_momentum < 0.55) { continue; }
if (fbeam_en=="4461" && el_momentum < 1.1) { continue; }
For pi -:
Bool_t Fiducial::PimiFiducialCutExtra(std::string beam_en, TVector3 momentum) {
} For electron, can I use the electron fit?
E truth cuts
11:45 https://github.com/adishka/e4nu/tree/e4v_truth 11:45 in the readme.txt file