CgemClusterCreate Class Reference

#include <CgemClusterCreate.h>

Inheritance diagram for CgemClusterCreate:

Public Member Functions
	CgemClusterCreate (const std::string &name, ISvcLocator *pSvcLocator)
	~CgemClusterCreate ()
StatusCode	initialize ()
StatusCode	execute ()
StatusCode	finalize ()

Detailed Description

Definition at line 37 of file CgemClusterCreate.h.

Constructor & Destructor Documentation

CgemClusterCreate()

CgemClusterCreate::CgemClusterCreate	(	const std::string &	name,
		ISvcLocator *	pSvcLocator )

Definition at line 73 of file CgemClusterCreate.cxx.

                                                                                   :
    Algorithm(name,pSvcLocator)
{ 
    declareProperty("PrintFlag",myPrintFlag=0);
    declareProperty("MotherParticleID",myMotherParticleID=443);
    declareProperty("Method",myMethod=2);
    declareProperty("ntuple",myNtuple=0);
    declareProperty("effCluster",myClusterEff=1.0);
    declareProperty("fillOpt",m_fillOption=-1);
    declareProperty("selGoodDigi",m_selGoodDigi=1);
    declareProperty("minDigiTime",m_minDigiTime=-8875);
    declareProperty("maxDigiTime",m_maxDigiTime=-8562);
    declareProperty("TPCFitMethod",m_selectTPC=1);
    declareProperty("minQDigi",myQMin=0.0);
    declareProperty("LUTfile", m_LUTfile = "LUTfile.root");
    declareProperty("clustering_gap", clustering_gap = 0);
    declareProperty("dropHitRandomly", myDropHit = 0);
    declareProperty("uTPC_correction", uTPC_correction = true);
 
    // --- for random noise (toy cluster)   
    declareProperty("addRandomToyCluster", myRandomCluster=0);
    declareProperty("NoiseRate_layer1",    myNoiseRate[0]=0.0);
    declareProperty("NoiseRate_layer2",    myNoiseRate[1]=0.0);
    declareProperty("NoiseRate_layer3",    myNoiseRate[2]=0.0);
 
 
    if(myMethod==0) reset(); 
    else if(myMethod==1||myMethod==3) resetFiredStripMap();
 
    m_totEvt = 0;
    m_nStripRead=0;
    m_iStripToDrop=-1;
}

~CgemClusterCreate()

CgemClusterCreate::~CgemClusterCreate

(

)

Definition at line 107 of file CgemClusterCreate.cxx.

{
    delete lutreader;
}

Member Function Documentation

execute()

StatusCode CgemClusterCreate::execute

(

)

Definition at line 201 of file CgemClusterCreate.cxx.

{
    MsgStream log(msgSvc(), name());
    log << MSG::INFO << "in execute()" << endreq;
 
    // check and register /Event/Recon
    DataObject *aReconEvent;
    eventSvc()->findObject("/Event/Recon",aReconEvent);
    if(aReconEvent==NULL) {
        aReconEvent = new ReconEvent();
        StatusCode sc = eventSvc()->registerObject("/Event/Recon",aReconEvent);
        if(sc!=StatusCode::SUCCESS) {
            log << MSG::FATAL << "Could not register ReconEvent" <<endreq;
            return( StatusCode::FAILURE);
        }
    }
 
 
    if(myMethod==0)      method0();
    else if(myMethod==1) method1();
    else if(myMethod==2) toyCluster();
    else if(myMethod==3) method2();
    m_totEvt++;
 
    return StatusCode::SUCCESS;
}

finalize()

StatusCode CgemClusterCreate::finalize

(

)

Definition at line 2132 of file CgemClusterCreate.cxx.

{
    MsgStream log(msgSvc(),name());
    log << MSG::INFO << "in finalize(): Number of events in CgemClusterCreate" << m_totEvt << endreq;
    //delete myNTupleHelper;
 
    return StatusCode::SUCCESS;
}

initialize()

StatusCode CgemClusterCreate::initialize

(

)

Definition at line 112 of file CgemClusterCreate.cxx.

{
    MsgStream log(msgSvc(), name());
    log << MSG::INFO << "in initialize()" << endreq;
 
    //CgemGeomSvc//
    //      StatusCode sc;
    //      ICgemGeomSvc* ISvcCgem;
    //      sc = service("CgemGeomSvc", ISvcCgem);
    //      if(sc != StatusCode::SUCCESS)
    //      {
    //    log << MSG::ERROR << "can not use CgemGeomSvc" << endreq;
    //    return StatusCode::FAILURE;
    //      }
 
    //BesTimerSvc//
    StatusCode tsc;
    IBesTimerSvc* m_timersvc;
    tsc = service( "BesTimerSvc", m_timersvc);
    if( tsc.isFailure() ) 
    {
        log << MSG::WARNING << name() << " Unable to locate BesTimerSvc" << endreq;
        return StatusCode::FAILURE;
    }
    m_timer = m_timersvc->addItem("Execution");
 
    if(myNtuple) myNTupleHelper=new NTupleHelper(ntupleSvc()->book("RecCgem/CgemCluster",CLID_ColumnWiseTuple,"CgemCluster"));
 
    if(myMethod==0||myMethod==2) {
        if(myNtuple) hist_def();
    }
 
    // --- get CgemGeomSvc ---
    ISvcLocator* svcLocator = Gaudi::svcLocator();
    ICgemGeomSvc* ISvc;
    StatusCode sc=svcLocator->service("CgemGeomSvc", ISvc);
    myCgemGeomSvc=dynamic_cast<CgemGeomSvc *>(ISvc);
    //StatusCode sc = service( "CgemGeomSvc", myCgemGeomSvc);
    if (!sc.isSuccess()) log<< MSG::INFO << "CgemClusterCreate::initialize(): Could not open CGEM geometry file" << endreq;
    myNCgemLayers = myCgemGeomSvc->getNumberOfCgemLayer();
    if(myPrintFlag) cout<<"CgemClusterCreate::initialize() "<<myNCgemLayers<<" Cgem layers"<<endl;
    for(int i=0; i<myNCgemLayers; i++)
    {
        myCgemLayer[i]=myCgemGeomSvc->getCgemLayer(i);
        myNSheets[i]=myCgemLayer[i]->getNumberOfSheet();
        bool IsReverse = myCgemLayer[i]->getOrientation();
        myRXstrips[i] = myCgemLayer[i]->getInnerROfAnodeCu2();
        myRVstrips[i] = myCgemLayer[i]->getInnerROfAnodeCu1();
        if(IsReverse) {
            myRXstrips[i] = myCgemLayer[i]->getOuterROfAnodeCu2();
            myRVstrips[i] = myCgemLayer[i]->getOuterROfAnodeCu1();
        }
        if(myRandomCluster>0) {
            myNoiseRate[i]*=myCgemLayer[i]->getNumberOfChannelX();
            cout<<"going to add toy cluster as noise with noise rate "<<myNoiseRate[i]<<" for CGEM layer "<<i<<endl;
        }
        if(myPrintFlag) cout<<"layer "<<i<<": "<<myNSheets[i]<<" sheets"<<", is reversed "<<IsReverse<<", RX="<<myRXstrips[i]<<", RY="<<myRVstrips[i]<<endl;
    }
 
    sc = service("CgemGeomSvc", m_SvcCgem);
    if(sc != StatusCode::SUCCESS)  {
      log << MSG::ERROR << "can not use CgemGeomSvc" << endreq;
      return StatusCode::FAILURE;
    }
 
    anode = m_SvcCgem->getReadoutPlane(0, 0);
    anode_mid_gap_radius[0] = anode->getMidRAtGap();
    anode = m_SvcCgem->getReadoutPlane(1, 0);
    anode_mid_gap_radius[1] = anode->getMidRAtGap();
    anode = m_SvcCgem->getReadoutPlane(2, 0);
    anode_mid_gap_radius[2] = anode->getMidRAtGap();
 
 
    // --- get CgemCalibFunSvc ---
    sc = service ("CgemCalibFunSvc", myCgemCalibSvc);
    if ( sc.isFailure() ){
        cout<< name() << "Could not load MdcCalibFunSvc!" << endl;
        return sc;
    }
    //cout<<"sigma from CgemCalibFunSvc: "<<myCgemCalibSvc->getSigma(0,0,0,0,0,0)<<endl;
 
 
    // LUT
    lutreader = new CgemLUTReader(m_LUTfile);
    lutreader->ReadLUT();
 
    return StatusCode::SUCCESS;
}

---

The documentation for this class was generated from the following files:

• CgemClusterCreate.h
• CgemClusterCreate.cxx