/data3/calcul/jacquem/working_dir/Micromegas/micromegasFrameWork/src/lcio/MLcio.cpp

Go to the documentation of this file.

/* @version $Revision: 1328 $ * @modifiedby $Author: jacquem $ * @lastmodified $Date: 2011-10-03 17:04:17 +0200 (Mon, 03 Oct 2011) $ */

#include "geometry/Detector.hh"
#include "geometry/Chamber.hh"
#include "geometry/Dif.hh"
#include "geometry/Chip.hh"
#include "geometry/Board.hh"

#include "event/Run.hh"

#include "tools/MicroException.hh"
#include "tools/Log.hh"
/*
#include "ChannelHit.hh"
#include "Channel.hh"

#include "SlowControlEntry.hh"
#include "ChamberSlowControl.hh"


#include "UTIL/LCTime.h"
#include "IMPL/LCFlagImpl.h"
#include "IMPL/LCCollectionVec.h"
#include "IMPL/RawCalorimeterHitImpl.h"
#include "IMPL/CalorimeterHitImpl.h"
#include "IMPL/LCTOOLS.h"

#include <string>

#include "IO/LCWriter.h"
#include "EVENT/LCIO.h"
#include "DATA/LCFloatVec.h"
#include "DATA/LCIntVec.h"


#include "IMPL/TrackerRawDataImpl.h"
#include "IMPL/TrackerDataImpl.h"
#include "IMPL/TrackerPulseImpl.h"

#include "UTIL/LCRelationNavigator.h"
#include "UTIL/LCTime.h"
#include "UTIL/CellIDEncoder.h"
#include "UTIL/LCTypedVector.h"

#include "UTIL/LCSplitWriter.h"

static const int NEVENT = 10 ; // events
static const int NMCPART = 10 ;  // mc particles per event
static const int NHITS = 5 ;  // calorimeter hits per event

#include "IMPL/LCTOOLS.h"
*/

#include "lcio.h"
#include "IMPL/LCRunHeaderImpl.h"
#include "IMPL/LCEventImpl.h"
#include "IMPL/LCCollectionVec.h"
#include "IMPL/RawCalorimeterHitImpl.h"
#include "IMPL/LCFlagImpl.h"

#include "UTIL/CellIDEncoder.h"

#include "lcio/MLcio.hh"
#include "event/Event.hh"


using namespace std;
using namespace lcio ;

/*------------------------------------------------------------------------------------*/
void MLcio::setLCEventImp(LCEventImpl& evt, Event& evtRef)
{
try {

        evt.setRunNumber(  1   ) ;
    evt.setDetectorName( "m2" ) ;
    evt.setWeight( 1 ) ;


    evt.setEventNumber( evtRef.getId() ) ;
    // Event timestamp express in millisecond
        // But lcio timestamp express in second
    evt.setTimeStamp( evtRef.getTimeStamp() * 1000000  ) ;


    evt.parameters().setValue("Dif_syncro",evtRef.getDifSynchro() ) ;
    evt.parameters().setValue("Temperature",evtRef.getTemperature() ) ;
    evt.parameters().setValue("Pressure",evtRef.getPressure()) ;
    evt.parameters().setValue("Global trigger counter",(int)evtRef.getGlobalTriggerCounter() ) ;
    evt.parameters().setValue("Dif trigger counter",(int)evtRef.getDifTriggerCounter() ) ;
    //evt.parameters().setValue("SetDriftVoltage",evtRef.difTriggerCounter ) ;
    //evt.parameters().setValue("SetMeshVoltage",evtRef.difTriggerCounter ) ;



    // now add some raw calorimeter hits
    LCCollectionVec* calVec = new LCCollectionVec( LCIO::RAWCALORIMETERHIT )  ;

    LCFlagImpl chFlag(0) ;
    //chFlag.setBit( LCIO::CHBIT_ID1 ) ;

    chFlag.setBit( LCIO::RCHBIT_LONG )   ;  // long(1) - short(0) , incl./excl. position
    chFlag.setBit( LCIO::RCHBIT_BARREL ) ;  // barrel(1) - endcap(0)
    //chFlag.setBit( LCIO::RCHBIT_TIME   ) ; // 1: time information stored
    //chFlag.setBit( LCIO::RCHBIT_ID1   )  ;  // cellid1 stored
    //chFlag.setBit( LCIO::RCHBIT_NO_PTR ) ; // 1: pointer tag not added
    //chFlag.setBit( LCIO::RCHBIT_ENERGY_ERROR )  ;   // 1: store energy error
    calVec->setFlag( chFlag.getFlag()) ;


//
//OFFSET_CHAMBER *NB_CHANNEL* GetChamberId()
//+ OFFSET_DIF *NB_CHANNEL*  GetDifId()
//+ OFFSET_CHIP*NB_CHANNEL * GetChipId ()
//+  channelId 
//
        // Chambre 1024 -> 2 exp 10
        // Dif     128  -> 2 exp 7
        // Chip    1024 -> 2 exp 10
        // Channel 1024 -> 2 exp 10
    std::string cellIDEncoding( "CHAMBER:10,DIF:7,CHIP:10,CHANNEL:10") ;
    CellIDEncoder<RawCalorimeterHitImpl> cellEncoder( cellIDEncoding , calVec ) ;


    ChamberMap_t chamberList = evtRef.getRun().getDetector().getChambers();
    for (ChamberMap_t::iterator iter=chamberList.begin(); iter!=chamberList.end(); iter++)
    {
      Chamber& chamber = *(iter->second);
      try
      {
          const  multiset<ChannelHit*,sorter>& list = evtRef.getChannelHitVector(chamber);
          int order = 0 ;
          for (set<ChannelHit*,sorter>::iterator iter=list.begin(); iter!=list.end(); iter++)
          {
            ChannelHit& channelHit = *(*iter);
            const Channel& channel = channelHit.getChannel();
                RawCalorimeterHitImpl* rawHit = NULL ;
            try
            {
                          // Create RawCalorimeterHit
                  rawHit = new RawCalorimeterHitImpl() ;

              if ( chamber.getType().compare(GASSIPLEXCHAMBER1)==0 ||  chamber.getType().compare(GASSIPLEXCHAMBER4)==0)
              {
                    channelHit.setTimestamp(evtRef.getTimeStamp());
                          }
                        
                        
                  cellEncoder["CHAMBER"] = channel.getChamber().getId(); 
                  cellEncoder["DIF"] = channel.getChip().getBoard().getDif()->getId(); 
                  cellEncoder["CHIP"] = channel.getChip().getId(); 
                  cellEncoder["CHANNEL"] = channel.getHardId(); 
                          cellEncoder.setCellID( rawHit ) ;
                        
                  rawHit->setAmplitude( channelHit.getDigitalValue()   )  ;
                  //rawHit->setTimeStamp( channelHit.getTimestamp() / 1000) ; 

              order++;  // because aEvent have already been  sorted
                  calVec->push_back( rawHit ) ;
            } // end try block
            catch (MicroException e)
            {
            }

          } // end for channelHits for this chamber
      }
      catch (MicroException e)
      {
                // get ChamberHitVector for next chamber
      }

    }
    //LCTOOLS::printRawCalorimeterHits(calVec);
    // add all collections to the event

    evt.addCollection( calVec , "m2" ) ;


  }
  catch( Exception& ex){

    cout << " an excpetion occured: " << endl ;
    cout << "  " << ex.what() << endl ;
  }



}

---