CalDetPIDModule.cxx

Go to the documentation of this file.

//$Id: CalDetPIDModule.cxx,v 1.3 2003/05/21 16:10:08 bv Exp $
//
//CalDetPIDModule.cxx
//
//Module for identifying particles in the caldet
//
//Author:  B. Rebel 10/2000
#include "TString.h"

#include "CalDetPIDSR/CalDetPIDModule.h"
#include "CalDetPIDSR/CalDetTOFId.h"
#include "CalDetPIDSR/CalDetCerenkovId.h"
#include "JobControl/JobCModuleRegistry.h"
#include "JobControl/JobCommand.h"
#include "MessageService/MsgService.h"
#include "MinosObjectMap/MomNavigator.h"
#include "CandData/CandRecord.h"

#include <cassert>

ClassImp(CalDetPIDModule)
CVSID("$Id: CalDetPIDModule.cxx,v 1.3 2003/05/21 16:10:08 bv Exp $");

// Declare this module to JobControll. Arguments are:
//  (1) The class name 
//  (2) The human-readable name 
//  (3) A short, human-readable description of what the module does
JOBMODULE(CalDetPIDModule, 
          "CalDetPIDModule",
          "A module used for particle id in the CalDet");

//......................................................................
CalDetPIDModule::CalDetPIDModule() :
  fBeamMomentum(1.0),
  fCkCollectionEfficiency(0.5),
  fCkDetectorEfficiency(0.27), //from pdg 2001
  fCkLength(6.0),
  fCkPedestal(0),
  fMassSq(-1.),
  fTDCPedestal1(0),
  fTDCPedestal2(0),
  fTDCGain1(0.1),
  fTDCGain2(0.1),
  fTOFSeparation(11.0)
{

  MSG("CalDPID", Msg::kDebug) << "CalDetPIDModule::Constructor" << endl;

}

//----------------------------------------------------------------------
CalDetPIDModule::~CalDetPIDModule()
{

  MSG("CalDPID", Msg::kInfo) << "CalDetPIDModule::Destructor" << endl;

}

//......................................................................

void CalDetPIDModule::BeginJob()
{
  //load in the appropriate root file with histogram here.    
  return;
}

//......................................................................

JobCResult CalDetPIDModule::Reco(MomNavigator *mom)
{
  MSG("CalDPID", Msg::kDebug) << "CalDetPIDModule::Reco\n";

  // Check that mom exists.
  assert(mom);

  //find the Raw record fragment in MOM - to be added.
  
  //get the raw tof digit
  RawTOFDigit rtofd;

  //instantiate the CandTOFId and CandCerenkovId objects
  CalDetTOFId tof = CalDetTOFId(fTOFSeparation, 
                                fTDCGain1, 
                                fTDCGain2, 
                                fTDCPedestal1, 
                                fTDCPedestal2, 
                                fFirstTOFPaddle, 
                                fSecondTOFPaddle, 
                                rtofd);
  //CalDetCerenkovId ck = CalDetCerenkovId(rcd);

  //Get the Mass Squared value
  fMassSq = tof.GetParticleMassSquared(fBeamMomentum);

  return JobCResult::kAOK; // All OK...
}

//......................................................................

void CalDetPIDModule::HandleCommand(JobCommand *command) 
{
  //
  //  Purpose:  Method to interpret module commands.
  //
  //  Arguments:
  //    command   in    Command to interpret.
  //
  //  Return:   n/a
  //
  //  Commands implemented:
  //    Set TOFSeparation  Set the separation to use for tof in meters.
  //    Set CerenkovPressure    Set the pressure of the cerenkov detector in bar.
  //    Set CerenkovPedestal     Set the pedestal value for cerenkov adc.
  //    Set CerenkovCollectionEfficiency  Set the collection efficiency of the cerenkov detector.
  //    Set CerenkovDetectorEfficiency  Set the integrated efficiency of the cerenkov detector.
  //    Set CerenkovLength      Set the length of the cerenkov detector in meters.
  //    Set BeamMomentum   Set the momentum of the beam, enter in GeV/c
  //    Set TDCPedestal Set the pedestal value for the TDC
  //    Set TDCGain  Set the gain for the TDC
  //    Set FirstTOFPaddle  Set which paddle is the first one used for tof info
  //    Set SecondTOFPaddel Set which paddle is the second one used for tof info
  //

   MSG("CalDPID", Msg::kDebug) << "CalDetPIDModule::HandleCommand" << endl;

   TString cmd = command->PopCmd();
   if(cmd == "Set"){   
     TString opt = command->PopOpt();
     if(opt == "TOFSeparation"){ fTOFSeparation = command->PopFloatOpt();}
     else if(opt == "CerenkovPressure"){ fCkPressure = command->PopFloatOpt();}
     else if (opt == "CerenkovPedestal"){ fCkPedestal = command->PopIntOpt();}
     else if (opt == "CerenkovCollectionEfficiency"){fCkCollectionEfficiency = command->PopFloatOpt();}
     else if (opt == "CerenkovDetectorEfficiency"){ fCkDetectorEfficiency = command->PopFloatOpt();}
     else if (opt == "CerenkovLength"){ fCkLength = command->PopFloatOpt();}
     else if (opt == "BeamMomentum"){ fBeamMomentum = command->PopFloatOpt();}
     else if (opt == "TDCPedestal1"){ fTDCPedestal1 = command->PopIntOpt();}
     else if (opt == "TDCGain1"){ fTDCGain1 = command->PopFloatOpt();}
     else if (opt == "TDCPedestal2"){ fTDCPedestal2 = command->PopIntOpt();}
     else if (opt == "TDCGain2"){ fTDCGain2 = command->PopFloatOpt();}
     else if (opt == "FirstTOFPaddle"){ fFirstTOFPaddle = command->PopIntOpt();}
     else if (opt == "SecondTOFPaddle"){ fSecondTOFPaddle = command->PopIntOpt();}
     else {
         MSG("CalDPID", Msg::kWarning)<< "CalDetPIDModule: Unrecognized option " << opt << endl;
      }
   }
   else {
     MSG("CalDPID", Msg::kWarning) << 
       "CalDetPIDModule: Unrecognized command " << cmd << endl;
   }
}

//......................................................................

void CalDetPIDModule::Help() 
{
  MSG("CalDPID", Msg::kInfo) 
    << "CalDetPIDModule::Help\n"
    <<"CalDetPIDModule is a module which identifies particles "
    <<"in the CalDet."
    << endl
    << "Commands implemented:"
    << "Set TOFSeparation  Set the separation to use for tof in meters."
    << "Set CerenkovPressure    Set the pressure of the cerenkov detector in bar."
    << "Set CerenkovPedestal     Set the pedestal value for cerenkov adc."
    << "Set CerenkovCollectionEfficiency  Set the collection efficiency of the cerenkov detector."
    << "Set CerenkovDetectorEfficiency  Set the integrated efficiency of the cerenkov detector."
    << "Set CerenkovLength      Set the length of the cerenkov detector in meters."
    << "Set BeamMomentum   Set the momentum of the beam, enter in GeV/c"
    << "Set TDCPedestal Set the pedestal value for the TDC"
    << "Set TDCGain  Set the gain for the TDC, in units of ns/tdc"
    << "Set FirstTOFPaddle  Set which paddle is the first one used for tof info"
    << "Set SecondTOFPaddel Set which paddle is the second one used for tof info"
    << endl;

  
}

---