BmnZDCEventData.cxx

Go to the documentation of this file.

/* 
 * File:   BmnZDCEventData.cxx
 * Author: pnaleks
 * 
 * Created on 5 июня 2020 г., 13:43
 */
 
#include "BmnZDCEventData.h"
 
#include <iostream>
using namespace std;
 
BmnZDCEventData::BmnZDCEventData() {
    fE = fEc = fEn = fEp = 0.;
    fH = fHc = fHn = fHp = 0;
    fX = fY = 0.;
    fAsymmetry = fMoment = 0.;
}
 
BmnZDCEventData::~BmnZDCEventData() {
}
 
void BmnZDCEventData::Set(TClonesArray * pArrayOfBmnZDCDigits, Float_t * pModuleScale, Float_t * pModuleThreshold) {
    fE = fEc = fEn = fEp = 0.;
    fH = fHc = fHn = fHp = 0;
    fX = fY = 0.;
    fAsymmetry = fMoment = 0.;
    
    Double_t SEx = 0, SEy = 0, SExx = 0, SEyy = 0;
    Int_t N = pArrayOfBmnZDCDigits->GetEntries();
    for (Int_t i = 0; i < N; i++) {
        BmnZDCDigit * p = (BmnZDCDigit *) pArrayOfBmnZDCDigits->At(i);
        
        Int_t ch = p->GetChannel();
        if (ch < 1 || ch > 104) continue;
        
        Float_t e = p->GetAmp() * 1e-3; // MeV -> GeV
        if (pModuleScale) e *= pModuleScale[ch];
        
        if (pModuleThreshold && e < pModuleThreshold[ch]) continue;
        
        fH++;
        fE += e;
        
        switch (ch) {
            case  3: case  4: case  5:
            case 10: case 11: case 12:
            case 17: case 18: case 19:
            case 24: case 25: case 26:
                fHn++;
                fEn += e;
                break;
            case 43: case 44: case 45:
            case 50: case 51: case 52:
            case 57: case 58: case 59:
            case 64: case 65: case 66:
                fHp++;
                fEp += e;
                break;
            default:
                if (ch > 68) {
                    fHc++;
                    fEc += e;
                }
        }
        
        SEx += e * p->GetX();
        SEy += e * p->GetY();
        SExx += e * p->GetX() * p->GetX();
        SEyy += e * p->GetY() * p->GetY();
    }
    
    if (fE > 0.) {
        fX = SEx / fE / 10.; // mm -> cm
        fY = SEy / fE / 10.; // mm -> cm
        fMoment = (SExx + SEyy) / fE / 100.; // mm^2 -> cm^2
        fAsymmetry = 2.*fEc/fE - 1.;
    }
}
 
Float_t BmnZDCEventData::GetEnergy(char op) {
    switch (op) {
        case 'c': case 'C': return fEc;
        case 'n': case 'N': return fEn;
        case 'p': case 'P': return fEp;
    }
    
    return fE;
}
 
Int_t BmnZDCEventData::GetNHits(char op) {
    switch (op) {
        case 'c': case 'C': return fHc;
        case 'n': case 'N': return fHn;
        case 'p': case 'P': return fHp;
    }
    
    return fH;
}
 
void BmnZDCEventData::Print(Option_t *option) const {
    cout << "BmnZDCEventData:"; 
    if (fH) {
        cout << endl;
        cout << "  Hits:   " << fH << endl;
        cout << "  Energy: " << fE << " GeV" << endl;
        cout << "  Asym.:  " << fAsymmetry << endl;
        cout << "  Moment: " << fMoment << " cm^2" << endl;
    } else {
        cout << " is empty" << endl;
    }
}