| 1 | 
/** | 
| 2 | 
 * \file CaloPreSampler.h | 
| 3 | 
 * \author Emiliano Mocchiutti (2007/07/18) | 
| 4 | 
 */ | 
| 5 | 
#ifndef calopresampler_h | 
| 6 | 
#define calopresampler_h | 
| 7 | 
 | 
| 8 | 
#include <PamLevel2.h> | 
| 9 | 
 | 
| 10 | 
#include <TTree.h> | 
| 11 | 
#include <TFriendElement.h> | 
| 12 | 
#include <TChain.h> | 
| 13 | 
#include <TFile.h> | 
| 14 | 
#include <TList.h> | 
| 15 | 
#include <TKey.h> | 
| 16 | 
#include <TSystemFile.h> | 
| 17 | 
#include <TSystemDirectory.h> | 
| 18 | 
#include <TSQLServer.h> | 
| 19 | 
 | 
| 20 | 
#include <iostream> | 
| 21 | 
 | 
| 22 | 
using namespace std; | 
| 23 | 
 | 
| 24 | 
/** | 
| 25 | 
 * | 
| 26 | 
 * Class to store and calculate variables useful for presampler analysis | 
| 27 | 
 */ | 
| 28 | 
class CaloPreSampler : public TObject { | 
| 29 | 
 | 
| 30 | 
 private: | 
| 31 | 
    // | 
| 32 | 
    PamLevel2 *L2; | 
| 33 | 
    Bool_t debug; | 
| 34 | 
    // | 
| 35 | 
    // needed to avoid reprocessing the same event over and over to obtain the variables; | 
| 36 | 
    // | 
| 37 | 
    UInt_t OBT; | 
| 38 | 
    UInt_t PKT; | 
| 39 | 
    UInt_t atime; | 
| 40 | 
    Bool_t ssel; | 
| 41 | 
    // | 
| 42 | 
    Bool_t usepl18x; | 
| 43 | 
    // | 
| 44 | 
    Bool_t simulation; ///< True when using simulated data, false by default; | 
| 45 | 
    // | 
| 46 | 
    CaloLevel0 *event; ///< CaloLevel0 event needed to process data from level0/1 to level2 | 
| 47 | 
    CaloStrip *cstrip; ///< CaloStrip object needed to easily recover parameter files | 
| 48 | 
    CaloLevel2 *pcalo; ///< CaloLevel2 object which contains variables calculated with a reduced calorimeter. | 
| 49 | 
    //   CaloLevel1 *c1; ///< CaloLevel1 object which contains variables calculated with a reduced calorimeter. | 
| 50 | 
    // | 
| 51 | 
    Int_t N; ///< Number of W planes to be used as presampler. NOTICE: none the silicon detectors attached to these planes will be used; | 
| 52 | 
    Int_t NC; ///< Number of W planes to be used as calorimeter. | 
| 53 | 
    Bool_t sel; ///< Selection mode: planes from 1 to 22-N are used, plane 18 - N is masked if "emulate18" variable is true; | 
| 54 | 
    Bool_t cont; ///< Contamination mode: planes from N to 22 are used; | 
| 55 | 
    Bool_t emulate18; ///< If true emulates the non-working plane 18; | 
| 56 | 
    Bool_t nox; ///< If true do not use X-view | 
| 57 | 
    Bool_t noy; ///< If true do not use Y-view | 
| 58 | 
    Bool_t forcecalo; | 
| 59 | 
    Int_t mask[2][22]; | 
| 60 | 
    // | 
| 61 | 
    Float_t rigdefault; ///< Rigidity to be used if no track is available or in standalone mode | 
| 62 | 
    Bool_t withtrk; ///< Use or not the tracker | 
| 63 | 
    Int_t forcefitmode; | 
| 64 | 
    // | 
| 65 | 
 | 
| 66 | 
 public: | 
| 67 | 
    // | 
| 68 | 
    // | 
| 69 | 
    // CaloLevel1* GetCaloLevel1(){ Process(); return c1;}; ///< Retrieve CaloLevel1 pointer which contains variables calculated with a reduced calorimeter. | 
| 70 | 
    CaloLevel2* GetCaloLevel2(){ Process(); return pcalo;}; ///< Retrieve CaloLevel2 pointer which contains variables calculated with a reduced calorimeter. | 
| 71 | 
    CaloLevel2* GetLevel2Pointer(){ printf(" OBSOLETE, use GetCaloLevel2Pointer()!\n"); return pcalo;}; ///< Retrieve CaloLevel2 pointer which contains variables calculated with a reduced calorimeter. | 
| 72 | 
    CaloLevel2* GetCaloLevel2Pointer(){ return pcalo;}; ///< Retrieve CaloLevel2 pointer which contains variables calculated with a reduced calorimeter. | 
| 73 | 
    // | 
| 74 | 
    CaloTrkVar* GetCaloTrack(Int_t t); | 
| 75 | 
    // | 
| 76 | 
    void SetNoWpreSampler(Int_t n); ///< Set the number of W planes to be used as presampler. NOTICE: none the silicon detectors attached to these planes will be used. Default: N = 4; | 
| 77 | 
    void SetNoWcalo(Int_t n); ///< Set the number of W planes to be used as calorimeter. Default: NC = 22-N (N presampler); | 
| 78 | 
    void SplitInto(Int_t NoWpreSampler, Int_t NoWcalo); | 
| 79 | 
    void Selection(){sel = true; cont = false;}; ///< Set selection mode: planes from 1 to 22-N are used, plane 18 - N is masked if "emulate18" variable is true (DEFAULT); | 
| 80 | 
    void Contamination(){sel = false; cont = true;}; ///< Set contamination mode: planes from N to 22 are used. | 
| 81 | 
    void EmulatePlane18(Bool_t emu){emulate18 = emu;} ///< Emulates the non-working plane 18; Default = true; | 
| 82 | 
    void Simulation(Bool_t sim){simulation = sim;}; ///< Set simulated data flag. Default FALSE; | 
| 83 | 
    // | 
| 84 | 
    CaloPreSampler(); | 
| 85 | 
    CaloPreSampler(PamLevel2 *L2); | 
| 86 | 
    ~CaloPreSampler(){ Delete(); }; | 
| 87 | 
    // | 
| 88 | 
    Int_t GetNoWpreSampler(){return N;}; ///< Get the number of W planes used as presampler.  | 
| 89 | 
    Int_t GetNoWcalo(){return NC;}; ///< Get the number of W planes used as calorimeter. | 
| 90 | 
    void SetDebug(Bool_t d){ debug=d; }; | 
| 91 | 
    void UseTracker(Bool_t ch){ withtrk = ch; }; | 
| 92 | 
    void SetDefaultRig(Float_t rig){ rigdefault = rig; }; | 
| 93 | 
    void SetForceFitMode(Int_t mode){ forcefitmode = mode; }; | 
| 94 | 
    void MaskX(Bool_t mskx){ nox = mskx; }; | 
| 95 | 
    void MaskY(Bool_t msky){ noy = msky; }; | 
| 96 | 
    void Mask(Int_t view, Int_t plane){ mask[view][plane] = 1;}; | 
| 97 | 
    void ForceCaloFit(){forcecalo=true;}; | 
| 98 | 
    void ForceNextProc(){atime=0;}; | 
| 99 | 
    // | 
| 100 | 
    void UsePlane18X(Bool_t use){usepl18x = use;}; | 
| 101 | 
    // | 
| 102 | 
    void Clear(); | 
| 103 | 
    void Clear(Option_t *option){Clear();}; | 
| 104 | 
    void Delete(); | 
| 105 | 
    void Delete(Option_t *option){Delete();}; | 
| 106 | 
    // | 
| 107 | 
    void Process(); ///< Process data  | 
| 108 | 
    void Print(); | 
| 109 | 
    void Print(Option_t *option){Print();}; | 
| 110 | 
    // | 
| 111 | 
    CaloTrkVar* AddCaloTrkVar(float *al,int trktag);//ELENA | 
| 112 | 
 | 
| 113 | 
    ClassDef(CaloPreSampler,3); | 
| 114 | 
}; | 
| 115 | 
 | 
| 116 | 
#endif | 
| 117 | 
 |