EventReducer.h

Go to the documentation of this file.

/*
EventReducer.h

This class is intended to handle all the hit list manipulation/reduction,
The reduction methods could initialized once for all, and then it will update
the hit list stored in SRawEvent. It is declared as a friend class of SRawEvent.

Author: Kun Liu, liuk@fnal.gov
Created: 06-17-2015
*/

#ifndef _EVENTREDUCER_H
#define _EVENTREDUCER_H

#include <GlobalConsts.h>

#include <list>
#include <map>
#include <TString.h>
#include <TRandom.h>

#include "SRawEvent.h"
#include "TriggerAnalyzer.h"

#include <phool/recoConsts.h>
#include <geom_svc/GeomSvc.h>

class EventReducer
{
public:
    EventReducer(TString options);
    ~EventReducer();

    //main external call
    int reduceEvent(SRawEvent* rawEvent);

    //external handle to set chamber efficiency/resolution
    void setChamEff(double val)   { chamEff = val;   }
    void setChamResol(double val) { chamResol = val; }

    //sagitta ratio reducer
    void sagittaReducer();

    //hough transform reducer
    void houghReducer();

    //hit cluster remover
    void deClusterize();
    void processCluster(std::vector<std::list<Hit>::iterator>& cluster);

    //hodosope maksing
    void initHodoMaskLUT();
    void hodoscopeMask(std::list<Hit>& chamberhits, std::list<Hit>& hodohits);
    bool lineCrossing(double x1, double y1, double x2, double y2,
                      double x3, double y3, double x4, double y4);

private:
    //pointer to geometry service, inited outside
    GeomSvc* p_geomSvc;

    //pointer to the reco configuration 
    recoConsts* rc;

    //pointer to trigger analyzer, inited inside
    TriggerAnalyzer* p_triggerAna;

    //Random number
    TRandom rndm;

    //temporary container for the hit list
    std::list<Hit> hitlist;
    std::list<Hit> hodohitlist;

    //loop-up table of hodoscope masking
    typedef std::map<int, std::vector<int> > LUT;
    LUT h2celementID_lo;
    LUT h2celementID_hi;
    LUT c2helementIDs;

    //flags of the hit manipulation method
    bool afterhit;            //after pulse removal
    bool hodomask;            //hodoscope masking
    bool outoftime;           //out of time hit removal
    bool decluster;           //remove hit clusters in chamber
    bool mergehodo;           //merge trigger hit with hit
    bool triggermask;         //use active trigger road for track masking
    bool sagitta;             //remove the hits which cannot form a sagitta triplet
    bool hough;               //remove the hits which cannot form a peak in hough space, will be implemented later
    bool externalpar;         //re-apply the alignment and calibration parameters
    bool realization;         //apply detector efficiency and resolution by dropping and smear
    bool difnim;              //treat the nim/FPGA triggered events differently, i.e. no trigger masking in NIM events

    //Adjustable parameters
    double timeOffset;        //timing correction
    double chamEff;           //chamber efficiency
    double chamResol;         //chamber resolution

    //SagittaReducer parameters
    double SAGITTA_DUMP_CENTER;
    double SAGITTA_DUMP_WIDTH;
    double SAGITTA_TARGET_CENTER;
    double SAGITTA_TARGET_WIDTH;
    double Z_TARGET;
    double Z_DUMP;

    //Hodo masking parameters
    double TX_MAX;
    double TY_MAX;
    bool USE_V1495_HIT;
    bool USE_TWTDC_HIT;
};

#endif