TTools.h

Go to the documentation of this file.

// Class TTools
//
// Author     : F. Cassol Brunner
// Created    : 12/9/2002
// Last update: $Date: 2005/06/28 21:19:21 $
//          by: $Author: sauvan $
// Comment: class for general tool routines

#ifndef __TTOOLS
#define __TTOOLS

#include <cstdio>
#include <cstdlib> 
#include <iostream> 
#include <fstream> 

#include "H1Skeleton/H1Tree.h"
#include "H1Pointers/H1BytePtr.h"
#include "H1Pointers/H1IntPtr.h"
#include "H1Pointers/H1ShortPtr.h"
#include "H1Pointers/H1FloatPtr.h"
#include "H1Mods/H1PartSelTrack.h"
#include "H1Mods/H1PartMuon.h"
#include <TKey.h>
#include <TVector3.h>

#define MYR2D 180.0/3.14159265359

// Definition of DATA/MC convention
#define DATAID 0
#define MCDJANGONCID 1
#define MCDJANGOCCID 2
#define MCGPRODDIRID 3
#define MCGPRODRESID 4
#define MCGPRODCCID 5
#define MCGPRODBBID 6
#define MCRAPGAPNCID 7
#define MCGRAPEID 8
#define MCEPVECID 9
#define MCCOMPTONID 10


// Define B Field constants for track extrapolation
#define BFIELDCONST 0.0003
#define BFIELDNOM   11.6083

// Define calorimeter radius
#define rCalo 100.5

// Define z of BBE, backward detector and IF1
#define zBACKDET -150.0
#define zEDBBE -218.5
#define zIF1 292.0


using namespace std;

class TTools {
public:
  Int_t num;                       // event counter
  Int_t numMax;                    // maximum number of events

  Int_t RunNumber;                 // Runnumber
  Int_t EventNumber;               // Eventnumber
  Bool_t NewRun;                   // Set to kTrue if new run started
  Float_t GenPl;
  Float_t GenPp;
  Short_t MCFlag;                  // Flag pour MC 

  TTools();
  TTools(Int_t Run, Int_t Event,Int_t MCFlag);

  virtual Bool_t Next();

//--- set an event list to the tree
  void SetEventList(Char_t *ListFileName);

//--- Which data set is being analysed
  Bool_t IsData2005(void);
  Bool_t IsData2004e(void);
  Bool_t IsData2004p(void);
  Bool_t IsData0304R1(void);
  Bool_t IsData0304L1(void);
  Bool_t IsData0304L2(void);
  Bool_t IsData0304R2(void);
  Bool_t IsData2003(void);
  Bool_t IsData2002(void);
  Bool_t IsDataHeraII(void);
  Bool_t IsData2000(void);
  Bool_t IsData1999p(void);
  Bool_t IsData1999e(void);
  Bool_t IsData1998(void);
  Bool_t IsData1997(void);
  Bool_t IsData1996(void);
  Bool_t IsData1995(void);
  Bool_t IsData1994(void);

  Bool_t IsData94(const UInt_t& irun);
  Bool_t IsData95S(const UInt_t& irun);
  Bool_t IsData95P(const UInt_t& irun);
  Bool_t IsData96P(const UInt_t& irun);
  Bool_t IsData97P(const UInt_t& irun);
  Bool_t IsData97S(const UInt_t& irun);
  Bool_t IsData98E(const UInt_t& irun);
  Bool_t IsData99S(const UInt_t& irun);
  Bool_t IsData99E(const UInt_t& irun);
  Bool_t IsData99P(const UInt_t& irun);
  Bool_t IsData00S(const UInt_t& irun);
  Bool_t IsData00P(const UInt_t& irun);


//--- Test trigger bits
  Bool_t IsL4Trigg(Int_t iTrigg);
  Bool_t IsL1Trigg(Int_t iTrigg);
  Bool_t IsL2TTTrigg(Int_t iTrigg);

  // Are HV bits on
  Bool_t IsHV(void);
  Bool_t IsHVCJCLoose(void);
  Bool_t IsCjcLooseOn();  

//--- Are CIP/COP on ?
  Bool_t IsCipCopOn();

//--- Is CIP on ?
  Bool_t IsCipOn();

  Bool_t IsBdcOn();
  Bool_t IsFmdOn();
  Bool_t IsPlugOn();

  // Is cosmic
  Bool_t IsCosmic(void) ;
  Bool_t IsIbgBit(Int_t iBit);

  // Is Special CC cosmic
  Bool_t IsCCCosmic(void);


  // Is NOISE or cosmic
  Bool_t IsNonEP();
  Bool_t IsNonEP03(void);

  // Is E-Pz physical (E-Pz < 75 GeV)
  Bool_t IsEPzOK(Float_t Epz);

  // Is Phi Crack
  Bool_t IsPhiCrack(Float_t Phi);
  Bool_t IsLArCrack(Float_t z,Float_t phi);
  Bool_t IsLArCrack(Float_t z,Float_t zBBE,Float_t phi);

  // Is z Crack
  Bool_t IsZCrack(Float_t z);

//--- Is timing ok
  Bool_t IsTimingOk();
  Bool_t Is2004CalibratedTimingOK(void);

  // Is CJC dead (Dead periond in 1999 data)
  Bool_t IsCJCDead(Float_t PhiTr);
  Bool_t IsCJCMCDead(Float_t PhiTr);

  // Is Q2 OK (Combination of diff. MCs)
  Bool_t IsMCQ2OK(void);

  Bool_t IsHVCaloCJC(void);


  // Calibrate jets with calib of qhadcor_98_00
  Float_t abscor(Int_t icalfac,Float_t theta_rad,Float_t ptj,Float_t qvar);

  // Get phi angle of 2 particles
  Double_t GetMod(Double_t x, Double_t y);
  Double_t GetRealPhi(Double_t PhiIn);
  Double_t GetMergedPhi(Double_t PhiPart1,Double_t PtPart1,Double_t PhiPart2,Double_t PtPart2);
  
  // Phi projection into (-22.5,22.5) to see cracks
  Double_t GetProjPhi(Double_t Phi);

//---- test grade of muon
  Int_t GetMuonGrade(H1PartMuon *Muon);

 //... Read in noise runs
  Bool_t ReadNoiseRuns(Char_t *InputFile);

 

  // QHADCOR calibration: NOT FULLY IMPLEMENTED
  Float_t qhadcor_98_00(Int_t ITRA,Int_t ITYPE,Float_t jetpos,Float_t ptin,Float_t qvar);
  Float_t qhadcor_98_00_trans(Int_t *itra,Int_t *itype,Float_t *jetpos,Float_t *ptin, Float_t *qvar);


//---- for track extrapolations
  void CalcExtraPolTrack(TVector3 XStart,H1PartSelTrack *Track,
                Float_t &XOut,Float_t &YOut,Float_t &ZOut,Float_t &ThTrOut,Float_t &PhTrOut);
  void OctgXY(Float_t *xl,Float_t phi,Int_t ic,Float_t rc,
                                Float_t roctog,Float_t *xout,Int_t &iflag);
  Int_t OctogonIntersection(TVector3 ClusStart,Float_t Theta, Float_t Phi,
        Float_t RCurv,Float_t Zmin,Float_t Zmax,Float_t RBarr,TVector3 &XOut,Float_t &PhiOut);
  void LinCir(Float_t *xl,Float_t phi,Float_t *xc,Float_t rc,
                        Float_t *xy1, Float_t *xy2, Int_t &iflag);
  

  ClassDef(TTools,1) // Make class known to root 
};

#endif

---