/data3/calcul/jacquem/working_dir/Micromegas/micromegasFrameWork/include/root/BeamProfile.hh

Go to the documentation of this file.

/* @version $Revision: 1566 $ * @modifiedby $Author: jacquem $ * @lastmodified $Date: 2012-03-08 14:00:10 +0100 (Thu, 08 Mar 2012) $ */
/* 
 * File:   BeamProfile.hh
 * Author: jacquem
 *
 * Created on May 6, 2010, 10:18 AM
 */
#include <string.h>
#include <map>
#include <string>
#include "TH1I.h"
#include "TH2I.h"
#include  "mTypes.h"

#include "root/MTChannel.hh"
//#include "geometry/Detector.hh"

class Detector;


#ifndef _COMMANDOPTION_HH
#define _COMMANDOPTION_HH


#define TIMECUT 1

class MTEvent;

class BeamProfile {
public:
    BeamProfile();
    BeamProfile(int argc, char **argv);
    BeamProfile(const BeamProfile& orig);

                void SetDetector( Detector* det) ;


    inline std::string  getFileName()              { return fileName; };
    inline std::string  getXMLFileName()           { return xmlFileName; };
    inline unsigned int          getTimeCutMin()   { return timeCutMin;};
    inline unsigned int          getTimeCutMax()   { return timeCutMax;};
    inline bool         getShowRunInfo()           { return showRunInfo;};
    inline bool         getShowXYDistribution()    { return showXYDistribution;};
    inline bool         getShowDeltaDistribution() { return showDeltaDistribution;};
    inline bool    getShowHitPerChipDistribution() { return showHitPerChipDistribution;};
    inline bool getShowHitPerChannelDistribution() { return showHitPerChannelDistribution;};
    inline int          getChamberId()             { return chamberId;};
    inline int          getDifId()                 { return difId;};
    inline int          getBoardId()               { return boardId;};
    inline int          getChipId()                { return chipId;};
        //std::map<ui16,std::string>&   getHotChip(void) const;

    bool FillXyDist( const  MTChannel &channel);
    bool FillHitPerChipDist(const   MTChannel &channel);
    bool FillHitPerChannel(const   MTChannel &channel);
    bool FillDeltaT(const   MTChannel &channel);
    bool FillHotChamber(const MTEvent& evt);
        bool FillHotChip( const MTChannel &channel);

    unsigned int GetHotChamber(const MTEvent& evt);

    void DrawXyDist(void);
    void DrawHitPerChipDist();
    void DrawHitPerChannelDist();
    void DrawChannel2D();
                                void DrawChip2D();
    void DrawDeltaT();
    void DrawHotChamber();

    void PrintHotChip() const;

    virtual ~BeamProfile();
private:
    void parseOptions();
    void askOptions();
    int getNbHist();
    bool filterOnId(const MTChannel &channel);
    bool filterOnCut(const MTChannel &channel);


private:
    int nbArgs;
    char **args;

    std::string fileName;
    std::string xmlFileName;
    unsigned int  timeCutMin;
    unsigned int  timeCutMax;
    int  analogCutMin;
    int  analogCutMax;
    int maxHistogram;
    //unsigned int cutType;
    bool analogCut;
    bool timeCut;
    bool showRunInfo;
    bool showXYDistribution;
    bool showDeltaDistribution;
    int nbhit;
    int adc;
    bool showHitPerChipDistribution;
    bool showHitPerChannelDistribution;
    bool showHotChamber;
    bool showHotChip;
    int chamberId;
    int difId;
    int boardId;
    int chipId;

    const Detector* detector;
    // Create map to store xy distribution histogram
    std::map<ui16,TH2I*> xyKeyChamber; // key is chamberId
    std::map<ui32,TH1I*> hitPerChipKeyChamberDif; // key is ChamberId and DifId to string;
    std::map<ui32,TH1I*> hitPerChannelKeyChip; // key is chip softId
    ui32 maxHitPerChannel;
    ui32 minHitPerChannel;
    std::map<ui32,TH1I*> deltaTKeyChip; // key is chip softId

    TH1I *chamberHotHist; //

    int minDeltaTValue;
    int maxDeltaTValue;
    int minHitPerChannelValue;
    int maxHitPerChannelValue;
                                int maxHitPerChip;
                                int minHitPerChip;
    std::map <std::string,int>  minHitPerChipValue;  // key i chamberd
    std::map <std::string,int>  maxHitPerChipValue;  // key i chamberd
    std::map<ui16,std::map<std::string,int>* >  nbHitPerChipPerChamber;


    std::map<std::string,std::map<Int_t ,Int_t>* > hitPerChipKeyChamberDifMap;
    std::map<std::string,std::map<Int_t ,Int_t>* > deltaTKeyChipMap;
    std::map<std::string,std::map<Int_t ,Int_t>* > hitPerChannelKeyChipMap;
};

#endif  /* _COMMANDOPTION_HH */

---