PndGemMagneticFieldVsTrackParameters.h

Go to the documentation of this file.

//* $Id: */

// -------------------------------------------------------------------------
// -----                    PndGemMagneticFieldVsTrackParameters header file                 -----
// -----                  Created 12/02/2009 by R. Karabowicz          -----
// -------------------------------------------------------------------------

#ifndef PNDGEMMAGNETICFIELDVSTRACKPARAMETERS_H
#define PNDGEMMAGNETICFIELDVSTRACKPARAMETERS_H 1

#include "FairTask.h"

#include "TStopwatch.h"

#include "TH1F.h"
#include "TH2F.h"
#include "TH3F.h"
#include "TF1.h"
#include "TF2.h"

#include <list>
#include <map>

class TClonesArray;
class PndGemDigiPar;

class PndGemMagneticFieldVsTrackParameters : public FairTask {

 public:
  PndGemMagneticFieldVsTrackParameters();

  PndGemMagneticFieldVsTrackParameters(Int_t iVerbose);

  PndGemMagneticFieldVsTrackParameters(const char *name, Int_t iVerbose);

  virtual ~PndGemMagneticFieldVsTrackParameters();

  virtual void Exec(Option_t *opt);

 private:
  //  FairGeoStsPar*     fGeoPar;       /** Geometry parameter container **/
  //  PndGemDigiPar*    fDigiPar;      /** Digitisation parameter container **/
  //  PndGemDigiScheme* fDigiScheme;   /** Digitisation scheme **/
  PndGemDigiPar *fDigiPar;

  Double_t fTrackFinderOnHits_ParThetaA;
  Double_t fTrackFinderOnHits_ParThetaB;

  Double_t fTrackFinderOnHits_ParTheta0;
  Double_t fTrackFinderOnHits_ParTheta1;
  Double_t fTrackFinderOnHits_ParTheta2;
  Double_t fTrackFinderOnHits_ParTheta3;

  Double_t fTrackFinderOnHits_ParRadPhi0;
  Double_t fTrackFinderOnHits_ParRadPhi2;

  Double_t fTrackFinderOnHits_ParMat0[3];
  Double_t fTrackFinderOnHits_ParMat1[3];

  TClonesArray *fMCTrackArray;
  TClonesArray *fMCPointArray;

  TStopwatch fTimer;

  Int_t fNofEvents; // event counter

  TList *fHistoList;
  // theta angle vs radius vs momentum
  TH3F *fhThetaVsRadiusVsMomentumAll;
  TH3F *fhThetaVsRadiusVsMomentum[10];
  TF2 *ffThetaVsRadiusVsMomentum[10];
  // radius vs angle
  TH2F *fhRadiusVsAngle;
  TH2F *fhRadiusVsAnglePair[20];
  TF1 *ffRadiusVsAngle[3];
  // momentum vs phi angle difference
  TH3F *fhMomentumVsPhiDiffVsRadius[10];
  TH2F *fhMomentumVsPhiDiffAll[10];
  TH2F *fhMomentumVsPhiDiff[10][20];
  TF2 *ffMomentumVsPhiDiffVsRadius[10];
  TF1 *ffMomentumVsPhiDiffAll[10];
  TF1 *ffMomentumVsPhiDiff[10][20];
  // transverse momentum vs trajectory helix radius
  TH2F *fhMomTransVsHRadius[10];
  // track phi versus hits phis
  TH3F *fhTrackPhiVsHitPhis[10];
  TH2F *fhTrackPhiVsCalcPhi[10];

  Int_t Fill1StationHistograms();
  Int_t Fill2StationsHistograms();

  void CreateHistos();

  virtual void SetParContainers();

  virtual void Finish();

  virtual InitStatus Init();

  virtual InitStatus ReInit();

  void Reset();

  ClassDef(PndGemMagneticFieldVsTrackParameters, 1);
};

#endif