Difference between revisions of "Analysis TTreeFormat"

//In plugin evnt()
vector<const DEventWriterROOT*> locEventWriterROOTVector;
locEventLoop->Get(locEventWriterROOTVector); //creates the TTrees for all DReactions upon first call
locEventWriterROOTVector[0]->Fill_Trees(locEventLoop);

//In plugin evnt()
vector<const DEventWriterROOT*> locEventWriterROOTVector;
locEventLoop->Get(locEventWriterROOTVector); //creates the TTrees for all DReactions upon first call
locEventWriterROOTVector[0]->Fill_Tree(locEventLoop, locReaction, locParticleCombos);

GetEntry(entry);
cout << ((TLorentzVector*)PiPlus2__P4_KinFit->At(0))->Px() << endl; //At(): particle combination 0, 1, 2, ...

MyTree->Draw("PiPlus2__P4_KinFit->Px()"); //draws all particle combinations

//Create the trees
vector<const DEventWriterROOT*> locEventWriterROOTVector;
locEventLoop->Get(locEventWriterROOTVector); //creates the TTrees for all DReactions upon first call
 
//Create the branch
double* locMyVariable = new double;
string locReactionName = locReaction->Get_ReactionName();
string locTreeName = locReactionName + string("_Tree");
japp->RootWriteLock(); //always acquire a lock before accessing the global ROOT file
{
  //get the tree
  gDirectory->cd("/");
  gDirectory->cd(locReactionName.c_str());
  TTree* locTree = (TTree*)gDirectory->Get(locTreeName.c_str());
  locTree->Branch("MyBranchName", locMyVariable, "D");
}
 
//Set the data in the branch and fill
*locMyVariable = 3.0; //or whatever you want to save
locEventWriterROOTVector[0]->Fill_Trees(locEventLoop);

//IDENTIFIERS
"ID": (unsigned int)[] //each particle has it's own #
"ParentID": (unsigned int)[] //the "Identifier" of the particle this particle decayed from (0 if none (e.g. photoproduced))
"PID_PDG": int[]
 
//KINEMATICS: THROWN  //At the production vertex 
"X4_Thrown": TClonesArray(TLorentzVector)
"P4_Thrown": TClonesArray(TLorentzVector)

//KINEMATICS: MEASURED //At the target center
"X4_Measured": TClonesArray(TLorentzVector)
"P4_Measured": TClonesArray(TLorentzVector)
 
//KINEMATICS: KINFIT //At the interaction vertex //only present if kinfit performed
"X4_KinFit": TClonesArray(TLorentzVector)
"P4_KinFit": TClonesArray(TLorentzVector)

//IDENTIFIERS
"ID": (unsigned int)[] //each physical particle has it's own # (to keep track of different pid hypotheses for the same particle)
 
//KINEMATICS: MEASURED //At the production vertex
"X4_Measured": TClonesArray(TLorentzVector) //t is the measured value in TOF/BCAL/FCAL projected back to Position_Measured
"P4_Measured": TClonesArray(TLorentzVector)
 
// KINEMATICS: END //at the reconstructed position of the BCAL/FCAL/TOF hit
"X4_End": TClonesArray(TLorentzVector)
"P4_End": TClonesArray(TLorentzVector)
 
//KINEMATICS: KINFIT //At the production vertex //only present if kinfit performed
"X4_KinFit": TClonesArray(TLorentzVector)
"P4_KinFit": TClonesArray(TLorentzVector)
 
// KINEMATICS: Other
"PathLength": double[] //from dPosition_KinFit/Measured (KinFit if performed) to dPosition_End
 
// PID QUALITY:
"NDF_Tracking": (unsigned int)[] //for charged only
"ChiSq_Tracking": double[] //for charged only
"NDF_Timing": (unsigned int)[]
"ChiSq_Timing": double[] //using kinematic fit data
"ChiSq_Timing_Measured": double[] //using measured data
"NDF_DCdEdx": (unsigned int)[] //for charged only
"ChiSq_DCdEdx": double[] //for charged only
 
// DEPOSITED ENERGY:
"dEdx_CDC": double[] //for charged only
"dEdx_FDC": double[] //for charged only
"dEdx_TOF": double[] //for charged only
"dEdx_ST": double[] //for charged only
"Energy_BCAL": double[]
"Energy_FCAL": double[]

//IDENTIFIERS
"ID": (unsigned int)[] //each physical particle has it's own # (to keep track of different pid hypotheses for the same particle)
"PID": (unsigned int)[] //Particle_t value
 
//KINEMATICS: MEASURED  //At the production vertex 
"P4_Measured": TClonesArray(TLorentzVector)
"X4_Measured": TClonesArray(TLorentzVector) //t is the measured value in TOF/BCAL/FCAL projected back to Position_Measured
 
// KINEMATICS: END //at the reconstructed position of the BCAL/FCAL/TOF hit
"P4_End": TClonesArray(TLorentzVector)
"X4_End": TClonesArray(TLorentzVector)
 
// KINEMATICS: Other
"PathLength": double[] //from Position_Measured to Position_End
 
// PID QUALITY:
"NDF_Tracking": (unsigned int)[] //for charged only
"ChiSq_Tracking": double[] //for charged only
"NDF_Timing": (unsigned int)[]
"ChiSq_Timing": double[] //using kinematic fit data
"ChiSq_Timing_Measured": double[] //using measured data
"NDF_DCdEdx": (unsigned int)[] //for charged only
"ChiSq_DCdEdx": double[] //for charged only
 
// DEPOSITED ENERGY:
"dEdx_CDC": double[] //for charged only
"dEdx_FDC": double[] //for charged only
"dEdx_TOF": double[] //for charged only
"dEdx_ST": double[] //for charged only
"Energy_BCAL": double[]
"Energy_FCAL": double[]

class DTreeParticle : public TObject
{
  public:
    // PID:
    Particle_t dPID;
 
    // KINEMATICS: Measured
    TVector3 dPosition_Measured; //the position where the particle is produced
    double dTime_Measured; //time of the track at dPosition_Measured, projected from the TOF/BCAL/FCAL hit
    TVector3 dMomentum_Measured; //momentum of the track at dPosition_Start
 
    // KINEMATICS: KinFit 
    TVector3 dPosition_KinFit; //the position where the particle is produced
    double dTime_KinFit; //time of the track at dPosition_KinFit
    TVector3 dMomentum_KinFit; //momentum of the track at dPosition_KinFit
 
    // KINEMATICS: End
    TVector3 dPosition_End; //detected particles: the reconstructed position of the BCAL/FCAL/TOF hit; decaying particles: the point where it decays
    double dTime_End; //time of the track at dPosition_End
    TVector3 dMomentum_End; //momentum of the track at dPosition_End
 
    // KINEMATICS: Other
    double dPathLength; //from dPosition_KinFit/Measured (if kinfit performed/not-performed) to dPosition_End
 
    // PID QUALITY:
    unsigned int dNDF_Tracking;
    double dChiSq_Tracking;
    unsigned int dNDF_Timing;
    double dChiSq_Timing;
    unsigned int dNDF_DCdEdx;
    double dChiSq_DCdEdx;
 
    // DEPOSITED ENERGY:
    map<DetectorSystem_t, double> dDepositedEnergies; //Is dE/dx for all but BCAL/FCAL (shower energies)
 
    // DTREESTEP POINTERS:
    const DTreeStep* dProductionStep; //the step object in which this DTreeParticle is produced (is a final-state particle)
    const DTreeStep* dDecayStep; //the step object in which this DTreeParticle decays (is an initial-state particle) (will be null if not a decaying particle!)
 
    // CUSTOM VARIABLES:
    map<string, double> dCustomVariables; //key is unique string, double is value
    map<string, const TObject*> dCustomObjects; //key is unique string, TObject* is object
 
  ClassDef(DTreeParticle, 1)
};

class DTreeStep : public TObject
{
  public:
    // INITIAL PARTICLES:
    const DTreeParticle* dInitialParticle; //if is null: decaying or beam particle not yet set!
    const DTreeParticle* dTargetParticle; //NULL for no target
 
    // FINAL PARTICLES:
    deque<DParticle_t> dFinalParticleIDs; //separate in case particle is NULL (e.g. decaying resonance)
    deque<const DTreeParticle*> dFinalParticles; //particle may be NULL if it is decaying or missing (especially if no kinematic fit was performed!!)
 
    // CUSTOM VARIABLES:
    map<string, double> dCustomVariables; //key is unique string, double is value
    map<string, const TObject*> dCustomObjects; //key is unique string, TObject* is object
 
  ClassDef(DTreeStep, 1)
};

class DTreeCombo : public TObject
{
  public:
    // STEPS:
    deque<const DTreeStep*> dTreeSteps;
 
    // RF:
    double dRFTime;
 
    // UNUSED PARTICLES:
    vector<const DTreeParticle*> dUnusedDetectedParticles;
    vector<const DTreeShower*> dUnusedDetectedShowers;
 
    // KINEMATIC FIT:
    DKinFitType dKinematicFitType; //Defined in DKinFitResults.h //d_NoFit if not performed
    double dChiSq_KinematicFit; //NaN if not performed
    unsigned int dNDF_KinematicFit; //0 if not performed
 
    // CUSTOM VARIABLES:
    map<string, double> dCustomVariables; //key is unique string, double is value
    map<string, const TObject*> dCustomObjects; //key is unique string, TObject* is object
 
  ClassDef(DTreeCombo, 1)
};

class DTreeEvent : public TObject
{
  public:
    // RUN, EVENT #'s:
    unsigned int dRunNumber;
    unsigned int dEventNumber;
 
    // DATA:
    map<string, deque<const DTreeCombo*> > dTreeCombos; //string key is (D)Reaction name, deque is the particle combos
    deque<const DTreeParticle*> dThrownParticles;
 
    // CUSTOM VARIABLES:
    map<string, double> dCustomVariables; //key is unique string, double is value
    map<string, const TObject*> dCustomObjects; //key is unique string, TObject* is object
 
  ClassDef(DTreeEvent, 1)
};