hmdccutstat.cc

Go to the documentation of this file.

//*-- AUTHOR : Jochen Markert

//_HADES_CLASS_DESCRIPTION 
///////////////////////////////////////////////////////////////
//HMdcCutStat
//
//  class for statistics on cut time1, time2 and time2-time1
//  and wires less than 4 per module
//////////////////////////////////////////////////////////////

#include "hmdccutstat.h"
#include "hades.h"
#include "hruntimedb.h"
#include "hmdcparasciifileio.h"
#include "hpario.h"
#include "hmdcdetector.h"
#include "hspectrometer.h"

HMdcCutStat* HMdcCutStat::fMdcCutStat=0;
HMdcCutStat::HMdcCutStat(const Char_t* name,const Char_t* title,
                         const Char_t* context,Int_t s,Int_t m)
    : HParSet(name,title,context), fSecs(s)
{
    fMdcCutStat=this;
    for (Int_t i=0;i<s; i++) fSecs.AddAt(new HMdcCutStatSec(m),i);
    setNameTitle();
    if (gHades) {
        fMdc = (HMdcDetector*)(((HSpectrometer*)(gHades->getSetup()))->getDetector("Mdc"));
    } else {
        fMdc = 0;
    }
}
HMdcCutStat::~HMdcCutStat()
{}
HMdcCutStat* HMdcCutStat::getExObject()
{
    return fMdcCutStat;
}
HMdcCutStat* HMdcCutStat::getObject()
{
    if(!fMdcCutStat) fMdcCutStat=new HMdcCutStat();
    return fMdcCutStat;
}
void HMdcCutStat::deleteHMdcCutStat()
{
    if(fMdcCutStat) delete fMdcCutStat;
    fMdcCutStat=0;
}
void HMdcCutStat::clear(void)
{
    for (Int_t i=0;i<getSize();i++)
    {
        HMdcCutStatSec &sector = (*this)[i];
        for (Int_t m=0; m<sector.getSize(); m++)
        {
            HMdcCutStatMod &mod = sector[m];
            mod.clear();
        }
    }
}
void HMdcCutStat::setNameTitle()
{
    strcpy(detName,"Mdc");
}
void HMdcCutStat::getCal1StatCut(Int_t s,Int_t mod, Int_t *cT1, Int_t *cT2,
                                 Int_t *cT12, Int_t *cWire,Int_t* cT1T12,Int_t* cBump,Int_t* cAll,Int_t* all)
{   // returns values of different cut statistics to Int_t*
    // cT1    : cut in T1
    // cT2    : cut in T2
    // cT12   : cut in T2-T1
    // cWire  : cut in number of Wires <4
    // cT1T12 : cut in T1 and T2-T1
    // cBump  : cut in T1 and T2-T1
    // cAll   : all cuts
    // all    : all hits

    HMdcCutStatMod &m = (*this)[s][mod];
    *cT1    = m.getCal1StatCutT1();
    *cT2    = m.getCal1StatCutT2();
    *cT12   = m.getCal1StatCutT12();
    *cWire  = m.getCal1StatCutWire();
    *cT1T12 = m.getCal1StatCutT1T12();
    *cBump  = m.getCal1StatCutBump();
    *cAll   = m.getCal1StatCutAll();
    *all    = m.getCal1StatAll();
}
void HMdcCutStat::setCal1StatCut(Int_t s,Int_t mod, Int_t cT1, Int_t cT2,
                                 Int_t cT12,Int_t cBump, Int_t cWire)
{
    // counts up values of different cut statistics
    // cT1    : cut in T1
    // cT2    : cut in T2
    // cT12   : cut in T2-T1
    // cWire  : cut in number of Wires <4
    // cT1T12 : cut in T1 and T2-T1 internaly filled
    // cBump  : cut in T1 and T2-T1 internaly filled
    // cAll   : all cuts            internaly filled
    // all    : all hits            internaly filled

    HMdcCutStatMod &m = (*this)[s][mod];
    if(cT1==1)    m.setCal1StatCutT1();
    if(cT2==1)    m.setCal1StatCutT2();
    if(cT12==1)   m.setCal1StatCutT12();
    if(cWire==1)  m.setCal1StatCutWire();
    if(cT1==1  ||
       cT12==1)   m.setCal1StatCutT1T12();
    if(cBump==1)  m.setCal1StatCutBump();
    if(cT1==1  ||
       cT2==1  ||
       cT12==1 ||
       cBump==1 ||
       cWire==1)  m.setCal1StatCutAll();
                  m.setCal1StatAll();
}
Bool_t HMdcCutStat::init(HParIo* inp,Int_t* set) {
    // intitializes the container from an input
    HDetParIo* input=inp->getDetParIo("HMdcParIo");
    if (input) return (input->init(this,set));
    return kFALSE;
}
Int_t HMdcCutStat::write(HParIo* output) {
    // writes the container to an output
    HDetParIo* out=output->getDetParIo("HMdcParIo");
    if (out) return out->write(this);
    return -1;
}
void HMdcCutStat::putAsciiHeader(TString& b) {
    b =
        "#######################################################################\n"
        "# Statistic onDrift time cuts for noise reduction of the MDC\n"
        "# Format:\n"
        "# sector mod ct1 ct2 ct12 cwires ct1t12 cbump call all \n"
        "#######################################################################\n";
}
Bool_t HMdcCutStat::writeline(Char_t*buf, Int_t s, Int_t m) {
    Bool_t r = kTRUE;
    if (fMdc) {
        if (fMdc->getModule(s,m) != 0) {
            if  (s>-1 && s<getSize()) {
                HMdcCutStatSec &sector = (*this)[s];
                if (m>-1 && m<sector.getSize()) {
                    sprintf(buf,"%i %i %i %i %i %i %i %i %i %i \n",
                            s,m,
                            (*this)[s][m].getCal1StatCutT1(),
                            (*this)[s][m].getCal1StatCutT2(),
                            (*this)[s][m].getCal1StatCutT12(),
                            (*this)[s][m].getCal1StatCutWire(),
                            (*this)[s][m].getCal1StatCutT1T12(),
                            (*this)[s][m].getCal1StatCutBump(),
                            (*this)[s][m].getCal1StatCutAll(),
                            (*this)[s][m].getCal1StatAll());

                } else r = kFALSE;
            } else r = kFALSE;
        } else { strcpy(buf,""); }
    }
    return r;
}
ClassImp(HMdcCutStat)
ClassImp(HMdcCutStatMod)
ClassImp(HMdcCutStatSec)