BeamSummaryAna.cxx

Go to the documentation of this file.

#include <iostream>
#include <cmath>

#include "BeamData/ana/Summary/BeamSummaryAna.h"
#include "MessageService/MsgService.h"

CVSID("$Id: BeamSummaryAna.cxx,v 1.4 2005/06/07 21:50:13 mdier Exp $");

using namespace std;

BeamSummaryAna::BeamSummaryAna()
{
    SetDefaults();
}

void BeamSummaryAna::SetDefaults()
{
    SetMinIntensCut(0.15); 
    SetMaxIntensCut(1000); 
    SetMinHorPosTgtCut(-3.0);
    SetMaxHorPosTgtCut(0.0);
    SetMinVerPosTgtCut(-1.0);
    SetMaxVerPosTgtCut(2.0);
    SetMinHorWidthCut(0.0);   
    SetMaxHorWidthCut(1.5);   
    SetMinVerWidthCut(0.0);   
    SetMaxVerWidthCut(1.5);   
    SetMinHornCurCut(190);  
    SetMaxHornCurCut(2000);  
    SelectTargetPosition("in");
}

void BeamSummaryAna::SetMinIntensCut(Double_t cutval){
    fIntensMin = cutval;
}
void BeamSummaryAna::SetMaxIntensCut(Double_t cutval){
    fIntensMax = cutval;
}
void BeamSummaryAna::SetMinHorPosTgtCut(Double_t cutval){
    fHorPosTgtMin = cutval;
}
void BeamSummaryAna::SetMaxHorPosTgtCut(Double_t cutval){
    fHorPosTgtMax = cutval;
}
void BeamSummaryAna::SetMinVerPosTgtCut(Double_t cutval){
    fVerPosTgtMin = cutval;
}
void BeamSummaryAna::SetMaxVerPosTgtCut(Double_t cutval){
    fVerPosTgtMax = cutval;
}
void BeamSummaryAna::SetMinHorWidthCut(Double_t cutval){
    fHorWidthMin = cutval;
}
void BeamSummaryAna::SetMaxHorWidthCut(Double_t cutval){
    fHorWidthMax = cutval;
}
void BeamSummaryAna::SetMinVerWidthCut(Double_t cutval){
    fVerWidthMin = cutval;
}
void BeamSummaryAna::SetMaxVerWidthCut(Double_t cutval){
    fVerWidthMax = cutval;
}
void BeamSummaryAna::SetMinHornCurCut(Double_t cutval){
    fHornCurMin = cutval;
}
void BeamSummaryAna::SetMaxHornCurCut(Double_t cutval){
    fHornCurMax = cutval;
}
void BeamSummaryAna::SetMinTargetPosCut(Double_t cutval){
    fTargetPosMin = cutval;
}
void BeamSummaryAna::SetMaxTargetPosCut(Double_t cutval){
    fTargetPosMax = cutval;
}

void BeamSummaryAna::SelectTargetPosition(const char* topt){
    string tchoice = topt;
    
    if (tchoice == "all") {
        SetMinTargetPosCut(-1e5);
        SetMaxTargetPosCut(1e6);
    } else if (tchoice == "out") {
        SetMinTargetPosCut(-1e5);
        SetMaxTargetPosCut(0.0);
    } else if (tchoice == "in") {
        SetMinTargetPosCut(0.0);
        SetMaxTargetPosCut(1e6);
    } else if (tchoice == "le") {
        SetMinTargetPosCut(0.0);
        SetMaxTargetPosCut(1000);
    } else if (tchoice == "ple") {
        SetMinTargetPosCut(3500);
        SetMaxTargetPosCut(4500);
    } else if (tchoice == "pme") {
        SetMinTargetPosCut(20000);
        SetMaxTargetPosCut(60000);
    } else if (tchoice == "phe") {
        SetMinTargetPosCut(60000);
        SetMaxTargetPosCut(1e6);
    } else {
        MSG("BeamSummaryAna",Msg::kWarning) << "Wrong option for SelectTargetPosition. "
                                            << "Will use default: target in" << endl;
        SetMinTargetPosCut(0.0);
        SetMaxTargetPosCut(1e6);
    }    
}

void BeamSummaryAna::PrintCuts()
{
   cout << "The cut values for beam spills are:" << endl;
   cout << "  " << fIntensMin << "E12 PoT < beam intensity < "
        << fIntensMax << "E12 PoT" << endl;
   cout << "  " << fHorPosTgtMin << "mm < horizontal position at target < "
        << fHorPosTgtMax << "mm " << endl;
   cout << "  " << fVerPosTgtMin << "mm < vertical position at target < "
        << fVerPosTgtMax << "mm " << endl; 
   cout << "     (note: target not at (0,0) in this coord system)" << endl; 
   cout << "  " << fHorWidthMin << "mm < horizontal width at target < "
        << fHorWidthMax << "mm " << endl;
   cout << "  " << fVerWidthMin << "mm < vertical width at target < "
        << fVerWidthMax << "mm " << endl;
   cout << "  " << fHornCurMin << "kA < absolute horn peak current < "
        << fHornCurMax << "kA " << endl;   
   cout << "  " << fTargetPosMin << "mils < target position < "
        << fTargetPosMax << "mils " << endl; 
}


Bool_t BeamSummaryAna::SimpleBeamCut(const BeamSummary& bsum)
{
    // Select spills with beam
    if (bsum.beamIntensity<fIntensMin ||
        bsum.beamIntensity>fIntensMax) return false;
    // Horizontal position
    if (bsum.hPosatTargetBPM<fHorPosTgtMin ||
        bsum.hPosatTargetBPM>fHorPosTgtMax) return false;
    // Vertical position 
    if (bsum.vPosatTargetBPM<fVerPosTgtMin ||
        bsum.vPosatTargetBPM>fVerPosTgtMax) return false;
    // Horizontal Beam width 
    if (bsum.hBeamWidth<fHorWidthMin ||
        bsum.hBeamWidth>fHorWidthMax) return false;
    // Vertical Beam width 
    if (bsum.vBeamWidth<fVerWidthMin ||
        bsum.vBeamWidth>fVerWidthMax) return false;
    // Horn current 
    Double_t hcur = bsum.hornPeakCurrent ;
    if ( fabs(hcur) < fHornCurMin ||
        fabs(hcur) > fHornCurMax) return false;
    // Target in
    if (bsum.nuTarZ<fTargetPosMin ||
        bsum.nuTarZ>fTargetPosMax) return false;
    
    return true;
}

---