FCPCFilterModule.h

Go to the documentation of this file.

//FCPCFilterModule.cxx
// A first pass Partially and fully contained Filter
// John Chapman chapman@hep.phy.cam.ac.uk
// (Original Code by Caius Howcroft)
#ifndef FCPCFILTERMODULE_H
#define FCPCFILTERMODULE_H
#ifndef JOBCMODULE_H
#include "JobControl/JobCModule.h"
#endif

#include <string>
#include "TString.h"

class TTree;
class TFile;

class FCPCFilterModule : public JobCModule
{
public:
  FCPCFilterModule();
  ~FCPCFilterModule();

public:
  // Handle job status changes
  void BeginJob();
  void EndJob();

  // Analysis and Reconstruction methods
  JobCResult Ana(const MomNavigator* mom);

  // Module configuration
  const Registry& DefaultConfig() const;
  void            Config(const Registry& r);


private:
  void ResetParams();
  Float_t TimeWalk(float adc);
  //Configuration data
  bool cDoPC;
  bool cDoFC;
  bool cDoUC;
  bool cDoUp;
  std::string cListIn;
  double cMinDmxWeight;
  double cMinDigitCharge;
  double cFidRadius;
  double cMinFidCharge;

  //Detector config
  int cSM1StartPlane;
  int cSM1StopPlane;
  int cSM2StartPlane;
  int cSM2StopPlane;


  //Cut config
  int cMinPlaneCut;        //Min required plane cut
  double cMinPlaneCharge;  
  int cEndPlaneVeto;       //The number of planes to veto on
  double cEndPlaneVetoCharge; //the cut on charge
  double cMaxFidPctChargeCut;  
  double cMaxRegionSep; //the maximum allowed distance between the mean hit postions in 2 adjacent regions
  
  //output variables
  int cMakeTree; //whether or not to make a tree of the results of FCPCFilter
  TFile* fFile;
  TTree* fTree;
  TString fFileName;
  //Event Info
  int fRun;
  int fSnarl;
  int fMinPlane;        //the max plane used in this module
  int fMaxPlane;        // the min plane used in this module
  double fCharge;       //total charge used (see above cuts)
  int fNPlane;          //Number of planes above some charge treshold
  double fFidCharge[8]; //the charge next to each side
  double fFidU[12];     //charge weighted mean U position in each side/end region
  double fFidV[12];     //charge weighted mean V position in each side/end region
  double fFidZ[12];     //charge weighted mean Z position in each side/end region
  double fSumFidCharge; //the total charge on all sides
  double fChargeInSM1;  //the charge within the SM1 fiducial volume
  double fChargeInSM2;  //the charge within the SM2 fiducial volum
  double fSM1begPlaneQ; //the low plane charge (SM1 beg)
  double fSM1endPlaneQ; //the SM1 end plane charge
  double fSM2begPlaneQ; //the SM2 beg charge
  double fSM2endPlaneQ; //the high plane charge (SM2 end)
  int fEventClass;      //the classification of the event 0=UC, 1=PC, 2=FC

  //some work areas
  double wPlaneChargePE[500];
  double wPlaneMeanTPos[500];
  double wPlaneZPos[500];
  int wPlaneContainment[2][2][2][2];
  double wMeanTime[12]; //upward going mus

};
#endif // FCPCFILTERMODULE_H

---